Question

基于matlab运动模糊图像处理的源代码

求一段matlab程序，把下面这个运动模糊图像变清楚，最好求出运动模糊方向和运动模糊像素的个数，十万火急！

Answer 1

等一会的，我来帮你

华东师范大学？？？

tuxianghuanyuan('3.jpg', 80, 8, 0.02);

function tuxianghuanyuan(im, a, b, NSPR)

i = imread(im);

f = im2double(i);

PSF = fspecial('motion', a, b);

frest1 = deconvwnr(f, PSF, NSPR);

subplot(221),imshow(f); title('原图像');

subplot(222),imshow(frest1); title('维纳滤波处理后图像'); 

end

Answer 2

PSF = fspecial('motion',len,ang); %建立扩散子，其中len是模糊长度，ang是模糊角度
img2=deconvlucy(img,PSF,n); %用lucy-richardson方法复原图像，其中img是运动模糊图像，PSF是扩散子，n是迭代次数，img2是复原图像

追问

我是要把模糊的图像变清楚，不是把清楚的图像用PSF变模糊再用lucy-richardson
变清楚。现在就是这幅图的模糊长度和模糊角度不知道，这是我的毕业设计题目，不可能这么简单，我已经做十多天了，试过很多方法，我需要从头到尾完整的源代码，就是copy后放在matlab软件上能立即出结果的程序

我是要把模糊的图像变清楚，不是把清楚的图像用PSF变模糊再用lucy-richardson变清楚。现在就是这幅图的模糊长度和模糊角度不知道，这是我的毕业设计题目，不可能这么简单，我已经做十多天了，试过很多方法，我需要从头到尾完整的源代码，就是copy后放在matlab软件上能立即出结果的程序

追答

一般模糊图像的模糊距离和模糊角度不好判断，但也可以通过查看图像的频谱图和实部二维图来大致估计。

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

% 求角度

img_fft=fftshift(fft2(img_gray));

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(img_gray);

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

IP=abs(fft2(J));

S=fftshift(real(ifft2(IP)));

figure;plot(S);

上面的代码可以求出图像实部二维图并查看

运动距离可以用两个最小值之间的距离估计

一般运动模糊是无法精确得到角度和距离的，所以都是估计方法。希望可以帮到你。

Answer 3

可用于求模糊方向。。可参考网页链接    网页链接

%% 读入并显示图像
filename = '112.jpg';
I = imread(filename);
% figure
% imshow(uint8(I));
% title('原图');
%% 生成运动模糊图像
PSF = fspecial('motion',20,10);
g = imfilter(I, PSF, 'circular');
% figure
% imshow(uint8(g));
% title('运动模糊图');
%% 对运动模糊图像进行灰度化，并进行二维快速傅里叶变换，生成其频谱图
gb = rgb2gray(g);
I=im2double(gb);
F=fft2(I);
F=fftshift(F);
imshow(F);
F=abs(F);
t=log(F+1);
figure;
imshow(t,[]);
%% 用canny算子将压缩居中后的频谱图进行边缘检测，二值化
T = graythresh(t);
bw=edge(t, 'canny', T);
imshow(bw);
title('边缘检测后的图');
%% 对二值化后的频谱图进行radon变换
theta = 1:180;
[R,xp]= radon(bw, theta);
%% 计算出通过radon变换求出的模糊角度
MAX = max(max(R));
[m, n] = find(R == MAX);
[M,N] = size(F);
beita =atan(tan(n*pi/180)*M/N)*180/pi; %beita为模糊方向

Answer 4

我也想知道怎么求PSF 你知道怎么弄了吗？求讲解