CN103236045A - 基于焦平面探测器mtf的微扫描图像重构方法 - Google Patents

Abstract

本发明公开了一种基于焦平面探测器MTF的微扫描图像重构方法，该方法基于焦平面探测器调制传递函数MTF构建Wiener滤波器M_W[u，v]，直流分量归一化后，对过采样图像g_os[m，n]进行Wiener滤波，复原出具有更高图像质量的重构图像；本发明在构建Wiener滤波器的过程中，考虑了焦平面探测器的占空比，克服了现有亚像元重构算法中假设焦平面探测器的占空比为100％的限制；另外利用焦平面探测器调制传递函数MTF构建的Wiener滤波器进行滤波，抑制了焦平面探测器的积分退化效应，在帧间过采样图像重构技术的基础上进一步提高图像质量，从而进一步提高微扫描成像系统的成像性能。

Description

基于焦平面探测器MTF的微扫描图像重构方法

技术领域

本发明涉及一种基于焦平面探测器MTF的微扫描图像重构方法，属于数字图像处理中的图像重构领域。

背景技术

随着科学技术的发展，CCD(电荷耦合器件)、CMOS(互补金属-氧化物-半导体)以及红外焦平面阵列等面阵探测器越来越多地应用在工业生产军事侦察、卫星遥感等诸多领域。随着加工工艺的不断发展，CCD等面阵探测器的像元尺寸逐渐变小，面阵规模逐渐变大。像元尺寸的越小，越接近脉冲采样，然而其集光能力越差，散粒噪声越严重，这已经成为成像系统分辨力的瓶颈。此外，相对于CCD和CMOS，红外焦平面阵列由于材料的限制，像元尺寸更大，所以其积分效应更为明显。由于焦平面探测器的像元具有有限大小，因此对应像素值与辐射到整个像元上的能量和成正比，这使像元无法区分更精细的局部图像，即像元对局部能量分布的积分过程使得采集图像发生模糊退化，我们将这种退化现象称之为积分退化效应。

为了克服探测器离散采样的限制，国内外先后出现了多种多帧图像超分辨率重建技术，例如频域法、凸集投影法、最大后验概率法、非均匀插值法以及学习法。然而，这些算法复杂性高，难以实现实时处理。微扫描成像系统是一种利用微扫描机构获得多帧具有帧间微位移的成像装置，它在不改变原有成像系统基本组成的情况下，只在光路中添加精密的微扫描装置即可用最低的成本提高系统分辨率。

针对微扫描成像系统获得的图像，已经提出了多种亚像元重构算法。大多数亚像元重构算法都假设焦平面探测器的占空比为100％，或者假设局部图像灰度相等。然而，这些假设都与实际情况不符。由于焦平面探测器占空比的存在，使得上述假设出现大量噪声，甚至使得重构算法失效。为此，考虑占空比条件下的微扫描图像重构是亟待解决的问题。

发明内容

本发明提出了一种基于焦平面探测器MTF的微扫描图像重构方法，基于焦平面探测器调制传递函数MTF构建Wiener滤波器M_W[u，v]，直流分量归一化后，对过采样图像g_os[m，n]进行Wiener滤波，复原出具有更高图像质量的重构图像，抑制了焦平面探测器的积分退化效应。

本发明的目的是通过下述技术方案实现的：

步骤一：利用基于离散傅里叶变换DFT的亚像素图像配准方法，计算同一场景多帧图像间的帧间位移量；

步骤二：在步骤一获得的帧间位移量的基础上，利用帧间差分过采样方法，将多帧图像重构成一帧过采样图像g_os[m，n]；

步骤三：基于焦平面探测器的调制传递函数MTF构造Wiener滤波器M_W[u，v]，直流分量归一化后得到Wiener滤波器M′_W[u，v]，用其对过采样图像g_os[m，n]进行Wiener滤波，复原出重构图像完成微扫描图像的重构；具体过程为：

步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H[u，v]，其中[u，v]为频域坐标，u为频域上的水平坐标，v为频域上的竖直坐标：

假设f(x，y)为经过光学系统成像在焦平面探测器上的模拟图像，g[m，n]为经焦平面探测器积分和采样后得到的退化图像，则

$\begin{array}{c}g[m,n]=\frac{1}{\mathrm{pq}}{\∫}_{\mathrm{nb}-q/2}^{\mathrm{nb}+q/2}{\∫}_{\mathrm{ma}-p/2}^{\mathrm{ma}+p/2}f(x,y)\mathrm{dxdy}\\ =\frac{1}{\mathrm{pq}}{\∫}_{-q/2}^{q/2}{\∫}_{-p/2}^{p/2}f(x+\mathrm{ma},y+\mathrm{nb})\mathrm{dxdy}\end{array}m=\mathrm{1,2},...,M,n=\mathrm{1,2},...,N---\left(1\right)$

