python求周长_Python如何使用osgeo.ogr.Geometry对象计算多边形周长

First of all, I apologize to post this easy question. I need to compute a certain number of gemotrical attributes (area, perimeters, Roundess, major and minor axis, etc). I am using GDAL/OGR to read a shapefile format of my polygon. What i wish to ask is:

is there a method to compute the perimeter using osgeo.ogr.Geometry?

is there a module build to compute metrics on polygon?

thanks in advance

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()

解决方案poly = [(0,10),(10,10),(10,0),(0,0)]

def segments(poly):

"""A sequence of (x,y) numeric coordinates pairs """

return zip(poly, poly[1:] + [poly[0]])

def area(poly):

"""A sequence of (x,y) numeric coordinates pairs """

return 0.5 * abs(sum(x0*y1 - x1*y0

for ((x0, y0), (x1, y1)) in segments(poly)))

def perimeter(poly):

"""A sequence of (x,y) numeric coordinates pairs """

return abs(sum(math.hypot(x0-x1,y0-y1) for ((x0, y0), (x1, y1)) in segments(poly)))