基于塔形PCA的乳腺钼靶X线图像块特征提取方法.pdf

摘要
申请专利号：	CN201410437632.4	申请日：	2014.08.30
公开号：	CN104182755A	公开日：	2014.12.03
当前法律状态：	授权	有效性：	有权
法律详情：	授权\|\|\|实质审查的生效IPC(主分类):G06K 9/46申请日:20140830\|\|\|公开
IPC分类号：	G06K9/46; G06T7/00; A61B6/03	主分类号：	G06K9/46
申请人：	西安电子科技大学
发明人：	李洁; 王颖; 刘璐; 高锐; 逄敏; 焦志成; 王斌; 路文; 李圣喜; 张琪
地址：	710071 陕西省西安市太白南路2号
优先权：
专利代理机构：	陕西电子工业专利中心 61205	代理人：	田文英;王品华
PDF下载：	PDF下载

内容摘要

本发明公开了一种基于塔形PCA的乳腺钼靶X线图像块特征提取方法，主要解决现有技术提取的特征中不包含肿块中间密度大边缘密度小这一特征的缺点。其实现步骤是：(1)预处理；(2)构成塔形结构；(3)获得各图像层的灰度特征向量；(4)训练各图像层灰度特征的特征空间；(5)获得各图像层的主成分特征；(6)获得基于塔形PCA的乳腺钼靶X线图像块特征。本发明能使乳腺钼靶X线图像块的特征表示更为鲁棒，更为有效的表示图像特征，提高乳腺钼靶X线摄影图像中肿块区域检测的准确率，从而辅助放射科医生进行临床诊断。

权利要求书

权利要求书
1.  一种基于塔形PCA的乳腺钼靶X线图像块特征提取方法，包括如下步骤：
(1)预处理：
(1a)采用中值滤波的方法，对乳腺钼靶X线摄影图像进行去噪处理；
(1b)对去噪后的乳腺钼靶X线摄影图像进行下5采样，得到下采样后的乳腺钼靶X线图像；
(1c)对下采样后的乳腺钼靶X线图像，从其上边缘向下裁剪40行，从其下边缘向上裁剪40行，从其左边缘向右裁剪10列，从其右边缘向左裁剪10列，得到剪裁后的乳腺钼靶X线摄影图像；
(1d)采用最大类间方差法，对裁剪后的乳腺钼靶X线图像进行二值化处理，获得乳腺组织区域；
(1e)设置滑动窗窗口的大小为100×100个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为100×100个像素乳腺钼靶X线图像块；
(2)构成塔形结构：
(2a)将乳腺钼靶X线图像块的第34行、第65行、第34列及第65列所围成的正方形区域作为第一层图像层；
(2b)由第一层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第二层图像层；
(2c)由第二层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第三层图像层；
(2d)由第三层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第四层图像层；
(2e)由第四层图像层向上、下、左、右四个方向延伸至填充整个乳腺钼靶X线图像块，将延伸的区域作为第五层图像层；
(3)获得各图像层的灰度特征向量：
采用灰度特征提取方法，逐层提取乳腺钼靶X线图像块的灰度特征，获得乳腺钼靶X线图像块各图像层的灰度特征向量；
(4)训练各图像层灰度特征的特征空间：
从乳腺影像数字数据库DDSM数据库中，选取500幅乳腺肿块图像作为训练图像，采用主成分分析PCA方法，训练得到各图像层灰度特征的特征空间：
(4a)按照步骤(2)的塔形结构划分方法，分别将500幅训练图像中每幅训练图像划分为5层图像层；
(4b)采用灰度特征提取方法，逐层提取500幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的500个灰度特征向量按列组成第一层层灰度特征矩阵G1，将第二层所提取的500个灰度特征向量按列组成第二层灰度特征矩阵G2，将第三层所提取的500个灰度特征向量按列组成第三层灰度特征矩阵G3，将第四层所提取的500个灰度特征向量按列组成第四层灰度特征矩阵G4，将第五层所提取的500个灰度特征向量按列组成第五层灰度特征矩阵G5；
(4c)采用主成分分析PCA方法，分别提取各层灰度特征矩阵的主成分，得到由20个主成分组成的第一层灰度特征的特征空间D1，由16个主成分组成的第二层灰度特征的特征空间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成分组成的第四层灰度特征的特征空间D4，由2个主成分组成的第五层灰度特征的特征空间D5；
(5)获得各图像层的主成分特征：
(5a)用乳腺钼靶X线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间D1的转置，得到乳腺钼靶X线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至20维；
(5b)用乳腺钼靶X线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间D2的转置，得到乳腺钼靶X线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至16维；
(5c)用乳腺钼靶X线图像块第三层图像层的灰度特征向量左乘第三层灰度特征的特征空间D3的转置，得到乳腺钼靶X线图像块第三个图像层的主成分特征，将第三个图像层灰度特征维数降低至8维；
(5d)用乳腺钼靶X线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间D4的转置，得到乳腺钼靶X线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至4维；
(5e)用乳腺钼靶X线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间D5的转置，得到乳腺钼靶X线图像块第五个图像层的主成分特征，将第五个图像层灰度特征维数降低至2维；
(6)获得基于塔形PCA的乳腺钼靶X线图像块特征：
将获得的乳腺钼靶X线图像块5个图像层的主成分特征，按第一层至第五层的顺序依次首尾相连，获得基于塔形主成分分析PCA的乳腺钼靶X线图像块特征向量。

2.  根据权利要求1所述的基于塔形PCA的乳腺钼靶X线图像块特征提取方法，其特征在于，步骤(1a)所述的中值滤波方法的步骤如下：
第一步，将中值滤波器的滑动窗设置为3×3个像素的正方形窗；
第二步，用正方形窗沿着乳腺钼靶X线摄影图像行的方向逐像素滑动，在每一次滑动期间内，将正方形窗中的所有像素的灰度值，按照由小到大的顺序进行排序，选取排序结果的中间值，替代正方形窗中心位置像素的灰度值。

3.  根据权利要求1所述的基于塔形PCA的乳腺钼靶X线图像块特征提取方法，其特征在于，步骤(1b)所述的下5采样的步骤如下：
第一步，将下采样的采样间隔设置为5；
第二步，在去噪后的乳腺钼靶X线摄影图像中，每隔5个像素保留一个像素，将保留的所有像素构成下采样后的乳腺钼靶X线图像。

4.  根据权利要求1所述的基于塔形PCA的乳腺钼靶X线图像块特征提取方法，其特征在于，步骤(1d)所述的最大类间方差法，按如下步骤进行：
第一步，对裁剪后的乳腺钼靶X线图像中所有像素的灰度值求平均值，得到裁剪后的乳腺钼靶X线图像的平均灰度值u；
第二步，在裁剪后的乳腺钼靶X线图像中所有像素的灰度值中，选取最小值与最大值范围之间的任意一个灰度值，作为目标与背景的分割阈值t；
第三步，按照下式，计算目标与背景的类间方差：
G＝w1×(u1-u)2+w2×(u2-u)2
其中，G表示裁剪后的乳腺钼靶X线图像目标与背景的类间方差，w1表示裁剪后的乳腺钼靶X线图像中灰度值大于分割阈值t的目标像素个数与乳腺钼靶X线图像总像素个数之比，u1表示裁剪后的乳腺钼靶X线图像目标像素的平均灰度值，u表示裁剪后的乳腺钼靶X线图像的平均灰度值，w2表示裁剪后的乳腺钼靶X线图像中灰度值小于或等于分割阈值t的背景像素个数与乳腺钼靶X线图像总像素个数之比，u2表示裁剪后的乳腺钼靶X线图像背景像素的平均灰度值；
第四步，遍历目标与背景的分割阈值t的所有取值，找寻类间方差G最大时，分割阈值t所对应的取值，将该值作为最佳分割阈值；
第五步，提取裁剪后的乳腺钼靶X线图像所有灰度值大于最佳分割阈值t的像素，构成乳腺组织区域。

5.  根据权利要求1所述的基于塔形PCA的乳腺钼靶X线图像块特征提取方法，其特征在于，步骤(3)、步骤(4b)所述的灰度特征提取方法的步骤如下：
第一步，逐列提取乳腺图像区域的灰度值；
第二步，组合提取的各列灰度值，得到乳腺图像区域的灰度特征列向量。

