简介这是一份面向土地复垦方案编制人员、国土空间规划与测绘技术人员及高校相关专业师生的专业课件围绕土地复垦方案图件的制作规范与制图要点展开可用于岗位培训、课程教学与方案编制自查。整套内容以18页PPT形式呈现先讲图面配置、文字说明与图号顺序等一般规定再逐类拆解标准分幅土地利用现状图、工程总平面布置图、矿产资源赋存分布及开采工艺流程图、破坏土地现状及预测分析图、预测剖面图、土地复垦规划图与工程单体设计图的制图要求如比例尺不小于1:10000、底图叠加与灰化处理、剖面线垂直地层走向并标注开采移动角等最后给出制图流程与图式图例。压缩包共1个pptx文件约222KB轻量便于下载后直接授课或自学。目前已有154人浏览学习适合作为编制复垦方案时的规范对照与培训讲义。1. 图件不是配图是复垦方案的证据链评审会上被退回的复垦方案相当一部分问题出在图件上损毁范围与现状地类图斑对不上复垦单元面积表和图上量算差了几个百分点图廓外缺坐标系与高程基准说明。文字报告写得再细这几处对不上方案就要返工。土地复垦方案图件制作规范PPT教案这个题目真正要处理两件事把分散在编制规程、制图规范和省级报件要求里的图件条款收敛成一套能照着做的作业流程再把流程拆成两三个课时讲得完、讲完就能上手的教案切片。适合两类人看刚接手复垦方案、出图总被打回的新手和要给团队做内训、手里只有零散条文加几份旧图的技术负责人。下面按图件清单、底层参数、成图流程、教案组织四段推进每一步都给到能直接抄的参数和命令。2. 复垦方案要交哪些图清单、规范依据与比例尺对照2.1 一套完整方案的图件构成复垦方案的图件按回答的问题分成三组。区位类回答项目在哪包括项目区位置图和交通位置图底图常用行政区划图叠加项目边界尺度可以粗但要能看出与周边城镇、道路、水系的相对关系。现状类回答现在是什么样包括土地利用现状图和土地损毁现状图前者按现行土地利用分类落到三级地类后者按损毁类型挖损、塌陷、压占、污染和损毁程度分色表达。复垦类回答以后做成什么样包括复垦规划图、复垦前后对比图、典型工程设计图与剖面图规模大的矿区还要加复垦单元划分图、适宜性评价图和监测点位布置图。判断一份方案图件齐不齐不看张数看每一张能不能独立回答一个问题。位置图看不出项目与县城的方位等于白画损毁现状图只有边界没有损毁类型评审时无法核对损毁面积的计算口径规划图有复垦方向却没有工程布局施工方拿不到可用的信息。2.2 规范依据怎么落到具体图件上常见做法是三层依据各管一段。复垦方案编制规程的系列文件负责图件种类、比例尺和必须表达的内容通则部分给出框图分矿种的部分再细化土地整治项目制图规范一类文件负责图式、图例、图廓整饰和注记规则省级自然资源主管部门的报件要求负责份数、装订方式、电子文件格式和坐标文件随附要求。三者出现冲突时优先级是报件要求 制图规范 行业规程因为审的是报件。提示把省级报件要求中关于图件的条目摘成一页纸贴在制图模板目录里。它比任何一本规程都更直接决定你的图能不能被收下。还有一个容易忽略的点图例必须与图上实际出现的符号一一对应图上没出现的地类不要在图例里留空位图例里有而图上没有的符号会被当成漏画。图例排列顺序按土地利用分类的层级排不要按画图先后排。2.3 比例尺、图幅与出图范围的对照比例尺不是随手定的它由最小上图图斑倒推。常用做法是图上最小图斑面积不低于 4 mm²换算到实地1:2000 时约 16 m²1:5000 时约 100 m²1:10000 时约 400 m²。小于这个面积的图斑要么归并到相邻图斑要么作为线状要素处理。图件名称常用比例尺成图图幅必备专题要素项目区位置图1:250001:50000A3 横版行政区界、项目区范围、交通、指北针土地利用现状图1:20001:10000A1/A0地类图斑、权属界线、地类符号、图斑面积表土地损毁现状图1:20001:5000A1损毁单元、损毁类型、损毁程度、现状地类叠加复垦规划图1:20001:5000A1复垦单元、复垦方向、工程布局、配套设施复垦前后对比图1:20001:5000A1前地类、后地类、对比统计表典型工程设计图1:2001:1000A3平面、剖面、尺寸标注、材料表项目区一旦超过一张图幅能覆盖的范围就要分幅出图并附接图表。接图表上标注每幅的图号与四角坐标图号命名建议用项目代号 图件代码 两位序号例如 FK-04-03这样批量出图时文件名天然有序。下面这份配置把图件清单写成代码后面驱动批量出图直接读它避免每次手改模板。# sheets.py —— 图件清单配置供批量出图脚本读取 MAP_SHEETS [ # code 图件代码, scale 比例尺分母, paper 出图纸型, must 必备专题要素 dict(code01, name项目区位置图, scale25000, paperA3, must[行政区界, 项目区范围, 交通]), dict(code02, name土地利用现状图, scale5000, paperA1, must[地类图斑, 权属界线, 地类符号]), dict(code03, name土地损毁现状图, scale5000, paperA1, must[损毁单元, 损毁类型, 损毁程度]), dict(code04, name复垦规划图, scale5000, paperA1, must[复垦单元, 复垦方向, 工程布局]), dict(code05, name复垦前后对比图, scale5000, paperA1, must[前地类, 后地类, 对比统计表]), dict(code06, name典型工程设计图, scale500, paperA3, must[平面, 剖面, 尺寸标注]), ] def paper_mm(paper): 返回纸型横版的宽高毫米 return {A3: (420, 297), A1: (841, 594), A0: (1189, 841)}[paper] for s in MAP_SHEETS: w, h paper_mm(s[paper]) # 图上宽度换算为实地宽度纸上米数 × 比例尺分母 cover_km w / 1000.0 * s[scale] / 1000.0 print(f{s[code]} {s[name]}: 1:{s[scale]} {s[paper]} 横版覆盖约 {cover_km:.2f} km)这段配置里scale决定成图精度和最小上图图斑paper决定一张图能装下多大的项目区A1 幅面 841 mm 宽在 1:5000 下只能覆盖约 4.2 km 宽项目区东西跨度超过这个数就必须分幅。must字段是给后面的自检环节用的出图前逐项对照缺一项就不允许导出。3. 底层参数先定死坐标系、分带与面积口径3.1 CGCS2000 与高斯-克吕格分带的选法平面坐标一律用 2000 国家大地坐标系高程用 1985 国家高程基准这两条写在图廓外的坐标系说明里不是可选项。投影用高斯-克吕格1:10000 及更小比例尺用 6 度带1:5000 及更大比例尺用 3 度带。3 度带的带号与中央经线关系是带号乘 338 带对应中央经线 114°E6 度带是带号乘 6 减 3。投影坐标有两种写法带带号的和不带带号的。带号版的 Y 值前面缀着带号例如 38 带写成 38500000 这样的形式不带带号的版本 Y 值从 500000 起算靠坐标系统定义本身来标识带号。麻烦往往出在这里——CAD 里画图时用不带带号导入 GIS 时按带带号解释两种数据叠加直接偏出去三百多公里。分带方式带号与中央经线CGCS2000 下的 EPSG 示例38 带适用比例尺3 度带带带号带号 n中央经线 3nEPSG:45261:5001:100003 度带不带带号CM 114°E无前缀Y 加 500 km 偏移EPSG:45471:5001:100006 度带带号 n中央经线 6n-3按实际带号查询1:100001:500003.2 面积量算的三套口径从一张图上能读出三种面积椭球面积大地坐标直接计算的曲面面积、投影平面面积高斯投影后的平面多边形面积、图上量算面积在成图上量取图斑得到的面积。三者数值不等投影平面面积相对椭球面积有变形离中央经线越远偏差越大且面积变形恒为正也就是投影面积一般略大于椭球面积。口径计算方式主要误差来源常用于椭球面积大地坐标直接积分边界采集精度面积汇总表、上报台账投影平面面积投影后多边形面积投影长度比变形GIS 图斑字段、与 CAD 对量图上量算面积成图上量取图纸伸缩、人工读数纸质图核验、现场比对同一份方案里损毁面积、复垦面积、耕地面积必须统一口径并且图上标注的图斑面积与汇总表、报告正文三处数字要完全一致。常见做法是图斑字段存投影平面面积汇总表另算一遍椭球面积两个数都留档报告里注明用的是哪一个。3.3 用 Python 自检坐标与面积出图前跑一遍下面的脚本坐标带号、面积口径和偏差量都能看住。import math from pyproj import Geod, Transformer from shapely.geometry import Polygon def gk_zone(lon, width3): 按经度返回高斯-克吕格带号与中央经线 if width 3: no round(lon / 3) return no, no * 3 no int(lon // 6) 1 return no, no * 6 - 3 print(gk_zone(113.4, 3)) # (38, 114) - 3 度带 38 带中央经线 114°E # CGCS2000 椭球参数与 GRS80 一致pyproj 可直接复用 geod Geod(ellpsGRS80) lons [113.4000, 113.4020, 113.4020, 113.4000] lats [34.1000, 34.1000, 34.1010, 34.1010] # 闭合环最后一个点与第一个点相同 lons_c, lats_c lons [lons[0]], lats [lats[0]] area_ell, peri geod.polygon_area_perimeter(lons_c, lats_c) print(f椭球面积 {abs(area_ell):.2f} m² 周长 {peri:.2f} m) # 投影到 CGCS2000 / 3 度带 38 带带带号再算平面面积 tr Transformer.from_crs(EPSG:4490, EPSG:4526, always_xyTrue) xs, ys tr.transform(lons_c, lats_c) area_plane Polygon(zip(xs, ys)).area print(f平面面积 {area_plane:.2f} m² 相对偏差 {(area_plane - abs(area_ell)) / abs(area_ell) * 100:.4f}%)脚本里gk_zone用来确认项目区落在哪个带、中央经线是多少与图廓外标注的坐标系统核对polygon_area_perimeter算椭球面积Transformer的两个参数分别是源坐标系 EPSG:4490CGCS2000 地理坐标和目标投影always_xyTrue保证按经纬度顺序传入。两位小数足够看出偏差量级同一个图斑两套口径的相对偏差通常在万分之几如果算出百分之几先查是不是坐标带号弄错了。3.4 参数环节的常见坑坐标 Y 值少一位、丢了带号前缀是最常见的一类导入 GIS 后整幅图跑到项目区以外把经纬度直接当平面坐标用画出来的图形完全对不上底图用 WGS84 采集的边界点没做坐标系转换就投影位置偏几十米跨两个 3 度带的项目区拼接时各自带了不同带号拼接后出现断裂。这些问题都能用gdalsrsinfo和ogrinfo -al -so快速确认出图前把每个图层的坐标系和图层范围打印一遍比事后返工便宜得多。4. 从底图到成图可复现的复垦图件制作流程4.1 图层组织与字段命名约定图层命名不统一是后期批量出图和改图最大的阻力。建议按业务对象 几何类型命名几何类型用点线面缩写字段名用中文或固定英文缩写二选一之后全项目不再变。字段里至少保留单元编号、地类代码、面积、复垦方向四项其中单元编号是图斑与报告表格之间的唯一关联键。图层名几何类型关键字段数据来源项目区范围面项目编号、面积批复文件附图、勘测定界权属界线线界线类型、权属单位不动产登记、村界现状地类图斑面地类代码、图斑面积年度国土变更调查损毁单元面损毁类型、损毁程度、面积实地调查、影像判读复垦单元面复垦方向、复垦面积、单元编号方案设计工程布局线/面工程类型、规格设计图4.2 图斑几何的硬规则与检查叠置关系上损毁单元之间、复垦单元之间不能有面积重叠也不能留下无归属的空隙复垦单元边界必须与权属界线严丝合缝不能让权属线从图斑中间穿过图斑共享边要完全重合不能出现肉眼看不见的细缝。这些在成图上不易察觉在面积汇总时立刻暴露。-- 检查复垦单元之间的重叠1 m² 以内视为容差 SELECT a.单元编号 AS 单元A, b.单元编号 AS 单元B, ST_Area(ST_Intersection(a.geom, b.geom)) AS 重叠面积 FROM 复垦单元 a JOIN 复垦单元 b ON ST_Intersects(a.geom, b.geom) AND a.单元编号 b.单元编号 WHERE ST_Area(ST_Intersection(a.geom, b.geom)) 1; -- 小于最小上图面积的图斑1:5000 下图上 4 mm² 约合实地 100 m² SELECT 单元编号, ST_Area(geom) AS 面积 FROM 复垦单元 WHERE ST_Area(geom) 100;第一条查询找重叠a.单元编号 b.单元编号避免同一个组合被统计两次容差 1 m² 是为了放过节点捕捉造成的浮点残差。第二条查询按最小上图面积筛图斑阈值随比例尺变1:2000 取 161:5000 取 1001:10000 取 400。筛出来的图斑要么归并要么在报告里单独说明保留理由。4.3 图例、注记与图廓整饰的参数化图例不要手工拼按符号库生成顺序与图上图层顺序一致。注记字号按比例尺换算图上字高不小于 1.8 mm1:5000 下对应实地 9 m 左右密集区可以降到 1.5 mm 但不能再小。图名用黑体地类与单元注记用宋体数字用等宽字体避免面积数字长短不一。图斑面积用表达式自动生成不手工敲。# ArcGIS Pro 计算字段表达式Python 3生成标准格式的面积注记 def fmt(area_m2, direction): 面积小于 1 hm² 时用 m²其余用 hm²保留四位小数 if area_m2 is None or direction is None: return None if area_m2 10000: return f{direction} {int(area_m2)} m² return f{direction} {area_m2 / 10000:.4f} hm² # 在计算字段工具中调用fmt(!Shape_Area!, !复垦方向!)Shape_Area是要素类的投影平面面积字段复垦方向取自属性表。写成函数而不是在表达式框里堆三元判断是为了让面积口径、单位切换规则和保留位数都能被审阅改规则时只改一处。图廓外必须齐备的要素包括图名、图号、数字比例尺与直线比例尺、指北针、图例、编制单位与日期、坐标系统与高程基准、分幅接图表缺一项在报件审查里都算形式不合格。4.4 批量出图与导出 PDF一张一张排版效率极低用图集Atlas按图幅网格批量出图是标准做法。下面的脚本在 QGIS 控制台里运行覆盖层一旦设定出图页数由网格要素数量决定。# QGIS 3.x 控制台按图幅网格批量导出 PDF from qgis.core import QgsProject, QgsLayoutExporter layout QgsProject.instance().layoutManager().layoutByName(复垦规划图) cov QgsProject.instance().mapLayersByName(图幅网格)[0] atlas layout.atlas() atlas.setCoverageLayer(cov) # 覆盖层决定出图页数与范围 atlas.setEnabled(True) atlas.setFilenameExpression(复垦规划图_ || \图幅号\) # 每页文件名 exp QgsLayoutExporter(layout) settings QgsLayoutExporter.PdfExportSettings() settings.dpi 300 # 报件与印刷都够用的分辨率 settings.exportMetadata False exp.exportToPdfs(layout, /data/out/复垦方案图件, settings)setCoverageLayer是图集的核心它把每一页的地图范围绑定到一个网格要素上setFilenameExpression用图幅号拼文件名导出的 PDF 顺序天然与分幅一致dpi设 300 是因为报件审查经常放大看注记150 会糊。若底图和矢量还没统一坐标系先用命令行统一避免在出图环节临场处理。# 栅格底图统一到 CGCS2000 / 3 度带 38 带带带号 gdalwarp -t_srs EPSG:4526 -tr 1 1 -r bilinear -co COMPRESSDEFLATE in.tif base_4526.tif # 矢量图层转坐标系并强制 UTF-8 属性防止注记乱码 ogr2ogr -f ESRI Shapefile -t_srs EPSG:4526 -lco ENCODINGUTF-8 \ -nln 复垦单元_4526 out/ 复垦单元.shp # 出图前核对成果的坐标系与图层范围 gdalsrsinfo -o proj4 base_4526.tif ogrinfo -al -so 复垦单元_4526.shp | head -n 20-tr 1 1把栅格重采样到 1 m 分辨率复垦图件在这个尺度上足够-r bilinear适合影像类底图若是分类栅格要换成-r near以保住类别值-lco ENCODINGUTF-8解决中文注记在部分软件里显示为乱码的问题。ogrinfo -al -so只看摘要不读要素几百兆的图层也能秒回适合放在出图流程的最后一道自检里。5. 把规范做成教案PPT 切片与出图前自检5.1 一页一图一参数的教案切片教案最忌讳按规范条文顺序讲。学员听完通则第几条之后还是不知道第一张图从哪下手。可行的组织方式是每页只讲一张图左侧放成图缩略图右上放这张图的必备要素清单右下放参数卡坐标系、比例尺、图幅、导出格式、命名规则。学员翻到哪一页就能照着把那张图做出来翻完六页等于拿到六张图的作业说明书。开场也不要用规范概述起手直接放一张被退回的图把批注圈出来——损毁范围压到了权属界线上、面积表与图斑差 3.2%、图廓外没有高程基准。让人先看到代价再讲规则。课时安排上一课时讲清单与参数一课时讲成图流程并演示批量出图第三课时留给演练。5.2 出图前自检清单与课堂演练自检清单是教案里最该被带走的产出物。课堂上让学员拿自己的项目跑一遍比讲十遍规范都有效。自检项检查方法判据坐标系与带号ogrinfo -al -so看图层范围X 值落在 38 带前缀范围内Y 值 6 位数加带号面积三处一致图斑字段 vs 汇总表 vs 报告正文三者完全一致口径在图中注明图斑重叠与空隙空间叠加查询重叠面积不超过 1 m²无归属空隙为零最小上图面积按比例尺阈值筛选1:5000 下不低于 100 m²图廓整饰要素对照要素清单逐项打勾八项齐全缺一不可图例与符号一致图例与图层清单双向核对无多余项、无缺漏项注记压盖成图缩放至 100% 目视关键注记不被符号压住导出文件打开 PDF 检查字体嵌入中文不变形线条不断线演练环节准备一套埋了错的样例数据一个图层丢了带号一个单元与相邻单元重叠了 500 m²一份汇总表用了椭球面积而图上字段是投影面积PDF 里图廓外少写了高程基准。限时二十分钟找出三处找出后当场用上面的命令和脚本验证一遍。教案的最后一页别放感谢页放这份自检清单的截图学员回去出图时第一件事就是对着它过一遍。本文还有配套的精品资源点击获取
简介这是一份面向土地复垦方案编制人员、国土空间规划与测绘技术人员及高校相关专业师生的专业课件围绕土地复垦方案图件的制作规范与制图要点展开可用于岗位培训、课程教学与方案编制自查。整套内容以18页PPT形式呈现先讲图面配置、文字说明与图号顺序等一般规定再逐类拆解标准分幅土地利用现状图、工程总平面布置图、矿产资源赋存分布及开采工艺流程图、破坏土地现状及预测分析图、预测剖面图、土地复垦规划图与工程单体设计图的制图要求如比例尺不小于1:10000、底图叠加与灰化处理、剖面线垂直地层走向并标注开采移动角等最后给出制图流程与图式图例。压缩包共1个pptx文件约222KB轻量便于下载后直接授课或自学。目前已有154人浏览学习适合作为编制复垦方案时的规范对照与培训讲义。1. 图件不是配图是复垦方案的证据链评审会上被退回的复垦方案相当一部分问题出在图件上损毁范围与现状地类图斑对不上复垦单元面积表和图上量算差了几个百分点图廓外缺坐标系与高程基准说明。文字报告写得再细这几处对不上方案就要返工。土地复垦方案图件制作规范PPT教案这个题目真正要处理两件事把分散在编制规程、制图规范和省级报件要求里的图件条款收敛成一套能照着做的作业流程再把流程拆成两三个课时讲得完、讲完就能上手的教案切片。适合两类人看刚接手复垦方案、出图总被打回的新手和要给团队做内训、手里只有零散条文加几份旧图的技术负责人。下面按图件清单、底层参数、成图流程、教案组织四段推进每一步都给到能直接抄的参数和命令。2. 复垦方案要交哪些图清单、规范依据与比例尺对照2.1 一套完整方案的图件构成复垦方案的图件按回答的问题分成三组。区位类回答项目在哪包括项目区位置图和交通位置图底图常用行政区划图叠加项目边界尺度可以粗但要能看出与周边城镇、道路、水系的相对关系。现状类回答现在是什么样包括土地利用现状图和土地损毁现状图前者按现行土地利用分类落到三级地类后者按损毁类型挖损、塌陷、压占、污染和损毁程度分色表达。复垦类回答以后做成什么样包括复垦规划图、复垦前后对比图、典型工程设计图与剖面图规模大的矿区还要加复垦单元划分图、适宜性评价图和监测点位布置图。判断一份方案图件齐不齐不看张数看每一张能不能独立回答一个问题。位置图看不出项目与县城的方位等于白画损毁现状图只有边界没有损毁类型评审时无法核对损毁面积的计算口径规划图有复垦方向却没有工程布局施工方拿不到可用的信息。2.2 规范依据怎么落到具体图件上常见做法是三层依据各管一段。复垦方案编制规程的系列文件负责图件种类、比例尺和必须表达的内容通则部分给出框图分矿种的部分再细化土地整治项目制图规范一类文件负责图式、图例、图廓整饰和注记规则省级自然资源主管部门的报件要求负责份数、装订方式、电子文件格式和坐标文件随附要求。三者出现冲突时优先级是报件要求 制图规范 行业规程因为审的是报件。提示把省级报件要求中关于图件的条目摘成一页纸贴在制图模板目录里。它比任何一本规程都更直接决定你的图能不能被收下。还有一个容易忽略的点图例必须与图上实际出现的符号一一对应图上没出现的地类不要在图例里留空位图例里有而图上没有的符号会被当成漏画。图例排列顺序按土地利用分类的层级排不要按画图先后排。2.3 比例尺、图幅与出图范围的对照比例尺不是随手定的它由最小上图图斑倒推。常用做法是图上最小图斑面积不低于 4 mm²换算到实地1:2000 时约 16 m²1:5000 时约 100 m²1:10000 时约 400 m²。小于这个面积的图斑要么归并到相邻图斑要么作为线状要素处理。图件名称常用比例尺成图图幅必备专题要素项目区位置图1:250001:50000A3 横版行政区界、项目区范围、交通、指北针土地利用现状图1:20001:10000A1/A0地类图斑、权属界线、地类符号、图斑面积表土地损毁现状图1:20001:5000A1损毁单元、损毁类型、损毁程度、现状地类叠加复垦规划图1:20001:5000A1复垦单元、复垦方向、工程布局、配套设施复垦前后对比图1:20001:5000A1前地类、后地类、对比统计表典型工程设计图1:2001:1000A3平面、剖面、尺寸标注、材料表项目区一旦超过一张图幅能覆盖的范围就要分幅出图并附接图表。接图表上标注每幅的图号与四角坐标图号命名建议用项目代号 图件代码 两位序号例如 FK-04-03这样批量出图时文件名天然有序。下面这份配置把图件清单写成代码后面驱动批量出图直接读它避免每次手改模板。# sheets.py —— 图件清单配置供批量出图脚本读取 MAP_SHEETS [ # code 图件代码, scale 比例尺分母, paper 出图纸型, must 必备专题要素 dict(code01, name项目区位置图, scale25000, paperA3, must[行政区界, 项目区范围, 交通]), dict(code02, name土地利用现状图, scale5000, paperA1, must[地类图斑, 权属界线, 地类符号]), dict(code03, name土地损毁现状图, scale5000, paperA1, must[损毁单元, 损毁类型, 损毁程度]), dict(code04, name复垦规划图, scale5000, paperA1, must[复垦单元, 复垦方向, 工程布局]), dict(code05, name复垦前后对比图, scale5000, paperA1, must[前地类, 后地类, 对比统计表]), dict(code06, name典型工程设计图, scale500, paperA3, must[平面, 剖面, 尺寸标注]), ] def paper_mm(paper): 返回纸型横版的宽高毫米 return {A3: (420, 297), A1: (841, 594), A0: (1189, 841)}[paper] for s in MAP_SHEETS: w, h paper_mm(s[paper]) # 图上宽度换算为实地宽度纸上米数 × 比例尺分母 cover_km w / 1000.0 * s[scale] / 1000.0 print(f{s[code]} {s[name]}: 1:{s[scale]} {s[paper]} 横版覆盖约 {cover_km:.2f} km)这段配置里scale决定成图精度和最小上图图斑paper决定一张图能装下多大的项目区A1 幅面 841 mm 宽在 1:5000 下只能覆盖约 4.2 km 宽项目区东西跨度超过这个数就必须分幅。must字段是给后面的自检环节用的出图前逐项对照缺一项就不允许导出。3. 底层参数先定死坐标系、分带与面积口径3.1 CGCS2000 与高斯-克吕格分带的选法平面坐标一律用 2000 国家大地坐标系高程用 1985 国家高程基准这两条写在图廓外的坐标系说明里不是可选项。投影用高斯-克吕格1:10000 及更小比例尺用 6 度带1:5000 及更大比例尺用 3 度带。3 度带的带号与中央经线关系是带号乘 338 带对应中央经线 114°E6 度带是带号乘 6 减 3。投影坐标有两种写法带带号的和不带带号的。带号版的 Y 值前面缀着带号例如 38 带写成 38500000 这样的形式不带带号的版本 Y 值从 500000 起算靠坐标系统定义本身来标识带号。麻烦往往出在这里——CAD 里画图时用不带带号导入 GIS 时按带带号解释两种数据叠加直接偏出去三百多公里。分带方式带号与中央经线CGCS2000 下的 EPSG 示例38 带适用比例尺3 度带带带号带号 n中央经线 3nEPSG:45261:5001:100003 度带不带带号CM 114°E无前缀Y 加 500 km 偏移EPSG:45471:5001:100006 度带带号 n中央经线 6n-3按实际带号查询1:100001:500003.2 面积量算的三套口径从一张图上能读出三种面积椭球面积大地坐标直接计算的曲面面积、投影平面面积高斯投影后的平面多边形面积、图上量算面积在成图上量取图斑得到的面积。三者数值不等投影平面面积相对椭球面积有变形离中央经线越远偏差越大且面积变形恒为正也就是投影面积一般略大于椭球面积。口径计算方式主要误差来源常用于椭球面积大地坐标直接积分边界采集精度面积汇总表、上报台账投影平面面积投影后多边形面积投影长度比变形GIS 图斑字段、与 CAD 对量图上量算面积成图上量取图纸伸缩、人工读数纸质图核验、现场比对同一份方案里损毁面积、复垦面积、耕地面积必须统一口径并且图上标注的图斑面积与汇总表、报告正文三处数字要完全一致。常见做法是图斑字段存投影平面面积汇总表另算一遍椭球面积两个数都留档报告里注明用的是哪一个。3.3 用 Python 自检坐标与面积出图前跑一遍下面的脚本坐标带号、面积口径和偏差量都能看住。import math from pyproj import Geod, Transformer from shapely.geometry import Polygon def gk_zone(lon, width3): 按经度返回高斯-克吕格带号与中央经线 if width 3: no round(lon / 3) return no, no * 3 no int(lon // 6) 1 return no, no * 6 - 3 print(gk_zone(113.4, 3)) # (38, 114) - 3 度带 38 带中央经线 114°E # CGCS2000 椭球参数与 GRS80 一致pyproj 可直接复用 geod Geod(ellpsGRS80) lons [113.4000, 113.4020, 113.4020, 113.4000] lats [34.1000, 34.1000, 34.1010, 34.1010] # 闭合环最后一个点与第一个点相同 lons_c, lats_c lons [lons[0]], lats [lats[0]] area_ell, peri geod.polygon_area_perimeter(lons_c, lats_c) print(f椭球面积 {abs(area_ell):.2f} m² 周长 {peri:.2f} m) # 投影到 CGCS2000 / 3 度带 38 带带带号再算平面面积 tr Transformer.from_crs(EPSG:4490, EPSG:4526, always_xyTrue) xs, ys tr.transform(lons_c, lats_c) area_plane Polygon(zip(xs, ys)).area print(f平面面积 {area_plane:.2f} m² 相对偏差 {(area_plane - abs(area_ell)) / abs(area_ell) * 100:.4f}%)脚本里gk_zone用来确认项目区落在哪个带、中央经线是多少与图廓外标注的坐标系统核对polygon_area_perimeter算椭球面积Transformer的两个参数分别是源坐标系 EPSG:4490CGCS2000 地理坐标和目标投影always_xyTrue保证按经纬度顺序传入。两位小数足够看出偏差量级同一个图斑两套口径的相对偏差通常在万分之几如果算出百分之几先查是不是坐标带号弄错了。3.4 参数环节的常见坑坐标 Y 值少一位、丢了带号前缀是最常见的一类导入 GIS 后整幅图跑到项目区以外把经纬度直接当平面坐标用画出来的图形完全对不上底图用 WGS84 采集的边界点没做坐标系转换就投影位置偏几十米跨两个 3 度带的项目区拼接时各自带了不同带号拼接后出现断裂。这些问题都能用gdalsrsinfo和ogrinfo -al -so快速确认出图前把每个图层的坐标系和图层范围打印一遍比事后返工便宜得多。4. 从底图到成图可复现的复垦图件制作流程4.1 图层组织与字段命名约定图层命名不统一是后期批量出图和改图最大的阻力。建议按业务对象 几何类型命名几何类型用点线面缩写字段名用中文或固定英文缩写二选一之后全项目不再变。字段里至少保留单元编号、地类代码、面积、复垦方向四项其中单元编号是图斑与报告表格之间的唯一关联键。图层名几何类型关键字段数据来源项目区范围面项目编号、面积批复文件附图、勘测定界权属界线线界线类型、权属单位不动产登记、村界现状地类图斑面地类代码、图斑面积年度国土变更调查损毁单元面损毁类型、损毁程度、面积实地调查、影像判读复垦单元面复垦方向、复垦面积、单元编号方案设计工程布局线/面工程类型、规格设计图4.2 图斑几何的硬规则与检查叠置关系上损毁单元之间、复垦单元之间不能有面积重叠也不能留下无归属的空隙复垦单元边界必须与权属界线严丝合缝不能让权属线从图斑中间穿过图斑共享边要完全重合不能出现肉眼看不见的细缝。这些在成图上不易察觉在面积汇总时立刻暴露。-- 检查复垦单元之间的重叠1 m² 以内视为容差 SELECT a.单元编号 AS 单元A, b.单元编号 AS 单元B, ST_Area(ST_Intersection(a.geom, b.geom)) AS 重叠面积 FROM 复垦单元 a JOIN 复垦单元 b ON ST_Intersects(a.geom, b.geom) AND a.单元编号 b.单元编号 WHERE ST_Area(ST_Intersection(a.geom, b.geom)) 1; -- 小于最小上图面积的图斑1:5000 下图上 4 mm² 约合实地 100 m² SELECT 单元编号, ST_Area(geom) AS 面积 FROM 复垦单元 WHERE ST_Area(geom) 100;第一条查询找重叠a.单元编号 b.单元编号避免同一个组合被统计两次容差 1 m² 是为了放过节点捕捉造成的浮点残差。第二条查询按最小上图面积筛图斑阈值随比例尺变1:2000 取 161:5000 取 1001:10000 取 400。筛出来的图斑要么归并要么在报告里单独说明保留理由。4.3 图例、注记与图廓整饰的参数化图例不要手工拼按符号库生成顺序与图上图层顺序一致。注记字号按比例尺换算图上字高不小于 1.8 mm1:5000 下对应实地 9 m 左右密集区可以降到 1.5 mm 但不能再小。图名用黑体地类与单元注记用宋体数字用等宽字体避免面积数字长短不一。图斑面积用表达式自动生成不手工敲。# ArcGIS Pro 计算字段表达式Python 3生成标准格式的面积注记 def fmt(area_m2, direction): 面积小于 1 hm² 时用 m²其余用 hm²保留四位小数 if area_m2 is None or direction is None: return None if area_m2 10000: return f{direction} {int(area_m2)} m² return f{direction} {area_m2 / 10000:.4f} hm² # 在计算字段工具中调用fmt(!Shape_Area!, !复垦方向!)Shape_Area是要素类的投影平面面积字段复垦方向取自属性表。写成函数而不是在表达式框里堆三元判断是为了让面积口径、单位切换规则和保留位数都能被审阅改规则时只改一处。图廓外必须齐备的要素包括图名、图号、数字比例尺与直线比例尺、指北针、图例、编制单位与日期、坐标系统与高程基准、分幅接图表缺一项在报件审查里都算形式不合格。4.4 批量出图与导出 PDF一张一张排版效率极低用图集Atlas按图幅网格批量出图是标准做法。下面的脚本在 QGIS 控制台里运行覆盖层一旦设定出图页数由网格要素数量决定。# QGIS 3.x 控制台按图幅网格批量导出 PDF from qgis.core import QgsProject, QgsLayoutExporter layout QgsProject.instance().layoutManager().layoutByName(复垦规划图) cov QgsProject.instance().mapLayersByName(图幅网格)[0] atlas layout.atlas() atlas.setCoverageLayer(cov) # 覆盖层决定出图页数与范围 atlas.setEnabled(True) atlas.setFilenameExpression(复垦规划图_ || \图幅号\) # 每页文件名 exp QgsLayoutExporter(layout) settings QgsLayoutExporter.PdfExportSettings() settings.dpi 300 # 报件与印刷都够用的分辨率 settings.exportMetadata False exp.exportToPdfs(layout, /data/out/复垦方案图件, settings)setCoverageLayer是图集的核心它把每一页的地图范围绑定到一个网格要素上setFilenameExpression用图幅号拼文件名导出的 PDF 顺序天然与分幅一致dpi设 300 是因为报件审查经常放大看注记150 会糊。若底图和矢量还没统一坐标系先用命令行统一避免在出图环节临场处理。# 栅格底图统一到 CGCS2000 / 3 度带 38 带带带号 gdalwarp -t_srs EPSG:4526 -tr 1 1 -r bilinear -co COMPRESSDEFLATE in.tif base_4526.tif # 矢量图层转坐标系并强制 UTF-8 属性防止注记乱码 ogr2ogr -f ESRI Shapefile -t_srs EPSG:4526 -lco ENCODINGUTF-8 \ -nln 复垦单元_4526 out/ 复垦单元.shp # 出图前核对成果的坐标系与图层范围 gdalsrsinfo -o proj4 base_4526.tif ogrinfo -al -so 复垦单元_4526.shp | head -n 20-tr 1 1把栅格重采样到 1 m 分辨率复垦图件在这个尺度上足够-r bilinear适合影像类底图若是分类栅格要换成-r near以保住类别值-lco ENCODINGUTF-8解决中文注记在部分软件里显示为乱码的问题。ogrinfo -al -so只看摘要不读要素几百兆的图层也能秒回适合放在出图流程的最后一道自检里。5. 把规范做成教案PPT 切片与出图前自检5.1 一页一图一参数的教案切片教案最忌讳按规范条文顺序讲。学员听完通则第几条之后还是不知道第一张图从哪下手。可行的组织方式是每页只讲一张图左侧放成图缩略图右上放这张图的必备要素清单右下放参数卡坐标系、比例尺、图幅、导出格式、命名规则。学员翻到哪一页就能照着把那张图做出来翻完六页等于拿到六张图的作业说明书。开场也不要用规范概述起手直接放一张被退回的图把批注圈出来——损毁范围压到了权属界线上、面积表与图斑差 3.2%、图廓外没有高程基准。让人先看到代价再讲规则。课时安排上一课时讲清单与参数一课时讲成图流程并演示批量出图第三课时留给演练。5.2 出图前自检清单与课堂演练自检清单是教案里最该被带走的产出物。课堂上让学员拿自己的项目跑一遍比讲十遍规范都有效。自检项检查方法判据坐标系与带号ogrinfo -al -so看图层范围X 值落在 38 带前缀范围内Y 值 6 位数加带号面积三处一致图斑字段 vs 汇总表 vs 报告正文三者完全一致口径在图中注明图斑重叠与空隙空间叠加查询重叠面积不超过 1 m²无归属空隙为零最小上图面积按比例尺阈值筛选1:5000 下不低于 100 m²图廓整饰要素对照要素清单逐项打勾八项齐全缺一不可图例与符号一致图例与图层清单双向核对无多余项、无缺漏项注记压盖成图缩放至 100% 目视关键注记不被符号压住导出文件打开 PDF 检查字体嵌入中文不变形线条不断线演练环节准备一套埋了错的样例数据一个图层丢了带号一个单元与相邻单元重叠了 500 m²一份汇总表用了椭球面积而图上字段是投影面积PDF 里图廓外少写了高程基准。限时二十分钟找出三处找出后当场用上面的命令和脚本验证一遍。教案的最后一页别放感谢页放这份自检清单的截图学员回去出图时第一件事就是对着它过一遍。本文还有配套的精品资源点击获取
