CN109143295B - 一种数字化几何定标场和gcp相结合的内方位元素定标方法 - Google Patents

一种数字化几何定标场和gcp相结合的内方位元素定标方法 Download PDF

Info

Publication number
CN109143295B
CN109143295B CN201811270813.7A CN201811270813A CN109143295B CN 109143295 B CN109143295 B CN 109143295B CN 201811270813 A CN201811270813 A CN 201811270813A CN 109143295 B CN109143295 B CN 109143295B
Authority
CN
China
Prior art keywords
gcp
calibration
internal
parameter
camera
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
Application number
CN201811270813.7A
Other languages
English (en)
Other versions
CN109143295A (zh
Inventor
龙小祥
王冰冰
李庆鹏
崔林
党安松
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Center for Resource Satellite Data and Applications CRESDA
Original Assignee
China Center for Resource Satellite Data and Applications CRESDA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Center for Resource Satellite Data and Applications CRESDA filed Critical China Center for Resource Satellite Data and Applications CRESDA
Priority to CN201811270813.7A priority Critical patent/CN109143295B/zh
Publication of CN109143295A publication Critical patent/CN109143295A/zh
Application granted granted Critical
Publication of CN109143295B publication Critical patent/CN109143295B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

一种数字化几何定标场和GCP相结合的内方位元素定标方法，包括高精度外业地面控制点布设、采集及人工内业立体环境刺点方案、利用数字化几何定标场数据采用密集匹配的方法进行控制点的自动量测、探元指向角相机内定校模型构建及模型参数求解策略和方法、相机内定标参数优化模型构建及模型参数求解。经实践验证表明该内定标方法可行、有效，结果稳定、可靠、内定标精度高，目前已工程化应用到商遥卫星地面数据处理系统中，运行稳定，取得了良好的工程化应用效果。

Description

一种数字化几何定标场和GCP相结合的内方位元素定标方法
技术领域
本发明涉及一种数字化几何定标场和GCP相结合的内方位元素定标方法，属于遥感影像在轨几何检校与处理技术领域。
背景技术
虽然卫星上的相机在卫星发射前会进行严格的实验室内标定，然而由于发射过程中的震动、材料放气、在轨运行时成像条件的改变以及器件的老化等因素的影响，使得相机内方位元素定标参数发生改变，因此需要对卫星进行在轨内方位元素定标，卫星在轨内方位元素定标是光学遥感卫星实现高精度几何定位的关键环节，直接影响卫星影像的内部几何精度和波段配准精度。
由于高分辨卫星地面像元分辨率(GSD)越来越高(注：高景卫星的GSD为0.5m、WorldViewⅢ更是达到了0.3m)、以及内方位元素定标精度要求越来越高(注：通常要求CCD探元的指向角标定精度优于0.3像元)，从而对数字化几何定标场数据(包括数字正射影像数据(DOM数据)和数字高程影像数据(DEM数据))的内部几何精度要求极高(注：DOM数据内部几何精度优于0.1m、DEM数据高程精度优于0.5m)，然而，目前已有的数字化几何检校场数据的内部几何精度为0.3m左右，采用传统的内方位元素定标方法根本无法满足要求，如果采用更好精度的数字化几何定标场数据将耗费巨大的经费，因此，必须针对GSD优于0.5米的高分辨率卫星的高精度内方位元素定标的特殊要求，采用适当的方法进行高精度内方位元素定标。
发明内容
本发明要解决的技术问题是：克服现有技术的不足，提供了一种数字化几何定标场和GCP相结合的内方位元素定标方法，包括高精度外业地面控制点布设、采集及人工内业立体环境刺点方案、利用数字化几何定标场数据采用密集匹配的方法进行控制点的自动量测、探元指向角相机内定校模型构建及模型参数求解策略和方法、相机内定标参数优化模型构建及模型参数求解。经实践验证表明该内定标方法可行、有效，结果稳定、可靠、内定标精度高，目前已工程化应用到商遥卫星地面数据处理系统中，运行稳定，取得了良好的工程化应用效果。
本发明目的通过以下技术方案予以实现：
一种数字化几何定标场和GCP相结合的内方位元素定标方法，包括如下步骤：
步骤一、在待定标影像的某一范围内进行外业地面控制点(GCP)采集并制成点之计，根据所述点之计在立体环境下生成像控点文件；
步骤二、将数字化几何定标场数据与步骤一中的待定标影像进行控制点自动量测，获取控制点坐标；
步骤三，利用卫星的姿态、轨道、行时和卫星的实验室定标参数，建立基于一元三次曲线的探元指向角在轨几何检校模型；
步骤四、利用所述步骤二中控制点坐标和所述步骤三中的探元指向角在轨几何检校模型，迭代求解相机内标定参数；
步骤五，利用所述步骤一中生成的像控点文件，对所述步骤四中相机内定标参数进行优化。
上述数字化几何定标场和GCP相结合的内方位元素定标方法，所述步骤二中控制点的数量不少于所述步骤一中待定标影像的某一范围的像素点的万分之五。
上述数字化几何定标场和GCP相结合的内方位元素定标方法，所述步骤二中，采用灰度匹配和线特征匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点初步自动量测，然后采用最小二乘匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点精确自动量测。
上述数字化几何定标场和GCP相结合的内方位元素定标方法，所述步骤三中基于一元三次曲线的探元指向角在轨几何检校模型为：
Figure BDA0001845935860000031
其中
Figure BDA0001845935860000032
式中，(Xg,Yg,Zg)与(Xgps,Ygps,Zgps)分别表示像点对应的物方点及GPS天线相位中心在WGS84坐标系下的坐标；λ为比例系数，
Figure BDA0001845935860000033
分别代表WGS84坐标系到J2000坐标系的旋转矩阵、J2000坐标系到卫星本体坐标系的旋转矩阵、卫星本体坐标系到相机坐标系的旋转矩阵；(BX,BY,BZ)body代表从传感器投影中心到GPS天线相位中心的偏心矢量在卫星本体坐标系下的坐标；(t)表示当前参数是一个随时间变化的量；(ψx(s),ψy(s))代表探元s在相机坐标系下的指向角，s代表探元列号；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数。
上述数字化几何定标场和GCP相结合的内方位元素定标方法，所述步骤四中利用步骤二中控制点坐标和步骤三中的探元指向角在轨几何检校模型，迭代求解相机内标定参数的具体方法包括如下步骤：
步骤(4a)、计算沿轨指向角残差Fx(XEC,XIC)与垂轨指向角残差Fy(XEC,XIC)：
Figure BDA0001845935860000034
式中，XEC为外定标参数；XIC为内定标参数；a1、b1、c1、a2、b2、c2、a3、b3、c3分别代表相机安装矩阵的9个元素；矢量
Figure BDA0001845935860000035
为物方矢量U，代表从相机投影中心到物方点的矢量在本体坐标系下的坐标；
步骤(4b)、建立每个外定标参数的误差方程：
ViEC＝AiX-Li，权为Pi，
其中
Figure BDA0001845935860000041
Figure BDA0001845935860000042
Figure BDA0001845935860000043
式中，ViEC为外定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000044
的误差向量，
Figure BDA0001845935860000045
为外定标参数的当前值；
Figure BDA0001845935860000046
为内定标参数的当前值；i为序号；Ai是外定标参数误差方程式的系数矩阵；X代表外定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；λ为比例系数，Δpitch(t)、Δroll(t)、Δyaw(t)分别为随时间自适应的俯仰、翻滚以及偏航方向夹角变化量；
步骤(4c)、计算外定标参数改正数X：
X＝(ATPA)-1(ATPL)
其中
Figure BDA0001845935860000047
Figure BDA0001845935860000048
Figure BDA0001845935860000049
式中，A为第一中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4d)、计算外定标参数XEC：
Figure BDA0001845935860000051
将计算所得的外定标参数XEC的值作为外定标参数的当前值
Figure BDA0001845935860000052
当外定标参数改正数X≤10-3时，转入步骤(4e)，否则转入步骤(4a)；
步骤(4e)、建立每个内定标参数的误差方程：
ViIC＝BiY-Li，权为Pi
其中
Figure BDA0001845935860000053
Y＝dXIC＝[dA0 dA1 dA2 dA3 dB0 dB1 dB2 dB3]T
Figure BDA0001845935860000054
式中，ViIC为内定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000055
的误差向量；Bi是内定标参数误差方程式的系数矩阵；Y代表内定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(4f)、计算外定标参数改正数Y：
Y＝(BTPB)-1(BTPL)
其中
Figure BDA0001845935860000056
Figure BDA0001845935860000057
Figure BDA0001845935860000061
式中，B为第四中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4g)、计算内定标参数XIC：
Figure BDA0001845935860000062
将计算所得的内定标参数XIC的值作为内定标参数的当前值
Figure BDA0001845935860000063
当内定标参数改正数Y≤10-3时，转入步骤(4h)，否则转入步骤(4a)；
步骤(4h)、利用步骤(4d)中计算的外定标参数XEC和步骤(4g)中计算内定标参数XIC，更新相机在轨几何检校参数文件。
上述数字化几何定标场和GCP相结合的内方位元素定标方法，所述利用步骤一中生成的像控点文件，对步骤四中相机内定标参数进行优化的方法包括如下步骤：
步骤(5a)、建立每个内定标参数优化的误差方程：
ViGCP＝BiGCPYGCP-LiGCP权为WiGCP
其中
Figure BDA0001845935860000064
YGCP＝dΔXIC＝[dΔA0 dΔA1 dΔA2 dΔA3 dΔB0 dΔB1 dΔB2 dΔB3]T
Figure BDA0001845935860000065
式中，ViGCP为内定标参数优化的误差，LiGCP是基于步骤四中内定标参数XIC和外定标参数XEC当前值
Figure BDA0001845935860000066
的误差向量；BiGCP是误差方程式的系数矩阵；YGCP代表内定标参数优化改正数dΔXIC，d代表改正数符号；WiGCP是当前GCP的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(5b)、计算内定标参数优化改正数YGCP：
YGCP＝(BGCP TWGCPBGCP)-1(BGCP TWGCPLGCP)
其中
Figure BDA0001845935860000071
Figure BDA0001845935860000072
Figure BDA0001845935860000073
式中，BGCP为第五中间矩阵，WGCP为第六中间矩阵，LGCP为第七中间矩阵；
步骤(5c)、根据步骤四中计算的相机内定标参数XIC的当前值
Figure BDA0001845935860000074
和步骤(5b)计算的内定标参数优化改正数YGCP，确定优化后的相机内定标参数XIC：
Figure BDA0001845935860000075
一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述方法的步骤。
本发明相比于现有技术具有如下有益效果：
一种数字化几何定标场和GCP相结合的内方位元素定标方法，包括高精度外业地面控制点布设、采集及人工内业立体环境刺点方案、利用数字化几何定标场数据采用密集匹配的方法进行控制点的自动量测、探元指向角相机内定校模型构建及模型参数求解策略和方法、相机内定标参数优化模型构建及模型参数求解；该方法首先，采用一元三次曲线探元指向角模型作为相机内定校模型，以及利用数字化几何定标场数据采用密集匹配的方法进行控制点的自动量测成果，然后，采用分步、分级迭代解求策略,解算相机内定标参数，最后，利用在待定标影像上某1000～2000行范围内相对均匀分布的外业地面控制点(GCP)采集和人工内业立体环境刺点成果，采用最小二乘平差策略对相机内定标参数进行优化，从而进一步提高相机内定标精度。经中国首颗0.5米分辨率商业遥感卫星及星座——高景01/02/03/04星在轨内定标的实践表明该内定标方法可行、有效，结果稳定、可靠、内定标精度高，目前已工程化应用到商遥卫星地面数据处理系统中，运行稳定，取得了良好的工程化应用效果。
附图说明
图1为本发明的步骤流程图；
图2为本发明实施例的影响控制点的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明的实施方式作进一步详细描述。
一种数字化几何定标场和GCP相结合的内方位元素定标方法，如图1所示，包括如下步骤：
步骤101、在待定标影像的某一范围内进行外业地面控制点(GCP)采集并制成点之计，根据所述点之计在立体环境下生成像控点文件。
步骤102、将数字化几何定标场数据与步骤101中的待定标影像进行控制点自动量测，获取控制点坐标；控制点的数量不少于所述步骤一中待定标影像的某一范围的像素点的万分之五。采用灰度匹配和线特征匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点初步自动量测，然后采用最小二乘匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点精确自动量测。
步骤103，利用卫星的姿态、轨道、行时和卫星的实验室定标参数，建立基于一元三次曲线的探元指向角在轨几何检校模型。
基于一元三次曲线的探元指向角在轨几何检校模型为：
Figure BDA0001845935860000091
其中
Figure BDA0001845935860000092
式中，(Xg,Yg,Zg)与(Xgps,Ygps,Zgps)分别表示像点对应的物方点及GPS天线相位中心在WGS84坐标系下的坐标；λ为比例系数，
Figure BDA0001845935860000093
分别代表WGS84坐标系到J2000坐标系的旋转矩阵、J2000坐标系到卫星本体坐标系的旋转矩阵、卫星本体坐标系到相机坐标系的旋转矩阵；(BX,BY,BZ)body代表从传感器投影中心到GPS天线相位中心的偏心矢量在卫星本体坐标系下的坐标；(t)表示当前参数是一个随时间变化的量；(ψx(s),ψy(s))代表探元s在相机坐标系下的指向角，s代表探元列号；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数。
步骤104、利用所述步骤102中控制点坐标和所述步骤103中的探元指向角在轨几何检校模型，迭代求解相机内标定参数。其具体步骤如下：
步骤(4a)、计算沿轨指向角残差Fx(XEC,XIC)与垂轨指向角残差Fy(XEC,XIC)：
Figure BDA0001845935860000094
式中，XEC为外定标参数；XIC为内定标参数；a1、b1、c1、a2、b2、c2、a3、b3、c3分别代表相机安装矩阵的9个元素；矢量
Figure BDA0001845935860000095
为物方矢量U，代表从相机投影中心到物方点的矢量在本体坐标系下的坐标；
步骤(4b)、建立每个外定标参数的误差方程：
ViEC＝AiX-Li，权为Pi，
其中
Figure BDA0001845935860000101
Figure BDA0001845935860000102
Figure BDA0001845935860000103
式中，ViEC为外定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000104
的误差向量，
Figure BDA0001845935860000105
为外定标参数的当前值；
Figure BDA0001845935860000106
为内定标参数的当前值；i为序号；Ai是外定标参数误差方程式的系数矩阵；X代表外定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；λ为比例系数，Δpitch(t)、Δroll(t)、Δyaw(t)分别为随时间自适应的俯仰、翻滚以及偏航方向夹角变化量；
步骤(4c)、计算外定标参数改正数X：
X＝(ATPA)-1(ATPL)
其中
Figure BDA0001845935860000107
Figure BDA0001845935860000108
Figure BDA0001845935860000109
式中，A为第一中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4d)、计算外定标参数XEC：
Figure BDA0001845935860000111
将计算所得的外定标参数XEC的值作为外定标参数的当前值
Figure BDA0001845935860000112
当外定标参数改正数X≤10-3时，转入步骤(4e)，否则转入步骤(4a)；
步骤(4e)、建立每个内定标参数的误差方程：
ViIC＝BiY-Li，权为Pi
其中
Figure BDA0001845935860000113
Y＝dXIC＝[dA0 dA1 dA2 dA3 dB0 dB1 dB2 dB3]T
Figure BDA0001845935860000114
式中，ViIC为内定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000115
的误差向量；Bi是内定标参数误差方程式的系数矩阵；Y代表内定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(4f)、计算外定标参数改正数Y：
Y＝(BTPB)-1(BTPL)
其中
Figure BDA0001845935860000116
Figure BDA0001845935860000117
Figure BDA0001845935860000121
式中，B为第四中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4g)、计算内定标参数XIC：
Figure BDA0001845935860000122
将计算所得的内定标参数XIC的值作为内定标参数的当前值
Figure BDA0001845935860000123
当内定标参数改正数Y≤10-3时，转入步骤(4h)，否则转入步骤(4a)；
步骤(4h)、利用步骤(4d)中计算的外定标参数XEC和步骤(4g)中计算内定标参数XIC，更新相机在轨几何检校参数文件。
步骤105，利用所述步骤101中生成的像控点文件，对所述步骤104中相机内定标参数进行优化。其具体步骤如下：
步骤(5a)、建立每个内定标参数优化的误差方程：
ViGCP＝BiGCPYGCP-LiGCP权为WiGCP
其中
Figure BDA0001845935860000124
YGCP＝dΔXIC＝[dΔA0 dΔA1 dΔA2 dΔA3 dΔB0 dΔB1 dΔB2 dΔB3]T
Figure BDA0001845935860000125
式中，ViGCP为内定标参数优化的误差，LiGCP是基于步骤四中内定标参数XIC和外定标参数XEC当前值
Figure BDA0001845935860000126
的误差向量；BiGCP是误差方程式的系数矩阵；YGCP代表内定标参数优化改正数dΔXIC，d代表改正数符号；WiGCP是当前GCP的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(5b)、计算内定标参数优化改正数YGCP：
YGCP＝(BGCP TWGCPBGCP)-1(BGCP TWGCPLGCP)
其中
Figure BDA0001845935860000131
Figure BDA0001845935860000132
Figure BDA0001845935860000133
式中，BGCP为第五中间矩阵，WGCP为第六中间矩阵，LGCP为第七中间矩阵；
步骤(5c)、根据步骤四中计算的相机内定标参数XIC的当前值
Figure BDA0001845935860000134
和步骤(5b)计算的内定标参数优化改正数YGCP，确定优化后的相机内定标参数XIC：
Figure BDA0001845935860000135
实施例：
一种数字化几何定标场和GCP相结合的内方位元素定标方法，包括如下步骤：
步骤1、在待定标影像上某1000～2000行范围内，进行外业控制点(GCP)采集并制成点之计，人工内业立体环境刺点，从而获取一定数量、相对均匀分布的像控点文件，通常位于影像中部左右、控制点数量约20个左右、影像两端需要有控制点，其分布如图2所示。
步骤2、利用数字化几何定标场数据(包括DOM数据和DSM数据)，采用密集匹配的方法进行控制点的自动量测，从而获取大量、密集的控制点坐标：首先，在最上层金字塔进行特征点提取(如：SUSAN算子、Harris算子、Forstner算子等)，其次，采用灰度相似性测度使匹配精度达到整像素，然后，利用匹配结果进行几何模型构建方法，消除视差，直到所有金字塔均匹配完成，最后，采用最小二乘匹配，使匹配精度达到亚像素级。
步骤3、利用卫星下传的姿态、轨道、行时等辅助数据以及实验室定标参数，构建基于一元三次曲线的探元指向角在轨几何检校模型：
Figure BDA0001845935860000141
Figure BDA0001845935860000142
式中，(Xg,Yg,Zg)与(Xgps,Ygps,Zgps)分别表示像点对应的物方点及GPS天线相位中心在WGS84坐标系下的坐标；λ为比例系数，
Figure BDA0001845935860000143
分别代表WGS84坐标系到J2000坐标系的旋转矩阵、J2000坐标系到卫星本体坐标系的旋转矩阵、卫星本体坐标系到相机坐标系的旋转矩阵；(BX,BY,BZ)body代表从传感器投影中心到GPS天线相位中心的偏心矢量在卫星本体坐标系下的坐标；(t)表示当前参数是一个随时间变化的量；(ψx(s),ψy(s))代表探元s在相机坐标系下的指向角，s代表探元列号；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数。
在以上几何定标模型中，待定标参数分为外定标参数XEC和内定标参数XIC，外定标参数
Figure BDA0001845935860000144
步骤4、利用步骤2自动匹配的密集控制点坐标和步骤3的探元指向角在轨几何检校模型，采用分步、分级迭代解求策略，求解内定标参数；实现方式为，首先基于最小二乘平差解算外定标参数，恢复相机坐标系在空间中的姿态；然后在此基础上，基于最小二乘平差解算相机内定标参数，确定相机CCD各探元在相机坐标系下的指向角。
步骤4.1、利用步骤2自动匹配的密集控制点信息，设在待定标影像上自动量测了N个均匀分布的地面控制点，各地面控制点对应的物方点和像方点分别记为GCPGi和GCPMi，物方点GCPGi的WGS84地心直角坐标为(Xi,Yi,Zi)，像方点GCPMi的影像坐标为(si,li)；
步骤4.2、令式(1)中：
Figure BDA0001845935860000151
式(1)转化为式(2)：
Figure BDA0001845935860000152
上式中，矢量
Figure BDA0001845935860000153
为物方矢量U，代表从相机投影中心到物方点的矢量在本体坐标系下的坐标；a1,b1,c1a2,b2,c2a3,b3,c3分别代表相机安装矩阵的9个元素；Fx(XEC,XIC)、Fy(XEC,XIC)分别为沿轨指向角残差与垂轨指向角残差；
步骤4.3、对外定标参数XEC和内定标参数XIC赋初值
Figure BDA0001845935860000154
步骤4.4、将内定标参数XIC的当前值视为真值，将外定标参数XEC视为待求的未知参数，将内定标参数XIC和外定标参数XEC的当前值
Figure BDA0001845935860000155
代入式(2)，对每个定向点，对式(2)进行线性化处理，建立误差方程(3)：
ViEC＝AiX-Li，权为Pi (3)
其中
Figure BDA0001845935860000161
ViEC为外定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000162
的误差向量，
Figure BDA0001845935860000163
为外定标参数的当前值；
Figure BDA0001845935860000164
为内定标参数的当前值；i为序号；Ai是外定标参数误差方程式的系数矩阵；X代表外定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；λ为比例系数，Δpitch(t)、Δroll(t)、Δyaw(t)分别为随时间自适应的俯仰、翻滚以及偏航方向夹角变化量；
按式(4)计算法方程系数矩阵，
Figure BDA0001845935860000165
Figure BDA0001845935860000166
上式中，矩阵
Figure BDA0001845935860000167
矩阵
Figure BDA0001845935860000168
矩阵
Figure BDA0001845935860000169
利用最小二乘平差计算X，如式(5)，
X＝(ATPA)-1(ATPL) (5)
利用式(6)更新外定标参数XE的当前值，当外定标参数改正数X≤10-3时，转入步骤4.5，否则转入步骤4.1；
Figure BDA00018459358600001610
需要说明的是，除了采用外定标参数改正数X≤10-3作为判定依据，也可以采用外定标参数的误差ViEC作为判定依据。即当外定标参数的误差ViEC不大于预设值时，转入步骤4.5，否则转入步骤4.1。
步骤4.5、解算内定标参数，将步骤4.4所得外定标参数XEC的当前值视为真值，而内定标参数XI则视为待求的未知参数，将内定标参数XIC和外定标参数XEC的当前值
Figure BDA0001845935860000171
代入公式(2)，对每个定向点，对式(2)进行线性化处理，建立误差方程(7)，
ViIC＝BiY-Li权为Pi (7)
其中，
Figure BDA0001845935860000172
Y＝dXIC＝[dA0 dA1 dA2 dA3 dB0 dB1 dB2 dB3]T
Figure BDA0001845935860000173
式中，ViIC为内定标参数的误差；Li是基于内外定标参数当前值
Figure BDA0001845935860000174
的误差向量；Bi是内定标参数误差方程式的系数矩阵；Y代表内定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
按式(8)计算法方程系数矩阵；
Figure BDA0001845935860000175
Figure BDA0001845935860000176
上式中
Figure BDA0001845935860000177
利用最小二乘平差计算Y，如式(9)；
Y＝(BTPB)-1(BTPL) (9)
利用式(10)更新内定标参数XI的当前值，当内定标参数改正数Y≤10-3时，转入步骤4.6，否则转入步骤4.1；
Figure BDA0001845935860000181
需要说明的是，除了采用内定标参数改正数Y≤10-3作为判定依据，也可以采用内定标参数的误差ViIC作为判定依据。即当外定标参数的误差ViEC不大于预设值时，转入步骤4.6，否则转入步骤4.1。
步骤4.6、根据步骤4.4和步骤4.5所得外定标参数的当前值XEC和内定标参数XIC的当前值，更新相机在轨几何检校参数文件。
步骤5、利用步骤1获取的高精度GCP数据，对步骤4所得的相机内定标参数进行优化。
将步骤4所得的相机内定标参数XIC和外定标参数XEC的当前值(XEC,XIC)为真值，对于步骤1获取的每个高精度GCP，将式(2)进行线性化处理，建立误差方程(11)，
ViGCP＝BiGCPYGCP-LiGCP权为WiGCP (11)
其中
Figure BDA0001845935860000182
YGCP＝dΔXIC＝[dΔA0 dΔA1 dΔA2 dΔA3 dΔB0 dΔB1 dΔB2 dΔB3]T
Figure BDA0001845935860000183
式中，ViGCP为内定标参数优化的误差，LiGCP是基于步骤4中内定标参数XIC和外定标参数XEC当前值
Figure BDA0001845935860000184
的误差向量；BiGCP是误差方程式的系数矩阵；YGCP代表内定标参数优化改正数dΔXIC，d代表改正数符号；WiGCP是当前GCP的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
按式(12)计算法方程系数矩阵；
Figure BDA0001845935860000191
Figure BDA0001845935860000192
上式中
Figure BDA0001845935860000193
利用最小二乘平差计算YGCP，如式(13)；
YGCP＝(BGCP TWGCPBGCP)-1(BGCP TWGCPLGCP) (13)
根据步骤4计算的相机内定标参数XIC和步骤5所得的相机内定标参数优化参数，更新相机在轨几何检校参数文件，如式(14)：
Figure BDA0001845935860000194
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

