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

2017-04-03 作者: bukun 浏览: 390 次 Help edit

摘要: 图层的概念 这里首先解释一下图层(Layer)的概念。 在GIS中,Layer是指是一个由同种要素(Feature,如点、线、多边形等)组合在一起的“层”。这个概念,ESRI在ArcGIS中定义的模型中的要素类(FeatureClass)。 在《Modelin....

图层的概念

这里首先解释一下图层(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()。

返回目录:Python与开源GIS

随机推荐

Copyright © 2014-2017 OSGeo中国中心 吉ICP备05002032号

Powered by TorCMS