CN100481133C - 基于有理多项式成像模型的正变换方法 - Google Patents
基于有理多项式成像模型的正变换方法 Download PDFInfo
- Publication number
- CN100481133C CN100481133C CNB2007100528099A CN200710052809A CN100481133C CN 100481133 C CN100481133 C CN 100481133C CN B2007100528099 A CNB2007100528099 A CN B2007100528099A CN 200710052809 A CN200710052809 A CN 200710052809A CN 100481133 C CN100481133 C CN 100481133C
- Authority
- CN
- China
- Prior art keywords
- lat
- lon
- scale
- latitude
- picpointed coordinate
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种RPC模型正变换方法,该方法将反变换结果引入正变换过程进行迭代变换,其过程包括以下步骤:设置高程、和RPC模型参数,设定地面矩形区域;计算地面矩形区域4角点和中心点的影像像素坐标;计算仿射变换系数;计算量测像点坐标对应的近似经纬度坐标;计算理论像点坐标;计算残差e,如果残差e满足设定的差值范围,输出近似经纬度坐标作为量测像点坐标对应的经纬度坐标,否则调整地面矩形区域,然后返回迭代,直到残差e满足差值范围。本发明的技术方案能够确保在影像供应商提供的有限条件下获得精确的地面点位置,在遥感影像应用领域具有重要意义。
Description
技术领域
本发明属于摄影成像技术领域,特别涉及基于有理多项式成像模型的正变换方法。
背景技术
由于严密成像几何模型具有较高的定位精度,因此一直是摄影测量学的首选,卫星遥感影像处理常常要利用严密成像几何模型。然而,一些高性能的传感器***虽然实现了商业化,但是由于卫星平台和有效载荷等关键技术参数是航天大国的技术机密,高分辨率传感器的核心信息和卫星轨道参数并未公开,如IKONOS等卫星,无法利用严格成像几何模型进行处理;同时传感器成像方式的多样化也对摄影测量软件提出了模型通用化的要求,这样在软件框架设计时就不必为现有的或将来可能会出现的各种类型的传感器一一确定其严格几何模型的形式,而可以统一采用通用模型进行处理,从而大大降低了程序设计的复杂性,更易于软件升级和维护,特别在同时处理多源卫星遥感影像数据时,通用成像几何模型更能显示出其优势。
常见的通用成像几何模型有以下几种:多项式、直接线性变换、仿射变换、有理函数模型等。RPC模型是卫星遥感影像的通用成像几何模型其中的一种,它适用于各类传感器包括最新的航空和航天传感器模型。基于RPC的成像模型并不要求了解传感器的实际特性和成像过程,是一种能获得和严格成像模型近似一致精度的形式简单的广义成像模型。
RPC模型具体定义关系式如下:
其中(X,Y)为正则化的影像坐标,(P,L,H)为正则化的地面点坐标。
NumL(P,L,H)、DenL(P,L,H)、Nums(P,L,H)、Dens(P,L,H)为三次多项式:
NumL(P,L,H)=a1+a2L+a3P+a4H+a5LP+a6LH+a7PH+a8L2+a9P2
+a10H2+a11PLH+a12L3+a13LP2+a14LH2+a15L2P+a16P3+a17PH2
+a18L2H+a19P2H+a20H3
DenL(P,L,H)=b1+b2L+b3P+b4H+b5LP+b6LH+b7PH+b8L2+b9P2
+b10H2+b11PLH+b12L3+b13LP2+b14LH2+b15L2P+b16P3+b17PH2
+b18L2H+b19P2H+b20H3
Nums(P,L,H)=c1+c2L+c3P+c4H+c5LP+c6LH+c7PH+c8L2+c9P2
+c10H2+c11PLH+c12L3+c13LP2+c14LH2+c15L2P+c16P3+c17PH2
+c18L2H+c19P2H+c20H3
Dens(P,L,H)=d1+d2L+d3P+d4H+d5LP+d6LH+d7PH+d8L2+d9P2
+d10H2+d11PLH+d12L3+d13LP2+d14LH2+d15L2P+d16P3+d17PH2
+d18L2H+d19P2H+d20H3
其中,三次多项式的系数a1,...,a20,b1,...,b20,c1,...,c20是RPC文件中提供的模型参数,b1和d1通常为1,l,s分别为影像列数值和行数值。
所谓正则化,是一项处理无限大、发散以及一些不合理表示式的方法,其方法透过引入一项辅助性的概念——正则化因子。在RPC模型中运用这种方法对地面点和影像点的坐标进行处理。
正则化地面点坐标定义公式为
其中,正则化因子LAT_OFF、LAT_SCALE、LONG_OFF、LONG_SCALE、HEIGHT_OFF、HEIGHT_SCALE是RPC文件中包含的地面点坐标正则化模型参数。Latitude表示经度、Longitude表示纬度、Height表示高程,此三项即可代表地面点的空间坐标。注:某点沿铅垂线方向到大地水准面的距离,通常称为绝对高程或海拔,本技术领域简称高程。
正则化影像坐标定义公式为
其中,正则化因子SAMP_OFF、SAMP_SCALE、LINE_OFF、LINE_SCALE也是RPC文件中包含的影像坐标正则化模型参数,sample代表影像列坐标,其数值即为s,line代表影像行坐标,其数值即为l。
在基于RPC模型的几何应用中首先需要基于RPC模型的正反变换。RPC正变换是指根据RPC模型的关系式,从地面控制点坐标到遥感影像平面坐标的变换过程。RPC反变换则是根据RPC模型的关系式,从原始遥感影像和影像中每个像素对应的高程数据到一定地图投影的变换过程。
但是在现有技术中,基于RPC模型的正变换都是利用最小二乘平差理论实现正变换,首先假定在地面的初始位置,然后迭代求解,由于假定的初始位置不一定准确,并且收敛的控制条件不一定是是最优的,有可能导致收敛位置的不准确。因此在RPC模型的实际应用中急切需要应用先进的科技手段,提高RPC模型正变换的可靠性,并提供一种新的RPC模型的正变换的方法以使基于RPC模型的几何处理获得更高的精度和可靠性。
发明内容
本发明目的在于解决现有技术不足,提出了一种全新的基于RPC模型反变换的正变换方法,解决了目前RPC模型正变换迭代不一定收敛的RPC模型应用中的瓶颈问题。
本发明的技术方案是:将有理多项式成像模型反变换结果引入有理多项式成像模型正变换过程进行迭代变换,即采用计算机进行以下步骤,
步骤1,设置高程数值Dhei和有理多项式成像模型参数,设定中心点经纬度坐标为(Dlat_off,Dlon_off)、边长为Dlat_scale和Dlon_scale的地面矩形区域,其中4个角点的经纬度坐标为
(Dlat_off+Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off+Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0);
步骤2,利用有理多项式成像模型的反变换模型计算地面矩形区域4个角点和中心点对应的影像像点坐标;
步骤3,利用上述5个点的经纬度坐标和影像像点坐标,根据仿射变换关系式 利用最小二乘原理计算影像的近似仿射变换参数f0,f1,f2,g0,g1,g2,x、y分别表示正则化的影像像点坐标的列数值和行数值,Dlat、Dlon分别表示正则化的地面点坐标的纬度数值和经度数值;
步骤6,计算残差e,计算公式为e=(|xa-xb|)2+(|ya-yb|)2,如果残差e满足设定的差值范围,直接输出作为量测像点坐标(xa,ya)对应的经纬度坐标否则以为中心点经纬度坐标缩小地面矩形区域范围,然后返回步骤2迭代运行,直到残差e满足设定的差值范围。
其中像点坐标(x1,y1)为地面点(Dhei)所对应的像点坐标,像点坐标(x2,y2)为地面点 所对应的像点坐标,θ为缩小取值系数;
步骤6.2,根据像素大小PS调整地面矩形区域范围,地面矩形区域4个角点的经纬度坐标为
其中es=|xa-xb|,表示列方向的残差,el=|ya-yb|,表示行方向的残差。
而且,缩小取值系数θ的取值为0.001~0.05。
本发明提出了一种全新的基于RPC模型反变换的正变换方法,将反变换结果引入进行迭代,克服了现有技术收敛难度,提高RPC模型正变换的可靠性,使基于RPC模型的几何处理获得更高的精度和可靠性。并且本发明的实现可在提供的现有模型参数下完成,实现方便。本发明提供的方案将可在卫星遥感影像处理的各方面工作中均到重要作用,例如影像的几何纠正、影像匹配和核线排列等方面。
附图说明
图1是本发明流程图。
具体实施方式
本发明将RPC模型反变换结果引入RPC模型正变换过程进行迭代变换,具体实施时可采用计算机程序完成,流程可参见附图。本发明的RPC模型正变化过程包括以下步骤,
步骤1,设置高程数值Dhei和RPC模型参数,设定中心点经纬度坐标为(Dlat_off,Dlon_off)、边长为Dlat_scale和Dlon_scale的地面矩形区域,其中4个角点的经纬度坐标为
(Dlat_off+Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off+Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0)。
经纬度坐标的含义:在大地坐标系中,某点大地经度,就是通过该点的子午面与起始子午面所构成的二面角。由起始子午面起算,向东为正,叫东经,向西为负,叫西经。这点的法线与赤道面的夹角,叫做大地纬度,由赤道面起算,向北为正,叫北纬,向南为负,叫南纬。由大地经度和大地纬度所表示的地面点坐标就成为经纬度坐标。经纬度坐标加上高程即可标注地面点的坐标。为了满足下列步骤的计算条件,本发明先设置高程数值Dhei,RPC模型参数则由影像供应商提供,具体实施时输入即可。设定的地面矩形区域,需要迭代时会进行调整。
步骤2,利用RPC反变换模型计算地面矩形区域4个角点和中心点对应的影像像点坐标。
步骤3,利用上述5个点的经纬度坐标和影像像点坐标,根据仿射变换关系式 利用最小二乘原理计算影像的近似仿射变换参数f0,f1,f2,g0,g1,g2,x、y分别表示正则化的影像像点坐标的列数值和行数值,Dlat、Dlon分别表示正则化的地面点坐标的纬度数值和经度数值;
仿射变换是平移、旋转、缩放等运算的组合,仿射变换关系式的系数即为仿射变换参数。
步骤5,根据RPC模型的反变换模型计算地面点坐标(Dhei)对应的理论像点坐标(xb,yb);
计算理论像点坐标(xa,ya),是用于将反变换结果引入正变换过程以便修正。
RPC模型的反变换模型即公式
在步骤4中获得的近似经纬度坐标结合高程数值Dhei,即可分别作为P、L、H带入上式求解X、Y值,获得对应的理论像点坐标(xb,yb)。步骤6,计算残差e,计算公式为e=(|xa-xb|)2+(|ya-yb|)2,如果残差e满足设定的差值范围,直接输出作为量测像点坐标(xa,ya)对应的经纬度坐标否则以为中心点经纬度坐标缩小地面矩形区域范围,然后返回步骤2迭代运行,直到残差e满足设定的差值范围。
本发明实施时,可根据精度需要设定差值范围,例如设定为残差e小于0.01像素,残差是指变量的真值与观测值之差。es=|xa-xb|,表示列方向的残差,el=|ya-yb|,表示行方向的残差,e=(|xa-xb|)2+(|ya-yb|)2,表示距离上的残差。
为了提高每次迭代运算结果精度,保证收敛过程,本发明还提供了迭代时调整地面矩形区域的具体方法:
其中像点坐标(x1,y1)为地面点(Dhei)所对应的像点坐标,并且数值为(xb,yb);像点坐标(x2,y2)为地面点 所对应的像点坐标,θ为缩小取值系数,sqrt表示开平方;缩小取值系数θ的取值范围可根据具体情况选取,例如0.001~0.05,本发明实施例选用0.01,则PS=0.01×Dlat_scale/sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))
具体实施时,对影像的像素大小PS的计算并不限于此种方式,例如调整纬度数值,计算 对应的像点坐标;如果改为调整经度,相应可改成计算 对应的像点坐标;或者经度纬度同时调整也可。相应的,公式中分子也需要相应的改变。
步骤6.2,根据像素大小PS调整地面矩形区域范围,地面矩形区域中心点和4个角点的经纬度坐标为
Claims (2)
1.基于有理多项式成像模型的正变换方法,其特征在于:将有理多项式成像模型反变换结果引入有理多项式成像模型正变换过程进行迭代变换,即采用计算机进行以下步骤,步骤1,设置高程数值Dhei和有理多项式成像模型参数,设定中心点经纬度坐标为(Dlat_off,Dlon_off)、边长为Dlat_scale和Dlon_scale的地面矩形区域,其中4个角点的经纬度坐标为
(Dlat_off+Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off+Dlon_scale/2.0);
(Dlat_off+Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0);
(Dlat_off-Dlat_scale/2.0,Dlon_off-Dlon_scale/2.0);
步骤2,利用有理多项式成像模型的反变换模型计算地面矩形区域4个角点和中心点对应的影像像点坐标;
步骤3,利用上述5个点的经纬度坐标和影像像点坐标,根据仿射变换关系式 利用最小二乘原理计算影像的近似仿射变换参数f0,f1,f2,g0,g1,g2,x、y分别表示正则化的影像像点坐标的列数值和行数值,Dlat、Dlon分别表示正则化的地面点坐标的纬度数值和经度数值;
步骤6,计算残差e,计算公式为e=(|xa-xb|)2+(|ya-yb|)2如果残差e满足设定的差值范围,直接输出作为量测像点坐标(xa,ya)对应的经纬度坐标否则以为中心点经纬度坐标缩小地面矩形区域范围,然后返回步骤2迭代运行,直到残差e满足设定的差值范围,
以为中心点经纬度坐标缩小地面矩形区域范围的具体步骤如下,步骤6.1,根据公式PS=θ×Dlat_scale/sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))计算 位置附近的影像的像素大小PS,
步骤6.2,根据像素大小PS调整地面矩形区域范围,地面矩形区域4个角点的经纬度坐标为
其中es=|xa-xb|,表示列方向的残差,el=|ya-yb|,表示行方向的残差。
2.根据权利要求1所述基于有理多项式成像模型的正变换方法,其特征在于:缩小取值系数θ的取值范围为0.001~0.05。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100528099A CN100481133C (zh) | 2007-07-23 | 2007-07-23 | 基于有理多项式成像模型的正变换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007100528099A CN100481133C (zh) | 2007-07-23 | 2007-07-23 | 基于有理多项式成像模型的正变换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101101674A CN101101674A (zh) | 2008-01-09 |
CN100481133C true CN100481133C (zh) | 2009-04-22 |
Family
ID=39035939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007100528099A Expired - Fee Related CN100481133C (zh) | 2007-07-23 | 2007-07-23 | 基于有理多项式成像模型的正变换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100481133C (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598785B (zh) * | 2009-06-30 | 2011-12-28 | 武汉大学 | 卫星影像各级产品的有理函数成像模型生成方法 |
CN101598786B (zh) * | 2009-06-30 | 2012-07-04 | 武汉大学 | 卫星影像二/三级产品的有理函数成像模型生成方法 |
CN102168972B (zh) * | 2010-12-15 | 2012-10-03 | 中国资源卫星应用中心 | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 |
CN102722250A (zh) * | 2012-06-07 | 2012-10-10 | 何潇 | 一种影像控制点交互编辑的方法及*** |
CN107316280B (zh) * | 2017-05-24 | 2019-10-11 | 三亚中科遥感研究所 | 离岛卫星影像rpc模型高精度几何定位方法 |
CN111121650B (zh) * | 2019-12-25 | 2021-08-31 | 中国民航科学技术研究院 | 一种利用dem辅助卫星影像测算建筑物高度的方法及*** |
CN113378104A (zh) * | 2021-06-11 | 2021-09-10 | 中国人民解放军战略支援部队信息工程大学 | 改进的自适应谱修正迭代法的rpc参数求解方法 |
CN117151991B (zh) * | 2023-07-06 | 2024-05-31 | 北京数慧时空信息技术有限公司 | 一种遥感影像rpc修复方法 |
-
2007
- 2007-07-23 CN CNB2007100528099A patent/CN100481133C/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
The Cubic Rational Polynomial Camera Model. Richard I. Hartley,Tushar Saxena.http://www.cs.albany.edu/~saxena/Papers/cubic.pdf. 2000 |
The Cubic Rational Polynomial Camera Model. Richard I. Hartley,Tushar Saxena.http://www.cs.albany.edu/~saxena/Papers/cubic.pdf. 2000 * |
无需初值的RPC模型参数求解算法研究. 秦绪文等.国土资源遥感,第4期. 2005 |
无需初值的RPC模型参数求解算法研究. 秦绪文等.国土资源遥感,第4期. 2005 * |
Also Published As
Publication number | Publication date |
---|---|
CN101101674A (zh) | 2008-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100481133C (zh) | 基于有理多项式成像模型的正变换方法 | |
CN103759714B (zh) | 一种三线阵卫星影像区域网平差方法 | |
Luhmann et al. | Close range photogrammetry | |
CN101598785B (zh) | 卫星影像各级产品的有理函数成像模型生成方法 | |
CN102866397B (zh) | 一种多源异构遥感影像联合定位方法 | |
CN101545776B (zh) | 基于数字地图的数码像片方位元素获取方法 | |
CN104457710B (zh) | 一种基于非量测数码相机的航空数字摄影测量方法 | |
CN102519433B (zh) | 一种利用rpc反演星载线阵传感器几何定标参数方法 | |
CN103823981A (zh) | 一种数字高程模型辅助的卫星影像区域网平差方法 | |
CN104406607A (zh) | 一种多视场复合光学敏感器的标定装置及方法 | |
CN103310487B (zh) | 一种基于时间变量的通用成像几何模型生成方法 | |
CN102778224B (zh) | 一种基于极坐标参数化的航空摄影测量光束法平差的方法 | |
CN103557841A (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
CN104820984B (zh) | 一种卫星遥感三线阵立体影像处理***及方法 | |
CN1996389A (zh) | 基于共线特征点的摄像机畸变快速校正方法 | |
CN103390102B (zh) | 一种卫星图像立体交会角的计算方法 | |
CN104123695A (zh) | 一种实现坐标转换方法 | |
CN101598786B (zh) | 卫星影像二/三级产品的有理函数成像模型生成方法 | |
CN106405581B (zh) | 卫星结构变形引起的多类型载荷间耦合指向精度评估方法 | |
CN101907705B (zh) | 通用的多源遥感影像几何校正模型联合平差方法 | |
CN104112078B (zh) | 一种自适应几何关系强弱的有理函数模型区域网平差方法 | |
CN101441075B (zh) | 采用遥感地图工具的地图投影方法 | |
CN104754323A (zh) | 摄像头光轴检测设备的标定方法 | |
CN103344252B (zh) | 一种航空高光谱成像***定位误差分析方法 | |
CN102346252A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090422 Termination date: 20170723 |
|
CF01 | Termination of patent right due to non-payment of annual fee |