首先,我很抱歉发布这个简单的问题。我需要计算一定数量的几何属性(面积、周长、圆度、长轴和短轴等)。我正在使用GDAL/OGR读取我的多边形的 shapefile 格式。我想问的是:
- 有没有使用 osgeo.ogr.Geometry 计算周长的方法?
- 是否有用于计算多边形指标的模块?
提前致谢
import osgeo.gdal, ogr
poly="C:\\\myshape.shp"
shp = osgeo.ogr.Open(poly)
layer = shp.GetLayer()
# For every polygon
for index in xrange(len(allFID)):
feature = layer.GetFeature(index)
# get "FID" (Feature ID)
FID = str(feature.GetFID())
geometry = feature.GetGeometryRef()
# get the area
Area = geometry.GetArea()