数字图像处理实验报告书.docVIP

数字图像处理实验报告 Page PAGE 3 of NUMPAGES 6 数字图像处理实验报告 一、实验内容及目的： 已知一幅图上有五个瓶子，装有液体。基于这幅图片，设计一个检测未完全装满的瓶子的解决方案。 二、实验思路分析： 形态学操作，使图像 形态学操作，使图像便于处理 二值化原图 找出不符合要求的标记块 显示原图 标记、计算 三、实验步骤 1、打开MATLAB软件，并将Current Folder改为所要处理的图片Fig1.jpg所在目录。我的图片所在目录为E:\12-13-2\数字图象处理\di： 2、读取和显示原始图像，显示原始图像的直方图： picOP = imread(Fig1.jpg); % 读入图像 figure,imshow(picOP),title(原始灰度图像); % 显示图像 figure, imhist(picOP), title(原始灰度图像直方图); % 显示原始图像的直方图 3、将图像进行二值化处理： 由于原始图像中白色比较明显，瓶身的颜色和背景比较相近，所以直接进行二值化处理，将未装满液体的信息提取出来。 picB = picOP; % 复制灰度图像到picB picBW = im2bw(picB,T/255); % 采用全局阈值进行灰度图像转变为二值图像 figure,imshow(picBW),title(全局阈值下二值图像); % 显示二值图像 4、进行形态学操作，转化为有利于处理的图像： picMORPHOLOGY = picBW ; % 复制准备形态学处理 se = strel(square,10); % 结构化元素 fo = imopen(picMORPHOLOGY,se); % 开操作 figure,imshow(fo),title(开运算之后图像); % 输出开运算之后图像 5、标记连通域，计算出各个白色连通区域（未装灌区域）面积和质心坐标： 图像处理后五个白色区域面积为： [2374; 2739; 8381; 2739; 1660] 五个质心的坐标数据： [20.79; 48.23; 138.8; 46.07; 256.7; 84.18; 376.8; 46.07; 484.7; 47.98] 6、根据面积和质心坐标计算结果判断那个装灌不合格： 通过没有装灌部分的面积大于2900或者质心大于50则判断为不合格，由步骤5的连通面积和质心坐标可得第三个瓶子不合格。 四、实验代码 clear all; close all; T = 200; % 全局阈值 200时效果要好一些，手动选出的值 %step1 读取和显示原始图像，显示原始图像的直方图 picOP = imread(Fig1.jpg); % 读入图像 figure,imshow(picOP),title(原始灰度图像); % 显示图像 figure, imhist(picOP), title(原始灰度图像直方图); % 显示原始图像的直方图 % step2 转化为二值图像 picB = picOP; % 复制灰度图像到picB picBW = im2bw(picB,T/255); % 采用全局阈值进行灰度图像转变为二值图像 figure,imshow(picBW),title(全局阈值下二值图像); % 显示二值图像 % step3 进行形态学操作，转化为有利于处理的图像 picMORPHOLOGY = picBW ; % 复制准备形态学处理 se = strel(square,10); % 结构化元素 fo = imopen(picMORPHOLOGY,se); % 开操作 figure,imshow(fo),title(开运算之后图像); % 输出开运算之后图像 % step4 计算出各个白色连通区域（未装灌区域）面积和质心坐标 [L,num] = bwlabel(fo,4); % 为白色连通区域做标签 stats = regionprops(L,basic); % 利用 regionprops（）函数计算面积、质心和边界等 s = [stats.Area]; % 矩阵S 中保存所有白色连通区域的面积 zx = [stats.Centroid]; % zx保存白色连通区域质心的坐标 % step5 根据面积和质心坐标计算结果判断那个装灌不合格 for indexnum = 1:num % 有实际几个连通区域就做几个循环,实际就是小于等于5个 if ( (s(indexnum) 2900) | (zx(2* indexnum) 50 ) ) % 如果面积大于2900或者质心大于50则为不合格 msgbox({不合格的瓶子