式中，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运算；

和

分别为向下取整和向上取整运算；sinc函数的定义为sinc(x)＝sin(πx)/(πx)；

由式(2)可以获得焦平面探测器的积分退化函数H[u，v]

$H[u,v]=\sin c\left(\frac{u}{M}\mathrm{\α}\right)\sin c\left(\frac{v}{N}\mathrm{\β}\right),---\left(3\right)$

步骤302、基于式(3)所表示的焦平面探测器的积分退化函数H[u，v]，构建Wiener滤波器M_W[u，v]

$M_W[u,v]=\frac{H^{*}[u,v]}{H^2[u,v]+\mathrm{\Γ}}---\left(4\right)$

式中，H[u，v]为积分退化函数；H^*[u，v]为H[u，v]的复共轭；Γ为常数，它并不随频率变化，而随着占空比、噪声等级以及图像内容变化；将M_W[u，v]用直流分量归一化后得到归一化后Wiener滤波器M′_W[u，v]

$M_W^{\′}[u,v]=\frac{M_W[u,v]}{M_W\left[\mathrm{0,0}\right]}---\left(5\right)$

式中，M_W[0，0]是M_W[u，v]在u＝0，v＝0情况下的值，即M_W[u，v]的直流分量；

步骤303、利用直流分量归一化后的Wiener滤波器M′_W[u，v]，对过采样图像g_os[m，n]进行滤波，获得重建图像

式中，

表示DFT运算，表示DFT逆运算。

本发明的有益效果：

本发明在基于焦平面探测器调制传递函数MTF构建Wiener滤波器M_W[u，v]的过程中，考虑了焦平面探测器的占空比，克服了现有亚像元重构算法中假设焦平面探测器的占空比为100％的限制；另外利用焦平面探测器调制传递函数MTF构建的Wiener滤波器进行滤波，抑制了焦平面探测器的积分退化效应，在帧间过采样图像重构技术的基础上进一步提高图像质量，从面进一步提高微扫描成像系统的成像性能。

附图说明

图1为本发明的流程图；

具体实施方式

为了更好地理解本发明的技术方案，下面结合附图及具体实施例对本发明做进一步详细描述。

如图1，下面以2×2微扫描图像重构为例，对本发明做进一步说明：

步骤一、利用基于DFT的亚像素图像配准算法，计算2、3、4帧微扫描图像相对于第1帧图像的帧间微位移量；

步骤二、在步骤一获得的帧间位移量的基础上，利用基于帧间差分的过采样重构算法，将4帧微扫描重建出一帧过采样图像g_os[m，n]；

步骤三、基于焦平面探测器调制传递函数MTF构建Wiener滤波器M_W[u，v]，经直流分量归一化后，对过采样图像g_os[m，n]进行Wiener滤波，抑制图像中因焦平面探测器有限像元尺寸引起的积分退化效应，复原出更高图像质量的高分辨率图像。

步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H[u，v]。

假设f(x，y)为经过光学系统成像在焦平面探测器上的模拟图像，gpm，n]为经焦平面探测器积分和采样后得到的退化图像，则

$\begin{array}{c}g[m,n]=\frac{1}{\mathrm{pq}}{\∫}_{\mathrm{nb}-q/2}^{\mathrm{nb}+q/2}{\∫}_{\mathrm{ma}-p/2}^{\mathrm{ma}+p/2}f(x,y)\mathrm{dxdy}\\ =\frac{1}{\mathrm{pq}}{\∫}_{-q/2}^{q/2}{\∫}_{-p/2}^{p/2}f(x+\mathrm{ma},y+\mathrm{nb})\mathrm{dxdy}\end{array}m=\mathrm{1,2},...,M,n=\mathrm{1,2},...,N---\left(1\right)$

式中，a和b为焦平面探测器的水平和竖直像元间隔，p×q为像元尺寸，M×N为焦平面探测器像面规格，(x，y)是模拟图像坐标系中的坐标，[m，n]是离散图像坐标系中的坐标。

对式(1)进行DFT，得

式中，F[u，v]和G[u，v]分别表示

和g[m，n]的DFT，[u，v]是频域坐标；α＝p/a和β＝q/b分别表示焦平面探测器在水平和竖直方向的占空比；

表示DFT运算；

和

