I=imread('D:/桌面/zxy.tif');
figure
subplot(1,3,1),imshow(I),title('原图像');
I=im2double(I);
F=fft2(I);
F=fftshift(F);
[M,N]=size(F);
[u,v]=meshgrid(1:M-1,1:N-1);
H=exp(-0.0025* ((u-M/2).^2+(v-N/2).^2).^(5/6));
FBlurred=F.*H;
IBlurred=ifftshift(FBlurred);
IBlurred=ifft2(IBlurred);
IBlurred=real(IBlurred);
figure
subplot(1,3,2),imshow(IBlurred,[]),title('无噪退化图像');
对于此运算,数组的大小不兼容,怎么解决
因为H的长宽比F都小1,而且H需要转置一下
I=imread('D:/桌面/zxy.tif');
I=im2gray(I);
figure(1)
subplot(1,3,1),imshow(I),title('原图像');
I=im2double(I);
F=fft2(I);
F=fftshift(F);
[M,N]=size(F);
[u,v]=meshgrid(1:M,1:N);
H=exp(-0.0025* ((u-M/2).^2+(v-N/2).^2).^(5/6));
FBlurred=F.*H';
IBlurred=ifftshift(FBlurred);
IBlurred=ifft2(IBlurred);
IBlurred=real(IBlurred);
figure(1)
subplot(1,3,2),imshow(IBlurred,[]),title('无噪退化图像');