一种基于R-D模型的侧视雷达影像严密定位方法技术

一种基于R-D模型严密定位多源数据联合平差方法，其特点是通过对R-D模型的改造，建立了侧视雷达影像像点坐标参数为因变量的显函数形式构像方程和误差方程，能够根据误差方程式直接计算像点量测坐标理论误差；通过建立的侧视雷达影像像点坐标观测值，联合POS观测值、传感器检校参数、地面点坐标等其它观测值的误差方程组成的误差方程组，实现侧视雷达影像观测值与非雷达影像观测值的联合平差，并根据相关观测值的先验精度对各观测值进行定权。本发明专利技术对于提高稀少控制点条件下R-D模型空中三角测量平差解算的稳定性、可靠性、精度、计算效率有重要作用。

【技术实现步骤摘要】
一种基于R-D模型的侧视雷达影像严密定位方法
本专利技术属于摄影测量领域，特别涉及到侧视雷达遥感影像的几何处理技术，具体的说是涉及一种基于距离一多普勒(R-D)模型的侧视雷达影像严密定位方法。
技术介绍
雷达影像包括真实孔径雷达影像和合成孔径雷达(SAR)影像。侧视雷达遥感影像的严格几何构像模型主要有基于传感器状态矢量、多普勒参数的模型和基于传感器状态矢量和姿态参数(外方位元素)的模型两类。在各式各样的雷达影像构像模型中，R-D模型应用得最为广泛。以R-D模型为基础，国内外对无控制点、稀少控制点及高山地合成孔径雷达影像的纠正、SAR影像的立体定位等各种几何处理均进行了较充分的研究。近10年来， SAR影像的区域网平差处理也逐渐得到重视，并对其进行了相关研究，虽然雷达影像定位技术得到了广泛的研究，但与光学遥感影像摄影测量软件相比，国内外仍缺少广泛应用的能够稳健、可靠地处理雷达遥感影像的摄影测量数据处理软件，满足稀少控制点条件下高分辨率机载SAR影像中大 比例尺地形测图的研究仍鲜有文献报导。稳健性好，适用性强的SAR 影像几何处理方法和系统仍然是未来一段时间摄影测量工作者重要的研究目标。SAR遥感影像的定位虽然已有较深入的研究，定位精度也逐渐达到了测绘中大比例尺地图精度的要求。R-D模型是SAR影像定位应用得最广泛的模型，其由距离等式和多普勒等式构成i (x-xv):-(；-;；r + (z-zv)2 = ^2iI Vxs (Xy-X) + VkS Fs-F) + VzXZs-Z) = -fMj 丨 2其中(X，Y，Z)为地面点坐标，(Xs, Ys, Zs)为传感器位置，(VXs, VYs, Vzs为传感器的速度，fD表示成像时刻多普勒中心频率，λ为发射脉冲波波长，&为像点对应的雷达测量距离。将上式改写成如下等式, Z1 = Cr- Xs)1 + (Y- Ys f + (Z- zv)2 -R = ο-^) + VM-O + VzXZs — ^ + Io^Rj/2 = 0像点坐标对应的地物点地心直角坐标系坐标(X, Y, Z)的计算,通过式上式联合考虑高程因子的精度无损地球椭球模型来进行，地球椭球模型为X2+ Y2 Z2 t-T +-- = I{β+Η)1其中，A为地球椭球的长半轴，B为地球椭球的短半轴，H为地面点的高程，从DEM 数据中提取。对于 WGS84 坐标系，A = 6378137. Om, B = 6356752. 314m。R-D模型误差方程的一般形式以传感器位置、速度、地面点坐标和像点坐标系统误差作为未知数，对变形后的误差方程进行线性化，得到相应像点坐标观测值的误差方程式IvX = fRX-Xs + frYs^S + fRZ^Zs + fgy^Vy + ZgnJFy + fRVjVZ + fgjflX jT fRydY + f- 1RI^= IdJXs+fDrsdYs+fDZsdZs+fDrJVx + fDVYdVY + fD1^Vz + fDXdX + fd￥dY + fDZdZ- Id当前，遥感影像一般都带有高精度POS数据和星历数据，影像数据与这些辅助数据进行联合处理，是提高稀少或无地面控制点条件下影像定向参数求解精度与稳定性的有效方法。R-D模型符合SAR影像的成像机理，其本身是严密的，但是，其误差方程直接根据 R-D模型建立的距离条件和多普勒条件来建立，不是针对像点坐标观测值。对于光学遥感影像，严密摄影测量根据共线方程展开，共线方程可表示成如下函数w) = fxXXs, Ys, Ζν,ω^φ,κ,尤 K Ζ)式中，X、y为像点坐标，X、Y、Z为地面点坐标，Xs、Ys、Zs)为传感器天线的空间坐标，ω、φ、K为姿态角。对误差方程进行线性化，得到相应像点坐标观测值的误差方程式[4 = f^dXs +fxZsdZs+fxJm + /邱却 + /Jk + 厶I+ fK…dx - 4U = UdXs+fyrsdYs+fyZsdZs++ /1ψφ + + fyXdX^ /^F+ /y^Z+ dy- Iy由于建立了像点坐标为显函数形式的误差方程，可以进行像点定向残差的快速计算，并根据像点量测先验精度对其进行定权。同时，光学影像的光束法平差根据像点量测坐标误差的最小二乘(残差平方和最小)来实现定向参数的解算，是一种最严密的影像空中三角测量方法。对于目前的R-D方程进行雷达影像的定位，由于像点坐标信息隐藏在R-D方程中， 没有建立起地面点坐标与像点坐标间明确的函数关系，所建立的误差方程不是像点量测坐标(即SAR影像测量观测值)的误差方程，意味着像点坐标的定向残差不能由误差方程直接计算得出，不利于像点量测或匹配坐标粗差的剔除；误差方程根据条件式来建立，没有直接的精度信息，雷达影像数据与其它数据联合平差时定权带来不便；根据定向参数计算得到的像点坐标定向残差(像点坐标改正数向量)直接体现了定位精度的高低，光学遥感影像的光束法平差正是利用这一特点，根据像点量测坐标误差的最小二乘(残差平方和最小)来实现定向参数的解算，是一种最严密的影像空中三角测量方法。而根据R-D条件式建立的虚拟观测值最小二乘求解的结果，理论上不等同定向参数的解算结果也能使得定像点坐标的定向残差也能实现最小二乘。
技术实现思路
本专利技术的目的在于提高R-D (距离_多普勒)模型严密定位时联合平差各观测值的先验定权精度、方便误差探测，提高稀少控制点条件下真实孔径雷达影像或合成孔径雷达(S A R)影像R-D模型严密定位的精度与平差解算的稳定性、可靠性。本专利技术根据给出的R-D误差方程，根据误差方程式直接计算像点量测坐标理论误差，并根据先验精度实现各观测值的定权，通过将原R-D方程转变为像点坐标为显函数形式的方程，使得像点坐标作为摄影测量观测值的属性得到体现，实现雷达影像观测值的先CN 102914771 A书明说3/7页验定权和像点定向残差的快速计算；联合机载全球定位系统(GPS)测量航迹和速度精化模型参数虚拟观测值和传感器系统差检校参数虚拟观测值的误差方程，实现高精度雷达影像的定位与定向。本专利技术采用计算机通过以下步骤进行(I)将R-D方程表达式转变为像点量测坐标参数为因变量的显函数形式方法为 对R-D模型的多普勒方程两边乘以系数Cx，距离变形为y为因变量的显函数形式，即将R-D 方程| 2[(,!--+ (； - +(Z-Zs)V^ilR -fn (X- As)2 + (,- }Wf + (Z-Zsf = {J% + Mr变形为雷达像点坐标显函数形式的R-D方程(x =Jcr-[(A^-As)Fli +{Y-Ys)Vri +{Z-Zs)Vzy{lR)-Cr-Zp^\ι-[ y = I 扎r-Xsf + (Y- Ysf + (Z- Zsf — Mr其中，Cx按下式计算Γ XiBn +Mr- r)Cx =-1~X 2\V\Ma 2\J \Mcx为非0比例系数，X、Y、Z为地面点坐标，Xs、Ys、Zs、VXs、VYs、Vzs为雷达传感器天线中心状态矢量，λ、fD、Rj为雷达波长、多普勒频率及相应像点雷达实测距离。Rtl为初始斜距，y为像点列坐标，χ为方位向像点坐标，Mr为斜距分辨率，Ma为方位向分辨率，|V|为传感器速度大小标量。像点本文档来自技高网...

【技术保护点】
一种基于R？D模型的侧视雷达影像严密定位方法，其特征在于：根据R？D模型建立了影像观测值的误差方程，按其误差方程式直接计算像点量测坐标理论误差，并根据先验精度实现各观测值的定权，包括以下步骤：（1）将R？D模型改变为像点坐标参数为因变量的显函数方程；（2）建立传感器位置、速度观测值的精化模型，并将该模型代入到步骤(1)建立的方程中,构建含有精化模型参数的显函数形式的构像方程；（3）以传感器位置、速度精化模型参数、地面点坐标参数和影像方位、距离向系统差等参数为未知数进行线性化，建立基于R？D模型的雷达影像像点坐标观测值误差方程；（4）将步骤（3）中像点坐标误差方程与POS观测值的精化模型参数、传感器检校参数、地面点坐标等观测值或虚拟观测值的误差方程一起组成误差方程组，并根据相关观测值的先验精度对各观测值进行定权；（5）迭代解算各种定向未知数、地面点坐标值，并根据误差方程计算像点定向残差，实现粗差的剔除；（6）根据定向参数实现雷达影像的几何处理与对地定位。

【技术特征摘要】

【专利技术属性】
技术研发人员：程春泉，张继贤，黄国满，赵争，杨书成，
申请(专利权)人：中国测绘科学研究院，
类型：发明
国别省市：