CN103890609A - 偶发同时计数推定方法以及偶发同时计数推定装置 - Google Patents
偶发同时计数推定方法以及偶发同时计数推定装置 Download PDFInfo
- Publication number
- CN103890609A CN103890609A CN201380003466.0A CN201380003466A CN103890609A CN 103890609 A CN103890609 A CN 103890609A CN 201380003466 A CN201380003466 A CN 201380003466A CN 103890609 A CN103890609 A CN 103890609A
- Authority
- CN
- China
- Prior art keywords
- tof
- mentioned
- event
- data
- enumeration data
- 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 37
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 230000015654 memory Effects 0.000 claims description 13
- 238000002600 positron emission tomography Methods 0.000 claims description 6
- 241001269238 Data Species 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000012636 positron electron tomography Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000030984 MIRAGE syndrome Diseases 0.000 description 2
- 208000035217 Ring chromosome 1 syndrome Diseases 0.000 description 2
- 208000032825 Ring chromosome 2 syndrome Diseases 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005251 gamma ray Effects 0.000 description 2
- TVLSRXXIMLFWEO-UHFFFAOYSA-N prochloraz Chemical compound C1=CN=CN1C(=O)N(CCC)CCOC1=C(Cl)C=C(Cl)C=C1Cl TVLSRXXIMLFWEO-UHFFFAOYSA-N 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 101000666896 Homo sapiens V-type immunoglobulin domain-containing suppressor of T-cell activation Proteins 0.000 description 1
- 238000012879 PET imaging Methods 0.000 description 1
- 102100038282 V-type immunoglobulin domain-containing suppressor of T-cell activation Human genes 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- IJJVMEJXYNJXOJ-UHFFFAOYSA-N fluquinconazole Chemical compound C=1C=C(Cl)C=C(Cl)C=1N1C(=O)C2=CC(F)=CC=C2N=C1N1C=NC=N1 IJJVMEJXYNJXOJ-UHFFFAOYSA-N 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Biology (AREA)
- Biotechnology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Artificial Intelligence (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Physiology (AREA)
- Bioethics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Evolutionary Computation (AREA)
- Nuclear Medicine (AREA)
Abstract
根据实施方式,偶发同时计数推定方法取得包含TOF信息的TOF列表模式计数数据,不使用TOF信息地将所取得的TOF列表模式计数数据转换成包含随机计数值的4D原始正弦图计数数据,为了生成4D插补正弦图计数数据,对4D原始正弦图计数数据进行插补,为了除去噪音,对4D插补正弦图计数数据使用低通滤波器,将使用低通滤波器得到的4D插补正弦图计数数据转换成使用滤波器得到的4D原始正弦图计数数据,对使用滤波器得到的4D原始正弦图计数数据有效地使用用于对湮没点进行滤波的TOF蒙片滤波器,从而根据使用滤波器得到的4D原始正弦图计数数据来生成5D TOF原始正弦图计数数据。
Description
技术领域
本说明书所述的实施方式一般涉及由γ射线检测系统取得的飞行时间数据中的随机事件的推定。
背景技术
正电子放射断层摄影(Positron Emission Tomography:PET)是将正电子放射型的放射性医药品导入被检体内的核医学的1个领域。当放射性医药品衰变时,生成正电子。具体而言,通过多个正电子分别在作为正电子湮没事件而熟知的现象中与电子进行反应,生成沿着同时计数线大致向相反方向移动的同时发生对的γ光子。在同时计数时间内检测到的γ光子的对通常作为湮没事件由PET扫描仪来记录。
在飞行时间(time-of-flight:TOF)成像中,还测量检测同时发生对的各γ光子的同时计数间隔内的时间。飞行时间信息表示沿着同时计数线检测到的事件的位置。为了重建或者生成扫描对象的被检体或者被检物的图像,使用多个湮没事件的数据。在此,该图像的重建或者生成一般使用统计学(反复的)或者分析学重建算法来进行。
图1表示在三维检测器中放射的正电子以及测量的同时计数线(line-of-response:LOR)的经轴坐标和轴向坐标。坐标(xe,ye,ze)或者(se,te,ze)表示放射出的正电子的图像坐标。所测定的LOR的投影坐标能够通过由z=(za+zb)/2表示的(s,φ,z,θ)表示,或者也可以包含TOF-LOR用的附加的维度t。
在PET中,根据同时计数窗口的有限宽度而产生偶发同时计数(随机),偶发同时计数用于检测真正的同时计数。当在同时计数窗口内检测两个不存在相关关系的单一事件时,有可能这两个被错误地识别并记录为真正的同时计数事件。随机事件的比例与各检测器上的单一事件的比例以及同时计数窗口的尺寸成比例,由式(1)表示。
Cij=2τrirj (1)
在此,Cij表示连结第i个检测器和第j个检测器的LOR上的偶发同时计数的计数率,ri、rj表示第i个检测器和第j个检测器各自的单一计数事件的比例,τ是同时计数窗口的尺寸。
特别地,在三维PET扫描仪以及高放射能浓度中,偶发同时计数可能形成所记录的瞬发同时计数事件(包含真正的同时计数事件、分散同时计数事件、偶发同时计数事件)的大部分。从而,当没有合适地进行校正时,由于偶发同时计数事件,在重建的图像中可能产生大量的定量误差。
因此,作为用于随机推定的最准确且一般使用的方法,能够列举出使用延迟同时计数窗口。延迟同时计数窗口能够从所记录的各个对事件中,除去2个单一事件的相关关系。延迟同时计数窗口通常使同时计数窗口的尺寸从数十到数百倍地延迟,因此2个所记录的单一事件来自1个湮灭的可能性极低。从而,在本方法中,只记录随机事件。
最近,为了通过由滤波器完全地除去无益于重建视野(field-of-view:FOV)的这些偶发同时计数,从而减少瞬发数据中的随机事件,设计有飞行时间蒙片(TOF蒙片)。在该处理中,通过在开始重建之前除去具有蒙片之外的TOF差的随机事件,从而能够通过列表模式重建以及更短的计算时间,得到更准确的重建图像。
在TOF列表模式重建中,随机事件与TOF信息一起由以下的方法推定。非TOF延迟列表模式数据最初进行重组插补在四维插补正弦图中,平滑化处理之后,反插补在四维原始正弦图中。在最终步骤中,通过将LOR计数均等地扩散到切线方向t的组距(bin),考虑沿着LOR的半径方向的t的范围的差以及倾斜角,从而生成五维TOF原始正弦图。
现有技术文献
专利文献
专利文献1:日本特开平09-184885号公报
发明内容
本发明要解决的问题在于,提供一种能够对随机事件的推定使用TOF蒙片的偶发同时计数推定方法以及偶发同时计数推定装置。
根据一个实施方式,在正电子放射断层摄影列表模式数据中推定偶发同时计数的偶发同时计数推定方法包括:取得包含飞行时间(time-of-flight:TOF)信息的TOF列表模式计数数据;不使用上述TOF信息地将上述取得的TOF列表模式计数数据转换成包含偶发同时计数值的四维原始正弦图计数数据;为了生成四维插补正弦图计数数据而对上述四维原始正弦图计数数据进行插补;为了除去噪音而对上述四维插补正弦图计数数据使用低通滤波器;将使用上述低通滤波器得到的四维插补正弦图计数数据转换成使用滤波器得到的四维原始正弦图计数数据;以及通过处理装置对使用上述滤波器得到的四维原始正弦图计数数据有效地使用用于对湮没点进行滤波的TOF蒙片滤波器,从而根据使用上述滤波器得到的四维正弦图计数数据来生成五维TOF原始正弦图计数数据。
附图说明
如果参照以下的详细的说明与附图建立关联来考虑,则能够易于更完全地理解本发明和其附随的大量的优点。
图1表示PET成像装置的几何学配置的例子。
图2表示新的TOF蒙片的横截面图以及矢状面图。
图3A表示没有600mm的FOV、TOF蒙片,而进行5D随机原始正弦图仿真得到的GATE IEC幻象数据的(rad,t)图。
图3B表示具有600mm的FOV、TOF蒙片以及3σTOF,进行5D随机原始正弦图仿真而得到的GATE IEC幻象数据的(rad,t)图。
图4表示基于一个实施方式的PET中的随机现象的推定方法的流程图。
图5表示基于一个实施方式的PET系统的硬件。
具体实施方式
以下,针对本发明所涉及的实施方式的细节进行说明。在本发明中,如上所示,提供一种能够在偶发同时计数的推定中使用TOF蒙片的偶发同时计数推定方法以及偶发同时计数推定装置。即,如果对上述的以往技术的偶发同时计数推定(随机推定)使用TOF蒙片,则随机的分布伴随着导入椭圆形的TOF蒙片会发生变化,由此在视野(field-of-view:FOV)的中心区域收集的随机事件增加,在FOV的周边部收集的随机事件减少。此时,不能使用没有使用TOF蒙片的上述步骤。其理由是因为:(1)通过随机的平滑化步骤,s维度中的随机的分布发生变化,由此在反插补步骤中会产生推定误差。另外,(2)当使LOR计数均等地向各切线方向t的组距(bin)扩散时,以随机事件数沿着t维度均匀地分布为前提,因此当使用TOF蒙片时,其不适合。这样,随机事件在FOV的中心更多地收集,但是在TOF蒙片的圆形的周边部更少地收集。
因此,根据本发明所涉及的一个实施方式,提供一种在正电子放射断层摄影列表模式数据中推定随机事件的方法,该方法包含:(1)取得包含飞行时间(time-of-flight:TOF)信息的TOF列表模式计数数据的步骤;(2)不使用TOF信息地将所取得的TOF列表模式计数数据转换成包含随机计数值的四维(4D)原始正弦图计数数据的步骤;(3)为了生成4D插补正弦图计数数据,对4D原始正弦图计数数据进行插补的步骤;(4)为了除去噪音,对4D插补正弦图计数数据使用低通滤波器的步骤;(5)将使用低通滤波器得到的4D插补正弦图计数数据转换成使用滤波器得到的4D原始正弦图计数数据的步骤;以及(6)通过由处理装置,对使用滤波器得到的4D原始正弦图计数数据有效地使用用于对湮没点进行滤波的TOF蒙片滤波器,从而根据使用滤波器得到的4D原始正弦图计数数据来生成五维(5D)TOF原始正弦图计数数据的步骤。
在另一个实施方式中,TOF蒙片滤波器沿着同时计数线(line-of-response:LOR),在重建视野外具有规定范围的滤波器区域。
在另一个实施方式中,生成步骤包含根据以下的公式对每个TOF组距(bin)t生成5D TOF原始正弦图计数数据 的步骤。
在此,tcoin是同时计数窗口的mm显示的尺寸,tofbinsz是mm显示的TOF组距(bin)的尺寸,r(rad,phi,slice)是使用滤波器得到的4D原始正弦图随机计数数据,(rad,phi,slice)是原始正弦图空间中的LOR的坐标,θ是LOR的倾斜角。另外,slice是与检测到湮没的γ射线的环(ring)相关的信息,包含{ring Diff,ring Sum}。“ring Diff”表示检测到的环间的差(|ring1-ring2|),“ring Sum”表示检测到的环的中点((ring1+ring2)/2)。
根据另一个实施方式,取得步骤包含对每个同时测量事件取得具有数据排列{xa,za,xb,zb,ea,eb,tof}的TOF列表模式数据的步骤,在此,xa、xb是事件a和事件b各自的横截方向的入射闪烁体的编号,za、zb是事件a和事件b各自的轴方向的入射闪烁体的编号,ea、eb是事件a和事件b各自的能级,tof是事件a和事件b的到达时间的差。
在另一个实施方式中,上述的方法包含使用延迟同时计数窗口法,在4D原始正弦图计数数据内推定随机计数值的步骤。
在此,非TOF列表模式数据一般对每个同时计数事件具有{xa,za,xb,zb,ea,eb}的数据排列。xa、xb是事件a和事件b各自的横截方向的入射闪烁体的编号,za、zb是事件a和事件b各自的轴方向的入射闪烁体的编号,ea、eb是事件a和事件b各自的能级。同样地,TOF列表模式数据具有数据排列{xa,za,xb,zb,ea,eb,tof},在此,tof是事件a与事件b的到达时间的差。
基于使用延迟同时计数窗口的TOF列表模式重建中的随机推定能够使用5D或者4D方法来执行。5D法包含以下的步骤。
(1)将TOF延迟窗口列表模式数据在5D原始正弦图(rad,phi,slice,t)数据中重组。
(2)进行向前插补(forward interpolation),将LOR正弦图(rad,phi,slice,t)数据转换成插补后的均匀的正弦图(s,φ、z、θ、t)数据。
(3)使用基于正弦图的5D低通滤波器,减少随机数据的误差。
(4)通过反插补(back interpolation)取得随机推定,将插补后的均匀的正弦图(s,φ、z、θ、t)数据重新转换成LOR正弦图(rad,phi,slice,t)数据。
同样地,4D法包含以下的步骤。
(1)将TOF延迟窗口列表模式数据在4D原始正弦图(rad,phi,slice)数据中重组。
(2)进行向前插补,将LOR正弦图(rad,phi,slice)数据转换成插补后的均匀的正弦图(s,φ、z、θ)数据。
(3)使用基于正弦图的4D低通滤波器,减少随机数据的误差。
(4)使用反插补,将插补后的均匀的正弦图(s,φ、z、θ)数据重新转换成LOR正弦图(rad,phi,slice)数据。
(5)将反插补的现象沿着t维度均等地分割成(rad,phi,slice)的切线方向t的组距(bin),生成(rad,phi,slice,t)中的随机推定值。
但是,上述的方法难以与在本实施方式中公开的新的飞行时间蒙片(TOF蒙片)一起使用。特别地,如果使用图2所示的椭圆形TOF蒙片,则随机的分布会发生变化。如图2所示,TOF蒙片的形状如以下的式(2)所示与半径方向距离s相关。
随机分布存在缓缓地变化的趋势,且遍布FOV整体大致均匀,因此随机分布在使用TOF蒙片的前后不同。因此,对应于新的TOF蒙片,需要新的随机校正法。
如图2所示,在大半径的位置(s),t在0的组距(bin)中半径最大且随机分布变得更多。从而,在上述方法的最终步骤中,不能使随机计数沿着t维度均匀地分布。即,需要考虑TOF蒙片的形状。
因此,在公开的实施方式中,在使用TOF蒙片的TOF列表模式重建中,使用用于推定各事件的随机的新的方法。代替利用使用TOF蒙片滤波器得到的随机列表模式数据,在发明的实施方式中利用不使用TOF蒙片的随机列表模式数据。以下,有时将其称为“原始随机数据”。
原始随机数据的分布(在使用TOF蒙片滤波器之前)不依存于TOF的大小,因此随机分布沿着切线方向的t维度是均匀的。如图3A所示,如果随机分布除去统计的不可靠性,则在半径方向的每个距离s沿着t维度大致均匀。在使用TOF蒙片并使用了滤波器之后,随机分布也沿着TOF蒙片内的t维度而均匀,但TOF蒙片的范围如图2的s坐标所示的那样沿着LOR的半径方向而变化。
如图3B所示,与周边部相比在半径方向中心部,t不是0的组距(bin)较多,另外如果LOR的倾斜角θ变大,则沿着t维度的随机的分布的范围减少。本发明发现在被TOF蒙片覆盖的区域中,随机的分布与有无使用TOF蒙片无关是完全相同的。
根据该发现,如以下那样公开了使用TOF蒙片的用于TOF列表模式重建的新的随机推定步骤。特别地,在一个实施方式中,如图4那样执行随机推定处理。
在步骤400中,得到包含飞行时间(time-of-flight:TOF)信息的TOF列表模式计数数据。如上所述,TOF列表模式数据具有数据排列{xa,za,xb,zb,ea,eb,tof},在此,tof是事件a与事件b的到达时间的差。
在步骤410中,非TOF延迟列表模式数据(没有TOF蒙片滤波器)最初在(rad,phi,slice)原始正弦图中重组。在该步骤中,忽视数据的TOF部分。
在步骤420中,(rad,phi,slice)原始正弦图数据在4D插补正弦图(s,φ、z、θ)数据中插补。
在步骤430中,为了除去数据内的统计噪音,使用基于正弦图的低通平滑化处理。
在步骤440中,平滑化处理后的正弦图(s,φ、z、θ)数据在4D(rad,phi,slice)原始正弦图数据r中反插补。
在步骤450中,将LOR计数均等地扩散到切线方向t的各个组距(bin),通过考虑沿着LOR的半径方向的t的范围的差以及倾斜角θ,从而生成5D TOF原始正弦图数据特别地,对t的每个组距(bin)推定出的随机计数根据以下的式(3)来确定。
在此,tcoin表示同时计数窗口的mm显示的尺寸,tbinsz表示mm显示的t的组距(bin)的尺寸,另外,trange(s)如式(2)那样进行定义,表示各s中的t的范围。另外,由式(3)生成的5D TOF原始正弦图具有使用图2所示的TOF蒙片的效果。
例如,当使用a±4.2纳秒的同时计数窗口时,成为tcoin=2×4200×0.15=1260mm,另外,当tbinsz=20mm以及θ=0时,各r(rad,phi,slice)以0.0159倍率进行变更,得到这样,使用式(3),实际上不用对使用TOF蒙片得到的数据使用滤波器,而能够预测使用TOF蒙片的效果。
在另一个实施方式中,可以省略步骤410和步骤440的双方或者一方。例如,在步骤400之后,能够执行步骤420,以使得从列表模式数据中直接取得4D插补正弦图数据。
另外,使用TOF蒙片的新的随机推定法不仅能够适用于列表模式重建,而且能够适用于基于正弦图的图像重建。此时,重建引擎能够使用TOF信息,或者也能够不使用TOF信息。当在重建中使用TOF信息中,推定出的随机在5D的TOF插补后的正弦图中使用。在重建中不使用TOF信息的情况下,也能够使用TOF蒙片不接受随机计数的一部分,但重建引擎本身不使用TOF信息,因此使用4D的没有TOF插补的正弦图。
用于上述的TOF列表模式重建的新的随机推定法具有几个优点。特别地,能够根据沿着t维度的随机现象的分布的均匀性,通过使用不使用TOF蒙片滤波器的非TOF延迟列表模式数据来推定瞬发数据中的随机现象。在重建间使用TOF蒙片的优点在此有效。另外,不会由于新的TOF蒙片的形状而增加随机校正的复杂性。如式(3)所示,推定出的随机计数的比例因子值不依存于半径方向距离s。
在图5中,示出能够与本实施方式一起使用的示例的PET硬件结构。在图5中,光电倍增管135以及140被配置于光导130上,闪烁体阵列105被配置于光导130下。第2闪烁体阵列125与闪烁体105对置,光导115和光电倍增管195、110重叠配置。
在图5中,如果从被检体放射γ射线,则γ射线在相互大约180度的相反方向前进。γ射线由闪烁体100以及120同时检测。并且,如果在既定限制时间内由闪烁体100以及120检测γ射线,则确定闪烁现象。这样,γ射线定时检测系统由闪烁体100以及120同时检测γ射线。但是,为了简化,在此只说明通过闪烁体100的γ检测。对于本领域的技术人员而言不言而喻,针对闪烁体100的说明也同样适用于通过闪烁体120的γ射线检测。
各光电倍增管110、135、140、195分别与数据取得装置150连接。数据取得装置150包含构成为对来自光电倍增管的信号进行处理的硬件。数据取得装置150测定γ射线的到达时间。数据取得装置150生成将对于系统时钟(未图示)的识别脉冲的时间进行编码的2个输出(1个用于光电倍增管(photomultiplier tube:PMT)135/140的组合,1个用于PMT110/195的组合)。在飞行时间型PET系统中,数据取得装置150典型的情况是以15~20皮秒的精度生成时间标记。数据取得装置测定各PMT的信号(来自数据取得装置150的4个输出)的大小。
数据取得装置的输出为了进行处理而发送至CPU170。该处理由根据数据取得装置的输出推定能量和位置的处理、以及根据对每个现象输出的时间标记推定到达时间的处理构成,为了改善能量、位置、以及时间的推定值的精度,还可以包含根据事先的校正,使用大量的校正步骤。
另外,CPU170构成为,按照图4所示的流程图以及上述内容执行随机现象的推定方法。
对于本领域的技术人员而言不言而喻,CPU170能够安装为特定用途集成电路(Application Specific Integrated Circuit:ASIC)、现场可编程门阵列(Field Programmable Gate Array:FPGA)、或者复杂可编程逻辑器件(Complex Programmable Logic Device:CPLD)等个别逻辑门。FPGA或者CPLD的安装也可以由VHDL、Verilog或者其他的硬件记述语言进行编码化,另外,也可以将该编码直接存储于FPGA或者CPLD内的电子存储器中,或者作为独立的电子存储器存储于电子存储器。另外,电子存储器也可以是ROM、电气PROM(electrically programmable read only memory:EPROM)、电可擦PROM(electrically erasable programmable read only memory:EEPROM)、闪存器(FLASH memory)等非易失性的。另外,电子存储器也可以是静态或者动态RAM等易失性的。另外,除了FPGA或者CPLD与电子存储器之间的对话之外,为了对电子存储器进行管理,还可以具备微控制器或微处理器等处理装置。
或者,CPU170也可以安装为,保存于上述的电子存储器和硬盘、CD、DVD,U盘(FLASH drive)等已知的存储介质的双方或者一方的任一方中的一连串的计算机可读命令。另外,计算机可读命令作为应用程序、后台程式、或者操作系统的构成要素、或者它们的组合来提供,与美国Intel公司制的Xeon处理器或者美国AMD公司制的Opteron处理器等处理装置、以及Microsoft VISTA、UNIX(注册商标)、Solaris、LINUX(注册商标)、以及Apple MAC-OS等本领域的技术人员所熟知的操作系统联动来执行。
如果由CPU170进行处理,则处理后的信号进行向电子存储装置180的存储和向显示器145的显示的双方或者其一方。对于本领域的技术人员而言不言而喻,电子存储装置180也可以是硬盘驱动器、CD-ROM驱动器、DVD驱动器、U盘、RAM、ROM等本技术领域所熟知的电子存储装置。显示器145也可以安装为液晶显示器(liquid crystaldisplay:LCD)、阴极射线管显示器(cathode ray tube:CRT)、等离子显示器、有机发光二极管(organic light emitting diode:OLED)、发光二级管(LED)等本技术领域所熟知的显示器。这样,在此说明的电子存储装置180以及显示器145的记载只是示例,绝没有限定本实施方式的进步的范围。
另外,图5包含将γ射线检测系统与其他的外部装置和用户的双方或者一方连接的接口175。例如,接口175也可以是通用串行总线(Universal Serial Bus:USB)接口、个人计算机存储卡国际协会(Personal Computer Memory Card International Association:PCMCIA)接口、以太网(注册商标)接口等本技术领域所熟知的接口。另外,接口175可以是有线的或者也可以是无线的,也可以包含与键盘和鼠标的双方或者一方等的用于与用户进行对话的人机接口。
虽然说明了本发明的几个实施方式,但这些实施方式是作为例子而提示的,并不意图限定本发明的范围。这些实施方式能够以其他的各种方式进行实施,在不脱离发明的要旨的范围内,能够进行各种的省略、置换、变更。这些实施方式或其变形与包含于发明的范围或要旨中一样,包含于权利要求书记载的发明及其均等的范围中。
Claims (10)
1.一种偶发同时计数推定方法,在正电子放射断层摄影列表模式数据中推定偶发同时计数,其中,包含:
取得包含飞行时间即TOF信息的TOF列表模式计数数据;
不使用上述TOF信息地将上述取得的TOF列表模式计数数据转换成包含偶发同时计数值的四维原始正弦图计数数据;
为了生成四维插补正弦图计数数据,对上述四维原始正弦图计数数据进行插补;
为了除去噪音,对上述四维插补正弦图计数数据使用低通滤波器;
将使用上述低通滤波器得到的四维插补正弦图计数数据转换成使用滤波器得到的四维原始正弦图计数数据;以及
通过处理装置,对使用上述滤波器得到的四维原始正弦图计数数据有效地使用用于对湮没点进行滤波的TOF蒙片滤波器,从而根据使用上述滤波器得到的四维正弦图计数数据来生成五维TOF原始正弦图计数数据。
2.根据权利要求1所述的偶发同时计数推定方法,其中,
上述TOF蒙片滤波器沿着同时计数线即LOR,在重建视野外具有规定范围的滤波器区域。
3.根据权利要求1或2所述的偶发同时计数推定方法,其中,
当生成上述五维TOF原始正弦图计数数据时,根据以下的公式在每个TOF组距t生成五维TOF原始正弦图计数数据
其中,上述tcoin是同时计数窗口的mm显示的尺寸,上述tbinsz是mm显示的TOF组距的尺寸,上述r(rad,phi,slice)是使用滤波器得到的四维原始正弦图随机计数数据,上述(rad,phi,slice)是原始正弦图空间中的LOR的坐标,上述θ是LOR的倾斜角。
4.根据权利要求1或2所述的偶发同时计数推定方法,其中,
在上述TOF列表模式计数数据的取得中,对每个同时测量事件取得具有数据排列{xa,za,xb,zb,ea,eb,tof}的TOF列表模式数据,xa、xb是事件a和事件b各自的横截方向的入射闪烁体的编号,za、zb是事件a和事件b各自的轴方向的入射闪烁体的编号,ea、eb是事件a和事件b各自的能级,tof是事件a和事件b的到达时间的差。
5.根据权利要求1或2所述的偶发同时计数推定方法,其中,
使用延迟同时计数窗口法在上述四维原始正弦图计数数据内推定偶发同时计数值。
6.一种偶发同时计数推定装置,在正电子放射断层摄影列表模式数据中推定偶发同时计数,其中,具备存储器和处理装置,
上述存储器存储包含飞行时间即TOF信息的TOF列表模式计数数据;
上述处理装置被构成为:
不使用TOF信息地将上述取得的TOF列表模式计数数据转换成包含偶发同时计数值的四维原始正弦图计数数据;
为了生成四维插补正弦图计数数据,对上述四维原始正弦图计数数据进行插补;
为了除去噪音,对上述四维插补正弦图计数数据使用低通滤波器,
将使用上述低通滤波器得到的四维插补正弦图计数数据转换成使用滤波器得到的四维原始正弦图计数数据;以及
对使用上述滤波器得到的四维原始正弦图计数数据,有效地使用用于对湮没点进行滤波的TOF蒙片滤波器,从而根据使用上述滤波器得到的四维正弦图计数数据来生成五维TOF原始正弦图计数数据。
7.根据权利要求6所述的偶发同时计数推定装置,其中,
上述TOF蒙片滤波器沿着同时计数线即LOR,在重建视野外具有规定范围的滤波器区域。
9.根据权利要求6或7所述的偶发同时计数推定装置,其中,
上述存储器对每个同时测量事件存储具有数据排列{xa,za,xb,zb,ea,eb,tof}的TOF列表模式数据,xa、xb是事件a和事件b各自的横截方向的入射闪烁体的编号,za、zb是事件a和事件b各自的轴方向的入射闪烁体的编号,ea、eb是事件a和事件b各自的能级,tof是事件a和事件b的到达时间的差。
10.根据权利要求6或7所述的偶发同时计数推定装置,其中,
上述处理装置还构成为使用延迟同时计数窗口法在上述四维原始正弦图计数数据内推定偶发同时计数值。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/467,090 | 2012-05-09 | ||
US13/467,090 US9241678B2 (en) | 2012-05-09 | 2012-05-09 | Random estimation in positron emission tomography with tangential time-of-flight mask |
JP2013-093004 | 2013-04-25 | ||
JP2013093004A JP6125309B2 (ja) | 2012-05-09 | 2013-04-25 | 偶発同時計数推定方法及び偶発同時計数推定装置 |
PCT/JP2013/063081 WO2013168778A1 (ja) | 2012-05-09 | 2013-05-09 | 偶発同時計数推定方法及び偶発同時計数推定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103890609A true CN103890609A (zh) | 2014-06-25 |
CN103890609B CN103890609B (zh) | 2016-03-23 |
Family
ID=49549306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380003466.0A Active CN103890609B (zh) | 2012-05-09 | 2013-05-09 | 偶发同时计数推定方法以及偶发同时计数推定装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9241678B2 (zh) |
JP (1) | JP6125309B2 (zh) |
CN (1) | CN103890609B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2011253779A1 (en) * | 2011-12-01 | 2013-06-20 | Canon Kabushiki Kaisha | Estimation of shift and small image distortion |
US11513243B2 (en) * | 2017-03-09 | 2022-11-29 | Shimadzu Corporation | Scatter estimation method, scatter estimation program, and positron CT device having same installed thereon |
JP6974159B2 (ja) * | 2017-12-26 | 2021-12-01 | 浜松ホトニクス株式会社 | 画像処理装置および画像処理方法 |
CN111684493B (zh) * | 2018-01-30 | 2024-08-20 | 皇家飞利浦有限公司 | 用于在列表模式重建中量化准确度提高的校正方法 |
CN109188489B (zh) * | 2018-06-20 | 2020-04-28 | 江苏赛诺格兰医疗科技有限公司 | 随机事例的估计方法、装置及计算机可读存储介质 |
CN110215227B (zh) * | 2019-06-05 | 2022-10-14 | 上海联影医疗科技股份有限公司 | 时间窗设置方法、装置、计算机设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07113873A (ja) * | 1993-10-14 | 1995-05-02 | Univ Tohoku | 陽電子断層撮影装置におけるγ線吸収体による散乱同時計数測定法及び陽電子断層撮影装置 |
CN101234024A (zh) * | 2007-01-30 | 2008-08-06 | 上海西门子医疗器械有限公司 | 一种计算机断层成像的扫描方法 |
WO2011117990A1 (ja) * | 2010-03-25 | 2011-09-29 | 独立行政法人放射線医学総合研究所 | Pet装置における同時計数判定方法及び装置 |
JP4933767B2 (ja) * | 2005-08-11 | 2012-05-16 | 株式会社島津製作所 | 放射線同時計数処理方法、放射線同時計数処理プログラムおよび放射線同時計数処理記憶媒体、並びに放射線同時計数装置およびそれを用いた核医学診断装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3610655B2 (ja) | 1995-12-31 | 2005-01-19 | 株式会社島津製作所 | ポジトロンect装置 |
US7649176B2 (en) * | 2003-11-14 | 2010-01-19 | Siemens Medical Solutions Usa, Inc. | Method for improving clinical data quality in positron emission tomography |
US20080317194A1 (en) | 2005-11-10 | 2008-12-25 | Koninklijke Philips Electronics N. V. | Pet Imaging Using Anatomic List Mode Mask |
US7414246B2 (en) * | 2006-01-03 | 2008-08-19 | Koninklijke Philips Electronics N.V. | Achieving accurate time-of-flight calibrations with a stationary coincidence point source |
US8110805B2 (en) | 2008-05-09 | 2012-02-07 | Siemens Medical Solutions Usa, Inc. | Iterative algorithms for crystal efficiencies estimations from TOF compressed normalization data |
-
2012
- 2012-05-09 US US13/467,090 patent/US9241678B2/en active Active
-
2013
- 2013-04-25 JP JP2013093004A patent/JP6125309B2/ja active Active
- 2013-05-09 CN CN201380003466.0A patent/CN103890609B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07113873A (ja) * | 1993-10-14 | 1995-05-02 | Univ Tohoku | 陽電子断層撮影装置におけるγ線吸収体による散乱同時計数測定法及び陽電子断層撮影装置 |
JP4933767B2 (ja) * | 2005-08-11 | 2012-05-16 | 株式会社島津製作所 | 放射線同時計数処理方法、放射線同時計数処理プログラムおよび放射線同時計数処理記憶媒体、並びに放射線同時計数装置およびそれを用いた核医学診断装置 |
CN101234024A (zh) * | 2007-01-30 | 2008-08-06 | 上海西门子医疗器械有限公司 | 一种计算机断层成像的扫描方法 |
WO2011117990A1 (ja) * | 2010-03-25 | 2011-09-29 | 独立行政法人放射線医学総合研究所 | Pet装置における同時計数判定方法及び装置 |
Non-Patent Citations (1)
Title |
---|
WENLI WANG ET AL.: "4-D Geometric Sensitivity for Multi-Headed Planar Detector PET Systems", 《2001 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD》, vol. 2, 10 November 2001 (2001-11-10), pages 1104 - 1108 * |
Also Published As
Publication number | Publication date |
---|---|
US9241678B2 (en) | 2016-01-26 |
JP2013234995A (ja) | 2013-11-21 |
CN103890609B (zh) | 2016-03-23 |
US20130304386A1 (en) | 2013-11-14 |
JP6125309B2 (ja) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103890609B (zh) | 偶发同时计数推定方法以及偶发同时计数推定装置 | |
EP3033637B1 (en) | A hybrid method based on simulation and experimental data to normalize pet data | |
US10126445B2 (en) | Radiation detecting element sensitivity correction method and radiation tomography device | |
RU2620863C2 (ru) | Способ энергетической калибровки цифровой позитронно-эмиссионной томографии (dpet) | |
JP2018506706A (ja) | Petシステム、デッドタイム補正因子の算出方法 | |
CN103547942B (zh) | 随机事件减少方法以及随机事件减少装置 | |
US10379228B2 (en) | Pet detector scintillator arrangement with light sharing and depth of interaction estimation | |
US20180114346A1 (en) | Reconstructing pet image | |
US10215864B2 (en) | System and method to improve image quality of emission tomography when using advanced radionuclides | |
US20170053423A1 (en) | Systems and methods for emission tomography quantitation | |
US20170192107A1 (en) | Apparatus, method and system for sparse detector | |
US10482634B2 (en) | Systems and methods for imaging with anisotropic voxels | |
CN105125231A (zh) | 一种pet图像环状伪影的去除方法和装置 | |
US11513243B2 (en) | Scatter estimation method, scatter estimation program, and positron CT device having same installed thereon | |
CN109887048B (zh) | Pet散射校正方法、图像重建方法、装置及电子设备 | |
CN108474862A (zh) | 具有lu谱减除的能量校准 | |
US7129497B2 (en) | Method and system for normalization of a positron emission tomography system | |
EP2910975B1 (en) | Arrangement determining method of detector module, gamma ray detector, and pet apparatus | |
CN104024885B (zh) | 核医学诊断装置及医用数据处理装置 | |
Gaens et al. | GPU-accelerated Monte Carlo based scatter correction in brain PET/MR | |
CN104146726B (zh) | Pet系统符合探测响应的生成方法 | |
US10359519B2 (en) | Nuclear medicine diagnostic apparatus and calibration method | |
KR20200086814A (ko) | 양전자 방출 단층 촬영기의 비행시간 정보를 이용하여 관측된 데이터에 포함된 반응깊이 정보를 보정하여 영상을 재구성하는 영상재구성방법 | |
CN112998737B (zh) | 一种扫描装置的时间偏移校正系统及其时间偏移校正方法 | |
US20200085397A1 (en) | Attenuation coefficient image estimation method, attenuation coefficient image estimation program, and positron ct apparatus equipped with the same |
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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160713 Address after: Japan Tochigi Tian Yuan City Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |