CN111324857B - 一种基于tdiccd推扫特性的快速反变换计算方法 - Google Patents
一种基于tdiccd推扫特性的快速反变换计算方法 Download PDFInfo
- Publication number
- CN111324857B CN111324857B CN202010197527.3A CN202010197527A CN111324857B CN 111324857 B CN111324857 B CN 111324857B CN 202010197527 A CN202010197527 A CN 202010197527A CN 111324857 B CN111324857 B CN 111324857B
- Authority
- CN
- China
- Prior art keywords
- imaging
- time
- ccd
- coordinate system
- satellite
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于TDICCD推扫特性的快速反变换计算方法,针对TDICCD推扫影像的延迟积分成像特征,提出了基于视场中心视轴地表轨迹与CCD阵列“近似垂直”约束的反变换算法。TDICCD正常成像条件下其视场中心视轴在地面的轨迹与CCD方向应该近似垂直的工作原理,推算出近似成像时间t,然后通过纠正镜头畸变和补偿偏流角残差,得到成像时刻,完成反算。可以达到提高计算效率和精度的技术效果。
Description
技术领域
本发明涉及航空航天技术领域,具体涉及一种基于TDICCD推扫特性的快速反变换计算方法。
背景技术
光学遥感卫星能够长时间、周期性地对地球成像,具备数据获取快速、不受区域限制的优势,已经成为对地观测的重要手段。目前,国内外光学遥感卫星主要采用TDICCD推扫成像体制来获取高分辨率和高信噪比的遥感影像,并通过对遥感影像进行几何校正来求解像素点的地理位置,且通过实现精确的偏流角调整来保障成像像质。当用TDICCD成像模式进行成像时,由地面坐标计算影像坐标的反变换需采用迭代求解算法来实现,过程复杂耗时,且存在不收敛情况。
本申请发明人在实施本发明的过程中,发现现有技术的方法,至少存在如下技术问题:
目前,反变换主要基于严格几何模型和有理多项式模型实现。张过等人在SPOT反变换算法基础上,依据稳态推扫成像特征将3次多项式简化为仿射模型,降低反变换计算量以提升效率。但该类算法由于求解格网多项式或仿射模型过程均涉及到大量的正变换计算,效率均较低。赵双明等人针对ADS40影像的几何预处理开展了研究,提出了一种基于共线方程的扫描行搜索算法来实现反变换。王密等、耿迅等学者针对基于共线方程的扫描行搜索算法计算量大的问题,研究并提出了基于物方的扫描行快速搜索算法。该算法可适用于稳态推扫成像的卫星,但对动中成像(尤其是曲线成像)卫星会因假设条件不成立而难以适用。基于有理多项式模型的反变换可直接根据定义式计算像方坐标,但要求首先生成能精确替代严格几何模型的RFM。金淑英等人分析了积分时间跳变对RFM替代精度的影响,提出积分时间相等的虚拟扫描行方法生成高替代精度的RFM,进而完成反变换。但该方法仅考虑了积分时间,忽略了姿态抖动对RFM替代精度的影响,适用性较低。
由此可知,现有的反变换方法存在计算复杂且耗时较长的技术问题。
发明内容
本发明提出一种基于TDICCD推扫特性的快速反变换计算方法,用于解决或者至少部分解决现有的反变换方法存在计算复杂且耗时较长的技术问题。
本发明的目的是通过下述技术方案实现的:
一种基于TDICCD推扫特性的快速反算计算方法,根据TDICCD正常成像条件下其视场中心视轴在地面的轨迹与CCD方向应该近似垂直的工作原理,推算出近似成像时间t,然后通过纠正镜头畸变和补偿偏流角残差,得到成像时刻,完成反算,得到影像坐标。
结合图1~图4,本发明提供的一种基于TDICCD推扫特性的快速反算计算方法,包括如下步骤:
S1:利用卫星下传的姿态、轨道数据、成像时间数据和相机参数,建立几何定位模型正变换;
S2:基于建立的几何定位模型的正变换计算CCD中心视轴的地表轨迹,其中,CCD中心视轴的地表轨迹与时间t具有函数关系;
S3:对于给定的地物点坐标,根据CCD中心视轴的地表轨迹与CCD方向之间关系、CCD中心视轴的地表轨迹与时间t具有的函数关系,获取地物点的近似成像时刻T0;
S4:根据S3中计算的成像时刻T0内插卫星成像位置、姿态,对于线阵CCD上的探元S,将相机内的成像指向分解为沿轨方向指向角和垂轨方向指向角,计算成像像点指向角根据设置的成像时刻增量计算沿轨指向角随成像时刻的增量关系;根据成像像点指向角计算T0时刻对应的成像像点的理论沿轨方向指向角根据沿轨指向角随成像时刻的增量关系以及更新成像时刻T0;
S5:重复执行步骤S4,直到满足预设条件,将此时更新后的时刻作为真实成像时刻;
S6:根据真实成像时刻内插影响行,根据T0时刻对应的成像像点的理论沿轨方向指向角和垂轨方向指向角内插影像列,进行反变换计算。
在一种实施方式中,S1包括建立几何定位模型正变换f,如下式所示:
其中,(XS YS ZS)T表示成像时卫星GPS测量的WGS84坐标系下的位置矢量,表示相机坐标系与卫星本体坐标系的转换矩阵,表示卫星本体坐标系与J2000坐标系的转换矩阵,表示J2000坐标系与WGS84坐标系的转换矩阵,m为卫星的成像比例,(ψx,ψy)为探元指向角,是相机内方位元素的综合表示,A=6378137.0m,B=6356752.3m,h为卫星成像时高度,(X Y Z)T为地面坐标,(x y)T为影像坐标。
在一种实施方式中,每片CCD像元数为Wid,影像总行数为Hei,计算各片CCD中心视轴的地表轨迹,S2包括:
在一种实施方式中,CCD中心视轴的地表轨迹与CCD方向之间关系为垂直关系,S3包括:
其中,式3表示矢量相乘,式4表示地物点的近似成像时刻的计算公式。
在一种实施方式中,S4包括:
本申请实施例中的上述一个或多个技术方案,至少具有如下一种或多种技术效果:
首先,利用卫星下传的姿态、轨道数据、成像时间数据和相机参数,建立几何定位模型正变换;然后基于建立的几何定位模型的正变换计算CCD中心视轴的地表轨迹,其中,CCD中心视轴的地表轨迹与时间t具有函数关系;接着对于给定的地物点坐标,根据CCD中心视轴的地表轨迹与CCD方向之间关系、CCD中心视轴的地表轨迹与时间t具有的函数关系,获取地物点的近似成像时刻T0;然后计算沿轨指向角随成像时刻的增量,根据ψy计算对应的理论更新步骤3中的成像时间,最后根据真实成像时刻内插影响行,根据T0时刻对应的成像像点的理论沿轨方向指向角和垂轨方向指向角内插影像列,进行反变换计算。
由于本发明提供的快速反变换计算方法针对TDICCD推扫影像的延迟积分成像特征,根据TDICCD正常成像条件下其视场中心视轴在地面的轨迹与CCD方向应之间的关系以及计算出的函数关系,推算出近似成像时间,然后通过纠正镜头畸变和补偿偏流角残差,得到真实成像时刻,完成反算,不需要进行反复迭代计算,可以解决现有技术中心计算复杂且耗时的技术问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中TDICCD工作原理图;
图2为本发明的具体实施方案步骤S2中基于垂直关系的成像行求解原理图,其中,O为地面点,P(t)为视轴的地表轨迹,通过O点做P(t)的垂线,垂足为S。
图3为本发明具体实施方案步骤S4中由时间误差引起的沿轨指向角偏差示意图,其中T0为成像时卫星所在的真实位置,T1为时间误差导致卫星所在的位置,ψx为T0时刻沿轨向指向角,ψy为T0时刻垂轨向指向角,ψx,为T1时刻沿轨向指向角,ψy,为T1时刻垂轨向指向角。则探元S的沿轨指向角会因成像时刻错误而造成偏差(ψx-ψ'x)。
图4为本发明具体实施方案中步骤1线阵推扫CCD排列特征示意图。
具体实施方式
本发明提出了一种基于TDICCD推扫特性的快速反算算法,是基于严格几何成像模型,从TDICCD推扫成像的本质特征出发,提出“近似垂直约束”的快速反投影算法,可以解决传统反变换算法计算复杂耗时,且不收敛的情况。为了实现上述目的,本发明的主要构思如下:
针对TDICCD推扫影像的延迟积分成像特征,提出了基于视场中心视轴地表轨迹与CCD阵列“近似垂直”约束的反变换算法。TDICCD正常成像条件下其视场中心视轴在地面的轨迹与CCD方向应该近似垂直的工作原理,推算出近似成像时间t,然后通过纠正镜头畸变和补偿偏流角残差,得到成像时刻,完成反算。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提供了一种基于TDICCD推扫特性的快速反算计算方法,该方法包括:
S1:利用卫星下传的姿态、轨道数据、成像时间数据和相机参数,建立几何定位模型正变换;
S2:基于建立的几何定位模型的正变换计算CCD中心视轴的地表轨迹,其中,CCD中心视轴的地表轨迹与时间t具有函数关系;
S3:对于给定的地物点坐标,根据CCD中心视轴的地表轨迹与CCD方向之间关系、CCD中心视轴的地表轨迹与时间t具有的函数关系,获取地物点的近似成像时刻T0;
S4:根据S3中计算的成像时刻T0内插卫星成像位置、姿态,对于线阵CCD上的探元S,将相机内的成像指向分解为沿轨方向指向角和垂轨方向指向角,计算成像像点指向角根据设置的成像时刻增量计算沿轨指向角随成像时刻的增量关系;根据成像像点指向角计算T0时刻对应的成像像点的理论沿轨方向指向角根据沿轨指向角随成像时刻的增量关系以及更新成像时刻T0;
S5:重复执行步骤S4,直到满足预设条件,将此时更新后的时刻作为真实成像时刻;
S6:根据真实成像时刻内插影响行,根据T0时刻对应的成像像点的理论沿轨方向指向角和垂轨方向指向角内插影像列,进行反变换计算。
在一种实施方式中,S1包括建立几何定位模型正变换f,如下式所示:
其中,(XS YS ZS)T表示成像时卫星GPS测量的WGS84坐标系下的位置矢量,表示相机坐标系与卫星本体坐标系的转换矩阵,表示卫星本体坐标系与J2000坐标系的转换矩阵,表示J2000坐标系与WGS84坐标系的转换矩阵,m为卫星的成像比例,(ψx,ψy)为探元指向角,是相机内方位元素的综合表示,A=6378137.0m,B=6356752.3m,h为卫星成像时高度,(X Y Z)T为地面坐标,(x y)T为影像坐标。
在一种实施方式中,每片CCD像元数为Wid,影像总行数为Hei,计算各片CCD中心视轴的地表轨迹,S2包括:
具体地,当地面点O坐标已知的情况下,通过地面点O作视轴轨迹P(t)的垂线,垂足S对应的时刻即为地物点O的近似成像时间,通过上述步骤可以得出视轴轨迹函数P(t)以及对应的系数。
在一种实施方式中,CCD中心视轴的地表轨迹与CCD方向之间关系为垂直关系,S3包括:
其中,式3表示矢量相乘,式4表示地物点的近似成像时刻的计算公式。
具体来说,根据视场中心视轴在地面的轨迹与CCD方向垂直的工作原理,利用矢量相乘,可以计算出第五点的近似成像时刻,近似成像时刻是指估算出的与真实时刻相近的时刻。
在一种实施方式中,S4包括:
总体来说,本发明相对于现有技术来说,具有如下优点或者有益技术效果,
1、本发明公开一种基于TDICCD推扫特性的快速反算计算方法,基于TDICCD推扫影像的延迟积分成像特性,提出了基于视场中心视轴地表轨迹与CCD阵列“近似垂直”约束的反变换算法,解决了传统方法反变换算法的复杂耗时和结果不收敛的情况。
2、本发明公开的一种基于TDICCD推扫特性的快速反算计算方法,提出了基于视场中心视轴地表轨迹与CCD阵列“近似垂直”约束的反变换算法,实现反变换的精度高,相较现有反变换算法具有跟高效率。
3、本发明公开的一种基于TDICCD推扫特性的快速反算计算方法,不仅适用于稳态推扫成像卫星,也可推广应用于未来超敏捷卫星的曲线成像。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于TDICCD推扫特性的快速反变换计算方法,其特征在于,包括:
S1:利用卫星下传的姿态、轨道数据、成像时间数据和相机参数,建立几何定位模型正变换;
S2:基于建立的几何定位模型的正变换计算CCD中心视轴的地表轨迹,其中,CCD中心视轴的地表轨迹与时间t具有函数关系;
S3:对于给定的地物点坐标,根据CCD中心视轴的地表轨迹与CCD方向之间关系、CCD中心视轴的地表轨迹与时间t具有的函数关系,获取地物点的近似成像时刻T0;
S4:根据S3中计算的成像时刻T0内插卫星成像位置、姿态,对于线阵CCD上的探元S,将相机内的成像指向分解为沿轨方向指向角和垂轨方向指向角,计算成像像点指向角根据设置的成像时刻增量计算沿轨指向角随成像时刻的增量关系;根据成像像点指向角计算T0时刻对应的成像像点的理论沿轨方向指向角根据沿轨指向角随成像时刻的增量关系以及更新成像时刻T0;
S5:重复执行步骤S4,直到满足预设条件,将此时更新后的时刻作为真实成像时刻;
S6:根据真实成像时刻内插影响行,根据T0时刻对应的成像像点的理论沿轨方向指向角和垂轨方向指向角内插影像列,进行反变换计算。
5.如权利要求4所述的方法,其特征在于,S4包括:
式中,为内插的卫星成像位置和姿态,(XS YS ZS)表示成像时卫星GPS测量的WGS84坐标系下的位置坐标,表示相机坐标系与卫星本体坐标系的转换矩阵,表示卫星本体坐标系与J2000坐标系的转换矩阵,表示J2000坐标系与WGS84坐标系的转换矩阵;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010197527.3A CN111324857B (zh) | 2020-03-19 | 2020-03-19 | 一种基于tdiccd推扫特性的快速反变换计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010197527.3A CN111324857B (zh) | 2020-03-19 | 2020-03-19 | 一种基于tdiccd推扫特性的快速反变换计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111324857A CN111324857A (zh) | 2020-06-23 |
CN111324857B true CN111324857B (zh) | 2022-03-04 |
Family
ID=71173462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010197527.3A Active CN111324857B (zh) | 2020-03-19 | 2020-03-19 | 一种基于tdiccd推扫特性的快速反变换计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111324857B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112184786B (zh) * | 2020-10-16 | 2024-03-26 | 西安应用光学研究所 | 一种基于合成视觉的目标定位方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226059A (zh) * | 2008-02-03 | 2008-07-23 | 航天东方红卫星有限公司 | 一种星载tdiccd相机积分时间计算及调整方法 |
CN101742050A (zh) * | 2009-12-03 | 2010-06-16 | 浙江大学 | 针对运动模糊核空间移变的tdiccd图像复原方法 |
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
CN105444781A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载自主引导成像地面验证方法 |
CN105528500A (zh) * | 2016-01-19 | 2016-04-27 | 国家测绘地理信息局卫星测绘应用中心 | 一种分米级星载tdi ccd立体测绘相机成像仿真方法和*** |
CN106525002A (zh) * | 2016-09-28 | 2017-03-22 | 北京空间机电研究所 | 一种tdiccd像移检测及补偿方法 |
CN106895851A (zh) * | 2016-12-21 | 2017-06-27 | 中国资源卫星应用中心 | 一种光学遥感卫星多ccd多相机统一处理的传感器校正方法 |
CN108344396A (zh) * | 2018-01-24 | 2018-07-31 | 浙江大学 | 一种敏捷卫星斜条带成像模式姿态计算方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103312994B (zh) * | 2013-05-17 | 2016-04-27 | 中国科学院长春光学精密机械与物理研究所 | 实现面阵cmos传感器双向扫描清晰成像的方法 |
CN103278140B (zh) * | 2013-05-27 | 2015-07-15 | 武汉大学 | 用于tdiccd线阵推扫传感器的坐标反算方法 |
CN103400345A (zh) * | 2013-07-18 | 2013-11-20 | 西南交通大学 | 一种降低tdi-ccd相机图像模糊度的方法 |
CN103986869A (zh) * | 2014-05-22 | 2014-08-13 | 中国科学院长春光学精密机械与物理研究所 | 一种高速tdiccd遥感相机图像采集与显示装置 |
CN104735321B (zh) * | 2015-03-25 | 2017-08-29 | 北京空间机电研究所 | 一种空间tdiccd遥感器光学双向扫描*** |
CN107966137B (zh) * | 2017-11-22 | 2019-05-31 | 安徽农业大学 | 一种基于tdiccd拼接区图像的卫星平台颤振探测方法 |
CN108898565B (zh) * | 2018-07-10 | 2021-09-24 | 中国科学院长春光学精密机械与物理研究所 | Tdi ccd相机摆扫成像影像几何畸变还原的反变换方法 |
-
2020
- 2020-03-19 CN CN202010197527.3A patent/CN111324857B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226059A (zh) * | 2008-02-03 | 2008-07-23 | 航天东方红卫星有限公司 | 一种星载tdiccd相机积分时间计算及调整方法 |
CN101742050A (zh) * | 2009-12-03 | 2010-06-16 | 浙江大学 | 针对运动模糊核空间移变的tdiccd图像复原方法 |
CN102229362A (zh) * | 2011-04-29 | 2011-11-02 | 航天东方红卫星有限公司 | 敏捷卫星成像反向推扫速度失配时的姿态补偿方法 |
CN105444781A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载自主引导成像地面验证方法 |
CN105528500A (zh) * | 2016-01-19 | 2016-04-27 | 国家测绘地理信息局卫星测绘应用中心 | 一种分米级星载tdi ccd立体测绘相机成像仿真方法和*** |
CN106525002A (zh) * | 2016-09-28 | 2017-03-22 | 北京空间机电研究所 | 一种tdiccd像移检测及补偿方法 |
CN106895851A (zh) * | 2016-12-21 | 2017-06-27 | 中国资源卫星应用中心 | 一种光学遥感卫星多ccd多相机统一处理的传感器校正方法 |
CN108344396A (zh) * | 2018-01-24 | 2018-07-31 | 浙江大学 | 一种敏捷卫星斜条带成像模式姿态计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111324857A (zh) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109903352B (zh) | 一种卫星遥感影像大区域无缝正射影像制作方法 | |
CN108830889B (zh) | 基于全局几何约束的遥感影像与基准影像的匹配方法 | |
CN103697864B (zh) | 一种基于大虚拟相机的窄视场双相机影像拼接方法 | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN107040695B (zh) | 基于rpc定位模型的星载视频稳像方法及*** | |
CN110006452B (zh) | 高分六号宽视场相机相对几何定标方法及*** | |
CN113793270A (zh) | 一种基于无人机姿态信息的航拍图像几何校正方法 | |
CN111324857B (zh) | 一种基于tdiccd推扫特性的快速反变换计算方法 | |
CN108447100B (zh) | 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法 | |
CN111538051B (zh) | 一种摆扫大幅宽光学卫星精准处理方法 | |
CN114972078A (zh) | 应用sar影像提升国产光学卫星影像无控几何质量方法及*** | |
CN110853140A (zh) | 一种dem辅助的光学视频卫星影像稳像方法 | |
CN107516291B (zh) | 夜景影像正射纠正处理方法 | |
CN111275773A (zh) | 无场几何定标方法及*** | |
CN114838739B (zh) | 一种考虑完整回归周期的卫星图像几何定标方法 | |
CN109579796B (zh) | 一种投影后影像的区域网平差方法 | |
CN114485574B (zh) | 基于卡尔曼滤波模型的三线阵影像pos辅助对地定位方法 | |
CN111161186A (zh) | 一种推扫式遥感器通道配准方法及装置 | |
CN113379648B (zh) | 一种高分七号和资源三号立体影像联合平差方法 | |
CN109493298B (zh) | 一种机载摆扫式高光谱数据快速几何校正方法 | |
CN111862332A (zh) | 一种卫星影像通用成像模型拟合误差的改正方法及*** | |
CN115164938B (zh) | 一种线阵摆扫式光学遥感卫星在轨几何定标方法 | |
CN117934346B (zh) | 一种无稳定平台的机载线阵高光谱遥感影像几何处理方法 | |
CN109741381B (zh) | 基于平行观测的星载推扫式光学传感器高频误差消除方法 | |
CN117008101A (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 |