CN103778612B - 一种基于全色影像的卫星颤振探测与补偿方法 - Google Patents
一种基于全色影像的卫星颤振探测与补偿方法 Download PDFInfo
- Publication number
- CN103778612B CN103778612B CN201410050336.9A CN201410050336A CN103778612B CN 103778612 B CN103778612 B CN 103778612B CN 201410050336 A CN201410050336 A CN 201410050336A CN 103778612 B CN103778612 B CN 103778612B
- Authority
- CN
- China
- Prior art keywords
- image
- flutter
- points
- point
- satellite
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000001514 detection method Methods 0.000 title claims abstract description 11
- 206010044565 Tremor Diseases 0.000 title abstract 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 abstract description 9
- 230000008878 coupling Effects 0.000 abstract 2
- 238000010168 coupling process Methods 0.000 abstract 2
- 238000005859 coupling reaction Methods 0.000 abstract 2
- 238000012892 rational function Methods 0.000 description 3
- 241000132092 Aster Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000011439 discrete element method Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 235000018936 Vitellaria paradoxa Nutrition 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004366 reverse phase liquid chromatography Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于全色影像的卫星颤振探测与补偿方法,包括以下步骤:1)SIFT特征检测器匹配少量均匀分布的同名点,计算像方仿射模型补偿参数改正数X1;2)利用基于几何约束的互相关方法在三视影像上进行同名点的匹配;3)使用最小二乘匹配方法对所有点进行精细匹配;4)剔除粗差点;5)将三线阵立体影像中的前视影像与下视影像、后视影像与下视影像,分别构成两组立体像对FN和BN,并分别建立误差方程,利用最小二乘方法解算正弦函数模型补偿参数改正数X2;6)判断X2是否大于阈值;7)输出仿射模型参数和正弦函数模型参数。与现有技术相比,本发明具有能够有效的消除非线性的、周期性的平台颤振对影像定位造成的影响。
Description
技术领域
本发明涉及一种卫星颤振探测与补偿方法,尤其是涉及一种基于全色影像的卫星颤振探测与补偿方法。
背景技术
高分辨率卫星在轨运行期间,由于星上转动部件的周期性运动或变轨、冷热交变等引发扰动,使星体产生一种幅值较小、具有一定周期性的颤振响应。NASA于上世纪八十年代深入研究了HST的反作用轮扰动影像,1990年Eyerman和Shea对卫星的颤振动力学环境进行了总结分析,认为最大的颤振源是反作用轮和热震颤。Iwasaki等利用ASTER卫星的多光谱影像波段间的配准,分析了ASTER卫星颤振响应。Lehner等结合地面DEM验证MOMS-2P影像中含有颤振分量。Schwind等利用ALOS-PRISM影像分析RPCs参数与相机严格模型的定位结果对比,发现pitch、yaw和roll三个方向上具有明显的颤振现象。这些研究从多重角度对卫星颤振进行了分析,但是这些分析中,利用了外部的DEM或者相机的严格模型,这些数据普通用户很难获得。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于全色影像的卫星颤振探测与补偿方法,在无控制点参与的条件下,精确探测高分辨率卫星的颤振特性,并对其进行有效的补偿,消除颤振影响。
本发明的目的可以通过以下技术方案来实现:
一种基于全色影像的卫星颤振探测与补偿方法,其特征在于,包括以下步骤:
1)SIFT特征检测器匹配少量均匀分布的同名点,计算像方仿射模型补偿参数改正数X1;
2)利用基于几何约束的互相关方法在三视影像上进行同名点的匹配,获取更多的同名点;
3)将步骤2)得到的同名点结果作为初始值,使用最小二乘匹配方法对所有点进行精细匹配;
4)匹配完成后,根据点的相关系数阈值以及同名像点前方交会的残差剔除粗差点;
5)将三线阵立体影像中的前视影像与下视影像、后视影像与下视影像,分别构成两组立体像对FN和BN,并分别建立误差方程,联立两组误差方程,其中下视影像的改正数在两组误差方程中相同,利用最小二乘方法解算正弦函数模型补偿参数改正数X2;
6)判断X2是否大于阈值,如果是,返回步骤5),如果否,则执行步骤7);
7)输出仿射模型参数和正弦函数模型参数。
所述的像方仿射模型补偿参数改正数X1计算如下:
101)建立像方平面误差的基本模型,并通过以下公式表示:
式中,(x,y)是像方同名点的列号和行号的观测值,(x′,y′)是由地面点坐标反投影到像方得到的像方坐标的列号和行号,ai、bi、ci为仿射模型补偿的系数,其中i=1、2,A1,w1,为正弦函数模型的系数;
102)将公式(1)中的仿射模型补偿部分整理为:
v1=B1X1-L1 (2)
其中,v1=[vx vy]T
X1=[a1 b1 c1 a2 b2 c2]T
通过三视影像的前方交会,计算得到对应的地面点坐标,将地面点坐标反投影至各影像上,量测的像方坐标与反投影计算的像方坐标之间存在相对偏差L1,利用公式(2)计算仿射模型系数X1。
所述的点的相关系数阈值为0.7.
所述的正弦函数模型补偿参数改正数X2具体计算如下:
采用正弦函数模型描述立体像对列方向上的颤振误差:
为量测的像方坐标经仿射模型改正后的值,为无颤振影响的同名像点列坐标的理论值,Ai、wi、分别为正弦函数的振幅、角频率和相位,ti为像点所在影像行的曝光时刻,其中i=1、2;
投影点p′1与量测点p2之间的偏差d表示为:
式中,A′1为投影点p′1的理论值和颤振影响的振幅,且与相等;
式(5)按泰勒公式展开至一次项,得误差方程:
其中,
X2为参数的改正数,d0为d的近似值;
式(6)可以整理为:
v2=B2X2-L2 (7)
与现有技术相比,本发明具有能够有效的消除非线性的、周期性的平台颤振对影像定位造成的影响。实验结果表明,在垂轨方向上存在约0.61~0.62HZ的颤振运动,振幅约为0.5-1个像素,通过像方补偿后,颤振引起的几何偏差从1个都像素减小至半个像素以内,并且周期性的颤振现象完全消除。
附图说明
图1为本发明的流程图;
图2为三线阵立体影像图;
图3为改正前的下视与前视影像的相对像方偏差在下视影像上的投影图;
图4为改正后的下视与前视影像的相对像方偏差在下视影像上的投影图;
图5为改正前的下视与后视影像的相对像方偏差在下视影像上的投影图;
图6为改正后的下视与后视影像的相对像方偏差在下视影像上的投影图;
图7为改正前的前视与后视影像的相对像方偏差在前视影像上的投影图;
图8为改正后的前视与后视影像的相对像方偏差在前视影像上的投影图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1所示,一种基于全色影像的卫星颤振探测与补偿方法,其特征在于,包括以下步骤:
1)SIFT特征检测器匹配少量均匀分布的同名点,计算像方仿射模型补偿参数改正数X1;
2)利用基于几何约束的互相关方法在三视影像上进行同名点的匹配,获取更多的同名点;
3)将步骤2)得到的同名点结果作为初始值,使用最小二乘匹配方法对所有点进行精细匹配;
4)匹配完成后,根据点的相关系数阈值以及同名像点前方交会的残差剔除粗差点;
5)将三线阵立体影像中的前视影像与下视影像、后视影像与下视影像,分别构成两组立体像对FN和BN,并分别建立误差方程,联立两组误差方程,其中下视影像的改正数在两组误差方程中相同,利用最小二乘方法解算正弦函数模型补偿参数改正数X2;
6)判断X2是否大于阈值,如果是,返回步骤5),如果否,则执行步骤7);
7)输出仿射模型参数和正弦函数模型参数。
1、有理函数模型(Rational function model,RMF)
有理函数模型(Rational function model,RMF)是一种独立于传感器平台物理意义、形式简单的通用传感器模型,广泛应用于目前商业卫星影像中。RMF是用地面点大地坐标作为自变量,用多项式的比值描述其对应的像点坐标。为了增强参数解的稳定性,对大地坐标和像点坐标标准化到-1到+1之间(式1)。
其中,l、s为标准化的影像行列坐标,B、P、H为标准化的大地坐标。
2、颤振影响建模
受卫星测轨精度以及卫星平台颤振运动的影响,卫星影像成像过程中存在线性的***性偏差和非线性的正弦颤振偏差,传统的相对定向过程仅考虑了***性线性偏差,而由于卫星平台颤振引起的周期性偏差仍然存在。针对该问题,本文创新性的提出一种无控制点参与的分步偏差补偿方案。在像方空间上,1)采用基于像方仿射模型的相对定向,补偿影像间相对的线性偏差;2)采用正弦函数模型来描述立体像对成像过程中由卫星平台颤振造成的周期性颤振偏差,并对其进行补偿。对于同轨前后视立体影像,行方向视差体现地形影响,因此主要在列方向上对颤振进行探测和建模。像方平面误差的基本模型可以表示为:
式中,(x,y)是像方同名点的列号和行号的观测值,(x′,y′)是由地面点坐标反投影到像方得到的像方坐标的列号和行号,ai,bi,ci(i=1,2)为仿射模型补偿的系数,A1,w1,为正弦函数模型的系数。
3、仿射模型补偿
式(2)中的仿射模型补偿部分可以整理为:
v1=B1X1-L1 (3)
其中,v1=[vx vy]T
X1=[a1 b1 c1 a2 b2 c2]T
通过三视影像的前方交会,计算得到对应的地面点坐标,将地面点坐标反投影至各影像上,量测的像方坐标与反投影计算的像方坐标之间存在相对偏差L1,利用式(3)计算仿射模型系数X1。
3、正弦函数模型补偿
在理想的立体像对模型中,由同名像点(p1,p2)定位得到的地面点坐标G,投影到影像上所得到的像点坐标(投影点p′1)与对应的量测像点坐标(量测点p2)相同。但是由于卫星平台的颤振影响,投影点与量测点之间存在一定的偏差d。为了解算影像各自的颤振值,在利用最小二乘估计地面点坐标G过程中,左影像给定较大的权重值10000,右影像给定较小权重值1,因此地面点坐标G的误差中主要包含了左影像的颤振分量,而不包含右影像的颤振分量。将地面点坐标G投影到右影像上得到的投影点p′1中也仅含有左影像的颤振分量,因此投影点p′1与量测点p2之间的偏差d,即为立体像对的两张影像在t1,t2时刻受到颤振影响的差值。
经过仿射模型补偿后,影像之间相对的线性偏差已经消除,可以看做理想的立体像对模型,剩余的偏差主要为卫星平台的非线性颤振引起。本文中采用正弦函数模型描述立体像对列方向上的颤振误差。
式中,为量测的像方坐标经仿射模型改正后的值,为无颤振影响的同名像点列坐标的理论值,Ai,wi,分别为正弦函数的振幅、角频率和相位,ti为像点所在影像行的曝光时刻。
投影点p′1与量测点p2之间的偏差d表示为:
式中,A′1为投影点p′1的理论值和颤振影响的振幅,且与相等。
式(5-2)按泰勒公式展开至一次项,得误差方程:
其中,
X2为参数的改正数,d0为d的近似值。
式(6-1)可以整理为:
v2=B2X2-L2 (6-2)
实验数据与结果
实验数据
如图2所示,实验数据采用国产某测绘卫星在山西地区的同轨三线阵立体影像,影像的前后视分辨率为3.5米,下视分辨率为2.1米,幅宽为52千米,前后视相机与下视相机的夹角为±22度。
实验结果
表1中为三视影像的颤振分量的正弦函数参数值,在垂轨方向上存在约0.61~0.62HZ的颤振运动,三视影像的颤振频率基本相同,颤振幅值从前视到后视逐渐减小。表2为颤振补偿后的两组立体模型,分别进行前方交会得到的地面点坐标,投影在下视影像上的列方向偏差,可以看出经过仿射模型和正弦函数模型补偿后的残差在半个像素以内。图3至图8为像方补偿前后在像方空间的偏差,图中显示在像方补偿前,三组立体模型在像方偏差均具有明显的非线性颤振现象,且呈现一定的周期性。经过仿射模型和正弦模型改正后,非线性的颤振现象基本消失,偏差呈线性分布,且集中在0值附近。
表1颤振探测结果
表2颤振补偿结果(单位:像素)
由于卫星平台颤振,卫星影像上会存在颤振引起的非线性的几何偏移,这种非线性的偏差很难通过常规的定向方法进行探测和补偿。而本文提出的不依赖控制点的,基于全色影像的卫星平台非线性颤振探测与补偿算法能够有效的消除非线性的、周期性的平台颤振对影像定位造成的影响。实验结果表明,在垂轨方向上存在约0.61~0.62HZ的颤振运动,振幅约为0.5-1个像素,通过像方补偿后,颤振引起的几何偏差从1个都像素减小至半个像素以内,并且周期性的颤振现象完全消除。
Claims (3)
1.一种基于全色影像的卫星颤振探测与补偿方法,其特征在于,包括以下步骤:
1)SIFT特征检测器匹配少量均匀分布的同名点,计算像方仿射模型补偿参数改正数X1;
2)利用基于几何约束的互相关方法在三视影像上进行同名点的匹配,获取更多的同名点;
3)将步骤2)得到的同名点结果作为初始值,使用最小二乘匹配方法对所有点进行精细匹配;
4)匹配完成后,根据点的相关系数阈值以及同名像点前方交会的残差剔除粗差点;
5)将三线阵立体影像中的前视影像与下视影像、后视影像与下视影像,分别构成两组立体像对FN和BN,并分别建立误差方程,联立两组误差方程,其中下视影像的改正数在两组误差方程中相同,利用最小二乘方法解算正弦函数模型补偿参数改正数X2;
6)判断X2是否大于阈值,如果是,返回步骤5),如果否,则执行步骤7);
7)输出仿射模型参数和正弦函数模型参数。
2.根据权利要求1所述的一种基于全色影像的卫星颤振探测与补偿方法,其特征在于,所述的点的相关系数阈值为0.7。
3.根据权利要求1所述的一种基于全色影像的卫星颤振探测与补偿方法,其特征在于,所述的正弦函数模型补偿参数改正数X2具体计算如下:
采用正弦函数模型描述立体像对列方向上的颤振误差:
为量测的像方坐标经仿射模型改正后的值,为无颤振影响的同名像点列坐标的理论值,Ai、wi、分别为正弦函数的振幅、角频率和相位,ti为像点所在影像行的曝光时刻,其中i=1、2;
投影点p′1与量测点p2之间的偏差d表示为:
式中,A′1为投影点p′1的理论值和颤振影响的振幅,且与相等;式(3)按泰勒公式展开至一次项,得误差方程:
其中,
X2为参数的改正数,d0为d的近似值;
式(4)可以整理为:
v2=B2X2-L2 (5)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410050336.9A CN103778612B (zh) | 2014-02-13 | 2014-02-13 | 一种基于全色影像的卫星颤振探测与补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410050336.9A CN103778612B (zh) | 2014-02-13 | 2014-02-13 | 一种基于全色影像的卫星颤振探测与补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103778612A CN103778612A (zh) | 2014-05-07 |
CN103778612B true CN103778612B (zh) | 2016-08-24 |
Family
ID=50570811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410050336.9A Active CN103778612B (zh) | 2014-02-13 | 2014-02-13 | 一种基于全色影像的卫星颤振探测与补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103778612B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983343B (zh) * | 2014-05-29 | 2016-05-11 | 武汉大学 | 一种基于多光谱影像的卫星平台震颤检测方法及*** |
CN104392457B (zh) * | 2014-12-11 | 2017-07-11 | 中国测绘科学研究院 | 倾斜影像的连接点自动匹配方法及装置 |
CN104864853B (zh) * | 2015-05-08 | 2017-06-13 | 同济大学 | 一种高分辨率三线阵卫星沿轨方向姿态颤振的探测方法 |
CN105469439B (zh) * | 2015-11-29 | 2018-01-19 | 北京工业大学 | 采用残差反馈的slm显微视觉数据重构方法 |
CN110796641B (zh) * | 2019-10-08 | 2022-02-01 | 武汉大学 | 基于连续快照模型的高分辨率卫星影像震颤检测方法 |
CN115773740B (zh) * | 2023-02-14 | 2023-04-21 | 中国船舶集团有限公司第七〇七研究所 | 一种谐振陀螺控制方法及控制*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4023544C2 (de) * | 1990-07-25 | 2002-09-19 | Eads Deutschland Gmbh | Anordnung zur Kompensation von Jitter eines digitalen Signals |
CN101521740A (zh) * | 2009-04-01 | 2009-09-02 | 北京航空航天大学 | 基于多尺度不变特征的实时运动估计方法 |
CN101750200B (zh) * | 2009-12-30 | 2011-06-15 | 航天东方红卫星有限公司 | 一种高分辨率小卫星颤振响应的确定方法 |
CN102289684A (zh) * | 2011-07-29 | 2011-12-21 | 北京航空航天大学 | 一种机载光学景像匹配实时图像自适应处理方法 |
CN102780846A (zh) * | 2012-07-11 | 2012-11-14 | 清华大学 | 一种基于惯导信息的电子稳像方法 |
CN103402045A (zh) * | 2013-08-20 | 2013-11-20 | 长沙超创电子科技有限公司 | 一种基于分区匹配和仿射模型相结合的图像去旋稳像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010109598A (ja) * | 2008-10-29 | 2010-05-13 | Canon Inc | 画像処理装置および画像処理方法 |
-
2014
- 2014-02-13 CN CN201410050336.9A patent/CN103778612B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4023544C2 (de) * | 1990-07-25 | 2002-09-19 | Eads Deutschland Gmbh | Anordnung zur Kompensation von Jitter eines digitalen Signals |
CN101521740A (zh) * | 2009-04-01 | 2009-09-02 | 北京航空航天大学 | 基于多尺度不变特征的实时运动估计方法 |
CN101750200B (zh) * | 2009-12-30 | 2011-06-15 | 航天东方红卫星有限公司 | 一种高分辨率小卫星颤振响应的确定方法 |
CN102289684A (zh) * | 2011-07-29 | 2011-12-21 | 北京航空航天大学 | 一种机载光学景像匹配实时图像自适应处理方法 |
CN102780846A (zh) * | 2012-07-11 | 2012-11-14 | 清华大学 | 一种基于惯导信息的电子稳像方法 |
CN103402045A (zh) * | 2013-08-20 | 2013-11-20 | 长沙超创电子科技有限公司 | 一种基于分区匹配和仿射模型相结合的图像去旋稳像方法 |
Non-Patent Citations (4)
Title |
---|
Geometric Processing of QuickBird Stereo Imageries for Urban Land Use Mapping: A Case Study in Shanghai, China;Xiaohua Tong 等;《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》;20090909;第2卷(第2期);第61-66页 * |
Influence Analysis of jitter on Image quality of remote sensing camera;Fan Chao 等;《Environmental Science and Information Application Technology》;20090705;第461-464页 * |
基于分时积分亚像元融合的地球静止轨道平台消颤振技术;陶小平 等;《光子学报》;20121130;第41卷(第11期);第1359-1364页 * |
复杂卫星抖动下的星敏感器姿态测量数据处理技术;王炯琦 等;《电子与信息学报》;20100831;第32卷(第8期);第1886-1892页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103778612A (zh) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103778612B (zh) | 一种基于全色影像的卫星颤振探测与补偿方法 | |
JP6483075B2 (ja) | 場面の3dパノラマモザイキングの方法 | |
Reigber et al. | Refined estimation of time-varying baseline errors in airborne SAR interferometry | |
CN107014399B (zh) | 一种星载光学相机-激光测距仪组合***联合检校方法 | |
JP4448187B2 (ja) | 映像の幾何補正方法およびその装置 | |
EP3041228A1 (en) | Camera calibration device, camera calibration system, and camera calibration method | |
CN102168972B (zh) | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 | |
CN104864853B (zh) | 一种高分辨率三线阵卫星沿轨方向姿态颤振的探测方法 | |
CN112597428B (zh) | 基于rfm模型光束法平差和影像重采样的颤振探测修正方法 | |
JP6413595B2 (ja) | 画像処理装置、システム、画像処理方法およびプログラム | |
KR20150112362A (ko) | 깊이 센서의 깊이 보정을 위한 이미지 처리방법 및 그 장치 | |
CN107967700A (zh) | 大视场宽工作距双目相机的在轨几何校正及精度验证方法 | |
de Franchis et al. | Automatic sensor orientation refinement of Pléiades stereo images | |
Zhao et al. | Direct georeferencing of oblique and vertical imagery in different coordinate systems | |
CN111156969A (zh) | 一种宽幅遥感影像立体测绘方法及*** | |
CN109029379B (zh) | 一种高精度小基高比立体测绘方法 | |
CN114998399A (zh) | 一种异源光学遥感卫星影像立体像对预处理方法 | |
JP5991821B2 (ja) | 写真測量装置 | |
Re et al. | Performance evaluation of 3DPD, the photogrammetric pipeline for the CaSSIS stereo images | |
KR100870894B1 (ko) | 선형 스캐닝 영상의 자동 영상기하 보정 방법 | |
JP5619866B2 (ja) | 対称な露光写真を利用した地球観測システムに係るアライメント誤差のキャリブレーション方法 | |
KR101694294B1 (ko) | 영상에서 평행사변형과 다중 뷰 제한을 이용한 카메라의 보정 장치 및 방법 | |
CN107146281B (zh) | 一种月球表面高分辨率dem提取方法 | |
JP5901379B2 (ja) | 撮像装置校正方法および画像合成装置 | |
KR101690950B1 (ko) | 최소 기준점을 이용한 푸시브롬 스캔 방식 스트립 위성영상의 기하보정 방법 및 이의 시스템 |
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 |