简介广东省shp地图是一套面向GIS从业者、城乡规划人员及区域研究学者的矢量地理数据包聚焦广东省省界、市级与县级行政区划并附带产业与人口专题图层。数据采用Shapefile标准格式可直接在ArcGIS等专业软件中打开。资源共30个文件以shp几何文件为核心配合dbf属性表、shx几何索引、sbn/sbx空间索引、prj投影定义及mxd工程文档压缩包仅1.25MB小巧但结构完整产业专题还提供jpg预览图便于直观核对。目前已有2657人学习下载。用户拿到后可快速提取广东各级行政边界叠加人口、产业属性进行地图制作、空间可视化和初步分析适用于区域经济对比、人口分布呈现、公共服务设施布点等场景。由于数据包含投影文件坐标定义明确也适合作为高校GIS课程练习的数据底图。对于需要开展广东地理信息基础分析或课题研究的人而言这套数据具备较好的实用价值。1. 广东省shp地图省级行政区划数据不该让你从画边界开始做地图可视化或空间分析的人最怕的不是没有数据而是数据拿回来打开一片空白、要素全挤在坐标零点、或者属性表全是乱码。广东省shp地图就是用来终结这类问题的它是一套覆盖全省的行政区划矢量底图包含省界、地级市界、区县级边界可用于数据可视化、选址分析、人口与经济数据空间化、Web GIS原型搭建等场景。适合GIS开发、数据分析师、前端可视化工程师和做量化研究的同学。你不用再从零描边界更不用求人手工画底图但前提是搞清楚它的坐标系、字段编码和几何完整性——这正是我接下来要拆的内容。2. 拿到shp先做三件事查字段、查投影、查几何完整性2.1 用QGIS快速体检看层级和属性字段我拿到任何shp数据的第一件事不是写代码而是先拖进QGIS看一眼。这能立刻暴露八成问题坐标对不对、边界全不全、属性字段有没有值。广东全省的地图数据通常按行政层级拆分成多个文件常见的是“省界、地级市、区县”三组。在QGIS左侧图层面板里把三个图层都加进来先看两点一是缩放比例尺看省界是否有闭合、各区县是否有明显缺块。二是打开属性表看字段结构。一般行政区划shp会带“名称”“行政区划代码”这类属性有的还有“级别”“面积”等派生字段。若属性表全空那它就只是“一张画了边界的图”后续做空间关联会很痛苦。一份能用的广东区县级shp属性表里至少要有名称字段、行政区划代码字段。代码字段尤其关键它是空间数据与统计年鉴数据做关联的主键。没有代码字段我只能靠名称模糊匹配翻车概率极高。提示在QGIS中快速检查要素数量右键图层→属性→信息会显示要素总数。地级市约20个区县在100个以上如果数量对不上说明数据有缺失或抓取不完整。2.2 用GeoPandas查投影和几何不能被“看起来正常”骗了QGIS能显示不代表坐标可靠。我见过太多shp在QGIS里显示正常导出给后端却发现坐标是经纬度和投影坐标混用的。所以第二步用GeoPandas做程序化检查把坐标系和几何类型打印出来。import geopandas as gpd gdf gpd.read_file(guangdong_county.shp) print(gdf.crs) # 查看坐标系 print(gdf.geometry.geom_type.unique()) # 几何类型 print(gdf.geometry.is_valid.sum(), /, len(gdf)) # 有效几何数量 print(gdf.bounds.head()) # 查看每个要素的经纬度范围这段代码做了四件事第一行打印坐标系定义正常情况应该能看到EPSG:4326或其他坐标系描述如果输出None说明shp缺投影文件第二行看几何类型行政区划应该是Polygon或MultiPolygon出现LineString说明数据有问题第三行统计有效几何数量能初步暴露自相交问题第四行打印要素边界范围广东的大致范围北纬20°到25.5°、东经109°到118°左右如果出现0或负数大值别急着分析先去查坐标定义。参数上要注意read_file默认按文件实际内容读取并不会自动修正坐标系。crs是None时后续to_crs做投影转换会直接报错。遇到这种情况我一般先查数据源说明文档确认原始坐标系是WGS84还是CGCS2000再手动指定。2.3 shp是一组文件不是单个文件缺组件会导致读不出来这是非常基础的常识但也是下载资源最容易翻车的点。shp格式是“多文件结构”一个完整的shapefile至少需要四个文件配合才能正常打开文件后缀作用缺失后果.shp存储几何坐标无法显示图形.shx几何索引部分软件打不开.dbf属性表数据属性全丢.prj投影坐标系定义坐标乱掉或无法转换.cpg属性表字符集说明中文乱码在命令行下可以直接用ls命令核对文件完整性ls -l guangdong_county.*输出里应该看到至少四个同主名文件。如果只有.shp和.dbf那必须回去补全组件。常见做法是用QGIS“另存为”重新生成一份完整的shapefile或者用GDAL的ogr2ogr命令做一次格式转换。ogr2ogr -f ESRI Shapefile guangdong_county_fixed.shp guangdong_county.shp这条命令会读取原始数据重新产出全套shp组件文件。–f参数指定输出格式为ESRI Shapefile输出文件名换一个避免覆盖原始数据。执行完后再跑一遍ls命令文件族就齐全了。这里有个容易误会的点很多人以为.shp一个文件就是全部这是不对的。你从网上下到的数据如果只有一个.shp那基本可以判断是某平台的数据导出工具偷懒只发了几何文件属性数据已经丢了。3. 把shp用起来区划合并、坐标转换与导出3.1 先用GBK或UTF-8把dbf属性读出来广东省的行政区划shp来自不同源头属性表编码很分裂有的是GBK有的是UTF-8。GeoPandas读取时如果编码不对中文名称会变成乱码。import geopandas as gpd gdf gpd.read_file( guangdong_county.shp, encodingutf-8 ) print(gdf[[name, adcode]].head())假如打印出来的名称字段是“锟斤拷”或“鏉窞”说明编码猜错了。把encoding参数改成gbk再读一次gdf gpd.read_file( guangdong_county.shp, encodinggbk )多数情况下老数据的dbf用GBK编码的概率更大因为早期ArcGIS的数据生产环境默认就是GBK。新导出的开源数据大多用UTF-8。判断依据可以看.shp同目录下有没有.cpg文件——它记录的就是编码信息。查看它的内容cat guangdong_county.cpg输出UTF-8就按UTF-8读输出GBK或GB2312就按GBK读。如果连.cpg都没有那就两个编码都试一遍这是最笨但也最有效的办法。注意不要用Excel直接打开dbf文件查看它不会告诉你真实编码改文件还可能损坏结构。用GeoPandas或QGIS读最稳。3.2 按地市级合并区县dissolve的字段保留策略拿到区县级shp后有时需要把边界合并成地级市级别比如做市级人均GDP分布图。GeoPandas里用dissolve方法import geopandas as gpd gdf gpd.read_file(guangdong_county.shp, encodingutf-8) city_gdf gdf.dissolve( bycity_name, # 按地市名字段分组 aggfuncsum, # 非几何字段聚合方式 as_indexFalse ) print(city_gdf.shape)dissolve的核心逻辑按分组字段合并几何把同一地市下所有区县边界融合为一个面。这里最容易犯的错是aggfunc参数设置。如果属性表里有面积字段用sum会重复累加——每个区县的面积加总成地市面积这没问题。但如果有“人口”字段sum也是对的。可如果字段是“人均GDP”这种派生指标sum就完全错了合并后必须重新计算。我一般会在dissolve之前先把不需要的字段全部drop掉只保留name、adcode这类后续用得上的标识字段然后aggfunc设为first避免数值字段被错误累加。3.3 坐标转换到CGCS2000七参数和四参数该信谁广东的shp数据最常见的两种坐标系是WGS84经纬度和CGCS2000经纬度。前者是GPS原始坐标基准后者是国家大地坐标系。很多制图规范要求输出CGCS2000所以坐标转换绕不开。gdf_cgcs2000 gdf.to_crs(EPSG:4490) print(gdf_cgcs2000.crs)EPSG:4490就是CGCS2000地理坐标系单位是度。如果只是经纬度层面转换因为WGS84与CGCS2000在民用精度下差异极小理论框架差异约几厘米到亚米级直接原地切换即可肉眼完全看不出差别。但如果你拿到的shp是投影坐标比如地图厂商导出过的高斯-克吕格投影那就要小心了。投影坐标转回地理坐标必须先确认数据落在哪个分带3度带还是6度带中央经线是多少。分带搞错图形会横向错位几十公里。遇到这种数据我一般优先用GDAL做转换ogr2ogr -t_srs EPSG:4490 output.shp input.shp -s_srs EPSG:4547这里-s_srs指定源坐标系EPSG:4547是CGCS2000 3度带高斯投影的某分带代称。实际分带要看数据范围广东区域常用120度中央经线那一带。参数不对转换结果就是整个边界变成一条斜线非常直观。3.4 导出成GeoJSON/PostGIS给上层应用用后端分析和前端可视化需要的格式通常不是shp。shp是传统GIS桌面软件的格式浏览器里没法直接用。所以最后一步通常是转成GeoJSON。gdf_cgcs2000.to_file( guangdong_city.geojson, driverGeoJSON, encodingutf-8 )导出后我习惯先打开GeoJSON文件看一眼头部确认坐标数组不是空壳属性字段是中文原文不乱码。如果数据量大超过10MB前端渲染会卡这时优先考虑转成TopoJSON或是接入PostGIS做空间查询只返回视野范围内的数据。另外如果项目是要部署到在线地图上务必确认坐标系与在线底图一致。绝大多数在线底图使用Web墨卡托投影对应EPSG:3857。直接从EPSG:4490转3857gdf_web gdf_cgcs2000.to_crs(EPSG:3857)但这只是坐标参考系的数学转换。另一个坑是地图厂商对国内坐标做了偏转处理直接用WGS84或CGCS2000的shp叠加在在线底图上会发现省界整体偏移几百米。这个问题不是转换能解决的后面避坑章节细说。4. 制图与可视化把边界落到一张能用的图上4.1 QGIS分级配色五分钟做出一张专题图数据检查无误后最快速的出图方式是QGIS。假设我要把广东各地市的某个统计指标做成等级图右键图层→符号系统→分级。在“分级”面板里选择要渲染的数值字段分级方法选“自然间断点”类别数量设为5左右颜色渐层用浅色到深色。分级符号化的核心参数是“模式”和“类别数”。自然间断点适合本身分布不均的数据能拉开层级差距等间距则适合数值均匀分布的指标。如果数据有极端值建议先看直方图把离群值单独处理。做行政区划专题图时还要注意边界线的视觉层级地市级边界线宽0.6到0.8线色灰色省界加宽到1.2县级边界只在需要局部展示时开启不然图面会显得很碎。提示出图前务必先用“检查几何有效性”工具跑一遍否则后续打印或导出PDF时边界可能会出现奇怪的破洞或重叠。4.2 在线渲染GeoJSON接入前端底图Web端展示行政区划数据最常见的选择就是Leaflet。把上面导出的GeoJSON直接加载进来fetch(guangdong_city.geojson) .then(res res.json()) .then(data { L.geoJSON(data, { style: function(feature) { return { color: #666, weight: 1, fillOpacity: 0.3 }; }, onEachFeature: function(feature, layer) { layer.bindTooltip(feature.properties.name); } }).addTo(map); });这段代码里有两个关键点。style回调控制样式web墨卡托底图下最好把fillOpacity调低避免遮挡底图底纹。onEachFeature中的bindTooltip把属性表中的名称字段挂到悬浮提示上前提是属性字段名与你GeoJSON里的properties键一致。如果数据量超过几万要素前端直接挂载GeoJSON会卡顿。这时候要把数据交给后端PostGIS做空间过滤前端只加载视野范围的数据。我用过Practical最顺的方式后端提供一个接口接收前端传过来的边界范围返回与范围相交的行政区要素。整套链路从shp到PostGIS的导入也很简单shp2pgsql -s 4490 -I guangdong_county.shp guangdong_county | psql -U postgres -d gisdb-s参数指定坐标系EPSG:4490-I参数会自动创建空间索引。导入后前端通过接口请求体验比直接加载几十MB的GeoJSON好一个量级。4.3 标注避让名称重叠的烦恼行政区划图的标注是个容易显穷的细节。地市级名称放在面中心用“地市质心”标注区县名称同理。但广东沿海区县密集小面积区县的中心点会互相重叠。在QGIS中解决标注重叠可以用“基于规则的标注”在放置设置里开启“禁止冲突”并把优先级设为按面积从大到小分配——面积大的先标面积小的自动让位或隐藏。这是出图时我每次都强制检查的选项。5. 广东省shp地图避坑实录乱码、边界漂移与拓扑错误5.1 属性表中文变成“锟斤拷”乱码现象在QGIS或GeoPandas中打开广东shp名称字段显示为“锟斤拷”“鏉窞”等奇怪字符。原因dbf属性表是GBK编码但读取时默认按UTF-8解码中文二进制字节被错误解析。这是shp数据历史遗留问题老数据生产环境多为GBK。解决先看同目录.cpg文件确认编码或在读取时指定encodinggbk。如果是QGIS打开在添加图层时选择“图层编码”为GBK。已经乱码的数据不要直接在Excel里改重新用正确编码读一遍最干净。5.2 加载后图形跑到海里或消失现象shp在QGIS里加载后图层显示完全空白或者边界跑到非洲几内亚湾附近。原因数据是投影坐标但软件按经纬度坐标渲染或者.prj文件缺失导致软件默认按EPSG:4326处理。投影坐标值通常在数百万级别与经纬度的百级数值差了几万倍。解决先打印数据范围的min/max值如果数值在几十到一百多大概率是经纬度如果是六位数甚至七位数就是投影坐标。用ogr2ogr指定正确的-s_srs和-t_srs重投影。我处理的模拟项目X里遇到过一次边界跑到大洋中央的原因就是源数据的投影分带参数写错换中央经线后恢复正常。5.3 shp只有一个文件缺了一整套组件现象下载的“广东省shp地图”解压后只有一个.shp文件拖进QGIS报“无效数据源”或打开后只有图形没有属性。原因资源发布者只上传了主文件没有把.shx、.dbf、.prj一起打包。很多网盘打包工具会自动过滤掉“看起来像隐藏文件”的同名文件导致组件缺失。解决优先联系发布者要完整包。如果拿不到完整包用ogr2ogr强制读取.shp并重新输出一套完整文件。如果.shp本身已经缺了.dbf那属性数据大概率已经彻底丢失只能靠图形做纯展示用途。5.4 合并后出现细长面或自相交现象用dissolve把区县合并成地市后生成的面出现自相交或在计算面积时报错甚至某些地市边界多出极细长的“缝”。原因原始区县边界在相邻处存在微小重叠或间隙或者个别要素是无效几何比如自相交多边形。dissolve时这些无效几何会被带进合并结果。解决在dissolve前跑一遍几何修复from shapely.validation import make_valid gdf[geometry] gdf.geometry.apply(make_valid)make_valid会重构无效几何自相交的面会被拆成合法多边形。修复后重新跑dissolve生成结果再做一次is_valid检查。如果几何数量不多也可以在QGIS里用“修复几何”工具手工处理。5.5 叠加在线地图总是偏几百米现象shp和在线地图都是“经纬度”但叠加后省界线与影像底图总是对不上整体偏移。原因国内在线底图厂商使用了偏转坐标算法把WGS84坐标加密到偏移后的坐标系。直接叠加未经转换的原始坐标必然偏移。解决两种办法。一是换底图使用无偏转的标准底图二是在应用中引入坐标纠偏库把WGS84的经纬度在客户端做一次偏转转换。这个问题本质上是坐标系基准不统一不是shp数据本身的问题。注意在处理专题图叠加影像做对比时最好确保两边坐标系基准一致否则“看着是同一个点实际隔了一条街”。6. 交付前做三个验证几何修复、面积自检与影像底图比对最后这步是我从多次翻车换来的习惯所有shp数据在用于正式项目前都要跑一遍几何有效性、面积量算和视觉比对三重验证。很多人拿到数据直接做分析最后出图全是破洞或边界错位回头改数据反倒浪费两倍时间。第一重验证是几何有效性。用GeoPandas跑is_valid统计如果无效几何比例超过千分之一就用make_valid修复再导出新文件。这一步能避免后续空间计算和出图时出现莫名其妙的异常。import geopandas as gpd gdf gpd.read_file(guangdong_county.shp, encodingutf-8) invalid gdf[~gdf.geometry.is_valid] print(f无效几何: {len(invalid)})第二重验证是面积自检。把地理坐标投影到适合广东区域的等积投影再计算各要素面积与官方统计年鉴上的行政区面积做对比。偏差超过2%就说明数据有问题可能是边界缺失或几何错误。广东海岸线曲折沿海区县如果漏了填海区域面积偏差会很明显。第三重验证是视觉比对。把shp叠加到无偏转的影像底图上透过半透明边界检查主要城镇轮廓是否吻合。行政边界与影像中明显的地物特征如河口、山脊线不应有肉眼可见的偏移。这一步靠代码无法替代必须人眼抽查三到五个区域。从那以后我每次拿到shp数据都强制走一遍这三个验证哪怕再急也至少跑完前两个。这套流程救过我不少次尤其是那些“看起来全对”的数据偏就偏在不起眼的多边形角落里。希望帮到你。本文还有配套的精品资源点击获取
