同态滤波及matlab实现

clearallclcI=imread('716.png');subplot(121),imshow(I);title('原图')I=double(rgb2gray(I));[M,N]=size(I);rL=0.3;rH=2.0;%可根据需要效果调整参数c=2;d0=10;I1=log(I+1);%取对数FI=fft2(I1);%傅里叶变换n1=floor(M/2);n2=floor(N/2);fori=1:Mforj=1:ND(i,j)=((i-n1).^2+(j-n2).^2);H(i,j)=(rH-rL).*(exp(c*(-D(i,j)./(d0^2))))+rL;%高斯同态滤波endendI2=ifft2(H.*FI);%傅里叶逆变换I3=real(exp(I2));subplot(122),imshow(I3,[]);title('同态滤波增强后')原图同态滤波增强后