CN111598810A - 旋转光学卫星影像退化信息感知与复原方法及系统 - Google Patents
旋转光学卫星影像退化信息感知与复原方法及系统 Download PDFInfo
- Publication number
- CN111598810A CN111598810A CN202010424134.1A CN202010424134A CN111598810A CN 111598810 A CN111598810 A CN 111598810A CN 202010424134 A CN202010424134 A CN 202010424134A CN 111598810 A CN111598810 A CN 111598810A
- Authority
- CN
- China
- Prior art keywords
- image
- degradation information
- frequency
- optical satellite
- information
- 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
- 230000003287 optical effect Effects 0.000 title claims abstract description 159
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000006731 degradation reaction Methods 0.000 claims abstract description 257
- 230000015556 catabolic process Effects 0.000 claims abstract description 250
- 238000003384 imaging method Methods 0.000 claims abstract description 36
- 230000033001 locomotion Effects 0.000 claims description 70
- 238000001914 filtration Methods 0.000 claims description 44
- 238000005259 measurement Methods 0.000 claims description 38
- 206010044565 Tremor Diseases 0.000 claims description 34
- 230000007246 mechanism Effects 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000000926 separation method Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 12
- 238000007500 overflow downdraw method Methods 0.000 claims description 9
- 230000002457 bidirectional effect Effects 0.000 claims description 8
- 230000004927 fusion Effects 0.000 claims description 8
- 238000010183 spectrum analysis Methods 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 claims description 5
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 239000000835 fiber Substances 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 206010020751 Hypersensitivity Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本发明提供了一种旋转光学卫星影像退化信息感知与复原方法及系统,所述旋转光学卫星影像退化信息感知与复原方法包括:旋转光学卫星进行旋转成像,形成旋转光学卫星影像,计算旋转光学卫星影像的退化信息;基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;对所述高频退化信息进行时空建模,形成高频退化信息时空模型;对所述低频退化信息进行时空建模,形成低频退化信息时空模型;基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
Description
技术领域
本发明涉及光学卫星技术领域,特别涉及一种旋转光学卫星影像退化信息感知与复原方法及系统。
背景技术
超高、超敏、超稳、超宽幅光学卫星是当前我国高分辨率对地观测系统的重点发展对象,动中成像是其重要成像模式。面对日益增长的区域大范围高时效对地观测需求,旋转对地成像逐渐成为一种新型模式并逐步得到研究与发展。旋转超大宽幅光学卫星采用面阵相机圆锥扫描对地观测技术体制,通过圆锥扫描方式实现超大区域观测,设计广域搜索成像与重点地区加密成像模式,具备广域、高分、多角度成像等特点。该卫星由平台和有效载荷组成,有效载荷包括:面阵光学相机分系统、载荷数据处理分系统、辅助载荷分系统和转台分系统。卫星平台包括:结构与机构分系统、热控分系统、电源与总体电路分系统、测控分系统、姿轨控分系统、星务分系统、数传与图像分发分系统。卫星采取天底指向,即转台旋转轴对地,转台保持正常工作转速,每旋转一周相机摆镜调整一档,摆镜共3个位置档位,旋转3圈后可覆盖天底方向5~60°范围,可实现千级公里直径区域覆盖。机动大范围扫描成像提高了影像数据获取效率,同时也带来了影像质量退化与畸变,如何解决旋转动中成像影像退化信息获取与复原,最终获得连续、高质的标准产品是本发明亟需解决的关键技术难题与重要瓶颈。
发明内容
本发明的目的在于提供一种旋转光学卫星影像退化信息感知与复原方法及系统,以解决现有的旋转大幅宽光学卫星成像机理复杂且影像质量退化模糊的问题。
为解决上述技术问题,本发明提供一种旋转光学卫星影像退化信息感知与复原方法,所述旋转光学卫星影像退化信息感知与复原方法包括:
旋转光学卫星进行旋转成像,形成旋转光学卫星影像,计算旋转光学卫星影像的退化信息;
基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;
对所述高频退化信息进行时空建模,形成高频退化信息时空模型;
对所述低频退化信息进行时空建模,形成低频退化信息时空模型;
基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,所述旋转光学卫星安装有光学相机面阵分系统和转台分系统,其中:
所述光学相机面阵分系统包括转动补偿机构,其中:
所述转动补偿机构包括摆镜档位机构、快速反射补偿镜及探测器转动机构;
所述旋转光学卫星进行旋转成像包括:所述旋转超大宽幅光学卫星采取天底指向并降低对地速度,所述转台分系统旋转轴对地,所述转台分系统保持正常工作转速旋转;
所述转台分系统每旋转一周,所述摆镜档位机构的位置挡位进行一次调整,以使光学相机面阵分系统完成对地观测的圆锥扫描。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,计算旋转光学卫星影像的退化信息包括:
所述旋转光学卫星采用快速反射补偿镜及探测器转动机构对平动像移与转动像移进行实时补偿,并采用绝对空间指向测量子系统以及相对空间指向测量子系统进行降质运动参数信息测量;
所述旋转光学卫星采用双向滤波整体加权融合方法执行绝对空间指向测量子系统以及相对空间指向测量子系统的高低频数据融合,计算旋转光学卫星影像的退化信息。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,计算旋转光学卫星影像的退化信息还包括:
对所述绝对空间指向测量子系统的观测数据进行融合,得到第一绝对指向;包括:提取光轴矢量、构建误差方程、计算绝对指向;
对所述相对空间指向测量子系统内部的摆镜档位机构、快速反射补偿镜及探测器转动机构、转台分系统、低频机械陀螺、光纤陀螺以及高频激光陀螺按照测量频率从小到大进行排序,得到第一相对指向、第二相对指向、第三相对指向,并以此类推;
将第一绝对指向与第一相对指向进行双向滤波整体加权融合方法得到第二绝对指向;将第二绝对指向与第二相对指向进行双向滤波整体加权融合方法得到第三绝对指向,并依次类推完成与每个相对指向测量子系统测量值的融合,实现旋转光学卫星影像的退化信息的计算;
双向滤波整体加权融合方法表示如下:
其中,根符号f表示正向滤波结果,根符号b表示反向滤波结果,根符号s表示整体加权平差结果;表示tk时刻正向滤波姿态最优估计值,表示tk时刻反向滤波姿态最优估计值的逆,表示tk时刻正向滤波与反向滤波结果解算出的四元数误差;表示系统状态变量误差矩阵,表示返回标量的正负号,表示tk时刻漂移最优估计值;与分别表示正向滤波与反向滤波误差协方差阵最优估计值的逆,表示整理加权的误差协方差阵估计值;表示系统状态变量误差矩阵的加权平均值,表示误差四元数整体加权平均矢量部分结果,表示误差漂移整体加权平均结果;表示误差四元数整体加权结果,表示误差漂移整体加权结果;表示tk时刻四元数整体加权结果,表示tk时刻漂移量整体加权结果。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,形成高频退化信息和低频退化信息包括:
对精确估计的旋转光学卫星的运动参数信号进行时-频-谱分析,采用快速傅里叶变换将所述旋转光学卫星影像的退化信息从时域信号转为频域信号,采用高通滤波器和低通滤波器对所述旋转光学卫星影像的退化信息进行信号分离,分别形成所述高频退化信息和所述低频退化信息。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,对所述高频退化信息进行时空建模,形成高频退化信息时空模型,对所述低频退化信息进行时空建模,形成低频退化信息时空模型包括:
采用以时间和姿态角共同作为变量的数学模型作为震颤扰动误差时序曲线的拟合模型,如式(2):
f(t)=F(ψ,ω,κ)W(t) (2)
其中,f(t)表示震颤运动特性参数,ψ,ω,κ分别表示卫星姿态的俯仰角、滚动角角和偏航角,F(ψ,ω,κ)表示姿态角为自变量的震颤扰动空间特性函数,W(t)表示以正弦函数作为核函数的震颤扰动时序特性函数。
可选的,在所述的旋转光学卫星影像退化信息感知与复原方法中,基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像包括:
旋转光学卫星成像过程中,通过快速反射补偿镜与探测器转动机构进行实时像移补偿,单帧影像降质退化信息建模为一个退化时空函数h(x,y)与一个加性噪声n(x,y),得到:
g(x,y)=h(x,y)*f(x,y)+n(x,y) (3)
式中,f(x,y)表示退化前影像,g(x,y)表示退化后影像,*表示卷积;对公式(3)进行傅里叶变换,得到频域中表示:
G(u,v)=H(u,v)*F(u,v)+N(u,v) (4)
式中,G(u,v),H(u,v),F(u,v),N(u,v)表示傅里叶频谱;根据得到的退化信息时空模型估计点扩散函数;将复杂运动过程简化成匀速直线运动,影像在x和y方向移动的距离分别为a和b,进一步得到降质传递函数:
式中,x(t)、y(t)表示运动参数信息,根据退化信息时空模型计算得到;确定退化时空函数后,采用逆滤波法、维纳滤波法以及约束最小二乘方滤波等方法实现单帧降质影像补偿复原,得到高质量影像产品。
本发明还提供一种旋转光学卫星影像退化信息感知与复原系统,所述旋转光学卫星影像退化信息感知与复原系统包括退化信息计算模块、退化信息分离模块、时空建模模块及单面阵影像补偿模块,其中:
旋转光学卫星进行旋转成像,形成旋转光学卫星影像,所述退化信息计算模块计算旋转光学卫星影像的退化信息;
所述退化信息分离模块基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;
所述时空建模模块对所述高频退化信息进行时空建模,形成高频退化信息时空模型;
所述时空建模模块对所述低频退化信息进行时空建模,形成低频退化信息时空模型;
所述单面阵影像补偿模块基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
在本发明提供的旋转光学卫星影像退化信息感知与复原方法及系统中,通过退化信息计算模块计算旋转光学卫星影像的退化信息,退化信息分离模块基于傅里叶模型检测旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息,时空建模模块对高频退化信息和低频退化信息进行时空建模,形成高低频退化信息时空模型,单面阵影像补偿模块基于高频退化信息时空模型和低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像,实现了旋转动中成像过程信息精确感知与建模,提升影像成像质量,通过本发明提供的方法可以解决旋转大幅宽卫星成像机理复杂且降质影像难以复原的关键问题,为高分辨率光学卫星高效、高质量后续应用奠定基础。
本发明提供了一种旋转光学卫星影像退化信息感知与复原方法,旋转宽幅光学卫星采用圆锥扫描结合面阵的模式进行对地成像获得大区域高时效覆盖,整个成像过程涉及的载荷与平台单机较多,其运动特性主要由相机载荷旋转、地速(卫星轨道运动加地球自转)、震颤扰动等复合而成。针对影像存在的模糊与像移,本发明设计采用快反镜摆动结合探测器转动方法对平动像移与转动像移进行实时补偿,并采用空间指向测量系统(震颤测量系统)进行降质运动参数信息测量。根据工作机理与误差特性进行广义分类分为:绝对空间指向测量子系统以及相对指向测量子系统。
进一步的,所述旋转光学卫星运动属于复合运动,包括旋转、地速(卫星轨道运动+地球自转)、扰动(震颤运动)。旋转与地速引起的像移通过快速反射补偿镜+探测器转动机构进行了实时补偿,影像的降质过程主要表现为震颤运动引起,故需要进一步对精确的估计的运动状态参数进行固有运动信号与震颤信号检测与分离。本发明对精确估计的运动参数信号进行时-频-谱分析,采用高通、低通滤波器进行信号分离。由于项目所涉及到的震颤时域信号属于离散的高采样率信号,故拟采用快速傅里叶变换(Fast Fourier Transform,FFT)将时域信号转到频域信号,并进行相应的频谱分析。
附图说明
图1~2是本发明一实施例的旋转光学卫星影像退化信息感知与复原方法示意图。
具体实施方式
以下结合附图和具体实施例对本发明提出的旋转光学卫星影像退化信息感知与复原方法及系统作进一步详细说明。根据下面说明和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
本发明的核心思想在于提供一种旋转光学卫星影像退化信息感知与复原方法及系统,以解决现有的旋转大幅宽光学卫星成像机理复杂且影像质量退化模糊的问题。
为实现上述思想,本发明提供了一种旋转光学卫星影像退化信息感知与复原方法及系统,所述旋转光学卫星影像退化信息感知与复原系统包括退化信息计算模块、退化信息分离模块、时空建模模块及单面阵影像补偿模块,其中:旋转光学卫星进行旋转成像,形成旋转光学卫星影像,所述退化信息计算模块计算旋转光学卫星影像的退化信息;所述退化信息分离模块基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;所述时空建模模块对所述高频退化信息进行时空建模,形成高频退化信息时空模型;所述时空建模模块对所述低频退化信息进行时空建模,形成低频退化信息时空模型;所述单面阵影像补偿模块基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
<实施例一>
本实施例提供一种旋转光学卫星影像退化信息感知与复原方法,如图1所示,所述旋转光学卫星影像退化信息感知与复原方法包括:旋转光学卫星进行旋转成像,形成旋转光学卫星影像,计算旋转光学卫星影像的退化信息;基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;对所述高频退化信息进行时空建模,形成高频退化信息时空模型;对所述低频退化信息进行时空建模,形成低频退化信息时空模型;基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
图2表示的是旋转大幅宽光学卫星影像退化信息感知与复原流程图,以下针对实施例流程中的各步骤,对本发明方法做进一步详细描述。
步骤1,旋转成像退化信息精确计算。
旋转宽幅光学卫星采用圆锥扫描+面阵的模式进行对地成像获得大区域高时效覆盖,整个成像过程涉及的载荷与平台单机较多,其运动特性主要由相机载荷旋转、地速(卫星轨道运动+地球自转)、震颤扰动等复合而成。针对影像存在的模糊与像移,卫星设计采用快反镜摆动+探测器转动方法对平动像移与转动像移进行实时补偿,并采用空间指向测量系统(震颤测量系统)进行降质运动参数信息测量。根据工作机理与误差特性进行广义分类分为:绝对空间指向测量子系统以及相对指向测量子系统。
本发明采用双向滤波整体加权的渐进式方法实现绝对与相对指向测量子系统高低频数据高精度精细化融合,实现旋转成像退化信息精确计算。拟采用的具体方法如下:首先对绝对指向测量子系统观测数据进行融合(光轴矢量提取、误差方程构建、绝对指向计算)得到绝对指向1,进一步对相对指向测量子系统内部的摆镜机构、快速反射补偿镜、探测器转动机构、转台、低频机械陀螺、光纤陀螺以及高频激光陀螺按照测量频率从小到大排序编号:相对指向1、相对指向2、相对指向3……,依次将绝对指向1与相对指向1进行双向滤波整体加权融合得到绝对指向2,将绝对指向2与相对指向2进行融合结果得到绝对指向3,依次类推完成与每个相对指向测量子系统测量值的融合,实现退化过程信息精确计算。双向滤波整体加权融合方法表示如下:
根符号f表示正向滤波结果,根符号b表示反向滤波结果,根符号s表示整体加权平差结果;表示tk时刻正向滤波姿态最优估计值,表示tk时刻反向滤波姿态最优估计值的逆,表示tk时刻正向滤波与反向滤波结果解算出的四元数误差;表示系统状态变量误差矩阵,表示返回标量的正负号,表示tk时刻漂移最优估计值;与分别表示正向滤波与反向滤波误差协方差阵最优估计值的逆,表示整理加权的误差协方差阵估计值;表示系统状态变量误差矩阵的加权平均值,表示误差四元数整体加权平均矢量部分结果,表示误差漂移整体加权平均结果;表示误差四元数整体加权结果,表示误差漂移整体加权结果;表示tk时刻四元数整体加权结果,表示tk时刻漂移量整体加权结果。
步骤2,基于傅立叶模型的高低频退化信息检测与分离。
旋转大宽幅光学卫星运动属于复合运动,包括旋转、地速(卫星轨道运动+地球自转)、扰动(震颤运动)。旋转与地速引起的像移通过快速反射补偿镜+探测器转动机构进行了实时补偿,影像的降质过程主要表现为震颤运动引起,故需要进一步对精确的估计的运动状态参数进行固有运动信号与震颤信号检测与分离。本发明对精确估计的运动参数信号进行时-频-谱分析,采用高通、低通滤波器进行信号分离。由于项目所涉及到的震颤时域信号属于离散的高采样率信号,故拟采用快速傅里叶变换(Fast Fourier Transform,FFT)将时域信号转到频域信号,并进行相应的频谱分析。
步骤3,高低频退化信息时空建模。
为实现后续降质退化模型构建以及补偿复原,须进一步对高低频退化运动参数信息进行时空建模。由于旋转动中成像时姿态角不断变化,震颤扰动对成像焦面的降质影响也随之变化,导致震颤扰动误差时空规律复杂,除了在时间维度上变化,在空间维度也在随卫星成像状态发生变化,采用传统仅考虑时序特性的单一模型难以准确描述震颤扰动误差变化规律。为此,本发明采用以时间和姿态角共同作为变量的数学模型作为震颤扰动误差时序曲线的拟合模型,如下式:
f(t)=F(ψ,ω,κ)W(t) (2)
其中,f(t)表示震颤运动特性参数,ψ,ω,κ分别表示卫星姿态的俯仰角、滚动角角和偏航角,F(ψ,ω,κ)表示姿态角为自变量的震颤扰动空间特性函数,W(t)表示以正弦函数作为核函数的震颤扰动时序特性函数。
步骤4,基于退化时空模型的单面阵影像补偿复原。
数字图像复原技术是数字图像处理技术的一个重要分支,利用退化现象的先验知识重建或恢复一幅退化图像。旋转宽幅光学卫星成像过程中,通过快速反射补偿镜+探测器转动机构进行了实时像移补偿,单帧影像降质机理主要由剩余扰动所造成,其退化过程被建模为一个退化时空函数h(x,y)与一个加性噪声n(x,y),得到:
g(x,y)=h(x,y)*f(x,y)+n(x,y) (3)
式中,f(x,y)表示退化前影像,g(x,y)表示退化后影像,*表示卷积。对上述模型进行傅里叶变换,得到频域中表示:
G(u,v)=H(u,v)*F(u,v)+N(u,v) (4)
式中,G(u,v),H(u,v),F(u,v),N(u,v)表示傅里叶频谱。影像补偿复原属于影像降质退化的逆过程,进一步根据研究内容二得到的退化信息时空模型估计点扩散函数(运动模糊核)。由于曝光时间T很短,将复杂运动过程简化成匀速直线运动,影像在x和y方向移动的距离分别为a和b,进一步得到降质传递函数:
式中,x(t)、y(t)表示运动参数信息,根据退化信息时空模型计算得到。确定退化时空函数后,进一步拟采用逆滤波法、维纳滤波法以及约束最小二乘方滤波等方法实现单帧降质影像补偿复原,得到高质量影像产品。
综上,上述实施例对旋转光学卫星影像退化信息感知与复原方法及系统的不同构型进行了详细说明,当然,本发明包括但不局限于上述实施中所列举的构型,任何在上述实施例提供的构型基础上进行变换的内容,均属于本发明所保护的范围。本领域技术人员可以根据上述实施例的内容举一反三。
<实施例二>
本实施例提供一种旋转光学卫星影像退化信息感知与复原系统,所述旋转光学卫星影像退化信息感知与复原系统包括退化信息计算模块、退化信息分离模块、时空建模模块及单面阵影像补偿模块,其中:旋转光学卫星进行旋转成像,形成旋转光学卫星影像,所述退化信息计算模块计算旋转光学卫星影像的退化信息;所述退化信息分离模块基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;所述时空建模模块对所述高频退化信息进行时空建模,形成高频退化信息时空模型;所述时空建模模块对所述低频退化信息进行时空建模,形成低频退化信息时空模型;所述单面阵影像补偿模块基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
在本发明提供的旋转光学卫星影像退化信息感知与复原方法及系统中,通过退化信息计算模块计算旋转光学卫星影像的退化信息,退化信息分离模块基于傅里叶模型检测旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息,时空建模模块对高频退化信息和低频退化信息进行时空建模,形成高低频退化信息时空模型,单面阵影像补偿模块基于高频退化信息时空模型和低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像,实现了旋转动中成像过程信息精确感知与建模,提升影像成像质量,通过本发明提供的方法可以解决旋转大幅宽卫星成像机理复杂且降质影像难以复原的关键问题,为高分辨率光学卫星高效、高质量后续应用奠定基础。
本发明提供了一种旋转光学卫星影像退化信息感知与复原方法,旋转宽幅光学卫星采用圆锥扫描结合面阵的模式进行对地成像获得大区域高时效覆盖,整个成像过程涉及的载荷与平台单机较多,其运动特性主要由相机载荷旋转、地速(卫星轨道运动加地球自转)、震颤扰动等复合而成。针对影像存在的模糊与像移,本发明设计采用快反镜摆动结合探测器转动方法对平动像移与转动像移进行实时补偿,并采用空间指向测量系统(震颤测量系统)进行降质运动参数信息测量。根据工作机理与误差特性进行广义分类分为:绝对空间指向测量子系统以及相对指向测量子系统。
进一步的,所述旋转光学卫星运动属于复合运动,包括旋转、地速(卫星轨道运动+地球自转)、扰动(震颤运动)。旋转与地速引起的像移通过快速反射补偿镜+探测器转动机构进行了实时补偿,影像的降质过程主要表现为震颤运动引起,故需要进一步对精确的估计的运动状态参数进行固有运动信号与震颤信号检测与分离。本发明对精确估计的运动参数信号进行时-频-谱分析,采用高通、低通滤波器进行信号分离。由于项目所涉及到的震颤时域信号属于离散的高采样率信号,故拟采用快速傅里叶变换(Fast Fourier Transform,FFT)将时域信号转到频域信号,并进行相应的频谱分析。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
上述描述仅是对本发明较佳实施例的描述,并非对本发明范围的任何限定,本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰,均属于权利要求书的保护范围。
Claims (8)
1.一种旋转光学卫星影像退化信息感知与复原方法,其特征在于,所述旋转光学卫星影像退化信息感知与复原方法包括:
旋转光学卫星进行旋转成像,形成旋转光学卫星影像,计算旋转光学卫星影像的退化信息;
基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;
对所述高频退化信息进行时空建模,形成高频退化信息时空模型;
对所述低频退化信息进行时空建模,形成低频退化信息时空模型;
基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
2.如权利要求1所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,所述旋转光学卫星安装有光学相机面阵分系统和转台分系统,其中:
所述光学相机面阵分系统包括转动补偿机构,其中:
所述转动补偿机构包括摆镜档位机构、快速反射补偿镜及探测器转动机构;
所述旋转光学卫星进行旋转成像包括:所述旋转超大宽幅光学卫星采取天底指向并降低对地速度,所述转台分系统旋转轴对地,所述转台分系统保持正常工作转速旋转;
所述转台分系统每旋转一周,所述摆镜档位机构的位置挡位进行一次调整,以使光学相机面阵分系统完成对地观测的圆锥扫描。
3.如权利要求2所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,计算旋转光学卫星影像的退化信息包括:
所述旋转光学卫星采用快速反射补偿镜及探测器转动机构对平动像移与转动像移进行实时补偿,并采用绝对空间指向测量子系统以及相对空间指向测量子系统进行降质运动参数信息测量;
所述旋转光学卫星采用双向滤波整体加权融合方法执行绝对空间指向测量子系统以及相对空间指向测量子系统的高低频数据融合,计算旋转光学卫星影像的退化信息。
4.如权利要求3所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,计算旋转光学卫星影像的退化信息还包括:
对所述绝对空间指向测量子系统的观测数据进行融合,得到第一绝对指向;包括:提取光轴矢量、构建误差方程、计算绝对指向;
对所述相对空间指向测量子系统内部的摆镜档位机构、快速反射补偿镜及探测器转动机构、转台分系统、低频机械陀螺、光纤陀螺以及高频激光陀螺按照测量频率从小到大进行排序,得到第一相对指向、第二相对指向、第三相对指向,并以此类推;
将第一绝对指向与第一相对指向进行双向滤波整体加权融合方法得到第二绝对指向;将第二绝对指向与第二相对指向进行双向滤波整体加权融合方法得到第三绝对指向,并依次类推完成与每个相对指向测量子系统测量值的融合,实现旋转光学卫星影像的退化信息的计算;
双向滤波整体加权融合方法表示如下:
其中,根符号f表示正向滤波结果,根符号b表示反向滤波结果,根符号s表示整体加权平差结果;表示tk时刻正向滤波姿态最优估计值,表示tk时刻反向滤波姿态最优估计值的逆,表示tk时刻正向滤波与反向滤波结果解算出的四元数误差;表示系统状态变量误差矩阵,表示返回标量的正负号,表示tk时刻漂移最优估计值;与分别表示正向滤波与反向滤波误差协方差阵最优估计值的逆,表示整理加权的误差协方差阵估计值;表示系统状态变量误差矩阵的加权平均值,表示误差四元数整体加权平均矢量部分结果,表示误差漂移整体加权平均结果;表示误差四元数整体加权结果,表示误差漂移整体加权结果;表示tk时刻四元数整体加权结果,表示tk时刻漂移量整体加权结果。
5.如权利要求4所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,形成高频退化信息和低频退化信息包括:
对精确估计的旋转光学卫星的运动参数信号进行时-频-谱分析,采用快速傅里叶变换将所述旋转光学卫星影像的退化信息从时域信号转为频域信号,采用高通滤波器和低通滤波器对所述旋转光学卫星影像的退化信息进行信号分离,分别形成所述高频退化信息和所述低频退化信息。
6.如权利要求5所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,对所述高频退化信息进行时空建模,形成高频退化信息时空模型,对所述低频退化信息进行时空建模,形成低频退化信息时空模型包括:
采用以时间和姿态角共同作为变量的数学模型作为震颤扰动误差时序曲线的拟合模型,如式(2):
f(t)=F(ψ,ω,κ)W(t) (2)
其中,f(t)表示震颤运动特性参数,ψ,ω,κ分别表示卫星姿态的俯仰角、滚动角角和偏航角,F(ψ,ω,κ)表示姿态角为自变量的震颤扰动空间特性函数,W(t)表示以正弦函数作为核函数的震颤扰动时序特性函数。
7.如权利要求6所述的旋转光学卫星影像退化信息感知与复原方法,其特征在于,基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像包括:
旋转光学卫星成像过程中,通过快速反射补偿镜与探测器转动机构进行实时像移补偿,单帧影像降质退化信息建模为一个退化时空函数h(x,y)与一个加性噪声n(x,y),得到:
g(x,y)=h(x,y)*f(x,y)+n(x,y) (3)
式中,f(x,y)表示退化前影像,g(x,y)表示退化后影像,*表示卷积;对公式(3)进行傅里叶变换,得到频域中表示:
G(u,v)=H(u,v)*F(u,v)+N(u,v) (4)
式中,G(u,v),H(u,v),F(u,v),N(u,v)表示傅里叶频谱;根据得到的退化信息时空模型估计点扩散函数;将复杂运动过程简化成匀速直线运动,影像在x和y方向移动的距离分别为a和b,进一步得到降质传递函数:
式中,x(t)、y(t)表示运动参数信息,根据退化信息时空模型计算得到;确定退化时空函数后,采用逆滤波法、维纳滤波法以及约束最小二乘方滤波等方法实现单帧降质影像补偿复原,得到高质量影像产品。
8.一种旋转光学卫星影像退化信息感知与复原系统,其特征在于,所述旋转光学卫星影像退化信息感知与复原系统包括退化信息计算模块、退化信息分离模块、时空建模模块及单面阵影像补偿模块,其中:
旋转光学卫星进行旋转成像,形成旋转光学卫星影像,所述退化信息计算模块计算旋转光学卫星影像的退化信息;
所述退化信息分离模块基于傅里叶模型检测所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息,将所述旋转光学卫星影像的退化信息中的高频退化信息和低频退化信息分离,形成高频退化信息和低频退化信息;
所述时空建模模块对所述高频退化信息进行时空建模,形成高频退化信息时空模型;
所述时空建模模块对所述低频退化信息进行时空建模,形成低频退化信息时空模型;
所述单面阵影像补偿模块基于所述高频退化信息时空模型和所述低频退化信息时空模型,补偿旋转光学卫星的单面阵影像,复原旋转光学卫星影像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010424134.1A CN111598810B (zh) | 2020-05-19 | 2020-05-19 | 旋转光学卫星影像退化信息感知与复原方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010424134.1A CN111598810B (zh) | 2020-05-19 | 2020-05-19 | 旋转光学卫星影像退化信息感知与复原方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111598810A true CN111598810A (zh) | 2020-08-28 |
CN111598810B CN111598810B (zh) | 2023-05-12 |
Family
ID=72183818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010424134.1A Active CN111598810B (zh) | 2020-05-19 | 2020-05-19 | 旋转光学卫星影像退化信息感知与复原方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111598810B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113163138A (zh) * | 2021-05-20 | 2021-07-23 | 苏州大学 | 基于双向循环网络的高分辨率视频复原系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020004691A1 (en) * | 2000-03-10 | 2002-01-10 | Yasuhiro Kinashi | Attitude determination and alignment using electro-optical sensors and global navigation satellites |
CN104820969A (zh) * | 2015-04-03 | 2015-08-05 | 西安交通大学 | 一种实时图像盲复原方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN109668579A (zh) * | 2019-01-23 | 2019-04-23 | 张过 | 基于角位移传感器的星载推扫光学载荷高频误差补偿方法 |
-
2020
- 2020-05-19 CN CN202010424134.1A patent/CN111598810B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020004691A1 (en) * | 2000-03-10 | 2002-01-10 | Yasuhiro Kinashi | Attitude determination and alignment using electro-optical sensors and global navigation satellites |
CN104820969A (zh) * | 2015-04-03 | 2015-08-05 | 西安交通大学 | 一种实时图像盲复原方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN109668579A (zh) * | 2019-01-23 | 2019-04-23 | 张过 | 基于角位移传感器的星载推扫光学载荷高频误差补偿方法 |
Non-Patent Citations (2)
Title |
---|
南一冰等: "卫星振动频率和振幅对高光谱成像质量影响的仿真分析和实验研究", 《光谱学与光谱分析》 * |
周春平等: "遥感图像MTF复原国内研究现状", 《航天返回与遥感》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113163138A (zh) * | 2021-05-20 | 2021-07-23 | 苏州大学 | 基于双向循环网络的高分辨率视频复原系统及方法 |
CN113163138B (zh) * | 2021-05-20 | 2023-01-17 | 苏州大学 | 基于双向循环网络的高分辨率视频复原系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111598810B (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7856154B2 (en) | System and method of super-resolution imaging from a sequence of translated and rotated low-resolution images | |
CN110296691A (zh) | 融合imu标定的双目立体视觉测量方法与系统 | |
EP3043160B1 (en) | Lenslet, beamwalk and tilt diversity for anisoplanatic imaging by large-aperture telescopes | |
Codona et al. | Focal plane wavefront sensing using residual adaptive optics speckles | |
CN111612693B (zh) | 一种旋转大幅宽光学卫星传感器校正方法 | |
CN113029132B (zh) | 一种结合地面影像与恒星光行差测量的航天器导航方法 | |
CN106768327A (zh) | 一种液晶可调谐滤波器成像光谱重建方法 | |
CN112598608A (zh) | 一种基于目标区域的光学卫星快速融合产品制作方法 | |
CN111561935B (zh) | 旋转大幅宽光学卫星在轨几何标定方法及系统 | |
CN111561936A (zh) | 旋转大幅宽光学卫星精准处理方法及系统 | |
CN111524196B (zh) | 一种摆扫大幅宽光学卫星在轨几何标定方法 | |
CN111598810A (zh) | 旋转光学卫星影像退化信息感知与复原方法及系统 | |
CN115752250A (zh) | 一种融合计算机视觉及加速度的桥梁高精度位移监测方法 | |
CN115469308A (zh) | 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质 | |
CN115086550A (zh) | 元成像方法与系统 | |
CN109029379B (zh) | 一种高精度小基高比立体测绘方法 | |
CN106558036B (zh) | 一种星载超分辨成像系统设计方法 | |
Jolissaint et al. | Adaptive optics point spread function reconstruction project at WM Keck Observatory: first results with faint natural guide stars | |
CN114170192A (zh) | 星载光学相机焦面震颤精细检测方法 | |
Westin et al. | Orthorectification of EROS A1 images | |
Fan et al. | Key techniques for virtual imaging model construction and correction of rotating imaging satellite | |
CN117934346B (zh) | 一种无稳定平台的机载线阵高光谱遥感影像几何处理方法 | |
Tao et al. | Image improvement based on sub-pixel post-integration for a staring imaging system in geostationary orbit | |
Lombardi et al. | Accuracy of high resolution CSK interferometric Digital Elevation Models | |
Zhu et al. | An integral time calculation model for agile satellite TDICCD camera |
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 |