Claims (7)

1.一种数字化几何定标场和GCP相结合的内方位元素定标方法，其特征在于：包括如下步骤：
步骤一、在待定标影像的某一范围内进行外业地面控制点(GCP)采集并制成点之计，根据所述点之计在立体环境下生成像控点文件；
步骤二、将数字化几何定标场数据与步骤一中的待定标影像进行控制点自动量测，获取控制点坐标；
步骤三，利用卫星的姿态、轨道、行时和卫星的实验室定标参数，建立基于一元三次曲线的探元指向角在轨几何检校模型；
步骤四、利用所述步骤二中控制点坐标和所述步骤三中的探元指向角在轨几何检校模型，迭代求解相机内定标参数；
步骤五，利用所述步骤一中生成的像控点文件，对所述步骤四中相机内定标参数进行优化；
所述步骤三中基于一元三次曲线的探元指向角在轨几何检校模型为：
Figure FDA0002720274020000011
其中
Figure FDA0002720274020000012
式中，(Xg,Yg,Zg)与(Xgps,Ygps,Zgps)分别表示像点对应的物方点及GPS天线相位中心在WGS84坐标系下的坐标；λ为比例系数，
Figure FDA0002720274020000013
分别代表WGS84坐标系到J2000坐标系的旋转矩阵、J2000坐标系到卫星本体坐标系的旋转矩阵、卫星本体坐标系到相机坐标系的旋转矩阵；(BX,BY,BZ)body代表从传感器投影中心到GPS天线相位中心的偏心矢量在卫星本体坐标系下的坐标；(t)表示当前参数是一个随时间变化的量；(ψx(s),ψy(s))代表探元s在相机坐标系下的指向角，s代表探元列号；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
所述利用步骤一中生成的像控点文件，对步骤四中相机内定标参数进行优化的方法包括如下步骤：
步骤(5a)、建立每个内定标参数优化的误差方程：
ViGCP＝BiGCPYGCP-LiGCP 权为WiGCP
其中
Figure FDA0002720274020000021
YGCP＝dΔXIC＝[dΔA0 dΔA1 dΔA2 dΔA3 dΔB0 dΔB1 dΔB2 dΔB3]T
Figure FDA0002720274020000022
式中，ViGCP为内定标参数优化的误差，LiGCP是基于步骤四中内定标参数XIC和外定标参数XEC当前值
Figure FDA0002720274020000023
的误差向量；BiGCP是误差方程式的系数矩阵；YGCP代表内定标参数优化改正数dΔXIC，d代表改正数符号；WiGCP是当前GCP的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(5b)、计算内定标参数优化改正数YGCP：
YGCP＝(BGCP TWGCPBGCP)-1(BGCP TWGCPLGCP)
其中
Figure FDA0002720274020000024
Figure FDA0002720274020000025
Figure FDA0002720274020000031
式中，BGCP为第五中间矩阵，WGCP为第六中间矩阵，LGCP为第七中间矩阵；
步骤(5c)、根据步骤四中计算的相机内定标参数XIC的当前值
Figure FDA0002720274020000032
和步骤(5b)计算的内定标参数优化改正数YGCP，确定优化后的相机内定标参数XIC：
Figure FDA0002720274020000033
2.根据权利要求1所述的一种数字化几何定标场和GCP相结合的内方位元素定标方法，其特征在于：所述步骤二中控制点的数量不少于所述步骤一中待定标影像的某一范围的像素点的万分之五。
3.根据权利要求1所述的一种数字化几何定标场和GCP相结合的内方位元素定标方法，其特征在于：所述步骤二中，采用灰度匹配和线特征匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点初步自动量测，然后采用最小二乘匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点精确自动量测。
4.根据权利要求1所述的一种数字化几何定标场和GCP相结合的内方位元素定标方法，其特征在于：所述步骤四中利用步骤二中控制点坐标和步骤三中的探元指向角在轨几何检校模型，迭代求解相机内定标参数的具体方法包括如下步骤：
步骤(4a)、计算沿轨指向角残差Fx(XEC,XIC)与垂轨指向角残差Fy(XEC,XIC)：
Figure FDA0002720274020000034
式中，XEC为外定标参数；XIC为内定标参数；a1、b1、c1、a2、b2、c2、a3、b3、c3分别代表相机安装矩阵的9个元素；矢量
Figure FDA0002720274020000041
为物方矢量U，代表从相机投影中心到物方点的矢量在本体坐标系下的坐标；
步骤(4b)、建立每个外定标参数的误差方程：
ViEC＝AiX-Li，权为Pi，
其中
Figure FDA0002720274020000042
Figure FDA0002720274020000043
Figure FDA0002720274020000044
式中，ViEC为外定标参数的误差；Li是基于内外定标参数当前值
Figure FDA0002720274020000045
的误差向量，
Figure FDA0002720274020000046
为外定标参数的当前值；
Figure FDA0002720274020000047
为内定标参数的当前值；i为序号；Ai是外定标参数误差方程式的系数矩阵；X代表外定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；λ为比例系数，Δpitch(t)、Δroll(t)、Δyaw(t)分别为随时间自适应的俯仰、翻滚以及偏航方向夹角变化量；
步骤(4c)、计算外定标参数改正数X：
X＝(ATPA)-1(ATPL)
其中
Figure FDA0002720274020000048
Figure FDA0002720274020000049
Figure FDA0002720274020000051
式中，A为第一中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4d)、计算外定标参数XEC：
Figure FDA0002720274020000052
将计算所得的外定标参数XEC的值作为外定标参数的当前值
Figure FDA0002720274020000053
当外定标参数改正数X≤10-3时，转入步骤(4e)，否则转入步骤(4a)；
步骤(4e)、建立每个内定标参数的误差方程：
ViIC＝BiY-Li，权为Pi
其中
Figure FDA0002720274020000054
Y＝dXIC＝[dA0 dA1 dA2 dA3 dB0 dB1 dB2 dB3]T
Figure FDA0002720274020000055
式中，ViIC为内定标参数的误差；Li是基于内外定标参数当前值
Figure FDA0002720274020000056
的误差向量；Bi是内定标参数误差方程式的系数矩阵；Y代表内定标参数改正数，d代表改正数符号；Pi是当前定向点RPi的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(4f)、计算外定标参数改正数Y：
Y＝(BTPB)-1(BTPL)
其中
Figure FDA0002720274020000061
Figure FDA0002720274020000062
Figure FDA0002720274020000063
式中，B为第四中间矩阵，P为第二中间矩阵，L为第三中间矩阵，K为控制点的个数；
步骤(4g)、计算内定标参数XIC：
Figure FDA0002720274020000064
将计算所得的内定标参数XIC的值作为内定标参数的当前值
Figure FDA0002720274020000065
当内定标参数改正数Y≤10-3时，转入步骤(4h)，否则转入步骤(4a)；
步骤(4h)、利用步骤(4d)中计算的外定标参数XEC和步骤(4g)中计算内定标参数XIC，更新相机在轨几何检校参数文件。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如下步骤：
步骤一、在待定标影像的某一范围内进行外业地面控制点(GCP)采集并制成点之计，根据所述点之计在立体环境下生成像控点文件；
步骤二、将数字化几何定标场数据与步骤一中的待定标影像进行控制点自动量测，获取控制点坐标；
步骤三，利用卫星的姿态、轨道、行时和卫星的实验室定标参数，建立基于一元三次曲线的探元指向角在轨几何检校模型；
步骤四、利用所述步骤二中控制点坐标和所述步骤三中的探元指向角在轨几何检校模型，迭代求解相机内定标参数；
步骤五，利用所述步骤一中生成的像控点文件，对所述步骤四中相机内定标参数进行优化；
所述步骤三中基于一元三次曲线的探元指向角在轨几何检校模型为：
Figure FDA0002720274020000071
其中
Figure FDA0002720274020000072
式中，(Xg,Yg,Zg)与(Xgps,Ygps,Zgps)分别表示像点对应的物方点及GPS天线相位中心在WGS84坐标系下的坐标；λ为比例系数，
Figure FDA0002720274020000073
分别代表WGS84坐标系到J2000坐标系的旋转矩阵、J2000坐标系到卫星本体坐标系的旋转矩阵、卫星本体坐标系到相机坐标系的旋转矩阵；(BX,BY,BZ)body代表从传感器投影中心到GPS天线相位中心的偏心矢量在卫星本体坐标系下的坐标；(t)表示当前参数是一个随时间变化的量；(ψx(s),ψy(s))代表探元s在相机坐标系下的指向角，s代表探元列号；pitch、roll、yaw分别为俯仰、翻滚以及偏航方向夹角；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
所述利用步骤一中生成的像控点文件，对步骤四中相机内定标参数进行优化的方法包括如下步骤：
步骤(5a)、建立每个内定标参数优化的误差方程：
ViGCP＝BiGCPYGCP-LiGCP 权为WiGCP
其中
Figure FDA0002720274020000074
YGCP＝dΔXIC＝[dΔA0 dΔA1 dΔA2 dΔA3 dΔB0 dΔB1 dΔB2 dΔB3]T
Figure FDA0002720274020000075
式中，ViGCP为内定标参数优化的误差，LiGCP是基于步骤四中内定标参数XIC和外定标参数XEC当前值
Figure FDA0002720274020000081
的误差向量；BiGCP是误差方程式的系数矩阵；YGCP代表内定标参数优化改正数dΔXIC，d代表改正数符号；WiGCP是当前GCP的像点量测精度对应的权；Fxi和Fyi分别为沿轨指向角残差Fx(XEC,XIC)、垂轨指向角残差Fy(XEC,XIC)的函数模型；A0、A1、A2、A3、B0、B1、B2、B3均一元三次曲线探元指向角内定标模型的系数；
步骤(5b)、计算内定标参数优化改正数YGCP：
YGCP＝(BGCP TWGCPBGCP)-1(BGCP TWGCPLGCP)
其中
Figure FDA0002720274020000082
Figure FDA0002720274020000083
Figure FDA0002720274020000084
式中，BGCP为第五中间矩阵，WGCP为第六中间矩阵，LGCP为第七中间矩阵；
步骤(5c)、根据步骤四中计算的相机内定标参数XIC的当前值
Figure FDA0002720274020000085
和步骤(5b)计算的内定标参数优化改正数YGCP，确定优化后的相机内定标参数XIC：
Figure FDA0002720274020000086
6.根据权利要求5所述的一种计算机可读存储介质，其特征在于：所述步骤二中控制点的数量不少于所述步骤一中待定标影像的某一范围的像素点的万分之五。
7.根据权利要求5所述的一种计算机可读存储介质，其特征在于：所述步骤二中，采用灰度匹配和线特征匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点初步自动量测，然后采用最小二乘匹配方法对数字化几何定标场数据与步骤一中的影像进行控制点精确自动量测。
CN201811270813.7A 2018-10-29 2018-10-29 一种数字化几何定标场和gcp相结合的内方位元素定标方法 Active CN109143295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811270813.7A CN109143295B (zh) 2018-10-29 2018-10-29 一种数字化几何定标场和gcp相结合的内方位元素定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811270813.7A CN109143295B (zh) 2018-10-29 2018-10-29 一种数字化几何定标场和gcp相结合的内方位元素定标方法