6.  根据权利要求1所述的基于塔形PCA的乳腺钼靶X线图像块特征提取方法，其特征在于，步骤(4c)所述的主成分分析PCA方法的步骤如下：
第一步，按照下式，构建第k层的灰度特征矩阵Gk：
Gk＝[gk1,gk2,…gki,…gk500]
其中，Gk表示第k层灰度特征矩阵，gki表示第i幅训练图像在第k层的灰度特征列向量；
第二步，按照下式，计算第k层的灰度特征向量平均值Ψk：
Ψk=1500Σi=1500gki]]>
其中，Ψk表示第k层的灰度特征向量平均值，i表示训练图像的标号，gki表示第i幅训练图像在第k层的灰度特征列向量；
第三步，按照下式，计算得到每一幅训练图像在第k层的灰度特征列向量与第k层的灰度特征向量平均值的差值：
dki＝gki-Ψk
其中，dki表示第i幅训练图像在第k层的灰度特征列向量与第k层的灰度特征向量平均值的差值，Ψk表示第k层的灰度特征向量平均值，i表示训练图像的标号，xi表示第i幅训练图像在该层的灰度特征列向量；
第四步，按照下式，构建第k层的差值矩阵Ak：
Ak＝[dk1,dk2,…dki,…dk500]
其中，Ak表示第k层的差值矩阵，i表示训练图像的标号，dki表示第i幅训练图像在第k层的灰度特征列向量与该层的灰度特征向量平均值的差值；
第五步，按照下式，构建第k层的协方差矩阵Ck：
Ck=1500AkTAk]]>
其中，Ck表示第k层的协方差矩阵，Ak表示第k层的差值矩阵，AkT表示第k层的差值矩阵Ak的转置；
第六步，计算第k层的协方差矩阵Ck的特征值λki和正交归一化特征向量uki，选择前pk个最大特征值对应的特征向量，按照下式，得到第k层灰度特征的特征空间：
Dk=[uk1,uk2,...,ukpk]]]>
其中，Dk表示第k层灰度特征的特征空间，表示前pk个最大特征值对应的特征向量。

说明书

