基于傅里叶变换的运动模糊图像恢复算法matlab仿真

1.算法运行效果图预览

2.算法运行软件版本

matlab2022a

3.算法理论概述

运动模糊是由于摄像设备和目标之间的相对运动导致的一种图像降质现象。这种模糊会对图像的视觉质量和识别性能产生严重影响。为了解决这个问题，研究者们提出了各种图像恢复算法。其中，基于傅里叶变换的方法因其计算效率和恢复效果受到广泛关注。

3.1、傅里叶变换与图像恢复

傅里叶变换是一种将信号从时域转换到频域的数学工具。在图像处理中，傅里叶变换可以帮助我们分析图像的频率成分，进而实现去噪、压缩等操作。对于运动模糊图像恢复，傅里叶变换的主要作用是在频域中对模糊进行建模和分析。

i 是虚数单位，p 和 j 是值范围从 0 到 m–1 的索引，q 和 k 是值范围从 0 到 n–1 的索引。在此公式中，X 和 Y 的索引平移 1 位，以反映 MATLAB® 中的矩阵索引。计算 X 的二维傅里叶变换等同于首先计算 X 每列的一维变换，然后获取每行结果的一维变换。换言之，命令 fft2(X) 等同于 Y = fft(fft(X).').'。

3.2、基于傅里叶变换的运动模糊图像恢复算法原理

假设我们有一个运动模糊图像g(x, y)和一个点扩散函数h(x, y)，那么模糊过程可以表示为卷积操作：

g(x, y) = f(x, y) * h(x, y) + n(x, y) （1）

其中f(x, y)是原始清晰图像，n(x, y)是加性噪声。我们的目标是从g(x, y)中恢复出f(x, y)。

在频域中，卷积操作变为乘法操作，因此我们对等式（1）两边进行傅里叶变换：

G(u, v) = F(u, v) * H(u, v) + N(u, v) （2）

其中G、F、H和N分别是g、f、h和n的傅里叶变换。通过解这个等式，我们可以得到原始清晰图像的估计值。具体步骤如下：

对模糊图像g进行傅里叶变换得到G。

估计点扩散函数h的傅里叶变换H。这通常通过分析模糊图像的频谱特性来实现。

根据等式（2）计算原始清晰图像的傅里叶变换估计值F'。

对F'进行逆傅里叶变换得到恢复的清晰图像f'。

4.部分核心程序

%获取角度
img_fft2 = img_fft2/max(max(img_fft2));% 归一化处理  
[Rr,Cc]  = size(img_fft2);% 获取频谱尺寸  
img_x    = zeros(1,Rr*Cc); % 初始化x坐标数组  
img_y    = zeros(1,Rr*Cc);% 初始化y坐标数组  
idx      = 0;
for i1 = 1:Rr% 遍历频谱的每一个像素  
    for j1 = 1:Cc
        if img_fft2(i1,j1) > 0.0005% 如果像素值大于阈值  
           idx        = idx+1;    % 计数器加1  
           img_x(idx) = i1+1-Rr/2;% 计算x坐标  
           img_y(idx) = j1+1-Cc/2;% 计算y坐标
        end
    end
end
  
% 对x和y进行线性拟合，获取拟合参数p和统计信息s 
[p,s] = polyfit(img_x,img_y,1);
% 如果拟合的残差平方和大于4500  
if s.normr>4500
    % 交换x和y再次进行拟合  
    [p,s] = polyfit(img_y,img_x,1);
    % 计算角度（弧度转角度）  
    Angles= atan(p(1))/pi*180;
    % 调整角度值 
    Angles= 90-Angles;
    % 如果角度小于0度  % 加上180度
    if Angles<0
       Angles = Angles+180;
    end
    % 使用拟合参数计算y的拟合值  
    y_fit = polyval(p,img_x);
    subplot(132)
    plot(-img_y,img_x,'g.',img_x,-y_fit,'b');
else
    Angles = atan(p(1))/pi*180;% 计算角度（弧度转角度） 
    if Angles < 0% 如果角度小于0度 
       Angles = Angles +180;% 加上180度
    end
    y_fit = polyval(p,img_x);%求拟合后的y值;
    subplot(132)
    plot(-img_y,img_x,'g.',-y_fit,img_x,'b');
end