一种基于超声图像特征相关性的乳腺肿瘤分类方法、装置与流程

本发明涉及一种基于超声图像特征相关性的乳腺肿瘤分类方法、装置。

背景技术：

乳腺癌已经成为一种常见的女性的恶性肿瘤之一，在女性癌症的发病率的25％，死亡率占15％。近年来，我国女性乳腺癌的发病率逐年上升，占恶性肿瘤的23％，严重威胁着女性的身体健康和生命安全。

随着医学影像技术的不断发展，钼靶、超声、mri、ct等成像设备已经被逐渐应用于乳腺癌的检测工作中。钼靶是一种常用的乳腺癌影像学分析方法，然而钼靶图像对软组织的分辨力差，在致密性乳腺癌诊断精度较低，仅能达到30％-48％，超声图像作为一种常见的补充辅助分析方式之一，价格较低，且能够对致密性乳腺肿瘤的分类具有更好的表现，有利于减少不必要的细胞活检，已经能逐渐引起研究学者的关注。

目前临床上，仅依靠临床医生对影像学图像的观察进行乳腺肿瘤的分析工作，大量的影像学设备及患者会产生大量的乳腺图像，极大的增加医生的工作量。同时分析的结果也受临床医生的经验和主观判断影响，图像中的阴影等噪声也会影响医生诊断的精度，因此难以较快获得客观的诊断结果，特别针对于经验不足的临床医生容易造成误判等现象。

为了提高乳腺肿瘤图像分析的精度和效率，克服临床医生的主观判断的影响等因素，基于计算机的乳腺癌计算机辅助诊断势在必行。当前的乳腺癌计算机诊断方法主要通过提取超声图像中的不同的特征进行分类，主要包括纹理特征，形态学特征，基于模型的方法等。不同的特征对乳腺超声图像具有不同的表现，不同特征之间的相关性在乳腺肿瘤良性恶性的区分上也具有一定研究价值，然而上述方法都没有对特征之间的相关性进行分析。图像的特征及特征之间的相关性都是乳腺肿瘤分类的有利依据。

技术实现要素：

针对现有技术中存在的不足，本发明提供了一种基于超声图像特征相关性的乳腺肿瘤分类方法，该方法提出了在超声图像上对乳腺肿瘤进行分类，主要包括以下三个关键部分：图像预处理，特征相关性提取和乳腺肿瘤的分类。上述乳腺癌的诊断方法，对于特征相关性在超声图像上乳腺癌的诊断上提出了新的观念并取得较好的效果。

本发明的技术方案为：

一种基于超声图像特征相关性的乳腺肿瘤分类方法，包括：

对乳腺超声图像进行预处理去噪，提取目标兴趣区域；

基于目标兴趣区域中的每个像素点，提取底层图像特征；

根据所述底层图像特征，利用协方差矩阵提取特征相关性，采用基于数据挖掘的分类器对所述特征相关性进行分类，进而得到乳腺肿瘤的分类结果。

进一步的，所述基于数据挖掘的分类器采用knn分类器。

进一步的，所述提取目标兴趣区域包括：

在所述乳腺超声图像上进行标记，根据标记提取包含有乳腺肿瘤区域的最小外接矩形，作为目标兴趣区域。

进一步的，所述提取底层图像特征包括：

根据每个像素点的灰度及其在x方向和y方向的一阶导数和二阶导数，提取底层图像特征。

进一步的，根据所述底层图像特征，利用协方差矩阵提取特征相关性包括，根据以下公式提取协方差矩阵：

其中，n表示图像中的像素个数，fk表示该像素的底层特征，μ表示所有像素的特征均值，k表示像素的序号；协方差矩阵的对角线元素表示的是不同特征自身的方差，非对角线元素表示的是不同特征之间的协方差，非对角线元素若值为正值，表示对应两个特征为正相关，反之，若值为负，表示对应两个特征为负相关，如果为0，表示两个特征相互独立。