说明书基于塔形PCA的乳腺钼靶X线图像块特征提取方法
技术领域
本发明属于图像处理技术领域。更进一步涉及在医学影像处理技术领域中一种基于塔形主成分分析(Principal Component Analysis，PCA)的乳腺钼靶X线图像块灰度特征提取方法。本发明是根据乳腺肿块图像中间密度大边缘密度小分布情况，对乳腺图像块进行分层特征提取，同时对分层提取的灰度特征分别进行主成分分析，从而提高乳腺钼靶X线摄影图像中肿块区域检测的准确率。本发明可用于临床医疗病变区域检测，提高检出率，降低检测假阳率，辅助放射科医生进行临床诊断。
背景技术
目前，在临床医疗诊断中运用的图像特征灰度直观特征、灰度统计特征、变换域特征、代数特征。在图像处理中，灰度特征是最直观也是最全面表达图像的特征，可以很好地反应图像的灰度分布信息。灰度特征是最容易提取的特征，被广泛运用到图像的特征提取中，很多算法都是对灰度进行抽取得到的。为了对图像有更深层次、更本质的描述，去除特征冗余性，很多特征是在对灰度进行主成分分析(PCA)后获得的。
浙江大学申请的专利“基于主元分析和支持向量机的人脸检测方法”(申请号：201110446113.0，公开号：CN102592145A)中公开了一种对输入图像区域进行主成分分析的特征提取方法。该方法先将彩色图像转换为大小统一的灰度图，再对处理后的图像进行主成分分析(PCA)，同时用支持向量机对分析的结果进行分类。该方法存在的不足是，直接将灰度图直接拉成列向量进行主成分分析(PCA)，没有将图像的分布信息和空间信息融合进去，不能表征图像之间空间与位置上的关系。
李男的学位论文“基于PCA的乳腺肿块辅助诊断方法研究”(沈阳工业大学，2013年)公开了一种降低图像特征维数的方法。该方法利用主成分分析(PCA)的方法降低医学图像特征维数，去除特征之间的相关性，使得处理速度得到提高而不影响实验精度。但是，该方法仍然存在的不足是，选取的特征是直接将感兴趣区域拉成列向量，没有将组织的分布特征和空间特征融合进去，由此造成组织信息表征不完全，影响了乳腺肿块的诊断准确性。
马萌的学位论文“基于Directionlet和稀疏表示的乳腺X线图像检测和增强”(硕士学位论文，西安电子科技大学，2012年)公开了一种将“本征脸”特征运用到乳腺图像中的特征提取方法。该方法先提取样本图像块的“本征脸”特征，构建字典，然后直接获取测试图像块的“本征脸”特征，对特征进行分类，以特征在字典中的稀疏表示系数作为分类评判标准。该方法存在的不足是：特征维数太高，存在信息冗余，没有突出整体图像块中间密度大边缘密度小的重要性，分类准确率不高。
中山大学申请的专利“一种基于结构主元分析的人脸识别方法”(申请号：200610124229.1，公开号：CN100423020C)，公开了一种用二维主成分分析(2DPCA)计算每个图像块的主元和主元特征的方法。该方法先将几何规整后的图像分割为30个互不重叠且大小一致的图像块，然后计算图像块的结构主元，再由图像块左乘各个结构主元得到图像块的结构主元特征。该方法的不足是：图像块的划分方式对图像块内包含的内容定义不明确，不能完全运用在乳腺特征中。
发明内容
本发明的目的在于克服上述现有技术的不足，提出一种基于塔形PCA的乳腺钼靶X线图像块特征提取方法，以使图像特征表示更为鲁棒，提高乳腺钼靶X线摄影图像的肿块区域检出率。
实现本发明的技术思路是，通过滑动窗提取乳腺钼靶X线摄影图像中乳腺部分的图像块，将乳腺钼靶X线图像块由内而外划分为5层，构成一个塔形结构，分层提取乳腺钼靶X线图像块的灰度特征，同时对乳腺钼靶X线图像块各层灰度特征进行分层降维。在分层降维的过程中，乳腺钼靶X线图像块第一层灰度特征降维后的特征维数最高，随着层数变高，降维后的特征维数变少，使符合乳腺钼靶X线图像块中间的灰度信息量多，边缘信息量少的特点，将乳腺钼靶X线图像块中间密度大边缘密度小的特征表征出来，达到更为准确地表征乳腺钼靶X线图像块的灰度特征，更为合理地加入乳腺肿块密度分布特征，更为鲁棒地表示乳腺钼靶X线图像块，提高乳腺钼靶X线摄影图像病变区域检出率的目的。
为实现上述目的，本发明包括如下主要步骤：
(1)预处理：
(1a)采用中值滤波的方法，对乳腺钼靶X线摄影图像进行去噪处理；
(1b)对去噪后的乳腺钼靶X线摄影图像进行下5采样，得到下采样后的乳腺钼靶X线图像；
(1c)对下采样后的乳腺钼靶X线图像，从其上边缘向下裁剪40行，从其下边缘向上裁剪40行，从其左边缘向右裁剪10列，从其右边缘向左裁剪10列，得到剪裁后的乳腺钼靶X线摄影图像；
(1d)采用最大类间方差法，对裁剪后的乳腺钼靶X线图像进行二值化处理，获得乳腺组织区域；
(1e)设置滑动窗窗口的大小为100×100个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为100×100个像素乳腺钼靶X线图像块；
(2)构成塔形结构：
(2a)将乳腺钼靶X线图像块的第34行、第65行、第34列及第65列所围成的正方形区域作为第一层图像层；
(2b)由第一层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第二层图像层；
(2c)由第二层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第三层图像层；
(2d)由第三层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第四层图像层；
(2e)由第四层图像层向上、下、左、右四个方向延伸至填充整个乳腺钼靶X线图像块，将延伸的区域作为第五层图像层；
(3)获得各图像层的灰度特征向量：
采用灰度特征提取方法，逐层提取乳腺钼靶X线图像块的灰度特征，获得乳腺钼靶X线图像块各图像层的灰度特征向量；
(4)训练各图像层灰度特征的特征空间：
从乳腺影像数字数据库DDSM数据库中，选取500幅乳腺肿块图像作为训练图像，采用主成分分析PCA方法，训练得到各图像层灰度特征的特征空间：
(4a)按照步骤(2)的塔形结构划分方法，分别将500幅训练图像中每幅训练图像划分为5层图像层；
(4b)采用灰度特征提取方法，逐层提取500幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的500个灰度特征向量按列组成第一层层灰度特征矩阵G1，将第二层所提取的500个灰度特征向量按列组成第二层灰度特征矩阵G2，将第三层所提取的500个灰度特征向量按列组成第三层灰度特征矩阵G3，将第四层所提取的500个灰度特征向量按列组成第四层灰度特征矩阵G4，将第五层所提取的500个灰度特征向量按列组成第五层灰度特征矩阵G5；
(4c)采用主成分分析PCA方法，分别提取各层灰度特征矩阵的主成分，得到由20个主成分组成的第一层灰度特征的特征空间D1，由16个主成分组成的第二层灰度特征的特征空间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成分组成的第四层灰度特征的特征空间D4，由2个主成分组成的第五层灰度特征的特征空间D5；
(5)获得各图像层的主成分特征：
(5a)用乳腺钼靶X线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间D1的转置，得到乳腺钼靶X线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至20维；
(5b)用乳腺钼靶X线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间D2的转置，得到乳腺钼靶X线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至16维；
(5c)用乳腺钼靶X线图像块第三层图像层的灰度特征向量左乘第三层灰度特征的特征空间D3的转置，得到乳腺钼靶X线图像块第三个图像层的主成分特征，将第三个图像层灰度特征维数降低至8维；
(5d)用乳腺钼靶X线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间D4的转置，得到乳腺钼靶X线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至4维；
(5e)用乳腺钼靶X线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间D5的转置，得到乳腺钼靶X线图像块第五个图像层的主成分特征，将第五个图像层灰度特征维数降低至2维；
(6)获得基于塔形PCA的乳腺钼靶X线图像块特征：
将获得的乳腺钼靶X线图像块5个图像层的主成分特征，按第一层至第五层的顺序依次首尾相连，获得基于塔形主成分分析PCA的乳腺钼靶X线图像块特征向量。
本发明与现有方法相比具有如下优点：
第一，由于本发明采用了按照构造的塔形结构进行分层降维的方式，克服了现有技术主成分分析方法直接提取特征主成分方法固定的缺点，使得本发明具有提取的特征符合肿块中间灰度信息量多，边缘信息量少的特点，获得的基于塔形PCA的灰度特征更为合理。
第二，由于本发明采用了分层提取乳腺钼靶X线图像块的灰度特征的方式，克服了现有技术灰度特征提取方法固定、不考虑肿块空间分布特征的缺点，使得本发明提取的乳腺钼靶X线图像块的灰度特征更为精确，提取的灰度特征包含乳腺肿块密度分布特征更为完整，提高了组织信息表征的完全性。
第三，由于本发明利用主成分分析PCA对提取的各层灰度特征进行分层降维，克服了现有技术提取的特征维数太高，使得本发明获得的基于塔形PCA的灰度特征去除了特征之间的相关性和冗余信息，突出了整体乳腺钼靶X线图像块中间密度大边缘密度小的重要性，提高了乳腺钼靶X线摄影图像中肿块的检出率。
第四，由于本发明构建塔形结构的方式与肿块延伸方式相关，克服了现有技术特征提取方法中不包含肿块由中间向外延伸生长特征的缺点，提高了特征提取方式对乳腺图像的适应性。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图1，对本发明实现的步骤作进一步的详细描述。
步骤1，预处理。
对乳腺钼靶X线摄影图像进行预处理，经滑动窗获取得到的乳腺钼靶X线图像块，乳腺钼靶X线图像块的宽度为100个像素，乳腺钼靶X线图像块的高度为100个像素。
对乳腺钼靶X线摄影图像进行预处理的方法按如下步骤进行：
第一步，采用中值滤波的方法，对乳腺钼靶X线摄影图像进行去噪处理：将中值滤波器的滑动窗设置为3×3个像素的正方形窗，用正方形窗沿着乳腺钼靶X线摄影图像行的方向逐像素滑动，在每一次滑动期间内，将正方形窗中的所有像素的灰度值，按照由小到大的顺序进行排序，选取排序结果的中间值，替代正方形窗中心位置像素的灰度值。
第二步，对去噪后的乳腺钼靶X线摄影图像进行下5采样，得到下采样后的乳腺钼靶X线图像：首先，将下采样的采样间隔设置为5；然后，在去噪后的乳腺钼靶X线摄影图像中每隔5个像素保留一个像素，将保留的所有像素构成下采样后的乳腺钼靶X线图像得到下采样后的乳腺钼靶X线图像。
第三步，对下采样后的乳腺钼靶X线图像，从其上边缘向下裁剪40行，从其下边缘向上裁剪40行，从其左边缘向右裁剪10列，从其右边缘向左裁剪10列，得到剪裁后的乳腺钼靶X线摄影图像。
第四步，采用最大类间方差法，对裁剪后的乳腺钼靶X线图像进行二值化处理，获得乳腺组织区域，具体实施方法如下：
对裁剪后的乳腺钼靶X线图像中所有像素的灰度值求平均值，得到裁剪后的乳腺钼靶X线图像的平均灰度值u。
在裁剪后的乳腺钼靶X线图像中所有像素的灰度值中，选取最小值与最大值范围之间的任意一个灰度值，作为目标与背景的分割阈值t。
按照下式，计算目标与背景的类间方差：
G＝w1×(u1-u)2+w2×(u2-u)2
其中，G表示裁剪后的乳腺钼靶X线图像目标与背景的类间方差，w1表示裁剪后的乳腺钼靶X线图像中灰度值大于分割阈值t的目标像素个数与乳腺钼靶X线图像总像素个数之比，u1表示裁剪后的乳腺钼靶X线图像目标像素的平均灰度值，u表示裁剪后的乳腺钼靶X线图像的平均灰度值，w2表示裁剪后的乳腺钼靶X线图像中灰度值小于或等于分割阈值t的背景像素个数与乳腺钼靶X线图像总像素个数之比，u2表示裁剪后的乳腺钼靶X线图像背景像素的平均灰度值。
遍历目标与背景的分割阈值t的所有取值，找寻类间方差G最大时，分割阈值t所对应的取值，将该值作为最佳分割阈值。
提取裁剪后的乳腺钼靶X线图像所有灰度值大于最佳分割阈值t的像素，构成乳腺组织区域。
第五步，设置滑动窗窗口的大小为100×100个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为100×100个像素乳腺钼靶X线图像块。
步骤2，构成塔形结构。
将乳腺钼靶X线图像块由中心向外划分为5层图像层，构成一个塔形结构：
第一步，将乳腺钼靶X线图像块的第34行、第65行、第34列及第65列所围成的正方形区域作为第一层图像层。
第二步，由第一层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第二层图像层。
第三步，由第二层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第三层图像层。
第四步，由第三层图像层向上、下、左、右四个方向各延伸8个像素，将延伸的区域作为第四层图像层。
第五步，由第四层图像层向上、下、左、右四个方向延伸至填充整个乳腺钼靶X线图像块，将延伸的区域作为第五层图像层。
步骤3，获得各图像层的灰度特征向量。
采用灰度特征提取方法，逐层提取乳腺钼靶X线图像块的灰度特征，获得乳腺钼靶X线图像块各图像层的灰度特征向量。
灰度特征提取方法按如下步骤进行：
第一步，逐列提取构成乳腺钼靶X线图像块第k层图像层的乳腺图像区域的灰度值。
第二步，组合提取的各列灰度值，得到构成乳腺钼靶X线图像块第k层图像层的乳腺图像区域的灰度特征列向量。
步骤4，训练各图像层灰度特征的特征空间。
从乳腺影像数字数据库DDSM数据库中，选取500幅乳腺肿块图像作为训练图像，采用主成分分析PCA方法，训练得到各图像层灰度特征的特征空间。
第一步，按照步骤2的塔形结构划分方法，分别将500幅训练图像中每幅训练图像划分为5层图像层。
第二步，采用灰度特征提取方法，逐层提取500幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的500个灰度特征向量按列组成第一层层灰度特征矩阵G1，将第二层所提取的500个灰度特征向量按列组成第二层灰度特征矩阵G2，将第三层所提取的500个灰度特征向量按列组成第三层灰度特征矩阵G3，将第四层所提取的500个灰度特征向量按列组成第四层灰度特征矩阵G4，将第五层所提取的500个灰度特征向量按列组成第五层灰度特征矩阵G5。
采用灰度特征提取方法逐层提取每幅训练图像各图像层的灰度特征向量的具体实施方法如下：
逐列提取构成第i幅训练图像第k层图像层的乳腺图像区域的灰度值。
组合提取的各列灰度值，得到构成第i幅训练图像第k层图像层的乳腺图像区域的灰度特征列向量。
第三步，采用主成分分析PCA方法，分别提取各层灰度特征矩阵的主成分，得到由20个主成分组成的第一层灰度特征的特征空间D1，由16个主成分组成的第二层灰度特征的特征空间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成分组成的第四层灰度特征的特征空间D4，由2个主成分组成的第五层灰度特征的特征空间D5。
采用主成分分析PCA方法分别提取各层灰度特征矩阵的主成分组成各层灰度特征的特征空间的具体实施方法如下：
本步骤中第二步得到的第k层灰度特征矩阵Gk的构成结构如下：
Gk＝[gk1,gk2,…gki,…gk500]
其中，Gk表示第k层灰度特征矩阵，gki表示第i幅训练图像在第k层的灰度特征列向量。
按照下式，计算第k层的灰度特征向量平均值Ψk：
Ψk=1500Σi=1500gki]]>
其中，Ψk表示第k层的灰度特征向量平均值，i表示训练图像的标号，gki表示第i幅训练图像在第k层的灰度特征列向量。
按照下式，计算得到每一幅训练图像在第k层的灰度特征列向量与第k层的灰度特征向量平均值的差值：
dki＝gki-Ψk
其中，dki表示第i幅训练图像在第k层的灰度特征列向量与第k层的灰度特征向量平均值的差值，Ψk表示第k层的灰度特征向量平均值，i表示训练图像的标号，xi表示第i幅训练图像在该层的灰度特征列向量。
按照下式，构建第k层的差值矩阵Ak：
Ak＝[dk1,dk2,…dki,…dk500]
其中，Ak表示第k层的差值矩阵，i表示训练图像的标号，dki表示第i幅训练图像在第k层的灰度特征列向量与该层的灰度特征向量平均值的差值。
按照下式，构建第k层的协方差矩阵Ck：
Ck=1500AkTAk]]>
其中，Ck表示第k层的协方差矩阵，Ak表示第k层的差值矩阵，AkT表示第k层的差值矩阵Ak的转置。
计算第k层的协方差矩阵Ck的特征值λki和正交归一化特征向量uki，根据贡献率选择前pk个最大特征值对应的特征向量，其中贡献率是指所选取的前pk个最大特征值之和与所有特征值之和的比值。将贡献率设定为90％，使满足前pk个最大特征值之和与所有特征值之和的比值大于等于90％的条件。按照下式，得到第k层灰度特征的特征空间：
Dk=[uk1,uk2,...,ukpk]]]>
其中，Dk表示第k层灰度特征的特征空间，表示前pk个最大特征值对应的特征向量，pk的值在第一层至第五层中分别对应于20、16、8、4、2。
步骤5，获得各图像层的主成分特征。
分别将乳腺钼靶X线图像块各图像层的灰度特征向量映射到相应图像层灰度特征的特征空间上，分层提取乳腺钼靶X线图像块各图像层灰度特征的主成分特征，达到分层降维的目的。
第一步，用乳腺钼靶X线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间D1的转置，得到乳腺钼靶X线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至20维。
第二步，用乳腺钼靶X线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间D2的转置，得到乳腺钼靶X线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至16维。
第三步，用乳腺钼靶X线图像块第三层图像层的灰度特征向量左乘第三层灰度特征的特征空间D3的转置，得到乳腺钼靶X线图像块第三个图像层的主成分特征，将第三个图像层灰度特征维数降低至8维。
第四步，用乳腺钼靶X线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间D4的转置，得到乳腺钼靶X线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至4维。
第五步，用乳腺钼靶X线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间D5的转置，得到乳腺钼靶X线图像块第五个图像层的主成分特征，将第五个图像层灰度特征维数降低至2维。
步骤6，获得基于塔形PCA的乳腺钼靶X线图像块特征。
将获得的乳腺钼靶X线图像块5个图像层的主成分特征，按第一层至第五层的顺序依次首尾相连构成一个50维的列向量，获得基于塔形主成分分析PCA的乳腺钼靶X线图像块特征向量。
本发明的效果可以通过以下仿真实验做进一步的说明。
1.仿真条件：
本发明是在中央处理器为Intel(R)Core i3-21003.10GHZ、内存4G、WINDOWS 7操作系统上，运用MATLAB软件进行的仿真。
2.仿真内容：
本发明对大量乳腺钼靶X线摄影图像进行了测试实验，实验图像来自于乳腺影像数字数据库DDSM数据库中任意选取的1191幅乳腺钼靶X线摄影图像。
本发明以检出率和假阳性率为指标对方法性能进行评测，仿真对比了不同特征提取方法对乳腺钼靶X线摄影图像进行肿块检测的检出率和假阳性率。仿真实验中，选取1191幅乳腺钼靶X线摄影图像中的200幅乳腺钼靶X线摄影图像作为肿块检测的训练集，其中包含肿块的乳腺钼靶X线摄影图像和正常的乳腺钼靶X线摄影图像各100幅，其余991幅乳腺钼靶X线摄影图像作为肿块检测的测试集，其中包含恶性肿块的乳腺钼靶X线摄影图像有508幅，包含良性肿块的乳腺钼靶X线摄影图像有483幅。对肿块检测的训练集图像进行特征提取形成字典，对肿块检测的测试图像进行滑动窗提取乳腺钼靶X线图像块而后进行特征提取，运用稀疏表示分类器得到第一次稀疏表示检测结果，再将标记出的肿块区域进行滑动窗合并，运用区域生长算法提取感兴趣区域，进行第二次稀疏表示检测，最终标记疑似肿块区域，达到检测的目的，得到检出率和假阳性率。检出率表征二次稀疏表示检测对肿块检测的测试集中正确检测出的肿块结果个数与肿块检测的测试集中所有肿块个数之比。假阳性率表征二次稀疏表示检测在肿块检测的测试集中平均每幅用作肿块检测的图像标记出的假阳个数。在保证检测模型相同的情况下，检出率和假阳性率直接反映了图像特征的优劣，因此检出率越高，假阳性率越低，表明图像特征越好。对于辅助放射科医生进行医学诊断的检测技术而言，两种指标中提高检出率更为重要。
本发明的仿真实验是对991幅乳腺钼靶X线摄影图像中，包含恶性肿块的乳腺钼靶X线摄影图像和包含良性肿块的乳腺钼靶X线摄影图像，进行肿块检测的检出率和假阳性率的进行了对比。对比实验结果如表1所示。
表1中的第一列中cancer表示数据库中包含恶性肿块的乳腺钼靶X线摄影图像的不同case，第三列中benign表示包含良性肿块的乳腺钼靶X线摄影图像的不同case，第二列和第五列中的分数表示检出率，其中，分子表示检测出的肿块数，分母表示该case中所有肿块数之和，第三列和第六列中的数字表示假阳性率。
由表1可见，恶性肿块的检出率更高，比良性肿块的检出率高8个百分点，同时恶性肿块的假阳性率相对良性肿块较低，平均每幅乳腺钼靶X线摄影图像中少0.7个。在对实验结果的分析中，发现本发明对于中心亮度区域较为敏感，这是因为本发明突出了乳腺肿块区域中间密度大边缘密度小的特点，所以能够更好地描述恶性乳腺肿块特征，由此获得在检出率上优于良性肿块的效果，验证了算法的特性。
表1本发明特征提取方法的乳腺肿块检测结果

