GDAL/OGR空间过滤之SetSpatialFilterRect()

• SetSpatialFilterRect函数，在osgeo.ogr模块下的Layer类中。也即，在图层（Layer）层面上执行空间筛选。
• 通过设置一个矩形范围，再进行空间筛选，这通过函数名称中Rect就可以看出来。
• SetSpatialFilterRect(double minx, double miny, double maxx, double maxy)，即分别设置矩形的最小x，最小y，最大x，最大y。x、y的值应当与图层在相同的坐标系下设定。
• 返回与给定矩形在几何上相交的要素。
• 示例程序如下：

from osgeo import ogr

g = ogr.Open('world.shp')
lyr = g.GetLyaer(0)
lyr.SetSpatialFilterRect(-14, 37, 38, 72)
for feat in lyr:
  print(feat.GetFieldAsString('NAME') + '------'+feat.GetFieldAsString('CAPITAL'))
lyr.ResetReading()

ALGERIA ---- ALGIERS
BELGIUM ---- BRUSSELS
LUXEMBOURG ---- LUXEMBOURG
SAN MARINO ---- SAN MARINO
AUSTRIA ---- VIENNA
CZECH REPUBLIC ---- PRAGUE
SLOVENIA ---- LJUBLJANA
HUNGARY ---- BUDAPEST
SLOVAKIA ---- BRATISLAVA
YUGOSLAVIA ---- BELGRADE [BEOGRADE]
BOSNIA AND HERZEGOVINA ---- SARAJEVO
ALBANIA ---- TIRANE
MACEDONIA, THE FORMER YUGOSLAV REPUBLIC ---- SKOPJE
LITHUANIA ---- VILNIUS
LATVIA ---- RIGA
BULGARIA ---- SOFIA
BELARUS ---- MINSK
MOLDOVA, REPUBLIC OF ---- KISHINEV
IRELAND ---- DUBLIN
ICELAND ---- REYKJAVIK
SPAIN ---- MADRID
SWEDEN ---- STOCKHOLM
FINLAND ---- HELSINKI
TURKEY ---- ANKARA
RUSSIAN FEDERATION ---- MOSCOW
GREECE ---- ATHENS
PORTUGAL ---- LISBON
POLAND ---- WARSAW
NORWAY ---- OSLO
GERMANY ---- BERLIN
ESTONIA ---- TALLINN
TUNISIA ---- TUNIS
CROATIA ---- ZAGREB
ROMANIA ---- BUCURESTI
UKRAINE ---- KIEV
NETHERLANDS ---- AMSTERDAM
JERSEY ---- SAINT HELIER
GUERNSEY ---- SAINT PETER PORT
FAROE ISLANDS ---- TORSHAVN
DENMARK ---- COPENHAGEN
MONACO ---- MONACO
ANDORRA ---- ANDORRA LA VELLA
LIECHTENSTEIN ---- VADUZ
SWITZERLAND ---- BERN
ISLE OF MAN ---- DOUGLAS
UNITED KINGDOM ---- LONDON
FRANCE ---- PARIS
VATICAN CITY (HOLY SEE) ---- VATICAN CITY
ITALY ---- ROME

• 要清除空间过滤集的唯一方式是调用SetSpatialFilter()函数。

lyr.SetSpatialFilterRect(NULL)