### 使用Python创建Shapefile文件 在GIS二次开发过程中，创建Shapefile是一项常见需求。可以利用`osgeo.ogr`库来完成这项工作。下面展示了一个完整的例子，说明如何使用Python编写程序创建一个新的Shapefile文件。 #### 准备环境 为了能够顺利运行以下代码片段，需提前安装GDAL库，该库包含了处理矢量数据所需的工具集。可以通过pip命令轻松安装： ```bash pip install gdal ``` #### 创建Shapefile的具体过程 ```python from osgeo import ogr, osr # 设置输出路径及文件名 output_path = "path/to/output" shp_filename = f"{output_path}/new_shapefile.shp" # 注册所有的驱动 driver_name = "ESRI Shapefile" drv = ogr.GetDriverByName(driver_name) # 如果存在同名文件则删除之 if drv is None: raise ValueError(f"无法找到名为 {driver_name} 的驱动") if os.path.exists(shp_filename): drv.DeleteDataSource(shp_filename) # 创建新的Shapefile数据源 data_source = drv.CreateDataSource(shp_filename) spatial_reference = osr.SpatialReference() spatial_reference.ImportFromEPSG(4326) # WGS84坐标系 # 建立图层（此处假设为点状要素） layer = data_source.CreateLayer("points", spatial_reference, geom_type=ogr.wkbPoint) # 添加字段定义 field_defn = ogr.FieldDefn("Name", ogr.OFTString) field_defn.SetWidth(50) layer.CreateField(field_defn) # 向图层中加入新特征 feature_definition = layer.GetLayerDefn() for i in range(1, 4): # 插入三个样本点 point_x = float(i * 10 + 79) point_y = float(i * 10 + 20) # 构建几何对象 point = ogr.Geometry(ogr.wkbPoint) point.AddPoint(point_x, point_y) # 新建Feature实例并设置属性值 feature = ogr.Feature(feature_definition) feature.SetGeometry(point) feature.SetField("Name", f"point_{i}") # 将新建的Feature写入到图层里 layer.CreateFeature(feature) # 清理不再使用的资源 feature.Destroy() # 关闭数据源连接 data_source.Destroy() print(f'Shapefile已成功保存至{shp_filename}') ``` 上述代码实现了基于给定坐标的若干个点要素组成的简单Shapefile文件创建[^1]。