I=imread('roadmark_59.jpg');
F=fftshift(fft2(I));
[M,N]=size(F);
n=2; %对n赋初值
D0=120;
n1=floor(M/2); %对M/2进行取整
n2=floor(N/2); %对N/2进行取整
for u=1:M
for v=1:N
D(u,v)=sqrt((u-n1)^2+(v-n2)^2);
H=1*exp(-(D(u,v)/D0)^n); %BLPF滤波函数
F(u,v)=H*F(u,v); %BLPF滤波后的频域表示
end
end
F=ifftshift(F);
s=uint8(real(ifft2(F)));
subplot(1,2,1);imshow(I),title('原图像');
subplot(1,2,2);imshow(s),title('指数滤波后图像D0=10 ');