简介本资源为基于GS(2023)2767审图号制作的GIS地理信息数据集面向从事地图制图、空间分析与城市规划的GIS从业者及地理信息相关专业师生。包内共75个文件以shp、shx、dbf、prj等Shapefile核心格式为主配套sbn、sbx空间索引与cpg编码文件另含xml元数据及一份mxd工程文件压缩包约67.62MB。数据涵盖国界、省面、省界线、海岸、珊瑚礁、地名及0-15KM、15-30KM缓冲国界等图层并附南海附图框与南海诸岛小框可直接用于边界缓冲分析、区域制图与空间叠加等场景。借助prj投影信息与元数据读者能快速完成坐标校核与数据集成mxd工程文件则保留了图层符号化与版面布局便于直接复用或二次编辑。目前已有264人学习下载适合需要规范审图号底图与缓冲分析素材的用户参考使用。1. GS(2023)2767审图号GIS数据从合规底图到业务落地的完整路径如果你手里拿到一份标注了“GS(2023)2767”审图号的GIS数据第一反应大概率是这数据能用吗能直接进项目吗边界在哪我最初接触这类数据时也踩过坑——把审图号当成普通数据版本号结果在坐标系和属性字段上翻车。GS(2023)2767审图号本质上是自然资源主管部门对地图内容合规性的一个审核编号它意味着这份GIS数据的空间范围、要素表达、注记方式已经通过了官方审查。但“通过审查”不等于“拿来即用”它解决的是合规性问题不解决你的业务适配问题。这份数据通常以矢量图层形式提供包含行政区划、交通路网、水系、居民地等基础地理要素适合做区域规划、选址分析、专题图制作等场景。如果你需要一份能过审、能对外发布、能作为项目底图的数据源这类带审图号的GIS数据是绕不开的起点。但怎么用、怎么改、怎么保证改完还不违规才是真正花时间的地方。2. 审图号GIS数据的坐标系与属性结构先看懂再动手2.1 审图号背后的数据规格约束带审图号的GIS数据不是随便导出的Shapefile或GeoJSON。它的生产流程决定了几个硬约束坐标系通常采用2000国家大地坐标系CGCS2000高程基准为1985国家高程基准投影方式多为高斯-克吕格投影分带方式按3度带或6度带。你拿到数据后第一件事不是加载看效果而是确认坐标系定义是否完整。我见过不少案例数据本身坐标值是对的但.prj文件缺失或写错导致在QGIS里叠加其他图层时偏移几百米。审图号数据还有一个特点属性字段往往经过规范化处理字段名可能是中文拼音缩写或标准代码比如“XZQDM”代表行政区代码“DLMC”代表地类名称。这些字段在后续空间分析中直接决定你能不能做聚合统计。常见做法是先用ogrinfo或QGIS的图层属性面板查看元数据。如果字段名是代码需要对照《基础地理信息要素分类与代码》标准来映射。不要试图改字段名来“方便自己”改完再导出时如果忘了还原送审会被打回。2.2 用GDAL快速检查数据完整性的命令拿到数据后我一般会跑一套检查命令确认坐标系、范围、字段和几何类型。以下命令在GDAL 3.x环境下执行# 查看数据基本信息包括坐标系、要素数量、范围 ogrinfo -so -al GS20232767.gdb # 如果数据是Shapefile检查.prj文件是否存在且内容正确 ogrinfo -so GS20232767.shp | grep -A 5 Coordinate System # 导出坐标系定义到文本便于核对EPSG代码 gdalsrsinfo -o epsg GS20232767.prj # 检查几何有效性找出自相交或空几何 ogr2ogr -f ESRI Shapefile /dev/null GS20232767.shp -dialect sqlite -sql SELECT * FROM GS20232767 WHERE ST_IsValid(geometry) 0第一行ogrinfo -so -al输出图层摘要重点看Extent是否在你预期的地理范围内Feature Count是否和文档一致。第二行检查坐标系如果输出为空或报错说明.prj缺失。第三行gdalsrsinfo把投影信息转成EPSG代码CGCS2000对应的地理坐标系EPSG是4490投影坐标系常见4547到4554。第四行用SQL过滤无效几何如果返回结果不为空说明有要素需要修复。参数上-dialect sqlite启用SQLite方言才能用ST_IsValid函数/dev/null在Windows下换成NUL。2.3 属性字段映射与业务字段扩展审图号数据的属性表通常“够用但不富裕”。比如行政区划图层可能只有代码和名称没有面积、人口这些业务字段。你需要做的是在保留原始字段的前提下扩展。我习惯用GeoPandas做字段映射和扩展因为可以链式操作且容易回滚import geopandas as gpd import pandas as pd # 读取审图号数据指定编码避免中文乱码 gdf gpd.read_file(GS20232767.gdb, layerXZQ, encodingutf-8) # 查看原始字段 print(gdf.columns.tolist()) # 建立代码到名称的映射字典示例实际按标准文档 code_map { 110000: 北京市, 120000: 天津市, # ... 其他省级代码 } # 新增业务字段不覆盖原始字段 gdf[PROVINCE_NAME] gdf[XZQDM].map(code_map) # 计算面积注意投影坐标系下单位是米地理坐标系下是度 gdf_proj gdf.to_crs(epsg4547) # 转为CGCS2000 3度带投影 gdf[AREA_SQKM] gdf_proj.geometry.area / 1e6 # 导出时保留原始字段和新增字段 gdf.to_file(GS20232767_extended.gpkg, driverGPKG, encodingutf-8)这段代码的关键点encodingutf-8在读取File GDB时可能不生效因为GDB内部用UTF-16但导出为GPKG时必须指定UTF-8。to_crs(epsg4547)把地理坐标转成投影坐标4547对应CGCS2000 3度带第114度中央经线具体用哪个带号取决于你的数据范围。面积计算前一定要转投影否则算出来是平方度没有业务意义。导出为GPKG而不是Shapefile因为Shapefile字段名限制10个字符PROVINCE_NAME会被截断。注意扩展字段后如果数据要重新送审新增字段的内容不能与地图表达冲突。比如你加了一个“规划用地性质”字段但地图上没表达审查时可能被要求补充图例或删除字段。3. 从审图号数据到业务底图裁剪、投影与符号化3.1 按行政区边界裁剪的两种策略审图号数据往往是全省或全市范围你的项目可能只关心一个区。裁剪是最常见的操作但有两种做法按矩形范围裁剪和按行政区边界裁剪。矩形裁剪快但会保留边界外的要素碎片按行政区裁剪干净但需要额外的边界图层。我一般用GeoPandas的clip函数配合行政区边界import geopandas as gpd # 读取审图号数据和目标行政区边界 source gpd.read_file(GS20232767.gdb, layerDLTB) target_boundary gpd.read_file(target_district.shp) # 确保两者坐标系一致 if source.crs ! target_boundary.crs: target_boundary target_boundary.to_crs(source.crs) # 按边界裁剪 clipped gpd.clip(source, target_boundary) # 检查裁剪结果 print(f原始要素数: {len(source)}, 裁剪后: {len(clipped)}) print(f裁剪后范围: {clipped.total_bounds}) clipped.to_file(clipped_output.gpkg, driverGPKG)gpd.clip内部用的是空间索引加几何相交判断比手动循环快得多。但要注意如果目标边界本身有拓扑错误比如自相交裁剪结果可能丢失要素。裁剪前用target_boundary.is_valid.all()检查一下。另一个坑是裁剪后的要素属性表会保留原始字段但如果你只需要部分字段可以在裁剪后做列筛选减少数据量。3.2 投影转换的参数选择与验证审图号数据可能以地理坐标系EPSG:4490提供但做面积统计、距离分析时必须转投影。选哪个投影带原则是中央经线尽量靠近数据范围中心分带方式按数据跨度选。如果数据跨多个3度带要么分带处理要么用Albers等面积投影。我一般用以下代码做批量转换和验证import geopandas as gpd from pyproj import CRS # 定义目标投影CGCS2000 3度带中央经线117度对应EPSG:4548 target_crs CRS.from_epsg(4548) gdf gpd.read_file(GS20232767.gdb, layerXZQ) gdf_proj gdf.to_crs(target_crs) # 验证转换前后面积变化投影后面积应更合理 area_geo gdf.geometry.area.sum() area_proj gdf_proj.geometry.area.sum() print(f地理坐标系面积和: {area_geo:.6f} 平方度) print(f投影坐标系面积和: {area_proj/1e6:.2f} 平方公里) # 检查转换后坐标值是否在合理范围投影坐标通常是6-7位数 print(fX范围: {gdf_proj.total_bounds[0]:.0f} - {gdf_proj.total_bounds[2]:.0f}) print(fY范围: {gdf_proj.total_bounds[1]:.0f} - {gdf_proj.total_bounds[3]:.0f})EPSG:4548对应中央经线117度适用于中国东部地区。如果你不确定用哪个带号可以用pyproj的CRS.from_dict自定义中央经线。验证环节很重要投影后X坐标应该是6位数带号500公里偏移Y坐标7位数。如果X是8位数说明带号没去掉或者加了假东偏移。面积对比可以快速判断投影是否合适——如果投影后面积和地理坐标系面积差异过大说明中央经线偏离太远。3.3 符号化与审图号合规表达的边界符号化是审图号数据落地时最容易出问题的地方。审图号审查的是地图内容包括符号、颜色、注记。你可以在QGIS或ArcGIS里调整样式但有几条红线不能改变行政境界的等级表达比如把省界画成市界不能删除审图号注记不能改变重要地物的符号类型比如把高速公路画成普通公路。我一般会在QGIS里做样式然后导出为SLD或QML文件方便复用!-- QGIS样式文件片段行政区填充 -- qgis renderer-v2 typesingleSymbol symbols symbol typefill name0 layer classSimpleFill prop kcolor v255,255,255,255/ prop koutline_color v0,0,0,255/ prop koutline_width v0.26/ /layer /symbol /symbols /renderer-v2 /qgis这个SLD片段定义了白色填充、黑色边框的行政区样式。实际项目中我会把审图号注记单独放在一个图层确保导出图片时不会被其他图层遮挡。如果项目要求输出PDF或图片记得在布局里加上审图号文本位置通常在右下角或左下角字号不小于6pt。提示修改样式后如果数据要重新发布建议保留一份原始样式的备份。审图号审查时可能要求提供样式说明原始样式是最安全的参照。4. 审图号GIS数据落地避坑5个血泪教训4.1 坐标系“看起来对”但叠加偏移现象把审图号数据和其他来源的底图叠加发现整体偏移几十到几百米但坐标系定义显示都是CGCS2000。原因审图号数据可能用的是“CGCS2000地理坐标系”但实际坐标值是“投影坐标值”或者.prj文件写的是EPSG:4490但数据本身是EPSG:4548。这种“定义与值不匹配”的情况在数据流转中很常见。解决用gdalsrsinfo导出实际坐标系再用ogrinfo查看坐标值范围。如果X值是6-7位数说明是投影坐标需要把.prj改成对应的投影EPSG。如果无法确定用已知控制点做仿射变换校正。4.2 字段名截断导致属性丢失现象导出Shapefile后原本的“XZQDM”字段变成“XZQDM”但“PROVINCE_NAME”变成“PROVINCE_”数据内容被截断。原因Shapefile格式限制字段名最多10个字符超出部分自动截断。审图号数据原始字段可能刚好10个字符以内但你扩展的字段超了。解决中间过程用GPKG或File GDB只在最终交付时转Shapefile并提前把字段名缩到10字符以内。或者直接用GPKG交付现在主流GIS平台都支持。4.3 裁剪后面积统计对不上现象按行政区裁剪后统计各地类面积总和比原始数据小或者某些地类消失。原因裁剪时边界图层有缝隙或重叠导致部分区域被重复裁剪或漏裁。另一个原因是裁剪后几何被切割但属性表中的面积字段没有重新计算。解决裁剪前用ST_Union把边界图层合并成一个要素消除内部缝隙。裁剪后用geometry.area重新计算面积字段不要依赖原始面积字段。4.4 审图号注记被遮挡或丢失现象出图时审图号文本被其他图层盖住或者导出PDF后审图号不见了。原因审图号注记可能放在某个图层的标注里导出时被其他图层覆盖。或者布局中的文本元素没有锁定位置页面尺寸变化后跑出画布。解决把审图号注记单独放在一个点图层设置标注优先级最高。在QGIS布局中把审图号文本放在最上层并锁定位置和大小。导出前用PDF预览检查。4.5 数据更新后审图号失效现象基于审图号数据做了编辑比如新增道路重新出图时被告知审图号不适用。原因审图号对应的是特定版本的数据内容。任何对地图内容的修改包括新增、删除、修改要素都可能导致审图号失效。解决如果只是样式调整不影响内容审图号通常仍有效。如果涉及内容修改需要重新送审。我的习惯是把审图号数据作为“底图参考”业务数据单独成层出图时叠加。这样底图不变审图号持续有效。5. 审图号数据的进阶用法版本管理与自动化检查5.1 用Git LFS管理审图号数据版本审图号数据文件大直接放Git仓库会爆。我一般用Git LFS跟踪.gdb和.gpkg文件配合.gitattributes# 初始化Git LFS git lfs install # 跟踪GIS数据格式 git lfs track *.gdb/** git lfs track *.gpkg git lfs track *.shp git lfs track *.prj # 提交配置 git add .gitattributes git commit -m 配置Git LFS跟踪GIS数据这样每次数据更新都有版本记录出问题时可以回滚到上一个审图号对应的版本。注意.gdb是文件夹git lfs track *.gdb/**要写成这样才能跟踪内部文件。5.2 自动化合规检查脚本每次数据编辑后跑一遍合规检查避免送审被打回。以下脚本检查坐标系、字段完整性、几何有效性和审图号注记import geopandas as gpd import os def compliance_check(gdb_path, layer_name, expected_crs4490): issues [] # 检查文件是否存在 if not os.path.exists(gdb_path): issues.append(数据文件不存在) return issues # 读取数据 gdf gpd.read_file(gdb_path, layerlayer_name) # 检查坐标系 if gdf.crs is None: issues.append(坐标系未定义) elif gdf.crs.to_epsg() ! expected_crs: issues.append(f坐标系为{gdf.crs.to_epsg()}预期{expected_crs}) # 检查几何有效性 invalid_count (~gdf.geometry.is_valid).sum() if invalid_count 0: issues.append(f存在{invalid_count}个无效几何) # 检查空几何 empty_count gdf.geometry.is_empty.sum() if empty_count 0: issues.append(f存在{empty_count}个空几何) # 检查必填字段按实际标准调整 required_fields [XZQDM, XZQMC] for field in required_fields: if field not in gdf.columns: issues.append(f缺少必填字段: {field}) return issues # 执行检查 problems compliance_check(GS20232767.gdb, XZQ) if problems: print(发现合规问题:) for p in problems: print(f - {p}) else: print(合规检查通过)这个脚本可以集成到CI流程里每次数据提交自动跑。expected_crs参数按实际审图号要求设置CGCS2000地理坐标系是4490投影坐标系按带号调整。必填字段列表根据《基础地理信息要素分类与代码》和项目要求维护。5.3 审图号数据的长期维护习惯我带过的项目里审图号数据出问题最多的情况不是技术问题而是管理问题。我的习惯是每份审图号数据单独建目录目录名包含审图号、日期和版本比如GS20232767_20230601_v1。目录里放三个文件原始数据、编辑后的数据、变更说明。变更说明用Markdown写记录改了什么、为什么改、谁改的。这样半年后回头看能快速定位问题。另一个习惯是审图号数据永远保留一份只读副本所有编辑在副本上进行。只读副本用chmod 444锁定防止误操作。希望帮到你。本文还有配套的精品资源点击获取