weizhi.zip_blur image_距离 角度_距离模糊_运动 模糊_运动模糊

%将彩色图像转为灰度图： filename='a.jpg'; imgRgb=imread(filename); % 读入一幅彩色图像 imshow(imgRgb); % 显示彩色图像 imgGray=rgb2gray(imgRgb); % 转为灰度图像 figure(1) % 打开一个新的窗口显示灰度图像 imshow(imgGray);% 显示转化后的灰度图像 title('灰度图像'); imwrite(imgGray, 'gray.jpg'); % 将灰度图像保存到图像文件 %当PSF未知时，我们首先要估计点扩散函数： 计算频谱： I=imread('gray.jpg'); figure(2) imshow(I); title('原图像'); K=fft2(I);%傅里叶变换 M=fftshift(K);%直流分量移动到频谱中心 N=abs(M);%计算频谱幅值 P=(N-min(min(N)))/(max(max(N))-min(min(N)))*255;%归一化 figure(3) imshow(P); title('模糊图像傅里叶变换频谱'); %计算模糊距离： I=imread('gray.jpg'); f=im2double(I); h=fspecial('sobel'); %sobel算子 J=conv2(f,h,'same');%sobel算子微分 IP=abs(fft2(J)); %图像能量谱密度 S=fftshift(real(ifft2(IP))); figure(4) plot(S); title('自相关函数');