简介上海市县区乡镇村各级点位Shapefile数据包专注于提供全市范围内区县、乡镇、村三级行政单位位置点信息适用于GIS空间分析、区域对比、专题制图及智慧城市相关研究。面向地理信息开发人员、数据分析师及城市规划研究者可快速用于点位展示、空间查询与批量匹配等场景。资源共24个文件以.shp几何数据、.dbf属性表、.shx索引、.prj坐标参考为核心另含.sbn/.sbx空间索引及.xml元数据文件覆盖三类图层数据完整、结构清晰包体仅460KB便于下载与处理。目前已有173人学习浏览可作为中小尺度区域分析的基础底图数据。通过该数据包用户既能获得各级点位的地理坐标又能结合属性字段开展分级统计、热点识别等深入分析有效提升区域研究工作与项目管理效率。1. 为什么拿到这套上海市县区乡镇村点位 SHP先别急着拖进 ArcMap做上海范围内的点位分析与区域分析最怕的不是没有数据而是坐标系混乱、字段乱码、边界对不上。这套点位包把市级、区县级、乡镇级、村级四个层级的点位置放在同一套 Shapefile 文件里短文件名直接表明层级属性表结构相对一致适合直接用于选址、可达性分析、商圈划分和 POI 质量校验。但 SHP 不是单个文件而是一组配套文件.shp只是几何.dbf才是属性少一个.prj就不知道经纬度还是投影坐标。下面按文件结构、字段清洗、实际分析、批量转换的顺序把整套流程拆开。2. Shapefile 文件家族解密从 .shp 到 .sbx每个文件都别乱删2.1 一套完整 SHP 数据的文件构成与加载时依赖顺序解压后你会看到区县点.shp、乡镇点.shp、村点.shp以及对应的.shx、.dbf、.prj、.sbn、.sbx还有一个上海市.shp可能代表市级点位。这里有个经验如果只把.shp拷给别人对方看不到属性表只拷.dbf又看不到图形。.shp存坐标.shx是几何索引.dbf存属性.prj是坐标系描述.sbn/.sbx是 ArcGIS 的空间索引。.shp.xml是元数据有没有都不影响读取。ArcGIS 整理过的数据通常会有.sbn/.sbx这套资料里也带了说明生成时已经跑过空间索引。QGIS 加载时只需要选中.shp但如果你是做二次开发要注意OGR读取时如果同目录缺少.shx或.dbf会直接抛错。移动目录时我会保留全量文件而不是像网盘里那样只传.shp和.dbf。文件后缀作用丢失后果.shp要素几何坐标无法读取要素.shx几何位置索引部分软件加载慢或失败.dbf属性字段丢失所有非空间属性.prj坐标系描述坐标系未知叠加错位.sbn / .sbxArcGIS 空间索引不损坏数据可重建.shp.xml元数据信息不影响读取提示.sbn/.sbx不是必须的但当你用 ArcCatalog 移动或删除文件时不要只删这几个索引文件容易留下无效引用。2.2 区县、乡镇、村三套要素类的字段差异与用途分析这套数据的组织方式是把不同行政级别的点拆成独立要素类区县点.shp、乡镇点.shp、村点.shp外加上海市.shp。这种按层级拆分比单要素表更利于分析因为村点属性通常比区县点多比如可能包含行政村名、行政区划代码、人口或备注等而区县点层面可能只有名称和代码。实际处理时我会先看每个图层的字段结构用 GDAL 的工具能快速完成ogrinfo -al 区县点.shp | head -n 30说明-al表示列出所有要素head -n 30限制输出行数用来快速检查字段名和要素类型。不同图层的名称字段可能叫NAME、MC、XZQMC如果不统一后面关联就会出问题。我的习惯是第一遍编辑就统一成name和code两个字段并把字段类型设成字符串避免 join 的时候因为字符型与数值型不一致导致匹配失败。2.3 坐标参考系统.prj检查为什么必须先做一致性验证每个.prj文件都是纯文本直接用记事本打开能看到坐标系定义。上海这边的数据可能用到GCS_WGS_1984、CGCS2000也可能用到高斯投影或上海地方坐标。这里最容易踩的坑是区县点是 WGS84 经纬度乡镇点却是 CGCS2000 的投影坐标叠在一个工程里显示就会差几百米甚至更大。我拿到 zip 之后不做任何分析先跑一句ogrinfo -so 上海市.shp 上海市 | grep Crs\|EXTENT说明-so只输出摘要不展开每个属性grep过滤出坐标系和经纬度范围。如果四个图层输出的Crs不一致就需要用ogr2ogr -t_srs EPSG:4490统一到 CGCS2000 地理坐标系。以后再叠加路网或 POI 数据就能避免“点在海里”的经典问题。3. 用 QGIS / GDAL 把上海三级点位 SHP 变成可分析的数据集3.1 加载 SHP 并叠加底图投影坐标与经纬度的区分常见的错误是把所有图层直接拖进 QGIS然后叠加高德或 OSM 底图结果点位偏移。在线底图一般为 Web 墨卡托投影 EPSG:3857而原始 SHP 可能是经纬度坐标系必须先统一。我通常会在 QGIS 里把原始数据“导出 另存为”目标坐标系选择 EPSG:3857生成一个临时分析副本。这样底图、点位、边界都在同一投影下测量距离和生成缓冲区的结果也更直观。要注意的是如果你要计算面积EPSG:3857 在高纬度地区变形大上海处于中纬度勉强能用但严谨起见还是转成 UTM 或上海地方投影。坐标系EPSG 代码典型应用WGS84 经纬度4326原始 GPS 数据通用交换CGCS2000 经纬度4490国内测绘成果推荐使用Web 墨卡托3857切片底图、前端展示WGS84 UTM 51N32651上海所在带距离面积较准3.2 属性字段清洗中文编码问题与字段重命名用 QGIS 打开含中文属性的.dbf时经常出现乱码。原因是 shapefile 的 dbf 文件常用 GBK 编码而 QGIS 默认按 UTF-8 读取。解决办法是在“数据源管理器”里手动把编码改成GBK或GB2312。如果字段名本身是乱码就需要重命名。在 QGIS Python 控制台里可以这样处理layer QgsVectorLayer(区县点.shp, 区县, ogr) layer.dataProvider().renameAttributes({0: name, 1: code}) layer.updateFields() print([f.name() for f in layer.fields()])说明renameAttributes的参数是一个字典键是字段索引值是新字段名。执行完必须调用updateFields()才生效。我只建议在副本上做这一步因为原始字段名可能承载着来源信息直接覆盖后不容易反查数据责任方。3.3 点位数据空间关联把坐标点匹配到区域边界用 shp 转 txt 做数据交换在做数据交换时需要把村点坐标导出成文本文件“shp 转 txt”是最常见的形式之一。用 Python GDAL 能一次性输出属性、经度、纬度from osgeo import ogr ds ogr.Open(r村点.shp) layer ds.GetLayer() fields [f.GetName() for f in layer.schema] with open(村点_输出.txt, w, encodingutf-8) as f: f.write(,.join(fields) ,lng,lat\n) for feat in layer: geom feat.geometry() if geom is None: continue lng, lat geom.GetX(), geom.GetY() row [str(feat.GetField(fn)) for fn in fields] row.extend([str(lng), str(lat)]) f.write(,.join(row) \n)说明layer.schema拿到全部属性字段名geom.GetX()/GetY()取出要素坐标。输出结果是源文件坐标系下的坐标如果源文件是投影坐标应该先调用geom.TransformTo(target_srs)转成经纬度再写文本。这个脚本改一下就适用于乡镇点.shp和区县点.shp。3.4 点密度热力图与区域统计一个可复制的实现思路做分析时常需要统计每个乡镇范围内包含多少村点。QGIS 的“按位置连接”能完成但手动点菜单效率低。用处理框架写一小段脚本即可result processing.run(native:buffer, { INPUT: 村点.shp, DISTANCE: 500, OUTPUT: 内存:村点缓冲 }) count processing.run(native:countpointsinpolygon, { POLYGONS: result[OUTPUT], POINTS: 区县点.shp, OUTPUT: 内存:计数 })说明第一个工具为村点生成 500 米缓冲区第二个工具统计每个缓冲区内的区县点数量。DISTANCE的单位跟随源文件坐标系如果图层是经纬度这个数值会被当成度500 度显然不合理。所以必须先投影到以米为单位的坐标系再执行这个脚本。4. 进阶应用SHP 批量转换与压缩的实战套路4.1 shp 批量转 KML用 QGIS 批处理与 GDAL 命令用 QGIS“图层 另存为”可以转单个文件但如果需要把整个目录下的 shp 都转成 KML手工操作效率太低。我在 Windows CMD 下会写一个循环for %f in (*.shp) do ogr2ogr -f KML %~nf.kml %f -dsco NameFieldNAME说明%~nf提取当前文件的文件名前缀-f KML指定输出格式-dsco NameFieldNAME把属性字段 NAME 作为 KML 地标名称。如果你的字段名是中文建议先重命名字段再转 KML否则生成的图层名和要素名可能是空值。4.2 单个 SHP 批量压缩成 Shapefile 或 GeoPackage解决超大文件拖拽卡顿点数据超大时QGIS 初始化很慢拖动缩放手感也差。我一般分两步优化。第一步是简化几何适当降低精度。GDAL 自带简化参数例如ogr2ogr -simplify 0.0001 简化后.shp 原文件.shp说明0.0001是容差单位与源坐标系一致。如果源文件是经纬度这个数值大约是 10 米对点位数据影响不大但能显著减小文件体积。第二步是转成 GeoPackage。GeoPackage 读写比 shapefile 更高效属性字段类型支持更好。转换命令如下ogr2ogr -f GPKG 上海市点数据.gpkg 区县点.shp注意如果源 dbf 是 GBK 编码转换时建议加-lco ENCODINGGBK否则中文属性进入 GeoPackage 后可能变成乱码。4.3 从 Excel 经纬度生成 SHP一套可复用的脚本业务方经常会拿着 Excel 表里面有名称、经度、纬度要求转成 shp。环境里没有 geopandas 的时候可以用 GDAL/OGR 手动建点图层from osgeo import ogr, osr src_data [(113.921, 30.510), (121.473, 31.230)] driver ogr.GetDriverByName(ESRI Shapefile) ds driver.CreateDataSource(输出点.shp) layer ds.CreateLayer(输出点, srsosr.SpatialReference().ImportFromEPSG(4326), geom_typeogr.wkbPoint) layer.CreateField(ogr.FieldDefn(name, ogr.OFTString)) for i, (lng, lat) in enumerate(src_data): feat layer.CreateFeature(ogr.Feature(layer.GetLayerDefn())) feat.SetField(name, fpoi_{i}) wkt fPOINT ({lng} {lat}) feat.SetGeometry(ogr.CreateGeometryFromWkt(wkt)) layer.CreateFeature(feat) ds.Destroy()说明ImportFromEPSG(4326)把坐标系定义为 WGS84如果 Excel 里是经纬度必须使用这个。CreateFeature每循环一次生成一个新要素最后Destroy()释放句柄并落盘。如果输出文件已经存在CreateDataSource会报错所以脚本里最好先判断并清理同名文件。输出格式推荐命令适用场景KMLogr2ogr -f KML前端地图、Google EarthGeoPackageogr2ogr -f GPKGQGIS 分析、服务端发布CSV/TXTGDAL Python 脚本数据交换、机器学习导入5. 最后一步验证点位的边界与拓扑完整性避免“点不在面上”的坑5.1 只保留外边界线检查点是否落在行政区边界内不少人问“shp 有没有办法只保留外边界线”。如果你拿到的是包含很多碎面的区县边界先用 QGIS 的“融合”工具把碎面合成一个整体再对点位做空间过滤。命令行方式更直接ogr2ogr -clipsrc 上海市边界.shp 边界内村点.shp 村点.shp说明-clipsrc用上海市边界.shp的范围裁剪村点.shp输出到边界内村点.shp。执行前提是两个文件坐标系一致否则会出现裁剪结果偏移。如果裁剪后少了几百个点先去查原始 SHP 的.prj再做坐标转换。5.2 属性完整性核对将缺失值输出到 txt点位坐标正确不代表属性完整name或code字段为空会导致后续关联失败。我常用一个简短脚本把缺失字段的要素编号导出from osgeo import ogr ds ogr.Open(r乡镇点.shp) layer ds.GetLayer() check_fields [name, code] with open(缺失值.txt, w) as f: for feat in layer: for field in check_fields: value feat.GetField(field) if value is None or str(value).strip() : f.write(f{feat.GetFID()},{field}\n)说明feat.GetFID()返回要素编号回到 QGIS 中可以通过“按属性选择”快速定位到具体点。这个脚本对区县、乡镇、村级点位图层都能复用只需要替换文件路径和检查字段列表。距离和面积计算前先把所有点位图层统一重投影到同一个以米为单位的坐标系副本再开始分析。这是两分钟内能完成的动作但能省掉后续大量因单位不一致导致的返工。本文还有配套的精品资源点击获取 上海区县乡镇村SHP点位数据处理：坐标系校验与批量转换实战
简介上海市县区乡镇村各级点位Shapefile数据包专注于提供全市范围内区县、乡镇、村三级行政单位位置点信息适用于GIS空间分析、区域对比、专题制图及智慧城市相关研究。面向地理信息开发人员、数据分析师及城市规划研究者可快速用于点位展示、空间查询与批量匹配等场景。资源共24个文件以.shp几何数据、.dbf属性表、.shx索引、.prj坐标参考为核心另含.sbn/.sbx空间索引及.xml元数据文件覆盖三类图层数据完整、结构清晰包体仅460KB便于下载与处理。目前已有173人学习浏览可作为中小尺度区域分析的基础底图数据。通过该数据包用户既能获得各级点位的地理坐标又能结合属性字段开展分级统计、热点识别等深入分析有效提升区域研究工作与项目管理效率。1. 为什么拿到这套上海市县区乡镇村点位 SHP先别急着拖进 ArcMap做上海范围内的点位分析与区域分析最怕的不是没有数据而是坐标系混乱、字段乱码、边界对不上。这套点位包把市级、区县级、乡镇级、村级四个层级的点位置放在同一套 Shapefile 文件里短文件名直接表明层级属性表结构相对一致适合直接用于选址、可达性分析、商圈划分和 POI 质量校验。但 SHP 不是单个文件而是一组配套文件.shp只是几何.dbf才是属性少一个.prj就不知道经纬度还是投影坐标。下面按文件结构、字段清洗、实际分析、批量转换的顺序把整套流程拆开。2. Shapefile 文件家族解密从 .shp 到 .sbx每个文件都别乱删2.1 一套完整 SHP 数据的文件构成与加载时依赖顺序解压后你会看到区县点.shp、乡镇点.shp、村点.shp以及对应的.shx、.dbf、.prj、.sbn、.sbx还有一个上海市.shp可能代表市级点位。这里有个经验如果只把.shp拷给别人对方看不到属性表只拷.dbf又看不到图形。.shp存坐标.shx是几何索引.dbf存属性.prj是坐标系描述.sbn/.sbx是 ArcGIS 的空间索引。.shp.xml是元数据有没有都不影响读取。ArcGIS 整理过的数据通常会有.sbn/.sbx这套资料里也带了说明生成时已经跑过空间索引。QGIS 加载时只需要选中.shp但如果你是做二次开发要注意OGR读取时如果同目录缺少.shx或.dbf会直接抛错。移动目录时我会保留全量文件而不是像网盘里那样只传.shp和.dbf。文件后缀作用丢失后果.shp要素几何坐标无法读取要素.shx几何位置索引部分软件加载慢或失败.dbf属性字段丢失所有非空间属性.prj坐标系描述坐标系未知叠加错位.sbn / .sbxArcGIS 空间索引不损坏数据可重建.shp.xml元数据信息不影响读取提示.sbn/.sbx不是必须的但当你用 ArcCatalog 移动或删除文件时不要只删这几个索引文件容易留下无效引用。2.2 区县、乡镇、村三套要素类的字段差异与用途分析这套数据的组织方式是把不同行政级别的点拆成独立要素类区县点.shp、乡镇点.shp、村点.shp外加上海市.shp。这种按层级拆分比单要素表更利于分析因为村点属性通常比区县点多比如可能包含行政村名、行政区划代码、人口或备注等而区县点层面可能只有名称和代码。实际处理时我会先看每个图层的字段结构用 GDAL 的工具能快速完成ogrinfo -al 区县点.shp | head -n 30说明-al表示列出所有要素head -n 30限制输出行数用来快速检查字段名和要素类型。不同图层的名称字段可能叫NAME、MC、XZQMC如果不统一后面关联就会出问题。我的习惯是第一遍编辑就统一成name和code两个字段并把字段类型设成字符串避免 join 的时候因为字符型与数值型不一致导致匹配失败。2.3 坐标参考系统.prj检查为什么必须先做一致性验证每个.prj文件都是纯文本直接用记事本打开能看到坐标系定义。上海这边的数据可能用到GCS_WGS_1984、CGCS2000也可能用到高斯投影或上海地方坐标。这里最容易踩的坑是区县点是 WGS84 经纬度乡镇点却是 CGCS2000 的投影坐标叠在一个工程里显示就会差几百米甚至更大。我拿到 zip 之后不做任何分析先跑一句ogrinfo -so 上海市.shp 上海市 | grep Crs\|EXTENT说明-so只输出摘要不展开每个属性grep过滤出坐标系和经纬度范围。如果四个图层输出的Crs不一致就需要用ogr2ogr -t_srs EPSG:4490统一到 CGCS2000 地理坐标系。以后再叠加路网或 POI 数据就能避免“点在海里”的经典问题。3. 用 QGIS / GDAL 把上海三级点位 SHP 变成可分析的数据集3.1 加载 SHP 并叠加底图投影坐标与经纬度的区分常见的错误是把所有图层直接拖进 QGIS然后叠加高德或 OSM 底图结果点位偏移。在线底图一般为 Web 墨卡托投影 EPSG:3857而原始 SHP 可能是经纬度坐标系必须先统一。我通常会在 QGIS 里把原始数据“导出 另存为”目标坐标系选择 EPSG:3857生成一个临时分析副本。这样底图、点位、边界都在同一投影下测量距离和生成缓冲区的结果也更直观。要注意的是如果你要计算面积EPSG:3857 在高纬度地区变形大上海处于中纬度勉强能用但严谨起见还是转成 UTM 或上海地方投影。坐标系EPSG 代码典型应用WGS84 经纬度4326原始 GPS 数据通用交换CGCS2000 经纬度4490国内测绘成果推荐使用Web 墨卡托3857切片底图、前端展示WGS84 UTM 51N32651上海所在带距离面积较准3.2 属性字段清洗中文编码问题与字段重命名用 QGIS 打开含中文属性的.dbf时经常出现乱码。原因是 shapefile 的 dbf 文件常用 GBK 编码而 QGIS 默认按 UTF-8 读取。解决办法是在“数据源管理器”里手动把编码改成GBK或GB2312。如果字段名本身是乱码就需要重命名。在 QGIS Python 控制台里可以这样处理layer QgsVectorLayer(区县点.shp, 区县, ogr) layer.dataProvider().renameAttributes({0: name, 1: code}) layer.updateFields() print([f.name() for f in layer.fields()])说明renameAttributes的参数是一个字典键是字段索引值是新字段名。执行完必须调用updateFields()才生效。我只建议在副本上做这一步因为原始字段名可能承载着来源信息直接覆盖后不容易反查数据责任方。3.3 点位数据空间关联把坐标点匹配到区域边界用 shp 转 txt 做数据交换在做数据交换时需要把村点坐标导出成文本文件“shp 转 txt”是最常见的形式之一。用 Python GDAL 能一次性输出属性、经度、纬度from osgeo import ogr ds ogr.Open(r村点.shp) layer ds.GetLayer() fields [f.GetName() for f in layer.schema] with open(村点_输出.txt, w, encodingutf-8) as f: f.write(,.join(fields) ,lng,lat\n) for feat in layer: geom feat.geometry() if geom is None: continue lng, lat geom.GetX(), geom.GetY() row [str(feat.GetField(fn)) for fn in fields] row.extend([str(lng), str(lat)]) f.write(,.join(row) \n)说明layer.schema拿到全部属性字段名geom.GetX()/GetY()取出要素坐标。输出结果是源文件坐标系下的坐标如果源文件是投影坐标应该先调用geom.TransformTo(target_srs)转成经纬度再写文本。这个脚本改一下就适用于乡镇点.shp和区县点.shp。3.4 点密度热力图与区域统计一个可复制的实现思路做分析时常需要统计每个乡镇范围内包含多少村点。QGIS 的“按位置连接”能完成但手动点菜单效率低。用处理框架写一小段脚本即可result processing.run(native:buffer, { INPUT: 村点.shp, DISTANCE: 500, OUTPUT: 内存:村点缓冲 }) count processing.run(native:countpointsinpolygon, { POLYGONS: result[OUTPUT], POINTS: 区县点.shp, OUTPUT: 内存:计数 })说明第一个工具为村点生成 500 米缓冲区第二个工具统计每个缓冲区内的区县点数量。DISTANCE的单位跟随源文件坐标系如果图层是经纬度这个数值会被当成度500 度显然不合理。所以必须先投影到以米为单位的坐标系再执行这个脚本。4. 进阶应用SHP 批量转换与压缩的实战套路4.1 shp 批量转 KML用 QGIS 批处理与 GDAL 命令用 QGIS“图层 另存为”可以转单个文件但如果需要把整个目录下的 shp 都转成 KML手工操作效率太低。我在 Windows CMD 下会写一个循环for %f in (*.shp) do ogr2ogr -f KML %~nf.kml %f -dsco NameFieldNAME说明%~nf提取当前文件的文件名前缀-f KML指定输出格式-dsco NameFieldNAME把属性字段 NAME 作为 KML 地标名称。如果你的字段名是中文建议先重命名字段再转 KML否则生成的图层名和要素名可能是空值。4.2 单个 SHP 批量压缩成 Shapefile 或 GeoPackage解决超大文件拖拽卡顿点数据超大时QGIS 初始化很慢拖动缩放手感也差。我一般分两步优化。第一步是简化几何适当降低精度。GDAL 自带简化参数例如ogr2ogr -simplify 0.0001 简化后.shp 原文件.shp说明0.0001是容差单位与源坐标系一致。如果源文件是经纬度这个数值大约是 10 米对点位数据影响不大但能显著减小文件体积。第二步是转成 GeoPackage。GeoPackage 读写比 shapefile 更高效属性字段类型支持更好。转换命令如下ogr2ogr -f GPKG 上海市点数据.gpkg 区县点.shp注意如果源 dbf 是 GBK 编码转换时建议加-lco ENCODINGGBK否则中文属性进入 GeoPackage 后可能变成乱码。4.3 从 Excel 经纬度生成 SHP一套可复用的脚本业务方经常会拿着 Excel 表里面有名称、经度、纬度要求转成 shp。环境里没有 geopandas 的时候可以用 GDAL/OGR 手动建点图层from osgeo import ogr, osr src_data [(113.921, 30.510), (121.473, 31.230)] driver ogr.GetDriverByName(ESRI Shapefile) ds driver.CreateDataSource(输出点.shp) layer ds.CreateLayer(输出点, srsosr.SpatialReference().ImportFromEPSG(4326), geom_typeogr.wkbPoint) layer.CreateField(ogr.FieldDefn(name, ogr.OFTString)) for i, (lng, lat) in enumerate(src_data): feat layer.CreateFeature(ogr.Feature(layer.GetLayerDefn())) feat.SetField(name, fpoi_{i}) wkt fPOINT ({lng} {lat}) feat.SetGeometry(ogr.CreateGeometryFromWkt(wkt)) layer.CreateFeature(feat) ds.Destroy()说明ImportFromEPSG(4326)把坐标系定义为 WGS84如果 Excel 里是经纬度必须使用这个。CreateFeature每循环一次生成一个新要素最后Destroy()释放句柄并落盘。如果输出文件已经存在CreateDataSource会报错所以脚本里最好先判断并清理同名文件。输出格式推荐命令适用场景KMLogr2ogr -f KML前端地图、Google EarthGeoPackageogr2ogr -f GPKGQGIS 分析、服务端发布CSV/TXTGDAL Python 脚本数据交换、机器学习导入5. 最后一步验证点位的边界与拓扑完整性避免“点不在面上”的坑5.1 只保留外边界线检查点是否落在行政区边界内不少人问“shp 有没有办法只保留外边界线”。如果你拿到的是包含很多碎面的区县边界先用 QGIS 的“融合”工具把碎面合成一个整体再对点位做空间过滤。命令行方式更直接ogr2ogr -clipsrc 上海市边界.shp 边界内村点.shp 村点.shp说明-clipsrc用上海市边界.shp的范围裁剪村点.shp输出到边界内村点.shp。执行前提是两个文件坐标系一致否则会出现裁剪结果偏移。如果裁剪后少了几百个点先去查原始 SHP 的.prj再做坐标转换。5.2 属性完整性核对将缺失值输出到 txt点位坐标正确不代表属性完整name或code字段为空会导致后续关联失败。我常用一个简短脚本把缺失字段的要素编号导出from osgeo import ogr ds ogr.Open(r乡镇点.shp) layer ds.GetLayer() check_fields [name, code] with open(缺失值.txt, w) as f: for feat in layer: for field in check_fields: value feat.GetField(field) if value is None or str(value).strip() : f.write(f{feat.GetFID()},{field}\n)说明feat.GetFID()返回要素编号回到 QGIS 中可以通过“按属性选择”快速定位到具体点。这个脚本对区县、乡镇、村级点位图层都能复用只需要替换文件路径和检查字段列表。距离和面积计算前先把所有点位图层统一重投影到同一个以米为单位的坐标系副本再开始分析。这是两分钟内能完成的动作但能省掉后续大量因单位不一致导致的返工。本文还有配套的精品资源点击获取