基于主成分分析的ct图像去噪方法
【技术领域】
[0001] 本发明涉及医学技术领域，具体涉及一种基于主成分分析的CT图像去噪方法。
【背景技术】
[0002] 生物医学图像处理是数字图像处理的一个重要方面。很多生物医学信息是以图像 形式表现出来的，如CT图像，它使人类视觉从表面向内部延伸，人们可以通过它们来获取 人体内部器官在解剖形态、生物化学和生理功能上的有用信息。由于CT图像中的病灶部位 会在灰度和形状上相似，用肉眼不易分辨。不同能力和背景的人对同一幅医学图像往往会 得出不同的结果来。要相对图像有定量的评估就更是不可能的。因此，用计算机对CT图像 作后处理，首要的任务就是对获取的医学图像进行增强信噪比的工作。即滤处图像的噪声 和干扰，突出感兴趣区域或边缘从而为进一步分析（如图像分割、三维重建）和计算奠定基 础。对CT图像处理的主要目的是要提高CT图像的可懂度，即通过去除噪声，增强对比度， 根据医生的需要，尽可能突出有用的特征来，这样来改善图像的视觉效果。然而目前常见的 CT图像去噪方法存在的主要问题是，在滤除噪声的同时不能够较好的保存图像的边缘、细 节信息。
【发明内容】
[0003] 本发明实施例所要解决的技术问题在于，提供一种基于主成分分析的CT图像去 噪方法，能够准确的去除CT图像中的噪声，并且保留完整的图像边缘、细节信息。
[0004] 为解决上述技术问题，本发明实施例提供了一种基于主成分分析的CT图像去噪 方法，包括：
[0005] 构建一个可自适应调整大小的搜索窗口；
[0006] 以CT图像中一个包含噪声的像素为中心，构建矩形参考模块；
[0007] 在搜索窗口中，寻找与参考模块相似的测试模块，构成测试模块集；
[0008] 测试模块集映射到CT图像主成分的张量基底空间中；
[0009] 对包含噪声的张量基底进行削减；
[0010] 将测试模块映射回图像空间中，重建去除噪声后的CT图像。
[0011] 在本发明的其他方案中，搜索窗口为矩形窗口。
[0012] 在本发明的其他方案中，搜索窗口的大小可以根据不同窗口的特征自动调整。
[0013] 在本发明的其他方案中，参考模块是以CT图像中一个包含噪声的像素为中心的 矩形窗口，窗口的大小为固定值。
[0014] 在本发明的其他方案中，窗口大小为15x15。
[0015] 在本发明的其他方案中，两个模块相似的评测标准为两个模块中相对应位置的像 素的灰度值差异最小。
[0016] 在本发明的其他方案中，通过预设阈值，提取CT图像的主成分，构成张量基底空 间。
[0017] 在本发明的其他方案中，CT图像主成分的张量为三阶张量，参考模块的大小为一 阶和二阶，测试模块数目为三阶。
[0018] 在本发明的其他方案中，削减包含噪声的张量基底，保留包含CT图像主要信息的 成分。
[0019] 本发明能够准确的去除CT图像中的噪声，并且可以保留CT图像中的边缘和细节 信息。
【附图说明】
[0020] 图1是本发明所提出的基于主成分分析的CT图像去噪方法的流程图；
[0021] 图2是本发明所提出的基于主成分分析的CT图像去噪方法的搜索窗口构建示意 图。
【具体实施方式】
[0022] 下面结合附图和具体实施例对本发明作进一步详细描述，但不作为对本发明的限 定。
[0023] 图1是本发明的基于主成分分析的CT图像去噪方法的流程图，包括以下几个步 骤：
[0024] 步骤 S1，
[0025] 构建一个可自适应调整大小的搜索窗口。
[0026] 搜索窗口可以为矩形窗口，搜索窗口的大小可以根据不同窗口的特征自动调整。 例如，构建一个大小为LXL的搜索窗口，引入两个评估标准来描述不同大小搜索窗口的权 重：中位数绝对偏差和四分位间距绝对偏差，定义两个窗口 KXK和!^父匕且K > L，通过比 较[P+1，K]范围内所有窗口的权重可得最优搜索窗口 B。将所有比较窗口展成向量Id1 = [1^，132，···ΛΧι]，其中 i e (Ρ+1，Κ]，那么 MAD可表示为：
[0027] MAD (bj = median I bi-median (bj I
[0028] 其中，median(b,_)表示按窗口像素灰度平均值排序后的中位数。IQRAD可表示为：
[0030] 其中，％和&表示按窗口中像素灰度平均值排序后位于第一个和第三个四分位 数的灰度值。最优搜索窗口 B可通过下式得到：
[0032] 其中，median(X。）表示像素 X。灰度值，L = i为最优搜索窗口大小。
[0033] 步骤 S2，
[0034] 以CT图像中一个包含噪声的像素为中心，构建矩形参考模块。
[0035] 参考模块可以是以CT图像中一个包含噪声的像素为中心的矩形窗口，窗口的大 小为固定值，优选15x15。
[0036] 步骤 S3，
[0037] 在搜索窗口中，寻找与参考模块相似的测试模块，构成测试模块集。
[0038] 在大小为1^\1^的搜索窗口中包含（仏-?+1)2-1)个与参考模块\相同的观测模 块，采用块匹配法寻找与Xe相似的测试模块，定义与参考模块的相似性测度小于预设阈值 的观测模块为详细模块，其中相似性测度计算公式为：
[0040] 其中，PXP表示模块大小，X。(p)和X1 (p)表示X。和X i中第p个像素。
[0041] 上述中，两个模块相似的评测标准可以为两个模块中相对应位置的像素的灰度值 差异最小。当然，还可以采用其他评测标准。
[0042] 步骤 S4，
[0043] 测试模块集映射到CT图像主成分的张量基底空间中。
[0044] 在步骤S4,可以通过预设阈值，提取CT图像的主成分，构成张量基底空间。例如， CT图像主成分的张量为三阶张量，参考模块的大小为一阶和二阶，测试模块数目为三阶。
[0045] 将与X。相似性测度差异最小的（N-I)个模块[X η ...，Xn J映射到三阶张量 s 中，其中，模块大小为一阶和二阶，训练样本为三阶。
[0046] 在获得与参考模块X。相似的一定数量训练模块后，采用基于张量的主成
分分析法完成张量基底训练。定义％&为η阶张量基底展开矩阵： 中第in行第（i n+1_l) In i+in i列的元素 ，I i，12, I3表示展开矩阵中张量的维数，那么
根据三阶奇异值分解原理，X可分解为：
表示酉矩阵，f 表示核张 量。每一阶正交矩阵可表示为I
[0049] 步骤 S5，
[0050] 对包含噪声的张量的基底进行削减。
[0051] 例如，可以削减包含噪声的张量基底，保留包含CT图像主要信息的成分。
[0052] 将（N-I)个相似模块映射到张量基底空间中，那么不包含噪声的能量主要集中在 几个主要成分中，而噪声能量均匀的分布在所有成分中，本节节引入LMMSE完成张量基底 肖IJ减。定义η = 1，2是Yw的第i列，那么线性最小均方差可表示为：
[0056] 其中，
分别表不信号和噪声的方差，表不X(n)的 第i列。根据Daubechies二阶小波变换，噪声可表示为：
「00581 其中，HH表示高颇小波系数。伸用最太似然估计得到信号方差为：
[0060] 其中，M表示Yw的列数。图像中纹理单一的区域，
通常 小于零，那么wn(p) = 〇, yw中的噪声可被去除。
[0061] 步骤 S6，
[0062] 将测试模块映射回图像空间中，重建去除噪声后的CT图像。
[0063] 定义你? ei?Px/w每一列包含Wn个元素，那么收缩变换#可表示为：
[0064]
[0065] 其中，（·）表示向量内积。通过收缩变换f可将相似模块从张量基底空间映射回 图像空间，得到去除噪声的结果X，可表示为：
[0066]
[0067] 定义先为X。去除噪声的模块。对整幅CT图像；[e 7? 中的M1X M2个像素都采 用上述方法完成去噪，得到去噪后的图像。由于每个像素可包含在多个参考模块 而得到多组去除噪声后的值，本节通过计算其灰度平均值得到去噪后的图像，计算公式如 下：
[0069] 其中，E1表示笨的数目，xR表示X。的多组去噪后的值。
[0070] 当然，以上所述是本发明的优选实施方式。为方便说明起见，使用了步骤SI、S2等 序号，但是应该认识到的是，这些步骤本身还可以包括其他过程，这些步骤之间还可以有其 他步骤，这也在本发明的保护范围之内。应当指出，对于本技术领域的普通技术人员来说， 在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发 明的保护范围。
【主权项】
1. 一种基于主成分分析的CT图像去噪方法，其特征在于，包括： 步骤S1，构建一个可自适应调整大小的搜索窗口； 步骤S2,以CT图像中一个包含噪声的像素为中心，构建矩形参考模块； 步骤S3,在搜索窗口中，寻找与参考模块相似的测试模块，构成测试模块集； 步骤S4,测试模块集映射到CT图像主成分的张量基底空间中； 步骤S5,对包含噪声的张量基底进行削减； 步骤S6,将测试模块映射回图像空间中，重建去除噪声后的CT图像。2. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 Sl，搜索窗口为矩形窗口。3. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 Sl，搜索窗口的大小可以根据不同窗口的特征自动调整。4. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 52, 参考模块是以CT图像中一个包含噪声的像素为中心的矩形窗口，窗口的大小为固定 值。5. 根据权利要求4所述的基于主成分分析的CT图像去噪方法，其特征在于，窗口大小 为 15x15。6. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 53, 两个模块相似的评测标准为两个模块中相对应位置的像素的灰度值差异最小。7. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 54, 通过预设阈值，提取CT图像的主成分，构成张量基底空间。8. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 54, CT图像主成分的张量为三阶张量，参考模块的大小为一阶和二阶，测试模块数目为三 阶。9. 根据权利要求1所述的基于主成分分析的CT图像去噪方法，其特征在于，在步骤 55, 削减包含噪声的张量基底，保留包含CT图像主要信息的成分。
【专利摘要】本发明公开一种基于主成分分析的CT图像去噪方法。该方法包括：构建一个可自适应调整大小的搜索窗口；以CT图像中一个包含噪声的像素为中心，构建矩形参考模块；在搜索窗口中，寻找与参考模块相似的测试模块，构成测试模块集；测试模块集映射到CT图像主成分的张量基底空间中；对包含噪声的张量的基底进行削减；将测试模块映射回图像空间中，重建去除噪声后的CT图像。本发明能够准确的去除CT图像中的噪声，并且可以保留CT图像中的边缘和细节信息。
【IPC分类】G06T5/00
【公开号】CN105046662
【申请号】CN201510393426
【发明人】郑重, 吴文波, 杨文晖, 赖暖翔
【申请人】嘉恒医疗科技（上海）有限公司
【公开日】2015年11月11日
【申请日】2015年7月6日