CN117475087B - 一种星载tdi ccd影像震颤成像畸变建模方法及装置 - Google Patents
一种星载tdi ccd影像震颤成像畸变建模方法及装置 Download PDFInfo
- Publication number
- CN117475087B CN117475087B CN202311784666.6A CN202311784666A CN117475087B CN 117475087 B CN117475087 B CN 117475087B CN 202311784666 A CN202311784666 A CN 202311784666A CN 117475087 B CN117475087 B CN 117475087B
- Authority
- CN
- China
- Prior art keywords
- imaging
- tremor
- integral
- image
- distortion
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 591
- 206010044565 Tremor Diseases 0.000 title claims abstract description 547
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000010354 integration Effects 0.000 claims description 69
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000004075 alteration Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 abstract description 10
- 238000012937 correction Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 12
- 238000001228 spectrum Methods 0.000 description 9
- 230000033001 locomotion Effects 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000010363 phase shift Effects 0.000 description 5
- 238000003786 synthesis reaction Methods 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 238000012952 Resampling Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
- G06T2207/10036—Multispectral image; Hyperspectral image
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Color Television Image Signal Generators (AREA)
Abstract
本发明公开了一种星载TDI CCD影像震颤成像畸变建模方法及装置，涉及卫星遥感技术领域，包括对多级积分TDI CCD影像震颤成像绝对畸变分析，获取多级积分TDI CCD影像的震颤像移曲线；基于TDI CCD分时成像多光谱影像和震颤像移曲线，构建严密视差成像模型；基于该模型；进行多级积分TDI CCD影像震颤成像绝对畸变建模。本发明实现单级积分多光谱影像的震颤成像绝对畸变检测，全色影像和多光谱影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模，为多级积分TDI CCD影像震颤补偿和畸变校正提供可靠数据基础。
Description
技术领域
本发明涉及卫星遥感技术领域，尤其涉及一种星载TDI CCD影像震颤成像畸变建模方法及装置。
背景技术
TDI CCD（Time Delayed and Integration Charge Coupled Device，时间延迟积分电荷耦合器件）是一种基于时间延迟积分技术的线性CCD，采用同步多级时间累积积分来延长曝光时间。相对于普通线性CCD，TDI CCD可以在低照度条件下，获得高灵敏度、高输出速度、高空间分辨率和宽动态范围。因此，TDI CCD成为了高分辨率光学遥感卫星普遍采用的核心感光器件。TDI CCD通过多级积分成像，能够对同一目标物体进行多次曝光获得光通量。TDI CCD的每一积分级数依次扫描同一区域，电荷转移和累积速度与光学扫描图像同步，直到N行积分成像完成，累积的信号电荷都被转移到水平位移寄存器中并读出图像行。但是在卫星运行过程中，平台震颤会导致CCD运动速度和电荷传输速率不匹配，最终引起每一级成像偏移各不相同，导致不同区域中电荷的连续错位和混叠，造成图像模糊和失真。并且随着TDI CCD积分级数的增加，这种震颤像移对图像质量的影响就愈发严重。例如，高分某卫星配备了全色多光谱五合一的光学相机，全色谱段和多光谱波段共TDI CCD相机焦平面，以一定间隔平行排列在同一相机上。多光谱四个谱段影像的TDI积分级数各不相同，单级积分时间相同，在相同成像条件下，震颤引起多光谱影像在每一级的像移各不相同，导致每个谱段的多光谱影像变形不同；全色谱段和多光谱谱段的TDI积分级数也不相同，并且两者的单级积分时间不一样（一般全色影像空间分辨率是多光谱的四倍，全色影像单级积分时间是多光谱的四分之一），因而相同成像条件下，全色影像和多光谱影像因震颤引起的影像变形程度并不一致，取决于他们的TDI积分级数和单级积分时间。为了实现遥感影像后续的高精度处理与应用，多级积分TDI CCD的平台震颤检测是一个急需解决的关键性问题。
目前，对高分辨率光学卫星CCD影像进行震颤处理的研究主要基于单级CCD和基于多级TDI CCD两种，不论是单级CCD还是多级TDI CCD，都需先进行平台震颤检测，才能进行后续影像的几何定位、变化检测等应用。其中单级CCD震颤检测常基于分时成像特点，震颤检测技术手段已经非常成熟。学者们先后提出基于三线阵立体影像、多光谱影像以及正交图像的震颤检测和补偿的框架，并且建立了精确的单级CCD影像抖动位移估计模型，结果表明都能准确估计出卫星震颤并基于此进行震颤补偿，提高遥感影像的几何质量。而多级积分TDI CCD影像的震颤畸变分析与建模研究相对较少，尚不成熟。初期学者们的研究都停留在震颤与TDI CCD积分级数的定性关系上，发现低频振动产生空间变化的模糊和几何变形，TDI CCD积分级数的增加可以使图像变形程度变小，振动初始相位的改变会引起图像的变形在空域上发生相移。近期有学者通过分析平台震颤引起的多级积分TDI CCD影像变形规律，初次建立了多级积分影像平台震颤检测模型；并且提出卫星平台震颤引起的TDI CCD多级积分后的图像运动模型，实现了多光谱影像间的视差曲线到震颤曲线的转换；以及还有研究学者通过震颤成像仿真开展了基于深度学习的TDI CCD影像盲复原研究，有效提升了图像质量，但并未从退化机理的角度反演出震颤引起的影像几何变形的规律。
发明内容
本发明目的在于提供一种星载TDI CCD影像震颤成像畸变精确建模方法及装置，本发明首先基于TDI CCD多级积分成像特点，分析了多级积分TDI CCD影像震颤成像绝对畸变特性，利用因平台震颤引起的多级积分TDI CCD影像像移公式，仿真多级积分TDI CCD影像像移曲线，得到不同积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；其次再结合TDI CCD分时成像多光谱影像，构建严密视差成像模型，该模型建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系；对多级积分多光谱影像的两个波段影像进行高精度密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变曲线；基于多级积分多光谱影像间的震颤成像相对畸变曲线，对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位（于多光谱影像而言，其一行成像时间正好是一个单位成像时间，即单级积分时间）；最后根据单级积分多光谱影像的震颤成像绝对畸变，全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，对全色影像和多光谱各波段影像分别进行多级积分TDI CCD影像震颤成像绝对畸变精确建模，并对全色影像和多光谱各波段影像进行震颤补偿，校正几何畸变。为实现上述目的，本发明提供如下技术方案：
本发明提供了一种星载TDI CCD影像震颤成像畸变建模方法，所述方法包括，
对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
基于TDI CCD 分时成像多光谱影像和多级积分TDI CCD影像的震颤像移曲线，构建严密视差成像模型；
对所述严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；
基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模。
具体的，所述构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，包括：
仿真因平台震颤引起的多级积分TDI CCD影像的像移曲线；
基于所述像移曲线，得到某一单级积分时间下不同TDI CCD积分级数与多级积分TDI CCD影像震颤成像绝对畸变的振幅和相位的对应关系；
基于所述对应关系，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。
具体的，所述多级积分后TDI CCD影像的震颤像移曲线，表示为：
其中，是在积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子，是积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的相位偏移量，代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率。
具体的，所述严密视差成像模型表示为：
其中，是当前成像时间，是成像时间间隔，，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变；，和，分别是在积分级数成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
具体的，对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，包括：
选取多级积分多光谱影像的两个波段影像进行密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变曲线；
基于所述震颤成像相对畸变曲线对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位。
具体的，所述多级积分多光谱影像间的震颤成像相对畸变曲线公式为：
其中，代表多级积分多光谱影像间的震颤成像相对畸变振幅，是多级积分多光谱影像间的震颤成像相对畸变初相位，是多级积分多光谱影像间的震颤成像相对畸变频率。
具体的，所述单级积分多光谱影像的震颤成像绝对畸变频率与多级积分多光谱影像间的震颤成像相对畸变频率一致；
所述单级积分多光谱影像的震颤成像绝对畸变振幅表示为：
其中，代表单级积分多光谱影像的震颤成像绝对畸变振幅；是多级积分多光谱影像间的震颤成像相对畸变振幅；，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
所述单级积分多光谱影像的震颤成像绝对畸变初相位表示为：
其中，是单级积分多光谱影像的震颤成像绝对畸变初相位；是多级积分多光谱影像间的震颤成像相对畸变初相位。
具体的，多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
对全色影像的多级积分TDI CCD影像震颤成像绝对畸变建模和对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模。
具体的，所述对全色影像的多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系得到全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
根据单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDICCD影像的震颤像移曲线对全色影像进行多级积分TDI CCD影像震颤成像绝对畸变建模；或，
对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系得到多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
根据单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDI CCD影像的震颤像移曲线对多光谱各波段影像进行多级积分TDI CCD影像震颤成像绝对畸变建模。
本发明还提供了一种星载TDI CCD影像震颤成像畸变建模装置，所述装置包括，
第一构建模块，用于对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
获取模块，用于基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
第二构建模块，用于基于TDI CCD 分时成像多光谱影像和多级积分TDI CCD影像的震颤像移曲线，构建严密视差成像模型；
得到模块，用于对所述严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；
建模模块，用于基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模。
具体的，所述装置还包括：
补偿模块，用于根据多级积分TDI CCD影像震颤成像绝对畸变建模结果，对全色影像和多光谱各波段影像分别进行震颤补偿。
一种电子设备，所述电子设备包括至少一个处理器和至少一个存储器，所述存储器与所述处理器数据连接，其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行上述的方法。
本发明的技术效果和优点：
本发明针对TDI CCD多级积分成像特点，根据理论分析和仿真验证得到了不同积分成像时间下的多级积分TDI CCD影像震颤成像畸变振幅衰减因子和相位偏移量；基于TDICCD分时成像多光谱影像，构建严密视差成像模型，建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像震颤成像绝对畸变之间的关系；只要多级积分多光谱影像间的震颤成像相对畸变确定，就可以根据严密视差成像模型反演计算得到单级积分多光谱影像震颤成像绝对畸变，进而对多级积分TDI CCD影像震颤成像绝对畸变进行精确建模，最终进行多级积分TDI CCD影像震颤补偿，校正几何畸变。以真实的遥感卫星数据进行验证实验，实现单级积分多光谱影像震颤成像绝对畸变的检测，全色影像和多光谱影像的震颤成像绝对畸变精确建模，为多级积分TDI CCD影像震颤补偿和震颤几何畸变校正提供可靠数据基础。
该模型建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像震颤成像绝对畸变之间的关系；对多级积分多光谱影像的两个波段影像进行高精度密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变频率、振幅和相位；对严密视差成像模型进行理论推导，得到单级积分多光谱影像震颤成像绝对畸变的振幅和相位计算公式，结合多级积分多光谱影像间的震颤成像相对畸变频率、振幅和相位，计算得到以多光谱单级积分时间为单位成像时间的震颤成像绝对畸变频率、振幅和相位（于多光谱影像而言，其一行成像时间正好是一个单位成像时间）；最后根据单级积分多光谱影像的震颤成像绝对畸变，全色影像和多光谱各波段影像震颤成像畸变的振幅衰减因子和相位偏移量，对全色影像和多光谱各波段影像分别进行震颤成像绝对畸变精确建模，并对全色影像和多光谱各波段影像进行震颤补偿，校正几何畸变。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书以及附图中所指出的结构来实现和获得。
附图说明
图1a为本发明的星载TDI CCD影像震颤成像畸变精确建模方法的流程图;
图1b为本发明的星载TDI CCD影像震颤成像畸变精确建模方法的实施图;
图2a为本发明在两种不同震颤成像绝对畸变频率下多级积分TDI CCD影像震颤成像绝对畸变振幅与的关系图;
图2b为本发明在两种不同震颤成像绝对畸变频率下多级积分TDI CCD影像震颤成像绝对畸变相位与的关系图;
图3a为使用本发明提出的方法进行震颤补偿前的全色影像区域一；
图3b为使用本发明提出的方法进行震颤补偿前的全色影像区域二；
图3c为使用本发明提出的方法进行震颤补偿前的全色影像区域三；
图4a为使用本发明提出的方法进行震颤补偿后的全色影像区域一；
图4b为使用本发明提出的方法进行震颤补偿后的全色影像区域二；
图4c为使用本发明提出的方法进行震颤补偿后的全色影像区域三；
图5为本发明的星载TDI CCD影像震颤成像畸变精确建模与补偿系统图;
图6为本发明具体实施例中的电子设备结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为解决现有技术的不足，本发明公开了一种星载TDI CCD影像震颤成像畸变精确建模方法，如图1a所示，所述方法包括，
步骤1：对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
步骤2：基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
步骤3：基于TDI·CCD 分时成像多光谱影像和多级积分TDI CCD影像的震颤像移曲线，构建严密视差成像模型；该模型建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系；
步骤4：对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；
步骤5：基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模；其中，所述多级积分TDI CCD影像震颤成像绝对畸变建模，包括：对全色影像的多级积分TDICCD影像震颤成像绝对畸变建模和对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模；
结合图1b,本发明方法包括，步骤6：根据所述全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，对全色影像和多光谱各波段影像进行震颤补偿。
本发明技术方案通过理论分析和仿真验证多级积分TDI CCD影像震颤成像绝对畸变规律，得到不同积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；通过对某卫星的全色影像和多光谱影像进行震颤补偿验证了所构建的严密视差成像模型的正确性，以及星载TDI CCD影像震颤成像畸变精确建模方法的可行性。
在本发明的一个具体实施例中，步骤1：对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，具体为：
1）仿真因平台震颤引起的多级积分TDI CCD影像的像移曲线：TDI CCD多级积分过程中每一级成像受到平台震颤的像移随时间而变化，将第一级积分时间作为多级积分合成影像的成像时间，利用多级积分TDI CCD影像像移公式，开展因平台震颤引起的多级积分TDI CCD影像像移仿真实验，仿真因平台震颤引起的多级积分TDI CCD影像的像移曲线。
在本发明的一个具体实施例中，以第一级积分时间作为多级积分影像合成的影像成像时间，平台震颤引起的多级积分TDI CCD影像像移公式为：
其中，为TDI CCD积分级数，为影像行号，为单级积分时间，为每一积分级数下细分数量，代表震颤成像绝对畸变周期，代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，表示在细分数量为时的多级积分TDI CCD影像第行的像移值,k是一个在区间[1,mN-1]里变化的整数。
2）基于所述像移曲线，得到某一单级积分时间下不同TDI CCD积分级数与多级积分TDI CCD影像震颤成像绝对畸变振幅和相位的对应关系；
仿真在相同单级积分多光谱影像的震颤成像绝对畸变振幅，初相位以及单级积分成像时间，不同震颤成像绝对畸变频率、不同TDI CCD积分级数下的多级积分TDI CCD影像像移曲线，然后对像移曲线进行正弦函数拟合，获取不同震颤成像绝对畸变频率下不同TDICCD积分级数的多级积分TDI CCD影像震颤成像绝对畸变振幅和相位参数，从而得到某一单级积分时间下不同TDI CCD积分级数与多级积分TDI CCD影像震颤成像绝对畸变振幅和相位的对应关系；
3）基于所述对应关系，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子、相位偏移量的关系。
TDI CCD积分级数的变化，意味着积分成像时间的变化，同时平台震颤对不同积分级数TDI CCD影像的影响是与震颤成像绝对畸变频率相关的。因此进一步得到（TDICCD积分级数、单级积分时间和震颤成像绝对畸变频率的乘积）与多级积分TDI CCD影像震颤成像绝对畸变振幅和相位的关系；将多级积分TDI CCD影像震颤成像绝对畸变振幅和相位与设置的单级积分多光谱影像的震颤成像绝对畸变振幅和初相位进行比较，得到（TDI CCD积分级数、单级积分时间和震颤成像绝对畸变频率的乘积）与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。
在本发明的一个具体实施例中，步骤2：基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；具体如下：
单级积分多光谱影像的震颤成像绝对畸变曲线可以用正弦函数模型表示：
其中，代表代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率。
基于多级积分TDI CCD影像震颤成像绝对畸变分析，多级积分TDI CCD影像的震颤像移曲线可表示为：
其中，是在积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子，是积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的相位偏移量；代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率。
在本发明的一个具体实施例中，步骤3：基于TDI CCD 分时成像多光谱影像和多级积分TDI CCD影像的震颤像移曲线，构建严密视差成像模型，具体如下：
1）基于分时成像特点和TDI CCD多级积分成像特点，以及引入的震颤振幅衰减因子和相位偏移量，提出严密视差成像模型。严密视差成像模型表示了多级积分多光谱影像间的震颤相对误差与单级积分绝对震颤误差的关系，可用下式表达严密视差成像模型：
其中，是当前成像时间，是成像时间间隔，，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变。
结合多级积分TDI CCD影像的震颤像移曲线，严密视差成像模型可以进一步用如下公式表示：
其中，，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和像移相位偏移量。由上述公式可知，构建的严密视差成像模型建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系。
在本发明的一个具体实施例中，步骤4：对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，具体为：
1）基于TDI CCD影像中不同波段的多级积分多光谱影像，得到多级积分多光谱影像间的震颤成像相对畸变
首先选取TDI CCD影像中多级积分多光谱影像的两个波段影像进行高精度密集影像匹配得到垂轨和沿轨两个方向的视差图；在本发明的一个具体实施例中，选取一组多级积分多光谱影像的最佳波段影像组合进行高精度密集匹配，常用到的影像高精度密集匹配方法有最小二乘影像匹配和相位匹配。本实验采用最小二乘影像匹配方法，先确定多光谱影像重叠区，然后采用最小二乘影像匹配进行逐像素匹配获得重叠区内在垂轨和沿轨两个方向上所有同名点的坐标差，生成两个方向的视差图；
基于视差图，得到沿扫描行的多级积分多光谱影像间的震颤成像相对畸变曲线；根据起始成像时间和积分时间，将扫描行转换为时刻得到多级积分多光谱影像间的震颤成像相对畸变的时变曲线；利用傅里叶变换得到多级积分多光谱影像间的震颤成像相对畸变曲线的频率和振幅，将获得的频率和振幅作为初始值，采用正弦函数模型对多级积分多光谱影像间的震颤成像相对畸变曲线进行整体最小二乘拟合，得到多级积分多光谱影像间的震颤成像相对畸变的频率、振幅和初始相位。
其中，多级积分多光谱影像间的震颤成像相对畸变曲线可以用如下公式表示：
其中，代表多级积分多光谱影像间的震颤成像相对畸变振幅，是多级积分多光谱影像间的震颤成像相对畸变初相位，是多级积分多光谱影像间的震颤成像相对畸变震颤频率。
2）结合多级积分多光谱影像间的震颤成像相对畸变曲线，对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变振幅和相位计算公式。
根据平台周期运动和运动合成原理，单级积分多光谱影像的震颤成像绝对畸变频率应与多级积分多光谱影像间的震颤成像相对畸变频率一致。
单级积分多光谱影像的震颤成像绝对畸变振幅计算公式表示为：
其中，代表单级积分多光谱影像的震颤成像绝对畸变振幅；是多级积分多光谱影像间的震颤成像相对畸变振幅；，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
单级积分多光谱影像的震颤成像绝对畸变相位计算公式表示为：
其中，是单级积分多光谱影像的震颤成像绝对畸变初相位；是多级积分多光谱影像间的震颤成像相对畸变初相位。
在本发明的一个具体实施例中，步骤5：基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模，其中，多级积分TDI CCD影像震颤成像绝对畸变建模，包括：对全色影像的多级积分TDI CCD影像震颤成像绝对畸变建模和对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模，具体为：
基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，根据不同多级积分TDI CCD影像的单级积分时间，TDI CCD 积分级数和震颤成像绝对畸变频率，分别得到全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量、多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDICCD影像的震颤像移曲线对全色影像进行多级积分TDI CCD影像震颤成像绝对畸变精确建模；
基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDI CCD影像的震颤像移曲线分别对多光谱各波段影像进行多级积分TDI CCD影像震颤成像绝对畸变精确建模。
在本发明的一个具体实施例中，步骤6：根据所述全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，对全色影像和多光谱各波段影像进行震颤补偿，具体为：
基于全色影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，利用双线性重采样对全色影像进行震颤补偿；
基于多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，利用双线性重采样分别对多光谱各波段影像进行震颤补偿。
本发明还提供了一种星载TDI CCD影像震颤成像畸变精确建模与补偿装置，如图5所示，所述装置包括，
第一构建模块，用于对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
获取模块，用于基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
第二构建模块，用于基于TDI CCD 分时成像多光谱影像和多级积TDI CCD影像的震颤像移曲线，构建严密视差成像模型；
得到模块，用于对所述严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；
建模模块，用于基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模。
补偿模块，用于根据多级积分TDI CCD影像震颤成像绝对畸变建模结果，对全色影像和多光谱各波段影像分别进行震颤补偿。
所述严密视差成像模型，构建了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系，表示为：
其中，是当前成像时间，是成像时间间隔，，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和像移相位偏移量；代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是震颤成像绝对畸变频率；是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变。
所述对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变振幅和相位计算公式，包括：单级积分多光谱影像的震颤成像绝对畸变振幅关系式：
其中，代表单级积分多光谱影像的震颤成像绝对畸变振幅；是多级积分多光谱影像间的震颤成像相对畸变振幅；，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
单级积分多光谱影像的震颤成像绝对畸变相位关系式：
其中，是单级积分多光谱影像的震颤成像绝对畸变初相位；是多级积分多光谱影像间震颤成像相对畸变初相位。
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
下面将结合具体的实施例对本发明技术方案进一步进行说明。
本发明技术方案通过理论分析和仿真验证得到了不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，通过对某卫星的全色影像和多光谱影像进行震颤补偿验证了所构建的严密视差成像模型的正确性，以及星载TDI CCD影像震颤成像畸变精确建模方法的可行性。
以下给出一个具体的实施例：
图1a表示的是星载TDI CCD影像震颤成像畸变精确建模方法的流程示意图，主要包括以下步骤：
步骤1：多级积分TDI CCD影像震颤成像绝对畸变分析。
根据多级积分TDI CCD影像像移仿真分析，得到不同积分成像时间与多级积分TDICCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，，即与多级积分TDICCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。
1）仿真因平台震颤引起的多级积分TDI CCD影像的像移曲线。
对于多级积分TDI CCD影像，最终输出的扫描行是由多级积分影像合成的影像，TDI CCD多级积分过程中每一级成像因平台震颤引起的像移随时间而变化，故其受到的平台震颤影响是每一级影像受到影响的合成效果，合成后的影像受到的影响不再与单级积分影像相同。若TDI CCD积分级数为，以第一级积分时间作为多级积分影像合成的影像成像时间，平台震颤引起的多级积分TDI CCD影像像移公式如下：
(1)
其中，为TDI CCD积分级数，为影像行号，为单级积分时间，为每一积分级数下细分数量，代表震颤成像绝对畸变周期，代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，表示在细分数量为时的多级积分TDI CCD影像第行的像移值，本发明中最佳取值为10。
2）获取某一单级积分时间下不同TDI CCD积分级数与多级积分TDI CCD影像震颤成像绝对畸变的振幅和相位的对应关系
利用公式（1）仿真在相同单级积分多光谱影像的震颤成像绝对畸变振幅，初相位以及单级积分成像时间，不同震颤成像绝对畸变频率、不同TDI CCD积分级数下的多级积分TDI CCD影像像移曲线，然后对像移曲线进行正弦函数拟合，获取不同震颤成像绝对畸变频率下不同TDI CCD积分级数的多级积分TDI CCD影像震颤成像绝对畸变振幅和相位参数，从而得到某一单级积分时间下不同TDI CCD积分级数与多级积分TDI CCD影像震颤成像绝对畸变的振幅和相位的对应关系；
其中，假设单级积分多光谱影像的震颤成像绝对畸变振幅为1个像素，初始相位为0，单级积分时间为0.00036s，在震颤成像绝对畸变频率为50Hz和100Hz的震颤条件下，从单级积分到167级积分，仿真因平台震颤引起的多级积分TDI CCD影像像移曲线，得到单级积分时间为0.00036s时不同震颤成像绝对畸变频率下不同TDI CCD积分级数与多级积分TDICCD影像震颤成像绝对畸变的振幅和相位的对应关系。
3）获取不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。
由于平台震颤对不同TDI CCD积分级数的多级积分TDI CCD影像的影响是与积分时间和震颤成像绝对畸变频率相关的，进一步得到多级积分TDI CCD影像震颤成像绝对畸变振幅以及相位与（TDI CCD积分级数、单级积分时间和震颤成像绝对畸变频率的乘积）的关系曲线图，如图2a和2b所示。由图2a可知，当单级积分时间和震颤成像绝对畸变频率一定时，随着TDI CCD积分级数的增加，多级积分TDI CCD影像震颤成像绝对畸变振幅逐步衰减；同理，由图2b可知，随着积分级数的增加，多级积分TDI CCD影像震颤成像绝对畸变的相位存在偏移；相位的偏移速度也随着震颤成像绝对畸变频率变化，并且高频震颤引起的相位偏移存在周期性。
将多级积分TDI CCD影像震颤成像绝对畸变振幅和相位与设置的单级积分多光谱影像的震颤成像绝对畸变振幅和初相位进行相比较，得到不同积分成像时间下的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，即（TDI CCD积分级数、单级积分时间和震颤成像绝对畸变频率的乘积）与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。由于该实施例中单级积分多光谱影像震颤成像绝对畸变振幅为1pixel，初相位为0，因此多级积分TDI CCD影像震颤成像绝对畸变振幅和相位就是多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
因此当TDI CCD积分级数、单级积分时间和震颤成像绝对畸变频率都确定时，根据的值来获取对应多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
步骤2：基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
基于多级积分TDI CCD影像震颤成像绝对畸变分析，多级积分TDI CCD影像的震颤像移曲线可表示为：
（2）
其中，是在积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子，是积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的相位偏移量；代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率。
步骤3：构建严密视差成像模型。
基于TDI CCD分时成像多光谱影像，构建的严密视差成像模型表示为：
(3)
其中，是当前成像时间，是成像时间间隔，，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变。
结合公式（2）和（3），严密视差成像模型进一步表示了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系，可以用如下公式表示：
(4)
其中，，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
步骤4：基于严密视差成像模型的单级积分多光谱影像的震颤成像绝对畸变反演。先基于TDI CCD影像中不同波段的多级积分多光谱影像，得到多级积分多光谱影像间的震颤成像相对畸变曲线；基于多级积分多光谱影像间的震颤成像相对畸变曲线，对严密视差成像模型进行反演得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位。
1）基于TDI CCD影像中不同波段的多级积分多光谱影像，得到多级积分多光谱影像间的震颤成像相对畸变
高精度密集匹配。利用高分某卫星的B1和B3两个波段的多级积分多光谱影像进行高精密度匹配，本实验采用最小二乘影像匹配方法。确定多光谱影像重叠区，然后采用最小二乘影像匹配进行逐像素匹配获得重叠区内在垂轨和沿轨两个方向上所有同名点的坐标差，生成两个方向的视差图。
多级积分多光谱影像间的震颤成像相对畸变分析。基于视差图，得到沿扫描行的多级积分多光谱影像间的震颤成像相对畸变曲线；根据起始成像时间和积分时间，将扫描行转换为时刻得到多级积分多光谱影像间的震颤成像相对畸变的时变曲线；利用傅里叶变换得到多级积分多光谱影像间的震颤成像相对畸变曲线的频率和振幅，将获得的频率和振幅作为初始值，采用正弦函数模型对多级积分多光谱影像间的震颤成像相对畸变曲线进行整体最小二乘拟合，得到多级积分多光谱影像间的震颤成像相对畸变的频率、振幅和初始相位。
其中，多级积分多光谱影像间的震颤成像相对畸变曲线可以用如下公式表示：
(5)
其中，代表多级积分多光谱影像间的震颤成像相对畸变振幅，是多级积分多光谱影像间的震颤成像相对畸变初相位，是多级积分多光谱影像间的震颤成像相对畸变频率。
2）结合公式（4）和公式（5），对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变振幅和相位计算公式。
根据平台周期运动和运动合成原理，单级积分多光谱影像的震颤成像绝对畸变频率应与多级积分多光谱影像间的震颤成像相对畸变频率一致。
单级积分多光谱影像的震颤成像绝对畸变振幅计算公式表示为：
(6)
其中，代表单级积分多光谱影像的震颤成像绝对畸变振幅；是多级积分多光谱影像间的震颤成像相对畸变振幅；，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量。
单级积分多光谱影像的震颤成像绝对畸变相位计算公式表示为：
(7)
其中，是单级积分多光谱影像的震颤成像绝对畸变初相位；是多级积分多光谱影像间的震颤成像相对畸变初相位。
3）该实施例中使用的高分某卫星的多级积分多光谱影像的B1、B3波段的积分级数分别为48和24，单级积分时间为0.000362395311s，震颤成像绝对畸变频率均为100Hz；基于步骤1的分析，根据（TDI CCD积分级数、单级积分时间和震颤成像绝对畸变的乘积）的值来获取多级积分多光谱影像B1和B3波段对应的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；然后根据多级积分多光谱影像间的震颤成像相对畸变频率、振幅和初相位，结合公式（6）和公式（7）计算得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位。
步骤5：多级积分TDI CCD影像震颤成像绝对畸变精确建模。根据单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，以及多级积分TDI CCD影像（含全色影像和多光谱影像）震颤成像绝对畸变的振幅衰减因子和相位偏移量，根据公式（2）对多级积分TDI CCD影像进行震颤成像绝对畸变精确建模。
1）基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，根据不同多级积分TDI CCD影像（TDI CCD 积分级数、单级积分时间和震颤成像绝对畸变的乘积）的值，分别得到全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；其中该实施例中，全色影像积分级数为48，单级积分时间为0.00009060s，多光谱影像的B1、B2及B3波段的积分级数分别为48、48和24，单级积分时间为0.000362395311s，震颤成像绝对畸变频率均为100Hz；
2）基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用公式（2）对全色影像进行多级积分TDI CCD影像震颤成像绝对畸变精确建模；
3）基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用公式（2）分别对多光谱各波段影像进行多级积分TDI CCD影像震颤成像绝对畸变精确建模。
步骤6：多级积分TDI CCD影像震颤补偿。根据全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，对全色影像和多光谱各波段影像进行震颤补偿。
1）基于全色影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，利用双线性重采样对全色影像进行震颤补偿；
2）基于多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果，利用双线性重采样分别对多光谱各波段影像进行震颤补偿。
其中震颤补偿前的全色影像如图3a-3c所示，图3a-图3c为全色影像的三个不同区域，震颤补偿后的全色影像如图4a-图4c所示，图4a-图4c为图3a-3c中相对应的进行震颤补偿后的三个区域。
综上所述，本发明针对TDI CCD多级积分成像特点，根据理论分析和仿真验证得到了不同积分成像时间下的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；基于TDI CCD分时成像多光谱影像，构建严密视差成像模型，建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系；对不同波段的多级积分多光谱影像进行密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变曲线；基于多级积分多光谱影像间的震颤成像相对畸变曲线对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；进而对多级积分TDI CCD影像震颤成像绝对畸变进行精确建模，最终进行多级积分TDI CCD影像震颤补偿，校正几何畸变，提供了星载TDI CCD影像震颤成像畸变精确建模方法及系统。
本发明提供了一种星载TDI CCD影像震颤成像畸变精确建模方法及系统，主要包括六大步骤。第一步是多级积分TDI CCD影像震颤成像绝对畸变分析。利用多级积分TDICCD影像像移公式，仿真因平台震颤引起的多级积分TDI CCD影像像移曲线，得到不同积分成像时间（TDI CCD 积分级数、单级积分时间和震颤成像绝对畸变频率的乘积）与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。第二步是获取多级积分TDI CCD影像的震颤像移曲线。第三步是构建严密视差成像模型。基于TDI CCD分时成像多光谱影像，构建严密视差成像模型，该模型建立了多级积分多光谱影像间的震颤成像相对畸变与单级积分多光谱影像的震颤成像绝对畸变之间的关系。第四步是基于严密视差成像模型的单级积分多光谱影像的震颤成像绝对畸变反演。选取多级积分多光谱影像的两个波段影像进行密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变曲线；然后基于多级积分多光谱影像间的震颤成像相对畸变曲线，对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位。第五步是多级积分TDI CCD影像震颤成像绝对畸变精确建模。基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，以及全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，对全色影像和多光谱各波段影像进行多级积分TDI CCD影像震颤成像绝对畸变精确建模。最后是多级积分TDI CCD影像震颤补偿。利用全色影像和多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模结果分别对全色影像和多光谱影像进行震颤补偿，校正影像几何畸变。
同时，本发明以真实的遥感卫星数据进行验证实验，实现单级积分多光谱影像的震颤成像绝对畸变的检测及全色影像和多光谱影像的多级积分TDI CCD影像震颤成像绝对畸变精确建模，为多级积分TDI CCD影像震颤补偿提供可靠数据基础。
需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。
基于上述公开的内容，对应地本发明还提供了一种电子设备。如图6所示，本公开实施例的电子设备包括电气连接的至少一个处理器和只要一个存储器，所述存储器与所述处理器电气连接，其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行如上控制器执行的所述方法步骤。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (12)
1.一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，所述方法包括，
对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDICCD影像的震颤像移曲线；
基于TDI CCD 分时成像多光谱影像和多级积分TDI CCD影像的震颤像移曲线，构建严密视差成像模型，基于分时成像特点和TDI CCD多级积分成像特点，以及引入的震颤振幅衰减因子和相位偏移量，提出严密视差成像模型；严密视差成像模型表示多级积分多光谱影像间的震颤相对误差与单级积分绝对震颤误差的关系；
对所述严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；其中，对所述严密视差成像模型进行反演，包括：基于 TDICCD 影像中不同波段的多级积分多光谱影像，确定多级积分多光谱影像间的震颤成像相对畸变；根据多级积分多光谱影像间的震颤成像相对畸变与多级积分多光谱影像的震颤成像绝对畸变的关系，得到单级积分多光谱影像的震颤成像绝对畸变振幅和相位计算公式；
基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDICCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模。
2.根据权利要求1所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，所述构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，包括：
仿真因平台震颤引起的多级积分TDI CCD影像的像移曲线；
基于所述像移曲线，得到某一单级积分时间下不同TDI CCD积分级数与多级积分TDICCD影像震颤成像绝对畸变的振幅和相位的对应关系；
基于所述对应关系，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系。
3.根据权利要求1所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，
所述多级积分TDI CCD影像的震颤像移曲线，表示为：
其中，是在积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子，是积分级数为的积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的相位偏移量，代表单级积分多光谱影像的震颤成像绝对畸变振幅，是单级积分多光谱影像的震颤成像绝对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率，表示多级积分TDI CCD影像的震颤像移,t是当前成像时间。
4.根据权利要求1或3所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，所述严密视差成像模型表示为：
其中，t是当前成像时间，是成像时间间隔，，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变；，和，分别是在积分级数成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变。
5.根据权利要求1所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，包括：
选取多级积分多光谱影像的两个波段影像进行密集匹配，得到多级积分多光谱影像间的震颤成像相对畸变曲线；
基于所述震颤成像相对畸变曲线对严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位。
6.根据权利要求5所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，所述多级积分多光谱影像间的震颤成像相对畸变曲线公式为：
其中，代表多级积分多光谱影像间的震颤成像相对畸变振幅，是多级积分多光谱影像间的震颤成像相对畸变初相位，是单级积分多光谱影像的震颤成像绝对畸变频率，是因平台震颤引起的多级积分多光谱影像间的震颤成像相对畸变,t是当前成像时间。
7.根据权利要求5或6所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，
所述单级积分多光谱影像的震颤成像绝对畸变频率与多级积分多光谱影像间的震颤成像相对畸变频率一致；
所述单级积分多光谱影像的震颤成像绝对畸变振幅表示为：
其中，代表单级积分多光谱影像的震颤成像绝对畸变振幅；是多级积分多光谱影像间的震颤成像相对畸变振幅；，和，分别是在积分级数积分成像时间下多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量,是成像时间间隔,是单级积分多光谱影像的震颤成像绝对畸变频率；
所述单级积分多光谱影像的震颤成像绝对畸变初相位表示为：
其中，是单级积分多光谱影像的震颤成像绝对畸变初相位；是多级积分多光谱影像间的震颤成像相对畸变初相位。
8.根据权利要求1所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
对全色影像的多级积分TDI CCD影像震颤成像绝对畸变建模和对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模。
9.根据权利要求8所述的一种星载TDI CCD影像震颤成像畸变建模方法，其特征在于，所述对全色影像的多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系得到全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
根据单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，全色影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDI CCD影像的震颤像移曲线对全色影像进行多级积分TDI CCD影像震颤成像绝对畸变建模；或，
对多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变建模，包括：
基于不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系得到多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量；
根据单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位，多光谱各波段影像的多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量，利用多级积分TDI CCD影像的震颤像移曲线对多光谱各波段影像进行多级积分TDI CCD影像震颤成像绝对畸变建模。
10.一种星载TDI CCD影像震颤成像畸变建模装置，其特征在于，所述装置包括，
第一构建模块，用于对多级积分TDI CCD影像震颤成像绝对畸变分析，构建不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系；
获取模块，用于基于单级积分多光谱影像的震颤成像绝对畸变曲线和不同积分成像时间与多级积分TDI CCD影像震颤成像绝对畸变的振幅衰减因子和相位偏移量的关系，获取多级积分TDI CCD影像的震颤像移曲线；
第二构建模块，用于基于TDI CCD 分时成像多光谱影像和多级积TDI CCD影像的震颤像移曲线，构建严密视差成像模型；基于分时成像特点和TDI CCD多级积分成像特点，以及引入的震颤振幅衰减因子和相位偏移量，提出严密视差成像模型；严密视差成像模型表示多级积分多光谱影像间的震颤相对误差与单级积分绝对震颤误差的关系；
得到模块，用于对所述严密视差成像模型进行反演，得到单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位；其中，对所述严密视差成像模型进行反演，包括：基于TDICCD 影像中不同波段的多级积分多光谱影像，确定多级积分多光谱影像间的震颤成像相对畸变；根据多级积分多光谱影像间的震颤成像相对畸变与多级积分多光谱影像的震颤成像绝对畸变的关系，得到单级积分多光谱影像的震颤成像绝对畸变振幅和相位计算公式；
建模模块，用于基于单级积分多光谱影像的震颤成像绝对畸变频率、振幅和初相位与多级积分TDI CCD影像的震颤像移曲线，进行多级积分TDI CCD影像震颤成像绝对畸变建模。
11.根据权利要求10所述的一种星载TDI CCD影像震颤成像畸变建模装置，其特征在于，所述装置还包括：
补偿模块，用于根据多级积分TDI CCD影像震颤成像绝对畸变建模结果，对全色影像和多光谱各波段影像分别进行震颤补偿。
12.一种电子设备，所述电子设备包括至少一个处理器和至少一个存储器，所述存储器与所述处理器数据连接，其中，
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-9中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311784666.6A CN117475087B (zh) | 2023-12-22 | 2023-12-22 | 一种星载tdi ccd影像震颤成像畸变建模方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311784666.6A CN117475087B (zh) | 2023-12-22 | 2023-12-22 | 一种星载tdi ccd影像震颤成像畸变建模方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117475087A CN117475087A (zh) | 2024-01-30 |
| CN117475087B true CN117475087B (zh) | 2024-04-26 |
Family
ID=89631540
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311784666.6A Active CN117475087B (zh) | 2023-12-22 | 2023-12-22 | 一种星载tdi ccd影像震颤成像畸变建模方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117475087B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103983343A (zh) * | 2014-05-29 | 2014-08-13 | 武汉大学 | 一种基于多光谱影像的卫星平台震颤检测方法及系统 |
| JP2016140576A (ja) * | 2015-02-02 | 2016-08-08 | 株式会社ソシオネクスト | 情報処理装置、振戦情報表示方法、情報処理システム、及びプログラム |
| CN106959454A (zh) * | 2017-03-20 | 2017-07-18 | 上海航天控制技术研究所 | 一种基于数字域tdi及连续复线阵成像模式的颤振反演方法 |
| CN107301628A (zh) * | 2017-06-27 | 2017-10-27 | 武汉大学 | 基于震颤像移轨迹的卫星平台震颤影像去模糊方法 |
| CN113870321A (zh) * | 2021-09-08 | 2021-12-31 | 同济大学 | 一种遥感影像的颤振变化仿真方法和系统 |
| CN114170192A (zh) * | 2021-12-10 | 2022-03-11 | 武汉工程大学 | 星载光学相机焦面震颤精细检测方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101430963B1 (ko) * | 2007-12-06 | 2014-09-23 | 삼성전자주식회사 | 손떨림 보정장치를 구비한 촬영장치 및 광학계 |
| CN110796641B (zh) * | 2019-10-08 | 2022-02-01 | 武汉大学 | 基于连续快照模型的高分辨率卫星影像震颤检测方法 |
-
2023
- 2023-12-22 CN CN202311784666.6A patent/CN117475087B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103983343A (zh) * | 2014-05-29 | 2014-08-13 | 武汉大学 | 一种基于多光谱影像的卫星平台震颤检测方法及系统 |
| JP2016140576A (ja) * | 2015-02-02 | 2016-08-08 | 株式会社ソシオネクスト | 情報処理装置、振戦情報表示方法、情報処理システム、及びプログラム |
| CN106959454A (zh) * | 2017-03-20 | 2017-07-18 | 上海航天控制技术研究所 | 一种基于数字域tdi及连续复线阵成像模式的颤振反演方法 |
| CN107301628A (zh) * | 2017-06-27 | 2017-10-27 | 武汉大学 | 基于震颤像移轨迹的卫星平台震颤影像去模糊方法 |
| CN113870321A (zh) * | 2021-09-08 | 2021-12-31 | 同济大学 | 一种遥感影像的颤振变化仿真方法和系统 |
| CN114170192A (zh) * | 2021-12-10 | 2022-03-11 | 武汉工程大学 | 星载光学相机焦面震颤精细检测方法 |
Non-Patent Citations (6)
| Title |
|---|
| Jitter Detection and Image Restoration Based on Continue Dynamic Shooting Model for High-Resolution TDI CCD Satellite Images;Jun Pan et al;IEEE Transactions on Geoscience and Remote Sensing;20201203;第4915-4933页 * |
| Jitter Detection and Image Restoration Based on Generative Adversarial Networks in Satellite Images;zilin wang et al;Sensors;20210709;第1-14页 * |
| TDI遥感图像的颤振退化仿真和图像复原方法;苏慧;中国优秀硕士论文电子期刊网;20190415;第3-20页 * |
| 利用多光谱影像检测资源三号卫星平台震颤;朱映 等;测绘学报;20150415;44(04);第399-413页 * |
| 卫星平台振动对色散型光谱仪光谱成像质量的影响分析;王晓燕;唐义;唐秋艳;王静;倪国强;;光谱学与光谱分析;20111215(12);第3407-3410页 * |
| 高分辨率光学遥感卫星影像平台震颤处理方法研究;朱映;中国博士学位论文电子期刊网;20200115;第3-35页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117475087A (zh) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7764848B2 (en) | High resolution enabling apparatus and method | |
| JP3588290B2 (ja) | デジタル画像のサブピクセル・アライメント方法 | |
| JP4151793B2 (ja) | 撮像装置および画像の高解像化方法 | |
| EP1547012A1 (en) | Imaging and measurement system | |
| JP4941565B2 (ja) | 対応点探索装置および対応点探索方法 | |
| CN114170192A (zh) | 星载光学相机焦面震颤精细检测方法 | |
| CN115469308A (zh) | 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质 | |
| CN110160468B (zh) | 一种针对运动对象的散焦光栅投影三维测量方法 | |
| CN117475087B (zh) | 一种星载tdi ccd影像震颤成像畸变建模方法及装置 | |
| CN115086550A (zh) | 元成像方法与系统 | |
| CN109029379B (zh) | 一种高精度小基高比立体测绘方法 | |
| JP2017059998A (ja) | 画像処理装置およびその方法、並びに、撮像装置 | |
| CN111222544A (zh) | 一种卫星颤振对相机成像影响的地面模拟测试系统 | |
| CN113920046B (zh) | 一种多分片卫星影像拼接和几何模型构建方法 | |
| CN114485417B (zh) | 一种结构振动位移识别方法及系统 | |
| CN113870321A (zh) | 一种遥感影像的颤振变化仿真方法和系统 | |
| JP2013021636A (ja) | 画像処理装置および方法、学習装置および方法、プログラム、並びに記録媒体 | |
| CN112365400B (zh) | 一种快速光场角度超分辨重建方法 | |
| CN112348947B (zh) | 一种基于参考信息辅助的深度学习的三维重构方法 | |
| CN105374009A (zh) | 遥感影像拼接方法及装置 | |
| JP2022081296A (ja) | 座標変換システム、座標変換方法及びプログラム | |
| JP5528158B2 (ja) | 色ずれ補正装置 | |
| JP3693168B2 (ja) | 変位測定方法およびそれを用いた変位測定装置 | |
| CN118154411B (zh) | 数字自适应光学架构与系统 | |
| Heng et al. | High accuracy registration of translated and rotated images using hierarchical method |
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 |