简介面向GIS初学者、地理信息专业学生及区域规划研究者的江苏省矢量地理数据包以Shapefile标准格式整理可用于地图制作、空间分析与专题制图练习。压缩包共83个文件约974KB涵盖shp几何数据、dbf属性表、shx索引及prj投影等配套文件并附sbn、sbx空间索引与xml元数据保证在ArcGIS、QGIS等软件中直接读取。内容覆盖河流、湖泊、县域、地级市、高速公路、国道、铁路、长江、丘陵、湿地及经纬网等要素可支撑洪水风险叠加评估、地区发展差异分析、交通网络优化等典型场景。已有2386人学习下载说明其在教学与项目实践中具备一定参考价值。数据分层清晰、要素齐全便于读者快速搭建江苏地理底图开展空间查询、统计与可视化表达为区域研究与决策提供基础支撑。1. 江苏省河流湖泊县域高速等矢量数据的落地用法从拿到 shp 到能上图能算量手里有一份江苏省的河流、湖泊、县域、高速 shp 数据第一反应往往不是“太好了”而是“这堆文件到底怎么用”。我见过太多人拿到 shp 之后卡在第一步双击打不开拖进软件报坐标系错误图层叠上去位置偏到海里。这类基础地理矢量数据真正的价值不在于“有”而在于能不能稳定地做空间叠加、面积量算、缓冲区分析和出图。它适合做区域规划、水利分析、交通可达性评估、专题制图的人也适合需要给业务系统准备底图数据的开发者。这一篇就按我实际处理这类数据的顺序把坐标系、字段、拓扑、叠加和验证一条线讲透让你拿到同类数据能直接复现。2. 先搞清 shp 到底是什么五个文件缺一不可2.1 shp 不是单个文件是一组同名文件很多人以为 shp 就是一个文件复制的时候只拷了.shp结果到另一台机器上属性全没了。Shapefile 是一套同名不同后缀的文件集合少一个就残废。常见组成如下后缀作用缺失后果.shp存几何坐标直接打不开.shx几何索引部分软件无法定位要素.dbf属性表能显示图形但没有字段.prj坐标系定义坐标数值在但不知道是什么坐标系.cpg字符编码中文属性乱码我一般拿到数据先做一件事按文件名分组确认每个图层都是完整的一套。河流、湖泊、县域、高速通常是四个独立图层每个图层各自一套文件。如果.prj缺失不要急着猜坐标系先看坐标数值范围经纬度一般是小数点后带多位、范围在 70 到 140 之间投影坐标往往是六位或七位的大数。这个判断能帮你避开后面一大半的偏移问题。2.2 用命令行先做一次体检图形界面打开之前我习惯先用命令行把基本情况摸清楚速度快且不依赖软件。GDAL 的ogrinfo是最稳的入口# 查看单个图层的概要信息包括几何类型、要素数、坐标系 ogrinfo -so -al jiangsu_county.shp # 查看属性表前几条记录确认字段和编码 ogrinfo -al -geomNO jiangsu_county.shp | head -60-so表示只输出摘要不读几何-al表示所有图层-geomNO跳过几何只读属性。输出里重点看三样Geometry 是 Polygon 还是 LineStringFeature Count 是否和预期一致Layer SRS 是否带 EPSG 编号。如果 SRS 显示(unknown)说明.prj缺失或损坏后面必须手动指定。提示ogrinfo读中文属性乱码时加环境变量SHAPE_ENCODINGGBK再执行能解决大部分老数据集的编码问题。2.3 几何类型决定你能做什么分析河流和高速是线要素湖泊和县域是面要素这个区别不是学术问题直接决定分析手段。线要素能做长度量算、缓冲区、网络连通性面要素能做面积、叠加、包含判断。把线当面试图算面积或者把面当线做路径分析结果一定是错的。我在处理高速数据时踩过一次坑某份高速图层里混进了互通立交的面要素几何类型不纯导致后续按长度统计时多出一批零长度记录。所以体检时一定要确认几何类型是否单一混合类型要么拆分要么过滤。3. 坐标系这一步错了后面全白干3.1 地理坐标系和投影坐标系别混用江苏地处东经 116 到 121 度、北纬 30 到 35 度之间常见的地理坐标系是 WGS84EPSG:4326或 CGCS2000EPSG:4490单位是度。做面积和长度量算时必须转到投影坐标系单位才是米。江苏常用的投影有 UTM 50NEPSG:32650或高斯克吕格 3 度带。我一般统一转到 CGCS2000 的 3 度带投影EPSG 用 4547 到 4554 这一组里对应中央经线的那个具体看数据覆盖范围。判断当前坐标系最直接的办法是看.prj内容或者用ogrinfo输出的 SRS 字段。如果显示经纬度但你要算面积先投影如果显示投影坐标但你要和在线底图叠加先转回经纬度。混用的典型症状是图层叠上去整体偏移几百米到几公里或者面积算出来是零点零几明显不对。3.2 用 ogr2ogr 批量统一坐标系四个图层坐标系不一致是常态手工一个个转容易漏。我一般写一条循环命令统一处理# 把所有 shp 统一转到 CGCS2000 3度带投影以 EPSG:4547 为例 for f in *.shp; do base${f%.shp} ogr2ogr -f ESRI Shapefile -t_srs EPSG:4547 \ reproj/${base}.shp $f \ -lco ENCODINGUTF-8 done-t_srs是目标坐标系-lco ENCODINGUTF-8让输出的 dbf 用 UTF-8 编码避免中文乱码。源坐标系如果.prj存在会自动读取不存在就得加-s_srs手动指定。转完之后不要只看命令有没有报错一定要抽一个已知地物核对坐标数值量级投影坐标应该是几十万到几百万的米制数值不再是 118、32 这种度数。注意-s_srs指定错了不会报错只会默默把数据转到错误位置。不确定源坐标系时宁可先不转用已知控制点反推也不要凭感觉填。3.3 投影选择对面积量算的影响同一个湖泊用不同投影算出来的面积能差百分之一到百分之几。对一般专题图无所谓对水利统计、用地平衡这类要报数的场景就不能含糊。我的做法是全省尺度用等积投影或对应中央经线的高斯投影局部小区域用 UTM 对应带。判断标准很简单量算结果和官方公布的控制面积对一下偏差在可接受范围内就固定用这套参数别每次换。4. 字段、编码和拓扑属性对了图才可信4.1 中文属性乱码的三种修法dbf 的编码问题几乎是这类数据的标配。表现是属性表里中文变成问号或方块。修法分三层第一层看有没有.cpg文件里面写着UTF-8或GBK第二层用ogrinfo加SHAPE_ENCODING环境变量试读第三层实在不行用 Python 直接读 dbf 探测编码再重写。我一般优先用 GDAL 重写一遍带 UTF-8 的输出比在软件里手动改编码稳。from osgeo import ogr, gdal gdal.SetConfigOption(SHAPE_ENCODING, GBK) # 源数据按 GBK 读 src ogr.Open(jiangsu_river.shp) driver ogr.GetDriverByName(ESRI Shapefile) dst driver.CreateDataSource(river_utf8.shp) # 复制图层时显式指定 UTF-8 输出编码 dst_layer dst.CreateLayer(river_utf8, src.GetLayer(0).GetSpatialRef(), ogr.wkbLineString, [ENCODINGUTF-8])这段的关键是SHAPE_ENCODING控制读入解码ENCODINGUTF-8控制写出编码两头对齐才不会乱。参数上源编码不确定时可以先用chardet之类工具对 dbf 做一次探测但别完全信自动探测中文短字段容易误判最终还是要人工抽几条看。4.2 县域面数据的拓扑检查县域图层最常见的两个问题相邻县之间有缝隙或重叠边界不闭合。缝隙会导致按县汇总面积时总和小于全省面积重叠会导致重复计算。检查方法是用拓扑规则跑一遍或者简单点把所有县合并成一个面和全省边界对比面积差。差得离谱就说明拓扑有问题。# 合并所有县为一个面用于和全省边界做面积对比 ogr2ogr -f ESRI Shapefile county_union.shp jiangsu_county.shp \ -dialect SQLite -sql SELECT ST_Union(geometry) AS geometry FROM jiangsu_countyST_Union会把所有面合并如果存在缝隙合并结果里会保留内部空洞。参数上注意 SQLite 方言下几何字段名通常就是geometry不同数据源可能叫SHAPE报错时先ogrinfo看字段名。合并后算面积和官方全省面积对一下偏差超过千分之几就值得回头查边界。4.3 河流线要素的方向和连通性河流数据如果要做网络分析线的方向起点到终点和端点是否相接很关键。很多公开数据的河流线是分段采集的相邻段端点差几米没接上网络分析就会断。检查办法是统计悬挂端点数量或者直接建网络看连通分量。我一般先用缓冲区把端点吸附一下再建网容差取几米到十几米具体看数据精度。高速数据同理出入口和主线是否拓扑相连直接决定可达性分析能不能跑通。5. 叠加分析与出图把四类数据用起来5.1 用县域面裁剪河流和高速最常见的需求是“每个县有多少公里河流、多少公里高速”。做法是用县域面去裁剪线要素再按县分组统计长度。命令行用ogr2ogr配合 SQL 就能完成不必打开桌面软件。# 按县裁剪高速线输出带县名的分段线 ogr2ogr -f ESRI Shapefile highway_by_county.shp jiangsu_county.shp \ -dialect SQLite -sql SELECT c.name AS county, ST_Intersection(h.geometry, c.geometry) AS geometry FROM jiangsu_county c, jiangsu_highway.shp.jiangsu_highway h WHERE ST_Intersects(c.geometry, h.geometry) 这里ST_Intersection求交得到每县境内的线段ST_Intersects先做粗筛减少计算量。参数上要注意两个图层必须同一坐标系否则求交结果为空或错位。输出后按county字段分组求长度就得到每县高速里程。河流同理只是把高速图层换掉。5.2 湖泊与县域的包含关系判断湖泊跨县是常态需要知道每个湖涉及哪些县、各占多少面积。用ST_Intersection求交后按湖和县分组算面积即可。判断包含关系时注意用ST_Within还是ST_Intersects前者要求完全在内部后者只要沾边就算。跨县湖泊用ST_Intersects才不会漏。算完面积记得核对各分块之和是否等于湖泊总面积差太多说明几何有问题。5.3 出图前的样式和标注处理出图阶段最容易翻车的是标注压盖和线宽随缩放变化。我的习惯是河流按等级分线宽高速用双线或粗线加颜色区分县域只描边不填充或浅填充湖泊填充蓝色。标注用县域名称设置最小显示比例避免小比例尺下糊成一片。导出前把图层顺序理清面在下、线在上、标注最上。这些是制图常识但数据本身字段不全时得先补一个名称字段再标注否则只能显示 FID。6. 避坑与排查这几处翻车最多6.1 现象图层叠加整体偏移几百米原因两个图层坐标系不同一个是地理坐标一个是投影坐标或者投影带选错。解决统一转到同一投影坐标系用已知地物如某个湖心核对坐标数值量级确认不再是度数。6.2 现象面积算出来只有零点零几原因在地理坐标系下直接算面积单位是“平方度”没有物理意义。解决先投影到米制坐标系再算投影参数按区域选择算完和已知面积对一下。6.3 现象属性表中文全是问号原因dbf 编码和读取端不一致常见是源数据 GBK 被当 UTF-8 读。解决用SHAPE_ENCODING指定源编码输出时显式写 UTF-8重写一遍数据。6.4 现象按县汇总面积总和小于全省原因县域之间存在缝隙或者边界不闭合。解决做一次合并对比面积定位缝隙位置必要时用容差做一次吸附或手工修边。6.5 现象网络分析跑不通路径断裂原因河流或高速线端点未相接存在悬挂点。解决建网前做端点吸附容差按数据精度取几米吸附后重新检查连通分量。7. 一套可复用的验证脚本与参数习惯处理完数据我最后一定会跑一遍验证不靠肉眼。核心是三件事坐标系对不对、面积长度量级对不对、拓扑有没有明显断裂。下面这段 Python 把前两件做了第三件用前面的合并对比法补上。from osgeo import ogr import os def check_layer(path): ds ogr.Open(path) lyr ds.GetLayer(0) srs lyr.GetSpatialRef() print(f图层: {os.path.basename(path)}) print(f 要素数: {lyr.GetFeatureCount()}) print(f 几何类型: {ogr.GeometryTypeToName(lyr.GetGeomType())}) if srs: print(f 坐标系: {srs.GetAuthorityName(None)}:{srs.GetAuthorityCode(None)}) else: print( 坐标系: 缺失需手动指定) # 抽第一个要素看坐标量级判断是经纬度还是投影 feat lyr.GetNextFeature() if feat: geom feat.GetGeometryRef() x, y geom.GetX(), geom.GetY() scale 经纬度 if abs(x) 180 and abs(y) 90 else 投影坐标 print(f 首点坐标量级: {scale} ({x:.2f}, {y:.2f})) for f in [jiangsu_river.shp, jiangsu_lake.shp, jiangsu_county.shp, jiangsu_highway.shp]: check_layer(f)这段脚本的价值在于把“感觉没问题”变成“有输出可核对”。参数上GetAuthorityCode返回 None 就说明.prj里没有 EPSG 编号需要人工确认首点坐标量级判断只是粗筛最终还是要看 SRS 字段。我一般把这套检查固化成习惯每批数据进来先跑一遍比事后返工省时间。再补一个参数习惯所有中间成果单独放一个目录文件名带坐标系后缀比如_4547避免不同版本混在一起。投影参数、编码、容差这些值记在一个文本里随数据走换人接手也不用猜。这套做法不复杂但能挡住大部分“数据看着对、一算就错”的问题。希望帮到你。本文还有配套的精品资源点击获取 
长河君 · 2026/11/4 20:09:34 · 阅读 1.2k 江苏河流湖泊县域高速shp数据实战：坐标系、拓扑与叠加分析全解析
简介面向GIS初学者、地理信息专业学生及区域规划研究者的江苏省矢量地理数据包以Shapefile标准格式整理可用于地图制作、空间分析与专题制图练习。压缩包共83个文件约974KB涵盖shp几何数据、dbf属性表、shx索引及prj投影等配套文件并附sbn、sbx空间索引与xml元数据保证在ArcGIS、QGIS等软件中直接读取。内容覆盖河流、湖泊、县域、地级市、高速公路、国道、铁路、长江、丘陵、湿地及经纬网等要素可支撑洪水风险叠加评估、地区发展差异分析、交通网络优化等典型场景。已有2386人学习下载说明其在教学与项目实践中具备一定参考价值。数据分层清晰、要素齐全便于读者快速搭建江苏地理底图开展空间查询、统计与可视化表达为区域研究与决策提供基础支撑。1. 江苏省河流湖泊县域高速等矢量数据的落地用法从拿到 shp 到能上图能算量手里有一份江苏省的河流、湖泊、县域、高速 shp 数据第一反应往往不是“太好了”而是“这堆文件到底怎么用”。我见过太多人拿到 shp 之后卡在第一步双击打不开拖进软件报坐标系错误图层叠上去位置偏到海里。这类基础地理矢量数据真正的价值不在于“有”而在于能不能稳定地做空间叠加、面积量算、缓冲区分析和出图。它适合做区域规划、水利分析、交通可达性评估、专题制图的人也适合需要给业务系统准备底图数据的开发者。这一篇就按我实际处理这类数据的顺序把坐标系、字段、拓扑、叠加和验证一条线讲透让你拿到同类数据能直接复现。2. 先搞清 shp 到底是什么五个文件缺一不可2.1 shp 不是单个文件是一组同名文件很多人以为 shp 就是一个文件复制的时候只拷了.shp结果到另一台机器上属性全没了。Shapefile 是一套同名不同后缀的文件集合少一个就残废。常见组成如下后缀作用缺失后果.shp存几何坐标直接打不开.shx几何索引部分软件无法定位要素.dbf属性表能显示图形但没有字段.prj坐标系定义坐标数值在但不知道是什么坐标系.cpg字符编码中文属性乱码我一般拿到数据先做一件事按文件名分组确认每个图层都是完整的一套。河流、湖泊、县域、高速通常是四个独立图层每个图层各自一套文件。如果.prj缺失不要急着猜坐标系先看坐标数值范围经纬度一般是小数点后带多位、范围在 70 到 140 之间投影坐标往往是六位或七位的大数。这个判断能帮你避开后面一大半的偏移问题。2.2 用命令行先做一次体检图形界面打开之前我习惯先用命令行把基本情况摸清楚速度快且不依赖软件。GDAL 的ogrinfo是最稳的入口# 查看单个图层的概要信息包括几何类型、要素数、坐标系 ogrinfo -so -al jiangsu_county.shp # 查看属性表前几条记录确认字段和编码 ogrinfo -al -geomNO jiangsu_county.shp | head -60-so表示只输出摘要不读几何-al表示所有图层-geomNO跳过几何只读属性。输出里重点看三样Geometry 是 Polygon 还是 LineStringFeature Count 是否和预期一致Layer SRS 是否带 EPSG 编号。如果 SRS 显示(unknown)说明.prj缺失或损坏后面必须手动指定。提示ogrinfo读中文属性乱码时加环境变量SHAPE_ENCODINGGBK再执行能解决大部分老数据集的编码问题。2.3 几何类型决定你能做什么分析河流和高速是线要素湖泊和县域是面要素这个区别不是学术问题直接决定分析手段。线要素能做长度量算、缓冲区、网络连通性面要素能做面积、叠加、包含判断。把线当面试图算面积或者把面当线做路径分析结果一定是错的。我在处理高速数据时踩过一次坑某份高速图层里混进了互通立交的面要素几何类型不纯导致后续按长度统计时多出一批零长度记录。所以体检时一定要确认几何类型是否单一混合类型要么拆分要么过滤。3. 坐标系这一步错了后面全白干3.1 地理坐标系和投影坐标系别混用江苏地处东经 116 到 121 度、北纬 30 到 35 度之间常见的地理坐标系是 WGS84EPSG:4326或 CGCS2000EPSG:4490单位是度。做面积和长度量算时必须转到投影坐标系单位才是米。江苏常用的投影有 UTM 50NEPSG:32650或高斯克吕格 3 度带。我一般统一转到 CGCS2000 的 3 度带投影EPSG 用 4547 到 4554 这一组里对应中央经线的那个具体看数据覆盖范围。判断当前坐标系最直接的办法是看.prj内容或者用ogrinfo输出的 SRS 字段。如果显示经纬度但你要算面积先投影如果显示投影坐标但你要和在线底图叠加先转回经纬度。混用的典型症状是图层叠上去整体偏移几百米到几公里或者面积算出来是零点零几明显不对。3.2 用 ogr2ogr 批量统一坐标系四个图层坐标系不一致是常态手工一个个转容易漏。我一般写一条循环命令统一处理# 把所有 shp 统一转到 CGCS2000 3度带投影以 EPSG:4547 为例 for f in *.shp; do base${f%.shp} ogr2ogr -f ESRI Shapefile -t_srs EPSG:4547 \ reproj/${base}.shp $f \ -lco ENCODINGUTF-8 done-t_srs是目标坐标系-lco ENCODINGUTF-8让输出的 dbf 用 UTF-8 编码避免中文乱码。源坐标系如果.prj存在会自动读取不存在就得加-s_srs手动指定。转完之后不要只看命令有没有报错一定要抽一个已知地物核对坐标数值量级投影坐标应该是几十万到几百万的米制数值不再是 118、32 这种度数。注意-s_srs指定错了不会报错只会默默把数据转到错误位置。不确定源坐标系时宁可先不转用已知控制点反推也不要凭感觉填。3.3 投影选择对面积量算的影响同一个湖泊用不同投影算出来的面积能差百分之一到百分之几。对一般专题图无所谓对水利统计、用地平衡这类要报数的场景就不能含糊。我的做法是全省尺度用等积投影或对应中央经线的高斯投影局部小区域用 UTM 对应带。判断标准很简单量算结果和官方公布的控制面积对一下偏差在可接受范围内就固定用这套参数别每次换。4. 字段、编码和拓扑属性对了图才可信4.1 中文属性乱码的三种修法dbf 的编码问题几乎是这类数据的标配。表现是属性表里中文变成问号或方块。修法分三层第一层看有没有.cpg文件里面写着UTF-8或GBK第二层用ogrinfo加SHAPE_ENCODING环境变量试读第三层实在不行用 Python 直接读 dbf 探测编码再重写。我一般优先用 GDAL 重写一遍带 UTF-8 的输出比在软件里手动改编码稳。from osgeo import ogr, gdal gdal.SetConfigOption(SHAPE_ENCODING, GBK) # 源数据按 GBK 读 src ogr.Open(jiangsu_river.shp) driver ogr.GetDriverByName(ESRI Shapefile) dst driver.CreateDataSource(river_utf8.shp) # 复制图层时显式指定 UTF-8 输出编码 dst_layer dst.CreateLayer(river_utf8, src.GetLayer(0).GetSpatialRef(), ogr.wkbLineString, [ENCODINGUTF-8])这段的关键是SHAPE_ENCODING控制读入解码ENCODINGUTF-8控制写出编码两头对齐才不会乱。参数上源编码不确定时可以先用chardet之类工具对 dbf 做一次探测但别完全信自动探测中文短字段容易误判最终还是要人工抽几条看。4.2 县域面数据的拓扑检查县域图层最常见的两个问题相邻县之间有缝隙或重叠边界不闭合。缝隙会导致按县汇总面积时总和小于全省面积重叠会导致重复计算。检查方法是用拓扑规则跑一遍或者简单点把所有县合并成一个面和全省边界对比面积差。差得离谱就说明拓扑有问题。# 合并所有县为一个面用于和全省边界做面积对比 ogr2ogr -f ESRI Shapefile county_union.shp jiangsu_county.shp \ -dialect SQLite -sql SELECT ST_Union(geometry) AS geometry FROM jiangsu_countyST_Union会把所有面合并如果存在缝隙合并结果里会保留内部空洞。参数上注意 SQLite 方言下几何字段名通常就是geometry不同数据源可能叫SHAPE报错时先ogrinfo看字段名。合并后算面积和官方全省面积对一下偏差超过千分之几就值得回头查边界。4.3 河流线要素的方向和连通性河流数据如果要做网络分析线的方向起点到终点和端点是否相接很关键。很多公开数据的河流线是分段采集的相邻段端点差几米没接上网络分析就会断。检查办法是统计悬挂端点数量或者直接建网络看连通分量。我一般先用缓冲区把端点吸附一下再建网容差取几米到十几米具体看数据精度。高速数据同理出入口和主线是否拓扑相连直接决定可达性分析能不能跑通。5. 叠加分析与出图把四类数据用起来5.1 用县域面裁剪河流和高速最常见的需求是“每个县有多少公里河流、多少公里高速”。做法是用县域面去裁剪线要素再按县分组统计长度。命令行用ogr2ogr配合 SQL 就能完成不必打开桌面软件。# 按县裁剪高速线输出带县名的分段线 ogr2ogr -f ESRI Shapefile highway_by_county.shp jiangsu_county.shp \ -dialect SQLite -sql SELECT c.name AS county, ST_Intersection(h.geometry, c.geometry) AS geometry FROM jiangsu_county c, jiangsu_highway.shp.jiangsu_highway h WHERE ST_Intersects(c.geometry, h.geometry) 这里ST_Intersection求交得到每县境内的线段ST_Intersects先做粗筛减少计算量。参数上要注意两个图层必须同一坐标系否则求交结果为空或错位。输出后按county字段分组求长度就得到每县高速里程。河流同理只是把高速图层换掉。5.2 湖泊与县域的包含关系判断湖泊跨县是常态需要知道每个湖涉及哪些县、各占多少面积。用ST_Intersection求交后按湖和县分组算面积即可。判断包含关系时注意用ST_Within还是ST_Intersects前者要求完全在内部后者只要沾边就算。跨县湖泊用ST_Intersects才不会漏。算完面积记得核对各分块之和是否等于湖泊总面积差太多说明几何有问题。5.3 出图前的样式和标注处理出图阶段最容易翻车的是标注压盖和线宽随缩放变化。我的习惯是河流按等级分线宽高速用双线或粗线加颜色区分县域只描边不填充或浅填充湖泊填充蓝色。标注用县域名称设置最小显示比例避免小比例尺下糊成一片。导出前把图层顺序理清面在下、线在上、标注最上。这些是制图常识但数据本身字段不全时得先补一个名称字段再标注否则只能显示 FID。6. 避坑与排查这几处翻车最多6.1 现象图层叠加整体偏移几百米原因两个图层坐标系不同一个是地理坐标一个是投影坐标或者投影带选错。解决统一转到同一投影坐标系用已知地物如某个湖心核对坐标数值量级确认不再是度数。6.2 现象面积算出来只有零点零几原因在地理坐标系下直接算面积单位是“平方度”没有物理意义。解决先投影到米制坐标系再算投影参数按区域选择算完和已知面积对一下。6.3 现象属性表中文全是问号原因dbf 编码和读取端不一致常见是源数据 GBK 被当 UTF-8 读。解决用SHAPE_ENCODING指定源编码输出时显式写 UTF-8重写一遍数据。6.4 现象按县汇总面积总和小于全省原因县域之间存在缝隙或者边界不闭合。解决做一次合并对比面积定位缝隙位置必要时用容差做一次吸附或手工修边。6.5 现象网络分析跑不通路径断裂原因河流或高速线端点未相接存在悬挂点。解决建网前做端点吸附容差按数据精度取几米吸附后重新检查连通分量。7. 一套可复用的验证脚本与参数习惯处理完数据我最后一定会跑一遍验证不靠肉眼。核心是三件事坐标系对不对、面积长度量级对不对、拓扑有没有明显断裂。下面这段 Python 把前两件做了第三件用前面的合并对比法补上。from osgeo import ogr import os def check_layer(path): ds ogr.Open(path) lyr ds.GetLayer(0) srs lyr.GetSpatialRef() print(f图层: {os.path.basename(path)}) print(f 要素数: {lyr.GetFeatureCount()}) print(f 几何类型: {ogr.GeometryTypeToName(lyr.GetGeomType())}) if srs: print(f 坐标系: {srs.GetAuthorityName(None)}:{srs.GetAuthorityCode(None)}) else: print( 坐标系: 缺失需手动指定) # 抽第一个要素看坐标量级判断是经纬度还是投影 feat lyr.GetNextFeature() if feat: geom feat.GetGeometryRef() x, y geom.GetX(), geom.GetY() scale 经纬度 if abs(x) 180 and abs(y) 90 else 投影坐标 print(f 首点坐标量级: {scale} ({x:.2f}, {y:.2f})) for f in [jiangsu_river.shp, jiangsu_lake.shp, jiangsu_county.shp, jiangsu_highway.shp]: check_layer(f)这段脚本的价值在于把“感觉没问题”变成“有输出可核对”。参数上GetAuthorityCode返回 None 就说明.prj里没有 EPSG 编号需要人工确认首点坐标量级判断只是粗筛最终还是要看 SRS 字段。我一般把这套检查固化成习惯每批数据进来先跑一遍比事后返工省时间。再补一个参数习惯所有中间成果单独放一个目录文件名带坐标系后缀比如_4547避免不同版本混在一起。投影参数、编码、容差这些值记在一个文本里随数据走换人接手也不用猜。这套做法不复杂但能挡住大部分“数据看着对、一算就错”的问题。希望帮到你。本文还有配套的精品资源点击获取