一种骨髓细胞形态学特征的提取方法与系统

1.本发明涉及医学图像特征提取的技术领域，尤其是指一种骨髓细胞形态学特征的提取方法与系统。


背景技术：

2.骨髓细胞形态学检查是诊断血液系统疾病、判断疗效及检测病情变化的最简单也是最重要的手段之一，在临床检验工作中通过显微镜观察染色后的骨髓细胞，以直接了解骨髓中各种血细胞数量、比例、形态以及有无异常细胞等。迄今，人工镜检仍是骨髓形态学检查的主要方法，但其检验过程耗时耗力，效率低下。近年来人工智能技术在医学图像领域发展迅速，在骨髓细胞分割及识别等方面提供了许多新的思路。
3.骨髓细胞形态学提取一般包括如下步骤：1)对骨髓细胞图像进行图像预处理；2)对预处理后的骨髓细胞图像进行细胞分割，得到单一的骨髓细胞图像；3)对单一骨髓细胞图像进行细胞形态学特征提取。骨髓细胞的细胞区域和细胞核区域是计算骨髓细胞形态学特征的重要依据，如何得到精准的细胞区域和细胞核区域至关重要。人工智能技术背景下应运而生的图像分割算法如unet模型及其改进模型在医学图像领域应用十分成功，其中unet 模型使得unet模型的深度可以自适应于数据集规模，改进效果明显，能够精准得到细胞区域。传统的图像分割方法如基于区域的分水岭算法应用十分广泛，在各项任务中取得了不错的效果。其改进算法如基于标记的分水岭算法通过融入预处理步骤来限制允许存在的区域数目，优化了原始分水岭算法对于存在噪声以及梯度不规则图像极易过度分割的问题。但在骨髓细胞的核质分割任务中，仍存在分割区域多，细胞核区域难以确定的问题，难以得到精准的细胞核区域。


技术实现要素：

