s101.m

简单的数字图像处理程序

% 同态滤波器  s101
clear;
m1=imread('2.jpg');
J1=rgb2gray(m1);
subplot(1,2,1),imshow(J1),title('原图');
f=double(J1);
c=1/log10(50);
[m,n]=size(f);
for i=1:m
    for j=1:n
        r(i,j)=10*log10(1+49*f(i,j));
    end
end
p=fft2(r);
g=fftshift(p);
[N1,N2]=size(g);
n=3;
d0=5;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
    for j=1:N2
        d=sqrt((i-n1)^2+(j-n2)^2);
        if d<2
            h=1/5;
        else
            h=1/(1+(d0/d)^(2*n));
        end
        result(i,j)=h*g(i,j);
    end
end
result=ifftshift(result);
x2=ifft2(result);
x3=exp(x2-2.4);
x4=uint8(real(x3));
subplot(1,2,2),imshow(x4),title('同态滤波处理后')