Python与开源GIS:获取图层信息

图层的概念

这里首先解释一下图层(Layer)的概念。 在GIS中,Layer是指是一个由同种要素(Feature,如点、线、多边形等)组合在一起的“层”。这个概念,ESRI在ArcGIS中定义的模型中的要素类(FeatureClass)。 在《Modeling our World》中,对这两个概念有明确的阐述:

  1. 要素数据集(要素集):是具有相同坐标系统的要素类的集合。我们可以选择在要素集的内部或外部组织简单要素类, 但拓扑要素类只能在要素集内部组织,以确保它们具有相同的坐标系统。
  2. 要素类:要素类是具有相同几何形状的要素的集合:点、线或多边形。 我们最关心的两种要素类是简单要素类和拓扑要素类。简单要素类包括没有任何拓扑关系的点、线、多边形或注记。也就是说,一个要素类内的点与另一要素类中的线的终点可以是一致的,但它们是不同的。这些要素可以彼此独立地编辑。拓扑要素类局限在一定的图形范围内,它是一个由完整拓扑单元组成的一组要素类限定的对象。

Python中操作Layer的方法

下面是Layer的方法集合:

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

这里要注意 GetLayer 的参数是从 0 开始的。对于Shapefile,它只有一个图层,一般都是填 0 的,如果不填的话默认也是 0 。下面是layer的主要方法:

'CommitTransaction', 'CreateFeature', 'CreateField', 'DeleteFeature', 'Dereference', 'GetExtent', 'GetFeature', 'GetFeatureCount', 'GetFeaturesRead', 'GetLayerDefn', 'GetName', 'GetNextFeature', 'GetRefCount', 'GetSpatialFilter', 'GetSpatialRef', 'Reference', 'ResetReading', 'RollbackTransaction', 'SetAttributeFilter', 'SetFeature', 'SetNextByIndex', 'SetSpatialFilter', 'SetSpatialFilterRect', 'StartTransaction', 'SyncToDisk', 'TestCapability'。

下面是一些简单的例子:

n = layer.GetFeatureCount()
print 'feature count:', n
# 读出上下左右边界
extent = layer.GetExtent()
print 'extent:', extent
print 'ul:', extent[0], extent[3]
print 'lr:', extent[1], extent[2]

可以看到,Layer核心操作都是针对的是要素的操作(Feature),并且通过 GetFeatureCount()GetFeature() 或者 GetNextFeature() 可以获取层内所有的要素,通过 GetFeaturesRead() 可以知道现在已经读取了多少条 Feature 了(获取一次值少1,如果要从头开始,用 ResetReading()