本发明专利技术公开了一种基于焦平面探测器MTF的微扫描图像重构方法，该方法基于焦平面探测器调制传递函数MTF构建Wiener滤波器MW[u，v]，直流分量归一化后，对过采样图像gos[m，n]进行Wiener滤波，复原出具有更高图像质量的重构图像；本发明专利技术在构建Wiener滤波器的过程中，考虑了焦平面探测器的占空比，克服了现有亚像元重构算法中假设焦平面探测器的占空比为100％的限制；另外利用焦平面探测器调制传递函数MTF构建的Wiener滤波器进行滤波，抑制了焦平面探测器的积分退化效应，在帧间过采样图像重构技术的基础上进一步提高图像质量，从而进一步提高微扫描成像系统的成像性能。
【技术实现步骤摘要】
本专利技术涉及一种基于焦平面探测器MTF的微扫描图像重构方法，属于数字图像处理中的图像重构领域。
技术介绍
随着科学技术的发展，CXD (电荷耦合器件)、CMOS (互补金属_氧化物_半导体)以及红外焦平面阵列等面阵探测器越来越多地应用在工业生产军事侦察、卫星遥感等诸多领域。随着加工工艺的不断发展，CXD等面阵探测器的像元尺寸逐渐变小，面阵规模逐渐变大。像元尺寸的越小，越接近脉冲采样，然而其集光能力越差，散粒噪声越严重，这已经成为成像系统分辨力的瓶颈。此外，相对于CCD和CMOS，红外焦平面阵列由于材料的限制，像元尺寸更大，所以其积分效应更为明显。由于焦平面探测器的像元具有有限大小，因此对应像素值与辐射到整个像元上的能量和成正比，这使像元无法区分更精细的局部图像，即像元对局部能量分布的积分过程使得采集图像发生模糊退化，我们将这种退化现象称之为积分退化效应。为了克服探测器离散采样的限制，国内外先后出现了多种多帧图像超分辨率重建技术，例如频域法、凸集投影法、最大后验概率法、非均匀插值法以及学习法。然而，这些算法复杂性高，难以实现实时处理。微扫描成像系统是一种利用微扫描机构获得多帧具有帧间微位移的成像装置，它在不改变原有成像系统基本组成的情况下，只在光路中添加精密的微扫描装置即可用最低的成本提高系统分辨率。针对微扫描成像系统获得的图像，已经提出了多种亚像元重构算法。大多数亚像元重构算法都假设焦平面探测器的占空比为100%，或者假设局部图像灰度相等。然而，这些假设都与实际情况不符。由于焦平面探测器占空比的存在，使得上述假设出现大量噪声，甚至使得重构算法失效。为此，考虑占空比条件下的微扫描图像重构是亟待解决的问题。
技术实现思路
本专利技术提出了一种基于焦平面探测器MTF的微扫描图像重构方法，基于焦平面探测器调制传递函数MTF构建Wiener滤波器Mw，直流分量归一化后，对过采样图像gos进行Wiener滤波，复原出具有更高图像质量的重构图像，抑制了焦平面探测器的积分退化效应。本专利技术的目的是通过下述技术方案实现的:步骤一:利用基于离散傅里叶变换DFT的亚像素图像配准方法，计算同一场景多帧图像间的帧间位移量；步骤二:在步骤一获得的帧间位移量的基础上，利用帧间差分过采样方法，将多帧图像重构成一巾贞过采样图像gmDn, η]；步骤三:基于焦平面探 测器的调制传递函数MTF构造Wiener滤波器Mw，直流分量归一化后得到Wiener滤波器M' w,用其对过采样图像gmQii, η]进行Wiener滤波，复原出重构图像./>Μφ完成微扫描图像的重构；具体过程为:步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H，其中为频域坐标，U为频域上的水平坐标，V为频域上的竖直坐标:假设f (X，y)为经过光学系统成像在焦平面探测器上的模拟图像，g为经焦平面探测器积分和采样后得到的退化图像，则权利要求1.一种基于焦平面探测器MTF的微扫描图像重构方法，其特征在于， 步骤一:利用基于离散傅里叶变换DFT的亚像素图像配准方法，计算同一场景多帧图像间的帧间位移量； 步骤二:在步骤一获得的帧间位移量的基础上，利用帧间差分过采样方法，将多帧图像重构成一巾贞过采样图像gjm, η]； 步骤三:基于焦平面探测器的调制传递函数MTF构造Wiener滤波器Mw，直流分量归一化后得到Wiener滤波器M' w,用其对过采样图像gjm,]!]进行Wiener滤波，复原出重构图像/，完成微扫描图像的重构；具体过程为: 步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H，其中为频域坐标，U为频域上的水平坐标，V为频域上的竖直坐标: 假设f (X，y)为经过光学系统成像在焦平面探测器上的模拟图像，g为经焦平面探测器积分和采样后得到的退化图像，则全文摘要本专利技术公开了一种基于焦平面探测器MTF的微扫描图像重构方法，该方法基于焦平面探测器调制传递函数MTF构建Wiener滤波器MW，直流分量归一化后，对过采样图像gos进行Wiener滤波，复原出具有更高图像质量的重构图像；本专利技术在构建Wiener滤波器的过程中，考虑了焦平面探测器的占空比，克服了现有亚像元重构算法中假设焦平面探测器的占空比为100％的限制；另外利用焦平面探测器调制传递函数MTF构建的Wiener滤波器进行滤波，抑制了焦平面探测器的积分退化效应，在帧间过采样图像重构技术的基础上进一步提高图像质量，从而进一步提高微扫描成像系统的成像性能。文档编号G06T5/00GK103236045SQ201310157478公开日2013年8月7日 申请日期2013年5月2日 优先权日2013年5月2日专利技术者金伟其, 路陆, 顿雄, 王霞, 范永杰, 李力 申请人:北京理工大学本文档来自技高网...
【技术保护点】
一种基于焦平面探测器MTF的微扫描图像重构方法，其特征在于，步骤一：利用基于离散傅里叶变换DFT的亚像素图像配准方法，计算同一场景多帧图像间的帧间位移量；步骤二：在步骤一获得的帧间位移量的基础上，利用帧间差分过采样方法，将多帧图像重构成一帧过采样图像gos[m，n]；步骤三：基于焦平面探测器的调制传递函数MTF构造Wiener滤波器MW[u，v]，直流分量归一化后得到Wiener滤波器M′W[u，v]，用其对过采样图像gos[m，n]进行Wiener滤波，复原出重构图像完成微扫描图像的重构；具体过程为：步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H[u，v]，其中[u，v]为频域坐标，u为频域上的水平坐标，v为频域上的竖直坐标：假设f(x，y)为经过光学系统成像在焦平面探测器上的模拟图像，g[m，n]为经焦平面探测器积分和采样后得到的退化图像，则g[m,n]=1pq∫nb-q/2nb+q/2∫ma-p/2ma+p/2f(x,y)dxdy=1pq∫-q/2q/2∫-p/2p/2f(x+ma,y+nb)dxdym=1,2,...,M,n=1,2,...,N---(1)式中，a和b为焦平面探测器的水平和竖直像元间隔，p×q为像元尺寸，p为像元在水平方向的尺寸，q为像元在竖直方向的尺寸，M×N为焦平面探测器像面规格，M为焦平面探测器像面在水平方向的尺寸，N为焦平面探测器像面在竖直方向的尺寸，(x，y)是模拟图像坐标系中的坐标，[m，n]是离散图像坐标系中的坐标，m为离散图像坐标系中的水平坐标，最大取值范围为焦平面探测器像面在水平方向的尺寸，n为离散图像坐标系中的竖直坐标，最大取值范围为焦平面 探测器像面在竖直方向的尺寸；对式(1)进行DFT，得式中，F[u，v]和G[u，v]分别表示f[m，n]和g[m，n]的DFT，α＝p/a和β＝q/b分别表示焦平面探测器在水平和竖直方向的占空比；表示DFT运算；和分别为向下取整和向上取整运算；sine函数的定义为sinc(x)＝sin(πx)/(πx)；由式(2)可以获得焦平面探测器的积分退化函数H[u，v]H[u,v]=sinc(uMα)sinc(vNβ),---(3)步骤302、基于式(3)所表示的焦平面探测器的积分退化函数H[u，v]，构建Wiener滤波器MW[u，v]MW[u,v]=H*[u,v]H2[u,v]+Γ---(4)式中，H[u，v]为积分退化函数；H*[u，v]为H[u，v]的复共轭；Γ为常数，它并不随频率变化，而随着占空比、噪声等级以及图像内容变化；将MW[u，v]用直流分量归一化后得到归一化后Wiener滤波器M′W[u，v]MW′[u,v]=MW[u,v]MW[0,0]---(5)式中，MW[0，0]是MW[u，v]在u＝0，v＝0情况下的值，即MW[u，v]的直流分量；步骤303、利用直流分量归一化后的Wiener滤波器M′W[u，v]，对过采样图 像gos[m，n]进行滤波，获得重建图像式中，表示DFT运算，表示DFT逆运算。FDA00003132924400011.jpg,FDA00003132924400021.jpg,FDA00003132924400022.jpg,FDA00003132924400023.jpg,FDA00003132924400024.jpg,FDA00003132924400025.jpg,FDA00003132924400026.jpg,FDA00003132924400028.jpg,FDA00003132924400031.jpg,FDA00003132924400032.jpg,FDA00003132924400033.jpg,FDA00003132924400034.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员：金伟其，路陆，顿雄，王霞，范永杰，李力，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。