matlab starcat_维纳滤波和编码曝光PSF去除运动模糊【matlab】

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%% 作者：WWC %%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% 维纳滤波

clc;

clear;

close all;

%% 读取图片

I=imread('飞机.png');

I2=imfilter(I,[0 0 0 0 0 0 0 0 5/25 4/25 4/25 3/25 3/25 2/25 2/25 1/25 1/25]);

I2=imnoise(I2,'gaussian',0,0.001);

I1=imfilter(I2,[1/9 1/9 1/9;1/9 1/9 1/9;1/9 1/9 1/9;]);

I_hsv=rgb2hsv(I1);

I=I_hsv(:,:,3)*255;

%% 频谱 估计运动模糊方向(光斑的角度与真实图像运动模糊距离垂直)

img_fft=fftshift(fft2(I));

N=abs(img_fft);

P=(N-min(min(N)))/(max(max(N))-min(min(N)))*225;

figure;imshow(P);

%% 估计运动模糊长度(两个最小点之间的距离为近似运动模糊距离)

h=fspecial('sobel');

img_double=double(I);

J=conv2(img_double,h,'same');

IP=abs(fft2(J));

S