1:1000地形图无人机航测技术设计书要点解析：从像控点到空三加密

简介《2021年航测技术设计书范本》是一份面向无人机航测项目技术设计环节的PDF参考文档适合测绘工程师、项目负责人及测绘相关专业学生用于编写1:1000比例尺地形图航测方案。文档以YS200无人机搭配Sony A7R航摄仪为例给出了飞行平台参数、航摄仪焦距与像素选择、基于高斯克吕格3°带投影的坐标系统设置并结合CORS控制网和350个像控点布设方案完整呈现了从测区踏勘、布标到航空摄影、空三加密、数字产品制作与验收的技术流程可作为实际项目技术设计书的编写模板。资源为单个PDF文件大小237KB内容涵盖任务来源、作业依据、技术设计、实施计划与质量控制等章节结构清晰便于直接参考或调整后使用。目前已有489人学习浏览对需要快速搭建航测技术设计框架的读者具有较强参考价值。1. 航测技术设计书为什么决定项目成败拿到这份 2021 年 3 月编制的 1:1000 地形图航测技术设计书范本时我特意先跳过了任务来源直接翻到技术设计那几节。项目覆盖巴克什营与虎什哈两个测区共 260 平方公里飞行平台是 YS200 固定翼无人机搭配 Sony A7R 全画幅相机。看完全文最直接的感觉是航测项目的质量往往不取决于飞机有多贵、相机有多少像素而取决于设计书里的参数是否闭环——从基准面高程到像控点判刺精度任何一环松动后端的空三加密和 DEM 生产都会加倍偿还。这份范本适合两类人正在准备 1:500 到 1:2000 比例尺项目技术设计的测绘工程师以及要给无人机航测配置数据采集方案的 GIS 开发者。它把很多只能靠经验判断的东西比如航线旁向重叠、像控点密度、GPS 摄站坐标获取变成了可执行的参数表。2. 测区踏勘与YS200/SonyA7R选型逻辑2.1 任务来源如何影响设计书走向范本里任务来源是典型的模板写法但真正有信息量的是“1:1000 地形图”这几个字。比例尺一旦定下来平面中误差、高程中误差、像控点布设密度、航线高度就全都被框住了。按 GB/T 15967 的要求1:1000 比例尺地形图的平面位置中误差平地、丘陵地不大于图上 0.3mm也就是实地 0.3m高程中误差按照地形类别分别限定。这就是为什么后续像控点联测结果要给出平面中误差 ±2.6cm、高程中误差 ±1.7cm 这样的量化指标而不是写“精度满足要求”了事。在写自己的设计书时我一般会把比例尺对应的限差直接做成表头放在技术设计章节的第一页这样评审专家和后续作业人员一眼就能看到验收标准不用来回翻规范。2.2 YS200飞行平台参数与单架次覆盖能力YS200 作为固定翼平台翼展 2100mm机长 850mm续航 75 分钟巡航速度 76km/h弹射起飞、伞降回收。这套参数决定了它适合做中等面积的正射数据采集尤其在丘陵地区固定翼的巡航效率比多旋翼高得多。计算单架次覆盖能力时可以用巡航速度乘以有效曝光时间估算。航向覆盖长度等于像幅沿飞行方向的尺寸乘以比例尺分母以 A7R 全画幅短边 24mm 沿航向为例地面覆盖约 24mm × H/f。基线距离由航向覆盖长度乘上重叠率余量决定。参数数值影响翼展2100mm抗风性与滑翔稳定性巡航速度76km/h单架次覆盖面积的上限续航时间75min分区长度与架次规划起飞/回收弹射/伞降对场地有要求不适合城市密集区我做架次规划时通常先按“75 分钟续航减 10 分钟冗余”来预估再根据测区形状把航线长度控制在 30km 以内避免逆风回程时电量紧张。弹射起飞和伞降回收对场地要求比较严格测区踏勘时必须确认起降场与测区边缘的距离固定翼爬升阶段不能直接进入航线。2.3 航摄仪选型Sony A7R与地面分辨率计算设计书选 Sony A7R 作为航摄仪3700 万像素、焦距 35mm。航测选相机不是只看像素数关键指标是像元尺寸、机械快门的稳定性以及镜头畸变是否经过标定。用地面分辨率公式 GSD H × pixel / f可以在外业前就确认航高是否满足比例尺要求。# 地面分辨率GSD与覆盖宽度计算 f 0.035 # 焦距单位m35mm pixel 0.0000049 # 像元尺寸单位m全画幅约7360x49123700万像素 h 500.0 # 相对航高单位m gsd h * pixel / f print(f相对航高{h}m时地面分辨率: {gsd*100:.2f} cm) # 全画幅传感器 36mm x 24mm计算地面覆盖 cover_x 0.036 * h / f # 旁向地面覆盖 cover_y 0.024 * h / f # 航向地面覆盖 print(f航向覆盖 {cover_y:.2f} m, 旁向覆盖 {cover_x:.2f} m)计算逻辑很直接把像元尺寸和焦距换算成同一单位再乘上相对航高。按上述参数500m 航高时 GSD 约 7cm满足 1:1000 地形图对影像分辨率的要求。注意这里用的是相对航高而不是绝对飞行高度基准面高程不同实际飞行高度要相应调整。覆盖宽度决定了航线的间距和每条航线的像片张数在设计书中必须写清楚。相机镜头最好做一次实验室检校把内方位元素和畸变参数交给空三软件否则边缘像元的畸变会影响连接点精度。2.4 作业依据规范清单如何约束技术设计作业依据列了 CH/Z 30012、CH/Z 3002、CH/Z 3003、CH/Z 3005、CH/Z 8021以及 GB/T 18314、GB 50026、GB/T 20257.1、GB/T 18316、GB/T 15967 等一批规范。这些规范看似冗余实际各自管一段CH/Z 3005 管航摄设计GB/T 15967 管比例尺地形图航空摄影规范GB/T 18316 管验收。写设计书时我会在每条依据后面标注被引用到的具体章节比如“GB/T 15967 第 5.2.3 条 基准面高程”这样外业查验时能快速定位依据也方便评审专家核对。3. 高斯克吕格3度带投影与CORS控制网设计3.1 坐标基准1980西安坐标系与1985国家高程基准的配对设计书里明确规定用 1980 西安坐标系和 1985 国家高程基准。这两个基准虽然在日常应用中被 2000 国家大地坐标系逐步替代但在存量项目和西部偏远测区仍然常见。1980 西安坐标系的椭球是 IAG-75高程基准的原点在青岛水准原点。外业时最容易翻车的是把 WGS84 的经纬高直接当成地方坐标用。RTK 采集到的原始坐标属于 CGCS2000 或 WGS84必须通过转换参数转到 1980 西安坐标系。转换参数通常由控制点计算获得至少需要 3 个以上公共点做布尔莎七参数转换这个步骤要在技术设计书里明确写出起算点来源和转换方法。3.2 高斯克吕格3度带投影参数计算投影方法选高斯克吕格 3° 带中央子午线 117°。3° 带的带号从东经 1.5° 开始每 3° 一带计算中央子午线的公式是带号 n int((L 1.5) / 3)中央子午线 L0 3n。对于 117°E 经线n 等于 39中央子午线正好是 117°不用跨带。# 计算3度带带号和中央子午线 lon 117.0 # 测区中心经度 zone int((lon 1.5) / 3) central_meridian zone * 3 print(f3度带带号: {zone}, 中央子午线: {central_meridian}°) # 用pyproj定义高斯克吕格投影1980西安坐标系的IAG-75椭球 from pyproj import CRS crs_xian80 CRS.from_proj4( projtmerc lat_00 lon_0117 k1 x_0500000 y_00 ellpsIAU76 unitsm ) print(crs_xian80)第一段代码给的是带号计算公式外业布点位时可以快速判断测区是否跨带。第二段代码定义了一个面向 117° 中央子午线的高斯克吕格投影对象x_0 是 500km 的加常数避免西移出现负坐标。实际项目里要用与当地测绘主管部门一致的椭球参数和转换参数不能照搬通用参数否则平面坐标差几米的案例并不少见。另外1:1000 地形图分幅采用 50cm×50cm 矩形分幅图名用自然村或单位名称命名保证唯一性这些在图廓整饰阶段就要准备好。3.3 CORS控制网与像控点联测中的已知点利用控制网设计采用 CORS 系统利用已知控制点作为起算点。CORS 的优势是省去静态观测缺点是它对通信链路要求高山区信号遮挡时会断链。像控点联测采用 GPS RTK 技术平面中误差 ±2.6cm、高程中误差 ±1.7cm这个精度数据说明测区环境里 CORS 信号质量较好。RTK 联测时的关键是每一站开机后必须先用至少两个已知点检核平面校正残差不超过 2cm、高程残差不超过 3cm 再开始采集。观测历元数不能省固定解状态下建议单个像控点采集 10 个历元以上测回间平面互差小于 2cm。范本里没写这些细节但在实际生产中像控点成果不合格往往不是中误差超标而是粗差没被剔除。4. 航摄分区、航线布设与像控点判刺实操4.1 基准面高程与航摄分区划分的计算航摄分区划分依据 CH/Z 3005 的 5.2.3 条核心要求是分区内地形高差不应大于 1/4 相对航高基准面高度取高点平均高程和低点平均高程的均值。这里的难点是“代表性地物”怎么选。我通常先在测区踏勘时用手持 GPS 或已有的 DEM 粗查高程变化按 1km×1km 网格统计高点与低点的平均值再用航线设计软件复算。如果高差超过限制就要把分区界线沿地形突变处破图廓划分否则同一分区内不同位置的地面分辨率波动太大。# 基准面高程计算与分区高差校验 h_high_avg 680.0 # 分区内高点平均高程单位m h_low_avg 510.0 # 分区内低点平均高程单位m rel_h 500.0 # 拟定相对航高单位m h_base (h_high_avg h_low_avg) / 2 if h_high_avg - h_low_avg rel_h / 4: print(f基准面高程 {h_base:.1f} m分区内高差在允许范围内) else: print(高差超限需重新划分航摄分区)这段代码做了两件事计算基准面高程以及校验分区高差是否小于等于相对航高的四分之一。如果校验结果超限就把分区标高重新拟定比如把相对航高从 500m 拉到 600m或者把分区分开分别计算各自的基准面。设计书里要写清楚相对航高与基准面高程的关系否则飞行操控员会按绝对高度飞导致重叠度漂移。4.2 航线布设方向与覆盖保证航线布设按东西向直线飞行常规航摄航线与图廓线平行特定情况下沿图廓中心线敷设。水域测区要避免像主点落水确保岛屿完整覆盖并构成正常重叠的立体像对。这是固定翼项目最容易出问题的地方。航线方向定成东西向主要是考虑太阳角度和侧风影响避免南北向飞行时太阳光从侧面照进相机产生光晕。实际布设航线时还要把航摄分区边界向外扩展航向覆盖超出摄区边界不少于 2 条基线旁向覆盖超出摄区边界一般不少于像幅的 50%最少不小于 35%。我用航线设计软件时会检查首末航线和测区边缘的距离如果不足就手动把范围外扩。重叠度直接影响空三的构网质量。航向重叠度一般设置为 70% 到 80%旁向重叠度设置在 30% 到 50% 之间。植被茂密或者城区高楼较多的测区建议航向重叠度取 80%旁向重叠度取 50%用冗余换连接点数量减少后期空三出现空洞的风险。重叠度的计算与航高、地面覆盖宽度联动通过代码可以快速算出航线间距。# 根据旁向重叠度计算航线间距 overlap_side 0.35 # 旁向重叠度 35% cover_x 514.29 # 旁向地面覆盖米由航高与焦距算出 spacing cover_x * (1 - overlap_side) print(f旁向重叠度 {overlap_side:.0%} 时航线间距 {spacing:.2f} m) # 根据航向重叠度计算曝光基线和张数 overlap_along 0.75 # 航向重叠度 75% cover_y 342.86 # 航向地面覆盖米 baseline cover_y * (1 - overlap_along) route_len 25000.0 # 单条航线长度米 photo_count route_len / baseline 1 print(f基线距离 {baseline:.1f} m单航线像片数 {photo_count:.0f} 张)参数说明spacing 是相邻航线之间的距离设计航线时直接按这个值导入飞控地面站baseline 是两次曝光之间的水平距离它决定同一航线上的像片张数。像片张数再加 2 条基线作为航向覆盖冗余就是最终每条航线的曝光数。实际操作中还要根据风速调整旁向间距侧风大会导致飞机偏离航线间距太小会出现旁向重叠不足。4.3 像片重叠度与补摄/重摄规则设计书要求航向和旁向重叠度根据项目情况设定同时规定绝对漏洞、相对漏洞及其他严重缺陷必须及时补摄。补摄必须按原设计航线飞行两端超出漏洞外两条基线。补摄作业中常见的坑是飞控位置重置后飞行员凭感觉“补飞一条线”结果补摄航带和原航带在高程基准上对不齐。补摄时要把原始航线的 POS 数据导出对比确保高度差在 10m 以内再接续。项目推荐值说明航向重叠度70%80%高差大的区域取上限旁向重叠度30%50%城区或林地取上限航向覆盖超出摄区边界≥2 条基线防止边缘构网失败旁向覆盖超出摄区边界≥35%50%一般取 50%补摄时机的选择也影响效率。一般建议当天飞行结束后立即用快拼影像检查重叠度发现漏洞第二天上午补飞等整个架次压一周后再补摄地面植被和光照可能已经变了空三时会出现色彩不一致。4.4 像控点布测密度、判刺目标与整饰本测区共布设像控点 350 个其中巴克什营 100 平方公里布设 150 个虎什哈 160 平方公里布设 200 个像控点全部采用白色反光带。布设密度约 0.74 点/平方公里对 1:1000 比例尺的丘陵测区来说属于中等偏密的水平。像控点加密时选择线状地物交角良好的交点或者影像尺寸小于 0.2mm 的点状地物中心点位判刺精度为图上 0.1mm。这个 0.2mm 和 0.1mm 的区别很容易被忽略0.2mm 是选点目标在影像上的最大尺寸限制0.1mm 是判刺时刺孔中心与真实目标位置的偏差上限。实际操作中选点要尽量避开高差突变的地方否则 RTK 测的是地面点而航片刺的是房顶或植被顶比高误差会进入空三计算。像控点整饰分为正反两面。正面以刺孔为中心画直径 7mm 的圆平高点用红色高程点用绿色点号与高程用分数式标注。背面用铅笔在 2cm×2cm 方框内画略图在方框旁边注明点位说明保证刺孔影像、桩位和略图说明一致。每次看到设计书里“并注明了点号选刺者、检查者的签名”这句我都会提醒外业同事签名不只是流程而是责任追溯。空三发现问题需要返工重刺时有签名能快速找到当事人一起判断是判刺偏差还是航线位移。5. 空三加密精度验证与数字产品质检技巧5.1 从空三到DEM作业流程中的关键节点设计书的作业流程图把全数字摄影测量分成数据采集、空三加密、数字产品制作、精度分析与验收几个阶段。实际生产中用 UAS Master 7.0 做空三、CASS 9.0 做数字化制图是这个项目比较典型的软件组合。空三加密的关键是连接点的数量和分布。单张影像至少要有 6 到 9 个连接点且不能集中在影像中心。我在检查空三成果时先看区域网平差後的像点残差中误差再看每个像控点的检查点残差是否超过 0.3m。像控点在空三里被分成两类参与平差的控制点和纯检查点。技术设计书里应该写明检查点比例不低于总数量的 20%这样精度分析才可信。5.2 用中误差公式验证像控点成果设计书给出像控点平面中误差 ±2.6cm、高程中误差 ±1.7cm。这个值在外业联测阶段就可以算出来方法与空三检查点精度评定的原理一致核心公式是中误差等于较差平方和的均值开根号。import math # 像控点RTK检核较差单位cm检测点实测值 - 已知值 dx [2.1, -3.2, 1.5, -0.8, 2.8, -1.2, 0.5, -2.1] dy [1.4, 2.2, -1.8, 1.9, -2.5, 0.7, 2.6, -1.1] dh [1.7, -2.2, 1.5, -0.9, -1.6, 2.0, -1.1, 1.8] n len(dx) mxy math.sqrt(sum(x*x y*y for x, y in zip(dx, dy)) / n) mh math.sqrt(sum(h*h for h in dh) / n) print(f平面中误差: ±{mxy:.2f} cm) print(f高程中误差: ±{mh:.2f} cm) if mxy 5.0 and mh 5.0: print(满足1:1000像控点限差要求) else: print(超限需要重测)这段代码的计算逻辑是先把每个检测点的平面差分解到 X 和 Y 方向分别求平方和后累计除以点数得到方差最后开根号。平面中误差同时包含 X 和 Y 两个方向的分量高程单独计算。代码里的 5cm 是 1:1000 地形图像控点在常规测区下的参考限差实际限差以技术设计书为准。如果超限先检查 RTK 解状态是否固定解再检查已知点坐标框架是否一致。5.3 三个容易被忽略的质检细节第一快拼影像不能只用来检查重叠度。我一般会在当天飞行结束后用原始影像生成测区低分辨率快拼图叠加航线轨迹和像控点文件检查刺点位置是否落在目标区域内。设计书里写“判刺精度为图上 0.1mm”但在快拼图上先粗检一遍能发现因为导航坐标偏移导致的系统性刺点偏差。第二补摄航带重投影后要重新做一次重叠度检查重点看补摄区域和原有航带在接边位置的重叠比例是否达到航向 70%。第三检查 DEM 生成的等高线看有没有沿航带方向出现条带状抬高或者凹陷如果有要么是空三区域网平差时高程系统没稳住要么是像控点在航带接边处的比高记录有误。以上这些环节不是只在验收时做一次而是在空三格网平差、DEM 编辑、成果提交三个阶段各做一遍。前两遍用自动脚本批量出图最后一遍人工重点检查接边和像控点附近区域。这样处理完的成果再进 CASS 9.0 做地形图编辑返修率会明显低于直接导出就画图的项目。本文还有配套的精品资源点击获取