一、使用qgis開源工具構建兩個文件,分別是filegdb和shp,每個文件包含一個圖層,圖層內容只包含一個字段:id,有兩個數據行,圖層幾何為多邊形,圖層都是如下的效果。
二、使用rust讀取上述文件
- rust依賴
[dependencies]
gdal = "0.17.1"
- 讀取shp
let source = "E:/a.shp";
let data = gdal::Dataset::open(source).unwrap();
let layer = data.layer(0).unwrap();
println!("layer_count:{}",layer.feature_count());
let start = 0;
let end = layer.feature_count();
for i in start..end{match layer.feature(i){Some(feature) => {let filed_id= &feature.field_as_integer_by_name("id").unwrap();println!("++數據索引:{}++字段內容:{:?}++幾何:{:?}",i,filed_id,feature.geometry());},None => {println!("++{}++None",i);},};
}
讀取結果:
layer_count:2
++數據索引:0++字段內容:Some(1)++幾何:Some(POLYGON ((-0.518043358996392 0.521232057117657,4.73046834618934 0.672511281341224,5.00976830389323 -2.57346955999557,-0.808980814937951 -2.73621961134104,-0.518043358996392 0.521232057117657),(0.925006422473741 -0.235197478165903,0.622431468294519 -1.95716715064606,3.67145600656206 -1.66637545113517,3.63654351184907 -0.188647804914247,0.925006422473741 -0.235197478165903),(-0.180555910104184 -0.561039123478221,0.215119029976337 0.044101814636796,-0.192193408341846 -0.060635653828846,-0.180555910104184 -0.561039123478221),(-0.285293394243145 -1.67800799341694,-0.366755881906781 -2.29442185095562,0.57588147534387 -2.34093376312268,-0.285293394243145 -1.67800799341694),(4.34643090434648 -1.49187934036916,3.69473100303738 -2.35256150119173,4.43953089024778 -2.35256150119173,4.61409336381271 -1.86411893486817,4.742105844427 -1.44534463965808,4.47444338496076 -1.22429231197673,4.48608088319843 -0.700680602759594,4.71883084795167 -0.700680602759594,4.1835059290192 -0.654133864500525,4.34643090434648 -1.49187934036916)))
++數據索引:1++字段內容:Some(2)++幾何:Some(POLYGON ((0.749541894170621 -3.51691666514472,3.61879402562659 -5.06406242230235,0.81986670131415 -5.02186753801624,0.749541894170621 -3.51691666514472)))
從結果看到,讀取shp文件的數據行索引是從 0 開始的
- 讀取filegdb
let source = "E:/a_gdb.gdb";
let data = Dataset::open(source).unwrap();
let mut layer = data.layer_by_name("a").unwrap();
println!("layer_count:{}",layer.feature_count());let start = 1;
let end = layer.feature_count()+1;
for i in start..end{match layer.feature(i){Some(feature) => {let filed_id= &feature.field_as_integer_by_name("id").unwrap();println!("++數據索引:{}++字段內容:{:?}++幾何:{:?}",i,filed_id,feature.geometry());},None => {println!("++{}++None",i);},};
}
讀取結果:
layer_count:2
++數據索引:1++字段內容:Some(1)++幾何:Some(MULTIPOLYGON (((-0.518043358999989 0.52123205700002,4.73046834600001 0.672511281000027,5.00976830399998 -2.57346955999998,-0.808980814999984 -2.73621961100002,-0.518043358999989 0.52123205700002),(0.925006422000024 -0.235197477999975,0.622431468 -1.95716715100002,3.67145600700002 -1.66637545100002,3.636543512 -0.188647805000016,0.925006422000024 -0.235197477999975),(-0.18055591000001 -0.561039123,0.215119029999983 0.044101814999976,-0.19219340799998 -0.060635654,-0.18055591000001 -0.561039123),(-0.285293394 -1.67800799299999,-0.366755882 -2.29442185099998,0.575881475000017
-2.34093376300001,-0.285293394 -1.67800799299999),(4.34643090399999 -1.49187934000003,3.69473100300002 -2.35256150100003,4.43953089000001 -2.35256150100003,4.61409336399998 -1.86411893500002,4.74210584399998 -1.44534463999997,4.47444338499997 -1.22429231199999,4.486080883 -0.700680603000023,4.71883084799998 -0.700680603000023,4.18350592899998 -0.654133865,4.34643090399999
-1.49187934000003))))
++數據索引:2++字段內容:Some(2)++幾何:Some(MULTIPOLYGON (((0.749541894 -3.516916665,3.61879402599999 -5.06406242200001,0.819866701000024 -5.02186753799998,0.749541894 -3.516916665))))
從讀取結果中看出,filegdb的數據行索引是從 1 開始的
三、總結
讀取這兩種不同的數據源,方法總體上是差不多的,需要在遍歷數據的時候,注意數據行起始索引的區別。