本发明的第二个仿真实验是对多种特征提取方法，对钼靶X线摄影图像肿块检测的检出率进行对比，对比的多种特征提取方法包括对灰度特征直接降维的方法、提取方向梯度直方图(Histogram of Oriented Gradient，HOG)特征并进行降维的方法和本发明方法。对比实验结果如表2所示。
表2灰度特征直接降维以及HOG特征降维与本发明检出率对比
Cases直接降维HOG特征降维本发明cancer_0119/4134/4137/41cancer_0230/6248/6254/62cancer_0535/5043/5043/50cancer_0616/3123/3128/31平均检出率100/184＝0.54148/184＝0.80162/184＝0.88
表2中第一列中cancer表示数据库中包含恶性肿块的乳腺钼靶X线摄影图像的不同case，第二列中的分数表示利用灰度特征直接降维的方法进行肿块检测的检出率，第三列中的分数表示利用方向梯度直方图HOG特征降维的方法进行肿块检测的检出率，第四列中的分数表示用本发明的方法进行肿块检测的检出率。
由表2可见，本发明的检出率是三种方法中最高的，这是因为本发明突出了乳腺肿块区域中间密度大边缘密度小的特点，所以能够更好地描述乳腺肿块特征，由此获得在检出率上优于其他特征提取方法的效果，进一步验证了算法的先进性。
综上所述，本发明能有效地表示乳腺钼靶X线图像块的特征，提高了乳腺肿块的检出率，从而辅助放射科医生进行医学诊断。

资源描述

