在处理从 `.shp` 文件中按“宗地代码”导出单个 `.dwg` 文件，并以“宗地代码 + 权利人姓名”命名的任务中，可以使用 Python 结合 `geopandas` 和 `pyshp`（`shapefile`）等库读取 `.shp` 文件，再通过 `ezdxf` 库创建 `.dwg` 文件。以下是一个完整的实现思路和示例代码。 ### 实现步骤 1. **读取 `.shp` 文件**：使用 `geopandas` 或 `shapefile` 读取属性数据和几何信息。 2. **筛选数据**：根据“宗地代码”字段分组或筛选。 3. **创建 `.dwg` 文件**：使用 `ezdxf` 创建新的 DXF 文档，并将对应的几何图形写入。 4. **命名规则**：以“宗地代码 + 权利人姓名”作为文件名保存。 ### 示例代码 ```python import geopandas as gpd import os from shapely.geometry import shape import ezdxf # 加载 shp 文件 gdf = gpd.read_file("path/to/your/file.shp") # 确保字段名正确 gdf.columns = gdf.columns.str.lower() required_columns = ['宗地代码', '权利人姓名'] # 遍历每一行数据 for index, row in gdf.iterrows(): 宗地代码 = row['宗地代码'] 权利人姓名 = row['权利人姓名'] geometry = row['geometry'] # 创建一个新的 DXF 文档 doc = ezdxf.new(dxfversion='R2010') msp = doc.modelspace() # 判断几何类型并写入 if geometry.geom_type == 'Polygon': coords = list(geometry.exterior.coords) points = [(float(x), float(y)) for x, y in coords] msp.add_lwpolyline(points, close=True) elif geometry.geom_type == 'MultiPolygon': for poly in geometry.geoms: coords = list(poly.exterior.coords) points = [(float(x), float(y)) for x, y in coords] msp.add_lwpolyline(points, close=True) # 构建文件名 filename = f"{宗地代码}_{权利人姓名}.dwg" output_path = os.path.join("output_folder", filename) # 保存 DXF 文件 doc.saveas(output_path) ``` ### 注意事项 - 确保字段名正确，必要时进行字段重命名或转换为小写处理。 - 若 `.shp` 文件包含多个几何对象（如多边形、点等），需根据实际类型进行判断和写入。 - `ezdxf` 支持 `.dxf` 格式，若需 `.dwg` 格式输出，可在保存后使用第三方工具（如 AutoCAD）进行转换。 ###