基于数字高程模型高程快速迭代的航拍图像目标定位方法

A target location method for aerial images through fast iteration of elevation based on DEM

在大倾角航空相机对地面目标定位过程中,借助数字高程模型(DEM)可有效解决地球椭球模型定位存在的大地高误差影响。为获取地面坐标的准确信息特别是高程信息,首先,根据载机的位置姿态信息以及航空相机的框架角等信息利用齐次坐标变换求解出成像系统视轴在地理坐标系下的指向,再利用数字高程模型确定目标点的坐标。针对成像过程中目标点高程计算繁琐、容易不迭代等问题,提出了一种对目标高程值进行快速迭代的方法。通过对目标区域高程进行折半查找处理,计算该处视轴光线高程与地面高程差值。继续计算该高程差中值并继续迭代,直到小于一定阈值。最后使用蒙特卡洛分析法对整个成像过程存在的误差项进行分析。实验结果表明:采用快速迭代法进行计算,当收敛阈值为十分之一DEM网格精度时,迭代效率提升45.5%,收敛速度大大提高;且通过数字高程模型计算,在飞行高度为15409 m,相机框架角大于74°时,对于山地区域目标的圆概率误差小于200 m,可以满足实际工程需要。

In the positioning process of aerial cameras with large inclination angles,the influence of height error in the earth ellipsoid model can be effectively solved with the help of a digital elevation model(DEM).This is very important for obtaining accurate ground coordinates,especially elevation.Firstly,the orientation of the line-of-sight angle in the geographic coordinate system is solved by transforming homogeneous coordinates according to the position and attitude information of the carrier aircraft and the frame angle information of the aerial camera,and then the longitude and latitude of the target point are determined by a digital elevation model.To overcome the tedious nature of calculating target elevation and the non-convergence in the imaging process,a fast iterative method is proposed to iterate over the target elevation’s value.The difference between the light elevation of the visual axis and the ground elevation is calculated by halving the target elevation.The median elevation difference is calculated iteratively until it is less than a certain threshold.Finally,Monte Carlo analysis was used to analyze the error terms in the whole imaging process.When the convergence threshold is 1/10 DEM in grid accuracy,the iteration efficiency increases by 45.5%and the convergence speed is greatly improved.Through the calculation of the digital elevation model,when the flight height is 15,409 meters and the camera frame’s angle is greater than 74°,a mountainous area’s target circular error probability is less than 200 m which meets the real engineering needs.