法律状态公告日
法律状态信息
法律状态
2015-10-21
授权
授权
2013-10-30
实质审查的生效 IPC(主分类):G01C25/00 申请日:20130614
实质审查的生效
2013-09-25
公开
公开
技术领域
本发明属于遥感影像几何处理领域,涉及一种基于物方定位一致性的卫星多光谱影像配 准方法。
背景技术
多光谱相机由于能够获取多个波段的影像,通过后续的配准、融合等处理,能够生成各 种专题影像产品,极大地丰富了影像信息,提升了影像数据的应用潜力,已成为当前遥感卫 星搭载的重要成像载荷。各波段影像之间的配准是多光谱影像处理的重要环节,配准效率、 精度以及可靠性直接影响其后续处理和应用的质量。
目前,基于影像匹配的卫星多光谱影像配准方法已经开展了大量的研究与实践工作,这 类方法对匹配质量的依赖使其适用性受到一定限制,对于那些纹理特征不明显、辐射差异较 大的区域,此类方法的配准质量难以保证;另外,利用该类方法每次进行影像配准前,都必 须进行影像匹配,其配准效率大大降低,严重影响后续的处理与应用。基于此,在卫星影像 地面预处理中,有必要研究一种不依赖影像匹配的多光谱影像自动配准方法,避免上述问题 的出现。
发明内容
本发明所要解决的问题是提供一种无需影像匹配的卫星多光谱影像自动配准方法。
本发明的技术方案为一种基于物方定位一致性的卫星多光谱影像配准方法,包括以下步 骤:
步骤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,采用以探元号为自变量的三次多项式如下,
其中,为三次多项式的系数,s代表探元号;
将参考谱段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)之间的 如下关系,
其中,e代表地球椭球扁率,变量a代表地球椭球长半轴;
得到像点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)。
而且,步骤1.3解算参考谱段B2相对于非参考谱段B1的几何畸变参数的实现方式如 下,
令
上式中,矢量
根据非参考谱段B1的严格几何成像模型,得到下式
令
上式中,vxi和vyi分别代表沿轨和垂轨方向的像方残差;
将步骤1.2中所得物方点Pi坐标代入上式中,对每个物方点 Pi构建如如下误差方程,
Vi=AiX-Li,Wi,
其中,
上式中,Vi、Ai、Li分别是利用物方点Pi构建的误差方程,的残差向量、待解参数的系 数矩阵以及常向量;X代表参考谱段B2相对于非参考谱段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)所示:
其中,为三次多项式的系数,s代表探元号。
将B1谱段相对于B2谱段的几何畸变参数记为(上标1代表谱段B1),则:
利用卫星获取的姿态、轨道、时间、相机参数以及物方高程信息,基于严格几何成像模 型,能够实现像点坐标与物方点坐标之间的换算。以谱段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)令
上式中,矢量
式(2)可转化为式(10):
令
上式中,vxi和vyi分别代表沿轨和垂轨方向的像方残差。
2)将步骤1.2中解算的物方点Pi坐标(i=1,…,n)代入式(11)中,对每 个物方点Pi均可构建如式(12)所示的误差方程(下标i代表利用物方点Pi建立的误差方程):
Vi=AiX-Li,Wi(i=1,…,n) (12)
其中,
上式中,Vi、Ai、Li分别是利用物方点Pi构建的误差方程式的残差向量、待解参数的系 数矩阵以及常向量;X代表B1谱段相对于B2谱段的几何畸变参数
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两谱段影像的配准。其中, 灰度重采样可采用现有技术中的双线性内插算法实现,本发明不予赘述。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技 术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不 会偏离本发明的精神或者超越所附权利要求书所定义的范围。
机译: 基于通过卫星定位电路获得的另一方向信息的电子来校正通过传感器获得的传感器和方向信息的方法及其电子设备
机译: 卫星导航系统全球定位系统,一种减少接收机干扰的方法,涉及基于本征特征反射来影响残余误差,这些特征反射是由于检测器类型特征而被推断出来的
机译: 一种方法是在基于热塑聚合物的模制件之间获得牢固的接头,另一方则是包含羧基的EP(D)橡胶,以及根据此方法获得的复合材料。