| 查看: 1019 | 回复: 0 | |||
[交流]
冈萨雷斯《数字图像处理》中的例子
|
|
第二版书中第四章图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 ] |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
关于英语语言要求，到底什么是 what which how ？有没有写好的例子给看下？
已经有13人回复
- 求ANSYS用于焊接分析的生死单元的例子
已经有4人回复
- William T. Freeman的基于例子超分辨率matlab程序跑不起来啊！
已经有6人回复
- 【分享】数字图像处理第二版、数字图像处理Matlab版（冈萨雷斯）【已搜索无重复】
已经有198人回复
- 【求助】Orca 2.8的一个计算赝势的例子报错
已经有12人回复
- 不予资助的委婉说法
已经有13人回复
- 【求助】按照计算Dmol分子轨道的例子计算的结果不同？
已经有16人回复
- 【求助】求GROMACS模拟表面活性剂在油水表面的行为的例子
已经有7人回复
- 【求助】求gamess 输入使用赝势基组与分裂基组的例子
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
5