简介这份合肥市乡镇街道行政区划shp文件面向GIS从业者、城市规划人员及地理信息相关专业师生提供乡镇街道级别的矢量边界数据可用于空间分析、专题制图与城市精细化管理研究。资源包共14个文件约445KB以shp主文件为核心配套shx索引、dbf属性表、prj投影信息及sbx、sbn、cpg等辅助文件另含xml元数据与lock锁定文件构成完整可用的地理数据集。目前已有1262人学习下载说明其在同类数据中具备一定参考价值。通过该数据可获取各街道乡镇的精确边界与属性字段便于导入ArcGIS、QGIS等软件开展人口分布、交通流量或公共服务设施布局分析也可作为Web GIS在线地图服务的基础底图适合需要合肥本地精细化行政区划底图的中高级用户直接调用。1. 合肥市乡镇街道行政区划 shp 文件一份能直接进 GIS 的矢量底图做合肥本地项目的人多半遇到过这个场景领导丢过来一句“把合肥各乡镇的边界画出来”或者要按街道统计人口、门店、管网、地块结果手头只有一张模糊的截图或者一份属性表里写着“XX街道”却没有几何。这时候真正卡住你的不是分析模型而是缺一份能直接拖进 ArcGIS、QGIS 或 PostGIS 的合肥市乡镇街道行政区划 shp 文件。这份资源就是干这个的——它把合肥市下辖各区县的乡镇、街道级行政边界整理成标准 shp 矢量数据每个面要素带名称、编码等属性字段坐标系通常是地理坐标CGCS2000 或 WGS84可以直接用于制图、空间连接、裁剪和统计。适合做本地政务、规划、物流、商业选址、数据可视化的从业者也适合刚接触 GIS 矢量数据、想拿一份真实行政区划练手的新手。下面我按“拿到文件先看什么、怎么用、哪里会翻车”的顺序拆一遍。2. 拿到 shp 先别急着打开字段、坐标系与拓扑自检shp 不是单个文件而是一组同名文件捆在一起工作的。很多人第一次拿到压缩包解压后只看到.shp就双击结果要么打不开要么打开是空白问题基本都出在配套文件缺失或坐标系没对上。这一章先把“这份数据到底能不能用”判断清楚再谈后面的操作。2.1 shp 文件组的构成与最小可用集合一个完整的 shapefile 至少包含三个文件.shp存几何、.shx存索引、.dbf存属性表。少任何一个ArcGIS 会直接报错QGIS 可能能读几何但属性全空。实际分发时还常带.prj坐标系定义、.cpg字符编码、.sbn/.sbx空间索引。合肥这份乡镇街道数据你解压后应该能看到一整套同名文件文件名一致、只差扩展名。文件作用缺失后果.shp几何坐标无法显示图形.shx几何索引部分软件报错或读取异常.dbf属性字段有图无属性无法按名称筛选.prj坐标系定义坐标数值在但不知道是什么坐标系.cpg字符编码中文名称乱码判断能不能用的第一步确认这五个文件齐全且文件名完全一致。常见坑是解压工具把中文名截断导致.dbf和.shp名字对不上软件就认不出来。2.2 坐标系识别为什么你的图和底图对不上坐标系是 shp 使用里最容易翻车的一环。合肥市位于东经约 117°、北纬约 31.8°如果数据是地理坐标系CGCS2000 或 WGS84坐标值应该是117.xx, 31.xx这种小数如果看到的是500000, 3400000这类大数说明它是投影坐标系常见 CGCS2000 3 度带中央经线 117°E。在 QGIS 里查看方式右键图层 → 属性 → 信息看 CRS 一栏。在 ArcGIS 里图层右键 → 属性 → 源看“空间参考”。如果.prj缺失软件会提示“未知坐标系”此时图形能显示但一旦和在线底图叠加就会偏移几百米甚至几公里。提示不确定坐标系时先加载一份已知坐标系的底图如天地图或本地正射影像把 shp 叠上去看是否吻合。吻合说明坐标系一致整体平移说明是基准差异缩放不对说明是投影 vs 地理坐标混用。2.3 拓扑自检乡镇边界该不该严丝合缝行政区划面数据理论上应该无缝无重叠——相邻乡镇共享一条边界全市拼起来正好是合肥市域。但实际拿到的数据常有缝隙、重叠、悬挂节点。做统计汇总时重叠会导致同一块地被算两次缝隙会导致漏算。快速自检方法QGIS矢量 → 几何工具 → 检查有效性能查出无效几何再用“拓扑检查器”插件查重叠和缝隙。ArcGIS 用“拓扑”工具集建立“不能有重叠”“不能有缝隙”规则。如果只是做展示轻微缝隙不影响如果要做面积统计或空间连接建议先修复。修复思路用“融合”按区县字段合并再“拆分多部件”或者用“消除”工具处理小缝隙。这一步不做后面按乡镇统计人口或 POI 时数据对不上排查起来非常痛苦。3. 从 shp 到可用成果筛选、裁剪与空间连接实操数据自检通过后真正的活是把它变成你要的结果。这一章覆盖三个最高频操作按名称筛出目标乡镇、按范围裁剪出研究区、把业务数据挂到行政区上做统计。每个操作给可复现的步骤和参数说明。3.1 按属性筛选目标乡镇街道假设你只要合肥市包河区的乡镇街道。在 QGIS 里右键图层 → 筛选 → 输入表达式区县 包河区字段名要换成你数据里实际的区县字段可能是NAME、区县名称、COUNTY等。先打开属性表看一眼字段名再写别凭猜。ArcGIS 里图层右键 → 属性 → 定义查询 → 输入同样的 SQL 表达式。如果字段里存的是编码而不是名称比如340111代表包河区那就用编码筛选。筛选完建议“导出 → 保存选定要素”生成一个新的 shp避免后续操作误改原数据。参数说明表达式里的字符串用单引号字段名用双引号QGIS或不加引号ArcGIS 视版本而定。中文乱码时先确认.cpg文件内容是否为UTF-8不是就改成 UTF-8 再重新加载。3.2 按研究区范围裁剪如果你手上有一个更大的范围比如整个安徽省的乡镇数据只想保留合肥市部分用“裁剪”工具。QGIS矢量 → 地理处理工具 → 裁剪输入图层选安徽数据叠加图层选合肥市界输出即合肥范围。ArcGIS分析工具 → 提取 → 裁剪注意“XY 容差”保持默认勾选“使用输入要素裁剪几何”。# 用 geopandas 做同样的裁剪适合批量处理 import geopandas as gpd anhui gpd.read_file(anhui_towns.shp) hefei gpd.read_file(hefei_boundary.shp) # 统一坐标系避免裁剪结果错位 anhui anhui.to_crs(hefei.crs) hefei_towns gpd.clip(anhui, hefei) hefei_towns.to_file(hefei_towns_clip.shp, encodingutf-8)逻辑说明to_crs保证两个图层坐标系一致这是裁剪前必须做的一步否则结果可能为空或错位。gpd.clip按叠加图层的外框裁剪输入图层保留相交部分。encodingutf-8保证中文属性写出不乱码。参数说明如果裁剪后要素数量异常少检查两个图层的坐标系是否一致如果属性丢失检查字段名是否超过 dbf 的 10 字符限制shapefile 格式限制长字段名会被截断。3.3 空间连接把业务数据挂到乡镇街道上最常见的需求有一批门店点位或人口表格想知道每个乡镇街道各有多少。点位用“空间连接”表格用“按位置选择”或“连接”。QGIS 空间连接矢量 → 数据管理工具 → 按位置连接属性输入图层选乡镇面连接图层选点位几何谓词选“包含”输出即每个面带上落在其中的点属性。# geopandas 空间连接示例 import geopandas as gpd towns gpd.read_file(hefei_towns.shp) points gpd.read_file(stores.shp) # 坐标系必须一致 points points.to_crs(towns.crs) joined gpd.sjoin(points, towns, howleft, predicatewithin) # 按乡镇统计点位数量 count_by_town joined.groupby(乡镇名称).size().reset_index(name门店数) print(count_by_town.head())逻辑说明sjoin把每个点匹配到包含它的面predicatewithin表示点在面内。howleft保留所有点没匹配上的点乡镇字段为空方便排查。最后按乡镇名称分组计数。参数说明predicate可选within、intersects、contains点位统计一般用within。如果点位正好落在边界上within可能不匹配改用intersects更稳妥。统计前先检查有没有点位落在合肥市域外那部分会被归为空值。4. 避坑与排查shp 使用中最容易翻车的五件事这一章是我自己踩过和帮人排查过的真实问题按“现象 → 原因 → 解决”写。每一条都对应一个具体操作场景不是泛泛而谈。4.1 打开是空白或图形跑到地图外现象shp 加载后看不到图形或者图形出现在地图角落、和底图完全对不上。原因坐标系缺失或错误。.prj文件丢失时软件按默认坐标系处理导致位置偏移或者数据本身是投影坐标但被当成地理坐标显示。解决先看坐标数值范围。如果是117, 31级别是地理坐标如果是500000, 3400000级别是投影坐标。在 QGIS 里用“设置图层 CRS”手动指定正确坐标系再“导出 → 另存为”写入.prj。ArcGIS 用“定义投影”工具。指定后仍偏移说明基准不同如 WGS84 vs CGCS2000需要做基准转换。4.2 中文属性乱码现象属性表里乡镇名称显示成????或锟斤拷。原因.cpg文件缺失或编码声明与实际不符。shapefile 的 dbf 不强制声明编码软件只能猜。解决新建一个同名.cpg文件内容写UTF-8重新加载。如果还是乱码用 QGIS 的“重新加载编码”试 GBK。ArcGIS 里可以在目录窗口中右键 dbf → 属性 → 更改编码。批量处理时用 Python 读 dbf 指定 encoding 参数。4.3 面积统计结果偏大或偏小现象按乡镇统计面积结果和官方公布数据差很多。原因坐标系是地理坐标度直接算面积得到的是“平方度”没有物理意义。必须投影到等面积坐标系再算。解决QGIS 用“重投影图层”转到 CGCS2000 3 度带合肥用中央经线 117°EEPSG 视具体带号或用 Albers 等面积投影。ArcGIS 用“投影”工具。投影后再计算几何面积。# 投影后计算面积单位平方米 towns gpd.read_file(hefei_towns.shp) towns_proj towns.to_crs(EPSG:4547) # CGCS2000 3度带 117E towns_proj[area_m2] towns_proj.geometry.area towns_proj[area_km2] towns_proj[area_m2] / 1e6参数说明EPSG:4547是 CGCS2000 3 度带中央经线 117°E 的投影代码合肥适用。不同带号代码不同选错会导致面积变形。算完和官方数据对比偏差在合理范围内即可。4.4 空间连接后统计数量对不上现象明明有 1000 个点连接后按乡镇汇总只有 800 个。原因部分点落在市域外、落在边界上未被within匹配、或坐标系不一致导致匹配失败。解决先检查未匹配的点乡镇字段为空落在哪里用“按位置选择”看它们和市域的关系。边界上的点改用intersects。坐标系不一致的先统一。另外注意重复匹配一个点落在两个面重叠区会被算两次需要先修复面的重叠。4.5 导出 GeoJSON 或转其他格式后属性丢失现象shp 转 GeoJSON、转 Excel、转 CAD 后部分字段没了或名称被截断。原因shapefile 的 dbf 字段名限制 10 个字符超长会被截断转 GeoJSON 时字段类型不兼容如日期格式会丢失。解决导出前先把字段名改成 10 字符以内或改用 GeoPackage.gpkg作为中间格式它没有字段名长度限制。转 CAD 时注意 CAD 只认几何和少量属性复杂属性表建议单独导出表格再关联。5. 进阶把乡镇 shp 用进自动化流程与在线地图前面讲的都是单机操作这一章说两个进阶方向一是把这份数据接进自动化脚本批量出图或批量统计二是转成在线地图能用的格式做 Web 端展示。最后给一个我常用的验证习惯。5.1 批量出图按乡镇循环导出专题图做汇报时经常要给每个乡镇单独出一张图。手动一个个筛太慢用 Python 循环import geopandas as gpd import matplotlib.pyplot as plt towns gpd.read_file(hefei_towns.shp) towns_proj towns.to_crs(EPSG:4547) for name, group in towns_proj.groupby(乡镇名称): fig, ax plt.subplots(figsize(8, 8)) towns_proj.plot(axax, color#e8e8e8, edgecolorwhite) group.plot(axax, color#c0392b, edgecolorwhite) ax.set_title(name) ax.axis(off) plt.savefig(foutput/{name}.png, dpi150, bbox_inchestight) plt.close()逻辑说明先画全部乡镇作为底再高亮当前乡镇保证每张图的空间参照一致。groupby按名称分组循环导出。bbox_inchestight去掉多余白边。参数说明dpi控制清晰度汇报用 150 够印刷用 300。中文字体需要提前设置plt.rcParams[font.sans-serif] [SimHei]否则标题乱码。输出目录要先建好否则报错。5.2 转 GeoJSON 与在线地图展示Web 端展示常用 GeoJSON 或矢量瓦片。转 GeoJSON 最简单ogr2ogr -f GeoJSON hefei_towns.geojson hefei_towns.shp -lco ENCODINGUTF-8逻辑说明ogr2ogr是 GDAL 自带工具-f GeoJSON指定输出格式-lco ENCODINGUTF-8保证中文编码。转出来的 GeoJSON 可以直接被 Leaflet、Mapbox、OpenLayers 加载。参数说明如果文件很大乡镇级一般不大几百 KB 到几 MB可以用-simplify简化几何减少体积。在线展示时注意坐标系用 WGS84EPSG:4326GeoJSON 标准要求如此投影坐标需要先转回地理坐标。5.3 一个我每次都走的验证习惯不管拿到什么来源的行政区划 shp我固定走三步验证第一步看属性表字段和记录数和官方公布的乡镇街道数量对一遍数量对不上说明有缺失或重复第二步叠加已知底图看位置确认没有整体偏移第三步随便挑一个乡镇算面积和公开数据比偏差超过 10% 就回去查坐标系。这个习惯帮我拦下过好几次“看起来能用、一算就错”的数据。从那以后我每次拿到新的 shp 都强制走一遍这三步宁可多花十分钟也不在错误数据上做半天分析。希望帮到你。本文还有配套的精品资源点击获取
简介这份合肥市乡镇街道行政区划shp文件面向GIS从业者、城市规划人员及地理信息相关专业师生提供乡镇街道级别的矢量边界数据可用于空间分析、专题制图与城市精细化管理研究。资源包共14个文件约445KB以shp主文件为核心配套shx索引、dbf属性表、prj投影信息及sbx、sbn、cpg等辅助文件另含xml元数据与lock锁定文件构成完整可用的地理数据集。目前已有1262人学习下载说明其在同类数据中具备一定参考价值。通过该数据可获取各街道乡镇的精确边界与属性字段便于导入ArcGIS、QGIS等软件开展人口分布、交通流量或公共服务设施布局分析也可作为Web GIS在线地图服务的基础底图适合需要合肥本地精细化行政区划底图的中高级用户直接调用。1. 合肥市乡镇街道行政区划 shp 文件一份能直接进 GIS 的矢量底图做合肥本地项目的人多半遇到过这个场景领导丢过来一句“把合肥各乡镇的边界画出来”或者要按街道统计人口、门店、管网、地块结果手头只有一张模糊的截图或者一份属性表里写着“XX街道”却没有几何。这时候真正卡住你的不是分析模型而是缺一份能直接拖进 ArcGIS、QGIS 或 PostGIS 的合肥市乡镇街道行政区划 shp 文件。这份资源就是干这个的——它把合肥市下辖各区县的乡镇、街道级行政边界整理成标准 shp 矢量数据每个面要素带名称、编码等属性字段坐标系通常是地理坐标CGCS2000 或 WGS84可以直接用于制图、空间连接、裁剪和统计。适合做本地政务、规划、物流、商业选址、数据可视化的从业者也适合刚接触 GIS 矢量数据、想拿一份真实行政区划练手的新手。下面我按“拿到文件先看什么、怎么用、哪里会翻车”的顺序拆一遍。2. 拿到 shp 先别急着打开字段、坐标系与拓扑自检shp 不是单个文件而是一组同名文件捆在一起工作的。很多人第一次拿到压缩包解压后只看到.shp就双击结果要么打不开要么打开是空白问题基本都出在配套文件缺失或坐标系没对上。这一章先把“这份数据到底能不能用”判断清楚再谈后面的操作。2.1 shp 文件组的构成与最小可用集合一个完整的 shapefile 至少包含三个文件.shp存几何、.shx存索引、.dbf存属性表。少任何一个ArcGIS 会直接报错QGIS 可能能读几何但属性全空。实际分发时还常带.prj坐标系定义、.cpg字符编码、.sbn/.sbx空间索引。合肥这份乡镇街道数据你解压后应该能看到一整套同名文件文件名一致、只差扩展名。文件作用缺失后果.shp几何坐标无法显示图形.shx几何索引部分软件报错或读取异常.dbf属性字段有图无属性无法按名称筛选.prj坐标系定义坐标数值在但不知道是什么坐标系.cpg字符编码中文名称乱码判断能不能用的第一步确认这五个文件齐全且文件名完全一致。常见坑是解压工具把中文名截断导致.dbf和.shp名字对不上软件就认不出来。2.2 坐标系识别为什么你的图和底图对不上坐标系是 shp 使用里最容易翻车的一环。合肥市位于东经约 117°、北纬约 31.8°如果数据是地理坐标系CGCS2000 或 WGS84坐标值应该是117.xx, 31.xx这种小数如果看到的是500000, 3400000这类大数说明它是投影坐标系常见 CGCS2000 3 度带中央经线 117°E。在 QGIS 里查看方式右键图层 → 属性 → 信息看 CRS 一栏。在 ArcGIS 里图层右键 → 属性 → 源看“空间参考”。如果.prj缺失软件会提示“未知坐标系”此时图形能显示但一旦和在线底图叠加就会偏移几百米甚至几公里。提示不确定坐标系时先加载一份已知坐标系的底图如天地图或本地正射影像把 shp 叠上去看是否吻合。吻合说明坐标系一致整体平移说明是基准差异缩放不对说明是投影 vs 地理坐标混用。2.3 拓扑自检乡镇边界该不该严丝合缝行政区划面数据理论上应该无缝无重叠——相邻乡镇共享一条边界全市拼起来正好是合肥市域。但实际拿到的数据常有缝隙、重叠、悬挂节点。做统计汇总时重叠会导致同一块地被算两次缝隙会导致漏算。快速自检方法QGIS矢量 → 几何工具 → 检查有效性能查出无效几何再用“拓扑检查器”插件查重叠和缝隙。ArcGIS 用“拓扑”工具集建立“不能有重叠”“不能有缝隙”规则。如果只是做展示轻微缝隙不影响如果要做面积统计或空间连接建议先修复。修复思路用“融合”按区县字段合并再“拆分多部件”或者用“消除”工具处理小缝隙。这一步不做后面按乡镇统计人口或 POI 时数据对不上排查起来非常痛苦。3. 从 shp 到可用成果筛选、裁剪与空间连接实操数据自检通过后真正的活是把它变成你要的结果。这一章覆盖三个最高频操作按名称筛出目标乡镇、按范围裁剪出研究区、把业务数据挂到行政区上做统计。每个操作给可复现的步骤和参数说明。3.1 按属性筛选目标乡镇街道假设你只要合肥市包河区的乡镇街道。在 QGIS 里右键图层 → 筛选 → 输入表达式区县 包河区字段名要换成你数据里实际的区县字段可能是NAME、区县名称、COUNTY等。先打开属性表看一眼字段名再写别凭猜。ArcGIS 里图层右键 → 属性 → 定义查询 → 输入同样的 SQL 表达式。如果字段里存的是编码而不是名称比如340111代表包河区那就用编码筛选。筛选完建议“导出 → 保存选定要素”生成一个新的 shp避免后续操作误改原数据。参数说明表达式里的字符串用单引号字段名用双引号QGIS或不加引号ArcGIS 视版本而定。中文乱码时先确认.cpg文件内容是否为UTF-8不是就改成 UTF-8 再重新加载。3.2 按研究区范围裁剪如果你手上有一个更大的范围比如整个安徽省的乡镇数据只想保留合肥市部分用“裁剪”工具。QGIS矢量 → 地理处理工具 → 裁剪输入图层选安徽数据叠加图层选合肥市界输出即合肥范围。ArcGIS分析工具 → 提取 → 裁剪注意“XY 容差”保持默认勾选“使用输入要素裁剪几何”。# 用 geopandas 做同样的裁剪适合批量处理 import geopandas as gpd anhui gpd.read_file(anhui_towns.shp) hefei gpd.read_file(hefei_boundary.shp) # 统一坐标系避免裁剪结果错位 anhui anhui.to_crs(hefei.crs) hefei_towns gpd.clip(anhui, hefei) hefei_towns.to_file(hefei_towns_clip.shp, encodingutf-8)逻辑说明to_crs保证两个图层坐标系一致这是裁剪前必须做的一步否则结果可能为空或错位。gpd.clip按叠加图层的外框裁剪输入图层保留相交部分。encodingutf-8保证中文属性写出不乱码。参数说明如果裁剪后要素数量异常少检查两个图层的坐标系是否一致如果属性丢失检查字段名是否超过 dbf 的 10 字符限制shapefile 格式限制长字段名会被截断。3.3 空间连接把业务数据挂到乡镇街道上最常见的需求有一批门店点位或人口表格想知道每个乡镇街道各有多少。点位用“空间连接”表格用“按位置选择”或“连接”。QGIS 空间连接矢量 → 数据管理工具 → 按位置连接属性输入图层选乡镇面连接图层选点位几何谓词选“包含”输出即每个面带上落在其中的点属性。# geopandas 空间连接示例 import geopandas as gpd towns gpd.read_file(hefei_towns.shp) points gpd.read_file(stores.shp) # 坐标系必须一致 points points.to_crs(towns.crs) joined gpd.sjoin(points, towns, howleft, predicatewithin) # 按乡镇统计点位数量 count_by_town joined.groupby(乡镇名称).size().reset_index(name门店数) print(count_by_town.head())逻辑说明sjoin把每个点匹配到包含它的面predicatewithin表示点在面内。howleft保留所有点没匹配上的点乡镇字段为空方便排查。最后按乡镇名称分组计数。参数说明predicate可选within、intersects、contains点位统计一般用within。如果点位正好落在边界上within可能不匹配改用intersects更稳妥。统计前先检查有没有点位落在合肥市域外那部分会被归为空值。4. 避坑与排查shp 使用中最容易翻车的五件事这一章是我自己踩过和帮人排查过的真实问题按“现象 → 原因 → 解决”写。每一条都对应一个具体操作场景不是泛泛而谈。4.1 打开是空白或图形跑到地图外现象shp 加载后看不到图形或者图形出现在地图角落、和底图完全对不上。原因坐标系缺失或错误。.prj文件丢失时软件按默认坐标系处理导致位置偏移或者数据本身是投影坐标但被当成地理坐标显示。解决先看坐标数值范围。如果是117, 31级别是地理坐标如果是500000, 3400000级别是投影坐标。在 QGIS 里用“设置图层 CRS”手动指定正确坐标系再“导出 → 另存为”写入.prj。ArcGIS 用“定义投影”工具。指定后仍偏移说明基准不同如 WGS84 vs CGCS2000需要做基准转换。4.2 中文属性乱码现象属性表里乡镇名称显示成????或锟斤拷。原因.cpg文件缺失或编码声明与实际不符。shapefile 的 dbf 不强制声明编码软件只能猜。解决新建一个同名.cpg文件内容写UTF-8重新加载。如果还是乱码用 QGIS 的“重新加载编码”试 GBK。ArcGIS 里可以在目录窗口中右键 dbf → 属性 → 更改编码。批量处理时用 Python 读 dbf 指定 encoding 参数。4.3 面积统计结果偏大或偏小现象按乡镇统计面积结果和官方公布数据差很多。原因坐标系是地理坐标度直接算面积得到的是“平方度”没有物理意义。必须投影到等面积坐标系再算。解决QGIS 用“重投影图层”转到 CGCS2000 3 度带合肥用中央经线 117°EEPSG 视具体带号或用 Albers 等面积投影。ArcGIS 用“投影”工具。投影后再计算几何面积。# 投影后计算面积单位平方米 towns gpd.read_file(hefei_towns.shp) towns_proj towns.to_crs(EPSG:4547) # CGCS2000 3度带 117E towns_proj[area_m2] towns_proj.geometry.area towns_proj[area_km2] towns_proj[area_m2] / 1e6参数说明EPSG:4547是 CGCS2000 3 度带中央经线 117°E 的投影代码合肥适用。不同带号代码不同选错会导致面积变形。算完和官方数据对比偏差在合理范围内即可。4.4 空间连接后统计数量对不上现象明明有 1000 个点连接后按乡镇汇总只有 800 个。原因部分点落在市域外、落在边界上未被within匹配、或坐标系不一致导致匹配失败。解决先检查未匹配的点乡镇字段为空落在哪里用“按位置选择”看它们和市域的关系。边界上的点改用intersects。坐标系不一致的先统一。另外注意重复匹配一个点落在两个面重叠区会被算两次需要先修复面的重叠。4.5 导出 GeoJSON 或转其他格式后属性丢失现象shp 转 GeoJSON、转 Excel、转 CAD 后部分字段没了或名称被截断。原因shapefile 的 dbf 字段名限制 10 个字符超长会被截断转 GeoJSON 时字段类型不兼容如日期格式会丢失。解决导出前先把字段名改成 10 字符以内或改用 GeoPackage.gpkg作为中间格式它没有字段名长度限制。转 CAD 时注意 CAD 只认几何和少量属性复杂属性表建议单独导出表格再关联。5. 进阶把乡镇 shp 用进自动化流程与在线地图前面讲的都是单机操作这一章说两个进阶方向一是把这份数据接进自动化脚本批量出图或批量统计二是转成在线地图能用的格式做 Web 端展示。最后给一个我常用的验证习惯。5.1 批量出图按乡镇循环导出专题图做汇报时经常要给每个乡镇单独出一张图。手动一个个筛太慢用 Python 循环import geopandas as gpd import matplotlib.pyplot as plt towns gpd.read_file(hefei_towns.shp) towns_proj towns.to_crs(EPSG:4547) for name, group in towns_proj.groupby(乡镇名称): fig, ax plt.subplots(figsize(8, 8)) towns_proj.plot(axax, color#e8e8e8, edgecolorwhite) group.plot(axax, color#c0392b, edgecolorwhite) ax.set_title(name) ax.axis(off) plt.savefig(foutput/{name}.png, dpi150, bbox_inchestight) plt.close()逻辑说明先画全部乡镇作为底再高亮当前乡镇保证每张图的空间参照一致。groupby按名称分组循环导出。bbox_inchestight去掉多余白边。参数说明dpi控制清晰度汇报用 150 够印刷用 300。中文字体需要提前设置plt.rcParams[font.sans-serif] [SimHei]否则标题乱码。输出目录要先建好否则报错。5.2 转 GeoJSON 与在线地图展示Web 端展示常用 GeoJSON 或矢量瓦片。转 GeoJSON 最简单ogr2ogr -f GeoJSON hefei_towns.geojson hefei_towns.shp -lco ENCODINGUTF-8逻辑说明ogr2ogr是 GDAL 自带工具-f GeoJSON指定输出格式-lco ENCODINGUTF-8保证中文编码。转出来的 GeoJSON 可以直接被 Leaflet、Mapbox、OpenLayers 加载。参数说明如果文件很大乡镇级一般不大几百 KB 到几 MB可以用-simplify简化几何减少体积。在线展示时注意坐标系用 WGS84EPSG:4326GeoJSON 标准要求如此投影坐标需要先转回地理坐标。5.3 一个我每次都走的验证习惯不管拿到什么来源的行政区划 shp我固定走三步验证第一步看属性表字段和记录数和官方公布的乡镇街道数量对一遍数量对不上说明有缺失或重复第二步叠加已知底图看位置确认没有整体偏移第三步随便挑一个乡镇算面积和公开数据比偏差超过 10% 就回去查坐标系。这个习惯帮我拦下过好几次“看起来能用、一算就错”的数据。从那以后我每次拿到新的 shp 都强制走一遍这三步宁可多花十分钟也不在错误数据上做半天分析。希望帮到你。本文还有配套的精品资源点击获取
单位体系与能量-力一致性新手完全指南)