为了满足机器学习算法的需求，减少计算复杂度我们采用了l²ecm算法对协方差矩阵进行向量化，将协方差矩阵中的值按照顺序构成一个新的特征向量，采用该特征向量表示特征相关性。

采用l2ecm算法对协方差矩阵进行向量化包括：将所述协方差矩阵用c＝u∑u^t表示，其中c表示该协方差矩阵，u表示正交阵，其中对角线元素为协方差矩阵c的特征值，计算得logc＝ulog(∑)u^t，将协方差矩阵中的值按照顺序构成一个新的特征向量，即vlog(c)＝[vlog(c1),vlog(c2),...,vlog(cn)]，采用该特征向量表示特征相关性。

本发明还提出了一种计算机可读存储介质，其中存储有多条指令，所述指令由处理器加载并执行以下处理：

对乳腺超声图像进行预处理去噪，提取目标兴趣区域；

基于目标兴趣区域中的每个像素点，提取底层图像特征；

根据所述底层图像特征，利用协方差矩阵提取特征相关性，采用基于数据挖掘的分类器对所述特征相关性进行分类，进而得到乳腺肿瘤的分类结果。

本发明又提出了一种基于超声图像特征相关性的乳腺肿瘤分类装置，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令由处理器加载并执行以下处理：

对乳腺超声图像进行预处理去噪，提取目标兴趣区域；

基于目标兴趣区域中的每个像素点，提取底层图像特征；

根据所述底层图像特征，利用协方差矩阵提取特征相关性，采用基于数据挖掘的分类器对所述特征相关性进行分类，进而得到乳腺肿瘤的分类结果。

本发明的有益效果：

(1)本发明利用了协方差矩阵提取图像底层特征之间的相关性，能够有效的降低图像特征的维度。对于一个m*n的图像，假设基于图像的像素点，提取s维的特征，则图像特征尺度高达m*n*s，而利用协方差矩阵提取上述s维特征的相关性，协方差矩阵的大小仅为s*s。利用l2ecm进行向量化后得到的图像特征仅为(s²+s)/2，极大的降低了特征维度。

(2)本发明探究了图像的底层特征之间的相关性，并将其用于解决乳腺肿瘤的分类问题中。本发明提出的基于特征相关性的乳腺肿瘤分类方法已经取得了较好的结果。本发明将协方差矩阵同其他的统计方法(如均值，方差，均值+方差等方法)进行了比较，本发明的准确度，精确性，敏感性及auc值都远高于pca等其他方法。实验表明，特征相关性在乳腺肿瘤的分类中能取得更好的结果，其准确性。

附图说明

图1是乳腺良性肿瘤超声图像；

图2是乳腺恶性肿瘤超声图像；

图4是本发明乳腺肿瘤分类方法的流程图；

具体实施方式：

下面结合附图与实施例对本发明作进一步说明：

应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

本发明的一种典型实施例是一种基于超声图像特征相关性的乳腺肿瘤分类方法，包括以下步骤：

对乳腺超声图像进行预处理，由于中值滤波能够有效的消除噪声点，本发明选择中值滤波进行去噪。中值滤波的具体方案为：依次扫描图像中的像素点，取该像素点3*3邻域内灰度均值，并将其赋值给该像素点。

提取感兴趣区域：肿瘤区域是区分肿瘤良性恶性的关键区域，故本发明首先提取肿瘤区域，而不是利用整幅图像进行后续工作。

在提取感兴趣区域时，本发明根据医生标记的肿瘤区域，提取肿瘤区域的最小外接矩形，并将其作为感兴趣区域，如此即保留了肿瘤区域，也为后续的特征提取等工作提供方便。

然后基于目标兴趣区域中的每个像素点，为提取特征间相关性，首先需要提取底层特征。

分别对图像感兴趣区域中的每个像素提取灰度，灰度在x方向和y方向的一阶导数和二阶导数等五维特征，对于每个乳腺肿瘤图像中的像素点，其底层特征f(x,y)可表示为

