基于DNST的医学CT图像去噪方法技术

基于DNST的医学CT图像去噪方法，包括如下步骤：步骤1)建立医学CT图像模型；步骤2)通过离散不可分离剪切波对CT图像进行分解；步骤3)对高频剪切波系数进行双变量收缩处理；步骤4)对低频剪切波系数进行旋转不变双边非局部均值滤波处理步骤4)对处理后的系数进行DNST逆变换。本发明专利技术通过实验分析与最新的去噪领域算法进行了对比，有效的应用在医学CT去噪领域；特别是针对CT低频采用旋转不变双边非局部均值滤波和对高频剪切波系数进行双变量收缩，能够更好的保护医学CT图像中软组织的细节信息。通过了大量的实验数据对比，提出了基于剪切波域双变量收缩和旋转不变双边非局部均值滤波结合的医学CT图像去噪方法，能够更好的有利于医师的分析诊断。

【技术实现步骤摘要】
基于DNST的医学CT图像去噪方法
本专利技术涉及于医学图像去噪领域，特别是涉及医学CT图像。设计一种基于剪切波域双变量收缩和旋转不变双边非局部均值滤波结合的医学CT图像去噪方法。
技术介绍
在人类发展史上，医学一直是人们极度重视和不断发展的学科。医学图像数字化的普及和发展极大的提高了医疗诊断的效率性和准确性。其中，计算机断层扫描(ComputedTomography，CT)在放射诊断方面具有很强的能力，它是第一种在不打开人体的情况下获得人体内部结构的轴向切片图，并且其具备组织结构密度分辨率高，对人体损害小，操作简便等优点广泛的应用在诊断诸如癌症，心血管疾病，传染病，阑尾炎，创伤和肌肉骨骼疾病等领域。医学CT的专利技术史造就了三位诺贝尔奖从侧面反映CT成像技术的重要性。CT是一种采用断层扫描的医学成像方法，能够很简单的获取物体的内部结构特征，如尺寸、形状、内部缺陷和密度，它是一种强大的非破坏性评估技术，广泛的用于生成人体的二维和三位横截面图像。CT成像运用物理技术，通过测量X射线在人体内的衰减系数为基础，充分利用X射线在人体内被吸收的特性。通过医学仪器中探测器将接收到射线变成电信号，这些电信号通过模拟数字转换得到射线衰减系数值。采用优秀的数学算法，根据得到的射线衰减系数值的二维分布矩阵，经过计算机求解出在人体某个剖面的图像画面的灰度值，从而得到人体某部分截面图。但由于CT定位、定性诊断的准确性仍然收到病变部位、大小、性质、病程长短等因素的影响，不能反映脏器的功能信息。图像也容易收到噪声干扰、被照物细节等因素的影响，从而观察到图像中有斑点和细粒，进一步影响到医生从病灶部位一定大小细节从背景中鉴别出来。图像中的噪声干扰必将影响图像的进一步分析，因此研究相关医学图像去噪算法将噪声从CT图像中滤除，可为观察图像提供好的视觉环境。本专利技术使用医学CT为研究对象。由于CT成像不免受到各种物理因素的影响，斑点噪声的存在严重影响了CT图像的质量，导致了医学图像质量较差。噪声的存在掩盖那些灰度差别很小的图像特征。对于临床医生而言，噪声对他们的准确诊断造成了很大的干扰，特别是对于经验不是很丰富的医生造成的影响更大。因此，从临床应用的角度出发，需要研究对CT医学图像去噪方法，为医生做出更准确的诊断提供技术支持，降低人工诊断的风险。综上所述，研究医学CT图像去噪方法具有广泛的应用。
技术实现思路
本专利技术要克服小波分析处理高维数据稀疏能力的不足和离散剪切波框架界效应差的不足，提供一种离散不可分离剪切波(DNST)医学CT图像去噪算法,用于解决医学CT图像去噪。现有技术中，小波变化能很好的用于图像去噪并且有效的抓住一维奇点，但不能反映直线和曲线的突变。脊波变换可以很好的捕捉线的奇异性，弥补小波的不足，但是仍然不能有效的捕捉曲线的奇异性。近些年来，通过离散剪切波算法对医学图像的处理，使得对医学图像去噪
有一定的突破。DNST相比较于离散可分离剪切波拥有更好的框架界和更好的定向选择，意味着DNST拥有更好的去噪效果。本专利技术上将DNST工具包用到医学CT噪声图像中，并且修改了传统的收缩算法使用双变量收缩算法，并且在低频部分结合CT图像中旋转不变双边非局部均值滤波。专利技术了具有速度快、去噪明显的设计一种基于剪切波域双变量收缩和旋转不变双边非局部均值滤波结合的医学CT图像去噪方法，最后通过仿真验证了方法的可行性与优化的效果。与现有技术相比，本专利技术的创造性与优点：提出了一种设计一种基于剪切波域双变量收缩和旋转不变双边非局部均值滤波结合的医学CT图像去噪方法，DNST变换克服了小波分析处理高维数据稀疏能力的不足，离散剪切波框架界的不足。通过双变量收缩算法比传统的硬阈值和软阈值更适合应用在CT去噪领域，并且在低频部分添加了适用与CT特征的旋转不变双边非局部均值滤波。并且此方法具有多分辨、多尺度、多方向性和时频局部性，将其应用于CT图像去噪能更好的保护图像边缘信息，给医师的诊断提供了方便。为使本专利技术的目的、技术方案和优点更加清晰，下面就对本专利技术的技术方案作进一步描述，分别分为以下四个步骤。基于DNST的医学CT图像去噪方法，包括以下步骤：步骤1)建立医学CT图像模型；为了解决CT噪声问题，我们不能凭借人的主观感受进行断定，通常噪声只能用概率统计的方法去认识。因此，最主要的是将CT图像中抽象的噪声进行原理具体化，建立一个符合基本特征的数学模型。X射线的能量以量子的单个能量块的方式进行传输，因此，这些无限量的X射线量子被X射线探测器检测到。由于统计错误，检测到的X射线量子可能因为再一次测量而不同。CT图像中的统计噪声由于在检测少量X射线量子点可能出现的波动而产生，因此统计噪声也可以称为量子噪声。从数学模型上来说这是一种高斯加性噪声，其数学模型如下：Inoisy(x,y)＝uorignal(x,y)+ηnoisy(x,y)(1)其中(x,y)表示图像中的像素点，Inosiy表示含有噪声CT图像即观察到的CT图像，uorignal表示无噪声CT图像，ηnoisy表示噪声，其分布符合高斯分布。步骤2)通过离散不可分离剪切波对CT图像进行分解；S1：输入一个二维CT图像信号f∈RX*Y，尺度参数J∈N，一个方向参数k∈NJ，以及选择方向滤波器DirectionFilter、低通滤波器QuadratureMirrorFilter。S2：转换CT图像的频率谱ffreq＝FFT(f)。S3：计算不同尺度不同方向下的子带i∈[0,nth]下的剪切波系数shearletCoeffs(i)∈RX*Y*nth，根据剪切波域的卷积理论可得:S4：通过计算得出离散不可分离剪切波系数shearletCoeffs(i)。其中第3步中nth代表了整个紧支撑DNST系统的冗余度，其计算如下:nth＝2*((2*2k[0]+1))+2*((2*2k[1]+1))+...+2*((2*2k[J]+1)(3)DNST能够得到更好的框架界。除此之外，由不可分离的剪切波发生器生成的剪切波的一个主要优点是扇形滤波器P在频域的每一个尺寸都提高了方向选择性。通过DNST分解，将医学CT图像在频域内分解成f1,f2,...,fnth-1张大小相等的高频CT图像和一张低频CT图像fnth。步骤3)对高频剪切波系数f1,f2,...,fnth-1进行双变量收缩处理；图像去噪最强大的技术之一就是小波域中常用的双变量阈值函数去噪方法。这里，利用这个想法被扩展到了剪切波域。通过利用非高斯双变量函数统计剪切波系数模型。首先定义了更精细尺度和粗糙尺度的关系，然后对剪切波系数进行有效的建模，每个系数取决于称为CS的较粗尺度中的相同空间位置处的系数，并且每个CS取决于称为FS的即时更精细尺度中的相同空间位置中的系数。因此，剪切波域中的每个CS有一个CS，每个CS有四个FS。该模型可以有效地捕捉框架系数与其CS之间的依赖关系。通常，对于给定的原始无噪声图像，如果它受到加性高斯白噪声的干扰，则降级后的图像在剪切波域表示如下所示：g＝f+ε(4)其中，g,f,ε分别表示观察到的剪切波系数，原始无噪声剪切波系数和噪声的系数。去噪的目的在于获得原始系数的估计使得估计值尽可能接近原始无噪声系数f。在所提出的方法中，最大后验概率本文档来自技高网...