《基于塔形PCA的乳腺钼靶X线图像块特征提取方法.pdf》由会员分享，可在线阅读，更多相关《基于塔形PCA的乳腺钼靶X线图像块特征提取方法.pdf（15页珍藏版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 104182755 A (43)申请公布日 2014.12.03 CN 104182755 A (21)申请号 201410437632.4 (22)申请日 2014.08.30 G06K 9/46(2006.01) G06T 7/00(2006.01) A61B 6/03(2006.01) (71)申请人西安电子科技大学地址 710071 陕西省西安市太白南路 2 号 (72)发明人李洁王颖刘璐高锐逄敏焦志成王斌路文李圣喜张琪 (74)专利代理机构陕西电子工业专利中心 61205 代理人田文英王品华 (54) 发明名称基于塔形 PCA。

2、的乳腺钼靶 X 线图像块特征提取方法 (57) 摘要本发明公开了一种基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，主要解决现有技术提取的特征中不包含肿块中间密度大边缘密度小这一特征的缺点。其实现步骤是： (1)预处理； (2)构成塔形结构； (3) 获得各图像层的灰度特征向量； (4) 训练各图像层灰度特征的特征空间； (5) 获得各图像层的主成分特征； (6)获得基于塔形PCA的乳腺钼靶X线图像块特征。本发明能使乳腺钼靶X 线图像块的特征表示更为鲁棒，更为有效的表示图像特征，提高乳腺钼靶 X 线摄影图像中肿块区域检测的准确率，从而辅助放射。

3、科医生进行临床诊断。 (51)Int.Cl. 权利要求书 4 页说明书 9 页附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请权利要求书4页说明书9页附图1页 (10)申请公布号 CN 104182755 A CN 104182755 A 1/4 页 2 1. 一种基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，包括如下步骤： (1) 预处理： (1a) 采用中值滤波的方法，对乳腺钼靶 X 线摄影图像进行去噪处理； (1b) 对去噪后的乳腺钼靶 X 线摄影图像进行下 5 采样，得到下采样后的乳腺钼靶 X 线图像； (1c) 对下采样后。

4、的乳腺钼靶 X 线图像，从其上边缘向下裁剪 40 行，从其下边缘向上裁剪 40 行，从其左边缘向右裁剪 10 列，从其右边缘向左裁剪 10 列，得到剪裁后的乳腺钼靶 X 线摄影图像； (1d) 采用最大类间方差法，对裁剪后的乳腺钼靶 X 线图像进行二值化处理，获得乳腺组织区域； (1e) 设置滑动窗窗口的大小为 100100 个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为 100100 个像素乳腺钼靶 X 线图像块； (2) 构成塔形结构： (2a) 将乳腺钼靶 X 线图像块的第 34 行、第 65 行、第 34 列及第 65。

5、列所围成的正方形区域作为第一层图像层； (2b) 由第一层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第二层图像层； (2c) 由第二层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第三层图像层； (2d) 由第三层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第四层图像层； (2e) 由第四层图像层向上、下、左、右四个方向延伸至填充整个乳腺钼靶 X 线图像块，将延伸的区域作为第五层图像层； (3) 获得各图像层的灰度特征向量：采用灰度特征提取方法，逐层提取乳腺钼靶X线图像块。

6、的灰度特征，获得乳腺钼靶X线图像块各图像层的灰度特征向量； (4) 训练各图像层灰度特征的特征空间：从乳腺影像数字数据库 DDSM 数据库中，选取 500 幅乳腺肿块图像作为训练图像，采用主成分分析 PCA 方法，训练得到各图像层灰度特征的特征空间： (4a) 按照步骤 (2) 的塔形结构划分方法，分别将 500 幅训练图像中每幅训练图像划分为 5 层图像层； (4b) 采用灰度特征提取方法，逐层提取 500 幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的500个灰度特征向量按列组成第一层层灰度特征矩阵G1，将第二层所提取的 500 个灰度。

7、特征向量按列组成第二层灰度特征矩阵 G2，将第三层所提取的 500 个灰度特征向量按列组成第三层灰度特征矩阵 G3，将第四层所提取的 500 个灰度特征向量按列组成第四层灰度特征矩阵 G4，将第五层所提取的 500 个灰度特征向量按列组成第五层灰度特征矩阵 G5； (4c) 采用主成分分析 PCA 方法，分别提取各层灰度特征矩阵的主成分，得到由 20 个主成分组成的第一层灰度特征的特征空间D1，由16个主成分组成的第二层灰度特征的特征空权利要求书 CN 104182755 A 2 2/4 页 3 间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成。

8、分组成的第四层灰度特征的特征空间 D4，由 2 个主成分组成的第五层灰度特征的特征空间 D5； (5) 获得各图像层的主成分特征： (5a) 用乳腺钼靶 X 线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间 D1的转置，得到乳腺钼靶 X 线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至 20 维； (5b) 用乳腺钼靶 X 线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间 D2的转置，得到乳腺钼靶 X 线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至 16 维； (5c) 用乳腺钼靶 X 线图像块第三层图。

9、像层的灰度特征向量左乘第三层灰度特征的特征空间 D3的转置，得到乳腺钼靶 X 线图像块第三个图像层的主成分特征，将第三个图像层灰度特征维数降低至 8 维； (5d) 用乳腺钼靶 X 线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间 D4的转置，得到乳腺钼靶 X 线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至 4 维； (5e) 用乳腺钼靶 X 线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间 D5的转置，得到乳腺钼靶 X 线图像块第五个图像层的主成分特征，将第五个图像层灰度特征维数降低至 2 维； (6) 获得基于塔。

10、形 PCA 的乳腺钼靶 X 线图像块特征：将获得的乳腺钼靶 X 线图像块 5 个图像层的主成分特征，按第一层至第五层的顺序依次首尾相连，获得基于塔形主成分分析 PCA 的乳腺钼靶 X 线图像块特征向量。 2. 根据权利要求 1 所述的基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，其特征在于，步骤 (1a) 所述的中值滤波方法的步骤如下：第一步，将中值滤波器的滑动窗设置为 33 个像素的正方形窗；第二步，用正方形窗沿着乳腺钼靶 X 线摄影图像行的方向逐像素滑动，在每一次滑动期间内，将正方形窗中的所有像素的灰度值，按照由小到大的顺序进行排序，选取排序。

11、结果的中间值，替代正方形窗中心位置像素的灰度值。 3. 根据权利要求 1 所述的基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，其特征在于，步骤 (1b) 所述的下 5 采样的步骤如下：第一步，将下采样的采样间隔设置为 5 ；第二步，在去噪后的乳腺钼靶 X 线摄影图像中，每隔 5 个像素保留一个像素，将保留的所有像素构成下采样后的乳腺钼靶 X 线图像。 4. 根据权利要求 1 所述的基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，其特征在于，步骤 (1d) 所述的最大类间方差法，按如下步骤进行：第一步，对裁剪后的乳腺钼靶 X 线图像中所有。

12、像素的灰度值求平均值，得到裁剪后的乳腺钼靶 X 线图像的平均灰度值 u ；第二步，在裁剪后的乳腺钼靶 X 线图像中所有像素的灰度值中，选取最小值与最大值范围之间的任意一个灰度值，作为目标与背景的分割阈值 t ；第三步，按照下式，计算目标与背景的类间方差：权利要求书 CN 104182755 A 3 3/4 页 4 G w1(u1-u)2+w2(u2-u)2 其中， G表示裁剪后的乳腺钼靶X线图像目标与背景的类间方差， w1表示裁剪后的乳腺钼靶 X 线图像中灰度值大于分割阈值 t 的目标像素个数与乳腺钼靶 X 线图像总像素个数之比， u1 表示裁剪后的乳腺钼靶。

13、 X 线图像目标像素的平均灰度值， u 表示裁剪后的乳腺钼靶 X 线图像的平均灰度值， w2 表示裁剪后的乳腺钼靶 X 线图像中灰度值小于或等于分割阈值 t 的背景像素个数与乳腺钼靶 X 线图像总像素个数之比， u2 表示裁剪后的乳腺钼靶 X 线图像背景像素的平均灰度值；第四步，遍历目标与背景的分割阈值 t 的所有取值，找寻类间方差 G 最大时，分割阈值 t 所对应的取值，将该值作为最佳分割阈值；第五步，提取裁剪后的乳腺钼靶X线图像所有灰度值大于最佳分割阈值t的像素，构成乳腺组织区域。 5. 根据权利要求 1 所述的基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法。

