ARTICLE DETAIL

资讯详情

深耕郑州网站建设与运营推广的一线实战洞察。

INPHO航测数据处理全流程解析:从空三平差到DOM生产

INPHO航测数据处理全流程解析:从空三平差到DOM生产 不用先把软件全家桶讲一遍再谈流程。直接说结论INPHO这套航测系统强项从头到尾都是“平差”和“匀色”这两件事。现在很多航测软件把空三做成一个黑盒点一下跑完就出报告但遇到大面积测区、山区或者高精度要求黑盒模型经常会给你挖坑。INPHO的数据处理流程之所以在测绘生产单位里用了这么多年还不淘汰核心就在于每个环节都能动手干预、能审计、能追溯——这才是航测数据处理的正常姿势。这篇文章就以我实际跑项目的顺序把INPHO从空三前准备到最终成果输出的整个流程拆开讲力求每个步骤都留下可复现的配置和判断逻辑。1. 先搞懂INPHO全家桶的分工再谈流程1.1 每个模块解决哪个环节的问题INPHO不是单体型软件是一整套航测解决方案。刚开始接触的人往往被ApplicationsMaster这个总控界面搞懵以为打开它就能跑所有流程。实际上ApplicationsMaster更像一个“调度台工程管理器”它把下面这些模块串在一起跑批处理任务模块核心职责对应流程节点ApplicationsMaster工程项目管理、任务流调度贯穿全流程Match-AT自动连接点提取、像点量测空三加密前期inBLOCK/inBAL区域网构建与光束法平差解算空三加密核心BINGOGPS辅助光束法区域网平差空三加密核心OrthoMaster单片正射纠正DOM生产OrthoVista镶嵌、匀色、分幅输出DOM生产DTMasterLiDAR点云分类、滤波、编辑点云数据处理SCOP不规则三角网构建、DEM内插高程模型生产这套分工逻辑和Pix4D或者CC那种“一键到底”的思路完全不同。INPHO是把测图生产的整个工艺拆成一个一个工段每一步生成中间成果你可以检查每一步的质量再进入下一步。从生产管理角度这反而是优点——出了问题定位快到底是匹配错了还是平差散了一目了然。1.2 INPHO与其他航测软件的流程差异在哪里如果只看官方宣传页会觉得INPHO好像只是多了一个“看起来很专业”的界面。真正把它用起来之后才能感受到差异。首先INPHO的自动空三不是只听你一句“开始匹配”就埋头干。Match-AT在提取连接点时有一个很关键的策略逻辑它会把影像分成不同金字塔层级从低分辨率到高分辨率逐层匹配每一层都在前一层结果的基础上加细细化像点坐标。这种“由粗到精”的做法意味着再弱的纹理区域也有可能被匹配到代价是运算次数多。相比之下很多全自动软件为了速度在低纹理区域直接放弃匹配最后造成区域网穿洞或者变形你得后期手动补点反而耽误时间。其次BINGO平差模块是INPHO的灵魂所在。它把GPS/IMU观测值、像控点坐标、像点量测值放在同一个平差模型里统一解算能处理带有系统误差的航飞数据。这个能力在用老飞机或者低成本POS做惯性导航数据时特别重要那些数据里往往带着漂移误差粗暴地把POS坐标当成固定真值用模型必变形。BINGO是真正“辅助”平差给POS观测值一个合理的权重而不是把它当成绝对约束。第三OrthoVista的匀色算法是我目前用过的软件里最稳定的。它做全局色彩平衡的时候是把整个测区的影像作为一个整体进行统计学优化不是简单地对单张影像做直方图拉伸也不会出现相邻影像色差一眼就能看出来的情况。这一点在后端交给甲方验收时太重要了——DOM漂亮愿意挑毛病的人就少一半。1.3 一条典型的INPHO生产链路长什么样给你一份我实际项目里跑的流程清单方便脑子里先搭起框架创建工程配置相机参数、坐标系统、投影参数导入影像、POS数据或GPS/IMU原始数据、像控点Match-AT提取连接点进行区域网构建和光束法平差迭代剔除粗差查看空三精度报告若不合格返回修正生成DSM或利用LiDAR点云生成DEM用OrthoMaster逐片生成正射影像用OrthoVista镶嵌匀色分幅输出DOM成果质检、格式转换、交付。后面所有章节都会围绕这条链路展开把每一环容易踩的坑单独拉出来说。2. 空三之前这些准备细节决定后面会不会返工2.1 相机参数和畸变模型的“坑”打开INPHO新建工程第一步不是导入影像而是配置相机参数。很多人图快直接从相机文件里读一个默认值就完事这是我最不推荐的做法。相机参数包括主距焦距对应的像素值、主点偏移x0y0、像元尺寸、径向畸变系数K1K2K3、切向畸变系数P1P2。注意INPHO和其他软件对畸变模型的符号约定不完全一样。比如Pix4D和ContextCapture里有自己的畸变参数体系直接迁移容易出错。建议把原始相机检校文件拿过来按照INPHO的参数定义逐一核对后再填入。一个非常实用的小技巧在导入影像时先在ApplicationsMaster里用“影像预览”刷新几张原片检查EXIF信息里有没有内方位元素。如果相机检校文件里有主点偏移值但EXIF里写的是零那就要以检校文件为准。这类问题不查出来空三跑完精度报告上控制点残差的分布会呈现一种奇怪的系统性偏移排查起来很麻烦。2.2 坐标系统和投影参数错了真的全白跑坐标系统设定是INPHO项目设置里最容易被忽略的地方也是返工率最高的地方。原因很简单航测原始数据里坐标经常是混合的。举个例子你从飞机上拿到的POS数据可能是WGS84经纬度或UTM投影但测区控制点是2000国家大地坐标系高斯投影下的平面坐标。INPHO允许你在工程里同时设定水平基准和垂直基准这一点相当灵活但灵活的另一面是容易混。我见过不止一次有人把POS数据坐标系选成了和像控点一样的坐标系结果投影参数差异导致空三解算直接发散或者精度报告看起来还行、实际成果叠到地形图上整体偏移几十米。所以我在新建工程时有一个强制习惯先做一个坐标系统清单把“影像POS原始坐标系”“控制点成果坐标系”“成果输出坐标系”三列写清楚再对照着设置。如果涉及地方独立坐标系注意核实中央子午线、投影面大地高、加常数这些参数。INPHO里提供了很多标准的坐标系统模板地方坐标系一般需要手动定义。2.3 控制点和检查点安排的门道控制点在INPHO空三中的作用不只是提供绝对位置约束还参与整个平差模型的形变控制。控制点分布不均匀——比如集中在测区中央四周没有——会导致区域网边缘精度不可控而且这类问题在精度报告里经常看不出来因为平差只会优化掉已知点附近的误差。按规范做就好测区四周要有控制点中间也要有控制点如果你用全野外布点按图幅间距布设如果是像控点加密至少保证每个测区外扩的点位覆盖。这里说一个我在实际生产中反复验证的经验检查点不要只放在平缓地区一定要有几个放在地形起伏较大的地方。因为航测平差对高程的约束和平面不完全一样起伏区域的高程误差更容易暴露模型变形。控制点坐标的格式需要注意。INPHO里导入控制点时接受文本格式的点名、XYZ或东、北、高程。Z值的高程基准要和DEM测区内一致如果控制点高程是正常高似大地水准面而POS高程是大地高你必须先做高程转换。不要指望在INPHO里一键转换航测生产里高程基准不一致是致命错误。2.4 POS数据导入前要做的清理工作POS数据导入INPHO之前我建议先做一次“预体检”。具体包括看格式常见的有SBET二进制、ASCII文本、以及飞控软件导出的自定义格式。INPHO虽然能识别多种格式但最好在Excel或者Python里先看一遍表头和数据预览。看时间标签POS数据每一行有时间戳影像也有采集时间。两者的时间系统必须一致GPS周秒、UTC或当地时间。如果时间基准不对影像外方位元素的初始值就全错了空三根本收敛不过去。看轨迹连续性把POS轨迹画出来看看有没有明显的跳变点。飞机转弯盘旋时的POS数据经常质量差导入前如果知道哪些航段数据有问题可以先标记出来在平差时给予低权重甚至剔除。这套“预体检”工作看起来不起眼实际上能省下后面排错的一个通宵。3. 自动空三的硬核过程Match-AT连接点与BINGO平差3.1 连接点提取从策略到参数连接点提取是所有后续操作的地基。Match-AT在提取连接点时不是单纯地在每一对相邻影像上找同名点而是构建影像间的拓扑关系图——它得知道哪些影像之间有重叠、重叠度多少然后在这个基础上做匹配并且会刻意提取“多度重叠点”也就是同一地面点出现在3张及以上影像上的像点。这类点多区域网的强度就高尤其是航线之间的连接靠的就是这些多度重叠点。实际设置时有几个参数值得认真调匹配策略/金字塔层级默认设置通常能从4层金字塔开始逐层加密。如果你的测区纹理强城区、农田可以适当减少高层级匹配直接在中层级开始如果测区是大面积水面、雪地、沙漠这类弱纹理区域就老老实实开着完整金字塔让低分辨率层先兜住全局。每张影像目标点数INPHO默认可能是几百到上千个点。我个人建议不要一上来就把数量拉太高先按默认或稍高的值跑一遍看平差结果再决定。点数太多会让粗差剔除变得困难反而影响收敛速度。多片点最低重合度建议至少要求每点被3张影像覆盖。如果测区航线重叠度低可能达不到这时要检查航飞设计——航线旁向重叠度低于30%时连接点的几何强度会严重不足这是后期无法修复的。3.2 区域网构建与平差方案选择连接点提取完接下来是构建区域网。这里的核心是把每个连接点、每条航线的GPS/IMU观测值、所有像控点组织成一个整体的误差方程然后解算。INPHO里通常用inBLOCK或inBAL来完成区域网构建与自由网平差。自由网平差的目的是在不引入控制点的情况下先看整个区域网的内部几何一致性如何——影像之间的相对关系是否严丝合缝。如果这一步通过再引入控制点和GPS/IMU辅助数据做绝对定向。这里要特别强调GPS辅助平差的意义。现在无人机航测基本都有RTK/PPK很多人觉得POS数据已经足够准了控制点放几个意思一下就行。实际上如果你把POS数据权重设得太高平差会把影像位置牢牢钉死在POS给出的坐标上像控点反而纠正不了系统偏差。正确做法是在BINGO里给POS观测值设定一个合理的先验中误差让它在平差中作为带权观测值参与修正而不是作为固定值。一套我常用的权重参考经验如果使用经过PPK解算的高精度POS其平面位置精度可标定为0.05m高程0.1m控制点坐标精度标定为0.02m这样赋予权重后才能得到整体最佳点位精度。当然这些数字要结合你的飞行平台和大地测量手段调整。3.3 粗差剔除与平差迭代的节奏把握不管匹配算法多先进连接点里总会有误匹配。BINGO平差的经典处理流程是“迭代-剔除-再平差”先进行第一次自由网平差得到像点残差设定残差阈值比如像点残差大于1.5~2个像素的点标记为粗差自动剔除粗差后重新平差重复上述过程直到没有新的粗差产生确认自由网精度合格后引入控制点和POS辅助数据做绝对平差。有些人希望一步到位直接加控制点平差发现精度不好就去手动改控制点这是很典型的错误思路。粗差问题必须在自由网阶段解决干净控制点阶段只负责绝对定向。粗差阈值怎么定一个像素还是两个像素取决于你的影像分辨率和地面分辨率。通常无人机低空影像地面分辨率在2~5cm之间允许像点残差1个像素相当于地面2~5cm偏差如果用大幅面航摄仪地面分辨率10cm以上一个像素就是10cm。所以不要盲目把阈值调到0.5个像素那会把正常点全部误杀。3.4 空三精度报告怎么看到位INPHO的精度报告是整个数据处理流程中最关键的输出物。你需要重点看这些指标像点RMSRMS of image coordinates所有连接点像点残差的均方根值一般要求小于1/3到1/2个像素。过大了说明匹配质量差或者相机参数有问题。控制点残差Control point residuals每个控制点在XYZ方向上的残差。要注意看最大残差和RMS残差前者能反映个别点位是否有粗差。检查点残差Check point residuals同样看XYZ向RMS。检查点是你的验收依据它的精度直接等于你最终成果的定位精度参考。这里特别提醒精度报告里控制点残差很好不代表检查点残差也好。如果控制点过多平差会把点位误差强推到控制点上检查点残差才会暴露真实问题。所以我在每个项目里都会专门保留一批点位不参与平差只作为检查点参与精度评价。4. DSM、DEM与3D点云数据处理流程4.1 从空三成果到密集点云空三平差结束后你已经有了每张影像的精确外方位元素。下一步是生成密集点云。INPHO在这一点上和主流三维建模软件思路类似——利用多视影像匹配生成高密度三维点云整个处理也是逐级匹配的策略通过多次迭代获取亚像素级的视差值然后通过前方交会生成三维坐标。IMPORTANT生成点云前一定要确认空三精度达标否则就是在垃圾数据上加工。点云的密度受地面分辨率和匹配质量双重影响地面分辨率高、纹理丰富点云密度就高。同样的数据如果你跑完发现点云有明显的“条带效应”带状高低错落大概率是空三没有完全消除航线系统误差回头检查平差比恼火地滤波更高效。4.2 点云滤波与地面点分类DTMaster扮演的角色如果项目需要生产DEM那么点云滤波环节就绕不开DTMaster。DTMaster的核心功能是对LiDAR点云进行交互式分类自动滤波分出地面点、植被点、建筑物点、噪点。你要把它理解成一个数据清洗车间而不是完全自动的黑盒。自动滤波算法对地形有一定假设——地面是一个连续、平滑变化的表面。在地形陡峭或植被密集的地方算法经常把真实地面点误分类为植被点或者在建筑物边缘把地面点错误保留。这时候就需要人工介入用剖面视图逐条航线查看分类结果。实际操作时我习惯先做一次自动分类然后在DTMaster里按“地面”和“非地面”两类显示。重点关注山脊线、陡坎、建筑边缘这些特征点沿等高线方向切剖面检查地面点是否贴合真实地形。如果有大面积错误可以先用“局部地面滤波”重跑绝大多数遗留问题都能通过局部滤波解决。这里插一句即使你的目的是做DOM而不是DEMDSM的质量也决定了正射纠正的质量。城市高层建筑的DSM如果出现大面积空洞纠正出来的建筑物就会变形。所以别把点云环节单独当成“做地形”的步骤它与DOM质量直接相关。4.3 从分类点云到DEM的生成逻辑点云分类完成后SCOP负责把地面点内插成DEM。重点说一下内插算法的选择。INPHO的SCOP提供多种内插算法常用的是线性预测least squares prediction和移动曲面拟合。线性预测本质上是带空间相关性的最小二乘拟合它对数据点分布要求不高能生成比较平滑的地表模型移动曲面拟合则更适合地形变化剧烈的区域能保留细节但容易出现过度拟合。实际选择标准是这样一句话总结的地形起伏大就选能保细节的算法地形平缓就选平滑的算法。比如城市里要保留道路和建筑的边界在山地里要注意避免把山谷填平。内插分辨率的选择要与最终成果对应——DEM栅格尺寸不能大于点云密度所对应的空间分辨率否则就是人为拉低成果精度。用2cm分辨率的点云去生成10m分辨率的DEM那是浪费反过来用稀疏点云生成0.5m分辨率的DEM中间全是空洞和插值假象。4.4 结合点云数据处理的完整闭环现在热词里高频出现“3d点云数据处理流程”其实在INPHO这套体系里点云数据的处理闭环基本如下原始点云导入影像匹配点云或LiDAR点云噪点剔除低点、孤立点点云分类地面点、非地面点地面点生成DEM全部点生成DSMDSM用于正射纠正DEM用于等高线、断面等测绘成果质检剖面检查、高程精度比对。这套流程在INPHO里主要由DTMaster和SCOP完成但它生成的DEM又可以反向辅助正射纠正形成数据闭环——这也是INPHO区别于其他只做影像流程的软件的地方。5. DOM生产正射纠正与镶嵌匀色的实战细节5.1 OrthoMaster单片正射纠正的正确打开方式OrthoMaster的原理不复杂利用DSM/DEM提供的高程信息把中心投影影像逐像素纠正为正射投影。关键在两点一是纠正时使用的高程模型二是纠正过程中的隐藏地物处理。如果你用DSM纠正城市里的每栋楼、每棵树都会按它的真实高度纠正位置效果最好。但是DSM如果有错误比如建筑物顶部纹理匹配错位纠正出的DOM上建筑物会出现拉花变形。如果直接用DEM纠正建筑物的顶点会因为高程不足而偏移纠正出来的建筑像被推倒了一样。所以实际生产中城市区域优先用DSM纠正平缓丘陵地区用DEM就够了。如果DSM质量不好宁可在DSM编辑阶段花时间修也不要用错误DSM产出劣质DOM。纠正时还需要为每张影像生成“纠正掩膜”OrthoMaster会把影像中哪些区域是有效的、哪些区域因为遮挡或高程异常需要从相邻影像补充标记出来。掩膜生成结果直接决定后续镶嵌时接缝处是否会有重影。5.2 OrthoVista镶嵌与匀色最容易拉开差距的环节OrthoVista是我心中INPHO最被低估的模块。它干三件事镶嵌线自动生成、色彩均衡、成果分幅输出。镶嵌线的质量决定DOM最后有没有“接缝感”。OrthoVista默认会自动搜索相邻影像之间的最优接缝位置它会优先让接缝绕过明显地物比如房屋、道路标线避免从地物中间切开。即使自动生成的结果也应该人工检查几处明显区域比如主干路、桥梁、大型广场这些地方接缝一旦穿帮整个图幅都要重新修。匀色参数上OrthoVista提供全局色彩平衡和局部匀色两种模式。我做城市DOM项目时倾向于先用全局平衡然后用局部匀色校正残留色差。全局平衡会把整个测区所有影像的亮度、色相映射到统一标准上防止出现“一块亮一块暗”。局部匀色是对相邻影像间的过度区域做处理但它有个副作用如果处理过度会出现色块晕染或者纹理模糊。我的经验是局部匀色的强度控制在刚好能消除接缝色差的范围内就好不要把整条接缝边缘的影像色彩都重算一遍。一句话总结——匀色是为了让色差不明显不是为了彻底重塑影像色彩。5.3 分幅输出与格式选择OrthoVista产出成果时支持按标准图幅或自定义规则分幅切割。常用的输出格式是GeoTIFF务必设置好压缩方式和色彩深度。一个容易出问题的地方TIFF内部带的坐标参考信息和直方图信息可能会和实际不符交付前最好用GIS软件将影像与矢量套合检查确认影像范围、坐标系和像素分辨率都与要求一致。如果是做真正射影像TrueDOM把DSM作为高程参考进行纠正时注意输出成果中建筑物底部的“空洞”区域——真正射会把高层建筑遮蔽的底部区域暴露出来这些区域在原始影像上是看不到的需要从相邻影像补充信息。这需要OrthoMaster性能与OrthoVista填补功能的配合真正常用的项目反而少但用到的项目都对它要求极高。6. 质检环节精度评定与返工判断6.1 空三质量与成果精度的一致性很多项目直到最后一步才做精度检查这是最被动的。我建议在空三阶段就把检查点精度指标与最终成果精度要求做一次对照评估。比如某个测图项目要求DOM平面中误差不超过0.5m而你空三检查点的平面中误差已经做到了0.3m同时考虑到正射纠正、镶嵌等环节还会引入少量误差那么这个空三精度是够用的。如果空三检查点平面中误差已经在0.45m左右那DOM最终成果大概率超限这时候就要回到空三去查原因而不是直接继续生产DOM。6.2 DOM成果检查的实操方法检查DOM不用打开专业遥感软件在常见的GIS平台里就能完成。我常用的两步检查法第一步把DOM和空三检查点叠加显示检查点位的影像坐标和实测坐标是否吻合。这一步主要查绝对定位精度。第二步沿航带接缝方向连续浏览DOM查接缝处地物是否错位、颜色是否突变。检查的重点是道路标线、田埂、房顶边缘这类线性特征。一个看起来简单却很有效的方法把DOM按RGB三通道单独显示如果某一通道在接缝处出现阶梯状跳变说明原始影像色差或纠正参数还有问题。6.3 哪些情况必须返工哪些可以修补结合这些年踩过的坑我把常见问题分为“必须返工”和“可修补”两类问题类型判断标准处理建议空三不收敛控制点/检查点残差远超规范必须返回空三重跑点云大面积空洞建筑、植被区域点云缺失严重必须重新生成点云局部接缝错位偏差在1~2个像素内可局部重新镶嵌局部色差仅亮度差异几何无错位OrthoVista局部匀色小范围噪点单个异常点或小簇DTMaster局部去噪控制点残差在边缘系统性偏移检查点残差整体偏大返回空三调整控制点权重这条判断表是我每次项目质检时的实际操作依据。记住一个原则精度不合格的成果修补不如重跑。因为所有修补都是在错误基础上打补丁只有源头纠正才能保证后续链条不再积累误差。7. 一点经验最容易被新手忽略的三件事最后写几句实操感受。INPHO这套流程我前前后后跑了三年如果说有什么经验教训值得单独拿出来讲就是以下三件最容易被忽略却影响最大的事第一坐标系规划永远最先做。INPHO里一切参数都可以调唯独坐标系统错了后面所有成果都无法挽回。我的做法是在项目文件夹里建一个“项目配置.txt”把POS坐标系、控制点坐标系、成果坐标系、投影参数、高程基准逐项写清楚每次开工前先读一遍。第二空三阶段多花一小时后期少花一个工作日。很多人急着出DOM空三精度报告大致看一眼就往下走。实际上空三阶段发现和修复问题的成本最低到了镶嵌阶段发现问题返工的就是全场数据。第三点云处理不是“影像处理的附属品”。DEM的质量直接决定正射纠正精度而正射纠正的质量直接决定DOM能否通过验收。把DSM/DEM当成一等公民对待你的DOM成果会明显比别人稳定。这套流程跑通之后你会发现INPHO不是最省事的航测软件但它是一套让你心里有底的系统——每一个环节的数据质量都能被量化检查每一个工序都能定点回溯。如果项目对成果精度要求高、测区情况复杂这套“看得见、管得着”的处理流程才是真正能打硬仗的工作方式。
返回列表