【技术保护点】
1.基于DNST的医学CT图像去噪方法，包括以下步骤：步骤1)建立医学CT图像模型；对CT中的噪声来源分析，对抽象的噪声进行统计，建立医学CT模型，具体包括：X射线的能量以量子的单个能量块的方式进行传输，因此，这些无限量的X射线量子被X射线探测器检测到；由于统计错误，检测到的X射线量子可能因为再一次测量而不同；CT图像中的统计噪声由于在检测少量X射线量子点可能出现的波动而产生，因此统计噪声也可以称为量子噪声；从数学模型上来说这是一种高斯加性噪声，其数学模型如下：Inoisy(x,y)＝uorignal(x,y)+ηnoisy(x,y)  (1)其中(x,y)表示图像中的像素点，Inosiy表示含有噪声CT图像即观察到的CT图像，uorignal表示无噪声CT图像，ηnoisy表示噪声，其分布符合高斯分布；步骤2)通过离散不可分离剪切波对CT图像进行分解；并计算分解后子带的剪切波系数，具体包括：S1：输入一个二维CT图像信号f∈RX*Y，尺度参数J∈N，一个方向参数k∈NJ，以及选择方向滤波器Direction Filter、低通滤波器Quadrature Mirror Filter；S2：转换CT图像的频率谱ffreq＝FFT(f)；S3：计算不同尺度不同方向下的子带i∈[0,nth]下的剪切波系数shearletCoeffs(i)∈RX*Y*nth，根据剪切波域的卷积理论可得:...