Publications (2)

Publication Number Publication Date
CN109143295A CN109143295A (zh) 2019-01-04
CN109143295B true CN109143295B (zh) 2021-04-02

Family

ID=64806383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811270813.7A Active CN109143295B (zh) 2018-10-29 2018-10-29 一种数字化几何定标场和gcp相结合的内方位元素定标方法

Country Status (1)

Country Link
CN (1) CN109143295B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114241064B (zh) * 2022-02-24 2022-05-17 中国科学院空天信息创新研究院 一种遥感卫星内外方位元素实时几何定标方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901519A (zh) * 2012-11-02 2013-01-30 武汉大学 一种基于探元指向角光学推扫卫星在轨分步几何定标方法
CN106873004A (zh) * 2016-12-21 2017-06-20 中国资源卫星应用中心 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法
US20170307763A1 (en) * 2016-04-26 2017-10-26 Uber Technologies, Inc. Road registration differential gps

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103673995B (zh) * 2013-11-29 2016-09-21 航天恒星科技有限公司 一种线阵推扫式相机在轨光学畸变参数标定方法
US10248124B2 (en) * 2016-07-21 2019-04-02 Mobileye Vision Technologies, Inc. Localizing vehicle navigation using lane measurements

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901519A (zh) * 2012-11-02 2013-01-30 武汉大学 一种基于探元指向角光学推扫卫星在轨分步几何定标方法
US20170307763A1 (en) * 2016-04-26 2017-10-26 Uber Technologies, Inc. Road registration differential gps
CN106873004A (zh) * 2016-12-21 2017-06-20 中国资源卫星应用中心 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
On-Orbit Geometric Calibration of ZY-3 Three-Line Array Imagery With Multistrip Data Sets;Yongjun Zhang et al.;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20140131;第52卷(第1期);第224-234页 *
线阵推扫式相机高精度在轨几何标定;孟伟灿等;《武汉大学学报·信息科学版》;20151031;第40卷(第10期);第1392-1399、1413页 *

