CN1903128A - 正电子ct设备的3d图像重构方法和正电子ct设备 - Google Patents
正电子ct设备的3d图像重构方法和正电子ct设备 Download PDFInfo
- Publication number
- CN1903128A CN1903128A CNA2006101078348A CN200610107834A CN1903128A CN 1903128 A CN1903128 A CN 1903128A CN A2006101078348 A CNA2006101078348 A CN A2006101078348A CN 200610107834 A CN200610107834 A CN 200610107834A CN 1903128 A CN1903128 A CN 1903128A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- munder
- data
- msubsup
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 230000008569 process Effects 0.000 claims abstract description 34
- 238000007792 addition Methods 0.000 claims description 24
- 230000014509 gene expression Effects 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000013179 statistical model Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 abstract description 11
- 230000006866 deterioration Effects 0.000 abstract description 8
- 238000006243 chemical reaction Methods 0.000 abstract description 6
- 238000004904 shortening Methods 0.000 abstract description 5
- 238000010348 incorporation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 14
- 230000005855 radiation Effects 0.000 description 14
- 230000009467 reduction Effects 0.000 description 11
- 238000013500 data storage Methods 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/1611—Applications in the field of nuclear medicine, e.g. in vivo counting using both transmission and emission sources sequentially
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
在降低3D数据获取处理中的存储数据量并且缩短从检查开始到成像结束的时间的同时,由于分辨率降低或者S/N比下降会引起图像质量恶化。在3D数据获取处理的同时,本发明的方法或设备执行正弦图加法、读取已求和的正弦图的子集以及图像重构。因此,减少了存储数据量,并且缩短了从检查开始到成像结束的时间。同时,因为3D迭代重构并不伴随着从3D数据到2D数据的转换,可以避免由于从3D数据转换到2D数据所发生的误差而引起的分辨率降低。3D迭代重构可以直接并入如衰减校正处理等处理。因此,还可以避免不直接并入这种处理而导致的S/N比下降。
Description
技术领域
本发明涉及正电子CT设备(下面适当地称为“PET设备”)的3D图像重构方法,并且涉及一种正电子CT设备,用于通过在沿多环检测器的轴向彼此相对地移动具有沿轴向方向排列的环式检测器的多环检测器和其上放置有被检查对象的床的同时,检测每一个移动位置中的环形对的重合计数数据,来获取环形对的每一个组合的正弦图(sinograms),并且根据这些正弦图来执行RI分布图像重构的3D图像重构。更具体地,本发明涉及一种改善RI分布图像的质量的技术。
背景技术
通常,根据通过3D数据获取而获取的正弦图来执行RI分布图像的图像重构处理的PET设备沿多环检测器的轴向,例如以环式检测器的间隔,彼此相对间断地移动具有沿轴向方向排列的环式检测器的多环检测器和其上放置有被检查对象的床。设备检测每一个移动位置中环形对的重合计数数据,用以收集环形对的每一个组合的正弦图,并且根据这些正弦图,来执行RI分布图像的图像重构处理。
传统的PET设备通常在通过3D数据获取而获取所有正弦图之后,执行用于成像的3D迭代重构。然而,通过3D数据获取而获取的正弦图组成了大量数据。此外,在完成正弦图获取之后开始伴随3D投影和反向投影处理的3D图像重构。存储的数据量必然会非常大,并且从检查开始到成像结束要花费较长时间。
另一种传统PET设备在将正弦图与在紧接前一位置获得的正弦图相加作为从对象的同一位置获取的每一个环形对的正弦图的同时,获取3D数据,并且与3D数据获取处理同时,顺序地读取已求和的正弦图。在将3D数据转换为2D数据之后,在3D数据获取处理的同时,执行2D迭代重构(参见日本未审专利申请No.2004-61113(4-8页,图1-4))。
发明人所提出的该现有设备在3D数据重构处理的同时,执行正弦图的加法和读取以及图像重构处理。因此,该设备可以减少存储的数据量并且缩短从检查开始到成像结束所花费的时间。
然而,由于分辨率降低和S/N比下降,以上PET设备易于引起图像质量恶化。因为该传统正电子设备将3D数据转换为2D数据,所以难以避免由于从3D数据转换到2D数据而发生的误差所引起的分辨率降低。此外,因为从2D数据中重构图像,该设备不能在图像重构时直接使用3D数据物理模型,例如检测器响应、衰减校正和统计属性。该设备在改进重构图像的S/N比上受限。
发明内容
考虑到上述现有问题而做出本发明,并且本发明的目的是提供一种正电子CT设备的3D图像重构方法,并且提供一种正电子CT设备,能够在减少在3D数据获取处理时存储的数据量、并且缩短从检查开始到成像结束的时间的同时,避免由于分辨率降低或S/N比下降而引起的图像质量恶化。
根据本发明,由正电子CT设备的3D重构方法实现以上目的,所述方法包括:3D数据获取步骤,用于在沿多环检测器的轴向,彼此相对地移动具有沿轴向方向排列的环式检测器的多环检测器和其上放置有被对象的床的同时,收集每一个移动位置中的环形对的重合计数数据,来获取环形对的每一个组合的正弦图;以及3D图像重构步骤,用于根据在3D数据获取步骤中所获取的正弦图,来执行RI分布图像重构的3D图像重构;其中,将在移动多环检测器时在每一个位置处获得的正弦图与紧接前一位置处获得的正弦图相比较,并且对从对象相同部位所获取的每一个环形对的正弦图求和,来执行3D数据获取步骤;并且,在3D数据获取步骤的同时,执行3D图像重构步骤,用以顺序地读取已求和的正弦图的子集,并且在每次读取子集时,根据读取的子集,来执行3D迭代重构。
利用根据本发明的方法,在3D数据获取时,通过在沿多环检测器的轴向彼此相对地移动多环检测器和其上放置有被检查对象的床的同时,收集在每一个移动位置中的环形对的重合计数数据,来获得环形对的每一个组合的RI分布图像重构的正弦图。在3D数据获取中,在将移动多环检测器时在每一个位置处所获得的正弦图与在紧接前一位置处所获得的正弦图相比较、并且对从对象相同部位所获取的每一个环形对的正弦图求和的同时,获取沿轴向方向具有较小灵敏度变化的正弦图。对于在多环检测器和支撑被检查对象的床的相对移动时在每一个位置中获得的正弦图,在紧接前一位置处所获得的正弦图与同每一个环形相邻的环形对的正弦图重叠(即从对象相同部位获取环形对的这些正弦图)。通过相加重叠的正弦图,获取沿轴向方向具有较少灵敏度变化的正弦图。
此外,根据本发明的方法在3D数据获取的同时,仅顺序地读取已求和的正弦图的子集。在每一次读取正弦图的子集时,根据读取的子集来执行3D迭代重构,以便完成最终的RI分布图像。
因此,根据本发明的方法在3D数据获取处理的同时执行正弦图的加法和子集的读取,并且在每次读取子集时执行3D图像重构处理。因此,减少了在数据获取处理时存储的数据量,并且缩短了从检查开始到成像结束的时间。
同时,上述方法使用3D数据并且并不伴随着从3D数据到2D数据的转换,来执行3D迭代重构。这种图像重构可以避免由于从3D数据转换到2D数据而发生误差所引起的分辨率降低。3D迭代重构可以直接并入如衰减校正处理等处理。因此,还可以避免由于不能直接并入这些处理而导致的S/N比下降。
如上所述,根据本发明的方法能够在减少3D数据获取处理时的数据存储量、并且缩短从检查开始到成像结束的时间的同时,避免由于分辨率降低或者S/N比下降而引起的图像质量恶化。
在本发明的另一个方案中,一种正电子CT设备包括:多环检测器,具有沿轴向方向排列的环式检测器;床,用于支撑被检查对象;移动设备,用于沿多环检测器的轴向彼此相对地移动多环检测器和床;3D数据获取设备,用于通过收集多环检测器和床的每一个相对移动位置中的环形对的重合计数数据,来获得环形对的每一个组合的正弦图,并且在将移动多环检测器时在每一个位置处获得的正弦图与在紧接前一位置处获得的正弦图相比较、并且对从对象相同部位所获取的每一个环形对的正弦图求和的同时,获取数据;读取设备,用于在3D数据获取处理的同时顺序地读取已求和的正弦图的子集;以及3D图像重构设备,用于在每一次读取设备读取子集时,根据读取的子集,来执行3D迭代重构。
利用根据本发明的设备,在3D数据获取时,在移动设备沿多环检测器的轴向彼此相对地移动多环检测器和其上放置有被检查对象的床的同时,3D数据获取设备通过收集每一个移动位置中的环形对的重合计数数据,获取环形对的每一个组合的RI分布图像重构的正弦图。在3D数据获取中,在将移动多环检测时在每一个位置处所获得的正弦图与在紧接前一位置处所获得的正弦图相比较、并且对从对象相同部位获取的每一个环形对的正弦图求和的同时,获取正弦图。对于在多环检测器和支撑被检查对象的床的相对移动时在每一个位置所获得的正弦图,在紧接前一位置处所获得的正弦图与同每一个环相邻的环形对的正弦图重叠(即,从对象相同部位获取环形对的这些正弦图)。通过相加重叠的正弦图,获取沿轴向方向具有较小灵敏度变化的正弦图。
此外,在根据本发明的设备中,与3D数据获取的同时,读取设备顺序地读取已求和的正弦图的子集。此外,在3D数据获取的同时,3D图像重构设备在每一次读取正弦图子集时,根据读取的字节来执行3D迭代重构,以便完成最终的RI分布图像。
因此,根据本发明的设备执行根据本发明的正电子CT设备的3D图像重构方法,用以针对RI分布图像重构,重构3D图像。
因此,根据本发明的正电子CT设备能够在减少3D数据获取处理时的数据存储量并且缩短从检查开始到成像结束的时间的同时,避免由于分辨率降低或者S/N比下降而引起的图像质量恶化。
在根据本发明的方法和设备中,优选地,如下执行3D迭代重构。
根据统计模型,其中xj是每一个重构图像的像素值,yi是获取的数据,并且假设获取的数据yi具有泊松分布,则重构图像的对数似然值L(x)由下面的运算表达式表示:
其中,
aij是用于匹配重构图像j和获取数据i的矩阵,并且const是常数;以及
通过确定关于x的最大似然值估计L(x)来重构图像,其中根据以下运算表达式来更新解,即,将获取的数据划分为正弦图的子集,并且更新每一个子集的解:
其中,
xj (k,q+1)是通过第k次迭代的子集Sq+1而获取的像素值,即最新的更新图像的像素值,
xj (k,q)是通过第k次迭代的子集Sq而获取的像素值,即紧接更新之前的图像的像素值,
λk是松弛因子,以及
在根据本发明的方法和设备中,优选地,在3D迭代重构之前或期间执行多次加法差校正处理,用于消除由于3D数据获取中的正弦图加法次数的差别而引起的数据不均匀。按照这种方式,可以避免由于正弦图加法次数的差别而引起的图像质量恶化。
附图说明
为了演示本发明,在附图中示出了多种当前优选的形式,然而,应该理解到,本发明不局限于所示的精确布置和手段。
图1是示出了根据本发明的PET设备的概要的方框图;
图2是示出了PET设备中被检查对象和多环检测器之间移动的位置关系的示意图;
图3A至3C是用于演示正弦图加法的示意图;
图4是示出了PET设备中衰减校正数据获取系统的构造的示意图;
图5是示出了PET设备中与组成被检查对象的给定移动位置的子集的正弦图相对应的环形对的示意图;
图6是示出了在PET设备的每一步骤中读取的正弦图子集的示意图;
图7是示出了在PET设备读取的正弦图子集中由于加法次数的差别而引起的数据不均匀的状态的示意图;
图8是示出了与PET设备读取的正弦图子集相对应的衰减校正数据的示意图;以及
图9是示出了从PET设备所获得的3D重构图像的示意图。
具体实施方式
下面参考附图来详细描述本发明的优选实施例。
现在参考附图来描述根据本发明的方法和设备。图1是示出了根据本发明的PET设备的概要的方框图。图2是示出了PET设备中被检查对象和多环检测器之间移动的位置关系的示意图。
根据本发明的PET设备包括:具有中心开口的托台1和以环形形式布置的辐射检测器形式的环形检测器3。多环检测器5具有例如环式辐射检测器3,所述环式辐射检测器组成沿托台1的开口的轴向布置的六个环(#1至#6)。
在托台1的前面布置的是可相对于托台1的开口前后(出入)移动的床7、以及用于沿图1中的两点式点划线所指示的方向前后移动床7的床移动设备9。由移动控制单元11控制床7的移动,使之沿与托台1的开口的中心对齐的身体轴MP移动,使得床7可以被移进或移出托台1的开口。
如上所述,多环检测器5具有组成第一至第六环#1至#6的环式检测器。在布置在第一至第六环#1至#6的每一个的周围的辐射检测器中,并且在组成第一至第六环#1至#6的辐射检测器中,执行重合计数。由重合数据获取单元13来执行该重合计数。
具体地,如图2中与环式检测器3相连的线所示,在环形对的各种组合中,例如在组成第一环#1的检测器中、在组成第一环#1和第二环#2的检测器中、在组成第一环#1和第三环#3的检测器中、在组成第一环#1和第四环#4的检测器中、在组成第一环#1和第五环#5的检测器中、以及在组成第一环#1和第六环#6的检测器中,执行重合计数。因为还在组成不同环的环式检测器3中执行重合计数,所以对于环形对的每一种组合,都存在重合计数线。
现在,再参考图3A至3C来进一步描述实施例。图3A至3C是用于演示在移位正弦图矩阵的同时执行的正弦图加法的示意图。图3A示出了在连续地获取数据时的相对移动。图3B示出了在移位所获取的矩阵的同时在各个移动位置处所获得的间接正弦图的加法。图3C是示出了此时所顺序地读取的一个矩阵的一部分(子集)的示意图。
如图3A所示,此处假设例如在相对于多环式检测器5以检测器间隔σ间断地移动床7的同时,在每一个移动位置处获取数据。在这种情况下,在图3B中,由重合数据获取单元13所收集的床7的每一个位置中的重合数据被表达为r1×r2的矩阵。组成矩阵的每一个小方块表示每一个环形对的正弦图。由到中轴(径向矢量s)的距离和使同时检测辐射的两个检测器相连接的线(重合线)的周向角度(方位角φ)来表达正弦图。将在每一个移动位置处所获得的数据作为正弦图传送到获取数据存储器15并且存储在其中。将按照这种方式所传送的正弦图与在紧接前一位置处所获取的正弦图相比较,并且对从患者M的相同部位所获取的每一个环形对的正弦图求和。因此,重合数据获取单元13和获取数据存储器15组成根据本发明的设备中的3D数据获取设备。
接下来,使用公认的运算表达式(变换),来具体地描述以上重合数据获取单元13、获取数据存储器15和处理单元17的协作处理的细节。
利用具有组成N个环的环式检测器3的PET设备,在床7的每一个位置k处获得N2个正弦图qk(s,φ,r,Δr)。此时,s和φ是正弦图的径向矢量方向和方位角,r=r1+r2是轴向片指数,并且Δr=r1-r2是环形差。r1和r2是经历了重合计数的环所属的环序号。为了便于理解,下面假设以与环式检测器3之间的间隔(宽度φ)相同的间隔在连续获取中采样数据。在这种情况下,获得图3B所示的数据空间。从图3B可见,重叠地获取相邻正弦图,因此在实时相加正弦图的同时,形成了一组新的正弦图ph(s,φ,Δr)。
此处,h=0,....,N+n-2是在以n步骤执行连续获取时所获得的虚拟环指数。此时,每一个正弦图由以下运算表达式(1)表达:
其中,Δr=(r1-r2)=0,±1,±2,...,±Δrmax,以及kmin=max{0,h-(N-1-|Δr|)}。kmax=min{h,n-1},K=kmax-kmin+1,并且Δrmax是用于重构的环差的最大数。
通过加法的次数K来归一化新的正弦图组ph(s,φ,Δr)。除了在相对端(h<N-1或h>N-1)处的(N-1)数据之外,加法次数K反比于环差Δr。在图3B中,小方块表示正弦图qk,其中根据数据加法次数K来以线变化施加阴影。对于每一个步骤都执行了加法的正弦图子集是如图3C所示的中间拐弯的字母L的形式的划分组。
利用根据本实施例的PET设备,在3D数据获取处理同时,在完成加法时仅读取正弦图子集。即,在每一个步骤中读取完成了加法的正弦图子集。在3D图像重构处理之前,设备执行加法差校正处理多次,以消除读取的子集中由于正弦图加法次数K的差别而引起的数据不均匀。在该实施例中,运算表达式(1)中乘法因子(1/K)与加法差校正处理的次数相对应。
另一方面,在本实施例的PET设备的情况下,如图4所示,多环检测器5包含用于获取衰减校正数据(传输数据)的环式辐射检测器20和设置在其开口端相邻处的环式外部辐射源21。利用环式辐射检测器20和环式外部辐射源21,可以由多环检测器5在3D数据获取之前从患者M的相同部位获取衰减校正数据。将衰减校正数据获取单元(未示出)根据从环式辐射检测器20输出的检测信号而获取的衰减校正数据以与正弦图子集相对应的形式存储在衰减校正数据存储器(未示出)中。
设置片隔板22,用于避免来自环式外部辐射源21的辐射向多环检测器5泄漏。
可以由可周向旋转的点外部辐射源来代替环式外部辐射源21。
PET设备不必具有获取衰减校正数据的功能。PET设备可以使用从例如X射线CT设备的独立设备所接收的衰减校正数据。
另一方面,在3D数据获取处理的同时,处理单元17根据通过3D数据获取所获取的正弦图,来执行RI分布图像重构的3D图像重构。具体地,在三维数据获取处理的同时,在每一次读取正弦图子集时,根据读取的子集来执行3D迭代重构。因此,处理单元17与用于执行RI分布图像重构的3D图像重构的3D图像重构设备相对应。
再参考图5和6。图5是示出了与给定移动位置的子集的正弦图相对应的环形对的示意图。图6是示出了在每一步骤中读取的正弦图子集的示意图。如图6所示,正弦图子集包括属于中间拐弯的字母L区域的多个正弦图,与N2正弦图矩阵的上侧部分和左手侧部分相对应。即,子集包括与图5中的线所连接的环形对相对应的多个正弦图。
子集固有地在数据中具有由于正弦图加法次数K的差别(如图6中线变化的阴影所示)而引起的不均匀性。然而,通过用于消除数据中由于正弦图加法次数K的差别而引起的不均匀的多次差校正处理,可以消除数据非均匀性,如图7所示。因此,在本实施例的PET设备中,因为在消除了由于正弦图加法次数K的差别而引起的数据中的不均匀性之后才执行3D迭代重构,所以避免了由于正弦图加法次数K的差别而引起的图像质量恶化。
由处理单元17所执行的3D迭代重构是根据下面给出的运算表达式而执行的所谓3D(三维)-DRAMA(动态行作用最大似然值算法(Dynamic Row Action Maximum likelihood Algorithm))方法的处理(例如,参见Atsusi Fukano和Hiroyuki Kudo的Conference Recordof 2004 IEEE Nuclear Science Symposium and Medical ImagingConference,Paper No.M4-1)。
根据统计模型,其中xj是用于组成3D重构图像fh(x,y,z)的每一个发射图像(重构图像)的像素值,yi是获取的数据(投影数据),并且假设获取的数据yi具有泊松分布,重构图像的对数似然值L(x)由以下运算表达式(2)表达:
其中,
aij是用于匹配重构图像j和获取数据i的矩阵,并且const是常数。
通过确定关于x的L(x)的最大似然值估计,可以重构图像。根据以下运算表达式(3)来更新3D-DRAMA方法的解:
其中
xj (k,q+1)是第k次迭代的子集Sq+1所获取的像素值,即最新更新图像的像素值,
xj (k,q)是第k次迭代的子集Sq所获取像素值,即紧接更新之前的图像的像素值,
λk是松弛因子,以及
即,运算表达式(3)将获取数据划分为正弦图子集,并且通过更新每一个子集的解,以使运算速度获得改进。
利用3D-DRAMA方法,如图8所示,以与正弦图子集相对应的矩阵形式获取并存储衰减校正数据,并且衰减校正数据是矩阵aij的确定因子之一。可以将该数据反映到使重构图像j和获取数据i相匹配的矩阵aij的值。结果,直接将用于抑制S/N比下降的例如衰减校正处理的处理并入3D图像重构处理中。还可以将抑制S/N比下降的检测器响应校正、统计波动校正等反映到矩阵aij的值来直接并入3D图像重构处理中。此外,可以将加法次数K反映到矩阵aij,用以校正在重构处理期间正弦图加法次数K的差别。
在本实施例的设备的3D图像重构处理中,通过在床7的每一个移动位置(每一个步骤)处获取虚拟环形指数的3D重构图像fh(x,y,z),完成RI分布图像,如图9所示。在相邻步骤之间,重构图像重叠。在3D图像重构处理中,前一个重叠图像被用作下一个重叠图像的初始图像。
通过利用重构处理的处理来更新和显示图像,可以检查期间在显示器19上实时持续地监控图像。
利用本实施例的设备所使用的3D-DRAMA方法,可以在按照正弦图的方位角φ来进一步划分每一个正弦图子集的数据之后,或者在按照正弦图的方位角φ和径向矢量方向s来进一步划分每一个正弦图子集的数据之后,执行重构处理。
在本实施例中,在3D数据获取处理的同时,该设备求和正弦图,读取子集,并且根据读取的子集来执行RI分布图像重构的图像重构处理。因此,该设备可以减少在3D数据获取处理期间的数据存储量,并且缩短从检查开始到成像结束的时间。
在本实施例中,该设备采用用于使用3D数据来执行图像重构、但并不伴随着从3D数据到2D数据的转换的3D迭代重构。因此,可以避免由于从3D数据到2D数据的转换所发生的误差而引起的分辨率降低。3D迭代重构可以直接并入如衰减校正处理等用于抑制S/N比下降的处理。与不直接并入用于抑制S/N比下降的处理的技术相比,这对于避免S/N比下降更为有效。
因此,根据本实施例中的方法和设备,可以在降低3D数据获取处理期间的数据存储量、并且缩短从检查开始到成像结束的时间的同时,避免由于分辨率降低和S/N比下降而引起的图像质量恶化。
本发明不局限于上述实施例,而可以如下进行修改:
(1)在上述实施例中,根据运算表达式(1)来执行3D数据获取。可以根据除了运算表达式(1)之外的其他表达式来执行本发明中的3D数据获取。
(2)在上述实施例中,根据运算表达式(2)和(3)来执行3D迭代重构。可以根据除了表达式(2)和(3)之外的其他表达式来执行本发明的3D迭代重构。
例如,作为表达式(3)的修改,引用以下三个表达式。
考虑以φ(x)=L(x)-βR(x)的最大化值代替表达式(2)中的最大化L(x),则可以通过以下修改的表达式来更新图像:
其中,R(x)是包括图像的先验信息的势函数,β是用于调节其度数的正则化参数,xj MLEM<k+1>是从表达式(3)、(3)’、(3)”或(3)推导出的像素值。
在不脱离本发明的精神或本质属性下,可以在其它具体形式中实现本发明,因此,除了上述说明书之外,还可以参考所附权利要求,作为本发明的范围的指示。
Claims (6)
1.一种正电子CT设备的3D图像重构方法,包括:
3D数据获取步骤,用于在沿多环检测器的轴向,彼此相对地移动具有沿轴向方向排列的环式检测器的多环检测器和其上放置有被检查对象的床的同时,收集每一个移动位置中的环形对的重合计数数据,来获取环形对的每一个组合的正弦图;以及
3D图像重构步骤,用于根据在所述3D数据获取步骤中所获取的正弦图,来执行RI分布图像重构的3D图像重构;
其中,将在移动多环检测器时在每一个位置处获得的正弦图与紧接前一位置处获得的正弦图相比较,并且对从对象相同部位所获取的每一个环形对的正弦图求和,来执行所述3D数据获取步骤;以及
在所述3D数据获取步骤的同时,执行所述3D图像重构步骤,以顺序地读取已求和的正弦图的子集,并且在每次读取子集时,根据读取的子集,来执行3D迭代重构。
2.根据权利要求1所述的正电子CT设备的3D图像重构方法,其中:
在执行3D迭代重构中,
根据统计模型,其中xj是每一个重构图像的像素值,yi是获取的数据,并且假设获取的数据yi具有泊松分布,则重构图像的对数似然值L(x)由下面的运算表达式表示:
其中,
aij是用于匹配重构图像j和获取数据i的矩阵,并且const是常数;以及
通过确定关于x的最大似然值估计L(x)来重构图像,其中根据以下运算表达式来更新解,即,将获取的数据划分为正弦图的子集,并且更新每一个子集的解:
其中,
xj (k,q+1)是通过第k次迭代的子集Sq+1而获取的像素值,即最新的更新图像的像素值,
xj (k,q)是通过第k次迭代的子集Sq而获取的像素值,即紧接更新之前的图像的像素值,
λk是松弛因子,以及
3.根据权利要求1所述的正电子CT设备的3D图像重构方法,其中,在3D迭代重构之前或期间执行多次加法差校正处理,用于消除由于3D数据获取中的正弦图加法次数的差别而引起的数据不均匀。
4.一种正电子CT设备,包括:
多环检测器,具有沿轴向方向排列的环式检测器;
床,用于支撑被检查对象;
移动设备,用于沿所述多环检测器的轴向彼此相对地移动所述多环检测器和所述床;
3D数据获取设备,用于通过收集所述多环检测器和所述床的每一个相对移动位置中的环形对的重合计数数据,来获得环形对的每一个组合的正弦图,并且在将移动多环检测器时在每一个位置处获得的正弦图与在紧接前一位置处获得的正弦图相比较、并且对从对象相同部位所获取的每一个环形对的正弦图求和的同时,获取数据;
读取设备,用于在所述3D数据获取步骤的同时顺序地读取已求和的正弦图的子集;以及
3D图像重构设备,用于在每一次所述读取设备读取子集时,根据读取的子集,来执行3D迭代重构。
5.根据权利要求4所述的正电子CT设备,其中,所述3D图像重构设备执行3D迭代重构,使得
根据统计模型,其中xj是每一个重构图像的像素值,yi是获取的数据,并且假设获取的数据yi具有泊松分布,则重构图像的对数似然值L(x)由下面的运算表达式表示:
其中,
aij是用于匹配重构图像j和获取数据i的矩阵,并且const是常数;以及
通过确定关于x的最大似然值估计L(x)来重构图像,其中根据以下运算表达式来更新解,即,将获取的数据划分为正弦图的子集,并且更新每一个子集的解:
其中,
xj (k,q+1)是通过第k次迭代的子集Sq+1而获取的像素值,即最新的更新图像的像素值,
xj (k,q)是通过第k次迭代的子集Sq而获取的像素值,即紧接更新之前的图像的像素值,
λk是松弛因子,以及
6.根据权利要求4所述的正电子CT设备,其中,在3D迭代重构之前或期间,所述3D图像重构设备执行多次加法差校正处理,用于消除由于3D数据获取中的正弦图加法次数的差别而引起的数据不均匀。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005214019 | 2005-07-25 | ||
JP2005214019A JP4371083B2 (ja) | 2005-07-25 | 2005-07-25 | ポジトロンct装置の3次元画像再構成処理方法、及び、その装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1903128A true CN1903128A (zh) | 2007-01-31 |
CN100484477C CN100484477C (zh) | 2009-05-06 |
Family
ID=37650552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101078348A Expired - Fee Related CN100484477C (zh) | 2005-07-25 | 2006-07-24 | 正电子ct设备的3d图像重构方法和正电子ct设备 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7391027B2 (zh) |
JP (1) | JP4371083B2 (zh) |
CN (1) | CN100484477C (zh) |
DE (1) | DE102006034356A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101953693A (zh) * | 2009-07-16 | 2011-01-26 | 西门子公司 | 完整或部分正电子发射断层造影衰减图的基于模型的估计 |
WO2013097390A1 (zh) * | 2011-12-30 | 2013-07-04 | 沈阳东软派斯通医疗系统有限公司 | Pet系统中图像的衰减校正方法及装置 |
CN105792754A (zh) * | 2013-07-18 | 2016-07-20 | 通用电气公司 | 用于轴向分段正电子发射断层造影数据的方法和系统 |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10307331B4 (de) * | 2003-02-17 | 2009-03-05 | BAM Bundesanstalt für Materialforschung und -prüfung | Bildgebendes Verfahren zur rechnergestützten Auswertung computer-tomographischer Messungen durch direkte iterative Rekonstruktion |
US7750304B2 (en) * | 2004-12-15 | 2010-07-06 | Koninklijke Philips Electronics N.V. | Concurrent reconstruction using multiple bed frames or continuous bed motion |
US7671331B2 (en) * | 2006-07-17 | 2010-03-02 | General Electric Company | Apparatus and methods for processing imaging data from multiple detectors |
US7890282B2 (en) * | 2007-09-17 | 2011-02-15 | Siemens Medical Solutions Usa, Inc. | Estimation of crystal efficiency with axially compressed sinogram |
WO2009054070A1 (ja) * | 2007-10-26 | 2009-04-30 | Shimadzu Corporation | 放射線検出器 |
EP2267483B1 (en) * | 2008-04-14 | 2017-06-14 | National Institutes for Quantum and Radiological Science and Technology | Pet device and method for determining arrangement of detectors |
JP5405187B2 (ja) * | 2009-05-08 | 2014-02-05 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
US8379948B2 (en) | 2010-12-21 | 2013-02-19 | General Electric Company | Methods and systems for fast iterative reconstruction using separable system models |
ES2651093T3 (es) * | 2011-07-13 | 2018-01-24 | Scivis Wissenschaftliche Bildverarbeitung Gmbh | Reconstrucción de imágenes de transmisión y obtención de imágenes utilizando datos de detector de Poisson |
US10575802B2 (en) | 2013-09-27 | 2020-03-03 | General Electric Company | Systems and methods for controlling motion of detectors having moving detector heads |
US10492745B2 (en) | 2013-09-27 | 2019-12-03 | Ge Medical Systems Israel, Ltd | Systems and methods for controlling motion of detectors having moving detector heads |
US9402595B2 (en) | 2014-01-29 | 2016-08-02 | General Electric Company | Systems and methods for positioning detector heads |
US9349495B2 (en) | 2014-04-22 | 2016-05-24 | General Electric Company | Systems and methods for improved collimation sensitivity |
US9763631B2 (en) | 2014-09-17 | 2017-09-19 | General Electric Company | Systems and methods for imaging plural axial locations |
US10143437B2 (en) | 2015-06-30 | 2018-12-04 | General Electric Company | Systems and methods for dynamic scanning with multi-head camera |
US11534128B2 (en) | 2017-05-23 | 2022-12-27 | General Electric Company | Systems and methods for image quality enhancement for multi-head camera |
US10213174B1 (en) | 2018-01-05 | 2019-02-26 | General Electric Company | Nuclear medicine imaging systems and methods having multiple detector assemblies |
US10481285B1 (en) | 2018-08-13 | 2019-11-19 | General Electric Company | Systems and methods for determination of depth of interaction |
US10976452B2 (en) | 2018-08-13 | 2021-04-13 | General Electric Medical Systems Israel, Ltd. (Il) | Systems and methods for improved medical imaging |
US11151758B2 (en) | 2018-11-07 | 2021-10-19 | General Electric Company | Systems and methods for image quality enhancement for multi-head camera |
CN109978831B (zh) * | 2019-02-28 | 2021-02-19 | 成都真实维度科技有限公司 | 一种关节非骨质建模分区体积比的正常参照值的获取方法 |
US11320545B2 (en) | 2020-07-07 | 2022-05-03 | GE Precision Healthcare LLC | Systems and methods for improved medical imaging |
US11092701B1 (en) | 2020-07-07 | 2021-08-17 | GE Precision Healthcare LLC | Systems and methods for improved medical imaging |
Family Cites Families (7)
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 |
US5990482A (en) * | 1997-05-30 | 1999-11-23 | Adac Laboratories | Method and apparatus for performing correction of transmission self-contamination in a medical imaging system |
US6332035B1 (en) * | 1999-06-23 | 2001-12-18 | The Board Of Trustees Of The University Of Illinois | Fast hierarchical reprojection algorithms for 3D radon transforms |
US6590213B2 (en) * | 2001-09-07 | 2003-07-08 | Ge Medical Systems Global Technology Company, Llc | Method and system for estimating scatter in a pet scanner |
US6915004B2 (en) * | 2002-02-28 | 2005-07-05 | Cti Pet Systems, Inc. | Continuous tomography bed motion data processing apparatus and method |
JP4238534B2 (ja) | 2002-07-24 | 2009-03-18 | 株式会社島津製作所 | ポジトロンct装置の3次元画像再構成方法及びその装置 |
US7173248B2 (en) * | 2004-10-20 | 2007-02-06 | General Electric Company | Methods and systems for positron emission tomography data correction |
-
2005
- 2005-07-25 JP JP2005214019A patent/JP4371083B2/ja not_active Expired - Fee Related
-
2006
- 2006-07-18 US US11/488,020 patent/US7391027B2/en not_active Expired - Fee Related
- 2006-07-24 CN CNB2006101078348A patent/CN100484477C/zh not_active Expired - Fee Related
- 2006-07-25 DE DE102006034356A patent/DE102006034356A1/de not_active Withdrawn
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101953693A (zh) * | 2009-07-16 | 2011-01-26 | 西门子公司 | 完整或部分正电子发射断层造影衰减图的基于模型的估计 |
CN101953693B (zh) * | 2009-07-16 | 2015-02-18 | 西门子公司 | 完整或部分正电子发射断层造影衰减图的基于模型的估计 |
WO2013097390A1 (zh) * | 2011-12-30 | 2013-07-04 | 沈阳东软派斯通医疗系统有限公司 | Pet系统中图像的衰减校正方法及装置 |
US9311723B2 (en) | 2011-12-30 | 2016-04-12 | Neusoft Positron Medical Systems Co., Ltd. Shenyang | Attenuation correction method and device for image in PET system |
CN105792754A (zh) * | 2013-07-18 | 2016-07-20 | 通用电气公司 | 用于轴向分段正电子发射断层造影数据的方法和系统 |
CN105792754B (zh) * | 2013-07-18 | 2020-01-21 | 通用电气公司 | 用于轴向分段正电子发射断层造影数据的方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
US7391027B2 (en) | 2008-06-24 |
JP4371083B2 (ja) | 2009-11-25 |
JP2007033118A (ja) | 2007-02-08 |
US20070018108A1 (en) | 2007-01-25 |
DE102006034356A1 (de) | 2007-02-01 |
CN100484477C (zh) | 2009-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1903128A (zh) | 正电子ct设备的3d图像重构方法和正电子ct设备 | |
CN107212899B (zh) | 医学成像方法及医学成像系统 | |
CN1292712C (zh) | 多模态成像方法和装置 | |
US8520974B2 (en) | Image processing method, an apparatus therefor and a tomographic apparatus for removing artifacts from a sectional image | |
JP5142664B2 (ja) | X線コンピュータ断層撮影装置 | |
CN109949411B (zh) | 一种基于三维加权滤波反投影和统计迭代的图像重建方法 | |
US8971599B2 (en) | Tomographic iterative reconstruction | |
CN1217625C (zh) | 用于具有任意节距的多切面ct成像的部分扫描加权 | |
JP5099750B2 (ja) | 断層撮影装置の画像再構成方法、故障診断方法、断層撮影装置、及び、システムマトリクスの管理プログラム | |
US20070274581A1 (en) | Methods and apparatus for BIS correction | |
CN108876730B (zh) | 校正运动伪影的方法、装置及设备和存储介质 | |
CN103054599B (zh) | X射线ct装置及x射线ct装置的动作方法 | |
CN111110260B (zh) | 一种图像重建方法、装置及终端设备 | |
CN1175783C (zh) | 正电子断层扫描中热源高分辨快速图像迭代重建方法 | |
JP6176828B2 (ja) | 画像再構成装置、画像再構成方法およびx線コンピュータ断層撮影装置 | |
Lee et al. | Pragmatic fully 3D image reconstruction for the MiCES mouse imaging PET scanner | |
CN105118039A (zh) | 实现锥束ct图像重建的方法及系统 | |
CN108447553A (zh) | 一种医学成像方法、装置、系统及计算机可读存储介质 | |
CN110264536B (zh) | 一种在平行束超分重建中计算高低分辨率投影关系的方法 | |
JP4238534B2 (ja) | ポジトロンct装置の3次元画像再構成方法及びその装置 | |
JP6615531B2 (ja) | X線コンピュータ断層撮影装置及び医用画像処理装置 | |
Farquhar et al. | An evaluation of exact and approximate 3-D reconstruction algorithms for a high-resolution, small-animal PET scanner | |
CN111080737B (zh) | 图像重建方法、装置及pet扫描系统 | |
US11468607B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
CN108242066B (zh) | Ct图像的空间分辨率增强装置和方法以及ct成像系统 |
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: 20090506 Termination date: 20160724 |