版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

sunal

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1547.2
帖子: 407
在线: 156小时
虫号: 1005201

[交流] 冈萨雷斯《数字图像处理》中的例子

第二版书中第四章图4.28运行结果和书中显示的不一样，请高手帮我看一下，程序是不是有问题？谢谢！
% 频域中拉普拉斯增强
clear all;
close all;
clc;
f=imread('Fig4.28(a).jpg');
subplot(221),imshow(f);
f1=im2double(f);
[M,N]=size(f1);
m=fix(M/2); n=fix(N/2);
for x=1:M
for y=1:N
      f2(x,y)=f1(x,y)*(-1)^(x+y);
end
end
f3=fft2(f2);
for u=1:M
for v=1:N
      h(u,v)=-((u-m/2)^2+(v-n/2)^2); % 拉普拉斯滤波器
      f4(u,v)=h(u,v)*f3(u,v);
end
end
f5=ifft2(f4);
r=real(f5);
for x=1:M
for y=1:N
      f6(x,y)=r(x,y)*(-1)^(x+y);
end
end
subplot(222),imshow(uint8(f6),[]);
% 对拉普拉斯变换后的图像重新标定
mi=min(min(f6));
ma=max(max(f6));
T=ma-mi;
for i=1:M
for j=1:N
      if f6(i,j)<0
         f7(i,j)=0;
      else
         f7(i,j)=(f6(i,j)-mi)/T*255;
      end
end
end
% 显示标定后的图像
subplot(223),imshow(uint8(f7),[]);
% 增强后的图像
f8=f1-f6;
subplot(224),imshow(f8,[]);
运行结果中的第三幅图不对，背景不应该是白色的，但不知错在哪了。

Fig4.28(a).jpg

运行结果

[ 来自科研家族控制&图像家族 ]

[ Last edited by sunal on 2012-5-7 at 21:28 ]

回复此楼

» 猜你喜欢

静等基金结果已经有16人回复
filecode 已经有6人回复
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人已经有11人回复
关于Filecode分析方法已经有3人回复
国自然结果已经有8人回复
FileCode能看出啥？已经有8人回复
应该是93bebmhtak前后十一个字符比较关键已经有22人回复
奇怪，两个人的filecode固定段从头到尾一模一样已经有4人回复
8月时间戳变的，举个手。玩一下，释放压力已经有13人回复
好奇怪的filecode 已经有5人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-05-07 20:56:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

关于英语语言要求，到底什么是 what which how ？有没有写好的例子给看下？已经有13人回复
求ANSYS用于焊接分析的生死单元的例子已经有4人回复
William T. Freeman的基于例子超分辨率matlab程序跑不起来啊！已经有6人回复
【分享】数字图像处理第二版、数字图像处理Matlab版（冈萨雷斯）【已搜索无重复】已经有198人回复
【求助】Orca 2.8的一个计算赝势的例子报错已经有12人回复
不予资助的委婉说法已经有13人回复
【求助】按照计算Dmol分子轨道的例子计算的结果不同？已经有16人回复
【求助】求GROMACS模拟表面活性剂在油水表面的行为的例子已经有7人回复
【求助】求gamess 输入使用赝势基组与分裂基组的例子已经有6人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 sunal 的主题更新