简介粤港澳大湾区“92”政区图是一套按省级权威底图整理的ArcGIS矢量数据覆盖广州、深圳、珠海、佛山、惠州、东莞、中山、江门、肇庆九座内地城市以及香港、澳门两个特别行政区。数据依据广东省自然资源厅发布的1:100万行政区划图重新制作作者并未采用现成转载版本而是花费大量时间处理海岸线和边界细节尤其对珠江口周边岸线进行了手工校正适合用于国土空间规划、人口经济数据关联、论文配图、课堂教学等场景。压缩包体积仅270KB共7个文件其中shp是几何主文件dbf保存行政区名称与代码等属性prj定义投影坐标系统sbn/sbx提供空间索引可提升编辑与查询效率xml则包含元数据信息解压后可直接在ArcGIS、QGIS等工具中加载、编辑或出图。当前已有4365人浏览学习无论是GIS初学者快速获取大湾区底图还是研究者在此基础上结合业务数据进行专题分析这套数据都能显著减少前期整理与校正工作量。1. 先看清“92”是什么一份政区图压缩包能省掉多少事做粤港澳大湾区相关的项目最磨人的往往不是分析本身而是找一张能用的底图。所谓“92”指的是珠三角九座城市加上香港、澳门两个特别行政区全称就是“粤港澳大湾区政区图92”。这一类压缩包解压以后一般是一整套行政区划矢量数据包含九市两区的边界、市名、区县名和常见属性字段。它的价值不在于“画好了一张图”而在于把分散各处的政区边界收敛成一个数据源不用再为广州找一个文件、深圳另找一个文件坐标系还不一致。适合做城乡规划、人口统计、物流可视化、高校课题区位图的开发者。但压缩包解开只是第一步投影、编码、合并、校验这些事一步都省不了。下文按“先检查数据、再导入软件、再统一坐标、最后合并出图”的顺序展开每一步都给出命令、参数和坑。2. 解压后的第一件事核对文件清单、字段和坐标系rar 解压之后别急着把 shp 拖进 GIS。先花十分钟搞清楚三件事文件是否完整、属性表里有什么、坐标系是什么。这三件事决定了后面每一步的操作方向跳过任何一个后面的报错都要回来补课。2.1 认清 shapefile 文件家族不是单个 shp 就有全部信息一个可用的 shapefile 是“一组文件”而不是“一个文件”。其中 .shp 存几何边界.dbf 存属性表.shx 是几何和属性之间的索引这三者缺一不可。很多政区图压缩包还会带 .prj 和 .cpg前者记录坐标系描述后者记录属性表编码。对做政区图的人来说.prj 几乎是整个包里最重要的文件没有它GIS 只能按默认坐标系去猜数据位置。以下是常见文件后缀及其缺失后果建议对照压缩包逐个检查。文件后缀作用缺失后果.shp几何数据面或线边界没有实际图形.dbf属性表存市名、区县名等图层开得出来但看不到名称.shx几何与属性间的索引部分软件直接拒绝打开.prj坐标系描述文本投影错乱、位置偏移.cpg属性表编码声明中文属性乱码注意一个常见的误判在 Windows 里直接双击 .shp 文件系统弹窗说“无法打开”不代表数据损坏。换个思路把整个文件夹拖进 QGIS或者用下面的命令行工具先读一遍再决定是否真的缺文件。2.2 用 ogrinfo 快速检查图层边界与坐标系我拿到任何一份政区图第一件事就是打开终端跑 ogrinfo。它是 GDAL 自带的命令行工具安装 QGIS 时一般会一并装好也可以在终端单独安装 GDAL 后使用。命令很简单进入解压目录指定要看的 shp 文件即可。cd /path/to/your/shp_dir ogrinfo -so -al guangzhou.shp输出里需要重点看的四行是Geometry 告诉你几何类型是 Polygon 还是 MultiPolygonFeature Count 是要素数量对应区县数量Extent 给出数据范围如果数值是 113 到 114 之间的经纬度量级说明是地理坐标系如果变成几十万上百万的数量级说明是投影坐标系Layer SRS 如果有一长串坐标系描述说明 .prj 存在。最怕的是输出里没有 SRS 这一行那意味着这份数据没有坐标系定义后面叠加在线底图时大概率要出问题。这条命令是只读操作不会改动原始文件可以反复执行。如果目录里有十个城市的 shp可以写一个 for 循环批量跑把每个文件的 Extent 和 SRS 列出来几分钟内就能完成整个压缩包的体检。2.3 从命名和属性字段判断这份数据是市界还是区县界政区图包常见两种组织方式一种是按市拆分文件比如 gz.shp、sz.shp、hk.shp 各管一片另一种是把九市两区全部揉进一个文件靠属性字段区分城市。后者的文件名常带“92”“dawanqu”一类的字样。判断粒度最直接的方法是看属性表有“地市名称”和“区县名称”两个字段的是区县级数据只有“地市”字段的通常只有市界线。这个判断直接决定了后面能不能做按市聚合的统计。如果拿到的只是市界级数据就别指望它包含乡镇街道边界如果要做九市的精细图还需要重新找更细粒度的来源。港澳数据在这个包里通常单独成文件命名常见的是 hk.shp、mo.shp 或 gangao.shp。它们的位置、范围、坐标系和内地文件往往不一致后面要单独处理。3. 装进 QGIS 和 ArcGIS空白、偏移和路径断链的三种处理文件体检通过之后进入软件才真正接触数据。最常见的异常无外乎三种图层放进去界面空白、位置偏移到远处、ArcGIS 工程文件再次打开数据源失效。下面逐个说。3.1 QGIS 打开后看不到图先指定图层 CRS别急着重投影在 QGIS 里拖入一个政区 shp结果什么都没有这是碰到的最高频问题。多数情况下数据本身没坏而是 QGIS 不知道它应该落在哪里。默认情况下 QGIS 会拿 WGS84 经纬度去解释所有坐标如果这个文件原本是投影坐标数值量级完全不同图形就会画到视野之外。处理的关键是分清两个概念指定 CRS 和重投影。指定 CRS 是告诉软件“这份数据原本在哪个坐标系下”好比告诉别人这张图纸用的是哪套尺子重投影是把坐标数值换算到另一个坐标系相当于把尺子换成另一套。很多人上来就点“导出-另存为-重投影”如果原始坐标系本来就没定义重投影只会错上加错。正确顺序是先右键图层 → 图层属性 → 源 → 关联 CRS选对原始坐标系再把图层缩放到全图范围然后导出重投影。如果一次要处理十几个城市文件用 PyQGIS 脚本批量赋 CRS 更省事。下面这段代码的作用是给单个图层声明坐标系不改坐标数值from qgis.core import QgsVectorLayer, QgsCoordinateReferenceSystem, QgsProject layer QgsVectorLayer(gz.shp, gz, ogr) crs QgsCoordinateReferenceSystem(EPSG:4490) # CGCS2000 地理坐标系 layer.setCrs(crs) QgsProject.instance().addMapLayer(layer)参数说明EPSG:4490 是 CGCS2000 大地坐标系的代码适合国内测绘数据如果确认数据是 WGS84就用 EPSG:4326搞不清楚时先别写去查这份数据的来源说明。setCrs 只是声明等于给数据贴上坐标系标签。声明完之后正常显示只是第一步真正的坐标换算还要靠下一章的 ogr2ogr 或导出操作。3.2 ArcGIS 工程文件打不开数据源路径断链的修复另一种常见事故发生在 ArcGIS 里第一次打开数据正常保存工程关闭第二天再打开提示数据源缺失。原因通常是工程文件记录的是绝对路径而数据所在目录有了变动或者路径里带了中文导致解析失败。这个问题和政区图本身无关但一旦发生在交付阶段非常尴尬。修复分两步。第一步打开“自定义 → ArcMap 选项 → 常规”勾选“存储数据源的相对路径名”这样以后移动工程文件夹时数据还能跟着走。第二步对已经断链的图层右键图层 → 数据 → 修复数据源重新指向对应的 shp 文件。习惯上我会把解压后的数据全部放到同一个根目录比如项目文件夹下的 data 子目录保证相对路径始终成立。值得留意的是ArcGIS 对 shapefile 的字段名有长度限制属性表里的中文字段名在导出时可能被截断或省略。做政区图交付时建议字段名统一用拼音或英文把中文名称放到字段值里这是规避一系列乱码和断链问题的经验做法。3.3 港澳两个图层为什么总对不上“92”数据包里香港和澳门两个文件往往是特殊的存在。香港政区数据常用香港 1980 网格坐标系澳门也有自己的独立大地基准跟内地的 CGCS2000 或 WGS84 不通用。直接把 hk.shp 拖进工程里和内地九市叠加结果经常是香港跑到深圳北边或珠江口中央澳门凭空悬在海面上——这就是典型的坐标系不统一的“玄学”。判断方法看范围单独加载 hk.shp如果状态栏显示的范围是几十万上百万量级而内地文件是 113.x、22.x 这种经纬度格式两者不在一个坐标系下。更隐蔽的情况是两者都声明为经纬度但使用的椭球参数不同肉眼看到的是几百米级别的漂移放大后沿着海岸线一看就能发现。处理港澳数据的正确姿势不是手动改数字也不是在图层属性里顺手换一个 CRS而是分别查明原始坐标系再统一重投影输出到同一目标坐标系。具体命令放在下一章这里先记住一个原则港澳图层一定要单独处理不能跟着内地文件一起批量重投影因为它们的源坐标系不同一起处理等于拿错参数硬套。4. 坐标转换与底图叠加GCJ-02、WGS84、CGCS2000 怎么选政区图值不值钱一半看边界精度另一半看坐标系是否统一。这一章把几套坐标系的关系、判断方式以及统一转换步骤讲清楚顺便给命令和代码。4.1 政区图坐标系乱的根源国内能见到的政区图数据坐标系大体分成四类。CGCS2000 是国家大地坐标系测绘成果的标准格式WGS84 是 GPS 原始坐标系也是很多国际数据源的默认值GCJ-02 是互联网地图常用的加密坐标再老一点还有西安 80、北京 54。CGCS2000 和 WGS84 在省级尺度上的差异通常在一米以内出图制图基本不需要处理。真正的坑集中在 GCJ-02它不是单纯平移而是经过非线性加密整个湾区如果混用边界偏移几百米很常见而且面积统计也会跟着失真。判断政区图是不是 GCJ-02最直接的办法是看 .prj 里的描述词。打开 .prj 文件如果出现“GCJ”“Mars”这类关键词多半就是加密坐标。还有一部分数据源会直接注明“火星坐标”或“高精度坐标”这类同样属于 GCJ-02 系统不要和 WGS84 混着算。坐标系常见来源与 WGS84 差异制图建议CGCS2000测绘成果、国土数据米级可用面积统计推荐WGS84GPS、国际数据源基准一致叠加国际底图首选GCJ-02互联网地图几百米级非线性偏移避免与其它数据混用西安80/北京54历史数据十米级以上需七参数转换后再用4.2 不读元数据也能判断坐标系的两个土办法当 .prj 被删掉、元数据也找不到的时候我用两个土办法兜底。第一个是地标点反推找一个位置已知的地标比如某城市电视塔查它的公开经纬度再在图层里找到同一个地标的坐标对比两者的差值。如果整体偏移零点几度基本可以断定是加密坐标或坐标系选错。第二个是叠加在线底图加载一份卫星底图把政区图叠上去看珠江口轮廓、岛屿、岸线能否对齐。对不齐的时候先怀疑坐标系再怀疑数据本身的精度。需要注意在线底图服务商使用的坐标系并不相同叠加时要先确认底图用的是什么坐标系。多数标准在线瓦片基于 WGS84 或它的投影形式而某些在线地图默认是 GCJ-02叠加之前如果不弄清楚反而会误判自己的数据。4.3 用 ogr2ogr 和 geopandas 统一到目标坐标系统一坐标系的落地方式取决于最终用途。叠加在线底图出图用 EPSG:4326 或 EPSG:3857做面积计算和空间统计用 CGCS2000 的分带投影坐标。下面这条命令把单个政区文件重投影到 WGS84 经纬度ogr2ogr -t_srs EPSG:4326 wgs84_dir/gz.shp raw_dir/gz.shp参数说明第一个参数是输出文件第二个是输入文件-t_srs 指定目标坐标系。这里最容易翻车的是把 -t_srs 和 -a_srs 搞混。-t_srs 是真正执行坐标换算-a_srs 只是把坐标系声明写进文件不移动任何坐标。用错 -a_srs 的结果是文件“看起来有了坐标系”但叠到底图照样漂移。批量处理整个压缩包时我一般用 geopandas 写循环比一条条敲命令直观还方便顺便处理字段。import geopandas as gpd cities [guangzhou, shenzhen, zhuhai, foshan, huizhou, dongguan, zhongshan, jiangmen, zhaoqing, hk, mo] for name in cities: gdf gpd.read_file(fraw/{name}.shp) gdf gdf.to_crs(EPSG:4326) gdf.to_file(fwgs84/{name}.shp, encodingutf-8)这段代码把十一个文件全部读入、转坐标系、写出。两个注意点其一to_crs 的前提是源文件坐标系定义正确如果源 .prj 本身就是错的转完照样是错的其二输出时显式指定 encodingutf-8能提前拦住后面的属性表乱码问题。港澳两个文件在循环里虽然一起转但源坐标系不同所以转换结果需要单独核对不要因为代码跑完了就默认没问题。5. 避坑政区图落地常见的 5 个翻车点这一章的每一条都是实际项目里撞过的墙统一按“现象 → 原因 → 解决”写。遇到对应问题可以直接跳到对应小节处理。5.1 属性表汉字乱码编码问题挡住九成新手现象打开属性表市名和区县名变成“”或一堆乱码英文字段正常。 原因shapefile 的属性表 dbf 大多以 GBK 或 ANSI 编码存储而 QGIS 默认按 UTF-8 读取。ArcGIS 在中文系统里通常正常一旦工程文件移到英文环境同样的问题也会出现。 解决在 QGIS 中打开数据源管理器选中图层后指定编码为 GBK重新加载用 geopandas 读数据时显式传 encodinggbk再转存成 UTF-8import geopandas as gpd gdf gpd.read_file(guangzhou.shp, encodinggbk) gdf.to_file(guangzhou_utf8.shp, encodingutf-8)改完编码后务必检查一个中文地名确认能正确显示再继续后面的合并和出图。这个检查只需要十秒能省掉后面所有属性关联环节的返工。5.2 九市合并出现缝隙与重叠拓扑修复不是可选项现象把九市 shp 合并成一个文件后相邻市边界处出现白色缝隙或者存在互相重叠的窄带。直接计算面积总和对不上湾区陆地总面积。 原因各市的边界线来自不同批次或不同数据源相邻市界上的节点坐标不一致面与面之间自然形成微小缝隙或重叠。 解决合并之后做拓扑修复。QGIS 里可以用“修复几何”配合“消除”处理缝隙ArcGIS 里用拓扑工具集规则指定为“不能有间隙”和“不能重叠”修复后再跑统计。如果只是出图展示白缝可以用地图背景色掩盖但以后一旦要做面积统计或空间查询这条不能省。5.3 港澳图层整体偏移坐标系没统一就叠加现象内地九市显示正常香港、澳门的边界单独看也没问题一旦叠加进同一工程就整体错位。 原因港澳文件与内地文件的坐标系不一致加载时又没有先核对坐标系就直接混放。不少压缩包里的港澳数据还是从在线底图勾出来的精度本身偏低转换后仍会差几十米。 解决按第 4 章的 ogr2ogr 命令对港澳文件分别查明原始坐标系后重投影到统一目标坐标系再做叠加。港澳数据转换完要放大到口岸、跨海通道近岸这些边界线密集的位置检查确认与九市边界衔接自然再进入下一步。5.4 压缩包缺 .prj 或缺字段拿到的数据不完整现象明明解压出一堆文件导入 GIS 后提示“未知的空间参考”或者图层能显示但属性表缺关键的市名、区县名字段。 原因部分政区图包在整理上传时为了压缩体积或表格处理方便删掉了 .prj、.shx 这类看似次要的文件导致坐标系信息丢失。 解决先看同包内其它城市是否有 .prj多数情况它们来自同一数据源坐标系一致把参考 .prj 复制过来改名即可。另一个更稳的做法是在 QGIS 里按已知城市中心点手动指定 CRS再另存为完整 shapefile。属性字段缺失只能回到压缩包元数据里找说明找不到就先按当前可用的字段推进出图时改用手工标注。5.5 Dissolve 后字段丢了大半合并前想清楚要留什么现象在 ArcGIS 里用 Dissolve 把区县面融合成市界面结果输出图层只剩下几何和统计字段城市名称、区县名称全没了。 原因Dissolve 默认只保留用于融合的分组字段和统计字段其它字段不会自动带过去。 解决在 Dissolve 对话框的“字段映射”区域里把需要保留的字段合并规则改成“第一条”或者按需选择“总和”“平均值”。QGIS 的“融合”工具也一样先在字段列表中勾选分组字段再配置输出字段。这一步做完再补一个“城市”字段用于配色和标注就顺手多了。6. 把九市两区合并成一个可用图层字段整理与出图前校验走到这一步数据已经能显示、坐标系也统一了。最后收个尾把九市两区合并成一个文件补上城市字段再做出图前校验。这也是交付给别人时最省心的形态。6.1 批量合并与统一坐标系一条命令跑完全部城市ogr2ogr 支持追加模式。先以第一个城市建目标文件后面每个城市用 -update -append 追加ogr2ogr -f ESRI Shapefile merged.shp guangzhou.shp -t_srs EPSG:4326 ogr2ogr -f ESRI Shapefile -update -append merged.shp shenzhen.shp -t_srs EPSG:4326 ogr2ogr -f ESRI Shapefile -update -append merged.shp hk.shp -t_srs EPSG:4326每个城市一条-t_srs EPSG:4326 保证追加进去的所有数据坐标系一致。全部跑完后用 ogrinfo 查看要素数量九市两区一共是几十个区县单元数量应该对上对不上就回头查是哪个文件没追加上别急着出图。6.2 补一个“城市”字段让配色和标注不再靠猜合并后的文件通常没有城市字段配色只能一个颜色到底。补字段用 geopandas 最直接按区县名称映射到城市再写回文件import geopandas as gpd gdf gpd.read_file(merged.shp, encodingutf-8) district_city {越秀区: 广州, 天河区: 广州} # 示意按实际属性替换 gdf[城市] gdf[区县名称].map(district_city) gdf.to_file(merged_city.shp, encodingutf-8)实际数据里的区县名跟示例未必对得上先跑一条 print 把不重复的区县名列出来再写映射表。漏写一个区县就会多出一个“其他”分类配色阶段才发现就要返工。6.3 出图前最后校验三处边界看完就能放心用个人经验是放大到三处跨市地物密集的位置检查一处跨市河口的入海口、一处跨海通道的东岸落点、一处陆地口岸两侧。这三个位置边界交错最容易暴露坐标系不一致、缝隙和缺边问题。每处都要对比九市与港澳的边界是否自然衔接再叠加一张在线底图确认岸线轮廓。我给自己定了个老规矩任何一份政区图先跑 ogrinfo、先看 .prj、先叠在线底图三件套做完才动手分析。这套流程不是流程本身而是给自己留后悔药——坐标系问题越早发现返工成本越低。希望你拿到的下一份“92”政区图能一次跑通、少踩两个坑。希望帮到你。本文还有配套的精品资源点击获取