Also Published As

Publication number Publication date
CN109143295A (zh) 2019-01-04

Similar Documents

Publication Publication Date Title
CN109903352B (zh) 一种卫星遥感影像大区域无缝正射影像制作方法
KR101965965B1 (ko) 위성영상과 제공 rpc로부터 제작된 수치표고모델의 자동 오차보정 방법
CN107644435B (zh) 顾及姿态校正的敏捷光学卫星无场地几何标定方法及系统
CN104897175B (zh) 多相机光学推扫卫星在轨几何定标方法及系统
CN112017224B (zh) Sar数据区域网平差处理方法和系统
CN107564057B (zh) 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法
CN102901519B (zh) 一种基于探元指向角光学推扫卫星在轨分步几何定标方法
CN111044037B (zh) 一种光学卫星影像的几何定位方法及装置
CN108226982B (zh) 单线阵卫星激光联合高精度定位处理方法
CN111174753B (zh) 基于有理函数模型的光学影像与激光测高数据平差方法
CN109506627B (zh) 一种基于rpc模型的卫星影像系统误差修正方法
CN110986998B (zh) 一种基于有理函数模型的卫星视频相机在轨几何定标方法
CN111486864A (zh) 基于立体正八边结构的多源传感器联合标定方法
CN113538595B (zh) 利用激光测高数据辅助提升遥感立体影像几何精度的方法
CN114972078A (zh) 应用sar影像提升国产光学卫星影像无控几何质量方法及系统
CN109188483B (zh) 一种时序化高精度外方位元素自动定标方法
CN108562900B (zh) 一种基于高程校正的sar图像几何配准方法
CN114241064B (zh) 一种遥感卫星内外方位元素实时几何定标方法
CN111508028A (zh) 光学立体测绘卫星相机的自主在轨几何定标方法及系统
CN109143295B (zh) 一种数字化几何定标场和gcp相结合的内方位元素定标方法
CN114838740A (zh) 一种考虑不同经纬度区域的卫星图像几何定标方法
CN105510901A (zh) 基于多定标场的光学卫星影像时变误差定标方法及系统
CN111275773A (zh) 无场几何定标方法及系统
CN111044076B (zh) 基于参考底图的高分一号b卫星几何检校方法
CN115311366A (zh) 基于rpc模型的星载分片线阵传感器几何定标方法及系统

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