提出了改进的自回归模型(AR)在SPECT图像减少噪音。一个AR滤波器应用于预滤器投影图像和预滤器有序子集期望最大化(OSEM)重建图像(AR-OSEM-AR方法)。这种方法的性能与过滤后投影(FBP)之前巴特沃斯滤波(BW-FBP方法)和OSEM重建方法之后,巴特沃斯滤波(OSEM-BW方法)。一个数学汽缸幻影被用于这项研究。它是热的和冷的对象。执行的测试是使用三个模拟SPECT数据集。图像质量评估通过百分比对比分辨率(CR %)和半宽度(应用)的直线气缸的传播功能。BW-FBP方法显示最高的CR %值和AR-OSEM-AR方法给冷堆CR %值最低。在热分析栈,BW-FBP方法有更高的CR %值比OSEM-BW方法。冷BW-FBP方法表现出最低的应用值栈和AR-OSEM-AR热栈的方法。 In conclusion, the AR-OSEM-AR method is a feasible way to remove noise from SPECT images. It has good spatial resolution for hot objects.
无数的方法去除噪声从SPECT图像提出了
在二维AR模型,每个值图像的退化在其周围像素值,称为预测区域。AR模型可以被看作是一个低通滤波器,将图像划分为两个附加组件,一个可预测的图像和预测误差图像。一个基于“增大化现实”技术的过程
在一个典型的scintigraphic形象,有大型本地空间图像的统计数字的变化。因此,相同的模型不能应用于整个图像,但是模型必须适应变化。在这个自适应方法,图像区域划分成更小的块和AR模型然后安装到每一块分别通过使用MATLAB子例程。最近,还引入了一个基于块去噪方法对三维超声图像(
自回归去噪过程的流程图。汪:自回归过滤;
数据模拟使用3 d-mac幻影(
Transaxial片使用过滤反射影重建方法(FBP) [
获得一个公平的比较的方法,相同数量的过滤应用在每个方法。这样做是通过一个圆形的(ROI) 150毫米直径的均匀的幻影和计算变异系数比例在ROI (x %),也就是说,标准差的比值平均值乘以100。这种表示确保过滤每个方法之间是平等的。
百分比对比分辨率(CR %)值计算活动在每个气缸和统一的活动。CR %可以由以下公式来表达(
空间分辨率是估计的半宽度(应用)的直线气缸的传播功能。一,二,四,six-pixel-thick概要文件被吸引在10 - 20 - 40 - 60-mm-wide圆柱体,分别。FWMH值计算使用爱马仕的质量控制软件(版本2.0)。
骨骼SPECT进行静脉注射后三个小时925兆贝可
数据分析使用WinSTAT Excel(版本2007.1;r .惠誉软件,Staufen,德国)。成对比较进行非参数Wilcoxon rank-sum测试。之间的比较是AR-OSEM-AR和BW-FBP方法,AR-OSEM-AR和OSEM-BW方法,BW-FBP和OSEM-BW方法。冷栈和数据池热堆分别进行分析。对于每个气缸,搭配一个变量的值是计算之间的区别。的值差异分类排序。最后,消极的平均等级差异比较与积极的差异。Wilcoxon rank-sum测试决定在多大程度上意味着等级的差异是显著的。一个
图像的均方误差改进不同的AR模型用于总结时误差项形象而非原始图像。预测地区的四个正交邻国的块大小
的块大小变化的影响总结与预测误差项图像区域3×3像素。可预测的预测区域形象的四个正交邻居5×5像素的块大小是使用。
| 总数量 | 块大小 | 均方误差 |
|---|---|---|
| 28705年 | 5×5 | 0.87 |
| 28705年 | 6×6 | 0.86 |
| 28705年 | 7×7 | 0.86 |
|
|
||
| 54469年 | 5×5 | 2.12 |
| 54469年 | 6×6 | 2.10 |
| 54469年 | 7×7 | 2.14 |
|
|
||
| 108938年 | 5×5 | 6.61 |
| 108938年 | 6×6 | 6.56 |
| 108938年 | 7×7 | 7.04 |
Transaxial Zubal幽灵的切片。(一)Poisson-noise-corrupted transaxial片。(b)迭代过滤可预测的图像。(c)过滤总结误差项的形象。(d)最终的图像。图像单独扩展自己的最大值。逆线性灰度用于比较原始的幻影。总数量水平是108791年Poisson-noise-corrupted形象,102237年迭代过滤可预测的图像,在过滤后的总结和6205年误差项的形象。
巴特沃斯滤波选择这样的方法有相同数量的统计波动统一幻影的一部分,由x %值(表确认
比例系数变化对不同重建技术。
| 方法 | 计数水平 | RelAct | CF(周期/厘米) | x % |
|---|---|---|---|---|
| AR-OSEM-AR | 50000年 | 0 | - - - - - - | 6.35 |
| BW-FBP | 50000年 | 2 | 0.83 | 6.37 |
| OSEM-BW | 50000年 | 4 | 0.84 | 6.36 |
|
|
||||
| AR-OSEM-AR | 100000年 | 0 | - - - - - - | 4.60 |
| BW-FBP | 100000年 | 2 | 0.80 | 4.56 |
| OSEM-BW | 100000年 | 4 | 0.84 | 4.65 |
|
|
||||
| AR-OSEM-AR | 150000年 | 0 | - - - - - - | 4.43 |
| BW-FBP | 150000年 | 2 | 0.89 | 4.50 |
| OSEM-BW | 150000年 | 4 | 0.86 | 4.43 |
AR-OSEM-AR:自回归滤波前后有序子集期望最大化算法;BW-FBP:巴特沃斯前置过滤,过滤后的投影;巴特沃斯OSEM-BW:有序子集期望最大化算法和预滤器;RelAct:活动相对于背景的活动1;CF:截止频率。过滤器的顺序是2;——:不是可定义。
百分比对比分辨率值不同的方法。
| 方法 | RelAct | 10Ø | 20Ø | 40Ø | 60Ø |
|---|---|---|---|---|---|
| 一个 | |||||
| AR-OSEM-AR | 0 | 18.4 | 56.3 | 74.6 | 84.2 |
| BW-FBP | 0 | 28.7 | 67.2 | 83.1 | 91.6 |
| OSEM-BW | 0 | 24.9 | 59.5 | 75.7 | 85.6 |
|
|
|||||
| AR-OSEM-AR | 2 | 9.0 | 86.1 | 84.3 | 103.5 |
| BW-FBP | 2 | 16.7 | 75.5 | 86.6 | 97.2 |
| OSEM-BW | 2 | 9.9 | 65.3 | 76.9 | 98.3 |
|
|
|||||
| AR-OSEM-AR | 4 | 136.5 | 253.0 | 270.4 | 296.5 |
| BW-FBP | 4 | 84.9 | 229.2 | 259.2 | 289.1 |
| OSEM-BW | 4 | 70.2 | 204.1 | 247.1 | 300.0 |
|
|
|||||
| B | |||||
| AR-OSEM-AR | 0 | 21.5 | 58.0 | 76.1 | 84.8 |
| BW-FBP | 0 | 20.4 | 61.6 | 83.2 | 91.8 |
| OSEM-BW | 0 | 32.9 | 61.9 | 78.0 | 86.6 |
|
|
|||||
| AR-OSEM-AR | 2 | 32.8 | 81.0 | 81.9 | 95.7 |
| BW-FBP | 2 | 35.1 | 80.2 | 87.4 | 102.7 |
| OSEM-BW | 2 | 29.7 | 75.4 | 78.8 | 98.3 |
|
|
|||||
| AR-OSEM-AR | 4 | 132.8 | 230.2 | 269.8 | 295.7 |
| BW-FBP | 4 | 94.6 | 236.9 | 270.3 | 316.2 |
| OSEM-BW | 4 | 85.6 | 224.6 | 255.1 | 308.5 |
|
|
|||||
| C | |||||
| AR-OSEM-AR | 0 | 24.0 | 55.9 | 75.0 | 85.5 |
| BW-FBP | 0 | 32.9 | 68.7 | 83.0 | 93.0 |
| OSEM-BW | 0 | 31.2 | 60.7 | 76.8 | 87.8 |
|
|
|||||
| AR-OSEM-AR | 2 | 41.1 | 76.0 | 86.9 | 94.9 |
| BW-FBP | 2 | 38.8 | 82.4 | 93.9 | 103.6 |
| OSEM-BW | 2 | 34.3 | 83.1 | 85.5 | 100.6 |
|
|
|||||
| AR-OSEM-AR | 4 | 119.4 | 237.1 | 266.3 | 297.7 |
| BW-FBP | 4 | 112.7 | 248.5 | 274.5 | 318.2 |
| OSEM-BW | 4 | 108.1 | 260.5 | 261.0 | 316.9 |
答:计数水平低;B:中间数水平;C:高数水平;AR-OSEM-AR:自回归滤波前后有序子集期望最大化算法;BW-FBP:巴特沃斯前置过滤,过滤后的投影;巴特沃斯OSEM-BW:有序子集期望最大化算法和预滤器;RelAct:活动相对于背景的活动1;Ø:直径。
在空间分辨率的分析,没有例外,BW-FBP和OSEM-BW方法表现出冷栈应用值低于AR-OSEM-AR方法(
半宽度值不同的方法。
| 方法 | RelAct | 10Ø | 20Ø | 40Ø | 60Ø |
|---|---|---|---|---|---|
| 一个 | |||||
| AR-OSEM-AR | 0 | 19.7 | 25.0 | 45.8 | 64.8 |
| BW-FBP | 0 | 17.9 | 23.2 | 43.5 | 63.8 |
| OSEM-BW | 0 | 19.4 | 23.2 | 43.6 | 62.8 |
|
|
|||||
| AR-OSEM-AR | 2 |
|
21.5 | 35.3 | 59.1 |
| BW-FBP | 2 |
|
23.2 | 36.6 | 59.6 |
| OSEM-BW | 2 |
|
22.2 | 35.8 | 59.1 |
|
|
|||||
| AR-OSEM-AR | 4 | 13.4 | 18.4 | 37.8 | 57.8 |
| BW-FBP | 4 | 15.8 | 20.4 | 37.9 | 57.9 |
| OSEM-BW | 4 | 15.2 | 19.9 | 37.7 | 57.9 |
|
|
|||||
| B | |||||
| AR-OSEM-AR | 0 | 23.7 | 23.7 | 43.0 | 62.3 |
| BW-FBP | 0 | 18.6 | 22.4 | 41.2 | 61.4 |
| OSEM-BW | 0 | 18.8 | 23.2 | 41.9 | 61.7 |
|
|
|||||
| AR-OSEM-AR | 2 | 19.4 | 19.8 | 38.9 | 59.0 |
| BW-FBP | 2 | 19.0 | 21.0 | 38.5 | 58.9 |
| OSEM-BW | 2 | 17.7 | 19.9 | 38.5 | 59.0 |
|
|
|||||
| AR-OSEM-AR | 4 | 13.0 | 19.7 | 37.8 | 58.2 |
| BW-FBP | 4 | 15.4 | 21.0 | 37.7 | 57.9 |
| OSEM-BW | 4 | 14.2 | 20.2 | 37.9 | 58.1 |
|
|
|||||
| C | |||||
| AR-OSEM-AR | 0 | 20.2 | 23.9 | 40.9 | 63.2 |
| BW-FBP | 0 | 18.5 | 21.1 | 40.6 | 61.7 |
| OSEM-BW | 0 | 16.8 | 23.2 | 40.3 | 62.3 |
|
|
|||||
| AR-OSEM-AR | 2 | 15.8 | 19.5 | 39.0 | 58.6 |
| BW-FBP | 2 | 16.9 | 20.8 | 38.6 | 58.7 |
| OSEM-BW | 2 | 16.5 | 20.3 | 39.1 | 58.6 |
|
|
|||||
| AR-OSEM-AR | 4 | 12.5 | 19.0 | 37.6 | 58.2 |
| BW-FBP | 4 | 14.0 | 19.6 | 36.8 | 58.0 |
| OSEM-BW | 4 | 13.8 | 19.3 | 37.4 | 57.9 |
答:计数水平低;B:中间数水平;C:高数水平;AR-OSEM-AR:自回归滤波前后有序子集期望最大化算法;BW-FBP:巴特沃斯前置过滤,过滤后的投影;巴特沃斯OSEM-BW:有序子集期望最大化算法和预滤器;RelAct:活动相对于背景的活动1;Ø:直径。
此外,OSEM-BW值较低的半最大值宽度比BW-FBP方法(
视觉上,三种方法产生的图像之间的差异是很小的(数字
Transaxial片幻影的圆柱体直径20毫米。(一)图像从无噪声的投影重建图像使用命令子集期望最大化重建。(b)自回归滤波前后有序子集期望最大化重建。(c)之前巴特沃斯滤波过滤反射影重建方法。期望最大化(d)有序子集重建巴特沃斯滤波紧随其后。中间数水平。图像单独扩展自己的最大值。
重新格式化冠状切片的幻影。(一)图像从无噪声的投影重建图像使用命令子集期望最大化重建。(b)自回归滤波前后有序子集期望最大化重建。(c)过滤反射影重建方法之前巴特沃斯。期望最大化(d)有序子集重建巴特沃斯滤波紧随其后。栈最高的活动。中间数水平。图像单独扩展自己的最大值。
通过幻影One-pixel-thick资料绘制。(a)概要文件的一个圆柱体直径20毫米。(b)档案的圆柱体,直径40毫米。理想OSEM:图像重建使用命令子集的期望最大化从无噪声的投影图像重建;AR-OSEM-AR:自回归滤波前后有序子集期望最大化重建;BW-FBP:过滤反射影重建方法之前巴特沃斯滤波;期望最大化OSEM-BW:命令子集重建巴特沃斯滤波紧随其后。中间数水平。概要文件是新,以便他们有相同数量的计数从无噪声的投影图像的轮廓重建的图像使用命令子集期望最大化重建。
Transaxial块骨骼SPECT。(一)自回归期望最大化过滤前后有序子集重建。(b)之前巴特沃斯滤波过滤反射影重建方法。期望最大化(c)有序子集重建巴特沃斯滤波紧随其后。
本文提出了一种改进的二维自适应AR滤波器,介绍了一个三维的自适应AR模型减少噪音SPECT图像。我们证明scintigraphic图像的质量可以提高时相同的基于“增大化现实”技术的过程并不应用于原始图像,总结误差项的形象。之前我们已经表明,如果预测地区的四个正交邻国的预测像素的块大小
SPECT过滤的目的是抑制统计噪声,同时保持对比度和空间分辨率(
BW-FBP方法产生更好的性能比其他两个栈的方法分析冷。之间的良好的性能和冷特性(之前已经注意到
FBP方法由投影数据的筛选和过滤数据的投影(
出口押汇的缺点是,它可以产生径向条纹文物因为过滤噪声投影资料不背投影相互抵消。在目前的研究中,这种现象被认为在临床数据。迭代重建算法也提供了一些其他的优势出口押汇。他们允许使用的几个重要的修正,如散射、衰减,和准直器响应修正,可以包含在图像重建过程。的解剖信息来源于磁共振成像或电脑断层摄影术也是可能的
巴特沃斯滤波器是由两个参数:截止频率和秩序(
在我们看来,AR-OSEM-AR方法的一个优点是它的简单性,但缺乏用户控制变量也可以被视为一个限制。有时,可调参数是必要的。没有特定的过滤器可以成为任何器官系统的最佳滤波器。然而,应该执行过滤本地空间域,不是全球在频域,因为正确的分辨率之间的权衡和在不同的点在图像平滑会有所不同。
因为AR-OSEM-AR方法只是略微比OSEM-BW方法,一个额外的研究需要确定图像质量会更好如果AR方法应用于迭代之间的中间结果。其次,AR-OSEM-AR方法尚未经过测试与正电子发射断层扫描(PET)数据,但该方法也应该适用于宠物数据。信噪比是相当高的比SPECT的宠物。因此,我们的模型可能会提供一个适合的宠物数据。
SPECT的AR-OSEM-AR方法是一种可行的去噪方法。具有良好的空间分辨率的热特性和简单的使用。它没有任何可调参数。
作者宣称没有利益冲突有关的出版。