CN114166243B - 一种带有片间约束的面阵成像卫星在轨几何定标方法 - Google Patents
一种带有片间约束的面阵成像卫星在轨几何定标方法 Download PDFInfo
- Publication number
- CN114166243B CN114166243B CN202111317290.9A CN202111317290A CN114166243B CN 114166243 B CN114166243 B CN 114166243B CN 202111317290 A CN202111317290 A CN 202111317290A CN 114166243 B CN114166243 B CN 114166243B
- Authority
- CN
- China
- Prior art keywords
- area array
- sliced
- slice
- imaging
- inter
- 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
Images
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种带有片间约束的面阵成像卫星在轨几何定标方法,包括以下步骤:步骤1,通过影像匹配方法,分别获得每一片分片面阵影像上的控制点和相邻分片面阵影像之间的连接点;步骤2,构建分片面阵成像器件在轨几何定标模型;步骤3,构建相邻分片面阵器件之间的片间约束模型;步骤4,整体优化求解所有分片面阵器件的定标参数。本发明可以在分片影像边缘控制点缺乏或者分布不均匀的情况下,实现所有分片面阵器件定标参数的整体优化求解,进而提升相邻分片影像的几何拼接精度。
Description
技术领域
本发明属于光学遥感卫星数据处理相关技术领域,特别是涉及一种带有片间约束的面阵成像卫星在轨几何定标方法。
背景技术
为了满足大幅宽成像的应用需求,光学遥感卫星通常采用多片线阵或面阵器件拼接成像的方式,获取具有一定重叠度的多片分片影像。在卫星数据地面处理***中,再对多片分片影像进行拼接处理,得到一景完整的大幅宽影像。
对于线阵成像卫星,多片线阵器件通常以共线或者品字型方式安置于焦平面上,分片器件之间的重叠关系比较简单。而对于面阵成像卫星,多片面阵器件往往以2×2或者3×3等方式安置于焦平面上,分片器件之间的重叠关系要复杂得多。如何实现面阵成像卫星分片影像精确几何拼接,一直是光学遥感卫星数据处理技术领域亟待解决的一个关键问题。
精确获知所有面阵器件的成像参数,是实现分片影像精确几何拼接的前提和基础。在轨几何定标是当前获取面阵卫星精确成像参数的最常用方法之一,已广泛应用于卫星数据地面处理***中,并取得了不错的效果。该方法根据面阵成像卫星的成像几何关系构建在轨几何定标模型,并利用地面控制点分别求解出每一片面阵器件的定标参数。
现有方法使用的地面控制点主要来源于高精度参考数据(数字正射影像和数字高程模型),这些控制点需要在所有分片影像上均匀分布,特别是在影像边缘区域,需要有足够数量均匀分布的控制点,才能够形成有效的边缘控制约束,从而保证相邻分片影像之间的几何拼接精度。然而,受参考数据内地物变化、参考数据和卫星影像之间辐射差异等因素的影响,往往很难在所有分片影像的边缘区域均能够获得均匀分布的控制点,这必然导致难以精确求解出每一片面阵器件的定标参数,以至于无法保证相邻分片影像之间的几何拼接精度,降低拼接后卫星影像产品的几何质量。
发明内容
本发明针对现有技术的不足,提供一种带有片间约束的面阵成像卫星在轨几何定标方法,以分片面阵影像上的地面控制点作为绝对控制约束,以相邻分片影像之间连接点形成的片间约束作为相对控制约束,整体优化求解所有分片面阵器件的定标参数,提高了相邻分片影像之间的几何拼接精度,解决了现有技术中分片面阵影像边缘区域内控制点缺乏或分布不均匀导致的定标参数难以精确求解的问题。
为了达到上述目的,本发明提供的技术方案是一种带有片间约束的面阵成像卫星在轨几何定标方法,包括以下步骤:
步骤1,对面阵影像进行影像匹配,分别获得每一片分片面阵影像上的控制点和相邻分片面阵影像之间的连接点;
步骤2,根据面阵成像卫星的成像几何关系,构建分片面阵器件在轨几何定标模型;
步骤3,根据面阵成像卫星面阵器件瞬时成像和多片面阵器件同时成像的特点,构建相邻分片面阵器件之间的片间约束模型;
步骤4,利用步骤1得到的所有分片面阵影像上的控制点和相邻分片面阵影像之间的连接点,整体优化求解所有分片面阵器件的定标参数;
步骤4.1,根据步骤2构建分片面阵器件在轨几何定标模型,针对每一片分片面阵影像上的每一个控制点构建误差方程;
步骤4.2,根据步骤3构建相邻分片面阵器件之间的片间约束模型,针对相邻分片影像之间的每一个连接点构建误差方程;
步骤4.3,联立步骤4.1构建的误差方程和步骤4.2构建的误差方程,根据最小二乘平差原理求解得到所有分片面阵器件的定标参数。
而且,所述步骤2中根据面阵成像卫星的成像几何关系,构建分片面阵器件在轨几何定标模型,如式(1)所示:
式中,(XGNSS,YGNSS,ZGNSS)WGS84为GNSS天线相位中心在WGS84坐标系下的空间直角坐标;(X,Y,Z)WGS84为地面点在WGS84坐标系下的空间直角坐标;为WGS84坐标系至J2000坐标系的旋转矩阵;为J2000坐标系至卫星姿态测量坐标系的旋转矩阵;λ为比例因子;(φx,φy)为地面点对应的成像探元在卫星姿态测量坐标系下的指向角;(s,l)为探元编号;(a0,a1,...,a9,b0,b1,...,b9)为定标参数。
而且,所述步骤3中根据面阵成像卫星面阵器件瞬时成像和多片面阵器件同时成像的特点,构建相邻分片面阵器件之间的片间约束模型,如式(3)所示:
式中,(s,l)为探元编号,(a0,p,a1,p,...,a9,p,b0,p,b1,p,...,b9,p)为第p片面阵成像器件的定标参数,(a0,q,a1,q,...,a9,q,b0,q,b1,q,...,b9,q)为第q片面阵成像器件的定标参数,第p片和第q片面阵器件为相机焦平面上相邻的两片器件。
而且,所述步骤4.1中针对每一片分片面阵影像上的每一个控制点,由式(1)求解得到(tanφx,g,m,ni,tanφy,g,m,ni),并根据式(2)构建误差方程:
Vg=AgX-Lg (4)
式中,
X=[a0,1 a1,1 … a9,1 b0,1 b1,1 … b9,1 … a0,m a1,m … a9,m b0,m b1,m … b9,m]T为未知数矩阵;下标g表示控制点;下标m为分片影像数量;下标ni=n1,n2,…,nm表示第i片分片影像上的控制点数量。
而且,所述步骤4.2中针对相邻分片影像之间的每一个连接点,根据式(3)构建误差方程:
Vt=AtX-Lt (5)
式中,
X=[… a0,p a1,p … a9,p b0,p b1,p … b9,p a0,q a1,q … a9,q b0,q b1,q … b9,q…]T
为未知数矩阵;下标t表示连接点,下标np表示第p片和第q片分片影像之间的连接点数量。
而且,所述步骤4.3中联立式(4)和式(5),根据最小二乘平差原理求解未知数矩阵X:
由此得到所有分片面阵器件的定标参数,面阵成像卫星的在轨几何定标工作完成。
与现有技术相比,本发明具有如下优点:本发明除了以分片面阵影像上的地面控制点作为绝对控制约束之外,还引入了相邻分片影像之间连接点形成的片间约束,将其作为相对控制约束,并整体优化求解所有分片面阵器件的定标参数,进而解决分片面阵影像边缘区域内控制点缺乏或分布不均匀导致的定标参数难以精确求解的问题,提高了相邻分片面阵影像的几何拼接精度,进而提升卫星影像产品的几何质量。
附图说明
图1为本发明实施例的流程框图。
具体实施方式
本发明提供一种带有片间约束的面阵成像卫星在轨几何定标方法,首先以分片面阵影像上的地面控制点作为绝对控制约束,然后以相邻分片影像之间连接点形成的片间约束作为相对控制约束,整体优化求解所有分片面阵器件的定标参数,以提高相邻分片影像之间的几何拼接精度。
下面结合附图和实施例对本发明的技术方案作进一步说明。
如图1所示,本发明实施例的流程包括以下步骤:
步骤1,对面阵影像进行影像匹配,分别获得每一片分片面阵影像上的控制点和相邻分片面阵影像之间的连接点。
步骤2,根据面阵成像卫星的成像几何关系,构建分片面阵器件在轨几何定标模型,如式(1)所示:
式中,(XGNSS,YGNSS,ZGNSS)WGS84为GNSS天线相位中心在WGS84坐标系下的空间直角坐标;(X,Y,Z)WGS84为地面点在WGS84坐标系下的空间直角坐标;为WGS84坐标系至J2000坐标系的旋转矩阵;为J2000坐标系至卫星姿态测量坐标系的旋转矩阵;λ为比例因子;(φx,φy)为地面点对应的成像探元在卫星姿态测量坐标系下的指向角;(s,l)为探元编号;(a0,a1,...,a9,b0,b1,...,b9)为定标参数。
步骤3,根据面阵成像卫星面阵器件瞬时成像和多片面阵器件同时成像的特点,构建相邻分片面阵器件之间的片间约束模型,如式(3)所示:
式中,(s,l)为探元编号,(a0,p,a1,p,...,a9,p,b0,p,b1,p,...,b9,p)为第p片面阵成像器件的定标参数;(a0,q,a1,q,...,a9,q,b0,q,b1,q,...,b9,q)为第q片面阵成像器件的定标参数;第p片和第q片面阵器件为相机焦平面上相邻的两片器件。
步骤4,利用步骤1得到的所有分片面阵影像上的控制点和相邻分片面阵影像之间的连接点,整体优化求解所有分片面阵器件的定标参数。
Vg=AgX-Lg (4)
式中,
X=[a0,1 a1,1 … a9,1 b0,1 b1,1 … b9,1 … a0,m a1,m … a9,m b0,m b1,m … b9,m]T为未知数矩阵;下标g表示控制点;下标m为分片影像数量;下标ni=n1,n2,…,nm表示第i片分片影像上的控制点数量。
步骤4.2,针对相邻分片影像之间的每一个连接点,根据式(3)构建误差方程:
Vt=AtX-Lt (5)
式中,
X=[… a0,p a1,p … a9,p b0,p b1,p … b9,p a0,q a1,q … a9,q b0,q b1,q … b9,q…]T为未知数矩阵;下标t表示连接点,下标np表示第p片和第q片分片影像之间的连接点数量。
步骤4.3,联立式(4)和式(5),根据最小二乘平差原理求解未知数矩阵X:
由此得到所有分片面阵器件的定标参数,面阵成像卫星的在轨几何定标工作完成。
具体实施时,以上流程可采用计算机软件技术实现自动运行流程。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种带有片间约束的面阵成像卫星在轨几何定标方法,其特征在于,包括如下步骤:
步骤1,对面阵影像进行影像匹配,分别获得每一片分片面阵影像上的控制点和相邻分片面阵影像之间的连接点;
步骤2,根据面阵成像卫星的成像几何关系,构建分片面阵器件在轨几何定标模型;
步骤3,根据面阵成像卫星面阵器件瞬时成像和多片面阵器件同时成像的特点,构建相邻分片面阵器件之间的片间约束模型,即:
式中,(s,l)为探元编号,(a0,p,a1,p,...,a9,p,b0,p,b1,p,...,b9,p)为第p片面阵成像器件的定标参数,(a0,q,a1,q,...,a9,q,b0,q,b1,q,...,b9,q)为第q片面阵成像器件的定标参数,第p片和第q片面阵器件为相机焦平面上相邻的两片器件;
步骤4,利用步骤1得到的所有分片面阵影像上的控制点和相邻分片面阵影像之间的连接点,整体优化求解所有分片面阵器件的定标参数;
步骤4.1,根据步骤2构建分片面阵器件在轨几何定标模型,针对每一片分片面阵影像上的每一个控制点构建误差方程;
步骤4.2,根据步骤3构建相邻分片面阵器件之间的片间约束模型,针对相邻分片影像之间的每一个连接点构建误差方程;
根据式(3)构建误差方程:
Vt=AtX-Lt (5)
式中,
X=[…a0,p a1,p…a9,p b0,p b1,p…b9,p a0,q a1,q…a9,q b0,q b1,q…b9,q…]T
为未知数矩阵;下标t表示连接点,下标np表示第p片和第q片分片影像之间的连接点数量;
步骤4.3,联立步骤4.1构建的误差方程和步骤4.2构建的误差方程,根据最小二乘平差原理求解得到所有分片面阵器件的定标参数。
2.如权利要求1所述的一种带有片间约束的面阵成像卫星在轨几何定标方法,其特征在于:步骤2中根据面阵成像卫星的成像几何关系,构建分片面阵器件在轨几何定标模型,如式(1)所示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111317290.9A CN114166243B (zh) | 2021-11-09 | 2021-11-09 | 一种带有片间约束的面阵成像卫星在轨几何定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111317290.9A CN114166243B (zh) | 2021-11-09 | 2021-11-09 | 一种带有片间约束的面阵成像卫星在轨几何定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114166243A CN114166243A (zh) | 2022-03-11 |
CN114166243B true CN114166243B (zh) | 2023-04-18 |
Family
ID=80478329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111317290.9A Active CN114166243B (zh) | 2021-11-09 | 2021-11-09 | 一种带有片间约束的面阵成像卫星在轨几何定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114166243B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI124586B (fi) * | 2012-02-03 | 2014-10-31 | Suunto Oy | Menetelmä ja laite korkeuden määrittämiseksi |
CN111508029A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 星载分片线阵ccd光学相机整体几何定标方法及*** |
CN112802118B (zh) * | 2021-01-05 | 2022-04-08 | 湖北工业大学 | 一种光学卫星传感器在轨分时几何定标方法 |
CN113313769B (zh) * | 2021-06-11 | 2022-08-30 | 湖北工业大学 | 一种光学卫星多面阵传感器片间无缝几何定标方法 |
-
2021
- 2021-11-09 CN CN202111317290.9A patent/CN114166243B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114166243A (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110006452B (zh) | 高分六号宽视场相机相对几何定标方法及*** | |
CN113313769B (zh) | 一种光学卫星多面阵传感器片间无缝几何定标方法 | |
CN105606906B (zh) | 一种毫米波相控阵测试标定方法 | |
CN105444778B (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN106887016B (zh) | 一种gf-4卫星序列图像自动相对配准方法 | |
CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN100376883C (zh) | 一种基于像素频率的星敏感器校准方法 | |
CN110363758B (zh) | 一种光学遥感卫星成像质量确定方法及*** | |
CN110297220B (zh) | 一种双天线InSAR***基线矢量的测量方法 | |
CN109959898B (zh) | 一种座底式水声被动定位基阵自校准方法 | |
CN109191532B (zh) | 一种机载三线阵ccd相机标定方法 | |
CN109544642B (zh) | 一种基于n型靶标的tdi-ccd相机参数标定方法 | |
CN107316280B (zh) | 离岛卫星影像rpc模型高精度几何定位方法 | |
US5550376A (en) | Method of calibration of imaging devices | |
CN108447100B (zh) | 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法 | |
CN110793542A (zh) | 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法 | |
CN105571598B (zh) | 一种卫星激光高度计足印相机姿态的测定方法 | |
CN114166243B (zh) | 一种带有片间约束的面阵成像卫星在轨几何定标方法 | |
CN104537614B (zh) | 一种环境一号卫星ccd影像正射校正方法 | |
CN115865142B (zh) | 一种雷达天线波束指向校准方法、***、装置及存储介质 | |
CN107941201B (zh) | 光线同姿约束的零交会光学卫星影像联合平差方法及*** | |
CN109696155B (zh) | 光线共面约束的弱交会光学卫星影像联合平差方法及*** | |
CN114594435A (zh) | 一种国产民用sar卫星几何检校与定位精度改进方法 | |
CN110044379B (zh) | 一种移动测量***标定方法 | |
CN115164938B (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 |