CN113900122B - 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 - Google Patents
星地联合的面阵成像遥感卫星全自主几何定标方法及系统 Download PDFInfo
- Publication number
- CN113900122B CN113900122B CN202111133176.0A CN202111133176A CN113900122B CN 113900122 B CN113900122 B CN 113900122B CN 202111133176 A CN202111133176 A CN 202111133176A CN 113900122 B CN113900122 B CN 113900122B
- Authority
- CN
- China
- Prior art keywords
- satellite
- calibration
- model
- cmos
- star
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 10
- 239000012482 calibration solution Substances 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 57
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000012937 correction Methods 0.000 claims description 24
- 239000000243 solution Substances 0.000 claims description 23
- 239000000523 sample Substances 0.000 claims description 18
- 230000003287 optical effect Effects 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 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 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000012634 fragment Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012876 topography Methods 0.000 abstract description 2
- 230000009897 systematic effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
Abstract
本发明借助高分辨率面阵成像卫星的敏捷机动特性,提出了一种星地联合的面阵成像遥感卫星全自主几何定标方法及系统,实现从数据获取到内外定标解算的全链路处理。本发明基于卫星的超高姿态稳定度和敏捷机动能力,在同一圈的太阳阴影区和阳照区通过多次机动获取对天和对地观测数据,通过阴影区的恒星观测数据完成卫星的外定标,通阳照区获取的地面重叠影像,在地形辅助下完成载荷的内定标,进而完全摆脱对地面定标场参考数据的依赖,不但可降低定标成本,且几乎可随时随地定标,保证定标的时效性和精度。
Description
技术领域
本发明属于光学卫星遥感影像几何处理领域,涉及一种星地联合的面阵成像遥感卫星全自主几何定标方法及系统。
背景技术
几何定标是改正光学遥感卫星成像模型中系统性几何误差的必要手段,虽然卫星载荷在发射前已经进行了严格的实验室标定,但由于卫星发射过程中成像环境变化以及相机调焦等因素,载荷的真实成像参数势必发生改变,需要在轨重新进行几何定标才能保证其影像的几何精度。传统在轨几何定标方法利用地面控制场的高精度参考影像(数字高程模型和数字正射影像)来补偿载荷几何成像模型中的系统误差,但随着卫星影像分辨率的不断提高,以及各行业对卫星影像处理的精度和时效性要求日益提高,这种方法在实际处理中的弊端日益凸显:1)定标场昂贵的建设和维护成本和较短的有效周期(地物变化)提高了在轨几何定标的代价;2)卫星较长的重访周期严重限制其在轨获取定标数据的时间窗口,且过境定标影像的采集受天气影响较大,造成传统定标处理所需时间过长,甚至影响卫星的业务化使用。3)采用航空摄影测量的方式生产的定标场数据,已越来越无法满足超高分辨率光学遥感卫星在轨几何定标的分辨率和精度需求。
发明内容
本发明要解决的问题是面阵成像高分辨率光学遥感卫星高精度定位的关键问题。本发明借助高分辨率面阵成像卫星的敏捷机动特性,提出一种星地联合的面阵成像遥感卫星全自主几何定标方法及系统,实现从数据获取到内外定标解算的全链路处理。本发明基于卫星的超高姿态稳定度和敏捷机动能力,在同一圈的太阳阴影区和阳照区通过多次机动获取对天和对地观测数据,通过阴影区的恒星观测数据完成卫星的外定标,通阳照区获取的地面重叠影像,在地形辅助下完成载荷的内定标,进而完全摆脱对地面定标场参考数据的依赖,不但可降低定标成本,且几乎可随时随地定标,保证定标的时效性和精度。
本发明的技术方案为一种星地联合的面阵成像遥感卫星全自主几何定标方法,包含以下步骤:
步骤1,当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
步骤2,基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
步骤3,基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
步骤4,当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
步骤5,基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
步骤6,在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
步骤7,在内定标解算参数的基础上,重复步骤2-3,重新进行外定标解算,进而得到内外定标结果。
进一步的,步骤1中当卫星运动到阴影区时,执行地相机的对天观测,由于不同的天区包含的星数以及恒星星分布情况均有所不同,在进行对天观测时需要考虑天区内恒星的观测情况;对全天区依据相机视场角划分网格,根据相机的敏感度选择天区内恒星成像情况和分布情况均最佳的天区作为观测天区,调整卫星姿态,对该天区成像,获取序列观星影像数据。
在导航星表的辅助下,通过影像去噪、二值化、边缘提取与质心拟合等图像处理方法提取精确的恒星观测值,得到每个恒星的赤经α和赤纬δ。
进一步的,步骤2的具体实现方式如下;
恒星观测值在天球坐标系下的观测向量为vstar=(cosαcosδ,sinαcosδ,sinδ)T,进而基于恒星观测向量和卫星相机视向量的共线关系,建立面阵光学卫星影像对天成像的严格几何成像模型,如下式:
其中,α和δ为恒星点的赤经和赤纬,和为像点在行列两个方向的指向角,μ为缩放系数;为J2000坐标系到卫星本体坐标系的旋转矩阵,由星敏、陀螺通过组合定姿得到;RAber为光行差改正矩阵,该矩阵可根据光行差θ和旋转向量n计算的变换四元数P=(p0,p1,p2,p3)确定:
其中,θ=v·sinβ/c,v为卫星运行速度,c为光速,β为恒星方向向量vstar与卫星运行方向向量vsat的夹角。
构建的外定标模型如公式(4)所示,代表从卫星本体坐标系到相机坐标系的广义旋转矩阵,将所有的外方位元素误差统一纳入到该矩阵中进行补偿,该矩阵由三个相机安装角确定,具体如下式:
进一步的,步骤3的具体实现方式如下;
31)构建外定标平差模型
采用最小二乘平差方法解算外定标参数,首先基于外定标模型(4)建立用于参数解算的平差模型,令式(6)中:
构建用于最小二乘平差解算的平差方程:
32)外定标参数解算
针对每个恒星观测值,根据其成像时间、姿态和轨道建立用于外参数解算的平差方程,并在此基础上,通过模型线性化得到误差方程:
其中,为三个外参数的改正向量;Ai是个恒星控制点的误差方程系数矩阵;Li为恒星控制点误差方程的常数向量,等于该点对应的平差方程的当前值,具体如下:
利用最小二乘平差计算x,如式(10);
其中,m为恒星观测值数量;根据解算的改正值更新外参数的当前值,当两次平差解算的结果小于限差时,迭代结束。
进一步的,步骤4中面阵卫星载荷是由多片CMOS组成,且各CMOS的畸变情况也不尽相同,因此在内定标中各片CMOS具有不同的待定标参数,在自主定标中对各片CMOS分开处理,因此需要每片CMOS获取的影像,在行列两个方向具有一定的重叠度Tcmos-l和Tcmos-s,要求两个方向的重叠度在45%到70%之间,由此可得到整景影像间的重叠度Tl和Ts与各片CMOS影像间重叠度的关系:
然后,在三景影像的两个重叠区域采用成熟的基于SIFT算子的匹配方法匹配密集同名点数据,作为内定标解算的观测值。
进一步的,步骤5中的具体实现方式如下;
51)面阵光学卫星对地成像模型构建
基于影像的成像时间,得到影像的姿轨辅助参数,然后基于像点、物点和投影中心三点共线的几何关系建立面阵光学卫星对地观测的严格几何成像模型,如下式:
其中:仍为卫星本体坐标系到相机坐标系的旋转矩阵,由上述外定标确定,为WGS84坐标系到J2000坐标系的旋转矩阵,由成像时刻的星历参数确定,其他矩阵与对天成像模型相同,(Xgps,Ygps,Zgps)表示GPS天线相位中心在WGS84坐标系下的坐标,由卫星上搭载的GPS获取;(Xg,Yg,Zg)表示像点对应的物方点在WGS84坐标系下的直角坐标,与地理坐标(Lat,Lon,Hei)(维度,精度,高程)之间的变换关系如下:
其中,N是地球卯酉圈的曲率半径,e是地球椭球体的第一个偏心率;
52)基于探元指向角的内定标模型
对于面阵相机的每个CMOS探元的指向,可利用其在相机坐标系下的指向角精确表示,但计算所有探元的指向角是不现实的,根据卫星载荷的畸变特性,通常利用两个二元三次多项式对其各探元在相机坐标系下的指向角进行拟合,进而得到基于二元三次多项式拟合的指向角内定标模型的:
其中,(s,l)为CMOS探元号,(a0,…a9,b0,…b9)为三次多项式系数,即为在轨几何内定标参数;
将构建的各片CMOS的内定标模型引入到严格几何成像模型(12)中,可得到其在轨几何定标模型,其中每片CMOS均具有其对应的一套待解算的内定标模型参数。
进一步的,步骤6的具体实现方式如下;
61)构建自主内定标平差模型
采用最小二乘平差方法解算内定标参数,因此首先基于几何定标模型建立用于参数解算的平差模型,令式(12)中:
构建用于最小二乘平差解算的平差方程:
其中,不同片CMOS的模型有平差模型中各自的指向角模型决定;
62)逐CMOS分片定标参数解算
对于每片CMOS,利用其获取的重叠影像对上的密集同名点整体解算其定标参数,由于指向角模型常数项(a0,b0)独立于影像间的共面约束条件,无法在内定标中精确解算,而该参数又与外定标参数完全相关,因此,在进行每个CMOS内定标解算时直接忽略常数项,只计算其他高次项参数;鉴于每片CMOS高阶内定标参数的解算过程相同,仅对一片的参数解算进行叙述,具体如下:
根据构建的平差模型(16),通过模型线性化处理,对每对同名像点(左右影像上各一个像点)可构建误差方程,如式(17)所示:
其中,和分别为左右影像上像点对应的改正数向量,y=[da1,…da9,db1,…db9]T为相机内定标参数改正向量;表示各同名像点物方平面坐标改正向量,物方坐标的高程直接从影像覆盖区域的DSM中插值得到,而不在内定标中解算;矩阵和分别为左右像点误差方程中定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量。以左影像上像点为例,其误差方程中各矩阵具体形式如下:
最后,利用最小二乘平差计算y,如式(18)所示:
其中,k代表该CMOS对应的重叠影像上同名点的数量;
同样地,基于最小二乘的内定标解算也为迭代计算,根据每次迭代解算的结果更新当前的内定标参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算;
63)常数项整体解算
为了确保各片CMOS间的拼接精度,需要根据相邻CMOS重叠影像间的同名像点解算各片CMOS的常数项,这里选择某一片CMOS为基准片,基准片的常数项仍不计算,以基准片为参考,计算所有非基准片的常数项;
基于上述构建的平差方程对同名点分别构建误差方程,方程的具体形式与(17)类似,不同之处在于这里的待解算参数为所有非基准CMOS指向角模型的常数项,具体如下式:
其中,z=[dz1,dz2…dzn]T为所有非基准CMOS的常数项改正数,其中dzi=[da0,db0]i,n为非基准CMOS个数;ti仍为各同名像点物方平面坐标改正向量,矩阵和分别为左右像点误差方程中常数项定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,同样是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量。
最后,利用最小二乘平差计算z,如式(18)所示:
其中,λ代表用于计算的同名点数量;
同样地,该解算为迭代计算,根据每次迭代解算的结果更新当前的内定标常数项参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算。
本发明还提供一种星地联合的面阵成像遥感卫星全自主几何定标系统,包括如下模块:
外定标观测值提取模块,用于当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
外定标模型构建模块,用于基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
外定标参数解算模块,用于基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
内定标观测值提取模块,用于当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
内定标模型构建模块,用于基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
内定标参数计算模块,用于在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
定标结果获取模块,用于在内定标解算参数的基础上,重复外定标模型构建模块和外定标参数解算模块,重新进行外定标解算,进而得到内外定标结果。
本发明的优点在于:通过基于卫星平台敏捷机动能力的星地联合处理模式,在完全无需地面控制点数据的条件下即可实现面阵高分辨率光学遥感卫星的全自主在轨几何定标,完全摆脱传统方法对地面定标场的依赖,该在轨几何定标方法成本低、时效性高、精度高。
附图说明
图1为本发明实施例的流程示意图;
图2为本发明中利用卫星机动能力的对天观测示意图;
图3为本发明中利用卫星机动能力的对地观测获取重叠影像示意图;
图4为本发明中地面重叠影像及其同名点分布图。
具体实施方式
以下结合附图和实施例详细说明本发明具体实施方式。本发明基于卫星的机动能力,通过基于恒星的外定标与基于地面重叠影像自约束的内定标相结合的星地联合方法实现面阵光学遥感卫星成像链路内外部系统性几何误差的补偿。参见图1,本发明的实施流程可以分为如下7个步骤:
步骤1,当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
步骤2,基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
步骤3,基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
步骤4,当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
步骤5,基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
步骤6,在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
步骤7,在内定标解算参数的基础上,重复步骤2-3,重新进行外定标解算,进而得到内外定标结果。
具体的,每个步骤的具体方法、公式以及流程如下:
步骤1,基于卫星敏捷机动的太阳阴影区恒星观测与恒星提取
当卫星运动到阴影区时,执行地相机的对天观测(如附图2)。由于不同的天区包含的星数以及恒星星等分布情况均有所不同,在进行对天观测时需要考虑天区内恒星的观测情况。对全天区依据相机视场角划分网格,根据相机的敏感度选择天区内恒星成像情况和分布情况均最佳的天区作为观测天区。调整卫星姿态,对该天区成像,获取序列观星影像数据。
在导航星表的辅助下,通过影像去噪、二值化、边缘提取与质心拟合等图像处理方法提取精确的恒星观测值,得到每个恒星的赤经α和赤纬δ。
步骤2,构建基于广义安装角的在轨几何外定标模型
21)面阵光学卫星对天成像模型构建
恒星观测值在天球坐标系下的观测向量为vstar=(cosαcosδ,sinαcosδ,sinδ)T,进而基于恒星观测向量和卫星相机视向量的共线关系,建立面阵光学卫星影像对天成像的严格几何成像模型,如下式:
其中,和为像点在行列两个方向的指向角,μ为缩放系数;为J2000坐标系到卫星本体坐标系的旋转矩阵,由星敏、陀螺通过组合定姿得到;RAber为光行差改正矩阵,该矩阵可根据光行差θ和旋转向量n计算的变换四元数P=(p0,p1,p2,p3)确定:
其中,θ=v·sinβ/c,v为卫星运行速度,c为光速,β为恒星方向向量vstar与卫星运行方向向量vsat的夹角。
22)外定标模型
本发明构建的外定标模型如公式(4)所示,代表从卫星本体坐标系到相机坐标系的广义旋转矩阵,本发明将所有的外方位元素误差统一纳入到该矩阵中进行补偿,该矩阵由三个相机安装角确定,具体如下式:
步骤3,基于序列恒星观测值的外定标参数解算
31)构建外定标平差模型
采用最小二乘平差方法解算外定标参数,首先基于几何外定标模型(4)建立用于参数解算的平差模型,令式(6)中:
构建用于最小二乘平差解算的平差方程:
32)外定标参数解算
针对每个恒星观测值,根据其成像时间、姿态和轨道建立用于外参数解算的平差方程,并在此基础上,通过模型线性化得到误差方程:
其中,为三个外参数的改正向量;Ai是个恒星控制点的误差方程系数矩阵;Li为恒星控制点误差方程的常数向量,等于该点对应的平差方程的当前值,具体如下:
利用最小二乘平差计算x,如式(10);
其中,m为恒星观测值数量。根据解算的改正值更新外参数的当前值,当两次平差解算的结果小于限差时,迭代结束。
步骤4,基于卫星敏捷机动的太阳阳照区地面重叠影像获取与同名点匹配
当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的卫星影像,如附图3,得到两次重叠影像如附图4所示,由于面阵卫星载荷可能是由多片CMOS组成,且各CMOS的畸变情况也不尽相同,因此在内定标中各片CMOS具有不同的待定标参数,在自主定标中也需要对各片CMOS分开处理,因此需要每片CMOS获取的影像,在行列两个方向具有一定的重叠度Tcmos-l和Tcmos-s,本发明要求两个方向的重叠度在45%到70%之间,最佳为55%,由此可得到整景影像间的重叠度Tl和Ts与各片CMOS影像间重叠度的关系:
然后,在三景影像的两个重叠区域采用成熟的基于SIFT算子的匹配方法匹配密集同名点数据,作为内定标解算的观测值。
步骤5,基于重叠影像共面约束的在轨几何内定标模型构建
51)面阵光学卫星对地成像模型构建
基于影像的成像时间,得到影像的姿轨辅助参数,然后基于像点、物点和投影中心三点共线的几何关系建立面阵光学卫星对地观测的严格几何成像模型,如下式:
其中:仍为卫星本体坐标系到相机坐标系的旋转矩阵,由上述外定标确定,为WGS84坐标系到J2000坐标系的旋转矩阵,由成像时刻的星历参数确定,其他矩阵与上述对天成像模型相同,(Xgps,Ygps,Zgps)表示GPS天线相位中心在WGS84坐标系下的坐标,由卫星上搭载的GPS获取;(Xg,Yg,Zg)表示像点对应的物方点在WGS84坐标系下的直角坐标,与地理坐标(Lat,Lon,Hei)(维度,精度,高程)之间的变换关系如下:
其中,N是地球卯酉圈的曲率半径,e是地球椭球体的第一个偏心率。
52)基于探元指向角的内定标模型
对于面阵相机的每个CMOS探元的指向,可利用其在相机坐标系下的指向角精确表示,但计算所有探元的指向角是不现实的,根据卫星载荷的畸变特性,通常利用两个二元三次多项式对其各探元在相机坐标系下的指向角进行拟合,进而得到基于二元三次多项式拟合的指向角内定标模型的:
其中,(s,l)为CMOS探元号,(a0,…a9,b0,…b9)为三次多项式系数,即为在轨几何内定标参数。
将构建的各片CMOS的内定标模型引入到严格几何成像模型(12)中,可得到其在轨几何定标模型,其中每片CMOS均具有其对应的一套待解算的内定标模型参数。
步骤6,基于重叠影像共面约束的内定标参数解算
在解算的外参数的基础上,进行内定标参数解算,在内定标中将上述解算的外定标参数视为真值。
61)构建自主内定标平差模型
采用最小二乘平差方法解算内定标参数,因此首先基于几何定标模型建立用于参数解算的平差模型,令式(12)中:
构建用于最小二乘平差解算的平差方程:
其中,不同片CMOS的模型有平差模型中各自的指向角模型决定。
62)逐CMOS分片定标参数解算
对于每片CMOS,利用其获取的重叠影像对(两组重叠影像对)上的密集同名点整体解算其定标参数,由于指向角模型常数项(a0,b0)独立于影像间的共面约束条件,无法在内定标中精确解算,而该参数又与外定标参数完全相关(可用外参数补偿其误差),因此,在进行每个CMOS内定标解算时直接忽略常数项,只计算其他高次项参数。鉴于每片CMOS高阶内定标参数的解算过程相同,仅对一片的参数解算进行叙述,具体如下:
根据构建的平差模型(16),通过模型线性化处理,对每对同名像点(左右影像上各一个像点)可构建误差方程,如式(17)所示:
其中,和分别为左右影像上像点对应的改正数向量,y=[da1,…da9,db1,…db9]T为相机内定标参数改正向量;表示各同名像点物方平面坐标改正向量,物方坐标的高程直接从影像覆盖区域的DSM中插值得到,而不在内定标中解算;矩阵和分别为左右像点误差方程中定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量。以左影像上像点为例,其误差方程中各矩阵具体形式如下:
最后,利用最小二乘平差计算y,如式(18)所示:
其中,k代表该CMOS对应的重叠影像上同名点的数量。
同样地,基于最小二乘的内定标解算也为迭代计算,根据每次迭代解算的结果更新当前的内定标参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算。每片CMOS解算过程相同,这里不再赘述。
63)常数项整体解算
为了确保各片CMOS间的拼接精度,需要根据相邻CMOS重叠影像间的同名像点解算各片CMOS的常数项,这里选择某一片CMOS为基准片,基准片的常数项仍不计算(采用其初值),以基准片为参考,计算所有非基准片的常数项。
基于上述构建的平差方程对同名点分别构建误差方程,方程的具体形式与(17)类似,不同之处在于这里的待解算参数为所有非基准CMOS指向角模型的常数项,具体如下式:
其中,z=[dz1,dz2…dzn]T为所有非基准CMOS的常数项改正数,其中dzi=[da0,db0]i,n为非基准CMOS个数;ti仍为各同名像点物方平面坐标改正向量,矩阵和分别为左右像点误差方程中常数项定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,同样是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量。
最后,利用最小二乘平差计算z,如式(18)所示:
其中,λ代表用于计算的同名点数量。
同样地,该解算为迭代计算,根据每次迭代解算的结果更新当前的内定标常数项参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算。
步骤7,在解算的内定标参数基础上,重复步骤2-3,重新进行外定标参数解算,得到新的外定标结果,即完成所有的外定标和内定标解算。
本发明实施例还提供一种星地联合的面阵成像遥感卫星全自主几何定标系统,包括如下模块:
外定标观测值提取模块,用于当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
外定标模型构建模块,用于基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
外定标参数解算模块,用于基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
内定标观测值提取模块,用于当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
内定标模型构建模块,用于基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
内定标参数计算模块,用于在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
定标结果获取模块,用于在内定标解算参数的基础上,重复外定标模型构建模块和外定标参数解算模块,重新进行外定标解算,进而得到内外定标结果。
各模块的具体实现方式和各步骤相应,本发明不予撰述。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于,包括如下步骤:
步骤1,当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
步骤2,基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
步骤3,基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
步骤4,当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
步骤5,基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
步骤6,在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
步骤7,在内定标解算参数的基础上,重复步骤2-3,重新进行外定标解算,进而得到内外定标结果。
2.如权利要求1所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤1中当卫星运动到阴影区时,执行地相机的对天观测,由于不同的天区包含的星数以及恒星星分布情况均有所不同,在进行对天观测时需要考虑天区内恒星的观测情况;对全天区依据相机视场角划分网格,根据相机的敏感度选择天区内恒星成像情况和分布情况均最佳的天区作为观测天区,调整卫星姿态,对该天区成像,获取序列观星影像数据;
在导航星表的辅助下,通过影像去噪、二值化、边缘提取与质心拟合图像处理方法提取精确的恒星观测值,得到每个恒星的赤经α和赤纬δ。
3.如权利要求2所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤2的具体实现方式如下;
恒星观测值在天球坐标系下的观测向量为vstar=(cosαcosδ,sinαcosδ,sinδ)T,进而基于恒星观测向量和卫星相机视向量的共线关系,建立面阵光学卫星影像对天成像的严格几何成像模型,如下式:
其中,α和δ为恒星点的赤经和赤纬,和为像点在行列两个方向的指向角,μ为缩放系数;为J2000坐标系到卫星本体坐标系的旋转矩阵,由星敏、陀螺通过组合定姿得到;RAber为光行差改正矩阵,该矩阵可根据光行差θ和旋转向量n计算的变换四元数P=(p0,p1,p2,p3)确定:
其中,θ=v·sinβ/c,v为卫星运行速度,c为光速,β为恒星方向向量vstar与卫星运行方向向量vsat的夹角;
构建的外定标模型如公式(4)所示,代表从卫星本体坐标系到相机坐标系的广义旋转矩阵,将所有的外方位元素误差统一纳入到该矩阵中进行补偿,该矩阵由三个相机安装角确定,具体如下式:
4.如权利要求3所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤3的具体实现方式如下;
31)构建外定标平差模型
采用最小二乘平差方法解算外定标参数,首先基于外定标模型(4)建立用于参数解算的平差模型,设置关系式如式(7):
构建用于最小二乘平差解算的平差方程:
32)外定标参数解算
针对每个恒星观测值,根据其成像时间、姿态和轨道建立用于外参数解算的平差方程,并在此基础上,通过模型线性化得到误差方程:
其中,为三个外参数的改正向量;Ai是个恒星控制点的误差方程系数矩阵;Li为恒星控制点误差方程的常数向量,等于该点对应的平差方程的当前值,具体如下:
利用最小二乘平差计算x,如式(10);
其中,m为恒星观测值数量;根据解算的改正值更新外参数的当前值,当两次平差解算的结果小于限差时,迭代结束。
5.如权利要求1所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤4中面阵卫星载荷是由多片CMOS组成,且各CMOS的畸变情况也不尽相同,因此在内定标中各片CMOS具有不同的待定标参数,在自主定标中对各片CMOS分开处理,因此需要每片CMOS获取的影像,在行列两个方向具有一定的重叠度Tcmos-l和Tcmos-s,要求两个方向的重叠度在45%到70%之间,由此可得到整景影像间的重叠度Tl和Ts与各片CMOS影像间重叠度的关系:
其中,ns,nl为两个方向的CMOS数量;然后,在三景影像的两个重叠区域采用成熟的基于SIFT算子的匹配方法匹配密集同名点数据,作为内定标解算的观测值。
6.如权利要求4所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤5中的具体实现方式如下;
51)面阵光学卫星对地成像模型构建
基于影像的成像时间,得到影像的姿轨辅助参数,然后基于像点、物点和投影中心三点共线的几何关系建立面阵光学卫星对地观测的严格几何成像模型,如下式:
其中:仍为卫星本体坐标系到相机坐标系的旋转矩阵,由上述外定标确定,为WGS84坐标系到J2000坐标系的旋转矩阵,由成像时刻的星历参数确定,其他矩阵与对天成像模型相同,(Xgps,Ygps,Zgps)表示GPS天线相位中心在WGS84坐标系下的坐标,由卫星上搭载的GPS获取;(Xg,Yg,Zg)表示像点对应的物方点在WGS84坐标系下的直角坐标,与地理坐标(Lat,Lon,Hei)(维度,精度,高程)之间的变换关系如下:
其中,N是地球卯酉圈的曲率半径,e是地球椭球体的第一个偏心率;
52)基于探元指向角的内定标模型
对于面阵相机的每个CMOS探元的指向,可利用其在相机坐标系下的指向角精确表示,但计算所有探元的指向角是不现实的,根据卫星载荷的畸变特性,通常利用两个二元三次多项式对其各探元在相机坐标系下的指向角进行拟合,进而得到基于二元三次多项式拟合的指向角内定标模型的:
其中,(s,l)为CMOS探元号,(a0,…a9,b0,…b9)为三次多项式系数,即为在轨几何内定标参数;
将构建的各片CMOS的内定标模型引入到严格几何成像模型(12)中,可得到其在轨几何定标模型,其中每片CMOS均具有其对应的一套待解算的内定标模型参数。
7.如权利要求6所述的星地联合的面阵成像遥感卫星全自主几何定标方法,其特征在于:步骤6的具体实现方式如下;
61)构建自主内定标平差模型
采用最小二乘平差方法解算内定标参数,因此首先基于几何定标模型建立用于参数解算的平差模型,设置关系式如式(15):
构建用于最小二乘平差解算的平差方程:
其中,不同片CMOS的模型有平差模型中各自的指向角模型决定;
62)逐CMOS分片定标参数解算
对于每片CMOS,利用其获取的重叠影像对上的密集同名点整体解算其定标参数,由于指向角模型常数项(a0,b0)独立于影像间的共面约束条件,无法在内定标中精确解算,而该参数又与外定标参数完全相关,因此,在进行每个CMOS内定标解算时直接忽略常数项,只计算其他高次项参数;鉴于每片CMOS高阶内定标参数的解算过程相同,仅对一片的参数解算进行叙述,具体如下:
根据构建的平差模型(16),通过模型线性化处理,对每对同名像点可构建误差方程,如式(17)所示:
其中,和分别为左右影像上像点对应的改正数向量,y=[da1,…da9,db1,…db9]T为相机内定标参数改正向量;表示各同名像点物方平面坐标改正向量,物方坐标的高程直接从影像覆盖区域的DSM中插值得到,而不在内定标中解算;矩阵和分别为左右像点误差方程中定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量;以左影像上像点为例,其误差方程中各矩阵具体形式如下:
最后,利用最小二乘平差计算y,如式(18)所示:
其中,k代表该CMOS对应的重叠影像上同名点的数量;
同样地,基于最小二乘的内定标解算也为迭代计算,根据每次迭代解算的结果更新当前的内定标参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算;
63)常数项整体解算
为了确保各片CMOS间的拼接精度,需要根据相邻CMOS重叠影像间的同名像点解算各片CMOS的常数项,这里选择某一片CMOS为基准片,基准片的常数项仍不计算,以基准片为参考,计算所有非基准片的常数项;
基于上述构建的平差方程对同名点分别构建误差方程,方程的具体形式与(17)类似,不同之处在于这里的待解算参数为所有非基准CMOS指向角模型的常数项,具体如下式:
其中,z=[dz1,dz2…dzn]T为所有非基准CMOS的常数项改正数,其中dzi=[da0,db0]i,n为非基准CMOS个数;ti仍为各同名像点物方平面坐标改正向量,矩阵和分别为左右像点误差方程中常数项定标参数对应的偏导数系数矩阵;矩阵和分别为左右像点误差方程中物方坐标对应的偏导数系数矩阵,同样是根据各自的影像建立的平差方程通过线性化得到的,和分别为左右像点误差方程中常数向量;
最后,利用最小二乘平差计算z,如式(20)所示:
其中,λ代表用于计算的同名点数量;
同样地,该解算为迭代计算,根据每次迭代解算的结果更新当前的内定标常数项参数,当连续两次计算的内定标参数改正数小于阈值时,停止迭代解算。
8.星地联合的面阵成像遥感卫星全自主几何定标系统,其特征在于,包括如下模块:
外定标观测值提取模块,用于当卫星运行到太阳阴影区时,利用卫星的敏捷机动能力,对合适的天区进行拍摄,获取恒星星图数据,并在导航星表的辅助下,提取恒星作为外定标观测值;
外定标模型构建模块,用于基于恒星观测向量和卫星载荷视线指向向量的共线关系建立面阵相机对天观测的几何成像模型,通过在模型中引入广义的相机安装角建立外定标模型;
外定标参数解算模块,用于基于序列恒星观测值,采用最小二乘方法迭代解算外定标参数;
内定标观测值提取模块,用于当卫星运行到同一圈的太阳阳照面时,利用卫星的机动能力在短时间内获取待标定相机拍摄的同一地区三次覆盖的影像,两次重叠中行列两个方向重叠度的比值应不同,并从影像的重叠区域匹配密集的同名点,作为内定标解算的观测值;
内定标模型构建模块,用于基于像点、物点和投影中心三点共线关系建立面阵相机对地观测的几何成像模型,通过在模型中引入二元多项式拟合的指向角模型建立内定标模型;
内定标参数计算模块,用于在外定标的基础上,基于同名像点的共面约束条件,在参考DSM的约束下,采用多片CMOS高阶模型参数和常数项分片优化的方法解算内定标参数;
定标结果获取模块,用于在内定标解算参数的基础上,重复外定标模型构建模块和外定标参数解算模块,重新进行外定标解算,进而得到内外定标结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111133176.0A CN113900122B (zh) | 2021-09-27 | 2021-09-27 | 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111133176.0A CN113900122B (zh) | 2021-09-27 | 2021-09-27 | 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113900122A CN113900122A (zh) | 2022-01-07 |
CN113900122B true CN113900122B (zh) | 2024-05-14 |
Family
ID=79029654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111133176.0A Active CN113900122B (zh) | 2021-09-27 | 2021-09-27 | 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113900122B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
CN107831502A (zh) * | 2017-10-17 | 2018-03-23 | 中国科学院电子学研究所 | 卫星指向夹角近似值的计算方法 |
KR20190026452A (ko) * | 2017-09-05 | 2019-03-13 | 순천대학교 산학협력단 | 위성영상과 제공 rpc로부터 제작된 수치표고모델의 자동 오차보정 방법 |
CN109828292A (zh) * | 2019-02-14 | 2019-05-31 | 上海卫星工程研究所 | 基于空间相机的可驱动天线标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10234533B2 (en) * | 2016-09-09 | 2019-03-19 | The Charles Stark Draper Laboratory, Inc. | Position determination by observing a celestial object transit the sun or moon |
-
2021
- 2021-09-27 CN CN202111133176.0A patent/CN113900122B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
KR20190026452A (ko) * | 2017-09-05 | 2019-03-13 | 순천대학교 산학협력단 | 위성영상과 제공 rpc로부터 제작된 수치표고모델의 자동 오차보정 방법 |
CN107831502A (zh) * | 2017-10-17 | 2018-03-23 | 中国科学院电子学研究所 | 卫星指向夹角近似值的计算方法 |
CN109828292A (zh) * | 2019-02-14 | 2019-05-31 | 上海卫星工程研究所 | 基于空间相机的可驱动天线标定方法 |
Non-Patent Citations (4)
Title |
---|
Robust camera distortion calibration via unified PRC model for optical remote sensing satellites;皮英冬等;《IEEE transactions on geoscience and remote sensing》;20220811;全文 * |
内外参数精确建模的太阳敏感器标定;樊巧云;张广军;魏新国;;北京航空航天大学学报;20111013(10);全文 * |
缺少地面控制点的光学卫星遥感影像几何精处理质量控制方法;皮英冬;《中国优秀博士学位论文全文数据库(基础科学辑)》;20220615;全文 * |
静止轨道高分辨率面阵相机几何处理关键技术研究;常学立;《中国优秀博士学位论文全文数据库(基础科学辑)》;20170315;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113900122A (zh) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113900125B (zh) | 星地联合的线阵成像遥感卫星全自主几何定标方法及系统 | |
Hu et al. | Understanding the rational function model: methods and applications | |
Wang et al. | On-orbit geometric calibration and geometric quality assessment for the high-resolution geostationary optical satellite GaoFen4 | |
CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及系统 | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN111508029A (zh) | 星载分片线阵ccd光学相机整体几何定标方法及系统 | |
CN110030978B (zh) | 一种全链路光学卫星几何成像模型构建方法及系统 | |
Burman | Calibration and orientation of airborne image and laser scanner data using GPS and INS | |
CN113514829B (zh) | 面向InSAR的初始DSM的区域网平差方法 | |
CN110006452A (zh) | 高分六号宽视场相机相对几何定标方法及系统 | |
CN111508028A (zh) | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 | |
CN113739767B (zh) | 针对国产面阵摆扫成像系统获取的影像生产正射影像的方法 | |
Yang et al. | Integrated geometric self-calibration of stereo cameras onboard the ZiYuan-3 satellite | |
CN114838740A (zh) | 一种考虑不同经纬度区域的卫星图像几何定标方法 | |
Ren et al. | A global adjustment method for photogrammetric processing of Chang’E-2 stereo images | |
CN117710479A (zh) | 基于空间约束的摆扫载荷在轨内参测量及标校方法 | |
CN115311366A (zh) | 基于rpc模型的星载分片线阵传感器几何定标方法及系统 | |
Pan et al. | Systematic geolocation errors of FengYun-3D MERSI-II | |
CN113900122B (zh) | 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 | |
Wolfe | MODIS geolocation | |
Cao | Geometric Calibration of the Whiskbroom DaQi-1 WSI Based on an Equivalent Linear-array CCD | |
CN111667533A (zh) | 一种基于有理函数模型获取卫星严格成像模型参数的方法 | |
Wang et al. | On-orbit calibration approach based on partial calibration-field coverage for the GF-1/WFV camera | |
Hwangbo et al. | Integration of orbital and ground image networks for the automation of rover localization | |
CN114838739A (zh) | 一种考虑完整回归周期的卫星图像几何定标方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |