数字图像处理实习报4-杨城-20131308077.docVIP

南京信息工程大学 实验（实习）报告 实验（实习）名称 图像锐化化算法实现 实验（实习）日期 2015.5.4 得分 指导教师 范春年 院系 *** 专业 **** 年级 **** 班次 * 姓名 ****** 学号 *********** 一、实验目的： 1．了解图像的目的和意义，巩固所学的图像的理论知识和相关算法； 2． 3．熟练掌握空域中常用的锐化滤波器；4．利用MATLAB程序进行图像Matlab自带的图像函数，但可以调用其它数学函数），观察图的效果。ameraman.tif）进行边缘检测，显示处理前图像和检测的边缘图 根据获得的梯度图，分别采用5种锐化输出处理方式，显示锐化输出前后的图像 三、实验要求： 用Matlab语言进行编程，实现上述功能，并尽量使得程序具有通用性，3种算子3个.m文件。 撰写实验报告并附上所用程序和结果。 四、实验代码及实验结果： Roberts算子： clear;clc; f=imread(cameraman.tif); [m,n]=size(f); f1=double(f); %Roberts算子 for i=1:m-1 for j=1:n-1 a=f1(i+1,j+1)-f1(i,j); b=f1(i+1,j)-f1(i,j+1); c(i,j)=abs(a)+abs(b); end end f2=uint8(c); %直接以梯度值输出 subplot(2,3,1),imshow(f),title(原图); subplot(2,3,2),imshow(f2),title(Roberts直接输出图); %辅以门槛判断 for i=1:m-1 for j=1:n-1 if (c(i,j)=125) f3(i,j)=c(i,j); else f3(i,j)=f1(i,j); end end end f3=uint8(f3); subplot(2,3,3),imshow(f3),title(Roberts门槛判断图); %给边缘规定一个特定灰度值 for i=1:m-1 for j=1:n-1 if(c(i,j)=125) f4(i,j)=255; else f4(i,j)=f1(i,j); end end end f4=uint8(f4); subplot(2,3,4),imshow(f4),title(Roberts边缘规定图); %给背景规定一个特定灰度级 for i=1:m-1 for j=1:n-1 if(c(i,j)=125) f5(i,j)=c(i,j); else f5(i,j)=0; end end end f5=uint8(f5); subplot(2,3,5),imshow(f5),title(Roberts背景规定图); %二值图像 for i=1:m-1 for j=1:n-1 if(c(i,j)=125) f6(i,j)=255; else f6(i,j)=0; end end end f6=uint8(f6); subplot(2,3,6),imshow(f6),title(Roberts二值图像); 效果图 Prewitt算子： clear; clc; f=imread(cameraman.tif); [m,n]=size(f); f1=double(f); for i=2:m-1 for j=2:n-1 a=f1(i+1,j-1)-f1(i-1,j-1)+f1(i+1,j)-f1(i-1,j)+f(i+1,j+1)-f1(i-1,j+1); b=f1(i-1,j+1)-f1(i-1,j-1)+f1(i,j+1)-f1(i,j-1)+f(i+1,j+1)-f1(i+1,j-1); c(i,j)=abs(a)+abs(b); end end f2=uint8(c); %直接以梯度值输出 subplot(2,3,1),imshow(f),title(原图); subplot(2,3,2),imshow(f2),title(Prewitt直接输出图);