降噪方法技术

在降低计算机断层摄影（CT）等的测定信号的噪声的方法中，抑制滤波器效果对于滤波器参数的依赖性。降噪方法包含：a）根据规定的噪声模型，决定测定信号各自的噪声方差的步骤；b）根据噪声方差的P次方，产生针对各测定信号的滤波器的离散核的步骤，该针对各测定信号的滤波器的离散核表示滤波器对于各测定信号的频率响应；以及c）应用分别对应的上述离散核，对测定信号进行滤波的步骤。

【技术实现步骤摘要】
降噪方法
本专利技术的实施方式涉及一种降噪方法。
技术介绍
受到弱光引起的噪声及条纹(streak)的影响，X射线CT图像的质量有时候会严重变差。这个问题虽然可以通过增加X射线的剂量得到缓解，但从患者安全性方面考虑，这种方法在临床上是不被允许的。为了在安全辐射剂量等级的前提下达到在诊断上有用的图像质量，在以往寻求用于大幅度地降低噪声及条纹的优选的解决策略，近数十年来进行了各种各样的尝试。受到关于低辐射剂量的意识提升的影响，所提出的上述研究在最近重要性显著增加，备受关注。为了改良结果，在以往的尝试中，利用自适应滤波器来代替固定滤波器。作为固定滤波器的若干示例有三角滤波器及双向滤波器。作为自适应滤波器的若干示例有自适应型高斯滤波器、自适应型裁剪均值滤波器。在以往技术的自适应滤波器中，其滤波器参数需要适当地进行调整。依照经验，对每个图像分别实现此调整的最佳化。在其他的例示的滤波器中，高斯核的方差和数据的噪声方差相同。同样地，由于是对图像分别进行调整，因此利用上述以往技术的自适应滤波器执行的噪声和/或条纹的降低并不充分。依然期望有一种有效的滤波器，该滤波器虽然使噪声及条纹的降低达到最大程度但其滤波器参数却不用复杂或特定的方法来进行调整。
技术实现思路
目的是在降低计算机断层摄影(CT)等的测定信号的噪声的方法中，抑制滤波器效果对于滤波器参数的依赖性。本实施方式涉及的降噪方法包含：a)根据规定的噪声模型，决定测定信号各自的噪声方差的步骤；b)根据上述噪声方差的P次方，产生针对各上述测定信号的滤波器的离散核的步骤，该针对各上述测定信号的滤波器的离散核表示针对各上述测定信号的滤波器的频率响应；以及c)应用分别对应的上述离散核，对上述测定信号进行滤波的步骤。附图说明图1是表示基于本实施方式的多切片X射线CT装置或扫描器的结构的图。图2是表示利用噪声模型预测的方差(y轴)和实测方差(x轴)的关系的图表。图3是表示在图2的计数范围、及利用噪声模型预测的方差(y轴)和实测方差(x轴)的关系变得格外非线性的计数较少的延伸范围内的上述关系的图表。图4是表示使用规定高斯滤波器的与维N、VarScale及VarPower相关的高斯滤波器的有效滤波器大小和计数的关系的图表。图5是表示使用基于本实施方式的具有第1组参数值的规定滤波器的噪声和/或条纹减少处理部的一实施方式中的滤波后的平均计数和对数转换后方差的关联的图表。图6是表示使用基于本实施方式的具有第2组参数值的规定滤波器的噪声和/或条纹减少处理部的一实施方式中的滤波后的平均计数和对数转换后方差的关联的图表。图7是表示使用基于本实施方式的具有第3组参数值的规定滤波器的噪声和/或条纹减少处理部的一实施方式中的滤波后的平均计数和对数转换后方差的关联的图表。图8是表示与使利用本实施方式重建CT图像前测定的数据的噪声和/或条纹大幅度地降低的一个例示性处理相关的步骤的流程图。图9A是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的肩部图像。图9B是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的肩部图像。图9C是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的肩部图像。图10A是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的其他肩部图像。图10B是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的其他肩部图像。图10C是作为基于本实施方式的自适应型高斯滤波器的应用结果而示出在临床上有意义的改善的其他肩部图像。图11A是基于图10A、10B及10C的同一数据而以肺窗(lungwindow)表示的缩放后的重建的图像。图11B是基于图10A、10B、及10C的同一数据而以肺窗表示的缩放后的重建的图像。图11C是基于图10A、10B、及10C的同一数据而以肺窗表示的缩放后的重建的图像。具体实施方式以下，一边参照附图，一边说明本实施方式。在此，参照附图，相同参照编号表示所有附图中对应的构造，尤其参照图1，图中表示包含台架100和其他处理部或单元的、基于本实施方式的多切片X射线CT装置或扫描器的一实施方式。如前视图所示，台架100还包含X射线管101、环状框架102、多列或二维排列型X射线检测器103。X射线管101及X射线检测器103是在直径方向上隔着被摄体(subject)S而安装于环状框架102，环状框架102绕着轴RA旋转。被摄体S沿着轴RA在图示的页面内外移动期间，旋转单元107使框架102以0.4秒/旋转等高速旋转。多切片X射线CT装置还包含电流调整器113和高电压产生装置109，该高电压产生装置109对X射线管101施加管电压，使X射线管101生成X射线。在一实施方式中，高电压产生装置109被安装于框架102。X射线向被摄体S辐射，用圆表示被摄体S的截面积。X射线检测器103为了检测透过被摄体S的辐射出的X射线，隔着被摄体S而位于X射线管101的相反侧。另外，参照图1，X射线CT装置或扫描器还包含用来处理来自X射线检测器103的被检测出的信号的其他处理部。数据收集电路或数据收集系统(DAS)104针对各通道将X射线检测器103输出的信号转换成电压信号，对其进行放大后再进一步转换成数字信号。X射线检测器103及DAS104构成为对每1旋转的规定的所有投影数(TPPR：totalnumberofprojectionsperrotation)进行操作。上述数据通过非接触式数据发送器105，被发送给台架100外部的控制中心内所容纳的前处理处理部106。前处理处理部106执行对原始数据的灵敏度校正等特定校正。接着，存储处理部112在即将重建处理前的阶段，存储还被称为投影数据的作为结果而生成的数据。存储处理部112和重建处理部114、显示处理部116、输入处理部115、及扫描计划帮助装置200一起，经由数据/控制总线而连接于系统控制器110。扫描计划帮助装置200包含为制定扫描计划而帮助摄像技师的功能。本实施方式的一实施方式还包含各种软件模块和硬件构成要素，使用规定滤波器使计算机断层摄影图像的条纹和/或噪声大幅度地降低。根据本实施方式的一方式，CT装置的降噪处理部117有利地执行噪声和/或条纹的减少。在一实施方式中，降噪处理部117经由数据/控制总线，动作性地连接于其他软件模块、和/或存储处理部112、重建处理部114、显示处理部116、及输入处理部115等系统构成要素。关于这一点，在本实施方式的其他实施方式中，未必由降噪处理部117单独地执行降噪功能和/或与此关联的任务。而且，降噪处理部117在本实施方式的代替实施方式中，通过任意选择而变成重建处理部114等其他处理部的一部分。一般来说，CT内的投影数据可以在规定的对数转换处理之后利用。在该对数转换处理中，将由于扫描对象而衰减的被测定后的X射线强度信号转换成线积分数据。然后，根据该线积分数据，利用已知的数学反演(mathematicalinversion)方法来重建CT图像。在基于本实施方式的噪声/条纹减少系统的例示性一实施方式中，降噪处理部117转换投影数据而使其返回至原始的X射线强度数据或光子计数测定值本文档来自技高网...

【技术保护点】
一种降噪方法，包含：a）根据规定的噪声模型，决定测定信号各自的噪声方差的步骤；b）根据上述噪声方差的P次方，产生针对各上述测定信号的滤波器的离散核的步骤，该针对各上述测定信号的滤波器的离散核表示滤波器针对各上述测定信号的频率响应；以及c）应用分别对应的上述离散核，对上述测定信号进行滤波的步骤。

【技术特征摘要】
2011.06.06 US 13/154,1241.一种降噪方法，包含：a)根据规定的噪声模型，决定针对测定信号各自的对数转换处理后的噪声方差的步骤；b)根据上述噪声方差的P次方，产生针对各上述测定信号的滤波器的离散核的步骤，该针对各上述测定信号的滤波器的离散核表示滤波器针对各上述测定信号的频率响应；以及c)应用分别对应的上述离散核，对上述测定信号进行滤波的步骤。2.根据权利要求1所述的降噪方法，其中上述规定的噪声模型是通过VarAL＝(Ve+WI)/I2定义的，上述I是上述测定信号，上述Ve是已知的电子噪声的值，上述W是信号增益。3.根据权利要求2所述的降噪方法，其中针对各上述测定信号的上述滤波器的上述离散核是通过对上述噪声方差的P次方乘以常数K而生成的。4.根据权利要求3所述的降噪方法，其中上述P是根据包含线性扩散式的尺度空间的导出而被决定的。5.根据权利要求3所述的降噪方法，其中上述P具有0.2至1的范围内的值，上述常数K具有1至10000的范围内的值。6.根据权利要求5所述的降噪方法，其中相对于上述测定信号，上述P为0.5，上述K为1。7.根据权利要求1所述的降噪方法，其中针对用于产生上述测定信号的检测器内的各检测器元件(i)的上述滤波器的上述离散核是由定义的，上述Δxi是1维图像中的第i个像素和规定的基准像素i0之间的距离，上述是上述基准像素i0中的上述滤波器的上述频率响应，上述VR是上述滤波器的参数。8.根据权利要求7所述的降噪方法，其中上述参数VR等于由VarAL＝(Ve+WI)/I2决定的对数转换后的方差，上述I是上述测定信号，上述Ve是已知的电子噪声的值。9.根据权利要求1所述的降噪方法，其中针对用于产生上述测定信号的检测器内的各检测器元件(i，j)的上述滤波器的上述离散核是由定义的，上述Δxi，j是2维图像中的(i，j)像素和规定的基准像素(i0，j0)之间的距离，其中，i和j是表示像素中的坐标，上述是上述基准像素(i0，j0)中的上述滤波器的频率响应，上述VR是上述滤波器的参数。10.根据权利要求9所述的降噪方法，其中上述参数VR等于由VarAL＝(Ve+WI)/I2决定的对数转换后的方差，上述I是上述测定信号，上述Ve是已知的电子噪声的值。11.根据权利要求1所述的降噪方法，其中...

【专利技术属性】
技术研发人员：杨智，A·察莫耶廷，邹宇，M·D·西尔弗，
申请(专利权)人：株式会社东芝，东芝医疗系统株式会社，
类型：发明
国别省市：