其中，i(x,y)表示像素点(x,y)的灰度，表示像素点(x,y)的灰度在x方向上的一阶导数，即表示像素点的灰度在y方向上的一阶导数，即表示像素点的灰度在x方向上的二阶导数，即表示像素点的灰度在y方向上的二阶导数，即

然后基于上述提取的底层特征，提取特征相关性，并将其作为后续分类的依据。

本文采用协方差矩阵来提取底层特征之间的相关性。基于上述提出的底层特征f(x,y)，本发明根据以下公式提取协方差矩阵，

其中，n表示图像中的像素个数，zk表示该像素的特征值，μ表示所有像素的特征均值。

同时，为了满足机器学习算法的需求，减少计算复杂度我们采用了l²ecm算法对特征矩阵进行向量化。l²ecm算法的基本方法如下：

由于协方差矩阵c是一个半正定矩阵，任意一个半正定矩阵都可以用c＝u∑u^t表示,其中u表示正交阵，∑＝diag(λ1,...,λn)表示对角阵，其中对角线元素为协方差矩阵c的特征值，那么对于矩阵c，可以计算得logc＝ulog(∑)u^t。

l²ecm就是基于logc对协方差矩阵c进行向量化，即对协方差矩阵计算logc，将协方差矩阵中的值按照顺序构成一个新的特征向量，由于协方差矩阵是对称矩阵，故在向量化过程中只需取协方差矩阵的上三角或下三角区域，即vlog(c)＝[vlog(c1),vlog(c2),...,vlog(cn)]。

其中vlogc指的是将logc矩阵取上三角区域或下三角区域的值按顺序组成的向量，如图3所示。

最后采用基于数据挖掘的分类器，优选knn算法进行分类，实现乳腺肿瘤的分类。

其中knn算法的思想如下：

测试集中的两个样本x和y，其中，x＝[x1,x2,...,xn]，y＝[y1,y2,...,yn]分别表示x和y的特征向量。则两个样本之间的距离可以表示为，

通过在训练集中，寻求与其距离较小的k个样本，采用最大投票准则，确定最后的分类结果。

此外，本发明还提出了一种计算机可读存储介质，其中存储有多条指令，所述指令由处理器加载并执行以下处理：

对乳腺超声图像进行预处理去噪，提取目标兴趣区域；

基于目标兴趣区域中的每个像素点，提取底层图像特征；

根据所述底层图像特征，利用协方差矩阵提取特征相关性，采用基于数据挖掘的分类器对所述特征相关性进行分类，进而得到乳腺肿瘤的分类结果。

在上述存储介质上，本发明又提出了一种基于超声图像特征相关性的乳腺肿瘤分类装置，包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令由处理器加载并执行以下处理：

对乳腺超声图像进行预处理去噪，提取目标兴趣区域；

基于目标兴趣区域中的每个像素点，提取底层图像特征；

根据所述底层图像特征，利用协方差矩阵提取特征相关性，采用基于数据挖掘的分类器对所述特征相关性进行分类，进而得到乳腺肿瘤的分类结果。

本发明公开了一种基于超声图像特征相关性的乳腺癌诊断方法，该方法主要包括以下三个关键部分：图像预处理，特征相关性提取和乳腺肿瘤的分类。上述乳腺癌的诊断方法，对于特征相关性在超声图像上乳腺癌的诊断上提出了新的观念并取得较好的效果。

本发明将协方差矩阵同其他的统计方法(如均值，方差，均值+方差等方法)进行了比较，如图1个图2所示，对图2的恶性肿瘤按照本发明的步骤进行分析，准确率可达77.96％，敏感性为62.75％，特意性为83.70％，roc曲线下面积为0.682。可以看到本发明分类方法的准确度，精确性，敏感性及auc值都远高于pca等其他方法。本方法分类的该实验结果表明，乳腺图像特征的相关性在乳腺肿瘤分类分析中具有非常重要的功能，是实现肿瘤类疾病辅助诊断的重要保障。

以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。