一种光学卫星传感器在轨分时几何定标方法技术

本发明专利技术公开了一种光学卫星传感器在轨分时几何定标方法，包括：以光学卫星传感器的中间片CCD/CMOS作为基准片，以其余片CCD/CMOS作为非基准片，依次对基准片进行外定标和内定标处理；在基准片外定标的基础上，分别对每一片非基准片进行内定标处理；从基准片开始，分别向左右两侧CCD/CMOS延伸，依次以相邻两片CCD/CMOS为一个定标单元，并以该单元中离基准片近的CCD/CMOS作为参考片、另一片作为非参考片，计算参考片和非参考片之间的片间错位量；根据片间错位量和成像探元大小，更新非参考片的内定标参数。本发明专利技术可以大大减小光学卫星传感器在轨几何定标对地面定标场幅宽的要求，从而减少定标场参考数据制作与更新过程中投入的人力、物力和财力，实现经济化在轨几何定标。实现经济化在轨几何定标。实现经济化在轨几何定标。

【技术实现步骤摘要】
一种光学卫星传感器在轨分时几何定标方法


[0001]本专利技术属于光学遥感卫星数据几何处理
，特别涉及一种光学卫星传感器在轨分时几何定标方法。

技术介绍

[0002]光学卫星传感器在轨几何定标旨在于获取传感器的精确成像参数，包括相机主点、主距、镜头畸变和相机安置角等，是光学卫星传感器实现多片/多相机无缝拼接、全色与多光谱影像精确配准、多光谱波段影像精确配准、高精度目标定位的前提和基础。
[0003]针对光学卫星传感器高精度在轨几何定标的应用需求，目前通常采用基于地面定标场的在轨几何定标方法，即：通过覆盖定标场的待定标影像和定标场参考DOM/DEM数据的密集匹配来获得大量控制点，并利用控制点精确求解出定标参数。该方法要求所有分片CCD/CMOS影像为同一时间段内采集的影像，且要求定标场参考数据在卫星影像幅宽方向上覆盖所有片CCD/CMOS影像。然而，随着卫星传感器技术的不断发展和卫星用户需求的不断提高，卫星影像幅宽已经有了大幅度增大。譬如，高分六号卫星宽幅相机的成像幅宽为800公里，海洋一号C/D卫星海岸带成像仪的成像幅宽达950公里。
[0004]在此情况下，现有定标方法必然面临两个主要问题：
[0005]1)现有定标场的参考数据无法满足大幅宽成像的在轨几何定标需求，而制作大幅宽定标场的参考数据需要投入大量的人力、物力和财力；
[0006]2)采集覆盖定标场的待定标影像受天气状况影响较大，难以保证同一时间段内所有片CCD/CMOS影像均无云或少云覆盖。

