Python与开源GIS：使用OGR获取要素信息

获取要素(Feature)信息

如果只从计算机角度来看，有人认为要素就是一些几何形状， 这种观点是欠妥的。几何形状，具体包括点、线、多边形、 弧段、向量、控制点等等非常多种，其实是要素的模型抽 象。 在Shapefile中，要素模型由点、线、面三种类型构成。 要素类都带有属性信息，一个要素一般对应属性表中的一行。

获取图层中的要素

下面看一下如何获取图层中的要素（feature）。

 >>> from osgeo import ogr >>> inshp = '/gdata/world_borders.shp' >>> datasource = ogr.Open(inshp) >>> layer = datasource.GetLayer(0) >>> feature = layer.GetFeature(0) 

另外还有按顺序读取feature，循环遍历所有的feature：

 >>> feat = layer.GetNextFeature() >>> while feat: > feat = layer.GetNextFeature() >>> layer.ResetReading() 

GetNextFeature()用来获取下一下要素，ResetReading()则 是进行重置，以便再次获取。

获取要素的属性

下面是一个简单的例子：

 >>> feat = layer.GetFeature(0) >>> feat.keys() ['CAT', 'FIPS_CNTRY', 'CNTRY_NAME', 'AREA', 'POP_CNTRY'] >>> fid = feat.GetField('AREA') >>> print(fid) 193.0 

下面是对所有的属性值进行遍历：

 >>> for i in range(feature.GetFieldCount()): > print feature.GetField(i) 1.0 AA Aruba 193.0 71218.0 

这会列出这个要素的所有属性值。 如果要看整个表的结构， 各个字段的名称等等信息，可以在layer的附加信息里看。

 >>> layerdef = layer.GetLayerDefn() >>> for i in range(layerdef.GetFieldCount()): > defn = layerdef.GetFieldDefn(i) > print(defn.GetName(),defn.GetWidth(),defn.GetType(),defn.GetPrecision()) ('CAT', 16, 2, 0) ('FIPS_CNTRY', 80, 4, 0) ('CNTRY_NAME', 80, 4, 0) ('AREA', 15, 2, 2) ('POP_CNTRY', 15, 2, 2) 

要素的形状(Geometry)

 >>> geom = feature.GetGeometryRef() >>> geom.GetGeometryName() 'POLYGON' >>> geom.GetGeometryCount() 1 >>> geom.GetPointCount() 0 >>> geom.GetX() 0.0 >>> geom.GetY() 0.0 >>> print(geom) # 会打印出所有的点 POLYGON ((-69.882232999999999 12.41111,-69.946944999999999 12.436666, ... 12.411664999999999,-69.882232999999999 12.41111)) >>> print(geom.ExportToWkt()) POLYGON ((-69.882232999999999 12.41111,-69.946944999999999 12.436666, ... 12.411664999999999,-69.882232999999999 12.41111)) >>> polygon = geom.GetGeometryRef(0) >>> polygon.GetGeometryName() 'LINEARRING' >>> polygon.GetGeometryCount() 0 >>> polygon.GetPointCount() 19 >>> polygon.GetX(0) -69.882232999999999 >>> polygon.GetY(0) 12.411110000000001 >>> polygon.GetZ(0) 0.0 >>> print(polygon.ExportToWkt()) LINEARRING (-69.882232999999999 12.41111,-69.946944999999999 12.436666, ... 12.411664999999999,-69.882232999999999 12.41111) 

返回目录：Python与开源GIS