14、，其特征在于，步骤 (3)、步骤 (4b) 所述的灰度特征提取方法的步骤如下：第一步，逐列提取乳腺图像区域的灰度值；第二步，组合提取的各列灰度值，得到乳腺图像区域的灰度特征列向量。 6. 根据权利要求 1 所述的基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，其特征在于，步骤 (4c) 所述的主成分分析 PCA 方法的步骤如下：第一步，按照下式，构建第 k 层的灰度特征矩阵 Gk： Gk gk1,gk2,gki,gk500 其中， Gk表示第 k 层灰度特征矩阵， gki表示第 i 幅训练图像在第 k 层的灰度特征列向量；第二步，按照下式，。

15、计算第 k 层的灰度特征向量平均值 k：其中， k表示第 k 层的灰度特征向量平均值， i 表示训练图像的标号， gki表示第 i 幅训练图像在第 k 层的灰度特征列向量；第三步，按照下式，计算得到每一幅训练图像在第k层的灰度特征列向量与第k层的灰度特征向量平均值的差值： dki gki-k 其中， dki表示第 i 幅训练图像在第 k 层的灰度特征列向量与第 k 层的灰度特征向量平均值的差值， k表示第 k 层的灰度特征向量平均值， i 表示训练图像的标号， xi表示第 i 幅训练图像在该层的灰度特征列向量；第四步，按照下式，构建第 k 层的差值矩阵 Ak： Ak。

16、 dk1,dk2,dki,dk500 其中， Ak表示第 k 层的差值矩阵， i 表示训练图像的标号， dki表示第 i 幅训练图像在第 k 层的灰度特征列向量与该层的灰度特征向量平均值的差值；第五步，按照下式，构建第 k 层的协方差矩阵 Ck：权利要求书 CN 104182755 A 4 4/4 页 5 其中， Ck表示第 k 层的协方差矩阵， Ak表示第 k 层的差值矩阵， AkT表示第 k 层的差值矩阵 Ak的转置；第六步，计算第 k 层的协方差矩阵 Ck的特征值 ki和正交归一化特征向量 uki，选择前 pk个最大特征值对应的特征向量，按照下式，得到第。

17、k 层灰度特征的特征空间：其中， Dk表示第 k 层灰度特征的特征空间，表示前 pk个最大特征值对应的特征向量。权利要求书 CN 104182755 A 5 1/9 页 6 基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法技术领域 0001 本发明属于图像处理技术领域。更进一步涉及在医学影像处理技术领域中一种基于塔形主成分分析 (Principal Component Analysis， PCA) 的乳腺钼靶 X 线图像块灰度特征提取方法。本发明是根据乳腺肿块图像中间密度大边缘密度小分布情况，对乳腺图像块进行分层特征提取，同时对分层提取的灰度特征分别进行主成。

18、分分析，从而提高乳腺钼靶 X 线摄影图像中肿块区域检测的准确率。本发明可用于临床医疗病变区域检测，提高检出率，降低检测假阳率，辅助放射科医生进行临床诊断。背景技术 0002 目前，在临床医疗诊断中运用的图像特征灰度直观特征、灰度统计特征、变换域特征、代数特征。在图像处理中，灰度特征是最直观也是最全面表达图像的特征，可以很好地反应图像的灰度分布信息。灰度特征是最容易提取的特征，被广泛运用到图像的特征提取中，很多算法都是对灰度进行抽取得到的。为了对图像有更深层次、更本质的描述，去除特征冗余性，很多特征是在对灰度进行主成分分析 (PCA) 后获得的。 000。

19、3 浙江大学申请的专利 “基于主元分析和支持向量机的人脸检测方法” ( 申请号： 201110446113.0，公开号： CN102592145A) 中公开了一种对输入图像区域进行主成分分析的特征提取方法。该方法先将彩色图像转换为大小统一的灰度图，再对处理后的图像进行主成分分析 (PCA)，同时用支持向量机对分析的结果进行分类。该方法存在的不足是，直接将灰度图直接拉成列向量进行主成分分析 (PCA)，没有将图像的分布信息和空间信息融合进去，不能表征图像之间空间与位置上的关系。 0004 李男的学位论文 “基于 PCA 的乳腺肿块辅助诊断方法研究” ( 沈阳工业大学， 2。

20、013 年 ) 公开了一种降低图像特征维数的方法。该方法利用主成分分析 (PCA) 的方法降低医学图像特征维数，去除特征之间的相关性，使得处理速度得到提高而不影响实验精度。但是，该方法仍然存在的不足是，选取的特征是直接将感兴趣区域拉成列向量，没有将组织的分布特征和空间特征融合进去，由此造成组织信息表征不完全，影响了乳腺肿块的诊断准确性。 0005 马萌的学位论文 “基于Directionlet和稀疏表示的乳腺X线图像检测和增强” (硕士学位论文，西安电子科技大学， 2012 年 ) 公开了一种将 “本征脸” 特征运用到乳腺图像中的特征提取方法。该方法先提取样本图像块的。

21、 “本征脸” 特征，构建字典，然后直接获取测试图像块的 “本征脸” 特征，对特征进行分类，以特征在字典中的稀疏表示系数作为分类评判标准。该方法存在的不足是：特征维数太高，存在信息冗余，没有突出整体图像块中间密度大边缘密度小的重要性，分类准确率不高。 0006 中山大学申请的专利 “一种基于结构主元分析的人脸识别方法” ( 申请号： 200610124229.1，公开号： CN100423020C)，公开了一种用二维主成分分析 (2DPCA) 计算每个图像块的主元和主元特征的方法。该方法先将几何规整后的图像分割为 30 个互不重叠且大小一致的图像块，然后计算。

22、图像块的结构主元，再由图像块左乘各个结构主元得到图说明书 CN 104182755 A 6 2/9 页 7 像块的结构主元特征。该方法的不足是：图像块的划分方式对图像块内包含的内容定义不明确，不能完全运用在乳腺特征中。发明内容 0007 本发明的目的在于克服上述现有技术的不足，提出一种基于塔形 PCA 的乳腺钼靶 X 线图像块特征提取方法，以使图像特征表示更为鲁棒，提高乳腺钼靶 X 线摄影图像的肿块区域检出率。 0008 实现本发明的技术思路是，通过滑动窗提取乳腺钼靶 X 线摄影图像中乳腺部分的图像块，将乳腺钼靶 X 线图像块由内而外划分为 5 层，构成一个塔。

23、形结构，分层提取乳腺钼靶 X 线图像块的灰度特征，同时对乳腺钼靶 X 线图像块各层灰度特征进行分层降维。在分层降维的过程中，乳腺钼靶 X 线图像块第一层灰度特征降维后的特征维数最高，随着层数变高，降维后的特征维数变少，使符合乳腺钼靶 X 线图像块中间的灰度信息量多，边缘信息量少的特点，将乳腺钼靶 X 线图像块中间密度大边缘密度小的特征表征出来，达到更为准确地表征乳腺钼靶 X 线图像块的灰度特征，更为合理地加入乳腺肿块密度分布特征，更为鲁棒地表示乳腺钼靶 X 线图像块，提高乳腺钼靶 X 线摄影图像病变区域检出率的目的。 0009 为实现上述目的，本发明包括如。

24、下主要步骤： 0010 (1) 预处理： 0011 (1a) 采用中值滤波的方法，对乳腺钼靶 X 线摄影图像进行去噪处理； 0012 (1b) 对去噪后的乳腺钼靶 X 线摄影图像进行下 5 采样，得到下采样后的乳腺钼靶 X 线图像； 0013 (1c) 对下采样后的乳腺钼靶 X 线图像，从其上边缘向下裁剪 40 行，从其下边缘向上裁剪 40 行，从其左边缘向右裁剪 10 列，从其右边缘向左裁剪 10 列，得到剪裁后的乳腺钼靶 X 线摄影图像； 0014 (1d) 采用最大类间方差法，对裁剪后的乳腺钼靶 X 线图像进行二值化处理，获得乳腺组织区域； 0015 。

25、(1e) 设置滑动窗窗口的大小为 100100 个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为 100100 个像素乳腺钼靶 X 线图像块； 0016 (2) 构成塔形结构： 0017 (2a) 将乳腺钼靶 X 线图像块的第 34 行、第 65 行、第 34 列及第 65 列所围成的正方形区域作为第一层图像层； 0018 (2b) 由第一层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第二层图像层； 0019 (2c) 由第二层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第三层图像。

