[发明专利]POCS超分辨率图像重建时边缘Halo效应的抑制方法

权利要求书

1.POCS超分辨率图像重建时边缘Halo效应的抑制方法，其特征在于，是在计算机中依次按以下步骤实现的：

步骤（1），计算机初始化：

设定：图像读取模块、高分辨率初始估计值获取模块、低分辨率图像配准模块以及凸集投影POCS方法重建模块；

步骤（2），依次按以下步骤抑制POCS超分辨率图像重建时边缘Halo效应，边缘Halo效应是指边缘晕圈效应：

步骤（2.1），所述图像读取模块从图像库读取在同一场景下的低分辨率图像序列g_i，1≤i≤n，i为低分辨率图像的序号，n为低分辨率图像的总数，n∈[6,36]，根据n确定迭代次数L∈[6,15]，所述低分辨率图像的分辨率为M×N，M∈[64,256]，N∈[64,256]，相邻低分辨率图像间有重叠部分，具有公共坐标系，用(x,y)表示坐标点；

步骤（2.2），所述高分辨率初始估计值获取模块，包含低分辨率参考帧获取子模块、小波分解子模块、双三次插值子模块以及小波逆变换子模块，对应的高分辨率图像的大小为qM×qN，q表示分辨率提高的倍数，q∈[2,4]，其中：

低分辨率参考帧获取子模块，任选所述低分辨率图像序列其中一帧低分辨率图像作为低分辨率参考帧g₀，并保存一份低分辨率参考帧；

小波分解子模块，对所述低分辨率参考帧进行q层离散小波分解变换，得到低频子图、水平方向的高频子图、垂直方向的高频子图、对角线方向的高频子图，所有子图的大小均为M/q×N/q；

双三次插值子模块，采用Matlab库中的双三次插值法分别对三个所述的高频子图和进行外推，得到分辨率都为M×N的高频子图g^H,g^V和g^D；

小波逆变换子模块，把所述低分辨率参考图像g₀作低频子图，三个M×N的所述图像g^H,g^V,g^D作高频子图进行小波逆变换，获得qM×qN的高分辨率图像初始估计值f₀，[f₀]＝IDWT[g₀,g^H,g^V,g^D,q]；

步骤（2.3），所述低分辨率图像配准模块，包括：预处理子模块、运动估计子模块以及运动向量补偿子模块，其中：

预处理子模块，利用低通滤波器F对包括保存的低分辨率参考帧在内的所有低分辨率图像进行滤波，F＝1/16×[1,2,1;2,4,2;1,2,1]，获得初步处理后的低分辨率图像；

运动估计子模块，按下式逐次求取所述各低分辨率观测帧g_i(x,y)与低分辨率参考帧g₀(x,y)之间的运动矢量Z＝[k₁,k₂,k₃,k4_]：

Z＝G^-1V，其中：

V=ΣR1(g0-gt)ΣR(g0-gt)Σ∂gt∂x(g0-gt)Σ∂gt∂y(g0-gt),]]>R=y∂gt∂x-x∂gt∂y,]]>R1=x∂gt∂x+y∂gt∂y,]]>

G=ΣR12ΣRR1ΣR1∂gt∂xΣR1∂gt∂yΣRR1ΣR2ΣR∂gt∂xΣR∂gt∂yΣR1∂gt∂xΣR∂gt∂xΣ(∂gt∂x)2Σ(∂gt∂x·∂gt∂y)ΣR1∂gt∂yΣR∂gt∂yΣ(∂gt∂x·∂gt∂y)Σ(∂gt∂y)2]]>

通过采取迭代方法逐步逼近真实值，s为迭代次数，当迭代误差满足一定精度Th≤10^-3时终止迭代；

运动向量补偿子模块，设定高分辨率图像补偿前的坐标为(x',y')，高分辨率图像补偿后的坐标为(x,y)，按下式把计算得到的所述运动矢量Z＝[k₁,k₂,k₃,k₄]补偿到高分辨率初始估计f₀(x,y)中，所述f₀是f₀(x,y)的简略表达式：

x=x′+k1x′+k2y′+k3y=y′+k1y′-k2x′+k4]]>

步骤（2.4），所述凸集投影POCS方法重建模块，包含：点扩散函数PSF获取子模块、松弛投影参数自适应调整子模块以及迭代投影修正子模块，依次按以下步骤降低重建结果中边缘Halo效应：

点扩散函数PSF获取子模块，依次按以下步骤求出具有边缘保持特性的点扩散函数PSF，其特征在于：

