[发明专利]基于多尺度几何分析的图像超分辨率重建方法

权利要求书

1.一种基于多尺度几何分析的图像超分辨率重建方法，其特征在于，具体按照以下步骤实施：

步骤1、利用降质模型生成低分辨率图像

从观测图像中选取连续三帧普通分辨率图像I₁、I₂和I₃，利用下列模型进行退化处理，得到相应的低分辨率图像T₁、T₂和T₃：

T_k＝D·B·M·I_k+n＝H·I_k+n，

式中，k＝1,2,3，D、B、M、n分别为下采样算子、混淆算子、运动模糊算子、噪声算子；

步骤2、图像配准

对于待配准的图像T_k，k＝1,2,3，选其中信噪比较大的一帧图像T₁作为参考帧，利用SIFT算法提取图像特征点，设T_k的特征点集合分别为：

P1={p1,1,p1,2,...,p1,m1},]]>

P2={p2,1,p2,2,...,p2,m2},]]>

P3={p3,1,p3,2,...,p3,m3},]]>

其中，m_k为对应于T_k的特征点数目，k＝1,2,3，j＝1,2,...,max(m₁,m₂,m₃)，

分别计算集合P₁与P₂、P₁与P₃中的特征点间的欧氏距离：

d(p1,j1,p2,j2)=Σi=1128(p1,j1i-p2,j2i)2,]]>j₁＝1,2,...,m₁，j₂＝1,2,...,m₂

d(p1,j1,p3,j3)=Σi=1128(p1,j1i-p3,j3i)2,]]>

其中，j₁＝1,2,...,m₁，j₂＝1,2,...,m₂，j₃＝1,2,...,m₃，

找出P₁与P₂中特征点间的最小欧氏距离d₁与次小欧氏距离d₂，定义距离比率为r＝d1/d2，定义配准准则为：

其中η为由实验得到的阈值，如果距离比率大于阈值即认为匹配成功，否则认为匹配失败；对T₃做同样的处理，得到配准后图像分别为：

T₂′＝M₁·T₂，

T₃′＝M₂·T₃，

其中，M₁、M₂均为运动补偿矩阵；

步骤3、图像融合

对配准后的图像进行Contourlet变换，即C_k＝C(T_k′)，其中，C(·)为轮廓波运算，假设匹配后待重建图像T_k′经过Contourlet变换后的低频系数记为如下式所示：

Xkg=xk,1,1gxk,1,2g...xk,1,jgxk,2,1gxk,2,2g...xk,2,jg............xk,i,1gxk,i,2g...xk,i,jg,]]>

其中k＝2,3，低频系数的个数为l＝i*j，记第s个尺度第τ个方向子带的高频系数为共有高频系数h＝p*q个，

Yks,τ=yk,1,1s,τyk,1,2s,τ...yk,1,qs,τyk,2,1s,τyk,2,2s,τ...yk,2,qs,τ............yk,p,1s,τyk,p,2s,τ...yk,p,qs,τ,s=1,2,...,g]]>

对于Contourlet变换后所得的低频系数，采用加权平均法进行处理，得到逆变换的低频系数，第k帧匹配后待重建图像T_k′经过Contourlet变换后的低频系数为ω_k为权重系数，且满足则最终选择作为Contourlet逆变换的低频系数为

对于Contourlet变换后所得的高频部分按以下原则进行融合：

分别比较相应的每个尺度每个方向子带上所对应高频系数方差的大小，取方差最小的系数作为相应尺度相应方向逆变换的高频系数，即采用表示第s个尺度第τ个方向子带高频系数的平均值，该平均值表达式为：

Yks,τ‾=1p·qΣi=1pΣj=1qYk,i,js,τ,]]>

采用表示第s个尺度第τ个方向子带高频系数的方差，该方差表达式为：Vks,τ=1p·qΣi=1pΣj=1q(Yk,i,js,τ-Yks,τ‾)2,]]>

即为第s个尺度第τ个方向子带Contourlet逆变换的高频系数，令：

G=G1,1G21,1000...0G1,2G2,2G3,2G22,200..................G1,gG2,gG3,gG4,gG5,gG2g,g,]]>

其中G^s,τ表示最终选择作为逆变换的第s个尺度第τ个方向子带的高频系数，G为各方向各尺度系数矩阵，其中行表示方向子带的方向索引，列表示尺度索引，

用低频系数L和每个尺度每个方向的高频系数G进行Contourlet逆变换，即I'＝C^-1(L,G)，C^-1(·)为轮廓波逆变换运算，I'即为融合图像；

步骤4、获得重建图像

对得到的融合图像I'进行双三次差值运算得到I''，对I''进行双边滤波处理，得到重建图像II。

2.根据权利要求1所述的基于多尺度几何分析的图像超分辨率重建方法，其特征在于：所述的步骤1中，利用点扩散函数PSF产生运动模糊函数和下采样算子，该函数有两个参数theta和len，表示摄像物体逆时针方向以theta角度运动了len个像素，通过调整这两个参数来改变图像的模糊程度。