【技术特征摘要】
1.基于DNST的医学CT图像去噪方法，包括以下步骤：步骤1)建立医学CT图像模型；对CT中的噪声来源分析，对抽象的噪声进行统计，建立医学CT模型，具体包括：X射线的能量以量子的单个能量块的方式进行传输，因此，这些无限量的X射线量子被X射线探测器检测到；由于统计错误，检测到的X射线量子可能因为再一次测量而不同；CT图像中的统计噪声由于在检测少量X射线量子点可能出现的波动而产生，因此统计噪声也可以称为量子噪声；从数学模型上来说这是一种高斯加性噪声，其数学模型如下：Inoisy(x,y)＝uorignal(x,y)+ηnoisy(x,y)(1)其中(x,y)表示图像中的像素点，Inosiy表示含有噪声CT图像即观察到的CT图像，uorignal表示无噪声CT图像，ηnoisy表示噪声，其分布符合高斯分布；步骤2)通过离散不可分离剪切波对CT图像进行分解；并计算分解后子带的剪切波系数，具体包括：S1：输入一个二维CT图像信号f∈RX*Y，尺度参数J∈N，一个方向参数k∈NJ，以及选择方向滤波器DirectionFilter、低通滤波器QuadratureMirrorFilter；S2：转换CT图像的频率谱ffreq＝FFT(f)；S3：计算不同尺度不同方向下的子带i∈[0,nth]下的剪切波系数shearletCoeffs(i)∈RX*Y*nth，根据剪切波域的卷积理论可得:S4：通过计算得出离散不可分离剪切波系数shearletCoeffs(i)；其中第3步中nth代表了整个紧支撑DNST系统的冗余度，其计算如下:nth＝2*((2*2k[0]+1))+2*((2*2k[1]+1))+...+2*((2*2k[J]+1)(3)DNST能够得到更好的框架界；除此之外，由不可分离的剪切波发生器生成的剪切波的一个主要优点是扇形滤波器P在频域的每一个尺寸都提高了方向选择性；通过DNST分解，将医学CT图像在频域内分解成f1,f2,...,fnth-1张大小相等的高频CT图像和一张低频CT图像fnth；步骤3)对高频剪切波系数f1,f2,...,fnth-1进行双变量收缩处理；通过剪切波系数模型之间系数的相互联系，通过双变量收缩函数进行处理，具体包括：利用非高斯双变量函数统计剪切波系数模型；首先定义了更精细尺度和粗糙尺度的关系，然后对剪切波系数进行有效的建模，每个系数取决于称为CS的较粗尺度中的相同空间位置处的系数，并且每个CS取决于称为FS的即时更精细尺度中的相同空间位置中的系数；因此，剪切波域中的每个CS有一个CS，每个CS有四个FS；该模型能有效地捕捉框架系数与其CS之间的依赖关系；通常，对于给定的原始无噪声图像，如果它受到加性高斯白噪声的干扰，则降级后的图像在剪切波域表示如下所示：g＝f+ε(4)其中，g,f,ε分别表示观察到的剪切波系数，原始无噪声剪...

【专利技术属性】
技术研发人员：程芸，周俊，田峥，赵恺伦，张聚，
申请(专利权)人：浙江医院，
类型：发明
国别省市：浙江,33