CN103180875A - 低剂量ct降噪 - Google Patents
低剂量ct降噪 Download PDFInfo
- Publication number
- CN103180875A CN103180875A CN2011800512979A CN201180051297A CN103180875A CN 103180875 A CN103180875 A CN 103180875A CN 2011800512979 A CN2011800512979 A CN 2011800512979A CN 201180051297 A CN201180051297 A CN 201180051297A CN 103180875 A CN103180875 A CN 103180875A
- Authority
- CN
- China
- Prior art keywords
- projection
- detection
- noise reduction
- photons
- estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000005855 radiation Effects 0.000 claims abstract description 27
- 230000009467 reduction Effects 0.000 claims description 63
- 238000001514 detection method Methods 0.000 claims description 46
- 238000000034 method Methods 0.000 claims description 18
- 238000009499 grossing Methods 0.000 claims description 13
- 238000007689 inspection Methods 0.000 claims description 13
- 238000003384 imaging method Methods 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000007476 Maximum Likelihood Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种系统包括围绕检查区域旋转并发射横穿所述检查区域的辐射的源、探测横穿所述检查区域的辐射并生成指示所探测到的辐射的投影数据的辐射敏感探测器阵列以及对所述投影数据降噪的投影数据降噪器,其中,所述降噪器基于所述投影的探测光子数量对所述投影降噪。
Description
技术领域
下文大体上涉及数据的降噪,其尤其适用于计算机断层摄影(CT),但是也同样适用于其他成像模态,例如,混合PET/CT系统、数字X射线系统、和/或其他成像模态。
背景技术
多切片计算机断层摄影(CT)扫描器包括安装在可旋转扫描架上的X射线管,所述可旋转扫描架绕纵轴或z轴围绕检查区域旋转。X射线管发射横穿检查区域以及其内的受试者或对象的辐射。二维探测器阵列限定了跨越检查区域与X射线管相对的具有一定角度的弧。所述探测器阵列包括多行探测器,它们相对于彼此对准并且沿z轴延伸。所述探测器探测横穿检查区域以及其内的受试者或对象的辐射,并生成指示其的投影数据。重建器对所述投影数据进行处理,并重建出指示其的三维(3D)体积图像数据。对所述体积图像数据进行处理,以生成所述检查区域的一幅或多幅图像,包括设置于所述检查区域中的对象和受试者的部分。
令人遗憾的是,CT扫描器发射电离辐射,因而使患者暴露于电离辐射之下,其可能提高患癌风险。一般而言,沉积在患者内的辐射剂量取决于多种因素,包括但不限于管电流(mAs)、管电压(kVp)、螺距/曝光时间(对于螺旋扫描而言)、切片厚度和间隔(对于轴向扫描而言)、研究当中的扫描次数以及患者体格(例如,较胖或较瘦)。可以通过降低管电流、管电压和/或扫描次数和/或提高螺距、切片厚度和/或切片间隔而降低所沉积的剂量。然而,图像噪声与辐射剂量成反比,因而降低辐射剂量不仅降低了患者体内沉积的剂量,还提高了所采集的数据内的图像噪声,在重建过程中,所述噪声将被传播至图像,从而降低了图像质量(即噪声更高的图像),由此可能降低所述过程的诊断价值。
基于图像的降噪算法已经得到了应用。然而,所述算法在处理“有条纹的”图像方面存在困难,在这样的图像中,噪声在图像数据的相邻体素之间存在很强的相关性。通常,投影测量当中的非常低水平的光子通量在重建图像当中生成条纹。而且,在所探测到的光子的平均数量非常低(例如,<10)时,所述算法操作将引入偏差,所述偏差将会作为偏移的平均CT值显现在所述图像内。迭代式重建,例如,基于最大似然(ML)的重建具有在这样的情况下获得改善的图像的潜力。然而,基于ML的重建计算成本极高,因而当前阻碍了其在例程实践中的应用。
发明内容
本申请的各个方面解决了上述和其他问题。
根据一个方面,一种系统包括围绕检查区域旋转并发射横穿所述检查区域的辐射的源、探测横穿所述检查区域的辐射并生成指示所探测到的辐射的投影数据的辐射敏感探测器阵列以及对所述投影数据降噪的投影数据降噪器,其中,所述降噪器基于所述投影的探测光子数量对所述投影降噪。
根据另一方面,一种方法包括获得成像系统生成的投影数据,针对所述投影数据的二维投影确定估计的探测光子数量并生成指示其的信号;以及基于所述信号对所述投影降噪。降噪程度与所述投影的估计的探测光子数量成比例,并且将在两个不同的程度上对至少两个投影降噪。
根据另一方面,一种方法包括对具有这样噪声的图像降噪:所述噪声在图像的相邻体素之间存在强相关性,其中,所述降噪不会向经降噪的图像中引入条纹伪像。
附图说明
本发明可以采取各种部件和部件布置的形式,以及各种步骤和步骤安排的形式。附图的作用在于对优选实施例进行图示,不应认为其对本发明构成限制。
图1结合投影数据降噪器示出了一种示范性成像系统。
图2示出了一种示范性投影数据降噪器。
图3示出了一种用于采用投影数据降噪器降低投影数据的噪声的示范性方法。
具体实施方式
图1示出了一种成像系统100,例如,其可以是计算机断层摄影(CT)扫描器。成像系统100包括基本上固定的扫描架102和旋转扫描架104。旋转扫描架104以可旋转的方式受到固定扫描架102的支持,并且绕纵轴或z轴108围绕检查区域106旋转。
诸如X射线管的辐射源110以可旋转的方式受到旋转扫描架104的支持。辐射源110与旋转扫描架104一起旋转,并发射横穿检查区域106的辐射。源准直器包括准直构件,其对辐射进行准直处理,从而使其形成大体为锥形、楔形、扇形或者其他形状的辐射束。
二维辐射敏感探测器阵列112界定了跨越检查区域106与辐射源110相对的呈一定角度的弧。探测器阵列112包括沿z轴108的方向延伸的多行探测器。探测器阵列112探测横穿检查区域106的辐射并生成指示其的投影数据。
投影数据降噪器114对投影数据降噪。就一个例子而言,在降噪器114采用的算法中,对应于较低光子数量的投影将比对应于较高光子数量的投影受到更为强烈的降噪,下文将对此给出更为详细的说明。这样的算法将能够减少重建图像中的由较低光子数量(例如,由于患者的身材、低剂量扫描等造成的)导致的条纹和/或偏差,同时保持所述数据中的强梯度(即,边缘),进而保持图像分辨率。
重建器116重建经降噪的投影数据,并生成指示其的三维(3D)体积图像数据。重建器116可以采用常规的3D滤波反投影重建、锥形射束算法、迭代算法和/或其他算法。
诸如卧榻的患者支撑118支持检查区域106内的诸如患者的对象或受试者。
通用计算系统或计算机起着操作员控制台120的作用。控制台120的处理器执行控制台126上的计算机可读指令,其允许操作员控制系统100的操作,例如,选择全剂量或低剂量扫描协议,激活投影数据降噪,启动扫描等。
在所示出的实施例中,将投影数据降噪器114示为单独的部件。在另一实施例中,投影数据降噪器114属于控制台120和/或其他计算装置的部分。
图2示出了投影数据降噪器114的例子。所示出的投影数据降噪器114包括对投影数据解记录(unlog)的数据解记录器200,其将衰减线积分转换为探测到的光子。
光子估计器202估计每一投影的探测光子数量,并生成指示其的信号。光子估计器202可以采用各种方案估计光子数量。例如,可以将探测到的光子数量估计为针对每一投影探测到的光子的平均数量。此外,可以(例如)采用移动平均值对这一均值进行平滑处理,由此可以促进缓解针对探测到的非常小的光子数量的泊松随机变量中的大的偏差。这里也可以设想其他技术,包括但不限于通常已知的用于估计探测光子数量的技术。
降噪器204基于估计的噪声对探测到的光子信号降噪。降噪器204能够为仅基于二维投影或者基于二维投影和一个或多个相邻二维投影对二维投影的探测光子进行降噪。如上所述,降噪器204可以采用这样的算法,即,对应于较低光子数量的投影将受到比对应于较高光子数量更为强烈的降噪,这将在保持边缘的同时缓解条纹。
数据记录器206记录经过降噪的探测光子,从而将经降噪的探测光子重新转换成衰减线积分,其能够通过重建器116得以重建。
评估器208可以用来基于预定的光子数量阈值确定是否应当对给定的投影降噪。在这种情况下,如果针对投影的估计光子数量表明该投影有足够数量的光子,那么不对该投影降噪。否则对该投影降噪。相对于对每一投影都降噪的做法而言,仅对被认定为没有足够的光子的投影降噪可以提高处理速度。在另一实施例中,可以省略评估器208。
下文提供了适当的降噪算法的非限制性例子,其以总变差最小化算法为基础,将每一2D投影作为图像处理。也可以采用具有适当的重新设计的类似的3D方法。
可以如方程1所示表示计算机断层摄影(CT)投影数据:
方程1:
其中,I表示测得的光子,I0表示输入光子,μx表示衰减函数,l表示采集线。
假设给定了光子测量结果的噪声特性,那么I表示平均光子测量结果,单个的实际测量结果f则是以I为均值,概率为的泊松随机变量的实现。
可以通过对代价函数求最小值而由f估计出平均光子测量结果I,其中,将所述代价函数表示为总变差项与加权最小二乘方项的和,如方程2所示:
方程2:
其中,v表示总统计加权,其将在噪声小时优选原始投影,在噪声大时优选具有小的总变差的投影。
一般而言,在估计的探测光子数量较低时,总变差项占据主导,而在估计的探测光子数量较高时,加权最小二乘方项则占据主导。
在方程2中,β是表示调节参数的任选变量,该参数控制进行平滑处理的算法的强烈程度,其中,较小的β值将导致总体更加强烈的平滑处理。
可以通过如方程3所示对欧拉-拉格朗日偏微分方程(PDE)离散化而解方程2:
方程3:
可以如方程4所示,通过每一投影的估计的探测光子的平均数量表示方程3:
方程4:
其中,ρsm(m,r)表示估计的探测光子的平均数量,m表示探测器的数量,r表示行数。通过取得初始测量结果f的平滑化版本导出方程4中的ρsm(m,r);然而未必使估计的光子平均数量平滑化,下文将对此予以更为详细的描述。由于所记录的测量结果的噪声方差与1/ρsm(m,r)成比例,因而采用ρsm(m,r)替代方程3中的v是合理的。
通过如方程4所示的换算(scaling),优选获得高噪声(低光子)区域内的小的总变差以及具有高光子计数值的区域内的原始图像接近度。注意,其他形式的借助β的探测光子换算也是有可能的。例如,在另一实施例中,β可以与观测角成比例。
探测光子的平均数量的估计可以如下文所述。
从具有ρ(m,r)的形式的输入2D投影开始,其中,所述数据处于记录衰减域内,与这一测量相关联的探测光子数量作为每一投影的探测光子平均数量ρsm(m,r)估计。
可以交替采用这一均值或平滑化的均值。平滑化均值可以促进缓解针对探测到的非常小的光子数量的泊松随机变量中的大的偏差。就后一种情况而言,采用(例如)对估计的探测光子数量求得的移动平均值使估计光子数量平滑化,如方程5所示:
方程5:
其中,j∈{m-n:m+n},k∈{r-n:r+n},nKernel表示平滑化内核,这里也可以设想其他平滑化方案。对于处于投影的边缘的探测器而言,可以外推出或者以其他方式确定填充必要的缓冲区的值。
可以在迭代固定点方案的基础上采用方程5解方程4,在所述方案中,根据方程6由前一更新图像i确定每一更新i+1:
方程6:
其中,ρeff(m,r)=β·ρsm(m,r),Wp表示权重,可以通过(例如)本领域已知的各种方式计算所述权重。
如方程7所示记录经降噪的探测光子:
方程7:
其中,ρ输出(m,r)表示所输出的记录降噪探测光子,ρdn(m,r)表示经降噪的探测光子,Nr(m,r)表示光子数量,s表示换算因子。方程7将降噪探测光子重新转换成供重建的衰减线积分。
上述降噪算法是一种在计算方面有效率的算法,例如,因为在采集时对投影数据实施一次降噪操作,继而进行单次滤波反投影重建,因而避免了最大似然迭代方案中的多重前投影和反投影的繁重计算负担。
此外,相对于不向投影应用降噪算法的配置而言,上述降噪算法为重建图像提供了显著的条纹降低,同时又不影响投影的良好的统计学特性。通过在解记录的投影域内应用所述降噪算法,能够缓解来自低光子计数的偏差。
图3示出了一种用于对投影数据降噪的示范性方法。
在302中,采集多个二维投影信号(投影数据)。可以通过系统100和/或其他成像系统生成所述投影。
在304中,针对每一投影估计探测光子的平均数量。
在306中,采用如文中所述的以下算法对所述投影的至少子集降噪:在噪声小时优选原始投影,在噪声大时优选具有小的总变差的投影。如上文指出的,没有必要对具有足够数量的光子的投影降噪。
在308中,将经降噪的投影转换为降噪衰减数据。
在310中,重建降噪衰减数据,以生成一幅或多幅图像。
可以通过计算机可读指令实现上述操作,计算机可读指令在被计算机处理器执行时,令处理器执行所描述的动作。在这种情况下,将所述指令存储在与相关计算机相关或者能够为其所访问的计算机可读存储介质中。
应当认识到,可以通过运行一个或多个计算机可读和/或可执行指令的一个或多个处理器实现所述投影数据降噪器114,其中,所述指令存储在或者被编码到诸如物理存储器的计算机可读存储介质上。或者或此外,可以将所述指令存储到诸如信号介质等的瞬时媒介上。
已经参考各种实施例描述了本发明。在阅读文中的描述的同时,本领域技术人员可以想到修改和变化。这意味着,应当将本发明推断为包括所有此类落在所附权利要求或者其等同要件的范围内的修改和变化。
Claims (20)
1.一种系统,包括:
源,其围绕检查区域旋转并且发射横穿所述检查区域的辐射;
辐射敏感探测器阵列,其探测横穿所述检查区域的辐射并生成指示所探测到的辐射的投影数据;以及
投影数据降噪器,其对所述投影数据降噪,其中,所述降噪器基于一投影的探测光子数量对所述投影进行降噪。
2.根据权利要求1所述的系统,所述降噪器包括:
光子估计器,其估计针对每一投影的探测光子数量,其中,所述投影数据降噪器基于所估计的探测光子数量对所述投影数据降噪。
3.根据权利要求2所述的系统,其中,所估计的探测光子数量表示对所估计的探测光子数量求得的平滑化移动平均值。
4.根据权利要求1所述的系统,所述投影数据降噪器还包括:
基于对代价函数求最小值以及所估计的探测光子数量对投影进行降噪的降噪器,所述代价函数包括至少两项:第一总变差项和第二加权最小二乘方项。
5.根据权利要求1所述的系统,所述投影数据降噪器还包括:
通过在所估计的探测光子数量大时优选所述投影,在所估计的探测光子数量小时优选具有小的总变差的投影而对投影进行降噪的降噪器。
6.根据权利要求5所述的系统,其中,降噪器采用迭代算法对所述代价函数的求最小值。
7.根据权利要求1所述的系统,还包括:
重建器,其重建经降噪的投影数据以生成体积图像数据。
8.根据权利要求1所述的系统,其中,相对于具有较高光子数量的投影而言,所述投影数据降噪器将在更大的程度对具有较低光子数量的投影降噪。
9.根据权利要求1所述的系统,所述投影数据降噪器还包括:
评估器,其基于对于投影的探测光子数量确定是否对所述投影进行降噪。
10.根据权利要求9所述的系统,其中,至少一个投影包含满足预定光子数量阈值的数量的探测光子,其不受到降噪,并且至少一个投影包含不满足所述预定光子数量阈值的数量的探测光子,其受到降噪。
11.根据权利要求1所述的系统,其中,所述投影数据降噪器在保持给定的图像分辨率的同时降低噪声。
12.一种方法,包括:
获得成像系统生成的投影数据;
针对所述投影数据的二维投影确定估计的探测光子数量,并生成指示估计的探测光子数量的信号;以及
基于所述信号对所述投影进行降噪。
13.根据权利要求12所述的方法,其中,降噪程度与所述投影的估计的探测光子数量成比例,并且在两个不同的程度上对至少两个投影进行降噪。
14.根据权利要求13所述的方法,其中,对所述投影降噪包括对代价函数求最小值,所述代价函数包括总变差项和加权最小二乘方项。
15.根据权利要求14所述的方法,其中,在所述估计的探测光子数量较低时,所述总变差项占据主导,在所述估计的探测光子数量较高时,所述加权最小二乘方项占据主导。
16.根据权利要求12所述的方法,其中,在所述估计的探测光子数量大时,所述代价函数优选所述投影,在所述估计的探测光子数量小时,所述代价函数优选具有小的总变差的投影。
17.根据权利要求12所述的方法,还包括:
仅对那些具有不满足预定阈值的估计光子数量的投影进行降噪,对其余的投影则基于相应的估计光子数量进行一定程度的降噪。
18.根据权利要求12所述的方法,其中,相对于具有较高的估计的探测光子数量的投影而言,将在更大的程度上对具有较低的估计的探测光子数量的投影进行降噪。
19.根据权利要求12所述的方法,还包括:
通过采用对所述估计的探测光子数量求得的移动平均值使所述估计的探测光子数量平滑化,其中,基于平滑化的估计的探测光子数量对所述投影进行降噪。
20.一种方法,包括:
对具有这样噪声的图像进行降噪:所述噪声在所述图像的相邻体素之间存在强相关性,其中,所述降噪不会向经降噪的图像内引入条纹伪像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US40704010P | 2010-10-27 | 2010-10-27 | |
US61/407,040 | 2010-10-27 | ||
PCT/IB2011/054588 WO2012056364A1 (en) | 2010-10-27 | 2011-10-17 | Low dose ct denoising |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103180875A true CN103180875A (zh) | 2013-06-26 |
CN103180875B CN103180875B (zh) | 2018-04-13 |
Family
ID=44993631
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180051297.9A Active CN103180875B (zh) | 2010-10-27 | 2011-10-17 | 低剂量ct降噪 |
Country Status (8)
Country | Link |
---|---|
US (1) | US9031299B2 (zh) |
EP (1) | EP2633494B1 (zh) |
JP (1) | JP5937093B2 (zh) |
CN (1) | CN103180875B (zh) |
BR (1) | BR112013009965A8 (zh) |
MX (1) | MX2013004540A (zh) |
RU (1) | RU2586968C2 (zh) |
WO (1) | WO2012056364A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104644203A (zh) * | 2014-09-02 | 2015-05-27 | 沈阳东软医疗系统有限公司 | 一种剂量调制扫描方法和装置 |
CN110574073A (zh) * | 2016-12-19 | 2019-12-13 | 皇家飞利浦有限公司 | 谱计算机断层摄影(ct)成像中的残余碘伪影的探测和/或校正 |
CN112656438A (zh) * | 2020-12-17 | 2021-04-16 | 中山大学 | 一种基于曲面全变差的低剂量ct投影域去噪及重建方法 |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538114B2 (en) * | 2011-06-06 | 2013-09-17 | Kabushiki Kaisha Toshiba | Method and system utilizing parameter-less filter for substantially reducing streak and or noise in computer tomography (CT) images |
DK2732437T3 (en) * | 2011-07-13 | 2017-12-11 | Scivis Wss Bildverarbeitung Gmbh | Reconstruction of transmission images as well as image reproduction using Poisson detector data |
EP2783344B1 (en) | 2011-11-23 | 2017-05-03 | Koninklijke Philips N.V. | Image domain de-noising |
WO2013088294A1 (en) | 2011-12-13 | 2013-06-20 | Koninklijke Philips Electronics N.V. | Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising |
US9105124B2 (en) * | 2012-06-21 | 2015-08-11 | General Electric Company | Methods and systems for reducing noise- related imaging artifacts |
CN104903933B (zh) * | 2012-11-26 | 2018-10-19 | 皇家飞利浦有限公司 | 投影数据降噪 |
KR101891833B1 (ko) * | 2012-12-31 | 2018-08-24 | 부산가톨릭대학교 산학협력단 | 전변분에 기초한 노이즈 제거 알고리즘을 이용한 유방 x선 조영장치에서의 선량 감소 방법 |
CN103226815A (zh) * | 2013-04-10 | 2013-07-31 | 东南大学 | 一种低剂量ct图像滤波方法 |
JP6214226B2 (ja) * | 2013-06-06 | 2017-10-18 | キヤノン株式会社 | 画像処理装置、断層撮影装置、画像処理方法およびプログラム |
KR101904493B1 (ko) * | 2016-11-22 | 2018-10-17 | 동서대학교산학협력단 | 저선량 촬영한 동적 3차원 의료영상을 시공간 특징 기반으로 개선하는 방법 |
CN108961237B (zh) * | 2018-06-28 | 2020-08-21 | 安徽工程大学 | 一种基于卷积神经网络的低剂量ct图像分解方法 |
EP3820372B1 (en) * | 2018-07-09 | 2023-09-06 | Koninklijke Philips N.V. | Imaging system detector clipping-induced bias correction |
US11054534B1 (en) | 2020-04-24 | 2021-07-06 | Ronald Nutt | Time-resolved positron emission tomography encoder system for producing real-time, high resolution, three dimensional positron emission tomographic image without the necessity of performing image reconstruction |
US11300695B2 (en) | 2020-04-24 | 2022-04-12 | Ronald Nutt | Time-resolved positron emission tomography encoder system for producing event-by-event, real-time, high resolution, three-dimensional positron emission tomographic image without the necessity of performing image reconstruction |
RU2744552C1 (ru) * | 2020-08-06 | 2021-03-11 | Государственное бюджетное учреждение здравоохранения города Москвы "Научно-практический клинический центр диагностики и телемедицинских технологий Департамента здравоохранения города Москвы" (ГБУЗ "НПКД ДиТ ДЗМ") | Способ исследования состояния легких при подозрении на COVID-19 с помощью низкодозной компьютерной томографии |
RU2753474C1 (ru) * | 2020-12-07 | 2021-08-17 | Государственное бюджетное учреждение здравоохранения города Москвы "Научно-практический клинический центр диагностики и телемедицинских технологий Департамента здравоохранения города Москвы" (ГБУЗ "НПКЦ ДиТ ДЗМ") | Способ низкодозного сканирования органов грудной клетки, адаптированный к массе тела пациента |
CN114793291B (zh) * | 2022-01-25 | 2023-11-07 | 深圳软牛科技有限公司 | Ios多设备实时投屏方法、系统、计算机设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6507633B1 (en) * | 2001-02-15 | 2003-01-14 | The Regents Of The University Of Michigan | Method for statistically reconstructing a polyenergetic X-ray computed tomography image and image reconstructor apparatus utilizing the method |
US20030219152A1 (en) * | 2002-02-14 | 2003-11-27 | Jonas August | Image smoothing with decoupled regularization |
CN101237820A (zh) * | 2005-04-25 | 2008-08-06 | 罗切斯特大学 | 用于ct成像的整体去噪的方法和装置 |
CN102203826A (zh) * | 2008-12-25 | 2011-09-28 | 梅迪奇视觉成像解决方案有限公司 | 医学图像的降噪 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4991093A (en) * | 1985-08-21 | 1991-02-05 | Exxon Research And Engineering Company | Method for producing tomographic images using direct Fourier inversion |
US5053615A (en) * | 1990-04-06 | 1991-10-01 | Samsung Electronics Co., Ltd. | Correction algorithm for contiguous CCD elements leakage |
US6016333A (en) * | 1998-05-15 | 2000-01-18 | International Business Machines Corporation | System and method for improving the quality of images produced by CT scanners operating on low power |
US6490476B1 (en) * | 1999-10-14 | 2002-12-03 | Cti Pet Systems, Inc. | Combined PET and X-ray CT tomograph and method for using same |
US7187794B2 (en) * | 2001-10-18 | 2007-03-06 | Research Foundation Of State University Of New York | Noise treatment of low-dose computed tomography projections and images |
US6493416B1 (en) * | 2001-11-21 | 2002-12-10 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for noise reduction in computed tomographic systems |
IL151634A0 (en) * | 2002-09-05 | 2003-04-10 | Real Time Radiography Ltd | Direct detection of high energy single photons |
US6950492B2 (en) * | 2003-06-25 | 2005-09-27 | Besson Guy M | Dynamic multi-spectral X-ray projection imaging |
US7317841B2 (en) * | 2003-12-22 | 2008-01-08 | Ge Medical Systems Global Technology Company, Llc | System and method for image noise reduction using a minimal error spatiotemporal recursive filter |
WO2006006090A1 (en) * | 2004-07-07 | 2006-01-19 | Koninklijke Philips Electronics, N.V. | Streak artifact reduction in cardiac cone beam ct reconstruction |
US7508968B2 (en) * | 2004-09-22 | 2009-03-24 | Siemens Medical Solutions Usa, Inc. | Image compounding based on independent noise constraint |
US7616841B2 (en) * | 2005-06-17 | 2009-11-10 | Ricoh Co., Ltd. | End-to-end design of electro-optic imaging systems |
JP2009513220A (ja) * | 2005-10-28 | 2009-04-02 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 分光コンピュータ断層撮影の方法および装置 |
WO2010011676A2 (en) * | 2008-07-23 | 2010-01-28 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
US7888651B2 (en) * | 2007-05-21 | 2011-02-15 | The Board Of Trustees Of The Leland Stanford Junior University | Method and system for using tissue-scattered coincidence photons for imaging |
US8014616B2 (en) | 2007-11-02 | 2011-09-06 | Siemens Aktiengesellschaft | System and method for fixed point continuation for total variation based compressed sensing imaging |
JP2010078338A (ja) * | 2008-09-24 | 2010-04-08 | Toshiba Corp | X線検出器 |
WO2010038536A1 (ja) * | 2008-09-30 | 2010-04-08 | 株式会社 日立メディコ | X線ct装置 |
WO2010096701A1 (en) * | 2009-02-20 | 2010-08-26 | Mayo Foundation For Medical Education And Research | Projection-space denoising with bilateral filtering in computed tomography |
US8031831B2 (en) * | 2009-05-28 | 2011-10-04 | Kabushiki Kaisha Toshiba | Voltage and or current modulation in dual energy computed tomography |
RU2543547C2 (ru) | 2009-10-06 | 2015-03-10 | Конинклейке Филипс Электроникс Н.В. | Способ уменьшения артефактов в изображениях конусно-лучевой компьютерной томографии |
BR112012012231A2 (pt) | 2009-11-25 | 2017-12-19 | Koninl Philips Electronics Nv | método e sistema |
-
2011
- 2011-10-17 MX MX2013004540A patent/MX2013004540A/es active IP Right Grant
- 2011-10-17 EP EP11784781.4A patent/EP2633494B1/en active Active
- 2011-10-17 CN CN201180051297.9A patent/CN103180875B/zh active Active
- 2011-10-17 BR BR112013009965A patent/BR112013009965A8/pt not_active IP Right Cessation
- 2011-10-17 WO PCT/IB2011/054588 patent/WO2012056364A1/en active Application Filing
- 2011-10-17 US US13/879,821 patent/US9031299B2/en active Active
- 2011-10-17 RU RU2013123937/14A patent/RU2586968C2/ru not_active IP Right Cessation
- 2011-10-17 JP JP2013535544A patent/JP5937093B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6507633B1 (en) * | 2001-02-15 | 2003-01-14 | The Regents Of The University Of Michigan | Method for statistically reconstructing a polyenergetic X-ray computed tomography image and image reconstructor apparatus utilizing the method |
US20030219152A1 (en) * | 2002-02-14 | 2003-11-27 | Jonas August | Image smoothing with decoupled regularization |
CN101237820A (zh) * | 2005-04-25 | 2008-08-06 | 罗切斯特大学 | 用于ct成像的整体去噪的方法和装置 |
CN102203826A (zh) * | 2008-12-25 | 2011-09-28 | 梅迪奇视觉成像解决方案有限公司 | 医学图像的降噪 |
Non-Patent Citations (2)
Title |
---|
LINDBLAD J. ET AL: "De-noising of SRμCT Fiber Images by Total Variation Minimization", 《2010 INTERNATIONAL CONFERENCE ON PATTERN RECOGNITION》, 26 August 2010 (2010-08-26), pages 4622 * |
张元科 等: "低剂量CT投影图像噪声分析及去噪算法研究", 《光电子•激光》, vol. 21, no. 7, 31 July 2010 (2010-07-31) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104644203A (zh) * | 2014-09-02 | 2015-05-27 | 沈阳东软医疗系统有限公司 | 一种剂量调制扫描方法和装置 |
CN110574073A (zh) * | 2016-12-19 | 2019-12-13 | 皇家飞利浦有限公司 | 谱计算机断层摄影(ct)成像中的残余碘伪影的探测和/或校正 |
CN112656438A (zh) * | 2020-12-17 | 2021-04-16 | 中山大学 | 一种基于曲面全变差的低剂量ct投影域去噪及重建方法 |
CN112656438B (zh) * | 2020-12-17 | 2023-02-21 | 中山大学 | 一种基于曲面全变差的低剂量ct投影域去噪及重建方法 |
Also Published As
Publication number | Publication date |
---|---|
RU2013123937A (ru) | 2014-12-10 |
BR112013009965A8 (pt) | 2016-09-13 |
BR112013009965A2 (pt) | 2016-08-02 |
EP2633494A1 (en) | 2013-09-04 |
JP5937093B2 (ja) | 2016-06-22 |
EP2633494B1 (en) | 2019-12-11 |
US20130208971A1 (en) | 2013-08-15 |
CN103180875B (zh) | 2018-04-13 |
WO2012056364A1 (en) | 2012-05-03 |
RU2586968C2 (ru) | 2016-06-10 |
MX2013004540A (es) | 2013-07-03 |
US9031299B2 (en) | 2015-05-12 |
JP2013542018A (ja) | 2013-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103180875A (zh) | 低剂量ct降噪 | |
US10395353B2 (en) | Model-based scatter in multi-modality multi-energy SPECT reconstruction | |
JP6100772B2 (ja) | 画像処理方法及びコンピューティング装置 | |
JP2023159080A (ja) | 医用画像処理装置及び医用画像処理システム | |
US8189735B2 (en) | System and method for reconstruction of X-ray images | |
JP5997180B2 (ja) | デュアルエネルギーct画像再構成のための方法及びシステム | |
JP6559678B2 (ja) | スペクトル画像データ処理 | |
EP2691932B1 (en) | Contrast-dependent resolution image | |
US10878544B2 (en) | Image data processing | |
CN102113020B (zh) | 用于谱x射线成像的系统和方法 | |
Remy-Jardin et al. | Thoracic helical CT: influence of subsecond scan time and thin collimation on evaluation of peripheral pulmonary arteries | |
CN103824266A (zh) | 具有增强噪声控制滤波的迭代重建 | |
CN102947861A (zh) | 用于在低剂量计算机断层摄影中降低噪声的方法和系统 | |
CN102048552A (zh) | 计算机断层造影灌注测量的射线硬化校正 | |
US20150043795A1 (en) | Image domain pansharpening method and system for spectral ct with large pixel energy discriminating detectors | |
CN103027705A (zh) | 产生运动补偿的ct 图像数据组的方法和系统 | |
CN103717137A (zh) | 图像生成装置 | |
JP6828061B2 (ja) | スペクトルctデータからのカルシウム含有量の特定 | |
Bailey et al. | Experience with fully 3D PET and implications for future high-resolution 3D tomographs | |
Watson | Count rate dependence of local signal-to-noise ratio in positron emission tomography | |
CN102096905A (zh) | 双源计算机断层成像拍摄中的降噪 | |
CN102670230B (zh) | 用于在双能计算机断层造影中减少运动伪影的方法 | |
Tadesse et al. | Techniques for generating attenuation map using cardiac SPECT emission data only: a systematic review | |
JP6425885B2 (ja) | 核医学診断装置、画像処理装置および画像処理プログラム | |
US11826190B2 (en) | System and method for quantitative blood volume imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |