|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
网上找的一段同态滤波的程序,中间有一段看不懂
clear all;
I = imread('rice.png');
[M N] = size(I);
subplot(1,2,1),imshow(I);
title('原图');
T = double(I);
L = log(T);
F = fft2(L);
for i = 1:M
for j = 1:N
D(i,j) = ((i-M/2)^2+(j-N/2)^2);
end
end
c = 1.05; %锐化参数,可调
Do = max(M,N);%本应该是方差,我简化了一下
H = (2.0-0.5)*(1 - exp(c*(-D/(Do^2)))) + 0.5;
F = F .* H;
F = ifft2(F);
Y = exp(F);
subplot(1,2,2),imshow(uint8(Y));
title('处理后的图');
运行成功,效果还不错,但就是这段程序不是很懂,不知道他的思路。
for i = 1:M
for j = 1:N
D(i,j) = ((i-M/2)^2+(j-N/2)^2);
end
end
c = 1.05; %锐化参数,可调
Do = max(M,N);%本应该是方差,我简化了一下
H = (2.0-0.5)*(1 - exp(c*(-D/(Do^2)))) + 0.5;
F = F .* H;
哪位高手看懂的话给我讲讲啊,谢谢! |
|