| 查看: 4929 | 回复: 43 | ||||||
| 【奖励】 本帖被评价39次，作者zhaoyuan0426增加金币 30.2 个 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
[资源]
利用Matlab模拟空间滤波实验
|
||||||
|
利用Matlab模拟空间滤波实验 阿贝成像原理是在透镜后焦面上得到光场空间频率分布的傅立叶变换,成像又是一次逆变换的过程,这种变换可由快速傅立叶变换（FFT）轻松实现。利用阿贝—波特实验装置和空间滤波系统,从改变频谱入手改造一幅光学图像,可以进行光学信息处理。本文在此基础上,用Matlab强大的计算及图像可视化功能完成阿贝—波特实验的物理模型的构建并进行计算机模拟,从而实现数字图像的处理。   二维光栅的频谱 x=ones(100,100) %创建矩阵 x(1:9:100,：)=0 %得到1 维光栅 y=x.*(x’) %得到2 维光栅 C=fft2(y,200,200) %傅立叶变换 p=abs(fftshift(z)) %变换像限并取模 imshow(0.01*p) %以一定比例显示图像 图11 即是二维光栅y,图12即是其频谱。可以看出,模拟的结果与实验中所观测到的图像是一致的。而且模拟图像更加清晰结果更加准确。  网格“光”字的滤波 滤波程序如下： clear f=imread('w1.bmp'); figure,imshow(f);title('原始图像'); PQ=paddedsize(size(f)); [U,V]=dftuv(PQ(1),PQ(2)); D0=0.03*PQ(2); F=fft2(f,PQ(1),PQ(2));%傅立叶变换 g=abs(fftshift(F)); figure,imshow(0.00001*g);title('点阵图'); x=1:15:1000;y=1:15:1000; figure,mesh(x,y,g(x,y));colormap(jet); H=lpfilter('gaussian',PQ(1),PQ(2),D0); g=dftfilt(f,H);H=fftshift(H); figure,mesh(H(1:10:1000,1:10:1000)); axis([0 100 0 100 0 1]) figure, imshow(H,[]) figure, imshow(g,[])   光栅为一周期性函数,其频谱是有规律排列的离散点阵,而字迹不是周期性函数,它的频谱是连续的,一般不容易看清楚,由于光字笔画较粗,其空间低频成分较多,因此频谱面的光轴附近只有光字信息而没有网格信息。经过低通滤波,滤掉了物的高频信息,处理图像显示了物的低频信息即光字。由于经过滤波后的图像能量有所损失导致输出图像比原图像模糊。 零频滤波 主程序如下: clear I1=imread('w3.bmp'); imshow(I1);title('原始图像'); figure f=double(I1); g=fft2(f); %傅立叶变换 g=fftshift(g); [result1,result2]=filt(g); result1=ifftshift(result1); result2=ifftshift(result2); X2=ifft2(result1); X3=uint8(real(X2)); imshow(X3); title('Butterworth滤波图像'); 显示出设计的一个高通滤波器程序如下： H=fftshift(hpfilter('gaussian',500,500,20)); mesh(H(1:10:500,1:10:500));colormap(jet); axis([0 50 0 50 0 1]) 其平面显示如图22所示,立体图如图23所示。   [ Last edited by zhaoyuan0426 on 2012-4-15 at 17:08 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 光学信息处理实验.doc
- 附件 2 : 空间滤波的计算机模拟.doc
2012-04-15 17:01:52, 493 K
2012-04-15 17:01:53, 1.62 M
» 收录本帖的淘帖专辑推荐
 matlab学习专辑 | 小木虫matlab全部资源 | 课题相关 | 好好学习天天向上 |
» 猜你喜欢
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有235人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助啊......matlab
已经有10人回复
- MATLAB用于三元回归参数
已经有3人回复
- matlab在概率统计中的计算
已经有76人回复
- Matlab图像处理问题
已经有10人回复
- matlab处于busy
已经有7人回复
- matlab反求未知输入参数
已经有4人回复
- 如何用matlab画出两组实验测定数据概率密度曲线，求两条曲线重叠部分面积
已经有8人回复
- 如何通过matlab把文件夹中的部分文件拷贝到另一个文件夹中?
已经有4人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- MATLAB 编程求解普通实验数据时间序列的最大李雅普诺夫指数！！！
已经有12人回复
- 利用MATLAB 分析图片颜色之间的差别
已经有4人回复
- 【求助】利用MATLAB求解的外围轮廓包络线为什么不与图形外围轮廓重合？怎样才可以精确
已经有8人回复
- 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码
已经有8人回复
9楼2012-05-25 21:38:27
2楼2012-04-15 17:12:14
8楼2012-05-25 21:38:04
12楼2012-09-22 08:15:43
简单回复
2012-04-15 21:33
回复
五星好评 顶一下，感谢分享！
2012-04-21 13:29
回复
五星好评 顶一下，感谢分享！
2012-05-13 22:49
回复
三星好评 顶一下，感谢分享！
diligent627811楼
2012-09-20 21:00
回复
五星好评 顶一下，感谢分享！
rabbitting13楼
2012-10-23 11:08
回复
五星好评 顶一下，感谢分享！