分别为向下取整和向上取整运算；sinc函数的定义为sinc(x)＝sin(πx)/(πx)。

由式(2)可以获得焦平面探测器的积分退化函数H[u，v]

$H[u,v]=\sin c\left(\frac{u}{M}\mathrm{\α}\right)\sin c\left(\frac{v}{N}\mathrm{\β}\right),---\left(3\right)$

可以看出式(3)就是焦平面探测器的MTF。

步骤302、基于式(3)描述的积分退化模型，构建Wiener滤波器M_W[u，v]

$M_W[u,v]=\frac{H^{*}[u,v]}{H^2[u,v]+\mathrm{\Γ}}---\left(4\right)$

式中，H[u，v]为退化函数；H^*[u，v]为H[u，v]的复共轭；Γ为常数，它并不随频率变化，而随着占空比、噪声等级以及图像内容变化；[u，v]是频域坐标。将M_W[u，v]用直流分量归化后得到归化后Wiener滤波器M′_W[u，v]

$M_W^{\′}=[u,v]=\frac{M_W[u,v]}{M_W\left[\mathrm{0,0}\right]}---\left(5\right)$

式中，M_W[0，0]是M_W[u，v]在u＝0，v＝0情况下的值，即M_W[u，v]的直流分量。

步骤303、利用直流分量归一化后的Wiener滤波器M′_W[u，v]，对过采样图像g_os[m，n]进行滤波，获得具有更好图像质量的重建图像

式中，

表示DFT运算，

表示DFT逆运算。由于过采样过程将采样间隔缩小1倍，Wiener滤波器M′_W[u，v]中的占空比应提高1倍，即α_os＝2α，β_os＝2β。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于焦平面探测器MTF的微扫描图像重构方法，其特征在于，

步骤一：利用基于离散傅里叶变换DFT的亚像素图像配准方法，计算同一场景多帧图像间的帧间位移量；

步骤二：在步骤一获得的帧间位移量的基础上，利用帧间差分过采样方法，将多帧图像重构成一帧过采样图像g_os[m，n]；

步骤三：基于焦平面探测器的调制传递函数MTF构造Wiener滤波器M_W[u，v]，直流分量归一化后得到Wiener滤波器M′_W[u，v]，用其对过采样图像g_os[m，n]进行Wiener滤波，复原出重构图像

完成微扫描图像的重构；具体过程为：

步骤301、建立焦平面探测器的积分退化模型，获得积分退化函数H[u，v]，其中[u，v]为频域坐标，u为频域上的水平坐标，v为频域上的竖直坐标：

假设f(x，y)为经过光学系统成像在焦平面探测器上的模拟图像，g[m，n]为经焦平面探测器积分和采样后得到的退化图像，则

$\begin{array}{c}g[m,n]=\frac{1}{\mathrm{pq}}{\∫}_{\mathrm{nb}-q/2}^{\mathrm{nb}+q/2}{\∫}_{\mathrm{ma}-p/2}^{\mathrm{ma}+p/2}f(x,y)\mathrm{dxdy}\\ =\frac{1}{\mathrm{pq}}{\∫}_{-q/2}^{q/2}{\∫}_{-p/2}^{p/2}f(x+\mathrm{ma},y+\mathrm{nb})\mathrm{dxdy}\end{array}m=\mathrm{1,2},...,M,n=\mathrm{1,2},...,N---\left(1\right)$

式中，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]=\sin c\left(\frac{u}{M}\mathrm{\α}\right)\sin c\left(\frac{v}{N}\mathrm{\β}\right),---\left(3\right)$

步骤302、基于式(3)所表示的焦平面探测器的积分退化函数H[u，v]，构建Wiener滤波器MW[u，v]

$M_W[u,v]=\frac{H^{*}[u,v]}{H^2[u,v]+\mathrm{\Γ}}---\left(4\right)$

式中，H[u，v]为积分退化函数；H^*[u，v]为H[u，v]的复共轭；Γ为常数，它并不随频率变化，而随着占空比、噪声等级以及图像内容变化；将M_W[u，v]用直流分量归一化后得到归一化后Wiener滤波器M′_W[u，v]

$M_W^{\′}[u,v]=\frac{M_W[u,v]}{M_W\left[\mathrm{0,0}\right]}---\left(5\right)$

式中，M_W[0，0]是M_W[u，v]在u＝0，v＝0情况下的值，即M_W[u，v]的直流分量；

步骤303、利用直流分量归一化后的Wiener滤波器M′_W[u，v]，对过采样图像g_os[m，n]进行滤波，获得重建图像

式中，

表示DFT运算，

表示DFT逆运算。

Images