瞳孔直径.ppt

《瞳孔直径.ppt》由会员分享，可在线阅读，更多相关《瞳孔直径.ppt（13页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、 1.图像的灰度变换图像的灰度变换 2.图像的灰度直方图图像的灰度直方图 3.图像的增强图像的增强 （一）图像的读入（一）图像的读入 在计算机中一幅二维数字图像可以用一个矩阵表示，该在计算机中一幅二维数字图像可以用一个矩阵表示，该矩阵中的元素是位于相应坐标位置的图像灰度值。对图像矩阵中的元素是位于相应坐标位置的图像灰度值。对图像处理，往往都是对像素灰度的操作，在图像处理中，直接处理，往往都是对像素灰度的操作，在图像处理中，直接对像素进行的操作称为对像素进行的操作称为空间域（简称空域）处理空间域（简称空域）处理。在。在matlab里读入一幅图片用函数里读入一幅图片用函数imread（）；例如读入

2、眼睛（）；例如读入眼睛图像语句为：图像语句为：im=imread（H：11.jpg）；）；figure，imshow 图片如下：图片如下：一、图像的灰度图像一、图像的灰度图像（二）图像的灰度变换二）图像的灰度变换 MATLAB中，用函数中，用函数imadjust可完成灰度图像的灰度变换。其基本语可完成灰度图像的灰度变换。其基本语法在格式为：法在格式为：gimadjust(p,low_in high_in,low_out high_out,gamma)其中，其中，p为读入到为读入到MATLAB中的待变换图像矩阵中的待变换图像矩阵;low_in high_in指定指定p图像中被执行变换操作的灰度范

3、围，图像中被执行变换操作的灰度范围，low_out high_out是是p图像中的像素变换后被映射到图像中的像素变换后被映射到low_outhigh_out的灰度级上。的灰度级上。默认状态下，默认状态下，low_out high_out及及low_in high_in为空，此时等价于为空，此时等价于0 1，表明输入与输出灰度级为，表明输入与输出灰度级为0255（uint8）或）或065535（uint16）。）。参数参数gamma是指变换映射的方式是指变换映射的方式，默认方式时，默认方式时gamma取作取作1，完，完成线性映射，这时变换前后的成线性映射，这时变换前后的灰度级没有被加权；灰度级没

4、有被加权；当当gamma小于小于1，则映射被加权，则映射被加权至更高的灰度级；至更高的灰度级；相反当相反当gamma大于大于1，输出则被，输出则被加权映射至较低的灰度级。加权映射至较低的灰度级。p=imread(football.jpg);g1=imadjust(p,0.1 0.6,);subplot(1,2,1);imshow(p);xlabel(a);subplot(1,2,2);imshow(g1);xlabel(b);例如例如：（三）彩色图像转换位灰度图像matlab中用函数rgb2gray（）来使彩色图像装化为灰度图像。语法为w=rgb2gray（q）;例如：代码为：A=imread

5、(H:11.jpg);subplot(1,2,1);imshow(A);B=rgb2gray(A);subplot(1,2,2),imshow(B)（四）对图像分为（四）对图像分为RGB三部分灰度函数为三部分灰度函数为 例如：A=imread(H:11.jpg);R=A(:,:,1);G=A(:,:,2);B=A(:,:,3);subplot(1,3,1);imshow(R);subplot(1,3,2),imshow(G);subplot(1,3,3),imshow(B)注：“：”表示方阵的大小，可省略。可用m n=size()二、图像的灰度直方图灰度直方图用imhist（）A=imread

6、(H:11.jpg);R=A(:,:,1);G=A(:,:,2);B=A(:,:,3);subplot(1,3,1);imhist(R);subplot(1,3,2),imhist(G);subplot(1,3,3),imhist(B)三、求眼睛直径代码：im=imread(1.jpg);im1=rgb2gray(im);im1=medfilt2(im1,10 10);medfilt2是一个二维中值滤波函数，其基本语法格式为：hp=medfilt2(p,m n,padopt)，m及n规定了邻域的大小，padopt指定了边界填充方式，默认方式为zerosBW=edge(im1,sobel);图像

7、边缘的检测 e,s=edge(p,method,para)method是边缘检测的类型，para为与edge对应的参数，e为与p同样大小的逻辑矩阵，在检测到边缘的位置时值为1，其它位置则为0，s为一可选参数。imx,imy=size(BW);msk=0 0 0 0 0;0 1 1 1 0;0 1 1 1 0;0 1 1 1 0;0 0 0 0 0;B=conv2(double(BW),double(msk);MATLAB图像处理工具箱中的函数bwlabel是把四连通或八连通的区域连接起来,调用语法为:L,num=bwlabel(BW,n)其中BW为输入图像;N可取值为4或8表示连接四连通或八连

8、通区域;NUM为找到的连通区域数目;L为输出图像矩阵,其元素值为整数,背景被标记为0,第一个连通区域被标记为1,第二个连通区域被标记为2,依此类推.l=bwlabel(B,8);mx=max(max(l);对灰度图片做圆形的边缘检测for i=1:mxr,c=find(l=i);rc=r c;sx sy=size(rc);n1=zeros(imx,imy);for j=1:sx x1=rc(j,1);y1=rc(j,2);n1(x1,y1)=1;end STATS=regionprops(n1,MinorAxisLength,MajorAxisLength,Area);ka=STATS.Min

9、orAxisLength/STATS.MajorAxisLength;if STATS.Area100&ka0.8 Centroid=regionprops(n1,Centroid);%中心 EquivDiameter=regionprops(n1,EquivDiameter);%直径figure,imshow(n1,);endend方法二 img=imread(E:m1.jpg);%读取图像读取图像1.jpgP=rgb2gray(img);%转换为灰度图像转换为灰度图像Pm,n=size(P);figure,imshow(P);%增加一个窗口，显示灰度图像增加一个窗口，显示灰度图像figure,imhist(P);%显示灰度图像直方图显示灰度图像直方图for ii=0:255 a(ii+1)=length(find(P=ii);end%读取灰度值所对应像素的个数读取灰度值所对应像素的个数s1=0;for ii=0:38 s1=s1+a(ii+1);end%计算灰度值为计算灰度值为0-38之间的像素的个数之间的像素的个数disp(s1=,num2str(s1);%输出输出s2=m*n;%计算图像总面积计算图像总面积pp=s1/s2;%0-38在总面积所占比例在总面积所占比例disp(pp=,num2str(pp);%输出输出