Python與開源GIS：根據(jù)空間過濾器選擇要素參考模板

1、Python與開源GIS：根據(jù)空間過濾器選擇要素2012-12-03 21:12 作者: bukunu 來源: 本站 瀏覽: 98 views 我要評論 字號: 大 中 小摘要: 空間過濾器Spatial filters 如果說按照屬性進行要素的選擇，還是帶有數(shù)據(jù)庫的特征的話，那么，根據(jù)空間位置進行選擇，則是地地道道的GIS功能。在OGR中，使用了Spatial filters（空間過濾）這一術(shù)語表征這一功能。 OGR提供的空間過濾功能.空間過濾器Spatial filters如果說按照屬性進行

2、要素的選擇，還是帶有數(shù)據(jù)庫的特征的話，那么，根據(jù)空間位置進行選擇，則是地地道道的GIS功能。在OGR中，使用了Spatial filters（空間過濾）這一術(shù)語表征這一功能。OGR提供的空間過濾功能有兩種，一種是SetSpatialFilter()，過濾某一類型的Feature，例如參數(shù)中使用Polygon，就是選出Layer中的所有Polygon所覆蓋的要素（注意，只要相交即可，不必完全包含）。下面這段代碼用了兩套數(shù)據(jù)。world_borders是全球國界數(shù)據(jù)，cover.shp則是覆蓋了非洲南部地區(qū)的一個多邊形。下面先定義一個根據(jù)圖層直接生成shape文件的函數(shù)，方便后面調(diào)用。123456

3、from osgeo import ogrdef create_shp_by_layer(shp, layer):    outputfile = shp    if os.access(outputfile, os.F_OK):    driver.DeleteDataSource(outputfile)    newds = driver. CreateDataSource ( outputfile )1 / 378 &#

4、160;  pt_layer = newds.CopyLayer ( layer, '')    newds.Destroy ()下面代碼是使用cover.shp中的多邊形來選擇全球國界數(shù)據(jù)：1234567891011121314driver = ogr.GetDriverByName("ESRI Shapefile")world_shp = '/gdata/world_borders.shp'cover_shp = '/gdata/cover.shp'world_d

5、s = ogr.Open(world_shp)cover_ds = ogr.Open(cover_shp)world_layer = world_ds.GetLayer(0)cover_layer = cover_ds.GetLayer(0)print(world_layer.GetFeatureCount()cover_feats = cover_layer.GetNextFeature() poly = cover_feats.GetGeometryRef()world_layer.SetSpatialFilter(poly)out_shp = '/gdata/world_cover.shp'create_shp_by_layer(out_shp, world_layer)結(jié)果可以通過下面的圖來看。另外還有SetSpatialFilterRect(minx, miny, maxx, maxy)，參數(shù)輸入四個坐標，可以選中矩形內(nèi)的要素。123world_layer.SetSpatialFilterRect(50, 60, 25, 35)out_shp = '/gdata/world_spatial_filter.shp'create_shp_by_layer(