基于扫描电镜的颗粒物形态识别方法
【专利摘要】一种基于扫描电镜的颗粒物形态识别方法，适用于对空气中悬浮的颗粒物质进行识别。利用扫描电镜获取被测空气中放大5000倍的悬浮颗粒图像，并对颗粒图像二值化，针对二值化图像中存在粘连颗粒图案，在二值化图像中找出的楔形像素，并作为分离点的始端楔形像素，寻找始端楔形像素对应的终端楔形像素，然后采用Bresenham算法画分离线，完成粘连颗粒图像的分割，之后通过等效直径和形状因子方法对已经完成粘连颗粒分割后的各个颗粒图案的形态特征进行识别。其识别率高、可以得到图像中颗粒的具体形态，便于分类。
【专利说明】基于扫描电镜的颗粒物形态识别方法
【技术领域】
[0001]本发明涉及一种颗粒物形态识别方法，尤其适用于对空气中的颗粒物质进行识别时使用的基于扫描电镜的颗粒物形态识别方法。
【背景技术】
[0002]空气中的“总悬浮颗粒”，为直径小于100微米的颗粒物。其中10微米以上的，几乎都可被鼻腔和咽喉所捕集；10微米以下的可经过呼吸道沉积于肺泡；而2.5微米以下能穿过人体内的末端细支气管到达气体交换区域。不同的环境中的颗粒物的具体形态特征不同，因而人们在不同环境中所需的防护措施也并不相同。
[0003]目前，空气污染中悬浮颗粒物的测量方法一般是通过颗粒采集器过滤膜的增重法来判断，通过将过滤膜上的颗粒收集并称重从而得到被测空气中存在的颗粒总重量，并通过颗粒总重量和被测量空气的体积从而判断出被测空气中的颗粒污染浓度，但是此测量方法不能获得空气中污染颗粒的具体形态，从而对空气污染的防护和治理产生影响。
【发明内容】
[0004]技术问题:本发明的目的是为克服现有技术的不足，提供了一种识别率高、可以得到图像中颗粒的具体形态的基于扫描电镜的颗粒物形态识别方法。
[0005]技术方案:本发明的基于扫描电镜的颗粒物的图像识别方法，包括以下步骤:
[0006]a.利用扫描电镜的低真空模式获取被测空气中放大的悬浮颗粒图像；
[0007]b.用矩阵实验室计算机对放大后的颗粒图像采用Otsu算法获取二值化后的颗粒图像；
[0008]c.对二值化后颗粒图像中存在相互粘连的颗粒图案进行分割，以提高颗粒图像中对颗粒的准确识别，颗粒图案分割过程为:
[0009](i)在二值化图像中，背景像素的灰度值为0，颗粒像素的灰度值为255，根据边界像素是由颗粒向背景转变的中间过程特性，确定二值化颗粒图像中的所有边界像素，对找到的每个边界像素，依次以边界像素为中心边界像素，取半径R为5个像素的圆形区域进行循环搜索，得到循环区域内的边界像素总数N，由公式:
Iv.1 V
[0010]Igmv Xi?，Jgrav
/V η^γIy
[0011]得到循环搜索区域中除中心边界像素外所有边界像素的重心agMV，jgMV);式中(i?, j?)为边界像素坐标值；
[0012]通过求无量纲偏心率E依次对每个中心边界像素进行判断是否为楔形像素，中心边界像素距离边界像素重心(igMV，Jgrav)的距离r为:
_3] r = ^i2grar+J2grar
[0014]无量纲偏心率￡ = i，式中半径R为5，将预设阈值CE与无量纲偏心率E进行比较，
如果E > CE，则判断此中心边界像素为待匹配的楔形像素，并将计算获得的所有待匹配楔形像素坐标按行依次写入数组ijwedge中，数组ijwedge中每一行代表一个楔形像素的坐标。
[0015](ii)从数组ijwedge中获取第一个的待匹配楔形点坐标值，将待匹配楔形像素作为线性分割的始端楔形像素，定义一个搜索区域，所述搜索区域以始端楔形像素为起始点，方向为与中心边界像素到重心的矢量方向相反，搜索半径为L，两边拓展的角度不超过η/3的扇形区域，其中搜索半径L为预设值，在搜索区域内通过双层循环逐行找寻楔形像素，如果找到多个楔形像素，按照与始端楔形像素欧氏距离最近的原则，从找到的楔形像素中确定终端楔形像素，将始端楔形像素与终端楔形像素视为分离点，采用Bresenham算法画分离线，完成粘连颗粒图像的分割，并将已经匹配的两个楔形像素在数组ijwedge中的坐标值标记；
[0016](iii)从数组ijwedge中找寻下一个未标记的楔形像素点坐标值，重复步骤ii中的操作，直到数组ijwedge中的楔形像素全部操作结束；
[0017]d.利用计算机通过等效直径和形状因子方法对已经完成粘连颗粒分割后的各个颗粒图案的形态特征进行识别，并将识别结果通过显示器显示出来。
[0018]所述对已经完成颗粒分割后的各个颗粒图案的形态特征的识别，包括:识别每个颗粒图案的等效直径de。和形状因子S ;先对得到颗粒分割后的二值化图像进行识别，找出像素灰度值为255的目标区域点的个数，结合单个像素点代表的实际长度求出颗粒的面积Area ;根据二值化图像使用Canny边缘检测算子得出其边缘图像，找出边缘图像中像素灰度值为255的边缘像素点的个数，继而求出颗粒的周长Perimeter。之后求得与颗粒面积Area相同的圆的直径，即为等效直径dec，然后根据自定义的比率公式S = (Area/ π )1/2/(Perimeter/2 η )求得形状因子S ;所述放大倍数为5000倍；所述数组ijwedge为一个10000行2列的数组，数组ijwedge中每一行代表一个楔形像素的坐标；所述标记在数组ijwedge中的坐标值的方法为，将已经匹配的楔形像素的坐标值置O。
[0019]有益效果:本发明提供了一种识别率高的基于扫描电镜的细颗粒物的形态识别方法，通过将拍摄后放大的空气图片进行二值化使图片中的颗粒图像凸显出来，并对二值化后颗粒图像中存在相互粘连的颗粒图案进行分割，防止了对颗粒图形识别时将相互粘连的数个颗粒识别成一个单独的颗粒，提高颗粒图像中对颗粒的准确，并通过求颗粒的等效直径和形状因子来获得单个颗粒的具体形态特征，这一系列操作为人们对颗粒物的防护起到有效的指导作用。
【专利附图】
【附图说明】
[0020]图1是本发明利用扫描电镜获取被测空气中放大的悬浮颗粒图像粘连示意图；
[0021]图2是本发明寻找待匹配楔形点的循环搜索区域示意图；
[0022]图3是本发明寻找待匹配楔形点的搜索区域示意图；
【具体实施方式】
[0023]下面结合附图对本发明的实施作进一步的描述，
[0024]本发明的基于扫描电镜的细颗粒物图像识别方法，其步骤如下:
[0025]a.如图1所示，利用扫描电镜的低真空模式获取被测空气中放大5000倍的悬浮颗粒图像，所述低真空模式为扫描电镜在10.6mm的工作距离和15kv的电压下获取颗粒图片，;
[0026]b.用矩阵实验室计算机对放大5000倍颗粒图像采用Otsu算法获取二值化后的颗粒图像；
[0027]c.对二值化后颗粒图像中存在相互粘连的颗粒图案进行分割，图1为粘连颗粒示意图，为了提高颗粒图像中对颗粒的准确识别，颗粒图案分割过程为:
[0028](i)在二值化图像中，背景像素的灰度值为0，颗粒像素的灰度值为255，根据边界像素是由颗粒向背景转变的中间过程特性，确定二值化后颗粒图像中的边界像素。对二值化后颗粒图像中的所有边界像素，依次以每个边界像素为中心边界像素，取半径R为5个像素的圆形区域进行循环搜索，得到循环区域内的边界像素总数N。如图2所示，标有数字2的为边界像素，标有数字4的像素为中心边界像素,标有数字3的为颗粒像素,标有数字I的为背景像素，计算循环搜索区域中所有标有数字2的边界像素的重心g。由公式:
【权利要求】
1.一种基于扫描电镜颗粒物形态识别方法，其特征在于包括以下步骤: a.利用扫描电镜的低真空模式获取被测空气中放大的悬浮颗粒图像； b.用矩阵实验室计算机对放大后的颗粒图像采用Otsu算法获取二值化后的颗粒图像； c.对二值化后颗粒图像中存在相互粘连的颗粒图案进行分割，以提高颗粒图像中对颗粒的准确识别，颗粒图案分割过程为: (i)在二值化图像中，背景像素的灰度值为0，颗粒像素的灰度值为255，根据边界像素是由颗粒向背景转变的中间过程特性，确定二值化颗粒图像中的所有边界像素，对找到的每个边界像素，依次以边界像素为中心边界像素，取半径R为5个像素的圆形区域进行循环搜索，得到循环区域内的边界像素总数N，由公式:
得到循环搜索区域中除中心边界像素外所有边界像素的重心Qgrav，jgrav);式中(i?, j?)为边界像素坐标值； 通过求无量纲偏心率E依次对每个中心边界像素进行判断是否为楔形像素，中心边界像素距离边界像素重心(igMV，Jgrav)的距离r为:
无量纲偏心率E = |，式中半径R为5，将预设阈值CE与无量纲偏心率E进行比较，如果E > CE，则判断此中心边界像素为待匹配的楔形像素，并将计算获得的所有待匹配楔形像素坐标按行依次写入数组ijwedge中，数组ijwedge中每一行代表一个楔形像素的坐标。 (ii)从数组ijwedge中获取第一个的待匹配楔形点坐标值，将待匹配楔形像素作为线性分割的始端楔形像素，定义一个搜索区域，所述搜索区域以始端楔形像素为起始点，方向为与中心边界像素到重心的矢量方向相反，搜索半径为L，两边拓展的角度不超过π/3的扇形区域，其中搜索半径L为预设值，在搜索区域内通过双层循环逐行找寻楔形像素，如果找到多个楔形像素，按照与始端楔形像素欧氏距离最近的原则，从找到的楔形像素中确定终端楔形像素，将始端楔形像素与终端楔形像素视为分离点，采用Bresenham算法画分离线，完成粘连颗粒图像的分割，并将已经匹配的两个楔形像素在数组ijwedge中的坐标值标记； (iii)从数组ijwedge中找寻下一个未标记的楔形像素点坐标值,重复步骤ii中的操作，直到数组ijwedge中的楔形像素全部操作结束； d.利用计算机通过等效直径和形状因子方法对已经完成粘连颗粒分割后的各个颗粒图案的形态特征进行识别，并将识别结果通过显示器显示出来。
2.根据权利要求1所述的基于扫描电镜的细颗粒物的形态识别方法，其特征在于:所述对已经完成颗粒分割后的各个颗粒图案的形态特征的识别，包括:识别每个颗粒图案的等效直径和形状因子S ;先对得到颗粒分割后的二值化图像进行识别，找出像素灰度值为255的目标区域点的个数，结合单个像素点代表的实际长度求出颗粒的面积Area ；根据二值化图像使用Canny边缘检测算子得出其边缘图像，找出边缘图像中像素灰度值为255的边缘像素点的个数，继而求出颗粒的周长Perimeter。之后求得与颗粒面积Area相同的圆的直径，即为等效直径dec，然后根据自定义的比率公式S = (Area/ π )1/2/(Perimeter/2 η )求得形状因子S。
3.根据权利要求1所述的基于扫描电镜的细颗粒物的形态识别方法，其特征在于:所述放大倍数为5000倍。
4.根据权利要求1所述的基于扫描电镜的细颗粒物的形态识别方法，其特征在于:所述数组ijwedge为一个10000行2列的数组，数组ijwedge中每一行代表一个楔形像素的坐标。
5.根据权利要求1所述的基于扫描电镜的细颗粒物的形态识别方法，其特征在于:所述标记在数组ijwedge中的坐标值的方法为，将已经匹配的楔形像素的坐标值置O。
【文档编号】G06K9/60GK104182759SQ201410413256
【公开日】2014年12月3日 申请日期:2014年8月20日 优先权日:2014年8月20日
【发明者】卢兆林, 冯瑾, 李世银, 陈燕 申请人:徐州坤泰电子科技有限公司