发明名称 单幅图像复原的解析迭代快速频谱外推方法
摘要 本发明公开了一种单幅图像复原的解析迭代快速频谱外推方法,包括迭代系统初始化过程和解析迭代过程。其中:初始化过程计算降质图像频谱、设置初始迭代解、系统初始参数;解析迭代过程通过不断精细的梯度信息阈值收缩估计图像的水平和垂直方向梯度信息,通过解析迭代频谱和梯度信息频谱解析预测得到预测频谱,通过图像频谱合成校正将降质图像频谱和预测频谱解析合成校正得到高分辨率图像。本发明方法基于快速傅立叶变换技术,通过截止频率以下频谱,外推截止频率以上频谱,以很小的复杂度达到图像快速复原。去模糊和抑制噪声能力明显优于常规方法,信噪比得到显著提高。该方法非常便于利用快速傅立叶变换DSP芯片组成的图像处理硬件实现。
申请公布号 CN101540043B 申请公布日期 2012.05.23
申请号 CN200910031270.8 申请日期 2009.04.30
申请人 南京理工大学 发明人 肖亮;黄丽丽;韦志辉;吴泽彬
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 唐代盛
主权项 1.一种单幅图像复原的解析迭代快速频谱外推方法,其特征在于包括迭代系统初始化过程和解析迭代过程,即1.1 所述的迭代系统初始化过程包括:输入一幅待复原的M×N大小的降质图像u<sub>0</sub>,计算降质图像频谱F(u<sup>0</sup>);设置初始迭代解u<sup>(0)</sup>=u<sub>0</sub>,初始频谱校正图像w<sup>(0)</sup>=u<sub>0</sub>;设置迭代系统参数初始值:梯度信息阈值α<sup>(0)</sup>>0,预测贡献度控制参数β<sup>(0)</sup>>0,合成贡献度参数γ<sup>(0)</sup>>1,图像频谱误差预设阈值ε≤10<sup>-4</sup>;并设置空间不变成像系统的系统调制传递函数F(h);1.2 所述的解析迭代过程包括如下步骤:步骤1:首先利用第k-1次的频谱校正图像w<sup>(k-1)</sup>,其中k≥1,分别计算其在水平方向和垂直方向的向前差分,并利用第k-1次的梯度信息阈值α<sup>(k-1)</sup>作阈值收缩判断,得到第k次的水平方向差分估计<img file="FSB00000519129900011.GIF" wi="100" he="63" />和垂直方向差分估计<img file="FSB00000519129900012.GIF" wi="130" he="62" />步骤2:将第k次的水平方向差分估计<img file="FSB00000519129900013.GIF" wi="101" he="63" />作快速傅立叶变换(FFT)得到水平方向高频频谱<img file="FSB00000519129900014.GIF" wi="209" he="89" />同样第k次的垂直方向差分估计<img file="FSB00000519129900015.GIF" wi="105" he="63" />作快速傅立叶变换(FFT)得到垂直方向高频频谱<img file="FSB00000519129900016.GIF" wi="210" he="88" />然后,利用计算得到的水平方向高频频谱<img file="FSB00000519129900017.GIF" wi="211" he="93" />垂直方向高频频谱<img file="FSB00000519129900018.GIF" wi="180" he="67" />以及第k-1次合成复原图像频谱F(u<sup>(k-1)</sup>)进行频谱预测得到第k次的预测频谱F(w<sup>(k)</sup>);步骤3:利用图像频谱合成校正,将第k次的预测频谱F(w<sup>k</sup>)和降质图像频谱F(u<sup>0</sup>)结合光学系统调制传递函数F(h)及其复共轭F(h)<sup>*</sup>通过矩阵元素逐点相乘、相加和矩阵元素逐点相除的方式进行合成校正,得到合成校正频谱F(u<sup>(k)</sup>);步骤4:采用2范数的平方<img file="FSB00000519129900019.GIF" wi="56" he="75" />计算第k次校正图像频谱F(u<sup>(k)</sup>)和第K-1次校正图像频谱F(u<sup>(k-1)</sup>)的误差,如果<img file="FSB000005191299000110.GIF" wi="632" he="68" />则转至步骤5,否则转至步骤6;步骤5:表明频谱外推没有收敛,则进行参数自适应校正,调整系统参数,转至步骤1,k=k+1;步骤6:判断是否达到收敛条件<img file="FSB00000519129900021.GIF" wi="634" he="67" />如果到达则算法结束,输出复原图像,对合成校正图像F(u<sup>(k)</sup>)进行快速逆傅立叶变换得到复原图像。
地址 210094 江苏省南京市孝陵卫200号