26、层； 0020 (2d) 由第三层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第四层图像层； 0021 (2e) 由第四层图像层向上、下、左、右四个方向延伸至填充整个乳腺钼靶 X 线图像块，将延伸的区域作为第五层图像层； 0022 (3) 获得各图像层的灰度特征向量：说明书 CN 104182755 A 7 3/9 页 8 0023 采用灰度特征提取方法，逐层提取乳腺钼靶 X 线图像块的灰度特征，获得乳腺钼靶 X 线图像块各图像层的灰度特征向量； 0024 (4) 训练各图像层灰度特征的特征空间： 0025 从乳腺影像数字数据库。

27、 DDSM 数据库中，选取 500 幅乳腺肿块图像作为训练图像，采用主成分分析 PCA 方法，训练得到各图像层灰度特征的特征空间： 0026 (4a) 按照步骤 (2) 的塔形结构划分方法，分别将 500 幅训练图像中每幅训练图像划分为 5 层图像层； 0027 (4b) 采用灰度特征提取方法，逐层提取 500 幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的 500 个灰度特征向量按列组成第一层层灰度特征矩阵 G1，将第二层所提取的 500 个灰度特征向量按列组成第二层灰度特征矩阵 G2，将第三层所提取的 500 个灰度特征向量按列组成第三层灰度特征。

28、矩阵 G3，将第四层所提取的 500 个灰度特征向量按列组成第四层灰度特征矩阵 G4，将第五层所提取的 500 个灰度特征向量按列组成第五层灰度特征矩阵 G5； 0028 (4c) 采用主成分分析 PCA 方法，分别提取各层灰度特征矩阵的主成分，得到由 20 个主成分组成的第一层灰度特征的特征空间D1，由16个主成分组成的第二层灰度特征的特征空间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成分组成的第四层灰度特征的特征空间 D4，由 2 个主成分组成的第五层灰度特征的特征空间 D5； 0029 (5) 获得各图像层的主成分特征： 0030 (5a) 。

29、用乳腺钼靶 X 线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间 D1的转置，得到乳腺钼靶 X 线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至 20 维； 0031 (5b) 用乳腺钼靶 X 线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间 D2的转置，得到乳腺钼靶 X 线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至 16 维； 0032 (5c) 用乳腺钼靶 X 线图像块第三层图像层的灰度特征向量左乘第三层灰度特征的特征空间 D3的转置，得到乳腺钼靶 X 线图像块第三个图像层的主成分特征，将第三个图。

30、像层灰度特征维数降低至 8 维； 0033 (5d) 用乳腺钼靶 X 线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间 D4的转置，得到乳腺钼靶 X 线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至 4 维； 0034 (5e) 用乳腺钼靶 X 线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间 D5的转置，得到乳腺钼靶 X 线图像块第五个图像层的主成分特征，将第五个图像层灰度特征维数降低至 2 维； 0035 (6) 获得基于塔形 PCA 的乳腺钼靶 X 线图像块特征： 0036 将获得的乳腺钼靶 X 线图像块 5 个图像。

31、层的主成分特征，按第一层至第五层的顺序依次首尾相连，获得基于塔形主成分分析 PCA 的乳腺钼靶 X 线图像块特征向量。 0037 本发明与现有方法相比具有如下优点： 0038 第一，由于本发明采用了按照构造的塔形结构进行分层降维的方式，克服了现有技术主成分分析方法直接提取特征主成分方法固定的缺点，使得本发明具有提取的特征符说明书 CN 104182755 A 8 4/9 页 9 合肿块中间灰度信息量多，边缘信息量少的特点，获得的基于塔形 PCA 的灰度特征更为合理。 0039 第二，由于本发明采用了分层提取乳腺钼靶 X 线图像块的灰度特征的方式，克服了现有技术。

32、灰度特征提取方法固定、不考虑肿块空间分布特征的缺点，使得本发明提取的乳腺钼靶 X 线图像块的灰度特征更为精确，提取的灰度特征包含乳腺肿块密度分布特征更为完整，提高了组织信息表征的完全性。 0040 第三，由于本发明利用主成分分析 PCA 对提取的各层灰度特征进行分层降维，克服了现有技术提取的特征维数太高，使得本发明获得的基于塔形 PCA 的灰度特征去除了特征之间的相关性和冗余信息，突出了整体乳腺钼靶 X 线图像块中间密度大边缘密度小的重要性，提高了乳腺钼靶 X 线摄影图像中肿块的检出率。 0041 第四，由于本发明构建塔形结构的方式与肿块延伸方式相关，克服了现有。

33、技术特征提取方法中不包含肿块由中间向外延伸生长特征的缺点，提高了特征提取方式对乳腺图像的适应性。附图说明 0042 图 1 为本发明的流程图。具体实施方式 0043 下面结合附图 1，对本发明实现的步骤作进一步的详细描述。 0044 步骤 1，预处理。 0045 对乳腺钼靶X线摄影图像进行预处理，经滑动窗获取得到的乳腺钼靶X线图像块，乳腺钼靶 X 线图像块的宽度为 100 个像素，乳腺钼靶 X 线图像块的高度为 100 个像素。 0046 对乳腺钼靶 X 线摄影图像进行预处理的方法按如下步骤进行： 0047 第一步，采用中值滤波的方法，对乳腺钼靶 X 线摄影图像进行去。

34、噪处理：将中值滤波器的滑动窗设置为33个像素的正方形窗，用正方形窗沿着乳腺钼靶X线摄影图像行的方向逐像素滑动，在每一次滑动期间内，将正方形窗中的所有像素的灰度值，按照由小到大的顺序进行排序，选取排序结果的中间值，替代正方形窗中心位置像素的灰度值。 0048 第二步，对去噪后的乳腺钼靶X线摄影图像进行下5采样，得到下采样后的乳腺钼靶X线图像：首先，将下采样的采样间隔设置为5 ；然后，在去噪后的乳腺钼靶X线摄影图像中每隔 5 个像素保留一个像素，将保留的所有像素构成下采样后的乳腺钼靶 X 线图像得到下采样后的乳腺钼靶 X 线图像。 0049 第三步，对。

35、下采样后的乳腺钼靶X线图像，从其上边缘向下裁剪40行，从其下边缘向上裁剪 40 行，从其左边缘向右裁剪 10 列，从其右边缘向左裁剪 10 列，得到剪裁后的乳腺钼靶 X 线摄影图像。 0050 第四步，采用最大类间方差法，对裁剪后的乳腺钼靶 X 线图像进行二值化处理，获得乳腺组织区域，具体实施方法如下： 0051 对裁剪后的乳腺钼靶 X 线图像中所有像素的灰度值求平均值，得到裁剪后的乳腺钼靶 X 线图像的平均灰度值 u。 0052 在裁剪后的乳腺钼靶 X 线图像中所有像素的灰度值中，选取最小值与最大值范围说明书 CN 104182755 A 9 5/9 页。

36、 10 之间的任意一个灰度值，作为目标与背景的分割阈值 t。 0053 按照下式，计算目标与背景的类间方差： 0054 G w1(u1-u)2+w2(u2-u)2 0055 其中， G表示裁剪后的乳腺钼靶X线图像目标与背景的类间方差， w1表示裁剪后的乳腺钼靶 X 线图像中灰度值大于分割阈值 t 的目标像素个数与乳腺钼靶 X 线图像总像素个数之比， u1表示裁剪后的乳腺钼靶X线图像目标像素的平均灰度值， u表示裁剪后的乳腺钼靶 X 线图像的平均灰度值， w2 表示裁剪后的乳腺钼靶 X 线图像中灰度值小于或等于分割阈值 t 的背景像素个数与乳腺钼靶 X 线图像总像素个数之比， u2。

37、表示裁剪后的乳腺钼靶 X 线图像背景像素的平均灰度值。 0056 遍历目标与背景的分割阈值 t 的所有取值，找寻类间方差 G 最大时，分割阈值 t 所对应的取值，将该值作为最佳分割阈值。 0057 提取裁剪后的乳腺钼靶 X 线图像所有灰度值大于最佳分割阈值 t 的像素，构成乳腺组织区域。 0058 第五步，设置滑动窗窗口的大小为 100100 个像素，在获得的乳腺组织区域中按行的方向滑动扫描，滑动窗每次滑动截取一个大小为 100100 个像素乳腺钼靶 X 线图像块。 0059 步骤 2，构成塔形结构。 0060 将乳腺钼靶 X 线图像块由中心向外划分为 5 层图像层。

38、，构成一个塔形结构： 0061 第一步，将乳腺钼靶 X 线图像块的第 34 行、第 65 行、第 34 列及第 65 列所围成的正方形区域作为第一层图像层。 0062 第二步，由第一层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第二层图像层。 0063 第三步，由第二层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第三层图像层。 0064 第四步，由第三层图像层向上、下、左、右四个方向各延伸 8 个像素，将延伸的区域作为第四层图像层。 0065 第五步，由第四层图像层向上、下、左、右四个方向延伸至。

