用于粒子检测仪中散点图的粒子分类方法
【专利摘要】本发明公开一种用于粒子检测仪中散点图的粒子分类方法，其包括步骤：A、将散点图转换成地形图，输出等高线图；B、按第一层间步长从灰度值最大的像素点所在的最高层开始查找每一层等高线图；C、判断第i层等高线图中的一个闭合等高线是否可以通过预设归类规则与第i+1层等高线图中类别P的像素点之间产生联系，如是，则把第i层等高线图中的该闭合等高线投影到散点图上所包围的像素点归入类别P，否则将不能归属于任何已知类别的新顶点标记为一个新类别，其中i为自然数；D、逐层归类直到达到第1层等高线图时判断是否满足终止条件。本发明可以提高散点图中粒子分类精度。
【专利说明】用于粒子检测仪中散点图的粒子分类方法
【技术领域】
[0001]本发明涉及一种医学图像处理技术，尤其是涉及一种用于粒子检测仪中对散点图通过自适应查找各个类别进行粒子分类的方法。
【背景技术】
[0002]粒子检测仪用于对各种粒子进行分类统计，例如血液细胞分析仪，其主要功能就是提供人体中各类细胞的数目，其中对各类白细胞的计数是基本的功能。使经过试剂处理后的、被液体包被的白细胞逐个通过检测区域(流动室)，将激光照射到单个的细胞上，通过收集不同角度的折射或散射光线(通常为前向散射光FSCdIU向散射光SSC)，再通过光电转换和AD转换，得到对应细胞的一组二维数据，将此二维数据映射到坐标系中，可以得到该细胞在二维坐标系中的位置，将样本中的所有白细胞都映射到坐标上可以得到二维的白细胞分布图(简称散点图)。在散点图中，同一类的白细胞聚集在一起，不同类的白细胞互相分离，如图1所示。
[0003]在散点图上划分多个区域，把落在统一区域的白细胞归为同一类，并统计落在这些类别内的粒子数据和百分比，用于分析被测样本的成分组成。
[0004]实际上，由于样本多种多样，散点图的形态也是多种多样，见错误！未找到引用源。A、图2B和图2C所示。各类粒子群落之间的分界线可能不明显；群落的形态可能不相同；某些群落粒子数有时较高而有时较低，甚至不存在；群落在散点图所处的位置会发生改变。以上这些因素都会影响到粒子分群，因此，简单地划分区域无法准确地提取出这些分
群信息。
【发明内容】
[0005]本发明提出一种用于粒子检测仪中散点图的粒子分类方法，以提高散点图中粒子分类精度。
[0006]本发明采用如下技术方案实现:一种用于粒子检测仪中散点图的粒子分类方法，其包括步骤:
[0007]A、将散点图转换成地形图，并将地形图以等高线图的形式输出；
[0008]B、按预设的第一层间步长，从灰度值最大的像素点所在的最高层开始查找每一层等闻线图；
[0009]C、将第i+Ι层等高线图中灰度值最大的像素点归入类别P，判断第i层等高线图中的一个闭合等高线是否可以通过预设归类规则与第i+Ι层等高线图中类别P的像素点之间产生联系，如是，则把第i层等高线图中的该闭合等高线投影到散点图上所包围的像素点归入类别P，否则将不能归属于任何已知类别的新顶点标记为一个新类别，其中i为自然数；
[0010]D、逐层归类直到达到第I层等高线图时判断是否满足终止条件，若否，则在第I层等高线图与像素灰度值为O的第O层等高线图之间采用预设的第二层间步长进行再分层，然后转入步骤C，否则结束；
[0011]其中，第二层间步长<第一层间步长/2。
[0012]其中，所述预设归类规则包括:类别P的像素点中灰度值最大的像素点或类别P的像素点在第i+ι层等高线图中形成的闭合等高线在第i层等高线图中的投影被包含在第i层等高线图中的一个闭合等高线内。
[0013]其中，所述预设归类规则包括:在第i层等高线图中的一个闭合等高线包围的像素点到类别P中灰度值最大的像素点在第i层等高线图的投影距离处于一个预设范围值之内。
[0014]其中，所述终止条件包括:所发现的类别数量超过设定阈值。
[0015]其中，假设散点图中所包含的粒子类别数量N，其中有η个类别的粒子无法识别，则设定阈值取Ν-η，其中N为整数，η小于N的整数。
[0016]与现有技术相比，本发明具有如下有益效果:
[0017]本发明是一个精度可变、可以适应不同用途的精度要求的散点图粒子分类方法，该方法较少地依赖于散点图中分类的形状、位置信息，即只要知道粒子的类别数量，即使粒子群出现的位置和形状发生改变，该方法也能通过查找新增类的方式，一个一个把粒子群找出来。另外，本发明通过查找闭合等高线的方式，能够智能地识别散点图中的各类边界，包括不规则的边界，而不是只能处理圆形、方形或椭圆形等规则边界，从而可以提高散点图中粒子分类精度。
【专利附图】
【附图说明】
[0018]图1是现有散点图的示意图
[0019]图2Α、图2Β和图2C分别是现有散点图的几种粒子群落形态分布不意图；
[0020]图3是本发明个优选实施例的流程示意图；
[0021]图4是将散点图转换生成地形图的示意图；
[0022]图5是将地形图转换生成等高线图的过程示意图；
[0023]图6是等高线图的示意图。
【具体实施方式】
[0024]为了弥补现有对散点图划分区域进行粒子分类方法存在的不足，本发明提出了一种提高散点图中粒子分类精度的方法，能够自适应查找各个类别，提供稳定、准确的分类结果，该方法尤其适用于粒子分析仪用，也可适用于所有存在粒子聚类的散点图的分类识别。
[0025]如图3所示，本发明的一个优先实施包括如下步骤:
[0026]步骤S1、以散点图的长和宽分别为X轴与Y轴、以散点图中像素点的灰度值为Z轴建立三维坐标，将散点图中每一个粒子在散点图中位置所在像素点的灰度值表示为该像素点的海拔高度，形成一张地形图，如图4所示。
[0027]步骤S2、将地形图以等高线图的形式输出。
[0028]在图4中，等高线是Z轴上高程相等(即灰度值相同)的各相邻点所连成的闭合曲线。以等高线在X-Y平面上投影后形成等高线图，如图5和图6所示。
[0029]步骤S3、按预设层间步长查找每一层等高线图。[0030]地形图上灰度值最大的像素点对应记为顶点，顶点为等高线的最高层。从最高层(顶点)开始，按预设层间步长在等高线图上逐层向下查找每一层等高线图。
[0031]层间步长也就是每一层等高线图包含的灰度等级，可以根据具体情况决定。需要精度越高时，层间步长越短，分的层数越多。层间步长与总灰度等级满足以下关系:分层数=总灰度等级/层间步长。
[0032]原则上可以先取一个较长的层间步长，如把整个地形图分成3-6层，可以有效提
高计算效率。
[0033]步骤S4、将第i+Ι层等高线图中灰度值最大的像素点归入类别P，判断下一层(第i层，i为自然数)等高线图中的一个闭合等高线是否可以通过预设归类规则与上一层(第i+Ι层)等高线图中某一类别p(p无实际含义，只是为了描述时区分的方便)的像素点之间产生联系，如是，则转入步骤S5，否则转入步骤S6。
[0034]把地形图上灰度值最大的像素点(顶点)记为第一类别的像素点。
[0035]如果第i层等高线图的一个闭合等高线通过预设归类规则能够与第i + Ι层等高线图中某个类别P的像素点之间产生联系:如果该类别P的像素点中顶点(即该类别的像素点中灰度值最大的像素点)或该类别P的像素点在第i+ι层等高线图中形成的闭合等高线在第i层等高线图中的投影被包含在第i层等高线图中的一个闭合等高线内，或者在第i层等高线图中的一个闭合等高线包围的像素点到该类别P中灰度值最大的像素点在第i层等高线图的投影距离处于一个预设范围值之内，则通过步骤S5把第i层等高线图中的该闭合等高线投影到散点图上所包围的像素点归入该类别P，而通过步骤S6将不能归属于任何已知类别的新顶点用来标记一个新类别。
[0036]步骤S7、逐层归类直到达到最下层(第I层)时判断是否满足终止条件(如所发现的类别数量超过设定阈值，比如该散点图中所包含的粒子类别数量N(N为整数)，或者如果有η个(η小于N的整数)类别不便于用这种方法识别却能通过其他方法识别，设定阈值也可以取Ν-η，这些可以根据需要分析的图像样本的具体情况而定)，若否，则转入步骤S8，否则结束
[0037]步骤S8、如果查找到最下层(由于存在预设层间步长，故最下层是灰度值不为O的底层)时未达到终止条件，在非O的最下层(第I层)与像素灰度值为O的第O层之间采用更较小的层间步长进行再分层，这时候分层数与这次所取的层间步长的关系是:
[0038]分层数=步骤S3的层间步长(或第一层间步长)/步骤S8的层间步长(或第二层间步长)。一般而言，第二层间步长<第一层间步长/2。
[0039]重新分层后，在第I层与第O层之间，按照新的分层重复步骤S4-步骤S7直到达到终止条件。
[0040]该方法介绍了一些在分层后归类的原则，但不限于这些方法。通过改变归类规则和终止条件可以调节分类识别的精度。
[0041]找到散点图上各种粒子的类别以后，可以通过所分析图像的特点，将找到的每一类别与实际需要解析的类别一一对应。对应规则需要通过所分析的具体样本确定，比如各类顶点的相对位置，或者各类粒子数。至此，通过这个方法就可以自动地从一个粒子散点图中把所需要查找的各类粒子提取出来并进行分析。以血液分析仪为例，这些粒子信息可以用作病人的疾病诊断。[0042]以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种用于粒子检测仪中散点图的粒子分类方法，其特征在于，所述方法包括步骤: A、将散点图转换成地形图，并将地形图以等高线图的形式输出； B、按预设的第一层间步长，从灰度值最大的像素点所在的最高层开始查找每一层等高线图； C、将第i+Ι层等高线图中灰度值最大的像素点归入类别P，判断第i层等高线图中的一个闭合等高线是否可以通过预设归类规则与第i+Ι层等高线图中类别P的像素点之间产生联系，如是，则把第i层等高线图中的该闭合等高线投影到散点图上所包围的像素点归入类别P，否则将不能归属于任何已知类别的新顶点标记为一个新类别，其中i为自然数； D、逐层归类直到达到第I层等高线图时判断是否满足终止条件，若否，则在第I层等高线图与像素灰度值为O的第O层等高线图之间采用预设的第二层间步长进行再分层，然后转入步骤C，否则结束； 其中，第二层间步长<第一层间步长/2。
2.根据权利要求1所述用于粒子检测仪中散点图的粒子分类方法，其特征在于，所述预设归类规则包括:类别P的像素点中灰度值最大的像素点或类别P的像素点在第i+Ι层等高线图中形成的闭合等高线在第i层等高线图中的投影被包含在第i层等高线图中的一个闭合等闻线内。
3.根据权利要求1所述用于粒子检测仪中散点图的粒子分类方法，其特征在于，所述预设归类规则包括:在第i层等高线图中的一个闭合等高线包围的像素点到类别P中灰度值最大的像素点在第i层等高线图的投影距离处于一个预设范围值之内。
4.根据权利要求1所述用于粒子检测仪中散点图的粒子分类方法，其特征在于，所述终止条件包括:所发现的类别数量超过设定阈值。
5.根据权利要求4所述用于粒子检测仪中散点图的粒子分类方法，其特征在于，假设散点图中所包含的粒子类别数量N，其中有η个类别的粒子无法识别，则设定阈值取Ν-η，其中N为整数，η小于N的整数。
【文档编号】G01N15/10GK103776751SQ201210407597
【公开日】2014年5月7日 申请日期:2012年10月23日 优先权日:2012年10月23日
【发明者】卓远, 易晗平 申请人:深圳市蓝韵实业有限公司