CN114838739B - 一种考虑完整回归周期的卫星图像几何定标方法 - Google Patents
一种考虑完整回归周期的卫星图像几何定标方法 Download PDFInfo
- Publication number
- CN114838739B CN114838739B CN202210557021.8A CN202210557021A CN114838739B CN 114838739 B CN114838739 B CN 114838739B CN 202210557021 A CN202210557021 A CN 202210557021A CN 114838739 B CN114838739 B CN 114838739B
- Authority
- CN
- China
- Prior art keywords
- satellite
- calibration
- image
- satellite image
- error
- 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 28
- 239000013598 vector Substances 0.000 claims description 21
- 238000003384 imaging method Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 9
- 239000000523 sample Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种考虑完整回归周期的卫星图像几何定标方法,包括:获取一个完整卫星回归周期内同一定标场的不同地方时的卫星图像集;基于该定标场控制数据,获取所述的卫星图像集中的沿卫星轨道方向多行像素内均匀分布的控制点;选取其中一个基准地方时图像中的控制点,基于第一几何定标模型,得到所述的地方时的相机定标参数;构建不同地方时的姿态误差补偿模型,建立考虑卫星完整回归周期的第二几何定标模型;利用其他地方时图像中的控制点,解算得到姿态补偿参数;根据所述的基准地方时的相机定标参数、不同地方时的姿态补偿参数,基于第二几何定标模型,更新定标文件。本发明可以实现同一区域不同地方时卫星图像的统一定标,提升定标精度。
Description
技术领域
本发明涉及卫星遥感技术领域,尤其涉及一种考虑完整回归周期的卫星图像几何定标方法。
背景技术
开展卫星图像高精度地面定标是高分辨率卫星图像高定位处理的前提条件。目前高分辨率卫星图像定标方法主要是采用常规几何定标模型,该模型主要基于某一定标场的单景卫星图像进行定标,即利用定标场的高精度参考数据和相应的某一时刻定标景卫星图像对相机参数进行在轨标定,然后该参数应用于整个卫星各个地方时图像处理。该定标方法在太阳同步轨道卫星在轨几何定标中效果较好,也是目前多数在轨卫星采用的定位方法,具有较好的普适性。
然而,对于目前逐步投入使用的高分辨率非太阳同步轨道卫星而言,由于其轨道的特殊性,在一个完整的卫星回归周期内,存在每次对同一地区成像都会有不同的地方时、不同的太阳辐照度,导致每次对同一地区成像时卫星器件受热情况差异较大,相机指向角存在随回归周期而发生周期性变化的现象。此时,若基于常规几何定标模型,并利用单景定标图像进行相机参数的标定,只能标定出该定标景图像成像时的相机参数,无法标定出整个回归周期内其他地方时下由于太阳辐射的差异引起的相机参数的变化,进而势必导致该定标参数无法很好的适用于卫星回归周期内不同地方时拍摄的所有图像的定位处理。因此,为保证不同地方时下的相机参数都能得到较好的标定,解算适于卫星常态化高精度定位处理的相机参数,急需开展针对非太阳同步轨道卫星的在轨几何定标方法研究。
由于常规几何定标方法主要采用常规几何定标模型,利用某一定标场区域的单景卫星图像进行在轨几何定标,并未考虑整个卫星整个回归周期内不同地方时的卫星影像情况,无法标定回归周期内同一地区不同地方时成像时由于太阳辐射差异较大而引起的相机参数变化情况,使得该定标参数对非太阳同步轨道卫星中不同地方时影像的适用性相对较差,制约了高分辨率卫星图像的在轨高定位处理,因此亟待解决非太阳同步轨道卫星回归周期内不同地方时的在轨几何定标技术难题。
发明内容
针对上述问题,本发明公开了一种考虑完整回归周期的卫星图像几何定标方法,包括:
S1,获取卫星在一个完整的卫星回归周期内同一定标场的不同地方时的卫星图像集;
S2,将步骤S1获取的卫星图像集中的各个卫星图像都分别与所述的定标场的正射影像和数字高程模型进行地面控制点匹配,获取所述的卫星图像集中的沿卫星轨道方向的多行连续像素内均匀分布的地面控制点的像方坐标和物方坐标;所述的地面控制点位于定标场中。
S3,建立第一几何定标模型,选取卫星图像集中某一地方时的卫星图像的地面控制点,该卫星图像对应的地方时称为基准地方时,并将该地面控制点的像方坐标和物方坐标代入第一几何定标模型,采用先解算外方位元素,后解算内方位元素的定标方法,解算第一几何定标模型,得到所述的基准地方时下的相机的内定标参数(m0,m1,m2,m3,...,mp)和(n0,n1,n2,n3,...,np)以及外定标参数p为相机的定标参数的阶数;
所述的第一几何定标模型,其表达式为:
式中,S为卫星图像中的地面控制点对应的相机的成像探元的探元编号,(XWGS84,YWGS84,ZWGS84)为地面控制点在WGS84坐标系下的物方空间坐标,(XGPS,YGPS,ZGPS)为卫星GPS天线相位中心在WGS84坐标系下的空间坐标,λ为比例因子,为J2000坐标系至WGS84坐标系的旋转矩阵,/>为卫星本体坐标系至J2000坐标系的旋转矩阵,(ψx,ψy)为卫星图像中的地面控制点对应的相机的成像探元在相机坐标系下的指向角。指向角指的是成像探元对应光线在相机坐标系下的方向。
S4,构建与地方时相关联的姿态误差补偿模型,建立考虑卫星完整回归周期的第二几何定标模型;
所述的步骤S4,包括,同一卫星回归周期内不同地方时的成像误差矩阵RT,即为绕X,Y,Z坐标轴旋转角度变化构成的矩阵;构建考虑卫星完整回归周期的在轨几何定标模型,其表达式为:
其中,Δωt为卫星的随地方时而变化的滚动角误差,为卫星的随地方时而变化的俯仰角误差,Δκt为卫星的随地方时而变化的偏航角误差;采用地方时的三次多项式模型进行姿态误差拟合,针对基准地方时,构建与地方时相关联的姿态误差补偿模型,其表达式为:
其中,Δt为卫星图像集中的其他图像对应的地方时相对于基准地方时的变化量,(a0,a1,a2,a3)、(b0,b1,b2,b3)、(c0,c1,c2,c3)分别为滚动角、俯仰角和偏航角的与地方时相关联的姿态误差补偿模型参数。
S5、利用除基准地方时外的其他地方时的卫星图像中的地面控制点的坐标,基于第二几何定标模型,解算与地方时相关联的姿态误差补偿模型参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3)。
所述的步骤S5,包括:
对在轨几何定标模型进行变形,引入两个中间矢量[X,Y,Z]和[x,y,z],其表达式为:
其中,根据像方坐标和物方坐标的转换关系,得到:
进一步得到两个中间矢量的关系式为:
对两个中间矢量的关系式在处进行泰勒级数展开,得到泰勒级数展开式,其表达式为:
对泰勒级数展开式建立误差方程,解算卫星在不同地方时的姿态误差对于步骤S2中获取的所述的卫星图像集中的沿卫星轨道方向的多行像素内均匀分布的地面控制点,基于所述的基准地方时的相机的内定标参数、外定标参数和两个中间矢量,对泰勒级数展开式建立误差方程,其表达式为:
Vi=AiΔX-Li,
通过对上述误差方程求解,得到ΔX,其中,ΔX为卫星在拍摄卫星图像集中的各个卫星图像所处的不同地方时的姿态误差补偿向量,Ai表示根据第i个地面控制点的物方坐标计算得到的误差方程的乘性系数,Li表示根据第i个地面控制点的像方坐标和物方坐标计算得到的误差方程的加性系数,(xi,yi,zi)和(Xi,Yi,Zi)分别表示第i个地面控制点的像方坐标和物方坐标,Vi表示根据第i个地面控制点的像方坐标和物方坐标计算得到的修正量;
根据与地方时相关联的姿态误差补偿模型,建立参数误差方程组:
其中,Vai、Vbi和Vci分别表示不同地方时下卫星的滚动角、俯仰角和偏航角的误差残差向量,为卫星的滚动角,ωi为卫星的俯仰角,κi为卫星的偏航角,Aai、Abi和Aci分别表示卫星的滚动角、俯仰角和偏航角的参数误差方程的参数矩阵,Xa、Xb和Xc分别表示与地方时相关联的姿态误差补偿模型的参数,Lbi、Lai和Lci分别表示卫星的滚动角、俯仰角和偏航角的误差向量,/>Δωi和Δκi分别表示不同地方时下卫星的姿态修正量,/>Δωi、Δκi分别表示第i个地方时下的卫星的姿态修正量,ti表示第i个地方时的时间;利用解算得到的卫星在拍摄不同地方时的卫星图像时的姿态误差,通过最小二乘法求解参数误差方程组,得到卫星图像集中各个卫星图像对应的与地方时相关联的姿态误差补偿模型的参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3)。
S6,根据基准地方时的相机的内定标参数、外定标参数以及与地方时相关联的姿态误差补偿模型参数,基于第二几何定标模型,更新相机的定标参数;利用更新后的定标参数对卫星图像集中的各个卫星图像进行定标处理,并利用地面控制点坐标,分别计算定标前和定标后的所述的卫星图像集的定位精度,当定标后的卫星图像集的定位精度满足预设值要求时,更新定标文件,完成考虑完整卫星回归周期的卫星图像几何定标。
所述的卫星图像集要求尽可能均匀分布于卫星完整回归周期内,确保周期内不同地方时的图像均有分布。其中地方时分别为t1、t2、...、tn(其中0:00:00≤n<24:00:00),对应卫星图像集中的图像分别为I1、I2、...、In。
所述的卫星图像集,是指卫星在一个完整的卫星回归周期内同一定标场的不同地方时进行拍摄所得到的卫星图像。
所述的地面控制点是经过实地测量或者通过其它方法获取到其精确位置坐标信息的地面点,用于对卫星图像的定位误差进行纠正。
本发明的有益效果为:
本发明可以顾及到卫星整个回归周期内不同地方时的系列图像,可以实现同一区域不同地方时卫星图像的统一定标,提升定标精度。
附图说明
图1为本发明方法的实施流程图。
具体实施方式
为了更好的了解本发明内容,这里给出一个实施例。
如图1所示,本发明公开了一种考虑完整回归周期的卫星图像几何定标方法,包括:
S1,获取卫星在至少一个完整的卫星回归周期内同一定标场的不同地方时的卫星图像集;
S2,将步骤S1获取的卫星图像集中的各个卫星图像都分别与所述的定标场的正射影像和数字高程模型进行地面控制点匹配,获取所述的卫星图像集中的沿卫星轨道方向的多行连续像素内均匀分布的地面控制点的像方坐标和物方坐标;所述的地面控制点位于定标场中。
S3,建立第一几何定标模型,选取卫星图像集中某一地方时的卫星图像的地面控制点,该卫星图像对应的地方时称为基准地方时,并将该地面控制点的像方坐标和物方坐标代入第一几何定标模型,采用先解算外方位元素,后解算内方位元素的定标方法,解算第一几何定标模型,得到所述的基准地方时下的相机的内定标参数(m0,m1,m2,m3,...,mp)和(n0,n1,n2,n3,...,np)以及外定标参数p为相机的定标参数的阶数;
所述的第一几何定标模型,其表达式为:
式中,S为卫星图像中的地面控制点对应的相机的成像探元的探元编号,(XWGS84,YWGS84,ZWGS84)为地面控制点在WGS84坐标系下的物方空间坐标,(XGPS,YGPS,ZGPS)为卫星GPS天线相位中心在WGS84坐标系下的空间坐标,λ为比例因子,为J2000坐标系至WGS84坐标系的旋转矩阵,/>为卫星本体坐标系至J2000坐标系的旋转矩阵,(ψx,ψy)为卫星图像中的地面控制点对应的相机的成像探元在相机坐标系下的指向角。指向角指的是成像探元对应光线在相机坐标系下的方向。
S4,构建与地方时相关联的姿态误差补偿模型,建立考虑卫星完整回归周期的第二几何定标模型;
所述的步骤S4,包括:同一卫星回归周期内不同地方时的成像误差矩阵RT,即为绕X,Y,Z坐标轴旋转角度变化构成的矩阵;构建考虑卫星完整回归周期的在轨几何定标模型,其表达式为:
其中,Δωt为卫星的随地方时而变化的滚动角误差,为卫星的随地方时而变化的俯仰角误差,Δκt为卫星的随地方时而变化的偏航角误差;采用地方时的三次多项式模型进行姿态误差拟合,针对基准地方时,构建与地方时相关联的姿态误差补偿模型为:
其中,Δt为卫星图像集中的其他图像对应的地方时相对于基准地方时的变化量,(a0,a1,a2,a3)、(b0,b1,b2,b3)、(c0,c1,c2,c3)分别为滚动角、俯仰角和偏航角的与地方时相关联的姿态误差补偿模型参数。
S5、将除基准地方时外的其他地方时的卫星图像中的地面控制点的坐标,代入第二几何定标模型,解算与地方时相关联的姿态误差补偿模型参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3)。
所述的步骤S5,包括:
对在轨几何定标模型进行变形,引入两个中间矢量[X,Y,Z]和[x,y,z],其表达式为:
其中,根据像方坐标和物方坐标的转换关系,得到:
进一步得到两个中间矢量的关系式为:
对两个中间矢量的关系式在处进行泰勒级数展开,得到泰勒级数展开式,其表达式为:
对泰勒级数展开式建立误差方程,解算卫星在不同地方时的姿态误差对于步骤S2中获取的所述的卫星图像集中的沿卫星轨道方向的多行像素内均匀分布的地面控制点,基于所述的基准地方时的相机的内定标参数、外定标参数和两个中间矢量,对泰勒级数展开式建立误差方程,其表达式为:
Vi=AiΔX-Li Pi
通过对上述误差方程求解,得到ΔX,其中,ΔX为卫星在拍摄卫星图像集中的各个卫星图像所处的不同地方时的姿态误差补偿向量,Ai表示根据第i个地面控制点的物方坐标计算得到的误差方程的乘性系数,Li表示根据第i个地面控制点的像方坐标和物方坐标计算得到的误差方程的加性系数,(xi,yi,zi)和(Xi,Yi,Zi)分别表示第i个地面控制点的像方坐标和物方坐标,Pi表示系数矩阵Ai的权值,Vi表示根据第i个地面控制点的像方坐标和物方坐标计算得到的修正量;
根据与地方时相关联的姿态误差补偿模型,建立参数误差方程组:
其中,Vai、Vbi和Vci分别表示不同地方时下卫星的滚动角、俯仰角和偏航角的误差残差向量,为卫星的滚动角,ωi为卫星的俯仰角,κi为卫星的偏航角,Aai、Abi和Aci分别表示卫星的滚动角、俯仰角和偏航角的参数误差方程的参数矩阵,Xa、Xb和Xc分别表示与地方时相关联的姿态误差补偿模型的参数,Lbi、Lai和Lci分别表示卫星的滚动角、俯仰角和偏航角的误差向量,/>Δωi和Δκi分别表示不同地方时下卫星的姿态修正量,/>Δωi、Δκi分别表示第i个地方时下的卫星的姿态修正量,ti表示第i个地方时的时间;利用解算得到的卫星在拍摄不同地方时的卫星图像时的姿态误差,通过最小二乘法求解参数误差方程组,得到卫星图像集中各个卫星图像对应的与地方时相关联的姿态误差补偿模型的参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3)。
S6,根据基准地方时的相机的内定标参数、外定标参数以及不同地方时的姿态误差补偿模型参数,基于第二几何定标模型,更新相机的定标参数;利用更新后的定标参数对卫星图像集中的各个卫星图像进行定标处理,并利用地面控制点坐标,分别计算定标前和定标后的所述的卫星图像集的定位精度,当定标后的卫星图像集的定位精度满足预设值要求时,更新定标文件,完成考虑完整卫星回归周期的卫星图像几何定标。
所述的卫星图像集要求尽可能均匀分布于卫星完整回归周期内,确保周期内不同地方时的图像均有分布。其中地方时分别为t1、t2、...、tn(其中0:00:00≤n<24:00:00),对应的卫星图像集中的原始图像分别为I1、I2、...、In。
所述的地面控制点是经过实地测量或者通过其它方法获取到其精确位置坐标信息的地面点,用于对卫星图像的定位误差进行纠正。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (4)
1.一种考虑完整回归周期的卫星图像几何定标方法,其特征在于,包括:
S1,获取卫星在至少一个完整的卫星回归周期内同一定标场的不同地方时的卫星图像集;
S2,将步骤S1获取的卫星图像集中的各个卫星图像都分别与所述的定标场的正射影像和数字高程模型进行地面控制点匹配,获取所述的卫星图像集中的沿卫星轨道方向的多行连续像素内均匀分布的地面控制点的像方坐标和物方坐标;
S3,建立第一几何定标模型,选取卫星图像集中某一地方时的卫星图像的地面控制点,该卫星图像对应的地方时称为基准地方时,并将该地面控制点的像方坐标和物方坐标代入第一几何定标模型,采用先解算外方位元素,后解算内方位元素的定标方法,解算第一几何定标模型,得到所述的基准地方时下的相机的内定标参数(m0,m1,m2,m3,...,mp)和(n0,n1,n2,n3,...,np)以及外定标参数和p为相机的定标参数的阶数;
S4,构建与地方时相关联的姿态误差补偿模型,建立考虑卫星完整回归周期的第二几何定标模型;
S5、利用除基准地方时外的其他地方时的卫星图像中的地面控制点的坐标,基于第二几何定标模型,解算与地方时相关联的姿态误差补偿模型参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3);
S6,根据基准地方时的相机的内定标参数、外定标参数以及与地方时相关联的姿态误差补偿模型参数,基于第二几何定标模型,更新相机的定标参数;利用更新后的定标参数对卫星图像集中的各个卫星图像进行定位处理,并利用地面控制点坐标,分别计算定标前和定标后的所述的卫星图像集的定位精度,当定标后的卫星图像集的定位精度满足预设值要求时,更新定标文件,完成考虑完整卫星回归周期的卫星图像几何定标。
2.如权利要求1所述的考虑完整回归周期的卫星图像几何定标方法,其特征在于,所述的第一几何定标模型,其表达式为:
式中,S为卫星图像中的地面控制点对应的相机的成像探元的探元编号,(XWGS84,YWGS84,ZWGS84)为地面控制点在WGS84坐标系下的物方空间坐标,(XGPS,YGPS,ZGPS)为卫星GPS天线相位中心在WGS84坐标系下的空间坐标,λ为比例因子,为J2000坐标系至WGS84坐标系的旋转矩阵,/>为卫星本体坐标系至J2000坐标系的旋转矩阵,(ψx,ψy)为卫星图像中的地面控制点对应的相机的成像探元在相机坐标系下的指向角。
3.如权利要求2所述的考虑完整回归周期的卫星图像几何定标方法,其特征在于,所述的步骤S4,包括:
同一卫星回归周期内不同地方时的成像误差矩阵RT,即为绕X,Y,Z坐标轴旋转角度变化构成的矩阵;构建考虑卫星完整回归周期的在轨几何定标模型,其表达式为:
其中,Δωt为卫星的随地方时而变化的滚动角误差,为卫星的随地方时而变化的俯仰角误差,Δκt为卫星的随地方时而变化的偏航角误差;采用地方时的三次多项式模型进行姿态误差拟合,针对基准地方时,构建与地方时相关联的姿态误差补偿模型,其表达式为:
其中,Δt为卫星图像集中的其他图像对应的地方时相对于基准地方时的变化量,(a0,a1,a2,a3)、(b0,b1,b2,b3)、(c0,c1,c2,c3)分别为滚动角、俯仰角和偏航角的与地方时相关联的姿态误差补偿模型参数。
4.如权利要求3所述的考虑完整回归周期的卫星图像几何定标方法,其特征在于,所述的步骤S5,包括:
对在轨几何定标模型进行变形,引入两个中间矢量[X,Y,Z]和[x,y,z],其表达式为:
其中,根据像方坐标和物方坐标的转换关系,得到:
进一步得到两个中间矢量的关系式为:
对两个中间矢量的关系式在处进行泰勒级数展开,得到泰勒级数展开式,其表达式为:
对泰勒级数展开式建立误差方程,解算卫星在不同地方时的卫星图像的姿态误差对于步骤S2中获取的所述的卫星图像集中的沿卫星轨道方向的多行像素内均匀分布的地面控制点,基于所述的基准地方时的相机的内定标参数、外定标参数和两个中间矢量,对泰勒级数展开式建立误差方程,其表达式为:
Vi=AiΔX-Li,
通过对上述误差方程求解,得到ΔX,其中,ΔX为卫星在拍摄卫星图像集中的各个卫星图像所处的不同地方时的姿态误差补偿向量,Ai表示根据第i个地面控制点的物方坐标计算得到的误差方程的乘性系数,Li表示根据第i个地面控制点的像方坐标和物方坐标计算得到的误差方程的加性系数,(xi,yi,zi)和(Xi,Yi,Zi)分别表示第i个地面控制点的像方坐标和物方坐标,Vi表示根据第i个地面控制点的像方坐标和物方坐标计算得到的修正量;
根据与地方时相关联的姿态误差补偿模型,建立参数误差方程组:
其中,Vai、Vbi和Vci分别表示不同地方时下卫星的滚动角、俯仰角和偏航角的误差残差向量,为卫星的滚动角,ωi为卫星的俯仰角,κi为卫星的偏航角,Aai、Abi和Aci分别表示卫星的滚动角、俯仰角和偏航角的参数误差方程的参数矩阵,Xa、Xb和Xc分别表示与地方时相关联的姿态误差补偿模型的参数,Lbi、Lai和Lci分别表示卫星的滚动角、俯仰角和偏航角的误差向量,/>Δωi和Δκi分别表示不同地方时下卫星的姿态修正量,/>Δωi、Δκi分别表示第i个地方时下的卫星的姿态修正量,ti表示第i个地方时的时间;利用解算得到的卫星在拍摄不同地方时的卫星图像时的姿态误差,通过最小二乘法求解参数误差方程组,得到卫星图像集中各个卫星图像对应的与地方时相关联的的姿态误差补偿模型的参数(a0,a1,a2,a3)、(b0,b1,b2,b3)和(c0,c1,c2,c3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210557021.8A CN114838739B (zh) | 2022-05-20 | 2022-05-20 | 一种考虑完整回归周期的卫星图像几何定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210557021.8A CN114838739B (zh) | 2022-05-20 | 2022-05-20 | 一种考虑完整回归周期的卫星图像几何定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114838739A CN114838739A (zh) | 2022-08-02 |
CN114838739B true CN114838739B (zh) | 2024-05-03 |
Family
ID=82573177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210557021.8A Active CN114838739B (zh) | 2022-05-20 | 2022-05-20 | 一种考虑完整回归周期的卫星图像几何定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114838739B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115563800B (zh) * | 2022-10-24 | 2024-09-03 | 中科南京移动通信与计算创新研究院 | 一种海量卫星仿真建模方法、系统、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007064836A (ja) * | 2005-08-31 | 2007-03-15 | Kyushu Institute Of Technology | カメラキャリブレーションの自動化アルゴリズム |
CN103274056A (zh) * | 2013-04-24 | 2013-09-04 | 上海卫星工程研究所 | 一种利用在轨数据对卫星降交点地方时漂移标定的方法 |
CN103940450A (zh) * | 2014-04-29 | 2014-07-23 | 上海卫星工程研究所 | 利用在轨数据对卫星降交点地方时的漂移进行标定的方法 |
CN111508028A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 |
-
2022
- 2022-05-20 CN CN202210557021.8A patent/CN114838739B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007064836A (ja) * | 2005-08-31 | 2007-03-15 | Kyushu Institute Of Technology | カメラキャリブレーションの自動化アルゴリズム |
CN103274056A (zh) * | 2013-04-24 | 2013-09-04 | 上海卫星工程研究所 | 一种利用在轨数据对卫星降交点地方时漂移标定的方法 |
CN103940450A (zh) * | 2014-04-29 | 2014-07-23 | 上海卫星工程研究所 | 利用在轨数据对卫星降交点地方时的漂移进行标定的方法 |
CN111508028A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 |
Non-Patent Citations (2)
Title |
---|
中高轨卫星海洋成像图像几何定位精度提升方法;刘晓磊;董小萌;潘忠石;张过;宋鹏飞;;航天器工程;20200415(02);全文 * |
高分辨率光学遥感卫星在轨几何定标现状与展望;王密;田原;程宇峰;;武汉大学学报(信息科学版);20171105(11);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114838739A (zh) | 2022-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114838740B (zh) | 一种考虑不同经纬度区域的卫星图像几何定标方法 | |
CN107644435B (zh) | 顾及姿态校正的敏捷光学卫星无场地几何标定方法及系统 | |
CN106403902B (zh) | 一种星地协同的光学卫星在轨实时几何定位方法及系统 | |
CN110388898B (zh) | 构建虚拟控制点约束的多源多重覆盖遥感影像平差方法 | |
CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及系统 | |
CN102901519B (zh) | 一种基于探元指向角光学推扫卫星在轨分步几何定标方法 | |
CN106885585B (zh) | 一种基于光束法平差的星载摄影测量系统一体化检校方法 | |
CN107194974B (zh) | 一种基于多次识别标定板图像的多目相机外参标定精度的提高方法 | |
CN105698764B (zh) | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 | |
CN102168972B (zh) | 一种基于rpc的三线阵立体卫星区域网平差改进及检校方法 | |
CN113902645B (zh) | 一种基于逆rd定位模型的星载sar图像rpc校正参数获取方法 | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN103697864B (zh) | 一种基于大虚拟相机的窄视场双相机影像拼接方法 | |
CN106873004B (zh) | 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法 | |
CN110986998B (zh) | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 | |
CN114838739B (zh) | 一种考虑完整回归周期的卫星图像几何定标方法 | |
CN110793542A (zh) | 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法 | |
CN107967700A (zh) | 大视场宽工作距双目相机的在轨几何校正及精度验证方法 | |
CN106504286B (zh) | 卫星影像定位方法和装置 | |
CN114972078A (zh) | 应用sar影像提升国产光学卫星影像无控几何质量方法及系统 | |
CN109115218A (zh) | 一种无地面控制点的遥感影像定位方法 | |
CN108447100A (zh) | 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法 | |
CN117710479A (zh) | 基于空间约束的摆扫载荷在轨内参测量及标校方法 | |
CN114993347B (zh) | 一种考虑不同太阳高度的卫星图像定位处理方法 | |
CN111275773A (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 |