[发明专利]结合倒易晶胞的平移不变小波遥感图像复原方法

权利要求书

1.结合倒易晶胞的平移不变小波遥感图像复原方法，其特征在于包括如下步骤：

(1)通过图像退化公式得到退化后的图像g：

g＝Δ_T·(H*f)+n

其中：f为进入图像处理系统前的理想图像；H为系统退化函数卷积核；n为噪声；T为传感器阵列的几何结构，Δ_T代表基于传感器阵列几何形状的采样单元；

(2)对退化后的图像g利用倒易晶胞进行频域内的图像复原，以去除混迭和噪声，具体方法如下：

为Voronoi倒易晶胞，满足如下凸集条件：

DVor*=Vor(Γ*)=={γ*∈Rs:∀γ0*∈Γ*,γ0*≠0⇒||γ*||<||γ*+γ0*||}]]>

其中：Г^*为对偶网格，为Г^*的覆盖，R^s为覆盖的半径，γ^*为倒易晶胞内的任意一个点到中心的距离，为倒易晶胞内的任意一个不等于0的数；

在倒易晶胞内进行逆滤波器卷积核的频域响应为：

k^=MTF‾P|MTF|2(P+αa2+βb2)]]>

=(MTFsenMTFmovMTFopt)‾P|MTFsenMTFmovMTFopt|2(P+αa2+βb2)]]>

=(4c2sinc(πξc)sinc(πηc)e-β′c|η|sinc(κπ<ξ,v>)e-α′cu2+v2)‾P|4c2sinc(πξc)sinc(πηc)e-β′c|η|sinc(κπ<ξ,v>)e-α′cu2+v2|2(P+αa2+βb2)]]>

其中：c为传感器的大小，P为图像信号的功率谱密度PSD，α′为光学系统滤波参数，β′为相邻传感器之间的传导率，a为混迭阈值，b为噪声阈值，α和β为调整因子，当α＝β＝1时，就是Weiner滤波，ξ为传感器的单位像元的长，η为传感器的单位像元的宽，κ为卫星平台在v方向上移动的距离，v为卫星平台的飞行方向，u为与v的垂直方向，

MTF为卷积核H的傅里叶变换，为MTF的共轭，MTF＝MTF_senMTF_movMTF_opt，MTF_sen为探测器MTF，MTF_mov为卫星平台运动MTF，MTF_opt为光学系统MTF；

在频域内采用维纳滤波得到复原后的退化图像g′，

g′＝k*g，其中k＝argminE{||·||²}表示约束最小；

(3)对倒易晶胞在频域内复原后的退化图像g′进行平移不变小波分解，并对分解后的小波系数进行软阈值去噪和去卷积，具体过程如下：

g→′=φjAj+Σd∈AΣj≥JψjdBjd]]>

其中：为退化图像g′在平移不变小波域的变换，φ_j为尺度函数，为二维平移不变小波的尺度函数矩阵，φ_jx为x方向的尺度函数，φ_jy为y方向的尺度函数，A_j为尺度系数，为小波函数，为小波系数，d为平移不变小波分解的6个方向，j为分解的层数；表示平移不变小波在6个方向上的高频子带总数；

a、对小波函数在6个方向上进行分解，得到下式：

6个方向上的二维平移不变小波函数为：

b、对小波系数进行紧缩软阈值去噪，得到去噪后的小波系数

Bjs=sign(Bjd)(|Bjd|-t),]]>|Bjd|≥t]]>

其中：t＝3σ_j以确保99％的高置信度，用以去除白噪声，σ_j为j尺度下的噪声方差；

c、对去噪后的小波系数去卷积，得到去卷积后的小波系数

B→jw=1Bj|Bjs|2|Bjs|2+λσj2]]>

其中：B_j为j尺度上的小波系数；

λ为噪声的权重调整因子；

(4)进行平移不变小波重构，得到复原后的图像

g→′′=φjAj+ΣψjdBjw.]]>

2.根据权利要求1所述的结合倒易晶胞的平移不变小波遥感图像复原方法，其特征在于：所述步骤(3)中分解的层数j为4。