测绘与空间地理信息
GEOMATICS&SPATIALINFORMATIONTECHNOLOGY
Vol.42ꎬNo.6Jun.ꎬ2019
基于机载GPS-IMU系统实现直接定向的研究
蒋勤玉1ꎬ丁兆连2
(1.广州南方卫星导航仪器有限公司ꎬ广东广州 510663ꎻ2.自然资源部第二航测遥感院ꎬ黑龙江哈尔滨150081)
摘要:主要研究利用高精度的机载GPS-IMU集成系统ꎬ还原航摄仪在摄影时像片的外方位元素ꎬ实现直接定
向ꎮ在航空摄影测量作业时ꎬ通过GPS差分定位和高精度IMU的内业数据处理、检校计算ꎬ获取到了位置及姿态信息ꎬ此方法简化了作业流程ꎬ节省了外业工作量ꎬ缩短了时间周期ꎮ文本对此过程的误差来源进行分析ꎬ提出优化方法ꎬ特别是通过检校计算补偿系统误差ꎬ对提高直接定向技术精度及应用效果有重大意义ꎬ尤其是在自然灾害应急保障方面的应用更加广泛[1]ꎮ关键词:GPS差分定位ꎻ直接定向ꎻ应急保障
中图分类号:P228.4 文献标识码:A 文章编号:1672-5867(2019)06-0139-04
ImageQualityInspectionSystemforFrameDigitalAerial
PhotographyInstrument
(1.GuangzhouSouthernSatelliteNavigationInstrumentCo.Ltd.ꎬGuangzhou510663ꎬChinaꎻ
JIANGQinyu1ꎬDINGZhaolian2
2.TheSecondInstituteofPhotoprammetryandRemoteSensingꎬMinistryofNaturalResourcesꎬHarbin150081ꎬChina)Abstract:ThispapermainlystudiestheairborneGPSandIMUsystembasedonthehighprecisionꎬandusestheGPSdifferentialposi ̄
tioningandhigh-precisionIMUdataprocessingintheaerialphotographyoperationtoobtainthepositionandattitudeinformationꎬandtorestoretheexteriororientationelementsoftheaerialcamerainthephotography.Thismethodsimplifiestheoperationprocessꎬsavestheworkloadinthefieldandshortensthetimeperiod.Thispaperanalyzestheerrorsourcesofthisprocessꎬwhichisofgreatsignifi ̄cancetoimprovetheaccuracyandapplicationeffectofdirectgeolocationtechnologyꎬespeciallyintheapplicationofnaturaldisastere ̄mergencysupport.
Keywords:GPSdifferentialpositioningꎻdirectgeolocationꎻemergencysupport
0 引 言
目前ꎬ航空摄影摄站点的定位、定姿ꎬ主要通过传统空中三角测量进行的ꎮ该方法需要外业进行像片控制测量ꎬ内业空三加密ꎬ需耗费大量的设备成本、人员成本和经济成本ꎬ最重要的是消耗了大量的时间ꎮ这种作业模式无法满足当下航空摄影测量的需要ꎬ尤其是自然灾害应急保障的工作ꎮ当发生地震、山洪、台风等自然灾害时ꎬ效率是第一位的ꎮ直接定向技术通过快速获取航空摄影时摄站的瞬间位置、姿态信息ꎬ从而快速生成所需要的数字测绘成果ꎮ直接定位具有人工投入少、效率高、自动化程度高等特点ꎬ它将是今后的发展方向ꎮ
本文根据摄影测量和GPS/IMU系统的基本原理ꎬ详细分析了直接定向的方法、原理ꎬ总结了一套高效的检校
收稿日期:2018-08-13
方案并验证了其可靠性ꎮ讨论影响直接定向误差来源及优化方法ꎬ尤其是检校计算ꎮ通过对POS飞控系统与航摄仪时间误差、偏心角、偏心分量的检校ꎬ实现了对系统各项误差补偿ꎬ提高了精度ꎮ
1 机载GPS/IMU的直接定向原理
直接定向DG(directgeoreferencing)ꎬ通过传感器直接获取6个外方位元素(位置、姿态)ꎬ在不进行地面布控的情况下ꎬ建立起影像和地面定向关系ꎮ其最主要的方式是通过航摄检校场的集成检校ꎬ获取GPS天线、IMU和像机之间位置角度关系ꎬ从而对外方位元素进行补偿ꎬ取代整个测区的空中三角测量ꎬ实现像片直接定向ꎮ
POS(PositionandOrientationSystem)系统ꎬ即机载
GPS/IMU组合系统ꎮ其中ꎬIMU(惯性测量单元)可以独
作者简介:蒋勤玉(1988-)ꎬ男ꎬ辽宁营口人ꎬ工程师ꎬ学士ꎬ主要从事测绘地理信息技术应用工作ꎮ
140
测绘与空间地理信息
2019年
立获取位置、姿态、加速度等数据ꎬ将其与GPS组合应用可以消除IMU累计误差ꎮ同时ꎬGPS在航测平台上高速运动可能出现无法跟踪信号的现象ꎬ两者正好相互补充ꎬ提高位置、姿态的获取能力ꎮ
IMU首先ꎬ采用差分GPS方式获取位置信息(XꎬYꎬZ)ꎬ
按照集成检校方案获取姿态信息(ꎬ获取位置和姿态系统误差的补偿数ϕꎬωꎬκ)ꎬ相机获取影像数据ꎻ然后ꎬ据ꎬ通过直接定向模型进行转换ꎬ恢复航摄仪曝光时6个外方位元素ꎮ转换模型通过已知的航摄仪中心、GPS相位中心和IMU之间的关系进行坐标转换ꎬ得到精准的像片6个外方位元素ꎻ最后ꎬ将影像数据和外方位元素数据用于测ꎮ
图或正摄影像制图ꎮ其整个作业流程如图1所示[2]图1 直接定向流程图Fig.1 Directgeolocationflowchart
1.1 位置、姿态信息获取
差分GPS技术是当下定位精度最高的定位方法ꎬ以采用POS/AV610系统进行比例尺为1∶5000的航摄为例ꎬ其高程精度可以达到10—15cmꎬ水平精度在10cm以内IMUꎮ其成果可以应用到正摄影像制作系统工作原理示意图ꎮ
ꎮ图2为机载GPS/图2 机载GPS/IMU系统工作原理示意图Fig.2 Schematicdiagramoftheworkingprincipleof
airborneGPS/IMUsystem
1.2 集成系统检校
检校场的布设方案可以参照规范GB/T27919—2011
«IMU包含在摄区内或靠近摄区/GPS辅助航空摄影技术ꎬ且相对航高与摄区保持一致»进行布设ꎮ其范围应尽量ꎮ
4其航线及控制点布设方案主要包括航线9控制点法两种ꎬ均能达到检校目的2航线ꎮ6采用直接定控制点法和向法测图ꎬ必须在每次航摄前进行一次集成系统检校ꎮ
为保证地面控制点的精度ꎬ可以在检校范围内布设人工标志ꎬ同时要确保人工标志坐标准确、保存完好ꎬ而且可以清晰辨别ꎮ
2 误差来源及优化方法
机载GPS/IMU系统进行直接定向的误差来源主要包括空间位置偏心差、视准轴误差、相机与POS系统的误差、大气折射、地面曲率等ꎮ
2.1 空间位置偏心差
摄影中心三者在空间位置上无法重合而且存在一个相对
GPS天线相位中心、惯性单元IMU中心和航摄相机
固定的空间关系ꎬ即偏心分量ꎬ如图3所示ꎮ偏心分量可以通过全站仪量测出来ꎬ但由于有人工参与ꎬ所以存在较大误差ꎮ
图3 偏心分量示意图
Fig.3 Eccentricitycomponentdiagram
2.2 视准轴误差
航摄仪的空间辅助坐标系与惯性单元IMU三轴陀螺坐标系在X、Y、Z3个方向上无法平行ꎬ有角度的差异ꎬ即偏心角ꎬ如图4所示ꎮ视准轴误差必须通过检校场检校去除ꎮ
图4 像机坐标系与IMU坐标系夹角示意图Fig.4 Diagramoftheanglebetweentheimagecoordinate
systemandtheIMUcoordinatesystem
2.3 时间同步误差
历元时刻GPS、IMUꎬIMU和像机三者都是独立工作采用曝光脉冲记录ꎬ像机采用曝光时刻ꎬGPS采用卫星
ꎮ三者时间上的微小差异都会导致位置上有很大的差异从而影响外方位元素ꎬ现有的POS系统会对这个误差进行补偿ꎮ
2.4 大气折射和地球曲率引起的误差
大气在不同的高度ꎬ密度和压力也不同ꎬ引起大气折射现象ꎬ像点的坐标也会因此产生变化ꎮ地球曲率变化会引起像点位移ꎬ这需要通过像点坐标改正来消除带来的影响ꎮ
第6期
蒋勤玉等:基于机载GPS-IMU系统实现直接定向的研究
解获取线元素、角元素的改正值ꎬ去除系统误差ꎮ
141
2.5 集成系统优化方法
本文通过对直接定向原理和误差的来源探讨ꎬ结合个人实践操作ꎬ总结出如下集成系统的优化方法:精确量取GPS、IMU偏心分量ꎬ最好采用全站仪量测3次GPS偏心分量ꎬ取平均数ꎻ布设基站点ꎬ采用差分GPS解算POSꎻ集成的传感要进行刚性牢固连接ꎻ飞行平台保持平稳飞行ꎬ从而确保飞行平台不发生突然抖动ꎻ当使用POS辅助空中三角测量时ꎬ设计构架航线连接区域网ꎬ有助于提高系统如IGI、Applanix具有良好的精度ꎬ建议使用ꎮ精度[3]ꎻ严格执行检校ꎬ进行系统误差补偿ꎮ国外的POS
3.2 航飞检校场的设计
航飞检校场的布设方案有多种ꎬ现有的POS生产商加拿大Applanix和德国IGI公司各自提出检校场的布设方案ꎮ出于航测成本和数据处理的便利性考虑ꎬ推荐以下设计方案ꎮ
1)按东西方向航测两条平行的航线ꎬ且两条航线对2)航向重叠度和旁向重叠确保超过60%ꎮ这种飞行3)在检校场区域内布设9控制点和2个检查点ꎬ控制
向飞行ꎬ每条航线定点曝光12次ꎬ这样便可以去除偏心分量和时间同步的误差ꎮ
3 检校场的布设及精度分析
机载GPS/IMU直接获取的位置及姿态数据存在较大的系统误差ꎬ为了实现直接定向法ꎬ获取符合测图规范要求的高精度外方位元素ꎬ便需要进行检校计算去除系统误差ꎮ即应用摄影测量原理进行航空摄影ꎬ通过像片控制测量所得的高精度控制点进行GPS/IMU辅助空中三角测量ꎬ反求出机载GPS/IMU各组件间的关系ꎮ为了验证其精度ꎬ我们通过布设检校试验场证明其可行性ꎮ
当前主流的检校方法有多种ꎬ包括塔吊式地面检校和机载式地面检校ꎬ前者模仿实测的过程ꎬ虽然节约经济成本ꎬ但存在环境搭建困难、影响效率、机动性不够的问题ꎮ后者提供了高航高、快速度ꎬ这更容易增大系统误差ꎬ有利于更正改正参数ꎮ所以ꎬ我们在实际应用的过程中应用机载式地面检校ꎮ此法主要对偏心分量、视准轴、时间同步及影响精度的各种系统误差进行补偿ꎬ完成外方位元素线元素、角元素的改正ꎮ
方式有助于提高外方位元素精度ꎬ减少外业布控高程点ꎮ点间的距离小于3条基线ꎮ
3.3 航飞检校实验
本文以大庆附近的一个区域作为检校场ꎬ采用运-5飞机作为航摄平台ꎬ航摄仪选择UCWAꎬPOS系统选择POSAV510ꎬGPS系统采用天宝双频接收机ꎬ基站采用Trimble5700型ꎮ此检校场飞行设计如下:飞行方向为东GSD大小为7cmꎬ图5为航飞检校场的轨迹图ꎮ地面控制点布设6个3°重叠控制点和2个检查点ꎮ
西向ꎬ相对航高500mꎬ航向重叠60%ꎬ旁向重叠60%ꎬ
图5 航飞检校场飞行轨迹图Fig.5 Flightflighttrajectorycontrolfield
3.1 具体流程
选定一块区域作为飞行检校场ꎬ设计航测飞行计划ꎬ同时在建立控制点布设方案ꎻ采用集成系统(机载GPS/IMU、航摄相机、座架平台)外业飞行ꎬ同时获取控制点数进行空中三角测量ꎬ获取摄站点定位、定姿数据ꎬ迭代求
利用pospac软件对pos系统进行解算ꎬ利用pospac检校模块获取外方位元素的线元素和角元素改正值并予以修正ꎬ获取准确外方位元素ꎮ通过photomod软件获取空三加密后的外方位元素ꎬ将两种外方位元素成果进行比对ꎬ位置的偏差结果见表1ꎬ表2为姿态角的偏差ꎮ
据ꎻ借助GPS/IMU辅助航空摄影测量ꎬ在数据后处理后ꎬ
表1 POS直接定向数据与空三后的外方位元素位置坐标对比
Tab.1 ComparisonofthepositioncoordinatesofthePOSdirectdirectionaldataandaerialtriangulation
POS直接定位坐标
Y32×××38.62132×××73.54732×××09.04632×××95.80732×××83.19632×××81.26732×××86.12532×××97.83232×××16.71432×××41.11332×××27.49
空三后坐标
Y32×××38.63232×××73.59632×××09.22232×××95.79132×××83.23932×××81.26932×××86.15332×××97.86932×××16.83332×××27.61232×××41.28
差值
像片号102103104105106107203204205206207
X58××79.77258××34.40758××90.42158××46.23258××01.87558××56.45159××35.02258××78.86958××23.52958××67.82358××24.35
Z4982.7964983.0964982.8594984.8034982.3564981.9014986.4534971.3594971.4874978.6124981.948
X58××24.32358××79.80858××34.73458××90.45458××46.24958××01.84958××56.51359××34.99358××78.81558××23.49358××67.82
Z4982.9554983.3464983.0684984.9824982.5744982.1094971.1284971.2514978.3514981.6794986.67
ΔX(M)-0.036-0.327-0.033-0.017-0.0620.0290.0540.0360.0030.0260.027
ΔY(M)-0.011-0.049-0.176-0.043-0.002-0.028-0.037-0.119-0.167-0.1220.016
ΔZ(M)-0.159-0.209-0.179-0.218-0.208-0.2170.2310.2360.2610.269-0.25
142
测绘与空间地理信息
2019年
表2 POS直接定向数据与空三后的外方位元素姿态角对比
Tab.2 ComparisonofattitudeanglesbetweenPOSdirectdirectionaldataandaerialtriangulation
像片号102103104105106107203204205206207
Omega(°)
POS直接定位后的姿态角
0.29367359.430930.286740.320850.29835359.716770.279000.28773Phi(°)
Kap(°)
Omega(°)-0.14587-0.12064-0.12689-0.122620.077020.086690.308940.12055-0.13060.51588-0.13031
空三后的姿态角
0.28886-0.569010.285510.316970.29265-0.26183-0.25346-0.247520.012350.273260.28569Phi(°)
Kap(°)
ΔOmega(°)0.016450.019080.014860.016010.006680.010020.007880.0070.01750.01570.01646
差值ΔPhi(°)0.002040.001230.003880.005740.002140.002150.002180.002410.00570.00040.00481
ΔKap(°)
359.88615359.87058359.89506359.88426359.88912359.894880.093690.318960.128430.08370.53496
91.3570291.7604989.3650589.9203590.24048269.55099269.07894269.55135271.4028390.4681890.4653
90.4743989.3729689.9272690.24656-90.43879-90.91168-90.43854-88.5877590.475
91.7708
91.3671-0.01008-0.01031-0.00909-0.00791-0.00691-0.00608-0.006820.010220.009380.010110.00942
359.72505359.99055359.72865
通过实际数据对比ꎬ证明了在进行完检校计算后ꎬ利用直接定向获取的外方位元素与到空三后获取的外方位差值很小ꎬ可以满足摄影测量要求ꎮ
0.105624mꎮ根据«数字航空摄影测量空中三角测量»规范的要求ꎬ在进行1∶500比例尺成图时ꎬ平地区域检查点的平面中误差需优于0.175mꎬ高程中误差需优于0.15mꎬ这一结果符合空三精度要求ꎮ实验结果证明了ꎬ利用直接定向获取的外方位元素可以达到空三的精度ꎬ利用检校计算获取修正参数的正确性与可靠性[4]ꎮ
3.4 直接定向模式下的空三精度
在本次实验测区内ꎬ实地量测30个平高点作为检查点ꎮ将采集的检查点导入Photomod软件的立体交互界面下进行刺点ꎬ采集对应点位的平面及高程数据ꎮ同时ꎬ实测点与立体测图点进行比对ꎬ其部分差值数据见表3所示ꎮ
表3 实测点与立体测图点差值
Tab.3 Pointdifferencebetweenthemeasuredpoint
andthestereoscopicmap
ΔX(m)-0.0266-0.01540.00350.00540.0842
ΔY(m)-0.0654-0.0811-0.0453-0.0980-0.0913-0.0396-0.0632-0.0956-0.0984-0.0983-0.0332-0.0652
点号p1p2p3p4p5p6p7p8p10p11p12p9
ΔZ(m)0.06360.03460.09460.01650.03470.09870.07450.03430.04930.38670.13650.168
4 结束语
本文探讨了直接定向的原理ꎬ总结了影响其精度的各方面误差ꎬ并提出优化方法ꎮ归纳了一套高效检校方案ꎬ特别指出要实现此项技术必须完成GPS、IMU、相机间系统误差的检校ꎮ直接定向为航空摄影测量提供了新的思路ꎬ机载GPS-IMU系统精确获取摄影时航片的外方位元素数据(3个位置参数和3个姿态参数)ꎬ从而可以大大减少外业控制点的布设ꎮ这样就可以实现在一定比例尺下无控地理定位ꎬ满足测图精度要求ꎬ节省了时间、人力、物力、财力ꎮ同时它也是处理自然灾害应急保障最有效、最可行的解决方法ꎮ
参考文献:
[1] 李学友.IMU/GPS辅助航空摄影测量综述[J].测绘科[2] 李学友.IMU/GPS辅助航空摄影测量原理、方法及实践[3] 王建超ꎬ郭大海ꎬ郑雄伟.机载POS直接地理定位软件[4] 杨娜.国产POS与数字航空相机集成检校的研究[D].青
岛:山东科技大学ꎬ2012.28-32.
[D].郑州:解放军信息工程大学ꎬ2005.
AeroDG的设计与实现[J].国土资源遥感ꎬ2007(4):学ꎬ2005ꎬ30(5):110-113.
-0.0124-0.0684-0.02940.09530.06540.0341
-0.0205
根据数据可得ꎬ在直接定向法模式下ꎬ检查点的平面位置中误差为0.09142mꎬ高程位置中误差为
[编辑:张 曦]
因篇幅问题不能全部显示,请点此查看更多更全内容