CN101061504B - 通过迭代约束去卷积进行核医学2d平面图像的恢复 - Google Patents
通过迭代约束去卷积进行核医学2d平面图像的恢复 Download PDFInfo
- Publication number
- CN101061504B CN101061504B CN2005800393974A CN200580039397A CN101061504B CN 101061504 B CN101061504 B CN 101061504B CN 2005800393974 A CN2005800393974 A CN 2005800393974A CN 200580039397 A CN200580039397 A CN 200580039397A CN 101061504 B CN101061504 B CN 101061504B
- Authority
- CN
- China
- Prior art keywords
- image
- initial
- function
- radiation
- algorithm
- 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.)
- Expired - Fee Related
Links
- 238000009206 nuclear medicine Methods 0.000 title description 4
- 230000005855 radiation Effects 0.000 claims abstract description 37
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000005316 response function Methods 0.000 claims abstract description 9
- 238000002059 diagnostic imaging Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 55
- 238000000034 method Methods 0.000 claims description 29
- 230000004044 response Effects 0.000 claims description 20
- 238000011084 recovery Methods 0.000 claims description 16
- 238000007476 Maximum Likelihood Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 7
- 201000004569 Blindness Diseases 0.000 claims description 6
- 238000009877 rendering Methods 0.000 claims description 5
- 230000000155 isotopic effect Effects 0.000 claims description 2
- 239000007858 starting material Substances 0.000 claims 3
- 238000003384 imaging method Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 210000000056 organ Anatomy 0.000 description 6
- 238000002603 single-photon emission computed tomography Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 229940121896 radiopharmaceutical Drugs 0.000 description 5
- 239000012217 radiopharmaceutical Substances 0.000 description 5
- 230000002799 radiopharmaceutical effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 230000002285 radioactive effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 239000013078 crystal Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012633 nuclear imaging Methods 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000009513 drug distribution Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 230000005610 quantum mechanics Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography 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
- G06T5/73—Deblurring; Sharpening
-
- 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
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Nuclear Medicine (AREA)
Abstract
医学成像系统(10)包括至少一个放置在对象接纳孔(18)附近以检测来自对象的辐射的辐射检测头(16)。将检测到的辐射重建成至少一个初始2D投影图像(μ)。利用拓展的迭代约束去卷积算法,通过合并对应于检测头和检测到的辐射发源点之间不同距离的系统响应函数的不同估计来恢复每个初始2D图像(μ)中的分辨率。所测得的响应函数用于恢复一系列图像。通过以品质因数进行自动搜索,通过用户观看,或者通过利用盲目去卷积进行系统响应函数的并行估计和更新原始图像来确定最佳图像。
Description
本发明涉及核医学诊断成像系统和方法。其尤其与单光子发射断层摄影(SPECT)系统联合应用,并将具体参照该系统进行描述。应当意识到,本发明还可应用于其它成像系统,如正电子发射断层摄影系统(PET)、计算机断层摄影系统(CT)和X射线成像等。
核成像采用放射性源对患者的解剖学结构进行成像。典型地,将放射性药物注入患者。放射性药物化合物包含以可预知速率和特征能量进行伽马辐射衰减的放射性同位素。一个或多个辐射检测器放置在患者附近用于监视和记录发射出的能量。有时,检测器绕患者旋转或转换角度以监视从多个方向发射的辐射。基于例如所检测到的位置和能量的信息,确定体内放射性药物分布并重建该分布的图像,以研究循环系统、选定器官或组织中的放射性药物吸收等。
在传统的闪烁检测器中,检测器具有由大闪烁晶体或较小闪烁晶体的矩阵组成的闪烁体。在无论哪种情况下,都是由传感器矩阵观察闪烁体。常用的传感器为光电倍增管(“PMT”)。包括栅格状或蜂窝状辐射吸收材料阵列的准直器位于闪烁体和受检查对象之间,以限制撞击到闪烁体上的辐射的接受角度。每个撞击到闪烁体上的辐射事件产生由PMT检测到的闪光(闪烁)。基于PMT的输出,伽马照相机映射辐射事件,即,其确定撞击到闪烁体上的射线的能量和位置。
SPECT图像的图像质量由检测器的计数灵敏度和准直器的几何形状确定。通常,很难获得高质量的SPECT图像,这是因为由包括系统参数在内的各种因素导致的有限空间分辨率,系统参数例如是准直器几何形状、非线性PMT响应和量子力学概率等。图像模糊或恶化通常表示为点扩散函数(PSF)。
已提出的用于提高平面(2D)图像分辨率的恢复技术采用经典的反向滤波器,如维纳滤波器、依赖于计数的梅斯滤波器、基于最大熵的滤波器和功率谱均衡滤波器等。然而,所提出的反向滤波技术假设点扩散函数是不依赖于深度的,例如,PSF是已知的并且定义在确定的深度处。已将恢复方法中的一种(迭代约束去卷积)扩大到采用未知点扩散函数或所谓的盲目去卷积,以通过采用投影信息获得投影之间的相关性来重建三维SPECT图像。然而,对于2D核医学平面图像,这种相关性信息不存在。
核医学图像的恢复是复杂的,主要由于图像模糊或点扩散函数取决于辐射衰减事件和照相机之间的距离这一事实,例如它取决于深度。所得到图像的准确点扩散函数通常是未知的,因为图像是来自不同深度的信息的汇编。此外,感兴趣器官的深度对于不同的研究会发生变化,从而在平面图像情况下难于预知给定研究的点扩散函数。
本发明提供一种克服上述问题和其它问题的新颖和改进的成像设备和方法。
根据本发明的一个方面,公开一种医学成像系统。至少一个辐射检测头放置在对象接纳孔附近以检测来自对象的辐射。一个装置将检测到的辐射重建成至少一个初始2D投影图像。迭代约束去卷积装置用多个系统响应函数迭代地恢复每个初始2D图像的分辨率,每个系统响应函数表示检测头和检测到的辐射的发源点之间的相应距离。
根据本发明的另一方面,公开一种医学成像方法。检测来自对象的辐射数据。将所检测到的辐射重建成初始2D图像。通过采用迭代约束去卷积迭代地恢复初始2D图像的分辨率,该迭代约束去卷积应用多个系统响应函数,每个系统响应函数对应于检测头和检测到的辐射的发源点之间的距离。
本发明的一个优点在于增加的空间分辨率。
另一个优点在于减少2D图像中的噪声。
本发明的另一个优点在于增加的重建图像的精度。
在阅读和理解下列优选实施例的详细说明后,对本领域技术人员而言,本发明的其它优点和有益效果将变得更为明显。
本发明可采用各种部件和部件设置,以及各种步骤和步骤设置形式。附图仅用于解释说明优选实施例的目的,而不应解释为限制本发明。
图1是成像系统的示意表示;
图2示意性示出三个点的点扩散函数;
图3是成像系统一部分的示意表示;和
图4是成像系统另一部分的示意表示。
参见图1,核成像系统10典型地包括支撑可旋转台架14的静止台架12。一个或多个检测器头16由旋转台架14承载以检测从感兴趣区域或检查区域18发射的辐射事件。每个检测器头包括检测器元件或闪烁体20的二维阵列。每个头16包括用于将每个辐射响应转换成表示其在检测器面上的位置(x,y)和其能量(z)的数字信号的电路22。在二维(2D)笛卡尔坐标系统内对闪烁体20上的事件位置以通称的x和y坐标进行分解和/或确定。然而,也可采用其它坐标系。准直器24控制方向和角度扩展,闪烁体20的每个检测器元件可在该方向和角度下接收辐射,即,闪烁体20可仅沿已知射线接收辐射。这样,所确定的在闪烁体20上检测到辐射的位置以及照相机16的角位置定义了每个辐射事件发生所沿的法线。由于准直器有限的高度、其孔径的物理尺寸、检测物理过程和其它系统参数,存在从辐射事件传播到检测器的实际射线的不确定性。潜在的射线随着离开检测器的距离而扩散,其中在由点扩散函数(PSF)给定的任何给定距离处具有位置不确定性。
典型地,将一种或多种放射性药物或放射性同位素注入待成像对象并将待成像对象放置在由床26支撑的检查区域18中。这种同位素的几个实例是Tc-99m、Ga-67和In-131。对象内存在的放射性药物从对象发出辐射。检测器头16检测该辐射,其能够绕检查区域18成角度地转换角度或旋转以采集一个或多个选定投影方向上的投影发射数据。投影发射数据,例如,位置(x,y)和/或能量(z),以及每个检测器头16围绕检查区域18的角位置(θ)(例如从角位置解算器28获得的)存储在数据存储器30内。来自每个选定投影方向的投影数据被传送至重建投影图像的图像处理器40,投影图像存储在投影图像存储器42内。
2D分辨率恢复处理器44恢复2D图像,如由全身骨扫描、肿瘤扫描和其它扫描生成的2D图像。如下面更为详细描述的,该图像通过去卷积而变得更清晰或再聚焦。模糊或点扩散函数取决于辐射事件和检测器20之间的距离。在图2中,三个点源a、b、c示出在不同深度d1、d2、d3处。在相同的检测器区采集来自点b和c的辐射数据。然而,点b的点扩散函数hb比点c的点扩散函数hc更宽。这表明点扩散函数依赖于深度,例如,从点源b到检测器20的距离d2比从点源c至检测器20的距离d3大。
再次参照图1,分辨率恢复2D图像存储在2D图像存储器48内。视频处理器50处理优化的2D图像以供显示在图像显示器52上。
在一个实施例中,3D重建处理器60将来自2D图像存储器48的2D图像处理成容积图像表示。这些图像表示存储在3D图像存储器62中以由视频处理器50操作并显示在图像显示器52上,如视频监视器、打印机等。
继续参照图1并进一步参照图3,迭代处理器或处理过程70通过利用具有已知深度的点扩散函数向原始图像μ施加迭代约束去卷积算法或处理过程。更具体地,函数装置72基于来自信息数据库74的输入为不同的深度确定一系列点扩散函数h1,h2,...,hn。优选地,基于给予对象的至少一种同位素和用于成像的准直器24的几何形状来确定点扩散函数h1,h2,...,hn。当然,还考虑到用户可通过操作者接口站78的输入装置76(如键盘和鼠标等)在工作中(on the fly)提供确定点扩散函数h1,h2,...,hn的信息。固定装置或处理器80针对每个点扩散函数(h1,h2,...,hn)向原始图像μ施加迭代约束去卷积算法或处理过程(ICD),以生成一系列恢复图像λ1,λ2,...,λμ。
优选地,固定处理器80实施最大似然(ML)算法:
其中h是已知的点扩散函数,
z,1表示图像的像素指数(index),
μ表示原始图像,
λ(k)是在第k次迭代时的恢复图像,和
i表示图像中第i个像素。
终止装置82基于预先指定的标准来终止迭代过程70,该预先指定的标准例如预先选定的迭代数、至少在一个选定的图像区域内预先定义的平均对比度、以及至少在一个选定的图像区域内预先定义的标准像素偏差等。
图像优化装置84,优选采用预先定义的优化分辨率标准自动地使优化图像或从恢复图像系列λ1,λ2,...,λμ选出最佳图像,例如通过使对比度噪声比、信噪比或某个频率分量最大化而使这些图像优化。恢复的2D图像存储在2D图像存储器48内。
当然,还考虑到由用户凭借视觉进行优化。在使初始图像对于每个点扩散函数h进行恢复或再聚焦后,每个恢复的图像存储在2D图像存储器48内。每个对应于不同深度最佳再聚焦的恢复图像系列显示在显示器52上。优选地,用户在该系列上滚动以选择他/她的观看偏好。在该过程期间的显示类似于基于需要观察的身体部分的深度而聚焦显微镜。例如,如果感兴趣器官位于靠近照相机16处,具有小深度的点扩散函数导致更佳质量的图像;而在其它情况下,当感兴趣器官远离照相机16时,具有大深度的点扩散函数导致更佳质量图像。
继续参照图1并进一步参照图4,迭代处理器70通过采用假设或未知的点扩散函数h1向原始图像μ施加迭代约束去卷积算法。更具体地,函数装置72估计点扩散函数h1的近似值。由于点扩散函数h1的全部知识是未知的,盲目去卷积装置或过程或处理器86从近似的点扩散函数h1开始对原始图像μ施加所谓的盲目去卷积。重新定义函数处理器或装置88估计点扩散函数,而重新定义图像处理器或装置90同时恢复图像。终止装置82通过监视盲目去卷积过程86确定恢复图像λ是否具有可接受的质量。在优选实施例中,终止装置82确定图像的平均对比度是否在预定义阈值以上、标准像素偏差是否在预定义标准像素偏差以上和是否已进行预选定数目迭代中的一个。最终的最佳恢复图像λ存储在2D图像存储器48中。优选地,盲目去卷积装置86在不放大太多噪声和/或不产生错误的伪影特征的情况下提高图像的对比度。这种盲目去卷积过程86是有利的,因为无需知道点扩散函数的准确信息。盲目恢复的点扩散函数比凭借经验测量的点扩散函数更精确,因为它没有噪声污染,并且是基于实际存在于数据集内的失真,而不是当测量点扩散函数的时候。结果,采用盲目恢复点扩散函数恢复的图像更为稳固和在统计学上更精确。
优先地,盲目去卷积过程86实施最大似然(ML)算法,其中图像λ和点扩散函数h在每次迭代中同时得到估计:
其中h(k)是第k次迭代时的估计点扩散函数;
z,1表示图像的像素索引,
μ表示原始图像,
λ(k)是第k次迭代的恢复图像,和
i表示图像中的第i个像素。
去卷积过程已参照整个投影图像进行简单描述。当然,可基于区域性使单个图像优化,其中使多个区域中的每一个对不同点扩散函数都是最佳的。这样,具有较浅感兴趣器官的区域和具有较深感兴趣器官的相同图像的另一区域都可以得到最佳的和较好分辨率的恢复。
迭代约束去卷积技术可用于恢复3D图像中的分辨率。在一个实施例中,通过采用3D容积响应函数在3D图像重建后对3D图像施加迭代约束去卷积。在另一个实施例中,在SPECT重建之前对每个2D投影施加迭代去卷积。在又一实施例中,迭代去卷积并入3D SPECT重建过程中。
本发明已参照优选实施例进行描述。显然,在阅读和理解前面详细描述后,其他人可进行修改和替换。本发明应当解释为包括全部这些修改和替换,只要它们落入所附的权利要求书及其等效物的范围内。
Claims (23)
1. 一种医学成像系统(10),包括:
至少一个辐射检测头(16),其放置在对象接纳孔(18)附近以检测来自对象的辐射;
用于将检测到的辐射重建成至少一个初始2D投影图像(μ)的装置(40);
迭代约束去卷积装置(70),其用多个系统响应函数(h)迭代地恢复每个初始2D图像的分辨率,每个系统响应函数表示检测头(16)和检测到的辐射的发源点之间的相应距离。
2. 根据权利要求1所述的系统,其中迭代装置(70)包括:
用于向初始2D图像(μ)施加盲目迭代约束去卷积算法的盲目装置(86),其中在图像分辨率恢复去卷积的迭代反复中调整初始系统响应函数(h)。
3. 根据权利要求2所述的系统,其中盲目装置(86)还包括:
用于迭代地重新定义2D图像的装置(90);和
用于迭代地重新定义系统响应函数的装置(88),使得图像重新定义装置(90)迭代地恢复2D图像,以及函数重新定义装置(88)同时迭代地重新定义系统响应函数。
4. 根据权利要求3所述的系统,还包括:
用于基于以下的至少一种终止由盲目装置(86)重新定义2D图像和系统响应函数的装置(82):
预先选定的迭代数,
至少在图像的一个选定区域中预先定义的平均对比度,和
至少在图像的一个选定区域中预先定义的标准像素偏差。
5. 根据权利要求2所述的系统,其中盲目去卷积施加统计估计算法以同时估计恢复的2D图像(λ)和响应函数(h),该统计估计算法是以下的至少一种:
最大似然,
最大后验,
共轭梯度,和
贝叶斯。
6. 根据权利要求1所述的系统,其中所述多个系统响应函数(h1,h2,...,hn)由函数装置(72)产生,其中每个系统响应函数(h1,h2,...,hn)对应于检测头(16)和检测到的辐射的发源点之间的多个距离的每一个,并且迭代装置(70)包括:
用于对于每个系统响应函数(h1,h2,...,hn)向初始2D图像(μ)施加固定迭代约束去卷积算法以产生初始2D图像(μ)的一系列恢复的2D图像(λ1,λ2,...,λμ)的装置(80)。
7. 根据权利要求6所述的系统,还包括:
用于基于以下的至少一种终止施加固定迭代约束去卷积算法步骤的装置(82):
预先选定的迭代数,
至少在图像的一个选定区域中预先定义的平均对比度,和
至少在图像的一个选定区域中预先定义的标准像素偏差。
8. 根据权利要求7所述的系统,还包括:
用于优化恢复的2D图像系列(λ1,λ2,...,λμ)以选定最佳分辨率恢复图像的装置(84);和
用于以用户可观看格式显示最佳分辨率恢复图像的监视器(52)。
9. 根据权利要求8所述的系统,其中优化装置(84)通过以下其中一种优化恢复的图像系列(λ1,λ2,...,λμ):
使对比度噪声比最大化,和
使信噪比最大化。
10. 根据权利要求7所述的系统,还包括:
用户输入装置(76)和显示器(52),用户通过这些装置在恢复的图像系列(λ1,λ2,...,λμ)上滚动。
11. 根据权利要求7所述的系统,还包括:
用于向函数装置(72)提供信息的信息数据库(74),函数装置(72)基于该信息产生多个系统响应函数(h1,h2,...,hn),并且其中该信息包括以下至少一种:
给予对象的同位素的名称,和
准直器(24)的几何形状。
12. 根据权利要求6所述的系统,其中固定迭代约束去卷积施加统计估计算法以估计恢复的2D图像(λ),该统计估计算法是以下的至少一种:
最大似然,
最大后验,
共轭梯度,和
贝叶斯。
13. 一种医学成像方法,包括:
检测来自对象的辐射数据;
将所检测到的辐射重建成初始2D图像(μ);和
通过采用迭代约束去卷积迭代地恢复初始2D图像(μ)中的分辨率,该迭代约束去卷积施加多个系统响应函数(h),这些函数对应于检测头和检测到的辐射的发源点之间的距离。
14. 根据权利要求13所述的方法,还包括任意指定初始系统响应函数,并且其中迭代恢复步骤包括:
从初始系统响应函数开始向初始2D图像(μ)施加盲目迭代约束去卷积算法;和
同时迭代地重新定义去卷积图像和系统响应函数。
15. 根据权利要求14所述的方法,还包括:
基于以下的至少一种终止重新定义2D图像和系统响应函数的步骤:
预先选定的迭代数,
至少在图像的一个选定区域中预先定义的平均对比度,和
至少在图像的一个选定区域中预先定义的标准像素偏差。
16. 根据权利要求14所述的方法,其中恢复步骤还包括:
施加统计估计算法,该算法包括最大似然、最大后验、共轭梯度和贝叶斯中的至少一种,以同时估计更新的2D图像(λ)和更新的响应函数(h)。
17. 根据权利要求13所述的方法,其中迭代分辨率恢复步骤包括:
对于多个系统响应函数(h1,h2,...,hn)中的每一个向初始2D图像(μ)施加固定迭代约束去卷积算法,以产生一系列恢复的2D图像(λ1,λ2,...,λn),其中每个系统响应函数对应于检测头和检测到的辐射的发源点之间的多个距离中的一个。
18. 根据权利要求17所述的方法,其中迭代分辨率恢复步骤施加统计估计算法,该算法包括最大似然、最大后验、共轭梯度和贝叶斯中的至少一种,以估计更新的2D图像(λ)。
19. 根据权利要求17所述的方法,还包括:
基于以下的至少一种终止施加固定迭代约束去卷积算法的步骤:
预先选定的迭代数,
至少在图像的一个选定区域中预先定义的平均对比度,和
至少在图像的一个选定区域中预先定义的标准像素偏差。
20. 根据权利要求19所述的方法,还包括:
通过以下至少一种来优化恢复的2D图像系列(λ1,λ2,...,λn):
使对比度噪声比最大化,和
使信噪比最大化。
21. 根据权利要求20所述的方法,还包括:
以人类可读格式显示恢复的2D图像系列。
22. 根据权利要求21所述的方法,其中显示包括:
按照最佳分辨率恢复的有效深度的顺序连续地显示图像;和
在操作者的控制下,在图像系列上前后移动以选定最佳图像。
23. 根据权利要求13所述的方法,还包括:
在对象周围的多个角位置处检测辐射;
重建对应于每个角位置的初始2D图像;
对每个初始2D图像重复迭代分辨率恢复步骤;
对每个角位置选定最佳分辨率恢复2D图像;和
由最佳分辨率恢复2D图像重建3D图像表示。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US62863804P | 2004-11-17 | 2004-11-17 | |
US60/628,638 | 2004-11-17 | ||
PCT/IB2005/053629 WO2006054192A1 (en) | 2004-11-17 | 2005-11-04 | Restoration of the nuclear medicine 2d planar image by iterative constrained deconvolution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101061504A CN101061504A (zh) | 2007-10-24 |
CN101061504B true CN101061504B (zh) | 2011-07-06 |
Family
ID=35841527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2005800393974A Expired - Fee Related CN101061504B (zh) | 2004-11-17 | 2005-11-04 | 通过迭代约束去卷积进行核医学2d平面图像的恢复 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7840052B2 (zh) |
EP (1) | EP1828977B1 (zh) |
JP (1) | JP5172347B2 (zh) |
CN (1) | CN101061504B (zh) |
WO (1) | WO2006054192A1 (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008073208A (ja) * | 2006-09-21 | 2008-04-03 | Konica Minolta Medical & Graphic Inc | 画像処理装置及び画像処理方法 |
US9575186B2 (en) * | 2008-01-18 | 2017-02-21 | Koninklijke Philips N.V. | Multi-segment reconstruction |
ATE531006T1 (de) * | 2008-01-29 | 2011-11-15 | Bergen Teknologioverforing As | Iterative blindentfaltung unter nebenbedingungen |
EP2347282A2 (en) * | 2008-10-10 | 2011-07-27 | Koninklijke Philips Electronics N.V. | Practical spect calibration method for quantification of nuclides with high-energy contributions |
JP2011002306A (ja) * | 2009-06-18 | 2011-01-06 | Shimadzu Corp | Pet装置の逐次近似画像再構成法 |
JP5398381B2 (ja) * | 2009-06-26 | 2014-01-29 | 株式会社東芝 | 核医学イメージング装置および画像処理プログラム |
CN102576467B (zh) * | 2009-07-14 | 2015-11-25 | 皇家飞利浦电子股份有限公司 | 包括移变模糊补偿的图像重建 |
JP5815573B2 (ja) | 2010-03-18 | 2015-11-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 機能的画像データ強調方法及びエンハンサ |
US8849002B2 (en) * | 2010-07-07 | 2014-09-30 | The Board Of Trustees Of The Leland Stanford Junior University | Noise robust decoder for multiplexing readout channels on an imaging sensor array |
US8971599B2 (en) * | 2010-12-20 | 2015-03-03 | General Electric Company | Tomographic iterative reconstruction |
EP2685284A4 (en) * | 2011-03-10 | 2014-11-05 | Shimadzu Corp | POSITRON TOMOGRAPHY DEVICE AND RECONSTRUCTION METHOD |
JP5883689B2 (ja) * | 2012-03-14 | 2016-03-15 | 株式会社日立製作所 | X線撮像装置およびx線撮像方法 |
CN103279926A (zh) * | 2013-05-15 | 2013-09-04 | 中国航空工业集团公司沈阳空气动力研究所 | 一种tsp/psp旋转部件测量的模糊修正方法 |
US10165263B2 (en) | 2013-09-30 | 2018-12-25 | Nikon Corporation | Point spread function estimation of optics blur |
WO2015070008A1 (en) | 2013-11-08 | 2015-05-14 | Schlumberger Canada Limited | Spectral analysis with spectrum deconvolution |
RU2632249C1 (ru) | 2013-11-08 | 2017-10-03 | Шлюмбергер Текнолоджи Б.В. | Определение режима течения для адаптации модели потока |
KR101850871B1 (ko) * | 2015-08-26 | 2018-04-23 | 주식회사 디알텍 | 방사선 영상의 처리방법 및 방사선 촬영시스템 |
US11576641B2 (en) | 2016-01-12 | 2023-02-14 | CFL Research and Development, L.L.C. | Device and method for performing nuclear imaging |
CN105913394B (zh) * | 2016-04-08 | 2019-04-12 | 南京理工大学 | 一种基于退化图像序列的图像盲复原方法 |
US9905044B1 (en) | 2016-08-25 | 2018-02-27 | General Electric Company | Systems and methods for functional imaging |
CN106530381B (zh) * | 2016-10-19 | 2019-01-29 | 浙江大学 | 一种基于gpu加速的三维荧光显微图像的去卷积算法 |
JP6981799B2 (ja) * | 2017-07-21 | 2021-12-17 | キヤノンメディカルシステムズ株式会社 | 画像データ復元装置、画像データ復元方法、及びプログラム |
US10803633B2 (en) | 2018-02-06 | 2020-10-13 | General Electric Company | Systems and methods for follow-up functional imaging |
US11402331B2 (en) | 2019-05-08 | 2022-08-02 | City University Of Hong Kong | Imaging and manipulation of biological sample |
US11309072B2 (en) | 2020-04-21 | 2022-04-19 | GE Precision Healthcare LLC | Systems and methods for functional imaging |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5331553A (en) * | 1992-04-15 | 1994-07-19 | Ugm Medical Systems, Inc. | Three dimensional image reconstruction for a positron emission tomograph |
US6171243B1 (en) * | 1997-05-30 | 2001-01-09 | Picker International, Inc. | Combination of collimated and coincidence information for positron imaging |
CN1489976A (zh) * | 2002-09-23 | 2004-04-21 | 西门子公司 | 运动检查对象的计算机断层造影中设置图像的方法及设备 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5047968A (en) * | 1988-03-04 | 1991-09-10 | University Of Massachusetts Medical Center | Iterative image restoration device |
US5596198A (en) * | 1994-04-22 | 1997-01-21 | The Regents, University Of California | Gamma ray camera |
US5532490A (en) * | 1994-12-27 | 1996-07-02 | The University Of Utah | Displaced center-of-rotation fan-beam tomography for cardiac imaging |
US5565684A (en) * | 1995-06-30 | 1996-10-15 | The University Of Utah | Three-dimensional SPECT reconstruction of combined cone-beam and fan-beam data |
US6381349B1 (en) * | 1997-11-12 | 2002-04-30 | The University Of Utah | Projector/backprojector with slice-to-slice blurring for efficient 3D scatter modeling |
US6539103B1 (en) * | 1997-11-12 | 2003-03-25 | The University Of Utah | Method and apparatus for image reconstruction using a knowledge set |
AU4900801A (en) * | 1999-10-22 | 2001-06-25 | Arch Development Corporation | Accurate image reconstruction for doi-pet systems |
US6388244B1 (en) * | 2000-03-20 | 2002-05-14 | Philips Medical Systems (Cleveland), Inc. | Virtual contouring for transmission scanning in spect and pet studies |
JP4389371B2 (ja) * | 2000-09-28 | 2009-12-24 | 株式会社ニコン | 画像修復装置および画像修復方法 |
US6670614B1 (en) * | 2001-06-01 | 2003-12-30 | Leonard F. Plut | Volume cone beam acquisition on a nuclear spect system using a digital flat panel |
WO2003046753A1 (de) * | 2001-11-28 | 2003-06-05 | Leica Microsystems Wetzlar Gmbh | Verfahren zur geführten 'blind deconvolution' mikroskopischer bilder und software |
SE0201249D0 (sv) | 2002-04-22 | 2002-04-22 | Karolinska Innovations Ab | Methods and Arrangements for biological imaging |
US7197193B2 (en) * | 2002-05-03 | 2007-03-27 | Creatv Microtech, Inc. | Apparatus and method for three dimensional image reconstruction |
JP4686147B2 (ja) * | 2003-07-31 | 2011-05-18 | 株式会社東芝 | 画像データ処理装置 |
US7139067B2 (en) * | 2003-09-12 | 2006-11-21 | Textron Systems Corporation | Three-dimensional imaging with multiframe blind deconvolution |
JP4519434B2 (ja) * | 2003-09-24 | 2010-08-04 | 株式会社東芝 | 超解像処理装置及び医用画像診断装置 |
US7602989B2 (en) * | 2004-05-26 | 2009-10-13 | Biggs David S C | Realtime 2D deconvolution system and method |
-
2005
- 2005-11-04 WO PCT/IB2005/053629 patent/WO2006054192A1/en active Application Filing
- 2005-11-04 US US11/719,420 patent/US7840052B2/en not_active Expired - Fee Related
- 2005-11-04 CN CN2005800393974A patent/CN101061504B/zh not_active Expired - Fee Related
- 2005-11-04 EP EP05800448.2A patent/EP1828977B1/en not_active Not-in-force
- 2005-11-04 JP JP2007540776A patent/JP5172347B2/ja not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5331553A (en) * | 1992-04-15 | 1994-07-19 | Ugm Medical Systems, Inc. | Three dimensional image reconstruction for a positron emission tomograph |
US6171243B1 (en) * | 1997-05-30 | 2001-01-09 | Picker International, Inc. | Combination of collimated and coincidence information for positron imaging |
CN1489976A (zh) * | 2002-09-23 | 2004-04-21 | 西门子公司 | 运动检查对象的计算机断层造影中设置图像的方法及设备 |
Non-Patent Citations (6)
Title |
---|
HATZINAKOSD.IEEESIGNALSPROCESSINGMAGAZINE.Blindimagedeconvolutionvol.13no.3.1996 vol.13(no.3) |
KUNDUR D |
KUNDUR D;HATZINAKOS D.IEEE SIGNALS PROCESSING MAGAZINE.Blind image deconvolutionvol.13 no.3.1996,vol.13(no.3),第43-64、50、表9. * |
MIGNOTTE M ET AL.COMPARISON OF DECONVOLUTION TECHNIQUESUSING A DISTRIBUTION MIXTURE PARAMETERESTIMATION APPLICATION IN SINGLE PHOTONEMISSION COMPUTED TOMOGRAPHY IMAGERY.JOURNAL OF ELECTRONIC IMAGING,SPIE/IS & T,USVOL.11 NO.1.2002,VOL.11(NO.1),PAGES 11-24,摘要,附图6,第0001、2.3.1、03.4段. |
MIGNOTTE M ET AL.COMPARISON OF DECONVOLUTION TECHNIQUESUSING A DISTRIBUTION MIXTURE PARAMETERESTIMATION APPLICATION IN SINGLE PHOTONEMISSION COMPUTED TOMOGRAPHY IMAGERY.JOURNAL OF ELECTRONIC IMAGING,SPIE/IS & T,USVOL.11 NO.1.2002,VOL.11(NO.1),PAGES 11-24,摘要,附图6,第0001、2.3.1、03.4段. * |
US 6171243 B1,全文. |
Also Published As
Publication number | Publication date |
---|---|
US20090202125A1 (en) | 2009-08-13 |
CN101061504A (zh) | 2007-10-24 |
EP1828977B1 (en) | 2013-09-11 |
EP1828977A1 (en) | 2007-09-05 |
US7840052B2 (en) | 2010-11-23 |
WO2006054192A1 (en) | 2006-05-26 |
JP5172347B2 (ja) | 2013-03-27 |
JP2008530517A (ja) | 2008-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101061504B (zh) | 通过迭代约束去卷积进行核医学2d平面图像的恢复 | |
CN1106579C (zh) | 医学成象系统中多形态数据的互补重合 | |
Dougherty | Digital image processing for medical applications | |
US6381349B1 (en) | Projector/backprojector with slice-to-slice blurring for efficient 3D scatter modeling | |
CN103824266B (zh) | 具有增强噪声控制滤波的迭代重建 | |
CN101305297B (zh) | 使用解剖列表模式掩模的pet成像 | |
US11179128B2 (en) | Methods and systems for motion detection in positron emission tomography | |
JP2000180550A (ja) | Ml―em画像再構成法及び医用画像形成装置 | |
JPH0679343B2 (ja) | 作像システムにおける図の間の動きを補正する装置及び方法 | |
CN101849247A (zh) | 用于生成衰减图的设备和方法 | |
CN110660111A (zh) | Pet散射校正及图像重建方法、装置及设备 | |
US8294109B2 (en) | Extracting location information using difference images from a non-parallel hole collimator | |
JP7254322B2 (ja) | 医用画像処理システムとその画像再構成方法のためのコリメータ | |
US5969358A (en) | Whole body scan coincidence imaging | |
WO1993016435A1 (en) | Method of recovering signals in tomographic measurement systems | |
US8359345B2 (en) | Iterative algorithms for variance reduction on compressed sinogram random coincidences in PET | |
EP3588145A1 (en) | Photon scatter imaging | |
KR101493683B1 (ko) | 콘-빔 기반 반응선 재구성을 이용한 초고해상도 pet 영상 재구성 장치 및 그 방법 | |
KR20230131823A (ko) | 방사선 검출 시스템 및 방법 | |
Mejia et al. | Methodological approaches to planar and volumetric scintigraphic imaging of small volume targets with high spatial resolution and sensitivity | |
Us | Reduction of Limited Angle Artifacts in Medical Tomography via Image Reconstruction | |
CN118177842A (zh) | 正电子发射断层显像(pet)系统中的列表模式图像重建的系统和方法 | |
CN116917773A (zh) | 辐射检测系统和方法 | |
Zhang et al. | Complete tomographic data for Compton imaging probes | |
Kulkarni | Image quality in MAP SPECT reconstructions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110706 Termination date: 20201104 |