数字图像处理实验报告.doc

数字图像处理实验报告数字图像处理实验报告一、实验内容及目的：已知一幅图上有五个瓶子，装有液体。基于这幅图片，设计一个检测未完全装满的瓶子的解决方案。二、实验思路分析：形态学操作，使图像便于处理二值化原图找出不符合要求的标记块显示原图标记、计算三、实验步骤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; % 复制灰度图像到picBpicBW = 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.6251053074979;48.0433866891323;138.604600219058;46.6575392479007;256.266674621167;84.9748240066818;376.604600219058;46.6575392479007;484.618072289157;47.90843373493986、根据面积和质心坐标计算结果判断那个装灌不合格：通过没有装灌部分的面积大于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; % 复制灰度图像到picBpicBW = 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* indexnu