拉普拉斯算子的应用.doc_第1页
拉普拉斯算子的应用.doc_第2页
拉普拉斯算子的应用.doc_第3页
拉普拉斯算子的应用.doc_第4页
免费预览已结束,剩余1页可下载查看

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

拉普拉斯算子的应用主要程序:1:function PQ=paddedsize(AB,CD,PARAM)if nargin=1 PQ=2*AB;elseif nargin=2&ischar(CD) PQ=AB+CD-1; PQ=2*ceil(PQ/2);elseif nargin=2 m=max(AB); P=2nextpow2(2*m); PQ=P,P;elseif nargin=3 m=max(AB CD); P=2nextpow2(2*m); PQ=p,p;else error(wrong number of inputs)end作用为提取填充参数。2:function g=dftfilt(f,H)F=fft2(f,size(H,1),size(H,2);g=real(ifft2(H.*F);g=g(1:size(f,1),1:size(f,2);此函数为滤波函数3:function u,v=dftuv(M,N)u=0:(M-1);v=0:(N-1);idx=find(uM/2);u(idx)=u(idx)-M;idy=find(vN/2);v(idy)=v(idy)-N;v,u=meshgrid(v,u);此函数的作用为计算任意点到指定点的距离。4:functionH,D=lpfilter(type,M,N,D0,n)u,v=dftuv(M,N)D=sqrt(u.2+v.2);switch type case ideal H=double(D=D0); casebtw if nargin=4 n=1; end H=1./(1+(D./D0).(2*n); casegaussian H=exp(-(D.2)./(2*(D02); otherwise error(unkown filter type.) end此函数为低通滤波器的函数5:function H=hpfilter(type,M,N,D0,n)if nargin=4 n=1;endHlp=lpfilter(type,M,N,D0,n);H=1-Hlp;此函数为高通滤波器的函数;主要过程:首先输入f=imread(D:照片收藏1.gif)将图像放进matlab环境中,然后使用p_down文件对图像进行低通处理。P_down文件:PQ=paddedsize(size(f);D0=0.05*PQ(1);H,D=lpfilter(ideal,PQ(1),PQ(2),D0);g=dftfilt(f,H);figure,imshow(g,)再使用p_high文件对图像进行高通处理P_high文件:PQ=paddedsize(size(f);D0=0.05*PQ(1);H=hpfilter(ideal,PQ(1),PQ(2),D0);g=dftfilt(f,H);figure,imshow(g,);处理的图像显示的结果为:第一个图为原图像。第二个

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论