同态滤波及matlab实现(共1页).docx

clear allclcI=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);for i=1:M for j=1:N D(i,j)=(i-n1).2+(j-n2).2); H(i,j)=(rH-rL).*(exp(c*(-D(i,j)./(d02)+rL; %高斯同态滤波 endendI2=ifft2(H.*FI);%傅里叶逆变换I3=real(exp(I2);subplot(122),imshow(I3,);title(同态滤波增强后)