技术实现思路
/>[0007]针对现有方法在光学卫星传感器在轨几何定标方面的不足，本专利技术提供了一种光学卫星传感器在轨分时几何定标方法，利用不同时间段采集的覆盖同一定标场或者不同定标场的分片CCD/CMOS影像分别开展几何定标处理，再通过相邻片CCD/CMOS影像之间的片间错位量和成像探元大小，更新非参考片的内定标参数，最终实现所有片CCD/CMOS的高精度定标。
[0008]为进一步实现上述目的，本专利技术采用以下技术方案：
[0009]本专利技术的技术构思是：在光学卫星传感器基准片内外定标和非基准片内定标的基础上，以相邻两片CCD/CMOS作为一个定标单元，通过定标单元中参考片与非参考片之间的片间错位量和成像探元大小，对非参考片的内定标参数进行更新处理，从而在保证每一片CCD/CMOS定标精度的同时，实现所有片CCD/CMOS之间的几何一致性。
[0010]一种光学卫星传感器在轨分时几何定标方法，包括以下步骤：
[0011]步骤1，以光学卫星传感器的中间片CCD/CMOS作为基准片，以其余片CCD/CMOS作为非基准片，利用基准片影像上的控制点依次对基准片进行外定标和内定标处理，采用的在轨几何定标模型如式(1)和式(2)所示：
[0012][0013][0014]式中：(X
GPS
,Y
GPS
,Z
GPS
)为GPS天线相位中心在WGS 84坐标系下的空间坐标；(X,Y,Z)为地面点在WGS 84坐标系下的物方空间坐标；为WGS84坐标系至J2000坐标系的旋转矩阵；为J2000坐标系至姿态测量坐标系的旋转矩阵；为传感器在姿态测量坐标系下的安置矩阵；λ为比例因子；(ψ
x
,ψ
y
)为地面点对应的成像探元在传感器坐标系下的指向角；S为成像探元编号；(a0,a1,a2,a3,b0,b1,b2,b3)为内定标参数；(p,r,y)为外定标参数；
[0015]步骤2，以步骤1获得的外定标参数作为非基准片的外定标参数，利用非基准片影像上的控制点，分别对每一片非基准片进行内定标处理；其中，每一片非基准片影像和基准片影像是分时成像的，即：非基准片影像和基准片影像为不同时间段采集的覆盖同一定标场或者不同定标场的影像，可以是同一轨道内的不同时间段成像，也可以是不同轨道之间的不同时间段成像；
[0016]步骤3，从基准片开始，分别向左右两侧CCD/CMOS延伸，依次以相邻两片CCD/CMOS为一个定标单元，并以该单元中离基准片近的CCD/CMOS作为参考片，另一片作为非参考片，计算参考片和非参考片之间的片间错位量，包括：
[0017]步骤3.1，分别计算参考片获取的影像和非参考片获取的影像的有理多项式系数，具体为：
[0018]步骤3.1.1，根据参考片和非参考片的定标参数，分别构建严密成像几何模型，除基准片所在定标单元外，当前定标单元中的参考片为前一个定标单元中的非参考片，因此当前定标单元中参考片的内定标参数为前一个定标单元中更新后的非参考片内定标参数；
[0019]步骤3.1.2，根据严密成像几何模型，分别生成虚拟控制点；
[0020]步骤3.1.3，根据虚拟控制点，分别求解参考影像和非参考影像的有理多项式系数；
[0021]步骤3.2，计算参考影像和非参考影像之间的重叠区域，具体为：
[0022]步骤3.2.1，根据非参考影像的有理多项式系数，将非参考影像的四个角点分别投影至物方，得到物方投影点；
[0023]步骤3.2.2，根据参考影像的有理多项式系数，分别将四个物方投影点进一步投影至参考影像上，得到像方投影点；
[0024]步骤3.2.3，计算参考影像四个角点构成的四边形和四个投影点构成的四边形之间的交集，得到参考影像上的重叠区域；
[0025]步骤3.3，参考影像和非参考影像之间匹配同名像点：
[0026]步骤3.3.1，分别在参考影像重叠区域内的行和列方向上，按照等像素间隔划分规则格网，得到若干格网点；
[0027]步骤3.3.2，在参考影像和非参考影像之间进行相关系数匹配，得到参考影像上的格网点在非参考影像上的初始同名点；
[0028]步骤3.3.3，在初始同名点的基础上，进行最小二乘影像匹配，得到参考影像和非参考影像之间的精确同名点；
[0029]步骤3.4，计算参考影像和非参考影像之间的片间错位量：
[0030]步骤3.4.1，根据参考影像的有理多项式系数，将参考影像上的每一个同名点分别投影至物方，得到物方投影点；
[0031]步骤3.4.2，根据非参考影像的有理多项式系数，进一步将物方投影点分别投影至非参考影像上，得到像方投影点；
[0032]步骤3.4.3，分别计算非参考影像上每个投影点的像方坐标与其对应匹配点的像方坐标的差值(Δx,Δy)，如式(3)所示：
[0033][0034]式中，(x
′
i
,y
′
i
)和(x
i
,y
i
)分别为投影点像方坐标及其对应匹配点像方坐标，i＝1,2,...,n，n为同名像点个数；
[0035]步骤3.4.4，计算步骤3.4.3中所有坐标差值的平均值，得到片间错位量；
[0036]步骤4，更新非参考片的内定标参数；
[0037]步骤4.本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种光学卫星传感器在轨分时几何定标方法，其特征在于，包括以下步骤：步骤1，以光学卫星传感器的中间片CCD/CMOS作为基准片，以其余片CCD/CMOS作为非基准片，利用基准片影像上的控制点依次对基准片进行外定标和内定标处理，采用的在轨几何定标模型如式(1)和式(2)所示：何定标模型如式(1)和式(2)所示：式中：(X
GPS
,Y
GPS
,Z
GPS
)为GPS天线相位中心在WGS 84坐标系下的空间坐标；(X,Y,Z)为地面点在WGS 84坐标系下的物方空间坐标；为WGS84坐标系至J2000坐标系的旋转矩阵；为J2000坐标系至姿态测量坐标系的旋转矩阵；为传感器在姿态测量坐标系下的安置矩阵；λ为比例因子；(ψ
x
,ψ
y
)为地面点对应的成像探元在传感器坐标系下的指向角；S为成像探元编号；(a0,a1,a2,a3,b0,b1,b2,b3)为内定标参数；(p,r,y)为外定标参数；步骤2，以步骤1获得的外定标参数作为非基准片的外定标参数，利用非基准片影像上的控制点，分别对每一片非基准片进行内定标处理；其中，每一片非基准片影像和基准片影像是分时成像的，即：非基准片影像和基准片影像为不同时间段采集的覆盖同一定标场或者不同定标场的影像，可以是同一轨道内的不同时间段成像，也可以是不同轨道之间的不同时间段成像；步骤3，从基准片开始，分别向左右两侧CCD/CMOS延伸，依次以相邻两片CCD/CMOS为一个定标单元，并以该单元中离基准片近的CCD/CMOS作为参考片，另一片作为非参考片，计算参考片和非参考片之间的片间错位量，包括：步骤3.1，分别计算参考片获取的影像和非参考片获取的影像的有理多项式系数，具体为：步骤3.1.1，根据参考片和非参考片的定标参数，分别构建严密成像几何模型，除基准片所在定标单元外，当前定标单元中的参考片为前一个定标单元中的非参考片，因此当前定标单元中参考片的内定标参数为前一个定标单元中更新后的非参考片内定标参数；步骤3.1.2，根据严密成像几何模型，分别生成虚拟控制点；步骤3.1.3，根据虚拟控制点，分别求解参考影像和非参考影像的有理多项式系数；步骤3.2，计算参考影像和非参考影像之间的重叠区域...

【专利技术属性】
技术研发人员：曹金山，叶志伟，杨航，王志林，刘熊峰，
申请(专利权)人：湖北工业大学，
类型：发明
国别省市：