CN110986998A - 一种基于有理函数模型的卫星视频相机在轨几何定标方法 - Google Patents
一种基于有理函数模型的卫星视频相机在轨几何定标方法 Download PDFInfo
- Publication number
- CN110986998A CN110986998A CN201911031231.8A CN201911031231A CN110986998A CN 110986998 A CN110986998 A CN 110986998A CN 201911031231 A CN201911031231 A CN 201911031231A CN 110986998 A CN110986998 A CN 110986998A
- Authority
- CN
- China
- Prior art keywords
- video camera
- image
- geometric calibration
- satellite video
- model
- 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.)
- Granted
Links
- 238000012892 rational function Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012937 correction Methods 0.000 claims description 18
- 238000003384 imaging method Methods 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 12
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 11
- 238000013519 translation Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 description 5
- 239000013598 vector Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Studio Devices (AREA)
Abstract
本发明公开了一种基于有理函数模型的卫星视频相机在轨几何定标方法,包括:步骤1,根据卫星视频影像的定向参数,构建有理函数模型;步骤2,根据卫星视频相机的物理结构特点,构建基于有理函数模型的在轨几何定标模型;步骤3,以视频影像覆盖地面范围内的数字正射影像和数字高程模型作为参考数据,利用影像匹配提取地面控制点;步骤4,根据卫星视频相机的在轨几何定标模型,利用地面控制点,按照最小二乘平差原理,求解卫星视频相机的几何定标参数,实现卫星视频相机在轨几何定标。本发明以有理函数模型为基础,构建卫星视频相机在轨几何定标模型,可以为用户提供一种简单、易行的卫星视频相机在轨几何定标方法。
Description
技术领域
本发明属于卫星视频影像几何处理技术领域,特别涉及一种基于有理函数模型的卫星视频相机在轨几何定标方法。
背景技术
视频卫星是近几年兴起的一种新型遥感卫星,星上搭载的视频相机不仅具有获取图像的能力,更大的优势是能够以视频录像的方式,获取比传统遥感卫星更多动态信息的视频数据,可以满足对兴趣目标实时、动态监测的需求。
卫星视频影像的高精度几何定位是其精确、可靠应用的前提。为了实现卫星视频影像的高精度几何定位,首先需求解决卫星视频相机的在轨几何定标问题。现有卫星相机在轨几何定标方法以严密传感器模型为基础,建立在轨几何定标模型,这涉及一系列坐标系的空间坐标转换,如卫星本体坐标系、星敏感器坐标系和J2000坐标系。由此可见,利用现有在轨几何定标方法需要很强的专业知识,不便于用户开展在轨几何定标工作。另一方面,为了建立基于严密传感器模型的在轨几何定标模型,必须获得相机成像过程中的各种特征参数,如卫星位置、卫星姿态及其物理特性参数等,而用户有时无法得到这些参数,以至于无法开展在轨几何定标工作。
发明内容
针对现有技术在卫星视频相机在轨几何定标方面存在的不足,本发明提供了一种基于有理函数模型的卫星视频相机在轨几何定标方法,本发明以有理函数模型为基础,构建在轨几何定标模型,可以为用户提供一种简单、易行的卫星视频相机在轨几何定标方法。
本发明提供的一种基于有理函数模型的卫星视频相机在轨几何定标方法,包括以下步骤:
步骤1,根据卫星视频影像的定向参数,构建有理函数模型;
步骤2,根据卫星视频相机的物理结构特点,构建基于有理函数模型的在轨几何定标模型,实现如下,
构建卫星视频相机的成像误差修正模型:
其中,(Δx1,Δy1)为卫星视频相机自身误差的修正值,(x,y)为像点坐标,(s0,s1,...,s7,t0,t1,...,t7)为卫星视频相机自身误差的修正参数;
构建卫星视频相机的在轨几何定标模型:
其中,(x,y)为像点坐标,(xoffset,yoffset)为像点坐标正则化平移参数,(xscale,yscale)为像点坐标正则化缩放参数,为地面点正则化坐标,p1,p2,p3,p4为有理函数模型中的三次多项式;(β0,β1,...,β7,θ0,θ1,...,θ7)为几何定标参数;
步骤3,以视频影像覆盖地面范围内的数字正射影像和数字高程模型作为参考数据,利用影像匹配提取地面控制点;
步骤4,根据卫星视频相机的在轨几何定标模型,利用地面控制点,按照最小二乘平差原理,求解卫星视频相机的几何定标参数,实现卫星视频相机在轨几何定标。
而且,步骤1中,所述有理函数模型为
其中,(x,y)为像点坐标;(xoffset,yoffset)为像点坐标正则化平移参数;(xscale,yscale)为像点坐标正则化缩放参数;为地面点正则化坐标;三次多项式p1,p2,p3,p4的具体形式为:
其中,(a1,a2,...,a20,b1,b2,...,b20,c1,c2,...,c20,d1,d2,...,d20)为卫星视频影像的定向参数。
而且,步骤3的实现方式包括以下步骤,
步骤3.1,利用SIFT算法在视频影像和数字正射影像上进行影像匹配,得到初始匹配点;
步骤3.2,利用仿射变换模型描述视频影像和数字正射影像之间的几何关系,并利用初始匹配点求解仿射变换系数;
步骤3.3,以仿射变换系数作为最小二乘影像匹配中几何纠正模型系数的初值,进行最小二乘影像匹配,得到精确匹配点;
步骤3.4,以视频影像上匹配点的像方坐标作为地面控制点的像方坐标,并根据数字正射影像上对应匹配点的像点坐标及数字高程模型计算其地理坐标,作为控制点的地理坐标。
而且,步骤4的实现方式包括以下步骤,
步骤4.1,针对每一个地面控制点,根据步骤3所得卫星视频相机在轨几何定标模型,建立误差方程式;
步骤4.2,按照最小二乘平差原理,求解几何定标参数;
步骤4.3,计算每一个地面控制点对应像点坐标的残差,并统计所有像点坐标残差的中误差;
步骤4.4,剔除像点坐标残差大于3倍中误差的地面控制点,并重复步骤4.1至4.3,直至所有像点坐标的残差小于3倍中误差时,将最新得到的几何定标参数作为视频相机的精确几何定标参数。
本发明以有理函数模型为基础,构建卫星视频相机在轨几何定标模型,可以为用户提供一种简单、易行的卫星视频相机在轨几何定标方法。
附图说明
图1为本发明实施例的具体流程图。
具体实施方式
为了更清楚地说明本发明实施例和/或现有技术中的技术方案,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
本发明实施例在卫星视频影像有理函数模型的基础上,构建一种基于有理函数模型的卫星视频相机在轨几何定标模型,并利用地面控制点精确求解几何定标参数,因而无需构建视频影像的严密成像几何模型,也无需用户具备很强的专业知识,非常便于用户开展卫星视频相机在轨几何定标工作。参见图1,实施例提供的基于有理函数模型的卫星视频相机在轨几何定标方法具体步骤如下:
步骤1,根据卫星视频影像的定向参数,构建有理函数模型,如式(1)所示:
其中,(x,y)为像点坐标;(xoffset,yoffset)为像点坐标正则化平移参数;(xscale,yscale)为像点坐标正则化缩放参数;为地面点正则化坐标;三次多项式p1,p2,p3,p4的具体形式为:
其中,(a1,a2,...,a20,b1,b2,...,b20,c1,c2,...,c20,d1,d2,...,d20)为卫星视频影像的定向参数。
步骤2,根据卫星视频相机物理结构特点,构建基于有理函数模型的在轨几何定标模型。
本步骤进一步包括:
2.1,构建卫星视频相机的成像误差修正模型。
卫星视频相机在成像过程中,主要受卫星位置与姿态观测误差、卫星相机与卫星定姿设备之间的集成误差以及卫星相机自身误差的影响。
卫星相机自身的误差主要包括主点、主距和镜头畸变误差,可由式(2)进行成像误差补偿:
其中,(Δx1,Δy1)为卫星视频相机自身误差的修正值;(x0,y0)为相机主点;f为相机主距;Δf为相机主距修正值;k1为镜头畸变修正参数。
本发明提出,将式(2)进行展开,可以得到:
其中,(s0,s1,...,s7,t0,t1,...,t7)为卫星视频相机自身误差的修正参数。
卫星位置与姿态观测误差、卫星相机与卫星定姿设备之间的集成误差可以通过像方仿射变换模型进行成像误差补偿,如式(4)所示:
其中,(Δx2,Δy2)为观测误差与集成误差修正值;(e1,e2,e3,f1,f2,f3)为观测误差与集成误差修正参数。
本发明提出,综合式(3)和(4),可以构建卫星视频相机的成像误差修正模型,如式(5)所示:
其中,(Δx,Δy)为成像误差修正值;(β0,β1,...,β7,θ0,θ1,...,θ7)为成像误差修正参数,这些参数综合描述了卫星位置与姿态观测误差、卫星相机与定姿设备集成误差以及卫星相机自身误差对卫星视频影像定位精度的影响。
2.2,构建卫星视频相机的在轨几何定标模型。
在有理函数模型(式(1))的基础上,引入卫星视频相机成像误差修正模型,即可构建卫星视频相机的在轨几何定标模型,如式(6)所示:
其中,(β0,β1,...,β7,θ0,θ1,...,θ7)即为几何定标参数。
步骤3,以视频影像覆盖地面范围内的数字正射影像和数字高程模型作为参考数据,利用影像匹配方法,提取地面控制点。
本步骤进一步包括:
3.1,利用SIFT算法在视频影像和数字正射影像上进行影像匹配,得到初始匹配点;
3.2,利用仿射变换模型描述视频影像和数字正射影像之间的几何关系,并利用初始匹配点求解仿射变换系数;
3.3,以仿射变换系数作为最小二乘影像匹配中几何纠正模型系数的初值,进行最小二乘影像匹配,得到精确匹配点;
3.4,以视频影像上匹配点的像方坐标作为地面控制点的像方坐标,并根据数字正射影像上对应匹配点的像点坐标及数字高程模型计算其地理坐标,作为控制点的地理坐标。
SIFT算法、仿射变换模型和最小二乘影像匹配具体实现为现有技术,本发明不予赘述,本发明利用这些方法得到视频影像和数字正射影像之间的匹配点。
步骤4,根据卫星视频相机的在轨几何定标模型,利用地面控制点,按照最小二乘平差原理,求解卫星视频相机的几何定标参数,实现卫星视频相机在轨几何定标。
本步骤进一步包括:
4.1,针对每一个地面控制点,根据如式(6)所示卫星视频相机在轨几何定标模型,建立误差方程式,如式(7)所示:
v=Bt-l (7)
t=[β0 β1 β2 β3 β4 β5 β6 β7 θ0 θ1 θ2 θ3 θ4 θ5 θ6 θ7]T为未知数向量。
4.2,按照最小二乘平差原理,求解几何定标参数,如式(8)所示:
t=(BTB)-1BTl (8)
4.3,计算每一个地面控制点对应像点坐标的残差,并统计所有像点坐标残差的中误差,如式(9)所示:
其中,vxi和vyi分别表示第i个像点在x和y方向上的坐标残差;σx和σy分别表示x和y方向上的像点坐标残差中误差;n表示控制点数量。
4.4,剔除像点坐标残差大于3倍中误差的地面控制点,并重复步骤4.1至4.3,直至所有像点坐标的残差小于3倍中误差,最近一次执行步骤4.2得到的几何定标参数即为视频相机的精确几何定标参数。
具体实施时,可采用软件方式实现流程的自动运行。运行流程的装置也应当在本发明的保护范围内。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明做出任何的修改和改变,都落入本发明的保护范围。
Claims (4)
1.一种基于有理函数模型的卫星视频相机在轨几何定标方法,其特征在于,包括以下步骤:
步骤1,根据卫星视频影像的定向参数,构建有理函数模型;
步骤2,根据卫星视频相机的物理结构特点,构建基于有理函数模型的在轨几何定标模型,实现如下,
构建卫星视频相机的成像误差修正模型:
其中,(Δx1,Δy1)为卫星视频相机自身误差的修正值,(x,y)为像点坐标,(s0,s1,...,s7,t0,t1,...,t7)为卫星视频相机自身误差的修正参数;
构建卫星视频相机的在轨几何定标模型:
其中,(x,y)为像点坐标,(xoffset,yoffset)为像点坐标正则化平移参数,(xscale,yscale)为像点坐标正则化缩放参数,为地面点正则化坐标,p1,p2,p3,p4为有理函数模型中的三次多项式;(β0,β1,...,β7,θ0,θ1,...,θ7)为几何定标参数;
步骤3,以视频影像覆盖地面范围内的数字正射影像和数字高程模型作为参考数据,利用影像匹配提取地面控制点;
步骤4,根据卫星视频相机的在轨几何定标模型,利用地面控制点,按照最小二乘平差原理,求解卫星视频相机的几何定标参数,实现卫星视频相机在轨几何定标。
3.如权利要求1或2所述的一种基于有理函数模型的卫星视频相机在轨几何定标方法,其特征在于:步骤3的实现方式包括以下步骤,
步骤3.1,利用SIFT算法在视频影像和数字正射影像上进行影像匹配,得到初始匹配点;
步骤3.2,利用仿射变换模型描述视频影像和数字正射影像之间的几何关系,并利用初始匹配点求解仿射变换系数;
步骤3.3,以仿射变换系数作为最小二乘影像匹配中几何纠正模型系数的初值,进行最小二乘影像匹配,得到精确匹配点;
步骤3.4,以视频影像上匹配点的像方坐标作为地面控制点的像方坐标,并根据数字正射影像上对应匹配点的像点坐标及数字高程模型计算其地理坐标,作为控制点的地理坐标。
4.如权利要求1或2所述的一种基于有理函数模型的卫星视频相机在轨几何定标方法,其特征在于:步骤4的实现方式包括以下步骤,
步骤4.1,针对每一个地面控制点,根据步骤3所得卫星视频相机在轨几何定标模型,建立误差方程式;
步骤4.2,按照最小二乘平差原理,求解几何定标参数;
步骤4.3,计算每一个地面控制点对应像点坐标的残差,并统计所有像点坐标残差的中误差;
步骤4.4,剔除像点坐标残差大于3倍中误差的地面控制点,并重复步骤4.1至4.3,直至所有像点坐标的残差小于3倍中误差时,将最新得到的几何定标参数作为视频相机的精确几何定标参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031231.8A CN110986998B (zh) | 2019-10-28 | 2019-10-28 | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031231.8A CN110986998B (zh) | 2019-10-28 | 2019-10-28 | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110986998A true CN110986998A (zh) | 2020-04-10 |
CN110986998B CN110986998B (zh) | 2021-09-14 |
Family
ID=70082560
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911031231.8A Active CN110986998B (zh) | 2019-10-28 | 2019-10-28 | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110986998B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111798523A (zh) * | 2020-06-19 | 2020-10-20 | 武汉大学 | 星相机在轨定标定姿及遥感影像几何定位方法、*** |
CN112082576A (zh) * | 2020-09-09 | 2020-12-15 | 桂林理工大学 | 一种“三无情况下”的卫星影像正射纠正方法 |
CN112924465A (zh) * | 2021-02-02 | 2021-06-08 | 长江空间信息技术工程有限公司(武汉) | 一种地下井室缺陷检测方法及其对应装置 |
CN114241064A (zh) * | 2022-02-24 | 2022-03-25 | 中国科学院空天信息创新研究院 | 一种遥感卫星内外方位元素实时几何定标方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104897175A (zh) * | 2015-06-23 | 2015-09-09 | 武汉大学 | 多相机光学推扫卫星在轨几何定标方法及*** |
CN105091906A (zh) * | 2015-06-30 | 2015-11-25 | 武汉大学 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及*** |
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及*** |
CN107144293A (zh) * | 2017-04-07 | 2017-09-08 | 武汉大学 | 一种视频卫星面阵相机的几何定标方法 |
CN108830889A (zh) * | 2018-05-24 | 2018-11-16 | 中国科学院遥感与数字地球研究所 | 基于全局几何约束的遥感影像与基准影像的匹配方法 |
-
2019
- 2019-10-28 CN CN201911031231.8A patent/CN110986998B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104897175A (zh) * | 2015-06-23 | 2015-09-09 | 武汉大学 | 多相机光学推扫卫星在轨几何定标方法及*** |
CN105091906A (zh) * | 2015-06-30 | 2015-11-25 | 武汉大学 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及*** |
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及*** |
CN107144293A (zh) * | 2017-04-07 | 2017-09-08 | 武汉大学 | 一种视频卫星面阵相机的几何定标方法 |
CN108830889A (zh) * | 2018-05-24 | 2018-11-16 | 中国科学院遥感与数字地球研究所 | 基于全局几何约束的遥感影像与基准影像的匹配方法 |
Non-Patent Citations (1)
Title |
---|
曹金山等: "资源三号影像对地目标定位的***误差补偿", 《测绘科学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111798523A (zh) * | 2020-06-19 | 2020-10-20 | 武汉大学 | 星相机在轨定标定姿及遥感影像几何定位方法、*** |
CN112082576A (zh) * | 2020-09-09 | 2020-12-15 | 桂林理工大学 | 一种“三无情况下”的卫星影像正射纠正方法 |
CN112924465A (zh) * | 2021-02-02 | 2021-06-08 | 长江空间信息技术工程有限公司(武汉) | 一种地下井室缺陷检测方法及其对应装置 |
CN112924465B (zh) * | 2021-02-02 | 2023-01-17 | 长江空间信息技术工程有限公司(武汉) | 一种地下井室缺陷检测方法及其对应装置 |
CN114241064A (zh) * | 2022-02-24 | 2022-03-25 | 中国科学院空天信息创新研究院 | 一种遥感卫星内外方位元素实时几何定标方法 |
CN114241064B (zh) * | 2022-02-24 | 2022-05-17 | 中国科学院空天信息创新研究院 | 一种遥感卫星内外方位元素实时几何定标方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110986998B (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110986998B (zh) | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 | |
CN105096329B (zh) | 一种精确校正超广角摄像头图像畸变的方法 | |
CN108830889B (zh) | 基于全局几何约束的遥感影像与基准影像的匹配方法 | |
CN107644435B (zh) | 顾及姿态校正的敏捷光学卫星无场地几何标定方法及*** | |
CN111524194B (zh) | 一种激光雷达和双目视觉相互融合的定位方法及终端 | |
CN112419380B (zh) | 一种基于云掩膜的静止轨道卫星序列影像高精度配准方法 | |
CN106885585B (zh) | 一种基于光束法平差的星载摄影测量***一体化检校方法 | |
CN115578448B (zh) | 一种基于ccd观测数据批处理的天文定位方法及*** | |
CN107040695B (zh) | 基于rpc定位模型的星载视频稳像方法及*** | |
CN109523585A (zh) | 一种基于方向相位一致性的多源遥感影像特征匹配方法 | |
CN114241064B (zh) | 一种遥感卫星内外方位元素实时几何定标方法 | |
CN112598608A (zh) | 一种基于目标区域的光学卫星快速融合产品制作方法 | |
CN114463442A (zh) | 一种非同轴相机的标定方法 | |
CN110006452A (zh) | 高分六号宽视场相机相对几何定标方法及*** | |
JP2017130067A (ja) | 衛星映像の位置正確度改善のための自動映像処理システム及びその方法 | |
CN114754798A (zh) | 一种陀螺误差特性参数在轨辨识与标定方法 | |
JP2020038550A (ja) | ラインセンサカメラのキャリブレーション装置及び方法 | |
CN111275773A (zh) | 无场几何定标方法及*** | |
CN112017108B (zh) | 基于独立模型法平差的卫星影像颜色相对校正方法 | |
CN111044076B (zh) | 基于参考底图的高分一号b卫星几何检校方法 | |
JP4523422B2 (ja) | 衛星画像の位置補正装置 | |
CN109143295B (zh) | 一种数字化几何定标场和gcp相结合的内方位元素定标方法 | |
CN107644444B (zh) | 一种基于压缩感知的单图像相机标定方法 | |
CN114549648B (zh) | 一种基于高分遥感卫星立体影像rpc模型定位*** | |
CN116091546B (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 |