CN103390102A - 一种卫星图像立体交会角的计算方法 - Google Patents
一种卫星图像立体交会角的计算方法 Download PDFInfo
- Publication number
- CN103390102A CN103390102A CN2013102969994A CN201310296999A CN103390102A CN 103390102 A CN103390102 A CN 103390102A CN 2013102969994 A CN2013102969994 A CN 2013102969994A CN 201310296999 A CN201310296999 A CN 201310296999A CN 103390102 A CN103390102 A CN 103390102A
- Authority
- CN
- China
- Prior art keywords
- partiald
- image
- coordinate
- point
- lon
- 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.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明属于空间对地观测技术领域,具体公开了一种卫星图像立体交会角的计算方法,其步骤:A、通过人工量测或者卫星图像自动匹配的方法获取卫星图像立体像对同名点的像点坐标。B、利用卫星图像同名点坐标,进行基于RFM模型的区域网平差方法,实现RPC参数的优化。C、通过投影轨迹法,计算同名点投影到不同高程投影面上的物方点所对应的空间直角坐标,完成同名光线的空间重建。D、利用重建同名光线的方向矢量,计算同名光线立体交会角。本发明适用于任何相互重叠卫星图像构成的立体模型立体交会角的计算,无需地面控制点,可用于立体模型测图精度评估,具有显著的经济和社会效益。
Description
技术领域
本发明属于空间对地观测技术领域,具体涉及一种卫星图像立体交会角的计算方法,该方法可以快速计算出相互重叠卫星图像所构成的“广义立体模型”的交会角,为卫星图像立体模型精度估算和选择利用提供了技术依据。
背景技术
近年来,随着对地观测技术的飞速发展,卫星遥感图像在时间分辨率、空间分辨率和光谱分辨率方面均得到了革命性的提升,尤其是空间分辨率。目前,高分辨率遥感卫星不断涌现,光学遥感卫星有Ikonos,Quickbird,GeoEye和WorldView等,SAR遥感卫星有TerraSAR-X,Cosmo-Skymed和Radarsat-2等,这使得空间三维信息获取手段从航空平台上升到航天平台。通常,商业立体卫星图像对拍摄方位、立体交会角都有严格的要求,能够保证所构成立体模型的测图精度,但是其价格比单景卫星图像昂贵许多。此外,在实际工程应用中,受云层覆盖等影响,有时无法及时获取高分辨率卫星立体图像,这将极大地延误公路勘察设计工期。从技术角度上说,任何从不同视角拍摄的相互重叠的两景卫星图像均可构成立体模型,又称为“广义立体模型”。利用广义立体模型获取高精度的地形数据,不仅可以降低数据获取的成本,而且可以降低对立体卫星图像的依赖,具有重要的实际意义。
立体交会角,是评估卫星图像立体模型精度的重要指标。对于航空影像,可直接利用基高比来估算立体交会角,但对于光学卫星图像,一般为多中心投影成像方式,且卫星图像公司往往不提供构建严格几何模型所需的各种轨道定向参数,因此,卫星图像立体交会角往往无法得到。
为此,本申请的发明人拟提出一种卫星图像立体交会角的计算方法,无需地面控制点,仅需利用卫星图像及其自带的有理函数模型(Rational Function Model,RFM),将产生显著的经济和社会效益。
发明内容
本发明的目的在于针对实际工程应用中利用单景卫星图像构成的“广义立体模型”空间信息获取的技术问题,提出了一种卫星图像立体交会角的计算方法。该方法直接利用卫星图像及其自带的RPC参数,无需任何辅助信息,可免除布设地面控制点工作,精确快速地计算出卫星图像的立体交会角,为“广义立体模型”测图精度估算提供有力依据,有效保证采集数据的质量。
为了达到上述目的,本发明采取了如下技术措施:
一种卫星图像立体交会角的计算方法,其步骤如下:
1、卫星图像同名点图像坐标获取:采用人工量测或者卫星图像自动匹配的方法得到;
2、卫星图像RPC参数优化:利用卫星图像同名点,进行基于RFM模型的区域网平差, 实现对卫星图像RPC参数的优化;
3、卫星图像同名光线空间重建:通过投影轨迹法,建立同名点在物方空间中同名光线的直线方程;
4、卫星图像立体交会角计算:利用重建的同名光线直线方程,计算出相应的立体交会角。
与现有的技术相比,本发明的显著优点和效果主要表现在:
1、直接利用卫星图像及其自带的支持文件信息,无需布设地面控制点,方法简单实用;
2、可在立体模型构建前计算出立体交会角,能对数据质量做出预估,方便卫星图像的合理选择和有效利用。
本发明创新性地基于RFM模型提出卫星图像“广义立体模型”立体交会角计算方法,可缩短数据获取周期、保障数据质量。本发明的技术方案已经成功应用于西藏1000多公里的公路勘察设计项目中,成功剔除测图精度不合要求的卫星图像立体模型,缩短数据获取周期1/3,展现出重大的工程应用价值。
附图说明
图1为本发明的一种卫星图像立体交会角的计算方法的流程示意图。
图2为基于RFM模型计算同名光线立体交会角的示意图。
具体实施方式
下面结合附图对本发明方法做进一步的详细描述。
实施例1:
一种卫星图像立体交会角的计算方法,如图1所示,各步骤详细阐述如下:
第1步:卫星图像同名点图像坐标获取,可以采用两种方法:
1)人工量测:人工判读出卫星图像间具有相同影像特征的点,并利用量测工具准确获取同名点真实图像坐标(sample,line)。
或2)自动匹配:直接利用卫星图像自动匹配软件,如ERDAS Imagine软件的Leica Photogrammetry Suite模块,获取同名点图像坐标,并进行检查,删除错误匹配点,最终获取卫星图像同名点真实图像坐标(sample,line)。
第2步:卫星图像RPC参数优化
利用所获取的同名点真实图像坐标,通过RFM模型区域网平差的方法,对RPC参数进行优化,从而实现同名光线对对相交。
对于卫星图像有理函数模型,图像坐标(x,y)表示为以对应地面点大地坐标(Lat,Lon,Height)为自变量的多项式的比值,即:
式中,为图像坐标(x,y)的规格化坐标,P,L,H为地面点大地坐标(Lat,Lon,Height)的规格化坐标,pi(P,L,H)(i=1,2,3,4)一般为三次有理函数多项式,其系数称为RPC(Rational Polynomial Coefficients)参数。
和P,L,H分别通过对图像坐标(x,y)和地面点大地坐标(Lat,Lon,Height)进行平移和缩放后得到,其变换关系为:
式中,LAT_SACLE,LON_SCALE,HEIGHT_SCALE表示物方大地坐标规格化缩放参数,LAT_OFF,LON_OFF,HEIGHT_OFF表示物方大地坐标规格化平移参数,SAMP_SCALE,LINE_SCALE表示像方坐标规格化缩放参数,SAMP_OFF,LINE_OFF表示像方坐标规格化平移参数。
由于卫星图像自带的RPC参数存在***误差,利用公式(1)、(2)计算得到的图像坐标(x,y)和真实图像坐标(sample,line)存在差异,采用像方仿射变换模型对其进行校正,从而建立地面点大地坐标(Lat,Lon,Height)和真实图像坐标(sample,line)间精确的空间几何关系,其表达式为:
令 对每个同名点的像点坐标,以其所在卫星图像对应的仿射变换系数和该点对应的物方坐标为未知数,对公式(3)线性化,可以建立基于RFM模型区域网平差的误差方程:
V=At+Bx-L,P (4)
式中,
x=[dLat dLon dHeight]T;
P—权矩阵。
式中,为Fx对地面点大地坐标的偏导数,为Fy对地面点大地坐标的偏导数,dLat、dLon、dHeight为同名点对应地面点大地坐标的改正数,(i=0,1,2)为Fx对图像仿射变换系数的偏导数,(i=0,1,2)为Fy对图像仿射变换系数的偏导数,(i=0,1,2)为图像仿射变换系数的改正数, 为未知参数取近似值时Fx、Fy的计算值,vx、vy为Fx、Fy虚拟观测值的改正数。
利用最小二乘平差原理,计算出每张卫星图像的像方仿射变换系数,从而完成RPC参数的优化。
第3步:卫星图像同名光线空间重建
采用投影轨迹法实现,如图2所示。图中,点(x1,y1)和点(x2,y2)是一个同名点对,首先,以测区平均高程为中心建立不同的高程平面其中△Z的取值可任意,本实施例中取为1000米。在优化PRC参数后,通过点(x1,y1)和点(x2,y2)可以确定两条同名光线L1、L2,同名光线L1、L2分别与两个高程平面相交于点(X1,Y1,Z1)、(X2,Y2,Z2)和(X3,Y3,Z3)、(X4,Y4,Z4)。通过这四个交点,可实现同名光线的空间重建。
四个交点坐标的计算方法如下,以交点(X1,Y1,Z1)的坐标计算为例:
V=Bx-L (6)
式中,Fx0,Fy0是指未知数Lat、Lon取初始值时的Fx,Fy的计算值,dLat,dLon分别代表Lat、Lon的改正数。通过迭代计算,可计算得到大地经纬度坐标(Lat,Lon)。
采用与交点(X1,Y1,Z1)坐标计算相同的方法,可以计算得到(X2,Y2,Z2)、(X3,Y3,Z3)和(X4,Y4,Z4)。利用坐标(X1,Y1,Z1)和(X2,Y2,Z2)、(X3,Y3,Z3)和(X4,Y4,Z4),可以分别完成同名光线L1、L2的空间重建。
同名光线L1的空间直线方程为:
同名光线L2的空间直线方程为:
第4步:卫星图像立体交会角计算
Claims (2)
1.一种卫星图像立体交会角的计算方法,其步骤如下:
(1)、卫星图像同名点图像坐标获取:采用人工量测或者卫星图像自动匹配的方法得到;
(2)、卫星图像RPC参数优化:利用卫星图像同名点,进行基于RFM模型的区域网平差,实现对卫星图像RPC参数的优化;
(3)、卫星图像同名光线空间重建:通过投影轨迹法,建立同名点在物方空间中同名光线的直线方程;
(4)、卫星图像立体交会角计算:利用重建的同名光线直线方程,计算出相应的立体交会角。
2.根据权利要求1所述的计算方法,其特征在于,步骤如下:
第1步:卫星图像同名点图像坐标获取,采用如下两种方法之一:
1)人工量测:人工判读出卫星图像间具有相同影像特征的点,并利用量测工具准确获取同名点真实图像坐标(sample,line);
或2)自动匹配:利用卫星图像自动匹配软件获取同名点图像坐标,并进行检查,删除错误匹配点,最终获取卫星图像同名点真实图像坐标(sample,line);
第2步:卫星图像RPC参数优化
利用所获取的同名点真实图像坐标,通过RFM模型区域网平差的方法,对RPC参数进行优化,从而实现同名光线对对相交;
对于卫星图像有理函数模型,图像坐标(x,y)表示为以对应地面点大地坐标(Lat,Lon,Height)为自变量的多项式的比值,即:
式中,为图像坐标(x,y)的规格化坐标,P,L,H为地面点大地坐标(Lat,Lon,Height)的规格化坐标,pi(P,L,H)(i=1,2,3,4)一般为三次有理函数多项式,其系数称为RPC参数;
式中,LAT_SACLE,LON_SCALE,HEIGHT_SCALE表示物方大地坐标规格化缩放参数,LAT_OFF,LON_OFF,HEIGHT_OFF表示物方大地坐标规格化平移参数,SAMP_SCALE,LINE_SCALE表示像方坐标规格化缩放参数,SAMP_OFF,LINE_OFF表示像方坐标规格化平移参数;
由于卫星图像自带的RPC参数存在***误差,利用公式(1)、(2)计算得到的图像坐标(x,y)和真实图像坐标(sample,line)存在差异,采用像方仿射变换模型对其进行校正,从而建立地面点大地坐标(Lat,Lon,Height)和真实图像坐标(sample,line)间精确的空间几何关系,其表达式为:
令 对每个同名点的像点坐标,以其所在卫星图像对应的仿射变换系数和该点对应的物方坐标为未知数,对公式(3)线性化,可以建立基于RFM模型区域网平差的误差方程:
V=At+Bx-L,P (4)
式中,
x=[dLat dLon aHeight]T;
P—权矩阵;
式中,为Fx对地面点大地坐标的偏导数,为Fy对地面点大地坐标的偏导数,dLat、dLon、dHeight为同名点对应地面点大地坐标的改正数,(i=0,1,2)为Fx对图像仿射变换系数的偏导数,(i=0,1,2)为Fy对图像仿射变换系数的偏导数,(i=0,1,2)为图像仿射变换系数的改正数,为未知参数取近似值时Fx、Fy的计算值,vx、vy为Fx、Fy虚拟观测值的改正数;
利用最小二乘平差原理,计算出每张卫星图像的像方仿射变换系数,从而完成RPC参数的优化;
第3步:卫星图像同名光线空间重建
采用投影轨迹法实现,点(x1,y1)和点(x2,y2)是一个同名点对,首先,以测区平均高程为中心建立不同的高程平面其中△Z的取值可任意;在优化PRC参数后,通过点(x1,y1)和点(x2,y2)可以确定两条同名光线L1、L2,同名光线L1、L2分别与两个高程平面相交于点(X1,Y1,Z1)、(X2,Y2,Z2)和(X3,Y3,Z3)、(X4,Y4,Z4);通过这四个交点,实现同名光线的空间重建;
四个交点坐标的计算方法如下,以交点(X1,Y1,Z1)的坐标计算为例:
V=Bx-L (6)
式中,Fx0,Fy0是指未知数Lat、Lon取初始值时的Fx,Fy的计算值,dLat,dLon分别代表Lat、Lon的改正数。通过迭代计算,可计算得到大地经纬度坐标(Lat,Lon);
利用大地测量学中大地坐标系与空间直角坐标系的关系式,可以获得点(x1,y1)投影到高程面的物方点在空间直角坐标系中的物方坐标(X1,Y1,Z1);
采用与交点(X1,Y1,Z1)坐标计算相同的方法,计算得到(X2,Y2,Z2)、(X3,Y3,Z3)和(X4,Y4,Z4),利用坐标(X1,Y1,Z1)和(X2,Y2,Z2)、(X3,Y3,Z3)和(X4,Y4,Z4),分别完成同名光线L1、L2的空间重建。
同名光线L1的空间直线方程为:
同名光线L2的空间直线方程为:
第4步:卫星图像立体交会角计算
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310296999.4A CN103390102B (zh) | 2013-07-16 | 2013-07-16 | 一种卫星图像立体交会角的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310296999.4A CN103390102B (zh) | 2013-07-16 | 2013-07-16 | 一种卫星图像立体交会角的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103390102A true CN103390102A (zh) | 2013-11-13 |
CN103390102B CN103390102B (zh) | 2018-11-30 |
Family
ID=49534371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310296999.4A Active CN103390102B (zh) | 2013-07-16 | 2013-07-16 | 一种卫星图像立体交会角的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103390102B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794490A (zh) * | 2015-04-28 | 2015-07-22 | 中测新图(北京)遥感技术有限责任公司 | 航空多视影像的倾斜影像同名点获取方法和装置 |
CN104978761A (zh) * | 2015-06-01 | 2015-10-14 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN105387846A (zh) * | 2015-10-26 | 2016-03-09 | 中国农业大学 | 卫星影像的正射校正方法和*** |
CN105466400A (zh) * | 2015-12-16 | 2016-04-06 | 贵州天地通科技有限公司 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
CN105740624A (zh) * | 2016-01-30 | 2016-07-06 | 武汉大学 | 复杂交会条件下区域网平差连接点自适应定权方法及*** |
CN106595600A (zh) * | 2016-12-23 | 2017-04-26 | 航天恒星科技有限公司 | 立体测绘卫星姿态测定***低频误差补偿方法及*** |
CN108388687A (zh) * | 2018-01-05 | 2018-08-10 | 合肥工业大学 | 考虑传感器侧摆的多成像卫星区域覆盖优化方法 |
CN109612439A (zh) * | 2018-12-13 | 2019-04-12 | 同济大学 | 基于有理函数模型的立体影像交会角和基线长度估计方法 |
CN114998397A (zh) * | 2022-05-20 | 2022-09-02 | 中国人民解放军61540部队 | 一种多视角卫星影像立体像对优化选择方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915913A (zh) * | 2010-07-30 | 2010-12-15 | 中交第二公路勘察设计研究院有限公司 | 一种稳健的高分辨率卫星影像连接点自动匹配方法 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
-
2013
- 2013-07-16 CN CN201310296999.4A patent/CN103390102B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915913A (zh) * | 2010-07-30 | 2010-12-15 | 中交第二公路勘察设计研究院有限公司 | 一种稳健的高分辨率卫星影像连接点自动匹配方法 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
Non-Patent Citations (4)
Title |
---|
刑帅等: "SPOT5与ERS-2卫星影像的联合立体定位", 《测绘科学》, vol. 34, no. 1, 20 January 2009 (2009-01-20), pages 46 - 47 * |
吕勇等: "三线阵测绘相机交会角在轨监测可行性分析", 《红外与激光工程》, vol. 41, no. 12, 25 December 2012 (2012-12-25), pages 3391 - 3393 * |
张过等: "卫星遥感影像的区域网平差成图精度", 《测绘科学技术学报》, vol. 23, no. 4, 30 August 2006 (2006-08-30), pages 239 - 241 * |
王伟玺等: "基于RFM的广义立体像对构建方法研究", 《测绘工程》, vol. 16, no. 2, 25 April 2007 (2007-04-25), pages 22 - 26 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794490A (zh) * | 2015-04-28 | 2015-07-22 | 中测新图(北京)遥感技术有限责任公司 | 航空多视影像的倾斜影像同名点获取方法和装置 |
CN104794490B (zh) * | 2015-04-28 | 2018-10-02 | 中测新图(北京)遥感技术有限责任公司 | 航空多视影像的倾斜影像同名点获取方法和装置 |
CN104978761A (zh) * | 2015-06-01 | 2015-10-14 | 北京航空航天大学 | 一种rpc模型校正系数获取方法 |
CN105387846B (zh) * | 2015-10-26 | 2017-11-10 | 中国农业大学 | 卫星影像的正射校正方法和*** |
CN105387846A (zh) * | 2015-10-26 | 2016-03-09 | 中国农业大学 | 卫星影像的正射校正方法和*** |
CN105466400A (zh) * | 2015-12-16 | 2016-04-06 | 贵州天地通科技有限公司 | 一种利用rpc探测多源卫星影像同名像点成像交会角方法 |
CN105740624A (zh) * | 2016-01-30 | 2016-07-06 | 武汉大学 | 复杂交会条件下区域网平差连接点自适应定权方法及*** |
CN105740624B (zh) * | 2016-01-30 | 2018-04-20 | 武汉大学 | 复杂交会条件下区域网平差连接点自适应定权方法及*** |
CN106595600A (zh) * | 2016-12-23 | 2017-04-26 | 航天恒星科技有限公司 | 立体测绘卫星姿态测定***低频误差补偿方法及*** |
CN106595600B (zh) * | 2016-12-23 | 2019-04-30 | 航天恒星科技有限公司 | 立体测绘卫星姿态测定***低频误差补偿方法及*** |
CN108388687A (zh) * | 2018-01-05 | 2018-08-10 | 合肥工业大学 | 考虑传感器侧摆的多成像卫星区域覆盖优化方法 |
CN109612439A (zh) * | 2018-12-13 | 2019-04-12 | 同济大学 | 基于有理函数模型的立体影像交会角和基线长度估计方法 |
CN114998397A (zh) * | 2022-05-20 | 2022-09-02 | 中国人民解放军61540部队 | 一种多视角卫星影像立体像对优化选择方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103390102B (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103390102A (zh) | 一种卫星图像立体交会角的计算方法 | |
CN103345737B (zh) | 一种基于误差补偿的uav高分辨率影像几何校正方法 | |
CN102506824B (zh) | 一种城市低空无人机***生成数字正射影像图的方法 | |
CN101750619B (zh) | 自检校pos直接对地目标定位方法 | |
CN102866397B (zh) | 一种多源异构遥感影像联合定位方法 | |
CN102645209B (zh) | 机载LiDAR点云和高分辨率影像进行空间点的联合定位方法 | |
CN107389029A (zh) | 一种基于多源监测技术融合的地面沉降集成监测方法 | |
CN101876701B (zh) | 一种侧视雷达遥感影像定位方法 | |
CN104567819B (zh) | 一种星载相机全视场偏流角确定与补偿方法 | |
CN109029368B (zh) | 像方补偿的遥感影像/sar影像高精度几何定位后处理方法 | |
CN106960174A (zh) | 高分影像激光雷达高程控制点提取及其辅助定位方法 | |
CN102426025B (zh) | 遥感卫星姿态机动时偏流修正角的仿真分析方法 | |
CN110046563B (zh) | 一种基于无人机点云的输电线路断面高程修正方法 | |
CN105004354A (zh) | 大斜视角下无人机可见光和红外图像目标定位方法 | |
CN105698766B (zh) | 考虑定向参数精度信息的卫星影像rfm模型区域网平差方法 | |
CN104422425A (zh) | 一种不规则外形物体空间姿态动态测量方法 | |
Liu et al. | Accurate mapping method for UAV photogrammetry without ground control points in the map projection frame | |
CN102243299B (zh) | 无人机载sar图像正射纠正装置 | |
CN111597496B (zh) | 基于有理多项式系数单幅星载遥感影像目标高度计算方法 | |
CN110986888A (zh) | 一种航空摄影一体化方法 | |
Zhao et al. | Direct georeferencing of oblique and vertical imagery in different coordinate systems | |
CN105004321B (zh) | 顾及曝光不同步的无人机gps辅助光束法平差方法 | |
CN105783879B (zh) | 一种条带约束下的同轨卫星连续景影像定位方法 | |
Wang et al. | Geometric calibration for the aerial line scanning camera GFXJ | |
Garcia et al. | The Influence of Ground Control Points Configuration and Camera Calibration for Dtm and Orthomosaic Generation Using Imagery Obtained from a Low-Cost Uav |
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 |