Python gdal チュートリアル: フィルター、単純な空間分析、関数とモジュール

Layer オブジェクトには SetAttributeFilter() というメソッドがあり、特定の条件を満たすレイヤー内のフィーチャをフィルターで除外できます。 Filter を設定した後、GetNextFeature() メソッドを使用して、条件を満たすフィーチャーを順番に取得できます。 SetAttributeFilter(None) はフィルターをクリアできます。たとえば

>>>layer.GetFeatureCount()

>>>layer.SetAttributeFilter("cover = 'shrubs'")

>>> ) 空間フィルターには 2 種類あります。 1 つは SetSpatialFilter() で、特定の種類の地物をフィルターします。たとえば、パラメーターに Polygon を入力すると、レイヤー内のすべての Polygon が選択されます

SetSpatialFilterRect(, ) もあります。 ;、< ;maxx>、)、パラメーターとして 4 つの座標を入力し、ボックスで [Feature

SetSpatialFilter(None)] を選択して空間属性フィルターをクリアできます。

たとえば、次のコードでは、layerAreas はポリゴンであり、layerSites はポイントです

>>> featAreas = LayerAreas.GetNextFeature()

>>> Poly = featAreas.GetGeometryRef()

> ;> & gt; レイヤーサイト.getfeatureCount ()

& gt; & gt; gt; & gt;layersites.getfeatureCount ()

33

>>>layerSites.SetSpatialFilterRect(460000, 4590000, 490000, 4600000)

>> layerSites.GetFeatureCount()

4

>>>layerSites.SetSpatialFilter(None)

>>>layerSites.GetFeatureCount()

42

SQL クエリ ステートメントの実行など、より複雑なフィルターもあります。 SQL の強力な機能を使用すると、より複雑なフィルターを実行できます。次のコードのようなタスクは、カバー タイプが草であるフィーチャを選択し、ID 番号の降順に並べ替えることです。

result = dsSites.ExecuteSQL("select * from sites where cover = 'grass' order by id desc")

resultFeat = result.GetNextFeature()

while resultFeat :

print resultFeat.GetField('id') print resultFeat.GetField('id')

resultFeat = result.GetNextFeature()

dsSites.ReleaseResultSet(result)

42