8_读取矢量数据

最新推荐文章于 2024-08-05 17:56:51 发布

原创

最新推荐文章于 2024-08-05 17:56:51 发布 · 619 阅读

7 ·

CC 4.0 BY-SA版权

文章标签：

#python #开发语言

本文详细介绍了如何使用Python的osgeo.ogr模块读取和操作ESRI Shapefile矢量数据，包括获取图层、图层属性、要素信息、几何形状以及矢量文件的投影参数。示例代码演示了从打开shapefile到获取特征的属性和几何信息的完整过程，并展示了如何遍历和提取属性表内容。最后，文章提到了释放资源的重要性。

读取矢量数据

读取矢量数据的方式就是上一章所讲的两种方法。

操作图层

现在要明白矢量数据的文件组织结构。

shp—>layer—>feature—>geometry—>wkt

一个shp文件中是只包含一个layer的；

layer就是同一个要素的集合，比如点集合、线集合、面集合，一个城市所有的学校为点集合；

feature指代单个的点线面之类的geometry，比如学校点集合中的某一个学校就是一个要素，在属性表中表示为一行数据；

geometry的组织方式就是wkt。

__author__ = 'He XK'
from osgeo import ogr

if __name__ == '__main__':
    shp_filename = 'shp/province_1.shp'
    driver = ogr.GetDriverByName('ESRI Shapefile')
    shp = driver.Open(shp_filename)
    print(shp.GetLayerCount())  # 1
    # 获取图层
    layer = shp.GetLayer(0)
    # 获取图层要素
    print(layer.GetFeatureCount())  # 37
    print(layer.GetGeomType())  # 3
    print(layer.GetExtent())  # (-3349547.371448596, 1834629.4342002147, 699750.9271820501, 5953688.75950143)