4.本发明的第一目的在于克服现有技术的缺点与不足，提供一种骨髓细胞形态学特征的提取方法，突破传统基于标记的分水岭算法在进行单一骨髓细胞图像进行核质分割时，细胞核区域判定困难的问题，根据骨髓细胞形态学检验知识预先得到不同类型骨髓细胞的细胞核容积占比信息，提高细胞核区域判定的准确性，从而提取更精准的形态学特征。
5.本发明的第二目的在于提供一种骨髓细胞形态学特征的提取系统。
6.本发明的第一目的通过下述技术方案实现：一种骨髓细胞形态学特征的提取方法，包括以下步骤：
7.1)获取经化学试剂染色后得到的骨髓细胞图像，并进行图像预处理；
8.2)利用unet 模型对骨髓细胞图像进行分割，得到若干个单一骨髓细胞图像；
9.3)结合医学诊断知识对基于标记的分水岭算法进行改进，根据细胞核数值描述对单一骨髓细胞图像进行细胞核分割，获取精准的细胞核区域；
10.4)提取单一骨髓细胞图像的形态学特征，并进行特征选择，筛选得到重要性高的形态学特征。
11.进一步，在步骤1)中，所述化学试剂为骨髓细胞瑞氏染色过程中所需的瑞氏染料；所述骨髓细胞图像为在骨髓涂片镜检后，显微镜摄像头拍摄而得的显微图像；所述图像预处理包括：
①
灰度化处理，将色彩图像转化为灰度图像；
②
滤波处理，通过滤波函数去除噪声；
③
边缘检测，通过边缘检测算子得到骨髓细胞图像中的边缘信息。
12.进一步，在步骤2)中，所述unet 模型是一种图像分割算法，其在原始unet模型基础上，提出自适应不同规模数据集的unet体系结构。
13.进一步，所述步骤3)包括如下步骤：
14.3.1)根据骨髓细胞形态学检验知识，获取骨髓细胞各系别、各生长阶段的不同类型细胞的细胞核在细胞中的容积占比信息，并转化为可量化计算的数值描述；数值描述是容积占比参数的集合，包含骨髓细胞类型和其容积占比参数的一一对应关系，表示为areas＝{type_i:areainfo_i,i＝1,
…
,n}，其中，areas表示数值描述，type_i表示某一种类型的骨髓细胞，areainfo_i表示type_i对应的容积占比参数的数值描述，n表示骨髓细胞类型的范围阈值；
15.3.2)根据步骤3.1)所得的数值描述改进基于标记的分水岭算法，对单一骨髓细胞图像进行细胞核区域判定；基于标记的分水岭算法在对单一骨髓细胞图像进行核质分割时，存在分割区域多的问题，且由于不同类型细胞形态差异大，细胞核区域的容积不尽相同，单一的容积占比参数无法适用于所有类型的骨髓细胞，使得细胞核区域的判定问题较困难；因此，在对不同类型的单一骨髓细胞图像进行核质分割时，使用对应的容积占比参数判定细胞核所在区域，从而提高细胞核区域判定的准确性；
16.3.3)输出单一骨髓细胞图像的细胞核像素总数和细胞像素总数。
17.进一步，在步骤4)中，所述形态学特征包括形状特征、纹理特征和颜色特征，所述特征选择为对所提取的形态学特征进行特征筛选操作。
18.本发明的第二目的通过下述技术方案实现：一种骨髓细胞形态学特征的提取系统，包括：
19.骨髓细胞图像获取和预处理模块，用于获取经化学试剂染色后得到的骨髓细胞图像，并进行图像预处理；
20.骨髓细胞图像分割模块，利用unet 模型对骨髓细胞图像进行分割，得到若干个单一骨髓细胞图像；
21.单一骨髓细胞图像细胞核分割模块，结合医学诊断知识对基于标记的分水岭算法进行改进，根据细胞核数值描述对单一骨髓细胞图像进行细胞核分割，获取更精准的细胞核区域；
22.形态学特征提取和筛选模块，提取单一骨髓细胞图像的形态学特征，并进行特征选择，筛选得到重要性高的形态学特征。
23.进一步，在骨髓细胞图像获取和预处理模块中，所述化学试剂为骨髓细胞瑞氏染色过程中所需的瑞氏染料；所述骨髓细胞图像为在骨髓涂片镜检后，显微镜摄像头拍摄而得的显微图像；所述图像预处理包括：
①
灰度化处理，将色彩图像转化为灰度图像；
②
滤波处理，通过滤波函数去除噪声；
③
边缘检测，通过边缘检测算子得到骨髓细胞图像中的边缘信息。
24.进一步，在骨髓细胞图像分割模块中，所述unet 模型是一种图像分割算法，其在
原始unet模型基础上，提出自适应不同规模数据集的unet体系结构。
25.进一步，所述单一骨髓细胞图像细胞核分割模块执行以下操作：
26.根据骨髓细胞形态学检验知识，获取骨髓细胞各系别、各生长阶段的不同类型细胞的细胞核在细胞中的容积占比信息，并转化为可量化计算的数值描述；数值描述是容积占比参数的集合，包含骨髓细胞类型和其容积占比参数的一一对应关系，表示为areas＝{type_i:areainfo_i,i＝1,
…
,n}，其中，areas表示数值描述，type_i表示某一种类型的骨髓细胞，areainfo_i表示type_i对应的容积占比参数的数值描述，n表示骨髓细胞类型的范围阈值；
27.根据所得的数值描述改进基于标记的分水岭算法，对单一骨髓细胞图像进行细胞核区域判定；基于标记的分水岭算法在对单一骨髓细胞图像进行核质分割时，存在分割区域多的问题，且由于不同类型细胞形态差异大，细胞核区域的容积不尽相同，单一的容积占比参数无法适用于所有类型的骨髓细胞，使得细胞核区域的判定问题较困难；因此，在对不同类型的单一骨髓细胞图像进行核质分割时，使用对应的容积占比参数判定细胞核所在区域，从而提高细胞核区域判定的准确性；
28.输出单一骨髓细胞图像的细胞核像素总数和细胞像素总数。
29.进一步，在形态学特征提取和筛选模块中，所述形态学特征包括形状特征、纹理特征和颜色特征，所述特征选择为对所提取的形态学特征进行特征筛选操作。
30.本发明与现有技术相比，具有如下优点与有益效果：
31.1、本发明中结合骨髓细胞形态学检验知识，将骨髓细胞各系别、各生长阶段的不同类型细胞的细胞核在细胞中的容积占比信息并转化为可量化计算的数值描述，作为确定细胞核区域的先验知识。
32.2、本发明对基于标记的分水岭算法进行改进，考虑不同类型的骨髓细胞的细胞核区域差异，从而得到更精准的细胞核区域。
附图说明
33.图1为本发明方法逻辑流程示意图。
34.图2为本发明系统的架构图。
具体实施方式
35.下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
36.实施例1
37.如图1所示，本实施例公开了一种骨髓细胞形态学特征的提取方法，其包括以下步骤：
38.1)获取经化学试剂染色后得到的骨髓细胞图像，并进行图像预处理。其中，所述化学试剂为骨髓细胞瑞氏染色过程中所需的瑞氏染料；所述骨髓细胞图像为在骨髓涂片镜检后，显微镜摄像头拍摄而得的显微图像；所述图像预处理包括：
①
灰度化处理，将色彩图像转化为灰度图像；
②
滤波处理，通过滤波函数去除噪声；
③
边缘检测，通过边缘检测算子得到骨髓细胞图像中的边缘信息。本实施例中共获取得到1000张染色后的骨髓细胞图像，均
为rgb图像格式。首先对1000张图像进行灰度化处理，得到1000张灰度图像；接着使用高斯滤波器对1000张灰度图像进行处理，高斯内核大小设置为7*7。最后利用canny算子进行边缘检测，边缘检测的最大最小阈值分别为200和60。
39.2)利用unet 模型对骨髓细胞图像进行分割，得到若干个单一骨髓细胞图像。其中，所述unet 模型是一种图像分割算法，其在原始unet模型基础上，提出自适应不同规模数据集的unet体系结构。因其结构简单，对小样本数据集友好，更适用于医学图像领域的图像分割任务。本实施例中采用基于pytorch框架的unet 模型，将1000张预处理后的图像按照8:2的比例划分为训练集和测试集，通过将训练集输入到模型中，得到训练后的unet 模型，接着用训练好的模型对测试集进行预测，得到若干个单一骨髓细胞图像。
40.3)结合医学诊断知识对基于标记的分水岭算法进行改进，根据细胞核数值描述对单一骨髓细胞图像进行细胞核分割，获取精准的细胞核区域，包括如下步骤：
41.3.1)根据骨髓细胞形态学检验知识，获取骨髓细胞各系别、各生长阶段的不同类型细胞的细胞核在细胞中的容积占比信息，并转化为可量化计算的数值描述。数值描述是容积占比参数的集合，包含骨髓细胞类型和其容积占比参数的一一对应关系，可表示为areas＝{type_i:areainfo_i,i＝1,
…
,n}，其中，areas表示数值描述，type_i表示某一种类型的骨髓细胞，areainfo_i表示type_i对应的容积占比参数的数值描述，n表示骨髓细胞类型的范围阈值。本实施例中骨髓细胞的形态学检验知识如表1所示：
42.表1：骨髓细胞形态学知识
43.序号骨髓细胞类型名称系别细胞核大小描述1早幼红细胞红系胞核大小一般大于细胞大小的3/42中幼红细胞红系胞核大小一般大于细胞大小的一半3晚幼红细胞红系胞核大小一般小于细胞大小的一半4原始粒细胞粒系胞质量较少5中性中幼粒细胞粒系胞核大小一般大于细胞大小的一半6淋巴细胞淋巴系胞核大小一般大于细胞大小的一半7单核细胞单核系胞核大小一般大于细胞大小的一半8浆细胞浆系胞核大小约等于细胞大小的一半
44.将表1转化为数值描述，即数值描述＝{早幼红细胞：0.75《x《1；中幼红细胞：0.5《x《1；晚幼红细胞：0《x《0.5；原始粒细胞：0.8《x《1；中性中幼粒细胞0.5《x《1；淋巴细胞：0.5《x《1；单核细胞：0.5《x《1；浆细胞：0.4《x《0.6}。其中，x表示容积占比参数。
45.3.2)根据步骤3.1)所得的数值描述改进基于标记的分水岭算法，对单一骨髓细胞图像进行细胞核区域判定。基于标记的分水岭算法在对单一骨髓细胞图像进行核质分割时，存在分割区域多的问题，且由于不同类型细胞形态差异大，细胞核区域的容积不尽相同，单一的容积占比参数无法适用于所有类型的骨髓细胞，使得细胞核区域的判定问题较困难。因此在对不同类型的单一骨髓细胞图像进行核质分割时，使用对应的容积占比参数判定细胞核所在区域，从而提高细胞核区域判定的准确性。
46.3.3)输出单一骨髓细胞图像的细胞核像素总数和细胞像素总数。
47.4)提取单一骨髓细胞图像的形态学特征，并进行特征选择，筛选得到重要性高的形态学特征。其中，所述形态学特征包括形状特征，纹理特征和颜色特征。所述特征选择为
对所提取的形状特征进行特征筛选操作。本实施例中共提取形态学特征143个，其中形状特征包括细胞容积(细胞像素总数)，细胞核容积(细胞核像素总数)，细胞核与细胞质的容积比，细胞核与细胞的容积比，细胞核最小外接矩形容积，细胞核圆度，细胞核矩形度等共14个；纹理特征共102个，颜色特征共27个。将提取得到的形态学特征使用主成分分析进行特征选择，挑选累计贡献率》0.99的主成分，共计16个，作为最终的特征。
48.实施例2
49.本实施例公开了一种骨髓细胞形态学特征的提取系统，如图2所示，其包括以下功能模块：
50.骨髓细胞图像获取和预处理模块，用于获取经化学试剂染色后得到的骨髓细胞图像，并进行图像预处理；
51.骨髓细胞图像分割模块，利用unet 模型对骨髓细胞图像进行分割，得到若干个单一骨髓细胞图像；
52.单一骨髓细胞图像细胞核分割模块，结合医学诊断知识对基于标记的分水岭算法进行改进，根据细胞核数值描述对单一骨髓细胞图像进行细胞核分割，获取更精准的细胞核区域；
53.形态学特征提取和筛选模块，提取单一骨髓细胞图像的形态学特征，并进行特征选择，筛选得到重要性较高的形态学特征。
54.进一步，在骨髓细胞图像获取和预处理模块中，所述化学试剂为骨髓细胞瑞氏染色过程中所需的瑞氏染料；所述骨髓细胞图像为在骨髓涂片镜检后，显微镜摄像头拍摄而得的显微图像；所述图像预处理包括：
①
灰度化处理，将色彩图像转化为灰度图像；
②
滤波处理，通过滤波函数去除噪声；
③
边缘检测，通过边缘检测算子得到骨髓细胞图像中的边缘信息。
55.进一步，在骨髓细胞图像分割模块中，所述unet 模型是一种图像分割算法，其在原始unet模型基础上，提出自适应不同规模数据集的unet体系结构，因其结构简单，对小样本数据集友好，更适用于医学图像领域的图像分割任务。
56.进一步，所述单一骨髓细胞图像细胞核分割模块执行以下操作：
57.根据骨髓细胞形态学检验知识，获取骨髓细胞各系别、各生长阶段的不同类型细胞的细胞核在细胞中的容积占比信息，并转化为可量化计算的数值描述；数值描述是容积占比参数的集合，包含骨髓细胞类型和其容积占比参数的一一对应关系，表示为areas＝{type_i:areainfo_i,i＝1,
…
,n}，其中，areas表示数值描述，type_i表示某一种类型的骨髓细胞，areainfo_i表示type_i对应的容积占比参数的数值描述，n表示骨髓细胞类型的范围阈值；
58.根据所得的数值描述改进基于标记的分水岭算法，对单一骨髓细胞图像进行细胞核区域判定；基于标记的分水岭算法在对单一骨髓细胞图像进行核质分割时，存在分割区域多的问题，且由于不同类型细胞形态差异大，细胞核区域的容积不尽相同，单一的容积占比参数无法适用于所有类型的骨髓细胞，使得细胞核区域的判定问题较困难；因此，在对不同类型的单一骨髓细胞图像进行核质分割时，使用对应的容积占比参数判定细胞核所在区域，从而提高细胞核区域判定的准确性；
59.输出单一骨髓细胞图像的细胞核像素总数和细胞像素总数。
60.进一步，在形态学特征提取和筛选模块中，所述形态学特征包括形状特征、纹理特征和颜色特征，所述特征选择为对所提取的形态学特征进行特征筛选操作。
61.综上所述，在采用以上方案后，本发明为骨髓细胞的形态特征提取提供了新的方法与系统，结合骨髓细胞形态学检验知识，改进了基于标记的分水岭算法，得到更为精准的骨髓细胞形态学特征，有效提高骨髓细胞形态学检查效率，具有实际利用价值，值得推广。
62.上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。