使用 gdal 读取 shp 文件

简介：
本教程详细介绍如何利用GDAL库在Python环境中高效地打开和处理SHP格式文件，助力地理数据科学入门者掌握基础操作。 使用Python编程读取shp文件需要借助gdal库。首先，在电脑上安装gdal。下面是一个带有注释的示例代码： ```python from osgeo import ogr # 打开.shp 文件 shapefile = path_to_your_shapefile.shp dataset = ogr.Open(shapefile) # 获取图层数量，一般shp文件只有一个图层 layer_count = dataset.GetLayerCount() for layer_index in range(layer_count): # 获取每个图层对象 layer = dataset.GetLayerByIndex(layer_index) # 打印当前处理的图层名称和要素数量 print(fProcessing Layer: {layer.GetName()}) feature_count = layer.GetFeatureCount() print(fNumber of features in this layer: {feature_count}) for i in range(feature_count): # 获取每个要素（即shp文件中的一个记录） feature = layer.GetNextFeature() # 打印要素属性 if feature: attribute_names = [field.name for field in feature.schema] print(fAttributes of Feature {i + 1}:) for name in attribute_names: print(f{name} : {feature[name]}) # 关闭数据集，释放资源 dataset.Destroy() ``` 以上代码展示了如何使用gdal库在Python中读取.shp文件中的图层和要素信息。请根据实际情况修改`path_to_your_shapefile.shp`为实际的shp文件路径。