资源描述：
《激光光斑中心位置判定及其圆拟合》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第一章彩色图像的二值化1图像的二值化原理2 图像的二值化的程序实现3二值化前后效果对比第二章去除噪声1去除噪声的原理2去除噪声的程序实现3去除噪声前后的图像对比第三章圆拟合1圆拟合原理2圆拟合的程序实现3拟合效果第四章求圆心及半径第五章完整的Matlab程序第一章彩色图像的二值化第一节图像的二值化原理图像的二值化处理就是将图像上的点的灰度置为0或255，也就是讲整个图像呈现出明显的黑白效果。即将256个亮度等级的灰度图像通过适当的阀值选取而获得仍然可以反映图像整体和局部特征的二值化图像。在数字图像处理中，二值图像占有非常重要的地位，特别是在实用的图像处理中，以二值图像处理实现而构成的系统
2、是很多的，要进行二值图像的处理与分析，首先要把灰度图像二值化，得到二值化图像，这样子有利于再对图像做进一步处理时，图像的集合性质只与像素值为0或255的点的位置有关，不再涉及像素的多级值，使处理变得简单，而且数据的处理和压缩量小。所有灰度大于或等于阀值的像素被判定为属于特定物体，其灰度值为255表示，否则这些像素点被排除在物体区域以外，灰度值为0，表示背景或者例外的物体区域。第二节 图像的二值化的程序实现方法一：首先将图像转变为灰度图像，再利用max，min等函数求阈值分割点，最后转化为二值化图像相关程序：J=imread('1.jpg');%读图像figure;imshow(J);%显
3、示原始图像P=rgb2gray(J);%转换为灰度图像[m,n]=size(P);%获取图像的行数和列数ma=max(max(P));%求最大值mi=min(min(P));%求最小值limen=(ma+mi)/2;%求分割阈值I=(P>limen);%二值化figure;imshow(I);%显示二值化图像方法二：首先将图像转变为灰度图像，再利用graythresh等函数求阈值分割点，最后用函数im2bw进行二值化。J=imread('1.jpg');%读图像P=rgb2gray(J);%转换为灰度图像level=graythresh(P);%求分割阈值I=im2bw(P,level)
4、;%二值化imshow(I);%显示二值化图像第三节二值化前后效果对比第二章去除噪声第一节去除噪声的原理噪声在理论上可以定义为“不可预测，只能用概率统计方法来认识的随机误差”。实际获得的图像一般都因受到某种干扰而含有噪声。引起噪声的原因有敏感元器件的内部噪声、相片底片上感光材料的颗粒、传输通道的干扰及量化噪声等。噪声产生的原因决定了噪声的分布特性及它和图像信号的关系。第二节去除噪声的程序实现方法一：首先将图像变成灰度图像，然后进行中值滤波。I=imread('1.JPG');I=rgb2gray(I);subplot(2,3,1),imshow(I);title('原始图像');k1=m
5、edfilt2(I);%进行3*3滤波k2=medfilt2(I,[5,5]);%进行5*5滤波subplot(2,3,3),imshow(k1);title('3*3中值滤波后的图像');subplot(2,3,4),imshow(k2);title('5*5中值滤波后的图像');方法二：先利用bwlabel函数标注连通区域，在通过regionprops函数找出其中的噪声区域，并置零，以此达到去除噪声的目的。L=bwlabel(I);%功能：标注二进制图像中已连接的部分。stats=regionprops(L,{'Area','ConvexHull','MajorAxisLength'
6、,...'MinorAxisLength','Eccentricity','Centroid'});%用来度量图像区域属性的函数。测量标注矩阵L中每一个标注区域的一系列属性。L中不同的正整数元素对应不同的区域。%A=[];%定义数组A%fori=1:length(stats)A=[Astats(i).Area];%'Area'是标量，计算出在图像各个区域中像素总个数。%End[mA,ind]=max(A);I1=I;I1(find(L~=ind))=0;Figure;imshow(I1);第三节去除噪声前后的图像对比第三章圆拟合第一节圆拟合原理首先通过regionprops函数找出与所求
7、区域具有相同标准二阶中心矩的椭圆的长轴长度，短轴长度，离心率等，并找出包含所求区域的最小凸多边形，的顶点坐标。通过一定的方法，用所求椭圆的相关参数确定一个正多边形。使该正多边形具有一个相当大的边数，是该多边形接近一个圆。这个圆就是所要拟合的圆。在指令t=linspace(0,2*pi，N)中，N-1表示所求正多边形的边数，N越大，所拟合的圆越精确。第二节圆拟合的程序实现stats=regionprops(L,{'Area','Con
