声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1080|回复: 1

[编程技巧] 谁能帮我看看这段程序!

[复制链接]
发表于 2007-7-24 11:03 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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;

哪位高手看懂的话给我讲讲啊,谢谢!
回复
分享到:

使用道具 举报

发表于 2007-7-24 12:36 | 显示全部楼层
高斯锐化滤波
请参考http://forum.vibunion.com/forum/thread-48079-1-1.html
先弄懂原理再看程序
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 08:04 , Processed in 0.061121 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表