CN105466400B - 一种利用rpc探测多源卫星影像同名像点成像交会角方法 - Google Patents
一种利用rpc探测多源卫星影像同名像点成像交会角方法 Download PDFInfo
- Publication number
- CN105466400B CN105466400B CN201510940824.1A CN201510940824A CN105466400B CN 105466400 B CN105466400 B CN 105466400B CN 201510940824 A CN201510940824 A CN 201510940824A CN 105466400 B CN105466400 B CN 105466400B
- Authority
- CN
- China
- Prior art keywords
- msub
- msup
- mrow
- mtr
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/30—Interpretation of pictures by triangulation
- G01C11/34—Aerial triangulation
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Image Processing (AREA)
Abstract
一种利用RPC探测多源卫星影像同名像点成像交会角方法,是一种基于有理多项式(RPC)参数和全球SRTM数据、适合同源或异源卫星遥感影像同名像点成像交会角的探测方法,在全球SRTM数据的支持下,通过高程的改变,利用RPC模型来计算地面点位置矢量的变化,进而计算地面目标的入射角,以及任意两景影像同名像点地面目标的成像交会角。本发明在不知晓传感器位置的情况下,利用RPC参数可以获得重叠影像同名点的成像交会角,为RPC模型定位成像交会角探测提供一种有效的方法,能够在影像几何处理中预防因弱交会成像而导致的地面点坐标求解病态问题提供基础信息。
Description
技术领域
本发明涉及卫星影像摄影测量方法,特别涉及到利用有理多项式模型进行空中三角测量或三维重建时探测多源卫星影像同名像点成像交会角探测方法。
背景技术
基高比或成像交会角反映了影像像点对应三维地面点交会求解的稳定性。基高比是摄影基线长度与摄站高度之比。摄影基线是多站摄影时,遥感器连续两次曝光瞬间镜头中心间的距离,一般由曝光时间间隔来控制;摄站高度系镜头中心沿铅垂线到地面的距离,即相对航高。基高比决定相邻影像在航向上的重叠度。基高比愈大,航向重叠愈小,反之重叠度愈大。而影像航向重叠愈大,立体观测和高程量测精度愈低。重叠度大,相邻两幅影像组成立体像对时构成的立体模型,其垂直比例尺夸大不明显,亦即立体效应差,影响观测。成像交会角是地面点与立体遥感影像成像时相应传感器位置间的夹角,交会角的大小决定了遥感影像目标对地三维定位时的精度与稳健性。基高比和成像交会角是判断经典光学遥感影像进行定位质量的重要参数,满足特定基高比或成像交会角要求也是经典影像进行摄影测量制图时需要满足的必须条件。
RPC(有理多项式参数)模型目前正成为星载遥感影像定位的主流方法,但RPC参数中不包含传感器的位置,使得传感器至地面目标的光线或电磁波入射向量难以使用传统的方法计算。因此,利用RPC参数进行成像交会角的探测,对RPC进行影像几何处理,无疑有着重要的作用。
卫星影像在成像的过程中,受到透视投影、摄影轴倾斜、大气折光、地球曲率及地形起伏等诸多因素影响,致使影像中各像点产生不同程度的几何变形而失真,无法直接与真实地图套合以进行后续应用。因此要精确地使用遥感卫星影像,首要的工作当先解决影像之传感器的几何模型问题。通常建立传感器的模型的方法主要可分成两大类,除严密模型外,通用模型在近几年受到前所未有的关注,其主要进展是有理多项式模型逐渐替代传统的严密模型,用于光学与雷达卫星遥感影像的几何处理中。
常用的多项式模型比较,RPC模型实际上是各种传感器几何模型的一种抽象的表达方式,它适用于各类传感器包括最新的航空和航天传感器,是多项式模型更精确的形式。它以像点坐标(R,C)表达为以相应地面点空间坐标(L,B,H)为自变量的多项式比值:
其中:
NumR(Ln,Bn,Hn)=a1+a2Bn+a3Ln+a4Hn+a5BnLn+a6BnHn+a7LnHn+a8Bn 2+a9Ln 2+a10Hn 2+a11BnLnHn+a12Bn 3+a13BnLn 2+a14BnHn 2+a15Bn 2Ln+a16Ln 3+a17LnHn 2+a18Bn 2Hn+a19Ln 2Hn+a20Hn 3
DenR(Ln,Bn,Hn)=b1+b2Bn+b3Ln+b4Hn+b5BnLn+b6BnHn+b7LnHn+b8Bn 2+b9Ln 2+b10Hn 2+b11BnLnHn+a12Bn 3+a13BnLn 2+b14BnHn 2+b15Bn 2Ln+b16Ln 3+b17LnHn 2+b18Bn 2Hn+b19Ln 2Hn+b20Hn 3
NumC(Ln,Bn,Hn)=c1+c2Bn+c3Ln+c4Hn+c5BnLn+c6BnHn+c7LnHn+c8Bn 2+c9Ln 2+c10Hn 2+c11BnLnZn+c12Bn 3+c13BnLn 2+c14BnHn 2+c15Bn 2Ln+c16Ln 3+c17LnHn 2+c18Bn 2Hn+c19Ln 2Hn+c20Hn 3
DenC(Ln,Bn,Hn)=d1+d2Bn+d3Ln+d4Hn+d5BnLn+d6BnHn+d7LnHn+d8Bn 2+d9Ln 2+d10Hn 2+d11BnLnHn+d12Bn 3+d13BnLn 2+d14BnHn 2+d15Bn 2Ln+d16Ln 3+d17LnHn 2+d18Bn 2Hn+d19Ln 2Hn+d20Hn 3
式中(Ln,Bn,Hn)为地面点(L,B,H)的归一化物方坐标系坐标,NumR(Ln,Bn,Hn),DenR(Ln,Bn,Hn),NumC(Ln,Bn,Hn),DenC(Ln,Bn,Hn)为三元三次多项式,ai,bi,ci,di(i=1….20)为影像附带的RPC参数,b1=1,d1=1;(Rn,Cn)为像点影像坐标(R,C)的归一化影像坐标,根据下式计算:
其中Roff,Coff,Loff,Boff,Hoff是像点的影像坐标和物方大地坐标的偏移参数,Rscale,Cscale,Lscale,Bscale,Hscale为像点和地面点坐标的比例系数。
与多项式模型一样,RPC模型参数的计算本质上是通过采样点来拟合。这里的采样点是计算RPC时已知的像点和地面点坐标的集合:
{Ri,Ci,Li,Bi,Hi}(i=1,2,3…为第i个已知点)
根据已知点的来源,RPC模型参数求解有与地形无关和与地形相关两种求解方式。在已有成像几何模型的情况下,采用与地形无关的求解方式,通过已有模型(通常是严密模型)产生虚拟的控制点;在严密模型未知的情况下,则采用与地形相关的求解方式,通过测量给定一定数目的控制点。
RPC模型在数据处理方式上虽然有与传感器和坐标系无关的特点。目前,已有大量文献对不同传感器影像RPC模型的精度进行了测试与应用,在一些高质量的高分影像如IKONOS和QuickBird中,通过平差定位精度可达到子像素级。
由于RPC模型属于传感器无关模型,它非常适合不同卫星影像的联合处理,在传感器摄影时摄站位置未知的情况下,探测地面目标的成像交会角和不同影像的弱交会条件,有重要的作用,但目前未出现有效的方法,本发明正是为解决这一难题而提出。
发明内容
本发明提供了一种利用RPC探测多源卫星影像同名像点成像交会角方法,是一种基于RPC参数和全球SRTM(航天飞机雷达地形测绘)数据的多源影像成像交会角的探测方法,在全球SRTM数据的支持下,通过地面点高程的变动,来计算地面点位置矢量的改变,进而计算任意影像像点对应的地面入射角,根据入射角计算任意两景影像同名像点对应的地面目标的成像交会角,并根据影像的性质判断多源影像的同名点是否属于弱交会成像。该方法通过以下步骤实现:
(1)利用卫星影像附带的RPC参数,高程采用全球SRTM数据,利用RPC模型计算各影像像点p(R,C)对应地面点P0的大地坐标经度、纬度和高程(L0,B0,H0),R、C为像点p的影像行、列坐标;
(2)依次取新的高程值H1=H0+1和H2=H0+1000,利用RPC模型重新计算得到p点对应地面点P1和P2的大地坐标经度、纬度和高程(L1,B1H1)、(L2,B2,H2);
(3)将P0、P1和P2的大地坐标转换为地心直角坐标系三维坐标,分别得到P0(X0,Y0,Z0)、P1(X1,Y1,Z1)、P2(X2,Y2,Z2),并计算地心直角坐标系中的向量和向量
(4)若p为光学卫星影像上的像点,对应地面点的入射向量若p为合成孔径雷达影像上的像点,入射向量为
(5)计算任意重叠影像1和2上的同名像点对应的成像交会角β=arccos[α1·α2/(|α1|·|α2|)],其中α1、α2分别按步骤(4)计算得到的影像1和2中同名像点对应地面点的入射向量;
(6)根据影像性质和成像交会角大小,判定同名点是否属于弱交会成像。
所述步骤(1)中,影像坐标为(R,C)的像点p对应地面点大地坐标P0(L0,B0,H0)的计算根据影像RPC参数和SRTM数据,利用RPC模型通过迭代的方法计算得到,包括以下步骤:
(1)给定影像点对应的地面点P0所在区域的近似高程h0作为初值,即令H0=h0;
(2)将像点坐标(R,C)、高程H0和RPC参数代入RPC模型,得到述方程组并求解地面点的归一化大地坐标经度和纬度(Ln,Bn):
其中,NumR(Ln,Bn,Hn),DenR(Ln,Bn,Hn),NumC(Ln,Bn,Hn),DenC(Ln,Bn,Hn)代表三元三次多项式,且有:
NumR(Ln,Bn,Hn)=a1+a2Bn+a3Ln+a4Hn+a5BnLn+a6BnHn+a7LnHn+a8Bn 2+a9Ln 2+a10Hn 2+a11BnLnHn+a12Bn 3+a13BnLn 2+a14BnHn 2+a15Bn 2Ln+a16Ln 3+a17LnHn 2+a18Bn 2Hn+a19Ln 2Hn+a20Hn 3
DenR(Ln,Bn,Hn)=b1+b2Bn+b3Ln+b4Hn+b5BnLn+b6BnHn+b7LnHn+b8Bn 2+b9Ln 2+b10Hn 2+b11BnLnHn+a12Bn 3+a13BnLn 2+b14BnHn 2+b15Bn 2Ln+b16Ln 3+b17LnHn 2+b18Bn 2Hn+b19Ln 2Hn+b20Hn 3
NumC(Ln,Bn,Hn)=c1+c2Bn+c3Ln+c4Hn+c5BnLn+c6BnHn+c7LnHn+c8Bn 2+c9Ln 2+c10Hn 2+c11BnLnZn+c12Bn 3+c13BnLn 2+c14BnHn 2+c15Bn 2Ln+c16Ln 3+c17LnHn 2+c18Bn 2Hn+c19Ln 2Hn+c20Hn 3
DenC(Ln,Bn,Hn)=d1+d2Bn+d3Ln+d4Hn+d5BnLn+d6BnHn+d7LnHn+d8Bn 2+d9Ln 2+d10Hn 2+d11BnLnHn+d12Bn 3+d13BnLn 2+d14BnHn 2+d15Bn 2Ln+d16Ln 3+d17LnHn 2+d18Bn 2Hn+d19Ln 2Hn+d20Hn 3
式中ai,bi,ci,di(i=1….20)为影像附带的RPC参数;Rn为影像行坐标R的归一化值,Cn为影像列坐标C的归一化值,Hn为地面点高程H0的归一化值,它们根据下式计算:
其中,Roff、Coff、Hoff为像点影像行向坐标、列向坐标和地面点高程的平移参数,Rscale、Cscale、Hscale为像点影像行向坐标、列向坐标和地面点高程的缩放系数,它们均取自影像附带的RPC参数值;
(3)在得到像点的归一化坐标(Ln,Bn)后,根据下式求解地面点的大地坐标(L0,B0):
上式中,Loff、Boff为地面点大地坐标经度和纬度的平移参数,Lscale、Bscale为地面点大地坐标经度和纬度的缩放系数,它们均取自影像附带的RPC参数;
(4)将大地坐标(L0,B0,H0)转换到横轴墨卡托地图投影坐标系中,得到坐标(x0,y0),并从SRTM中内插出该投影坐标对应的高程值h,给定阈值10,若|h-H0|<10,则(L0,B0,H0)即为P0点的大地坐标值,终止迭代,否则,取H0=h,重回步骤(2)。
所述步骤(3)中,P0、P1和P2的大地坐标转换到地心直角坐标系中,坐标转换时所采用的椭球参数为WGS-84椭球参数,转换按下述公式进行:
上式中,(X,Y,Z)为地面点在地心直角坐标系中的坐标,(L,B,H)为大地坐标的经度、纬度和高程,N为地球椭球的卯酉圈曲率半径,a,b,e分别为地球椭球的长半轴、短半轴和第一偏心率,其中a=6378137,b=6356752.3142;
所述步骤(6)中,弱交会判断方法为首先给定阈值σ,并根据影像性质进行判定:
(1)当同名点所在影像对同为光学影像,或同为合成孔径雷达影像时:若交会角|β|<σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像;
(2)当同名点所在影像对一景为光学影像,另一景为合成孔径雷达影像时:若交会角满足|β-90|≤σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像。
RPC模型定位成像交会角的探测方法,其特征是在RPC参数和SRTM数据的支持下,通过改变地面点的高程值,计算地面点的变化向量,并利用变化向量计算光学影像和雷达影像的对应地面点入射向量,以及任意同源或异源影像同名像点对应目标的成像交会角。
附图说明
图1是本发明提供的利用RPC探测多源卫星影像同名像点成像交会角方法的示意图。
具体实施方式
所述利用RPC探测多源卫星影像同名像点成像交会角方法是在全球SRTM数据的支持下,通过地面点高程的变动,计算地面点位置矢量的改变,进而计算任意影像像点对应的地面入射角,以及任意两景影像同名像点对应的地面目标的成像交会角。如图1所示,该方法通过以下步骤实现:
步骤1,利用卫星遥感影像附带的RPC参数,高程采用全球SRTM数据,利用RPC模型计算各像点p(R,C)对应地面点P0的大地坐标经度、纬度和高程(L0,B0,H0);大地坐标P0(L0,B0,H0)的计算根据RPC参数和SRTM数据,利用RPC模型通过迭代的方法计算得到,包括以下步骤:
1)在SRTM中,给定P0所在区域中附近的近似高程h0,并令H0=h0;
2)利用像点影像坐标(R,C)和RPC参数构建有理多项式模型,并按下式所构建的方程组求解物方坐标经纬度的归一化坐标(Ln,Bn):
其中NumR(Ln,Bn,Hn),DenR(Ln,Bn,Hn),NumC(Ln,Bn,Hn),DenC(Ln,Bn,Hn)为三元三次多项式,:
NumR(Ln,Bn,Hn)=a1+a2Bn+a3Ln+a4Hn+a5BnLn+a6BnHn+a7LnHn+a8Bn 2+a9Ln 2+a10Hn 2+a11BnLnHn+a12Bn 3+a13BnLn 2+a14BnHn 2+a15Bn 2Ln+a16Ln 3+a17LnHn 2+a18Bn 2Hn+a19Ln 2Hn+a20Hn 3
DenR(Ln,Bn,Hn)=b1+b2Bn+b3Ln+b4Hn+b5BnLn+b6BnHn+b7LnHn+b8Bn 2+b9Ln 2+b10Hn 2+b11BnLnHn+a12Bn 3+a13BnLn 2+b14BnHn 2+b15Bn 2Ln+b16Ln 3+b17LnHn 2+b18Bn 2Hn+b19Ln 2Hn+b20Hn 3
NumC(Ln,Bn,Hn)=c1+c2Bn+c3Ln+c4Hn+c5BnLn+c6BnHn+c7LnHn+c8Bn 2+c9Ln 2+c10Hn 2+c11BnLnZn+c12Bn 3+c13BnLn 2+c14BnHn 2+c15Bn 2Ln+c16Ln 3+c17LnHn 2+c18Bn 2Hn+c19Ln 2Hn+c20Hn 3
DenC(Ln,Bn,Hn)=d1+d2Bn+d3Ln+d4Hn+d5BnLn+d6BnHn+d7LnHn+d8Bn 2+d9Ln 2+d10Hn 2+d11BnLnHn+d12Bn 3+d13BnLn 2+d14BnHn 2+d15Bn 2Ln+d16Ln 3+d17LnHn 2+d18Bn 2Hn+d19Ln 2Hn+d20Hn 3
式中ai,bi,ci,di(i=1….20)为影像附带的RPC参数;(Rn,Cn)为像点影像坐标(R,C)的归一化坐标,Hn为地面点高程H的归一化坐标,它们根据下式计算:
上式中,(Roff,Coff,Hoff)为像点影像坐标和地面点高程的平移参数,(Rscale,Cscale,Hscale)为像点坐标和地面点高程的缩放系数,它们均取自影像附带的RPC参数;
3)在得到(Ln,Bn)后,根据下面公式求解地面点的大地坐标(L0,B0):
上式中,(Loff,Boff)为地面点大地坐标经度和纬度的平移参数,(Lscale,Bscale)为地面点大地坐标经度和纬度的缩放系数,它们均取自影像附带的RPC参数;
4)将大地坐标(L0,B0,H0)转换到UTM投影坐标系中,得到坐标(x0,y0),并从SRTM中内插出该平面坐标对应的高程值h,给定阈值10,若|h-H0|<10,坐标(L0,B0,H0)即为P0点的坐标,终止迭代,否则,取H0=h,重回步骤2)。
步骤2,分别取新的高程值H1=H0+1和H2=H0+1000,利用RPC重新计算像点p对应地面点的大地坐标P1(B1,L1,H1)和P2(B2,L2,H2);
步骤3,将P0、P1和P2的大地坐标转换到地心直角坐标中,坐标转换时所采用的椭球参数为WGS-84椭球参数,转换按下述公式进行:
上式中,(X,Y,Z)为地面点在地心直角坐标系中的坐标,(L,B,H)为大地经纬度和大地高,N为椭球的卯酉圈曲率半径,a,b,e分别为地球椭球的长半轴、短半轴和第一偏心率,其中a=6378137,b=6356752.3142;
利用上式分别计算P0、P1和P2的地心直角坐标P0(X0,Y0,Z0)、P1(X1,Y1,Z1)、P2(X2,Y2,Z2),并计算空间直角坐标系向量和向量
步骤4,计算像点p对应地面点入射向量,对于光学卫星影像,入射向量对于合成孔径雷达影像,入射向量
步骤5,计算任意两个影像1和影像2上的同名像点对应的成像交会角:β=arccos[α1·α2/(|α1|·|α2|)],其中α1、α2分别为影像1和影像2中同名像点对应地面点的入射向量;
步骤6,根据影像性质和成像交会角大小,判定同名点是否属于弱交会成像。弱交会成像判断方法为首先给定阈值σ,并根据影像性质进行判定:
(1)当同名点所在影像对同为光学影像,或同为合成孔径雷达影像时:若交会角|β|<σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像;
(2)当同名点所在影像对一景为光学影像,另一景为合成孔径雷达影像时:若交会角满足|β-90|≤σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像。
通过图1可以看出本发明主要包括以下步骤:(1)利用卫星遥感影像附带的RPC参数,高程采用全球SRTM数据,计算同名影像像点坐标对应物方坐标P0(L0,B0,H0);
(2)分别取新的高程值H1=H0+1和H2=H0+1000,利用RPC参数重新计算地面点P1和P2的大地坐标;
(3)将P0、P1和P2的大地坐标转换为地心直角坐标系中的坐标,并计算地心直角坐标系中空间向量和向量
(4)根据向量和向量分别针对光学传感器和雷达传感器,计算地面入射向量;
(5)根据任意两个重叠影像上的同名点对应的入射向量,计算重叠影像同名像点的成像交会角;
(6)根据立体影像的性质和成像交会角的大小,判定同名点是否属于弱交会成像。弱交会成像判断方法为首先给定阈值σ,并根据影像性质进行判定:
1)若立体影像同为光学影像,或同为合成孔径雷达影像时:在成像交会角满足|β|<σ时,认为是弱交会成像,否则认为不属于弱交会成像;
2)若立体影像一景为光学影像,另一景为合成孔径雷达影像时:若成像交会角满足|β-90|≤σ,则认为同名点属于弱交会成像,否则认为不属于弱交会成像。
Claims (3)
1.一种利用RPC探测多源卫星影像同名像点成像交会角方法,包括以下步骤:
(1)利用卫星影像附带的RPC参数,高程采用全球SRTM数据,利用RPC模型计算各影像像点p(R,C)对应地面点P0的大地坐标经度、纬度和高程(L0,B0,H0),R、C为像点p的影像行、列坐标;
(2)依次取新的高程值H1=H0+1和H2=H0+1000,利用RPC模型重新计算得到p点对应地面点P1和P2的大地坐标经度、纬度和高程(L1,B1H1)、(L2,B2,H2);
(3)将P0、P1和P2的大地坐标转换为地心直角坐标系三维坐标,分别得到P0(X0,Y0,Z0)、P1(X1,Y1,Z1)、P2(X2,Y2,Z2),并计算地心直角坐标系中的向量和向量
(4)若p为光学卫星影像上的像点,对应地面点的入射向量若p为合成孔径雷达影像上的像点,入射向量为
(5)计算任意重叠影像1和2上的同名像点对应的成像交会角β=arccos[α1·α2/(|α1|·|α2|)],其中α1、α2分别按步骤(4)计算得到的影像1和2中同名像点对应地面点的入射向量;
(6)根据影像性质和成像交会角大小,判定同名点是否属于弱交会成像,包括下列步骤:
弱交会成像判断方法为首先给定阈值σ,并根据影像性质进行判定:
(1)当同名点所在影像对同为光学影像,或同为合成孔径雷达影像时:若交会角|β|<σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像;
(2)当同名点所在影像对一景为光学影像,另一景为合成孔径雷达影像时:若交会角满足|β-90|≤σ时,认为同名点属于弱交会成像,否则认为不属于弱交会成像。
2.根据权利要求1所述利用RPC探测多源卫星影像同名像点成像交会角方法,其特征在于:所述步骤(1)中,影像坐标为(R,C)的像点p对应地面点大地坐标P0(L0,B0,H0)的计算根据影像RPC参数和SRTM数据,利用RPC模型通过迭代的方法计算得到,包括以下步骤:
(1)给定影像点对应的地面点P0所在区域的近似高程h0作为初值,即令H0=h0;
(2)将像点坐标(R,C)、高程H0和RPC参数代入有理多项式RPC模型,得到下述方程组并求解地面点的归一化大地坐标经度和纬度(Ln,Bn):
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>R</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>R</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,NumR(Ln,Bn,Hn),DenR(Ln,Bn,Hn),NumC(Ln,Bn,Hn),DenC(Ln,Bn,Hn)代表三元三次多项式,且有:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>R</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>7</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>8</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>9</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>10</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>13</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>14</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>15</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>16</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>17</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>18</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>19</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>20</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>R</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>4</mn>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>5</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>6</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>7</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>8</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>9</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>10</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mn>11</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>13</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>14</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>15</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>16</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>17</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>18</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>19</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>20</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>5</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>6</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>7</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>8</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>9</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>10</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>11</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>Z</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>12</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>13</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>14</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>15</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>16</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>17</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>18</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>19</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>20</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>3</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>5</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>6</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>7</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>8</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>9</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>10</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>11</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>12</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>13</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>14</mn>
</msub>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>15</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>16</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>17</mn>
</msub>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>18</mn>
</msub>
<msup>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>19</mn>
</msub>
<msup>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>20</mn>
</msub>
<msup>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中ai,bi,ci,di(i=1….20)为取自影像附带的RPC参数数据集;Rn为影像行坐标R的归一化值,Cn为影像列坐标C的归一化值,Hn为地面点高程H0的归一化值,它们根据下式计算:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mo>(</mo>
<mi>R</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Roff、Coff、Hoff为像点影像行向坐标、列向坐标和地面点高程的平移参数,Rscale、Cscale、Hscale为像点影像行向坐标、列向坐标和地面点高程的缩放系数,它们均取自影像附带的RPC参数数据集;
(3)在得到像点的归一化坐标(Ln,Bn)后,根据下式求解地面点的大地坐标(L0,B0):
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mi>n</mi>
</msub>
<msub>
<mi>B</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
上式中,Loff、Boff为地面点大地坐标经度和纬度的平移参数,Lscale、Bscale为地面点大地坐标经度和纬度的缩放系数,它们均取自影像附带的RPC参数;
(4)将大地坐标(L0,B0,H0)转换到横轴墨卡托地图投影坐标系中,得到坐标(x0,y0),并从SRTM中内插出该投影坐标对应的高程值h,给定阈值10,若|h-H0|<10,则(L0,B0,H0)即为P0点的大地坐标值,终止迭代,否则,取H0=h,重回步骤(2)。
3.根据权利要求1所述的利用RPC探测多源卫星影像同名像点成像交会角方法,其特征在于:所述步骤(3)中,P0、P1和P2的大地坐标转换到地心直角坐标系中,坐标转换时所采用的椭球参数为WGS-84椭球参数,转换按下述公式进行:
<mrow>
<mi>e</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
</mrow>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mfrac>
</msqrt>
</mrow>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mi>a</mi>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>B</mi>
</mrow>
</msqrt>
</mfrac>
</mrow>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Y</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Z</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>(</mo>
<mi>N</mi>
<mo>+</mo>
<mi>H</mi>
<mo>)</mo>
<mi>cos</mi>
<mi> </mi>
<mi>B</mi>
<mi> </mi>
<mi>cos</mi>
<mi> </mi>
<mi>L</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>(</mo>
<mi>N</mi>
<mo>+</mo>
<mi>H</mi>
<mo>)</mo>
<mi>cos</mi>
<mi> </mi>
<mi>B</mi>
<mi> </mi>
<mi>sin</mi>
<mi> </mi>
<mi>L</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>(</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>H</mi>
<mo>)</mo>
<mi>sin</mi>
<mi> </mi>
<mi>B</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
上式中,(X,Y,Z)为地面点在地心直角坐标系中的坐标,(L,B,H)为大地坐标的经度、纬度和高程,N为地球椭球的卯酉圈曲率半径,a,b,e分别为地球椭球的长半轴、短半轴和第一偏心率,其中a=6378137,b=6356752.3142。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510940824.1A CN105466400B (zh) | 2015-12-16 | 2015-12-16 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510940824.1A CN105466400B (zh) | 2015-12-16 | 2015-12-16 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105466400A CN105466400A (zh) | 2016-04-06 |
CN105466400B true CN105466400B (zh) | 2018-05-04 |
Family
ID=55604360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510940824.1A Active CN105466400B (zh) | 2015-12-16 | 2015-12-16 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105466400B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646562A (zh) * | 2016-09-09 | 2017-05-10 | 华东师范大学 | 高精度三维实景室内外一体化定位方法及装置 |
CN107192376B (zh) * | 2017-04-28 | 2019-05-24 | 北京航空航天大学 | 基于帧间连续性的无人机多帧图像目标定位校正方法 |
CN108761453B (zh) * | 2018-03-29 | 2021-02-05 | 航天东方红卫星有限公司 | 一种光学卫星与sar卫星图像融合的成像视角优化方法 |
CN109612439B (zh) * | 2018-12-13 | 2020-04-28 | 同济大学 | 基于有理函数模型的立体影像交会角和基线长度估计方法 |
CN109919835B (zh) * | 2019-03-20 | 2022-07-26 | 湖北省电力勘测设计院有限公司 | 基于多源卫星遥感影像联合平差的境外电力选线方法 |
CN111144350B (zh) * | 2019-12-30 | 2022-02-11 | 长光卫星技术有限公司 | 一种基于参考底图的遥感影像定位精度评价方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390102B (zh) * | 2013-07-16 | 2018-11-30 | 中交第二公路勘察设计研究院有限公司 | 一种卫星图像立体交会角的计算方法 |
CN105091906B (zh) * | 2015-06-30 | 2018-03-02 | 武汉大学 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
-
2015
- 2015-12-16 CN CN201510940824.1A patent/CN105466400B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105466400A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105466400B (zh) | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 | |
CN103558619B (zh) | 获取航拍像片的外方位元素的方法 | |
CN103557841B (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
KR100912715B1 (ko) | 이종 센서 통합 모델링에 의한 수치 사진 측량 방법 및장치 | |
CN103345737B (zh) | 一种基于误差补偿的uav高分辨率影像几何校正方法 | |
CN107316325A (zh) | 一种基于图像配准的机载激光点云与影像配准融合方法 | |
KR101346323B1 (ko) | 지상기준점과 부가매개변수를 이용한 비측량용 디지털 카메라 촬영 시스템의 자체검정 방법 | |
CN104363438B (zh) | 全景立体影像制作方法 | |
KR100529401B1 (ko) | 합성 개구 레이더 영상을 이용한 수치표고자료 제작 장치및 그 방법 | |
CN104764445B (zh) | 水下物点坐标确定方法及装置 | |
CN105716581B (zh) | 基于双介质摄影技术的水下物点坐标确定方法和装置 | |
CN105783875A (zh) | 一种集成非扫描激光雷达的航空摄影测量系统及方法 | |
CN113284171B (zh) | 一种基于卫星遥感立体成像对的植被高度分析方法及系统 | |
CN106408601A (zh) | 一种基于gps的双目融合定位方法及装置 | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
CN102721410B (zh) | 一种基于gps/imu定位定向技术的海岛空中三角测量方法 | |
Gonçalves | Automatic orientation and mosaicking of archived aerial photography using structure from motion | |
CN116625354A (zh) | 一种基于多源测绘数据的高精度地形图生成方法及系统 | |
CN108253942B (zh) | 一种提高倾斜摄影测量空三质量的方法 | |
RU2571300C2 (ru) | Способ дистанционного определения абсолютного азимута целевой точки | |
CN113340272B (zh) | 一种基于无人机微群的地面目标实时定位方法 | |
CN107063191B (zh) | 一种摄影测量区域网整体相对定向的方法 | |
CN114396920A (zh) | 基于倾斜摄影的地形图空间要素获取方法 | |
Ye et al. | Accuracy analysis of real-time object positioning without gcp for images from UAV oblique ultra-long focal small view field whiskbroom camera system | |
CN114140397A (zh) | 全数字地面成像法重力近区地形改正方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |