数字图像处理——计数题目对一幅具有一些特征点的图像得到灰.DOC

数字图像处理——计数 题目 对一幅具有一些特征点的图像 得到灰度图、二值图，最后利用获得的二值图像对原图像中的特征点进行计数。 实验方法和关键步骤 在助教的指导下我们利用matlab的一些基本图像数据变换函数很快就能对原先的rgb3色图做灰度图和二值图，例如a = imread(star.jpg)为图片读取函数；b = rgb2gray(a)为三色图转灰度图函数；level = graythresh(b)为灰度图平均灰度值提取函数；bw = im2bw(b,level)灰度图转二值图函数；imshow(bw) 显示图像；imwrite(b,out.jpg) 保存图像。 然而对特征点的计数比较棘手，首先我们要定义怎样的点在二值图里才能算特征点，如何判定特征，怎样从许多具有特征的点中找出最具特征的点。其次，必须想出一种合理的算法能够使既不遗漏主要特征点还能够不重复特征点的计数。 具体实验过程与结果分析 1找出原图 特地挑选了比较具有特征的星空图作为分析对象,由于此图中主要能分辨出的有12颗星星，其亮度变化比较明显，而其他背景基本可视为黑色，所以我认为这幅图片可以作为理想的特征原图。 图1 星系原图 2做灰度图变换 经对原图的读取和灰度图转换函数的运用，我立刻得到了星系图的灰度图像如下图所示，从图中我们可以看到虽然从rgb变为灰度图，但亮度变化仍很丰富：图2星系灰度图 3做二值图变换 为了能够利于之后的数数分析，我特地对im2bw(b,level)中的level进行了手工调整，以致左上角的那颗星星不至于太暗。下图即为二值变换结果，从中我们仍能数出12颗星。 图3星系二值图 4数数算法 得到自己认为比较理想的分析对象后，经二值变换现在的星系图中只含有0，1两种数据，而且1表示两点，按图片内容提取数组先从上到下后从左往右遍历图像上每一个点。定义标志位flag1，当在某一列中有连续几个像素点的值都为1时，flag1置1，并记录连续出现1的个数，当连续数达到number（number认为是图中最小星星的对应连续点数目）时，在新数组data[]中相应位置记录；当每一列都遍历过之后，对daga数组进行，同样，当有连续number个元素都为1时，就可以认为特征对象已被找到，累加所有特征点的个数并最后输出结果。 5结果分析与小结 经运行程序可得12个特征点记录结果（如右图所示），与实际要求相吻合，证明了这个通过二值图分析计数的方法是可行的，但我认为本程序中还有些不足之处，就是如果纵向排列的点数比较多的话可能就会发生记录错误。但总体来说通过本次实验我还是很好的感受到了数字图像处理方法在实际应用中有着很大的价值，在此感谢老师和助教在课程中给予的帮助。 图4星系数数结果 源代码 counter = 0; node = 0; %计数 flag1 = 1; flag2 = 1; allinone = 0; %把连续的标记位归为一个点 temp = 0; number = 3; %经观察最小的星星连续点数为3 rgb = imread(star1.jpg); %读文件 imshow(rgb); gray = rgb2gray(rgb); %转灰度图 figure; imshow(gray); imwrite(gray,star2.jpg); %保存 bw=im2bw(gray,0.3); %转二值图 figure; imshow(bw); imwrite(bw,star3.jpg); [m,n]=size(bw); data(n+1)=0; %用数组记录下连续点的信息 for x= 1:n for y = 1:m value = bw(y,x); if(value == 1) if(flag2 == 1