实验5图像锐化算法实现[精选].docVIP

南京信息工程大学 实验（实习）报告 实验（实习）名称 图像锐化算法实现 实验（实习）日期2013.12.2 得分 指导教师 范春年 专业软件工程年级2011 班次 1 姓名 陆超 学号 20112344015 实验目的： 1．了解图像的目的和意义，巩固所学的图像的理论知识和相关算法； 2． 3．熟练掌握空域中常用的锐化滤波器；4．利用MATLAB程序进行图像Matlab自带的图像函数，但可以调用其它数学函数），观察图的效果。close all; clear; clc; rice=imread(rice.png); imshow(rice); 分别利用Roberts、Prewitt 和Sobel 边缘检测算子，对一幅灰度数字图像（cameraman.tif）进行边缘检测，显示处理前图像和检测的边缘图 根据获得的梯度图，分别采用5种锐化输出处理方式，显示锐化输出前后的图像 close all; clear; clc; rice=imread(rice.png); imshow(rice); rice_noise=imnoise(rice,salt pepper,0.02); figure,subplot(231);imshow(rice_noise); rice_max=rice_noise; rice_min=rice_noise; rice_mid1=rice_noise;%中值 rice_mid2=rice_noise;%中点 rice_ave=rice_noise; [wide,high]=size(rice_noise); for i=1:wide; for j=1:high; temp=zeros(3,3); if(i==1j==1) temp(2:3,2:3)=rice_noise(i:i+1,j:j+1); continue; end if(i==1j==high) temp(2:3,1:2)=rice_noise(i:i+1,j-1:j); continue; end if(i==widej==1) temp(1:2,2:3)=rice_noise(i-1:i,j:j+1); continue; end if(i==widej==high) temp(1:2,1:2)=rice_noise(i-1:i,j-1:j); continue; end if(i==1) temp(2:3,1:3)=rice_noise(i:i+1,j-1:j+1); continue; end if(j==1) temp(1:3,2:3)=rice_noise(i-1:i+1,j:j+1); continue; end if(i==wide) temp(1:2,1:3)=rice_noise(i-1:i,j-1:j+1); continue; end if(j==high) temp(1:3,1:2)=rice_noise(i-1:i+1,j-1:j); continue; end if((i~=1)(j~=1)(i~=wide)(j~=high)) temp=rice_noise(i-1:i+1,j-1:j+1); continue; end rice_max(i,j)=max(max(temp)); rice_min(i,j)=min(min(temp)); rice_mid1(i,j)=round((rice_max(i,j)+rice_min(i,j))/2); temp=sort(temp); rice_mid2(i,j)=temp(2,2); rice_ave(i,j)=round(mean2(temp)); end end subplot(232);imshow(rice