数字图像处理大作业.pdfVIP

数字图像处理大作业 [HW5] [编号 51] 学号： S 刘志辉 （1）程序代码 %灰度图像转换 RGB=imread(lzh.jpg); subplot(3,3,1);imshow(RGB);title( 原始图像 ); Y=imresize(RGB,[480 640]); subplot(3,3,2);imshow(Y);title( 缩小后的图像 ); I=rgb2gray(Y); %I1=double(I); subplot(3,3,3);imshow(I);title( 灰度图像 ); %均值滤波 K1=filter2(fspecial(average,5),I)/255; subplot(3,3,4);imshow(K1);title(5*5 均值滤波图像 ); %锐化 a=1; W1=[0 -a 0;-a 1+4*a -a;0 -a 0]; rh=imfilter(I,W1,symmetric,conv); subplot(3,3,5);imshow(rh);title( 拉普拉斯锐化图像 ); %直方图均衡化 J=histeq(I); subplot(3,3,6);imshow(I);title( 原始灰度图像 ); subplot(3,3,7);imhist(I);title( 原始图像直方图 );axis([0 255 0 7000]); subplot(3,3,8);imshow(J);title( 均衡化后的图像 ); subplot(3,3,9);imhist(J);title( 均衡化后的直方图 );axis([0 255 0 7000]); 处理结果图： （2）程序代码 %边缘检测 BW1=edge(I,sobel); BW2=edge(I,prewitt); BW3=edge(I,roberts); subplot(2,2,1);imshow(I);title( 原始灰度图像 ); subplot(2,2,2);imshow(BW1);title(sobel 边缘检测图像 ); subplot(2,2,3);imshow(BW2);title(prewitt 边缘检测图像 ); subplot(2,2,4);imshow(BW3);title(roberts 边缘检测图像 ); 运行结果： （3）程序代码 %数字分割 RGB=imread(lzh.jpg); Y=imresize(RGB,[480 640]); I=rgb2gray(Y); figure;imshow(I);title( 原始灰度图像 ); I1=im2bw(I,0.7); figure;imshow(I1);title( 二值化图像 ); I2=edge(I1,roberts,0.09,both); figure;imshow(I2);title(roberts 边缘检测图像 ); se=[1;1;1]; % 线型结构元素 I3=imerode(I2,se); %腐蚀图像 figure,imshow(I3),title( 腐蚀后边缘图像 ); se=strel(rectangle,[25,25]); % 矩形结构元素 I4=imclose(I3,se);% 图像聚类、填充图像 figure;imshow(I4),title( 填充后图像 ); I5=bwareaopen(I4,2500);% 去除聚团灰度值小于