按下式计算所述高分辨率初始估计值域f₀(x，y)中各像素与支撑域中心像素之间的空域中的欧氏距离，用h(x,y)表示：

h(x,y)=exp{-(x-m)2+(y-n)22δs2},]]>其中，

δ_s为点扩散函数PSF的标准差，δ_s∈[0.5,2]，(m,n)为位于各相邻像素之间的5×5像素支撑域内中心像素的坐标，

按下式计算所述高分辨率初始估计值域f₀(x，y)中各像素之间的灰度关系

wδr(|f0(x,y)-f0(m,n)|)=exp{-|f0(x,y)-f0(m,n)|22δr2},]]>其中，

f₀(x,y),f₀(m,n)分别为坐标(x,y),(m,n)处的灰度值，δ_r为标准差，δ_r∈[0.8,15]；

按下式求出考虑了所述高分辨率初始估计值域中各像素之间的灰度关系的具有边缘保持特性的点扩散函数

h~(x,y)=U·h(x,y)·wδr(|f0(x,y)-f0(m,n)|),]]>

U为归一化系数，U∈(0.001,1)；

松弛投影参数自适应调整子模块，依次按以下步骤根据低分辨率观测帧和低分辨率参考帧之间的相关性自适应调整松弛投影参数，获得可变长的松弛投影参数：

按下式计算各低分辨率观测帧与低分辨率参考帧之间的相关性Rl：

Rl=Σx,y[gi(x,y)g0(x,y)]Σx,y[gi(x,y)]2Σx,y[g0(x,y)]2,]]>0＜Rl＜1，

g_i(x,y)为各低分辨率当前观测帧，g₀(x,y)为低分辨率参考帧，

设定：松弛投影参数λ，λ＝2Rl，所述λ是一个可变长的松弛投影参数，反映了低分辨率当前观测帧参与投影的程度，该λ值与Rl成正比；

迭代投影修正子模块，依次按以下步骤对所述高分辨率初始估计值域f₀(x,y)进行凸集迭代投影修正，以便通过设定的迭代次数获得高分辨率图像：

已知：先验信息及其先验界限值δ_i，所述先验信息包括正定、能量有限、数据一致性光滑在内的任何一种信息，根据数据一致性信息定义凸集C_m，n，i，(m,n)为低分辨率图像g_i(x,y)中坐标点(x,y)对应于高分辨率估计值中的坐标，i为低分辨率图像观测帧的序号,

设定高分辨图像迭代投影过程中的估计值域l∈[0,L)，L为迭代的总次数，L∈[6,15]，则任意一个高分辨图像迭代投影过程中的估计值域到凸集C_m，n，i上的投影为：

P(x,y)[f^l(m,n)]=f^l(m,n)+(λ·rif^l(m,n)(x,y)-δ0)·h~i(x,y;m,n),rif^l(m,n)(x,y)>δi0,|rif^l(m,n)(x,y)|≤δi(λ·rif^l(m,n)(x,y)+δ0)·h~i(x,y;m,n),rif^l(m,n)(x,y)<-δi,]]>

其中，δ_i＝cδ₀，c∈[2,5]，δ₀∈[0.1,3.8]，

Cm,n,i={f^l(m,n):|rif^l(m,n)(x,y)|≤δi},]]>

为残差函数，rif^l(m,n)(x,y)=gi(x,y)-Σm,nf^l(m,n)h~i(x,y;m,n),]]>

为高分辨率估计值中点(m,n)处灰度值与点(m,n)的5×5支撑域中各像素灰度值之间的约束关系，表达式为：

h~i(x,y;m,n)=U·h(x,y;m,n)·wδr(|f^l(x,y)-f^l(m,n)|),]]>U∈(0.001,1)

h(x,y;m,n)=exp{-(x-m)2+(y-n)22δs2},]]>δ_s∈[0.5，]，

wδr(|f^l(x,y)-f^l(m,n)|)=exp{-|f^l(x,y)-f^l(m,n)|22δr2},]]>δ_r∈[0.8,15]，

实质上是具有边缘保持特性的点扩散函数，通过把高分辨率图像估计值进行点扩散处理来获得低分辨率图像；

令由计算得到依次类推，逐次进行迭代；若未达到迭代次数，则把代入残差函数更新残差否则迭代结束，从中截断超出8比特图像取值范围[0，255]部分，得到重建的高分辨率图像。

2.根据权利要求1所述的POCS超分辨率图像重建时边缘Halo效应的抑制方法，其特征在于：选取低分辨率图像序列中的第一帧作为参考帧。