matlab自编sobel算子,prewitt算子,log算子锐化图像.doc_第1页
matlab自编sobel算子,prewitt算子,log算子锐化图像.doc_第2页
matlab自编sobel算子,prewitt算子,log算子锐化图像.doc_第3页
matlab自编sobel算子,prewitt算子,log算子锐化图像.doc_第4页
matlab自编sobel算子,prewitt算子,log算子锐化图像.doc_第5页
免费预览已结束,剩余2页可下载查看

下载本文档

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

文档简介

课程名称 图像处理与模式识别 实验项目名称 图象的锐化 学生姓名 专业班级 学号 实验成绩 指导老师(签名 ) 日期 一. 实验目的和要求二. 实验内容、原理及实验结果与分析I=imread(coins.png);%读入图像subplot(2,2,1),imshow(I);title(原图像);%显示原图像H=fspecial(sobel); %应用sobel算子锐化图像I2=filter2(H,I); %sobel算子滤波锐化subplot(2,2,2);imshow(I2); %显示sobel算子锐化图像title(sobel算子锐化图像);H=fspecial(prewitt);%应用prewitt算子锐化图像I3=filter2(H,I);%prewitt算子滤波锐化subplot(2,2,3);imshow(I3); %显示prewitt算子锐化图像title(prewitt算子锐化图像);H=fspecial(log); %应用log算子锐化图像I4=filter2(H,I); %log算子滤波锐化subplot(2,2,4);imshow(I4);%显示log算子锐化图像title(log算子锐化图像);A=imread(baboon.bmp); A=rgb2gray(A); subplot(2,2,1),imshow(A),title(原始图象); image=double(A); m,n=size(A); u=zeros(1,9); k=zeros(1,9); for i=2:m-1, for j=2:n-1, u(1)=0*image(i,j);u(2)=2*image(i,j+1);u(3)=1*image(i-1,j+1); u(4)=0*image(i-1,j);u(5)=-1*image(i-1,j-1);u(6)=-2*image(i,j-1); u(7)=-1*image(i+1,j-1);u(8)=0*image(i+1,j);u(9)=1*image(i+1,j+1); rimage1(i,j)=abs(sum(u); k(1)=0*image(i,j);k(2)=0*image(i,j+1);k(3)=1*image(i-1,j+1); k(4)=2*image(i-1,j);k(5)=1*image(i-1,j-1);k(6)=0*image(i,j-1); k(7)=-1*image(i+1,j-1);k(8)=-2*image(i+1,j);k(9)=-1*image(i+1,j+1); rimage2(i,j)=abs(sum(k); xuhua(i,j)=rimage1(i,j)+rimage2(i,j); end, end, subplot(2,2,2), imshow(xuhua,map), title(Sobel算子锐化图象); for i=2:m-1, for j=2:n-1, u(1)=0*image(i,j);u(2)=1*image(i,j+1);u(3)=1*image(i-1,j+1); u(4)=0*image(i-1,j);u(5)=-1*image(i-1,j-1);u(6)=-1*image(i,j-1); u(7)=-1*image(i+1,j-1);u(8)=0*image(i+1,j);u(9)=1*image(i+1,j+1); rimage1(i,j)=abs(sum(u); k(1)=0*image(i,j);k(2)=0*image(i,j+1);k(3)=1*image(i-1,j+1); k(4)=1*image(i-1,j);k(5)=1*image(i-1,j-1);k(6)=0*image(i,j-1); k(7)=-1*image(i+1,j-1);k(8)=-1*image(i+1,j);k(9)=-1*image(i+1,j+1); rimage2(i,j)=abs(sum(k); xuhua(i,j)=rimage1(i,j)+rimage2(i,j); end, end, subplot(2,2,3), imshow(xuhua,map), title(Prewitt算子锐化图象 ); u=zeros(1,25); for i=3:m-2, for j=3:n-2, u(1)=24*image(i,j);u(2)=-2*image(i-2,j-2);u(3)=-4*image(i-2,j-1); u(4)=-4*image(i-2,j);u(5)=-4*image(i-2,j+1);u(6)=-2*image(i-2,j+2); u(7)=-4*image(i-1,j-2);u(8)=0*image(i-1,j-1);u(9)=8*image(i-1,j); u(10)=0*image(i-1,j+1);u(11)=-4*image(i-1,j+2);u(12)=-4*image(i,j-2); u(13)=8*image(i,j-1);u(14)=8*image(i,j+1);u(15)=-4*image(i,j+2); u(16)=-4*image(i+1,j-2);u(17)=0*image(i+1,j-1);u(18)=8*image(i+1,j); u(19)=0*image(i+1,j+1);u(20)=-4*image(i+1,j+2);u(21)=-2*image(i+2,j-2); u(22)=-4*image(i+2,j-1);u(23)=-4*image(i+2,j);u(24)=-4*image(i+2,j+1); u(25)=-2*image(i+2,j+2); rimage1(i,j)=abs(sum(u); xuhua(i,j)=rimage1(i,j); end, end, subplot(2,2,3), imshow(xuhua,map), title(log算子锐化图象)伪彩色增强A=imread(baboon.bmp);A1=rgb2gray(A);Im=double(A1); Im1=Im/256;H V=size(Im);temp=zeros(H,V,3);for m=1:H for n=1:V if Im1(m,n)=0.25 temp(m,n,3)=256*4*Im1(m,n); elseif Im1(m,n)=0.5 temp(m,n,3)=256*(2-4*Im1(m,n); temp(m,n,2)=256*4*(Im1(m,n)-0.25); elseif Im1(m,n)=0.75 temp(m,n,1)=256*4*(Im1(m,n)-0.5); temp(m,n,2)=256; elseif Im1(m,n)=1 temp(m,n

温馨提示

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

评论

0/150

提交评论