CN103323028A - 一种基于物方定位一致性的卫星多光谱影像配准方法 - Google Patents
一种基于物方定位一致性的卫星多光谱影像配准方法 Download PDFInfo
- Publication number
- CN103323028A CN103323028A CN2013102369393A CN201310236939A CN103323028A CN 103323028 A CN103323028 A CN 103323028A CN 2013102369393 A CN2013102369393 A CN 2013102369393A CN 201310236939 A CN201310236939 A CN 201310236939A CN 103323028 A CN103323028 A CN 103323028A
- Authority
- CN
- China
- Prior art keywords
- spectral coverage
- coordinate
- object space
- centerdot
- wgs84
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000003384 imaging method Methods 0.000 claims abstract description 40
- 230000003595 spectral effect Effects 0.000 claims description 170
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000006467 substitution reaction Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 4
- 238000002271 resection Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000001228 spectrum Methods 0.000 abstract 9
- 230000003287 optical effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
一种基于物方定位一致性的卫星多光谱影像配准方法,包括选定多光谱相机中某一谱段作为参考谱段,其余谱段为非参考谱段,针对各非参考谱段与参考谱段之间的相对几何畸变分别进行在轨检校,并将检校结果保存;基于物方定位一致性,利用保存的校验结果建立各谱段的严格几何成像模型,将非参考谱段与参考谱段进行精确配准。本发明提供了一种真正几何意义上的光学卫星多光谱影像高精度自动配准方法,能够在无需影像匹配的情况下对卫星多光谱影像进行高精度自动配准,不仅提高了处理效率,并且配准质量与影像质量、地物类型无关;同时,波段间的几何纠正模型为严格几何成像模型,在理论上具有严密性。实践表明该方法可行、有效,配准质量稳定、精度高。
Description
技术领域
本发明属于遥感影像几何处理领域,涉及一种基于物方定位一致性的卫星多光谱影像配准方法。
背景技术
多光谱相机由于能够获取多个波段的影像,通过后续的配准、融合等处理,能够生成各种专题影像产品,极大地丰富了影像信息,提升了影像数据的应用潜力,已成为当前遥感卫星搭载的重要成像载荷。各波段影像之间的配准是多光谱影像处理的重要环节,配准效率、精度以及可靠性直接影响其后续处理和应用的质量。
目前,基于影像匹配的卫星多光谱影像配准方法已经开展了大量的研究与实践工作,这类方法对匹配质量的依赖使其适用性受到一定限制,对于那些纹理特征不明显、辐射差异较大的区域,此类方法的配准质量难以保证;另外,利用该类方法每次进行影像配准前,都必须进行影像匹配,其配准效率大大降低,严重影响后续的处理与应用。基于此,在卫星影像地面预处理中,有必要研究一种不依赖影像匹配的多光谱影像自动配准方法,避免上述问题的出现。
发明内容
本发明所要解决的问题是提供一种无需影像匹配的卫星多光谱影像自动配准方法。
本发明的技术方案为一种基于物方定位一致性的卫星多光谱影像配准方法,包括以下步骤:
步骤1,选定多光谱相机中某一谱段作为参考谱段,其余谱段为非参考谱段,针对各非参考谱段与参考谱段之间的相对几何畸变分别进行在轨检校,并将检校结果保存;
设参考谱段记为B2,任一非参考谱段记为B1,针对非参考谱段B1与参考谱段B2之间的相对几何畸变进行在轨检校实现方式如下,
步骤1.1,设为物方点P在非参考谱段B1与参考谱段B2影像上有同名像点p1和p2,
构建参考谱段B2的严格几何成像模型如下,
构建非参考谱段B1的严格几何成像模型如下,
上式中,(Xp,Yp,Zp)WGS84为物方点P的WGS84地心直角坐标;(xC1,yC1,-f)和(xC2,yC2,-f)分别为同名像点p1和p2在相机坐标系下的坐标,f代表相机主距;m1和m2为摄影比例尺因子;(XS1,YS1,ZS1)和(XS2,YS2,ZS2)为投影中心S1和S2的WGS84地心直角坐标;RBS为相机在卫星本体坐标系下的安装角矩阵;RBJ1和RT1分别为p1成像时卫星本体坐标系与地心惯性坐标系、地心惯性坐标系与WGS84地心直角坐标系之间的旋转矩阵;RBJ2和RT2则为p2成像时相应的旋转矩阵;对于附加参数项Δx1和Δy1,采用以探元号为自变量的三次多项式如下,
将参考谱段B2相对于非参考谱段B1的几何畸变参数记为,则
采用符号f1表示坐标正投影换算,建立将像点(x,y)正投影至物方获取其对应的物方点的WGS84地心直角坐标(X,Y,Z)WGS84的公式如下,
采用符号f2表示坐标反投影换算,建立将物方点的WGS84地心直角坐标(X,Y,Z)WGS84反投影至像方所得到的像点坐标(x,y)的公式如下,
步骤1.2,在非参考谱段B1与参考谱段B2影像上量测n对同名像点同名像点对应物方点Pi,n为同名像点对的总数,i=1,…,n;对每个同名点对根据参考谱段B2影像上的像点的坐标,利用参考谱段B2的严格几何成像模型及物方高程信息,执行坐标正投影换算f1,将像点投影至物方,获取相应物方点Pi的WGS84地心直角坐标
步骤1.3,利用步骤1.2所得物方点Pi的坐标基于非参考谱段B1的严格几何成像模型,利用空间后方交会的原理解算参考谱段B2相对于非参考谱段B1的几何畸变参数消除非参考谱段B1与参考谱段B2之间的相对几何畸变;
步骤2,基于物方定位一致性,利用各谱段的严格几何成像模型,将非参考谱段与参考谱段进行精确配准;针对非参考谱段B1与参考谱段B2进行精确配准实现方式如下,
步骤2.1,根据参考谱段影像像点坐标和物方高程信息获取物方点坐标,包括对参考谱段B2影像上的每个像元p2(xc2,yc2)执行坐标正投影换算f1,获取其物方点P的WGS84地心直角坐标(XP,YP,ZP)WGS84;
步骤2.2,获取非参考谱段影像像点坐标,包括利用步骤1所得参考谱段B2相对于非参考谱段B1的几何畸变参数构建非参考谱段B1的严格几何成像模型;对步骤2.1所得物方点P的坐标(XP,YP,ZP)WGS84执行坐标反投影计算f2,获得非参考谱段B1影像上对应的像点p1的坐标(xc1,yc1);
步骤2.3,根据步骤2.2所得非参考谱段B1影像上对应的像点p1的坐标(xc1,yc1)进行灰度重采样,完成非参考谱段B1与参考谱段B2的精确配准。
而且,坐标正投影换算和坐标反投影换算的实现方式如下,
设参考谱段B2影像上像点p2的坐标为(xc2,yc2),对应的物方点P的WGS84地心直角坐标为(Xp,Yp,Zp)WGS84,令
a1、a2、a3、b1、b2、b3、c1、c2、c3为矩阵的元素;
根据物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84与其大地坐标(Bp,Lp,Hp)之间的如下关系,
得到像点p2的坐标(xc2,yc2)与大地坐标(Bp,Lp,Hp)关系式如下,
坐标正投影换算时,由像点p2的坐标(xc2,yc2)以及给定的物方高程值Hp,利用上式解算物方点P的大地经纬度(Bp,Lp);再根据物方点P的大地坐标(Bp,Lp,Hp)获取相应WGS84地心直角坐标(Xp,Yp,Zp)WGS84;
坐标反投影换算时,基于参考谱段B2的严格几何成像模型,由物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84,利用下式解算,
得到对应的像点p2的坐标(xc2,yc2)。
令
上式中,矢量 代表从相机投影中心到物方点的矢量在相机坐标系下的坐标;
根据非参考谱段B1的严格几何成像模型,得到下式
令
上式中,vxi和vyi分别代表沿轨和垂轨方向的像方残差;
Vi=AiX-Li,Wi,
其中,
上式中,Vi、Ai、Li分别是利用物方点Pi构建的误差方程,的残差向量、待解参数的系数矩阵以及常向量;X代表参考谱段B2相对于非参考谱段B1的几何畸变参数 Wi是非参考B1谱段影像上的像点的量测精度对应的权;
基于最小二乘平差原理,利用下式计算参考谱段B2相对于非参考谱段B1的几何畸变参数,
记录计算所得几何畸变参数。
本发明选定多光谱相机中某一谱段作为参考谱段,其余谱段为非参考谱段。利用一景质量较优的影像,基于各谱段严格几何成像模型,仅利用谱段间的同名像点,对各非参考谱段其与参考谱段之间的相对几何畸变进行在轨检校;对非参考波段其与参考波段间的相对几何畸变进行在轨检校后,基于同名像元物方定位一致性的约束条件,利用各谱段严格几何成像模型,在无需影像匹配的情况下实现子像素级的多光谱影像自动配准,波段间的几何纠正模型为严格几何成像模型,理论上具有严密性。本发明的优点在于:1.影像配准之前无需进行影像匹配,一方面大大提升了数据处理效率,同时,配准质量与影像质量、地物类型无关,对于水域、沙漠以及山地等纹理特征不丰富、匹配质量难以保证的区域,本发明的配准质量也能得到保证。2.一种真正几何意义上的配准,各波段影像之间的几何纠正模型为严格几何成像模型,在理论上具有严密性,能够进一步保证配准精度。
附图说明
图1为本发明实施例的配准流程示意图;
图2为本发明基于物方定位一致性的卫星多光谱影像配准原理示意图。
具体实施方式
以下结合附图和实施例详细说明本发明具体实施方式。实施例的流程可以分为两个步骤,每个步骤实施的具体方法、公式以及流程如下:
步骤1,选定多光谱相机中某一谱段作为参考谱段,其余谱段为非参考谱段,利用一景质量较优的影像,针对各非参考谱段,对其与参考谱段之间的相对几何畸变进行在轨检校,并将检校结果保存,用于后续的影像配准。设参考谱段记为B2,任一非参考谱段记为B1,针对非参考谱段B1,其与参考谱段B2之间的相对几何畸变在轨检校的具体步骤及公式如下:
步骤1.1,利用姿态、轨道、时间等辅助数据以及相机参数,构建参考谱段B2的严格几何成像模型;如式(1);
同理,构建B1谱段的严格几何成像模型,并在其内定向参数模型中引入附加参数项Δx1和Δy1,用于补偿B1谱段相对于B2谱段的几何畸变,实现构建基于扩展共线条件方程的自检校平差模型,如式(2)。
上式中,(Xp,Yp,Zp)WGS84为物方点P的WGS84地心直角坐标;(xC1,yC1,-f)和(xC2,yC2,-f)分别为物方点P在B1谱段与B2谱段影像上同名像点p1和p2在相机坐标系下的坐标,f代表相机主距;m1和m2为摄影比例尺因子;(XS1,YS1,ZS1)和(XS2,YS2,ZS2)为投影中心S1和S2的WGS84地心直角坐标;RBS为相机在卫星本体坐标系下的安装角矩阵;RBJ1和RT1分别为p1成像时卫星本体坐标系与地心惯性坐标系、地心惯性坐标系与WGS84地心直角坐标系之间的旋转矩阵,地心惯性坐标系为J2000坐标系等;RBJ2和RT2则为p2成像时相应的旋转矩阵。对于B1谱段成像模型中的附加参数项Δx1和Δy1,采用以探元号为自变量的三次多项式,如式(3)所示:
利用卫星获取的姿态、轨道、时间、相机参数以及物方高程信息,基于严格几何成像模型,能够实现像点坐标与物方点坐标之间的换算。以谱段B2为例,对具体的换算过程及公式进行阐述。
假设谱段B2影像上某像点p2的坐标为(xc2,yc2),其对应的物方点P的WGS84地心直角坐标为(Xp,Yp,Zp)WGS84,令公式(1)中:
a1、a2、a3、b1、b2、b3、c1、c2、c3为矩阵的元素。
则利用公式(4)及物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84可解算像点p2的坐标(xc2,yc2):
根据大地测量学的基本原理可知,物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84与其大地坐标(Bp,Lp,Hp)(Bp、Lp分别为物方点P的纬度和经度,Hp为其椭球高)之间有如下关系:
其中,e代表地球椭球扁率,变量a代表地球椭球长半轴。将上式(5)代入公式(4)中有公式(6):
坐标正投影换算时,由像点p2的坐标(xc2,yc2)以及给定的物方高程值Hp,利用公式(6)即可解算物方点P的大地经纬度(Bp,Lp);将物方点P的大地坐标(Bp,Lp,Hp)代入公式(5)即可获取其WGS84地心直角坐标(Xp,Yp,Zp)WGS84。坐标反投影换算时,基于参考谱段B2的严格几何成像模型,由物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84,利用式(4)解算对应的像点p2的坐标(xc2,yc2)
为了便于描述,本文利用公式(7)和(8)简要表达上述的像点坐标与物方点坐标正反换算。其中,公式(7)中的符号f1表示将像点(x,y)正投影至物方获取其对应的物方点的WGS84地心直角坐标(X,Y,Z)WGS84,即坐标正投影换算;公式(8)中的符号f2则表示将物方点的WGS84地心直角坐标(X,Y,Z)WGS84反投影至像方所得到的像点坐标(x,y),即坐标反投影换算。
步骤1.2,在B1和B2两个谱段影像上量测n对同名像点表示B1和B2谱段影像上的一对同名像点),同名像点对应物方点Pi。n为同名像点对的总数,可由本领域技术人员根据具体情况设定。具体量测实现为现有技术,为了提高检校结果的精度以及可靠性,建议在实施时对像点在非参考谱段影像上的分布作如下要求(对应的同名像点在参考谱段影像上的分布形状与在非参考谱段上是一致的):(1)在影像行方向(即推扫方向)上尽量分布在较短的一段区域内;(2)在影像列方向(即沿CCD方向)上均匀覆盖整个线阵CCD。其中,第一点要求是为了降低外方位元素误差对解算结果的影响;第二点则是为了保证解算结果对线阵CCD所有探元均适用。对每个同名点对(i=1,…,n),将B2谱段影像上的像点的坐标代入公式(7)中的(x,y),利用式(1)所示B2谱段的严格几何成像模型及物方高程信息,执行坐标正投影换算f1,将像点投影至物方,获取其物方点Pi的WGS84地心直角坐标即公式(7)中的(X,Y,Z)WGS84;
步骤1.3,利用前述步骤1.2得到的物方点Pi的坐标(i=1,…,n),基于B1谱段的严格几何成像模型(式2),利用空间后方交会的原理解算,消除B1谱段与B2谱段影像之间的相对几何畸变。具体解算公式如下。
1)令
上式中,矢量 称为物方矢量U,代表从相机投影中心到物方点的矢量在相机坐标系下的坐标;
式(2)可转化为式(10):
令
上式中,vxi和vyi分别代表沿轨和垂轨方向的像方残差。
Vi=AiX-Li,Wi(i=1,…,n) (12)
其中,
上式中,Vi、Ai、Li分别是利用物方点Pi构建的误差方程式的残差向量、待解参数的系数矩阵以及常向量;X代表B1谱段相对于B2谱段的几何畸变参数 Wi是B1谱段影像上的像点的量测精度对应的权。
3)基于最小二乘平差原理,利用式(13)计算B1谱段相对于B2谱段的几何畸变参数,实现对B1谱段相对于B2谱段之间的相对几何畸变进行在轨检校,并将几何畸变参数用文件记录下来,用于后续多光谱影像波段配准。
步骤2,通过步骤1预先获取了各非参考谱段其相对于参考谱段的几何畸变参数后,基于物方定位一致性,利用各谱段影像的严格几何成像模型,对任意一景多光谱影像,将其各非参考谱段与参考谱段分别进行精确配准。仍设参考谱段记为B2,任一非参考谱段记为B1,结合图1和图2对本步骤中的具体原理和流程进行阐述如下。其中,图1为流程示意图,图2为原理示意图。图2中,p1和p2分别为B1和B2谱段影像上的一对同名像点,P则代表其对应的物方点;S1和S2分别为像点p1和p2成像时的相机投影中心;f1和f2则代表公式7和8中的像点坐标与物方点坐标之间的正反换算。
步骤2.1,根据参考谱段影像像点坐标和物方高程信息获取物方点坐标:对B2谱段影像上的每个像元p2(xc2,yc2),将其像点坐标(xc2,yc2)代入公式(7)中的(x,y),即利用B2谱段的严格几何成像模型(式1)及物方高程信息,执行坐标正投影换算f1,将像点p2(xc2,yc2)投影至物方,获取其物方点P的WGS84地心直角坐标(XP,YP,ZP)WGS84(即公式(7)中的(X,Y,Z)WGS84);
步骤2.2,获取非参考谱段影像像点坐标:利用步骤1获取的B1谱段其相对于B2谱段的几何畸变参数,构建式(2)所示的严格几何成像模型。将获取的物方点P的坐标(XP,YP,ZP)WGS84代入式(8)中的(X,Y,Z)WGS84,执行坐标反投影计算f2,获得B1谱段影像上对应的像点p1的坐标(xc1,yc1),即公式(8)中的(x,y);由于式(2)已经对B1和B2两谱段的相对几何畸变进行了补偿,基于同名像点物方定位一致性的约束关系,所得到的像点p1(s1,l1)即为B2谱段影像上像点p2(s,l)在B1谱段影像上对应的同名像点,由此可建立B1、B2两波段影像上同名像点之间的映射关系,实现多光谱影像的配准,这就是基于物方定位一致性的多光谱影像配准原理,如图2所示。
步骤2.3,灰度重采样:在步骤2.1获取B2谱段影像的每个像元p2(s,l)在B1谱段影像上对应的同名像点坐标p1(s1,l1)后,通过灰度重采样,即可实现B1与B2两谱段影像的配准。其中,灰度重采样可采用现有技术中的双线性内插算法实现,本发明不予赘述。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (3)
1.一种基于物方定位一致性的卫星多光谱影像配准方法,包括以下步骤:
步骤1,选定多光谱相机中某一谱段作为参考谱段,其余谱段为非参考谱段,针对各非参考谱段与参考谱段之间的相对几何畸变分别进行在轨检校,并将检校结果保存;
设参考谱段记为B2,任一非参考谱段记为B1,针对非参考谱段B1与参考谱段B2之间的相对几何畸变进行在轨检校实现方式如下,
步骤1.1,设为物方点P在非参考谱段B1与参考谱段B2影像上有同名像点p1和p2,
构建参考谱段B2的严格几何成像模型如下,
构建非参考谱段B1的严格几何成像模型如下,
上式中,(Xp,Yp,Zp)WGS84为物方点P的WGS84地心直角坐标;(xC1,yC1,-f)和(xC2,yC2,-f)分别为同名像点p1和p2在相机坐标系下的坐标,f代表相机主距;m1和m2为摄影比例尺因子;(XS1,YS1,ZS1)和(XS2,YS2,ZS2)为投影中心S1和S2的WGS84地心直角坐标;RBS为相机在卫星本体坐标系下的安装角矩阵;RBJ1和RT1分别为p1成像时卫星本体坐标系与地心惯性坐标系、地心惯性坐标系与WGS84地心直角坐标系之间的旋转矩阵;RBJ2和RT2则为p2成像时相应的旋转矩阵;对于附加参数项Δx1和Δy1,采用以探元号为自变量的三次多项式如下,
采用符号f1表示坐标正投影换算,建立将像点(x,y)正投影至物方获取其对应的物方点的WGS84地心直角坐标(X,Y,Z)WGS84的公式如下,
采用符号f2表示坐标反投影换算,建立将物方点的WGS84地心直角坐标(X,Y,Z)WGS84反投影至像方所得到的像点坐标(x,y)的公式如下,
步骤1.2,在非参考谱段B1与参考谱段B2影像上量测n对同名像点同名像点对应物方点Pi,n为同名像点对的总数,i=1,…,n;对每个同名点对根据参考谱段B2影像上的像点的坐标,利用参考谱段B2的严格几何成像模型及物方高程信息,执行坐标正投影换算f1,将像点投影至物方,获取相应物方点Pi的WGS84地心直角坐标
步骤1.3,利用步骤1.2所得物方点Pi的坐标基于非参考谱段B1的严格几何成像模型,利用空间后方交会的原理解算参考谱段B2相对于非参考谱段B1的几何畸变参数消除非参考谱段B1与参考谱段B2之间的相对几何畸变;
步骤2,基于物方定位一致性,利用各谱段的严格几何成像模型,将非参考谱段与参考谱段进行精确配准;针对非参考谱段B1与参考谱段B2进行精确配准实现方式如下,
步骤2.1,根据参考谱段影像像点坐标和物方高程信息获取物方点坐标,包括对参考谱段B2影像上的每个像元p2(xc2,yc2)执行坐标正投影换算f1,获取其物方点P的WGS84地心直角坐标(XP,YP,ZP)WGS84;
步骤2.2,获取非参考谱段影像像点坐标,包括利用步骤1所得参考谱段B2相对于非参考谱段B1的几何畸变参数构建非参考谱段B1的严格几何成像模型;对步骤2.1所得物方点P的坐标(XP,YP,ZP)WGS84执行坐标反投影计算f2,获得非参考谱段B1影像上对应的像点p1的坐标(xc1,yc1);
步骤2.3,根据步骤2.2所得非参考谱段B1影像上对应的像点p1的坐标(xc1,yc1)进行灰度重采样,完成非参考谱段B1与参考谱段B2的精确配准。
2.如权利要求1所述基于物方定位一致性的卫星多光谱影像配准方法,其特征在于:坐标正投影换算和坐标反投影换算的实现方式如下,
设参考谱段B2影像上像点p2的坐标为(xc2,yc2),对应的物方点P的WGS84地心直角坐标为(Xp,Yp,Zp)WGS84,令
a1、a2、a3、b1、b2、b3、c1、c2、c3为矩阵的元素;
根据物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84与其大地坐标(Bp,Lp,Hp)之间的如下关系,
得到像点p2的坐标(xc2,yc2)与大地坐标(Bp,Lp,Hp)关系式如下,
坐标正投影换算时,由像点p2的坐标(xc2,yc2)以及给定的物方高程值Hp,利用上式解算物方点P的大地经纬度(Bp,Lp);再根据物方点P的大地坐标(Bp,Lp,Hp)获取相应WGS84地心直角坐标(Xp,Yp,Zp)WGS84;
坐标反投影换算时,基于参考谱段B2的严格几何成像模型,由物方点P的WGS84地心直角坐标(Xp,Yp,Zp)WGS84,利用下式解算,
得到对应的像点p2的坐标(xc2,yc2)。
令
上式中,矢量 代表从相机投影中心到物方点的矢量在相机坐标系下的坐标;
根据非参考谱段B1的严格几何成像模型,得到下式
令
上式中,vxi和vyi分别代表沿轨和垂轨方向的像方残差;
Vi=AiX-Li,Wi,
其中,
上式中,Vi、Ai、Li分别是利用物方点Pi构建的误差方程,的残差向量、待解参数的系数矩阵以及常向量;X代表参考谱段B2相对于非参考谱段B1的几何畸变参数 Wi是非参考B1谱段影像上的像点的量测精度对应的权;
基于最小二乘平差原理,利用下式计算参考谱段B2相对于非参考谱段B1的几何畸变参数,
记录计算所得几何畸变参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310236939.3A CN103323028B (zh) | 2013-06-14 | 2013-06-14 | 一种基于物方定位一致性的卫星多光谱影像配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310236939.3A CN103323028B (zh) | 2013-06-14 | 2013-06-14 | 一种基于物方定位一致性的卫星多光谱影像配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103323028A true CN103323028A (zh) | 2013-09-25 |
CN103323028B CN103323028B (zh) | 2015-10-21 |
Family
ID=49191925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310236939.3A Active CN103323028B (zh) | 2013-06-14 | 2013-06-14 | 一种基于物方定位一致性的卫星多光谱影像配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103323028B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635241A (zh) * | 2015-03-05 | 2015-05-20 | 北京航空航天大学 | 一种艇载多光谱高光谱对地观测装置 |
CN105046679A (zh) * | 2014-11-19 | 2015-11-11 | 航天东方红卫星有限公司 | 遥感卫星图像多波段配准的方法及装置 |
CN105444778A (zh) * | 2015-11-10 | 2016-03-30 | 北京空间飞行器总体设计部 | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 |
CN105701830A (zh) * | 2016-01-18 | 2016-06-22 | 武汉大学 | 基于几何模型的lasis波段影像配准方法及*** |
CN107036629A (zh) * | 2017-04-20 | 2017-08-11 | 武汉大学 | 视频卫星在轨相对辐射定标方法及*** |
US10341565B2 (en) | 2016-05-10 | 2019-07-02 | Raytheon Company | Self correcting adaptive low light optical payload |
CN110006452A (zh) * | 2019-04-17 | 2019-07-12 | 武汉大学 | 高分六号宽视场相机相对几何定标方法及*** |
CN112066950A (zh) * | 2020-07-24 | 2020-12-11 | 北京空间机电研究所 | 一种多光轴平行的测绘相机单中心投影转换方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3422471B2 (ja) * | 1999-04-05 | 2003-06-30 | 日本電気株式会社 | 被写体監視方法、被写体追跡方法及び装置 |
CN101922930A (zh) * | 2010-07-08 | 2010-12-22 | 西北工业大学 | 一种航空偏振多光谱图像配准方法 |
CN102901519A (zh) * | 2012-11-02 | 2013-01-30 | 武汉大学 | 一种基于探元指向角光学推扫卫星在轨分步几何定标方法 |
-
2013
- 2013-06-14 CN CN201310236939.3A patent/CN103323028B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3422471B2 (ja) * | 1999-04-05 | 2003-06-30 | 日本電気株式会社 | 被写体監視方法、被写体追跡方法及び装置 |
CN101922930A (zh) * | 2010-07-08 | 2010-12-22 | 西北工业大学 | 一种航空偏振多光谱图像配准方法 |
CN102901519A (zh) * | 2012-11-02 | 2013-01-30 | 武汉大学 | 一种基于探元指向角光学推扫卫星在轨分步几何定标方法 |
Non-Patent Citations (2)
Title |
---|
HUISHENG ZHANG: ""Accurate Multi-spectral Image Registration Based on Scale Invariant Feature"", 《2012 2ND INTERNATIONAL CONFERENCE ON COMPUTER SCIENCE AND NETWORK TECHNOLOGY》, 31 December 2012 (2012-12-31), pages 847 - 852, XP032419976, DOI: doi:10.1109/ICCSNT.2012.6526062 * |
胡芬,等: ""基于投影基准面的线阵推扫式卫星立体像对近似核线影像生成方法"", 《测绘学报》, vol. 38, no. 5, 31 October 2009 (2009-10-31) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046679A (zh) * | 2014-11-19 | 2015-11-11 | 航天东方红卫星有限公司 | 遥感卫星图像多波段配准的方法及装置 |
CN104635241A (zh) * | 2015-03-05 | 2015-05-20 | 北京航空航天大学 | 一种艇载多光谱高光谱对地观测装置 |
CN105444778A (zh) * | 2015-11-10 | 2016-03-30 | 北京空间飞行器总体设计部 | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 |
CN105444778B (zh) * | 2015-11-10 | 2018-10-09 | 北京空间飞行器总体设计部 | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 |
CN105701830A (zh) * | 2016-01-18 | 2016-06-22 | 武汉大学 | 基于几何模型的lasis波段影像配准方法及*** |
CN105701830B (zh) * | 2016-01-18 | 2018-09-21 | 武汉大学 | 基于几何模型的lasis波段影像配准方法及*** |
US10341565B2 (en) | 2016-05-10 | 2019-07-02 | Raytheon Company | Self correcting adaptive low light optical payload |
CN107036629A (zh) * | 2017-04-20 | 2017-08-11 | 武汉大学 | 视频卫星在轨相对辐射定标方法及*** |
CN110006452A (zh) * | 2019-04-17 | 2019-07-12 | 武汉大学 | 高分六号宽视场相机相对几何定标方法及*** |
CN110006452B (zh) * | 2019-04-17 | 2023-06-23 | 武汉大学 | 高分六号宽视场相机相对几何定标方法及*** |
CN112066950A (zh) * | 2020-07-24 | 2020-12-11 | 北京空间机电研究所 | 一种多光轴平行的测绘相机单中心投影转换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103323028B (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103323028A (zh) | 一种基于物方定位一致性的卫星多光谱影像配准方法 | |
CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
CN110675450B (zh) | 基于slam技术的正射影像实时生成方法及*** | |
Teo et al. | DEM-aided block adjustment for satellite images with weak convergence geometry | |
CN103759714A (zh) | 一种三线阵卫星影像区域网平差方法 | |
CN102968631A (zh) | 山区多光谱遥感卫星影像的自动几何纠正与正射校正方法 | |
CN103129752A (zh) | 一种基于地面导航的光学遥感卫星姿态角误差动态补偿方法 | |
CN105551053A (zh) | 一种小面阵星载tdi ccd相机的快速几何精校正方法 | |
CN103093459A (zh) | 利用机载LiDAR点云数据辅助影像匹配的方法 | |
Cao et al. | Bundle adjustment of satellite images based on an equivalent geometric sensor model with digital elevation model | |
CN113514829B (zh) | 面向InSAR的初始DSM的区域网平差方法 | |
CN102901519A (zh) | 一种基于探元指向角光学推扫卫星在轨分步几何定标方法 | |
Jiang et al. | CCD distortion calibration without accurate ground control data for pushbroom satellites | |
Liu et al. | A new approach to fast mosaic UAV images | |
Fraser et al. | Precise georefrencing of long strips of ALOS imagery | |
CN109029379B (zh) | 一种高精度小基高比立体测绘方法 | |
Wang et al. | Geometric calibration for the aerial line scanning camera GFXJ | |
Tao et al. | On-orbit geometric calibration of the panchromatic/multispectral camera of the ZY-1 02C satellite based on public geographic data | |
CN116753916B (zh) | 多视角卫星影像区域网平差方法及*** | |
Zhao et al. | Digital Elevation Model‐Assisted Aerial Triangulation Method On An Unmanned Aerial Vehicle Sweeping Camera System | |
Liu et al. | High precision DTM and DOM generating using multi-source orbital data on Chang’e-4 landing site | |
Liu et al. | Development of an Attitude Transformation Method From the Navigation Coordinate System to the Projection Coordinate System | |
CN111044076B (zh) | 基于参考底图的高分一号b卫星几何检校方法 | |
Shen et al. | An improved method for transforming GPS/INS attitude to national map projection frame | |
Zhang et al. | Calibrating an airborne linear-array multi-camera system on the master focal plane with existing bundled images |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240104 Address after: No. 120, Haixiang Middle Road, Fengcheng Street, Haiyang City, Yantai City, Shandong Province, 265100 Patentee after: Land sea space (Yantai) Information Technology Co.,Ltd. Address before: 430072 Hubei Province, Wuhan city Wuchang District of Wuhan University Luojiashan Patentee before: WUHAN University |
|
TR01 | Transfer of patent right |