收集平时遇到的GDAL/OGR开发问题(不定时更新中……)
1、使用OGR的图层筛选Layer.SetAttributeFilter()和Layer.SetSpatialFilter(),筛选的结果使用
for feat in Layer遍历,不能直接使用Layer.GetFeature()获取筛选后的要素
2、计算Geometry面积和周长
若Geometry为Polygon,先需要获取边界,再计算周长Length();面积则直接计算Area()。
3、如果Geometry为MultiPolygon,需要注意对其polygon进行处理
for subGeometry in Geometry:
...
4、使用ospybook.vectorplotter.VectorPlotter绘图,绘制点的时候,注意传入点样式的参数,不然会报错。例如,点要素图层
from ospybook.vectorplotter import VectorPlotter
v