39、填充整个乳腺钼靶 X 线图像块，将延伸的区域作为第五层图像层。 0066 步骤 3，获得各图像层的灰度特征向量。 0067 采用灰度特征提取方法，逐层提取乳腺钼靶 X 线图像块的灰度特征，获得乳腺钼靶 X 线图像块各图像层的灰度特征向量。 0068 灰度特征提取方法按如下步骤进行： 0069 第一步，逐列提取构成乳腺钼靶 X 线图像块第 k 层图像层的乳腺图像区域的灰度值。 0070 第二步，组合提取的各列灰度值，得到构成乳腺钼靶X线图像块第k层图像层的乳腺图像区域的灰度特征列向量。 0071 步骤 4，训练各图像层灰度特征的特征空间。 0072 从乳腺影像数字数据库。

40、 DDSM 数据库中，选取 500 幅乳腺肿块图像作为训练图像，采用主成分分析 PCA 方法，训练得到各图像层灰度特征的特征空间。说明书 CN 104182755 A 10 6/9 页 11 0073 第一步，按照步骤2的塔形结构划分方法，分别将500幅训练图像中每幅训练图像划分为 5 层图像层。 0074 第二步，采用灰度特征提取方法，逐层提取 500 幅训练图像中每幅训练图像各图像层的灰度特征向量，将第一层所提取的 500 个灰度特征向量按列组成第一层层灰度特征矩阵G1，将第二层所提取的500个灰度特征向量按列组成第二层灰度特征矩阵G2，将第三层所提取的。

41、500 个灰度特征向量按列组成第三层灰度特征矩阵 G3，将第四层所提取的 500 个灰度特征向量按列组成第四层灰度特征矩阵 G4，将第五层所提取的 500 个灰度特征向量按列组成第五层灰度特征矩阵 G5。 0075 采用灰度特征提取方法逐层提取每幅训练图像各图像层的灰度特征向量的具体实施方法如下： 0076 逐列提取构成第 i 幅训练图像第 k 层图像层的乳腺图像区域的灰度值。 0077 组合提取的各列灰度值，得到构成第 i 幅训练图像第 k 层图像层的乳腺图像区域的灰度特征列向量。 0078 第三步，采用主成分分析 PCA 方法，分别提取各层灰度特征矩阵的主成分，得到由。

42、 20 个主成分组成的第一层灰度特征的特征空间 D1，由 16 个主成分组成的第二层灰度特征的特征空间D2，由8个主成分组成的第三层灰度特征的特征空间D3，由4个主成分组成的第四层灰度特征的特征空间 D4，由 2 个主成分组成的第五层灰度特征的特征空间 D5。 0079 采用主成分分析 PCA 方法分别提取各层灰度特征矩阵的主成分组成各层灰度特征的特征空间的具体实施方法如下： 0080 本步骤中第二步得到的第 k 层灰度特征矩阵 Gk的构成结构如下： 0081 Gk gk1,gk2,gki,gk500 0082 其中， Gk表示第 k 层灰度特征矩阵， gki表示第 i 幅训。

43、练图像在第 k 层的灰度特征列向量。 0083 按照下式，计算第 k 层的灰度特征向量平均值 k： 0084 0085 其中， k表示第 k 层的灰度特征向量平均值， i 表示训练图像的标号， gki表示第 i 幅训练图像在第 k 层的灰度特征列向量。 0086 按照下式，计算得到每一幅训练图像在第 k 层的灰度特征列向量与第 k 层的灰度特征向量平均值的差值： 0087 dki gki-k 0088 其中， dki表示第 i 幅训练图像在第 k 层的灰度特征列向量与第 k 层的灰度特征向量平均值的差值， k表示第 k 层的灰度特征向量平均值， i 表示训练图像的标号， xi表示第。

44、 i 幅训练图像在该层的灰度特征列向量。 0089 按照下式，构建第 k 层的差值矩阵 Ak： 0090 Ak dk1,dk2,dki,dk500 0091 其中， Ak表示第 k 层的差值矩阵， i 表示训练图像的标号， dki表示第 i 幅训练图像在第 k 层的灰度特征列向量与该层的灰度特征向量平均值的差值。 0092 按照下式，构建第 k 层的协方差矩阵 Ck：说明书 CN 104182755 A 11 7/9 页 12 0093 0094 其中， Ck表示第 k 层的协方差矩阵， Ak表示第 k 层的差值矩阵， AkT表示第 k 层的差值矩阵 Ak的转置。 0095 计算。

45、第k层的协方差矩阵Ck的特征值ki和正交归一化特征向量uki，根据贡献率选择前pk个最大特征值对应的特征向量，其中贡献率是指所选取的前pk个最大特征值之和与所有特征值之和的比值。将贡献率设定为 90，使满足前 pk个最大特征值之和与所有特征值之和的比值大于等于 90的条件。按照下式，得到第 k 层灰度特征的特征空间： 0096 0097 其中， Dk表示第 k 层灰度特征的特征空间，表示前 pk个最大特征值对应的特征向量， pk的值在第一层至第五层中分别对应于 20、 16、 8、 4、 2。 0098 步骤 5，获得各图像层的主成分特征。 0099 分别将乳腺钼靶 X 线。

46、图像块各图像层的灰度特征向量映射到相应图像层灰度特征的特征空间上，分层提取乳腺钼靶 X 线图像块各图像层灰度特征的主成分特征，达到分层降维的目的。 0100 第一步，用乳腺钼靶 X 线图像块第一层图像层的灰度特征向量左乘第一层灰度特征的特征空间 D1的转置，得到乳腺钼靶 X 线图像块第一个图像层的主成分特征，将第一个图像层灰度特征维数降低至 20 维。 0101 第二步，用乳腺钼靶 X 线图像块第二层图像层的灰度特征向量左乘第二层灰度特征的特征空间 D2的转置，得到乳腺钼靶 X 线图像块第二个图像层的主成分特征，将第二个图像层灰度特征维数降低至 16 维。 0102。

47、第三步，用乳腺钼靶 X 线图像块第三层图像层的灰度特征向量左乘第三层灰度特征的特征空间 D3的转置，得到乳腺钼靶 X 线图像块第三个图像层的主成分特征，将第三个图像层灰度特征维数降低至 8 维。 0103 第四步，用乳腺钼靶 X 线图像块第四层图像层的灰度特征向量左乘第四层灰度特征的特征空间 D4的转置，得到乳腺钼靶 X 线图像块第四个图像层的主成分特征，将第四个图像层灰度特征维数降低至 4 维。 0104 第五步，用乳腺钼靶 X 线图像块第五层图像层的灰度特征向量左乘第五层灰度特征的特征空间 D5的转置，得到乳腺钼靶 X 线图像块第五个图像层的主成分特征，将第五。

48、个图像层灰度特征维数降低至 2 维。 0105 步骤 6，获得基于塔形 PCA 的乳腺钼靶 X 线图像块特征。 0106 将获得的乳腺钼靶 X 线图像块 5 个图像层的主成分特征，按第一层至第五层的顺序依次首尾相连构成一个 50 维的列向量，获得基于塔形主成分分析 PCA 的乳腺钼靶 X 线图像块特征向量。 0107 本发明的效果可以通过以下仿真实验做进一步的说明。 0108 1. 仿真条件： 0109 本发明是在中央处理器为 Intel(R)Core i3-21003.10GHZ、内存 4G、 WINDOWS 7 操作系统上，运用 MATLAB 软件进行的仿真。说明。

49、书 CN 104182755 A 12 8/9 页 13 0110 2. 仿真内容： 0111 本发明对大量乳腺钼靶 X 线摄影图像进行了测试实验，实验图像来自于乳腺影像数字数据库 DDSM 数据库中任意选取的 1191 幅乳腺钼靶 X 线摄影图像。 0112 本发明以检出率和假阳性率为指标对方法性能进行评测，仿真对比了不同特征提取方法对乳腺钼靶 X 线摄影图像进行肿块检测的检出率和假阳性率。仿真实验中，选取 1191 幅乳腺钼靶 X 线摄影图像中的 200 幅乳腺钼靶 X 线摄影图像作为肿块检测的训练集，其中包含肿块的乳腺钼靶X线摄影图像和正常的乳腺钼靶X线摄影图像各100幅，其余991 幅乳腺钼靶 X 线摄影图像作为肿块检测。

展开阅读全文