CN103126701A - 正电子发射计算机断层摄影装置和图像处理装置 - Google Patents
正电子发射计算机断层摄影装置和图像处理装置 Download PDFInfo
- Publication number
- CN103126701A CN103126701A CN2012104967058A CN201210496705A CN103126701A CN 103126701 A CN103126701 A CN 103126701A CN 2012104967058 A CN2012104967058 A CN 2012104967058A CN 201210496705 A CN201210496705 A CN 201210496705A CN 103126701 A CN103126701 A CN 103126701A
- Authority
- CN
- China
- Prior art keywords
- count information
- subject
- moving
- positron emission
- over time
- 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
- 238000002591 computed tomography Methods 0.000 title abstract description 5
- 230000005855 radiation Effects 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 54
- 238000002600 positron emission tomography Methods 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 7
- 230000002123 temporal effect Effects 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 24
- 230000029058 respiratory gaseous exchange Effects 0.000 description 15
- 230000005251 gamma ray Effects 0.000 description 12
- FVAUCKIRQBBSSJ-UHFFFAOYSA-M sodium iodide Chemical compound [Na+].[I-] FVAUCKIRQBBSSJ-UHFFFAOYSA-M 0.000 description 9
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000000241 respiratory effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 235000009518 sodium iodide Nutrition 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 208000019901 Anxiety disease Diseases 0.000 description 1
- 229910052765 Lutetium Inorganic materials 0.000 description 1
- OSGMVZPLTVJAFX-UHFFFAOYSA-N [Gd].[Lu] Chemical compound [Gd].[Lu] OSGMVZPLTVJAFX-UHFFFAOYSA-N 0.000 description 1
- ANDNPYOOQLLLIU-UHFFFAOYSA-N [Y].[Lu] Chemical compound [Y].[Lu] ANDNPYOOQLLLIU-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000036506 anxiety Effects 0.000 description 1
- 229910052797 bismuth Inorganic materials 0.000 description 1
- JCXGWMGPZLAOME-UHFFFAOYSA-N bismuth atom Chemical compound [Bi] JCXGWMGPZLAOME-UHFFFAOYSA-N 0.000 description 1
- ORCSMBGZHYTXOV-UHFFFAOYSA-N bismuth;germanium;dodecahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.[Ge].[Ge].[Ge].[Bi].[Bi].[Bi].[Bi] ORCSMBGZHYTXOV-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- OHSVLFRHMCKCQY-UHFFFAOYSA-N lutetium atom Chemical compound [Lu] OHSVLFRHMCKCQY-UHFFFAOYSA-N 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001360 synchronised effect Effects 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)
-
- 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
- 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/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1647—Processing of scintigraphic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
实施方式涉及正电子发射计算机断层摄影装置和图像处理装置。提供一种能够简易地取得被检体的体动信息的正电子发射计算机断层摄影装置和图像处理装置。实施方式涉及的正电子发射计算机断层摄影装置具备检测器、同时计数信息生成部和体动检测部。检测器检测从被检体发射出的湮没放射线。同时计数信息生成部从根据检测器的输出信号而生成的计数信息的列表中，检索大致同时计数了一对湮没放射线而得到的计数信息的组，通过按照检索到的每个计数信息的组生成同时计数信息，来生成同时计数信息的时间序列列表。体动检测部基于同时计数信息的时间序列列表，检测被检体的体动随时间的变化。
Description
相关申请的交叉引用
本申请基于并要求2011年11月30日提出的日本专利申请号2011-262756的优先权，本申请中引用该日本专利申请的全部内容。
技术领域
本发明的实施方式涉及一种正电子发射计算机断层摄影装置和图像处理装置。
背景技术
以前，作为核医学成像装置，已知有正电子发射计算机断层摄影装置(以下适当称为“PET(Positron Emission computed Tomography)装置”)。在PET装置的摄影中，向被检体供给用正电子发射核素标识的放射性医药品。于是，被选择性地取入到被检体内生物组织中的正电子发射核素就发射出正电子，发射出的正电子与电子相结合而湮没。此时，正电子向大致相反方向发射一对湮没放射线(以下适当称为“湮没伽马(gamma)射线”)。另一方面，PET装置使用环形地配置在被检体周围的检测器来检测该湮没伽马射线，根据检测结果，生成同时计数信息的时间序列列表(也称作“符合(Coincidence)列表”。然后，PET装置使用同时计数信息的时间序列列表进行利用逆投影处理的重构，生成PET图像。
且说，在PET装置的摄影中存在与被检体的呼吸同步地进行摄影的情况。在该呼吸同步摄影中，除了PET装置主体以外，还使用用于取得被检体体动信息的外部装置。所述外部装置例如是呼吸位移监视器。呼吸位移监视器例如在摄影中向放置在被检体胸部上的红外线反射标记(marker)照射红外线，并利用反射照相机对其进行摄影，通过跟踪标记图像来取得呼吸位移的信息。可是，在使用这样的外部装置的情况下，必须个别地进行外部装置的设置或调整，此外，由于红外线反射标记在摄影中会产生偏移或脱落这样的故障，因此很麻烦。
发明内容
本发明所要解决的问题在于，提供一种能够简易地取得被检体的体动信息的正电子发射计算机断层摄影装置和图像处理装置。
实施方式涉及的正电子发射计算机断层摄影装置具备检测器、同时计数信息生成部和体动检测部。检测器检测从被检体发射出的湮没放射线。同时计数信息生成部从根据检测器的输出信号而生成的计数信息的列表中，检索大致同时计数了一对湮没放射线而得到的计数信息的组，通过按照检索到的每个计数信息的组生成同时计数信息，来生成同时计数信息的时间序列列表。体动检测部基于同时计数信息的时间序列列表，检测被检体的体动随时间的变化。
发明效果
根据实施方式的磁共振成像装置，能够简易地取得被检体的体动信息。
附图说明
图1是示出第一实施方式涉及的PET装置的结构的框图。
图2是用于说明第一实施方式涉及的检测器模块的图。
图3是用于说明第一实施方式涉及的数据存储部的图。
图4是用于说明第一实施方式中的计数信息的列表的图。
图5是用于说明第一实施方式中的同时计数信息的时间序列列表的图。
图6是用于说明第一实施方式涉及的体动检测部的图。
图7是用于说明第一实施方式涉及的体动检测部的图。
图8是用于说明第一实施方式的变形例涉及的体动检测部的图。
图9是用于说明第二实施方式涉及的体动检测部的图。
具体实施方式
以下说明实施方式涉及的正电子发射计算机断层摄影装置和图像处理装置。
(第一实施方式)
第一实施方式涉及的PET装置100不使用呼吸位移监视器等外部装置，而是从通过摄影收集到的扫描数据自身中取得被检体的体动信息。具体而言，第一实施方式涉及的PET装置100根据同时计数信息的时间序列列表中包含的各同时计数信息，并使用TOF(Time Of Flight：传播时间)信息，分别计算偶湮没坐标，基于该多个偶湮没坐标，检测被检体的体动随时间的变化。
所涉及的体动检测由后述的体动检测部26等来实现。以下首先说明第一实施方式涉及的PET装置100的结构，接着详细地说明体动检测部26进行的处理。
图1是示出第一实施方式涉及的PET装置100的结构的框图。如图1所示，第一实施方式涉及的PET装置100具备架台装置10和控制台装置20。
架台装置10利用环形地环绕被检体P的周围而配置的检测器，检测从被检体P内的正电子发射的一对湮没伽马射线，并根据检测器的输出信号，生成计数信息并对其进行收集。如图1所示，架台装置10具备顶板11、诊视床12、诊视床驱动部13、多个检测器模块14和计数信息收集部15。再有，架台装置10如图1所示地具有成为摄影口的空洞。
顶板11是载置被检体P的床，配置在诊视床12之上。诊视床驱动部13在后述的诊视床控制部23的控制下，使顶板11移动。例如，诊视床驱动部13通过使顶板11移动，来使被检体P移动到架台装置10的摄影口内。
检测器模块14检测从被检体P内的正电子发射出的湮没伽马射线。如图1所示，将多个检测器模块14配置成环形地环绕被检体P的周围。
图2是用于说明第一实施方式涉及的检测器模块14的图。如图2所示，检测器模块14是光子计数方式的安杰(Anger)型检测器，具有闪烁体141、光电倍增管(PMT(Photomultiplier Tube))142、和光波导143。
闪烁体141将从被检体P内的正电子发射后射入的湮没伽马射线变换成可见光，并输出变换后的可见光(以下适当称为“闪烁光”)。闪烁体141例如由NaI(Sodium Iodide：碘化钠)、BGO(Bismuth Germanate：锗酸铋)、LYSO(Lutetium Yttrium Oxyorthosilicate)、LSO(Lutetium Oxyorthosilicate)、LGSO(Lutetium Gadolinium Oxyorthosilicate)等闪烁晶体形成，如图2所示地二维排列。光电倍增管142将从闪烁体141输出的闪烁光倍增后变换成电信号。如图2所示，配置有多个光电倍增管142。光波导143将从闪烁体141输出的闪烁光传导到光电倍增管142。例如，由光透射性良好的塑料原料等形成光波导143。
再有，光电倍增管142具有接收闪烁光后使其产生光电子的光电阴极、提供用于加速所产生的光电子的电场的多段倍增极、以及作为电子的流出口的阳极。因光电效应而从光电阴极发射出的电子，朝向倍增极加速并与倍增极的表面碰撞，击出多个电子。通过遍及多段倍增极重复该现象，电子数雪崩式倍增，阳极上的电子数达到大约100万。在有关的例子中，光电倍增管142的增益率为100万倍。此外，为了利用雪崩现象进行放大，在倍增极与阳极之间通常施加1000伏以上的电压。
这样地，检测器模块14利用闪烁体141，将从被检体P内的正电子发射的湮没伽马射线变换成闪烁光，并且利用光电倍增管142，将变换后的闪烁光变换成电信号，由此检测从被检体P发射的湮没伽马射线。
回到图1，计数信息收集部15根据检测器模块14的输出信号，生成计数信息，并将生成的计数信息保存在后述的数据存储部24中。再有，图1中省略了图示，多个检测器模块14被区分成多个块，在每块中具备计数信息收集部15。例如，在第一实施方式中，1个检测器模块14是1个块，因此，在每个检测器模块14中都具备计数信息收集部15。
计数信息收集部15将检测器模块14的输出信号变换成数字数据，生成计数信息。该计数信息中包含湮没伽马射线的检测位置、能量值以及检测时间。例如，计数信息收集部15确定多个将闪烁光按相同定时变换成电信号的光电倍增管142。然后，计数信息收集部15使用所确定的各光电倍增管142的位置和电信号的强度，计算重心的位置，并确定示出湮没伽马射线射入的闪烁体141的位置的闪烁体号码(P)。再有，在光电倍增管142是位置检测型光电倍增管的情况下，光电倍增管142也可以确定位置。
此外，计数信息收集部15还通过对从各光电倍增管142输出的电信号的强度进行积分计算，来确定射入到检测器模块14中的湮没伽马射线的能量值(E)。此外，计数信息收集部15还确定由检测器模块14检测出湮没伽马射线的检测时间(T)。例如，计数信息收集部15以10-12秒(微微(pico)秒)为单位的精度下确定检测时间(T)。再有，检测时间(T)可以是绝对时刻，也可以是从摄影开始时刻起的经过时间。这样地，计数信息收集部15就生成包含闪烁体号码(P)、能量值(E)和检测时间(T)的计数信息。
控制台装置20受理操作者对PET装置100的操作来控制PET图像的摄影，并且，使用由架台装置10收集到的计数信息重构PET图像。如图1所示，控制台装置20具备输入部21、显示部22、诊视床控制部23、数据存储部24、同时计数信息生成部25、体动检测部26、图像重构部27和系统控制部28。再有，控制台装置20所具备的各部分经由总线(bus)相连接。
输入部21是由PET装置100的操作者在各种指示或各种设定的输入中使用的鼠标或键盘等，将输入的各种指示或各种设定传送到系统控制部28。例如，将输入部21用于摄影开始指示的输入。显示部22是由操作者参照的监视器等，在系统控制部28的控制下，显示被检体的呼吸波形或PET图像，或者显示用于从操作者受理各种指示或各种设定的GUI(Graphical User Interface：图形用户界面)。诊视床控制部23控制诊视床驱动部13。
数据存储部24对在PET装置100中使用的各种数据进行存储。图3是用于说明第一实施方式涉及的数据存储部24的图。如图3所示，数据存储部24具备计数信息存储部24a、同时计数信息存储部24b和PET图像存储部24c。再有，数据存储部24例如由RAM(Random Access Memory:随机存储器)、快闪存储器等半导体存储器元件或硬盘、光盘等来实现。
计数信息存储部24a存储由各计数信息收集部15收集到的计数信息的列表。此外，计数信息存储部24a所存储的计数信息的列表还用于同时计数信息生成部25的处理中。再有，计数信息存储部24a所存储的计数信息的列表，可以在用于同时计数信息生成部25的处理之后被删除，也可以存储规定期间。
图4是用于说明第一实施方式中的计数信息的列表的图。如图4所示，计数信息存储部24a与识别检测器模块14的模块ID相对应地存储包含闪烁体号码(P)、能量值(E)和检测时间(T)在内的计数信息。
同时计数信息存储部24b存储由同时计数信息生成部25生成的同时计数信息的时间序列列表。此外，同时计数信息存储部24b所存储的同时计数信息的时间序列列表用于体动检测部26或图像重构部27的处理中。再有，同时计数信息存储部24b所存储的同时计数信息的时间序列列表也可以在用于图像重构部27的处理之后删除，也可以存储规定期间。
图5是用于说明第一实施方式中的同时计数信息的时间序列列表的图。如图5所示，同时计数信息存储部24b与作为同时计数信息的编号的符合No.相对应地存储计数信息的组。再有，在第一实施方式中，同时计数信息的时间序列列表根据计数信息的检测时间(T)，大致按时间序列的顺序排列。
PET图像存储部24c存储由图像重构部27重构的PET图像。此外，PET图像存储部24c所存储的PET图像由系统控制部28显示在显示部22中。
回到图1，同时计数信息生成部25使用由计数信息收集部15收集到的计数信息的列表，生成同时计数信息的时间序列列表。具体而言，同时计数信息生成部25从计数信息存储部24a中存储的计数信息的列表中，基于计数信息的检测时间(T)，检索大致同时计数了一对湮没伽马射线而得到的计数信息的组。此外，同时计数信息生成部25按照检索到的每个计数信息的组生成同时计数信息，并将生成的同时计数信息按大致时间序列的顺序排列，同时保存在同时计数信息存储部24b中。
例如，同时计数信息生成部25基于操作者输入的同时计数信息生成条件，生成同时计数信息。在同时计数信息生成条件中指定了时间窗宽度。例如，同时计数信息生成部25基于时间窗宽度“600微微秒”，生成同时计数信息。
例如，同时计数信息生成部25参照计数信息存储部24a，在检测器模块14间检索检测时间(T)的时间差处于时间窗宽度“600微微秒”以内的计数信息的组。例如，同时计数信息生成部25在检索“P11、E11、T11”和“P22、E22、T22”的组作为满足同时计数信息生成条件的组时，生成该组作为同时计数信息，并存储在同时计数信息存储部24b中。再有，同时计数信息生成部25也可以与时间窗宽度一起使用能量窗宽度来生成同时计数信息。
体动检测部26基于由同时计数信息生成部25生成的同时计数信息的时间序列列表，检测被检体P的体动随时间的变化。再有，对于体动检测部26的处理，以后详细叙述。
图像重构部27重构PET图像。具体而言，图像重构部27读取同时计数信息存储部24b中存储的同时计数信息的时间序列列表，作为投影数据(正弦图(sinogram)数据)，通过将读取的投影数据进行逆投影处理来重构PET图像。此外，图像重构部27将重构的PET图像存储在PET图像存储部24c中。
在此，图像重构部27也可以使用由体动检测部26检测到的体动随时间的变化来重构PET图像。例如，图像重构部27根据由体动检测部26检测到的体动随时间的变化，分析被检体P的呼吸周期，按照该呼吸周期，将同时计数信息的时间序列列表分割成多个阶段(phase)(例如，将1个呼吸周期分割成6～15个阶段)。然后，图像重构部27使用各阶段中包含的同时计数信息，按照每个阶段重构PET图像。由于考虑到各阶段内的被检体的体动很小，因此，成为仅根据不怎么受体动影响的同时计数信息来重构1个PET图像，能够有助于PET图像的画质提高。
系统控制部28通过控制架台装置10和控制台装置20来进行PET装置100的整体控制。例如，系统控制部28控制PET装置100中的摄影。再有，上述的同时计数信息生成部25、体动检测部26、图像重构部27和系统控制部28等各部分，由ASIC(Application Specific Integrated Circuit：专用集成电路)或FPGA(Field Programmable Gate Array：现场可编程门阵列)等集成电路、CPU(Central Processing Unit：中央处理单元)或MPU(Micro Processing Unit：微处理单元)等电子电路来实现。
(第一实施方式中的体动检测)
接着，说明第一实施方式中的体动检测。如上所述，第一实施方式涉及的体动检测部26，基于同时计数信息的时间序列列表，检测被检体P的体动随时间的变化。具体而言，体动检测部26根据同时计数信息的时间序列列表中包含的各同时计数信息，分别计算LOR(Line Of Response：线响应)上的偶湮没坐标，并基于计算出的多个偶湮没坐标，检测被检体P的体动随时间的变化。
图6和图7是用于说明第一实施方式涉及的体动检测部26的图。在图6中概念性地示出大致按时间序列的顺序排列的同时计数信息。例如，体动检测部26如图6所示地按照每个规定时间区间(以下适当称为“帧(frame)”)分割同时计数信息的时间序列列表。在求出呼吸相位来作为体动随时间的变化的情况下，期望该帧比被检体P的呼吸周期短、且比心周期长，例如是1～2秒。
接着，体动检测部26对每个规定的帧，根据该帧中包含的全部同时计数信息的各个同时计数信息，分别计算LOR上的偶湮没坐标，并根据计算出的多个偶湮没坐标，计算偶湮没坐标的平均(偶湮没坐标分布的重心坐标)。再有，体动检测部26使用TOF信息计算偶湮没坐标。即，体动检测部26根据同时计数信息中包含的一组计数信息，计算检测时间差，并使用该检测时间差，计算LOR上的正电子的空间位置、即偶湮没坐标。LOR是连结与一对湮没伽马射线相对应的一组检测位置的线。
然后，体动检测部26针对各帧计算偶湮没坐标的重心坐标，检测重心坐标的随时间的变化。例如，体动检测部26选择偶湮没坐标的重心坐标(x坐标、y坐标、z坐标)中的随时间的变化最大的坐标轴，通过标绘所选择的坐标轴的重心坐标来求出呼吸相位。例如，体动检测部26选择z轴的重心坐标，如图7所示地在以横轴为时间、以纵轴为z轴的重心坐标的图表上，沿着时间轴标绘针对各帧计算出的z轴的重心坐标。于是，就如图7所示地在图表上描画出被检体P的呼吸波形。再有，在图7中，图表上部记载的模式图是被检体P的胸部剖面的模式图。记号×表示最大呼气时的z轴的重心坐标和最大吸气时的z轴的重心坐标。此外，例如也可以将图7所示的图表显示在显示部22中，或者也可以将其作为内部计算处理加以执行。
再有，在上述中，体动检测部26选择了偶湮没坐标的重心坐标(x坐标、y坐标、z坐标)中的随时间的变化最大的坐标轴，但实施方式不限于此。例如，体动检测部26也可以确定随时间的变化最大的方向，设定向所确定的方向旋转后的新的坐标轴，通过标绘该新坐标轴的重心坐标来求出呼吸相位。
如上所述，根据第一实施方式，由于从通过摄影收集到的扫描数据自身中取得被检体的体动信息，因此，能够简易地取得被检体的体动信息。此外，由于体动信息包含在扫描数据自身中，因此能够根据需要反复使用。例如能够在多次的分析中使用。
(第一实施方式的变形例)
在第一实施方式中说明了体动检测部26根据多个偶湮没坐标计算偶湮没坐标的平均(重心坐标)的例子，但实施方式不限于此。例如，体动检测部26也可以计算偶湮没坐标的标准偏差作为统计量。
图8是用于说明第一实施方式的变形例涉及的体动检测部的图。体动检测部26对每个规定的帧，根据该帧中包含的全部同时计数信息的各个同时计数信息，分别计算LOR上的偶湮没坐标，根据计算出的偶湮没坐标，计算偶湮没坐标的分布的标准偏差。标准偏差是指某个帧中包含的多个偶湮没坐标的扩大(图8中参照箭头)。并且，体动检测部26针对全部帧计算偶湮没坐标的标准偏差，通过沿着时间轴标绘计算出的标准偏差，就可以求出呼吸相位。例如，体动检测部26如图8所示地在以横轴为时间、以纵轴为标准偏差的图表上标绘针对各帧计算出的z轴的标准偏差。于是，就如图8所示地在图表上描绘出被检体P的呼吸波形。
(第二实施方式)
下面说明第二实施方式。在第一实施方式中说明了使用同时计数信息的时间序列列表中包含的全部同时计数信息来检测体动随时间的变化的方式，但实施方式不限于此。在第二实施方式中，对选择一部分同时计数信息来检测体动随时间的变化的方式进行说明。
图9是用于说明第二实施方式涉及的体动检测部26的图。图9是被检测器环绕的被检体P的胸部剖面的模式图，黑圆a表示被检体P的心脏。此外，在第二实施方式中，求出心相位来作为体动随时间的变化。
例如，体动检测部26从同时计数信息的时间序列列表中，仅选择LOR经过被检体P内心脏的多个同时计数信息。在此，可以根据顶板11上载置的被检体P的姿势或被检体P的体型等，对经过心脏的LOR(即、检测位置的组)大致预先确定。因此，体动检测部26从同时计数信息的时间序列列表中，仅选择与预先确定的检测位置的组相对应的同时计数信息。例如，图9中，a1～a4的LOR是经过心脏的LOR。
然后，体动检测部26根据所选择的各同时计数信息，分别计算LOR上的偶湮没坐标，与第一实施方式同样，基于计算出的多个偶湮没坐标，检测被检体P的体动随时间的变化。
(第二实施方式的变形例)
此外，作为第二实施方式的变形例，例如，说明求出呼吸相位来作为体动随时间的变化但去掉心相位的影响的方式。在第一实施方式中说明了能够求出呼吸相位来作为体动随时间的变化的这点，但在放射性医药品集聚到心脏的情况下等，具有会在体动随时间的变化中反映心相位的影响之虞。该情况下，就会在示出呼吸相位的呼吸波形上，作为噪声而附有示出比呼吸周期短的周期的心相位的波形。
因此，例如，体动检测部26从同时计数信息的时间序列列表中，仅选择LOR不经过被检体P内心脏的多个同时计数信息。即，体动检测部26从同时计数信息的时间序列列表中，除去与作为经过心脏的LOR而预先确定的检测位置的组相对应的同时计数信息，而选择剩余的同时计数信息。例如，在图9中，除a1～a4以外的LOR是不经过心脏的LOR。
然后，体动检测部26根据所选择的各同时计数信息，分别计算LOR上的偶湮没坐标，与第一实施方式同样，基于计算出的多个偶湮没坐标，检测被检体P的体动随时间的变化。
如上所述，根据第二实施方式，由于选择一部分同时计数信息来检测体动随时间的变化，因此，能够着眼于特定部位来检测体动随时间的变化，或者去掉确定部位所产生的影响来检测体动随时间的变化。
(其他实施方式)
再有，在上述实施方式中，作为PET装置100的结构，例示了图1，但实施方式不限于此。例如，可以将计数信息收集部15设置在控制台装置20一侧，反之，也可以将同时计数信息生成部25设置在架台装置10一侧。此外，也可以将数据存储部24中存储的各种数据设置在架台装置10一侧，或者设置在控制台装置20一侧。各个数据保持在PET装置100中的期间也任意。
此外，例如也可以将体动检测部26设置在与PET装置100主体不同的图像处理装置中。例如，该情况下，图像处理装置具备存储同时计数信息的时间序列列表的同时计数信息存储部、基于同时计数信息的时间序列列表，检测被检体的体动随时间的变化的体动检测部、以及使用被检体的体动随时间的变化，从同时计数信息的时间序列列表重构被检体的图像的图像重构部。再有，也可以将被检体的体动随时间的变化用于除图像的重构以外的其他目的。该情况下，也不需要图像重构部。
根据以上所述的至少一个实施方式的正电子发射计算机断层摄影装置和图像处理装置，能够简易地取得被检体的体动信息。
已经说明了本发明的几个实施方式，但是这些实施方式是作为例子而提出的，并不是想限定发明范围。这些实施方式可以以其他各种各样的方式进行实施，可以在不脱离发明主旨的范围内进行各种各样的省略、置换和变更。这些实施方式或其变形包含在发明范围或主旨内，并且也同样包含在权利要求书记载的发明及其均等的范围内。
Claims (15)
1.一种正电子发射计算机断层摄影装置，具备：
检测器，检测从被检体发射出的湮没放射线；
同时计数信息生成部，从根据所述检测器的输出信号而生成的计数信息的列表中，检索大致同时计数了一对湮没放射线而得到的计数信息的组，通过按照检索到的每个计数信息的组生成同时计数信息，来生成同时计数信息的时间序列列表；和
体动检测部，基于所述同时计数信息的时间序列列表，检测所述被检体的体动随时间的变化。
2.根据权利要求1所述的正电子发射计算机断层摄影装置，所述体动检测部根据所述同时计数信息的时间序列列表中包含的各同时计数信息，分别计算线响应LOR上的偶湮没坐标，基于计算出的多个偶湮没坐标，检测所述被检体的体动随时间的变化。
3.根据权利要求1所述的正电子发射计算机断层摄影装置，所述体动检测部从所述同时计数信息的时间序列列表中，选择线响应LOR经过所述被检体内的规定部位的多个同时计数信息，根据所选择的各同时计数信息，分别计算线响应LOR上的偶湮没坐标，基于计算出的多个偶湮没坐标，检测所述被检体的体动随时间的变化。
4.根据权利要求1所述的正电子发射计算机断层摄影装置，所述体动检测部从所述同时计数信息的时间序列列表中，选择线响应LOR不经过所述被检体内的规定部位的多个同时计数信息，根据所选择的各同时计数信息，分别计算线响应LOR上的偶湮没坐标，基于计算出的多个偶湮没坐标，检测所述被检体的体动随时间的变化。
5.根据权利要求2所述的正电子发射计算机断层摄影装置，所述体动检测部按照每个规定时间区间，计算根据各时间区间中包含的同时计数信息而计算出的偶湮没坐标的统计量，基于按照每个时间区间计算出的统计量随时间的变化，检测所述被检体的体动随时间的变化。
6.根据权利要求3所述的正电子发射计算机断层摄影装置，所述体动检测部按照每个规定时间区间，计算根据各时间区间中包含的同时计数信息而计算出的偶湮没坐标的统计量，基于按照每个时间区间计算出的统计量随时间的变化，检测所述被检体的体动随时间的变化。
7.根据权利要求4所述的正电子发射计算机断层摄影装置，所述体动检测部按照每个规定时间区间，计算根据各时间区间中包含的同时计数信息而计算出的偶湮没坐标的统计量，基于按照每个时间区间计算出的统计量随时间的变化，检测所述被检体的体动随时间的变化。
8.根据权利要求5所述的正电子发射计算机断层摄影装置，所述体动检测部计算每个时间区间的平均或者标准偏差，作为所述统计量。
9.根据权利要求6所述的正电子发射计算机断层摄影装置，所述体动检测部计算每个时间区间的平均或者标准偏差，作为所述统计量。
10.根据权利要求7所述的正电子发射计算机断层摄影装置，所述体动检测部计算每个时间区间的平均或者标准偏差，作为所述统计量。
11.根据权利要求1所述的正电子发射计算机断层摄影装置，还具备图像重构部，所述图像重构部使用所述被检体的体动随时间的变化，从所述同时计数信息的时间序列列表重构所述被检体的图像。
12.根据权利要求2所述的正电子发射计算机断层摄影装置，还具备图像重构部，所述图像重构部使用所述被检体的体动随时间的变化，从所述同时计数信息的时间序列列表重构所述被检体的图像。
13.根据权利要求3所述的正电子发射计算机断层摄影装置，还具备图像重构部，所述图像重构部使用所述被检体的体动随时间的变化，从所述同时计数信息的时间序列列表重构所述被检体的图像。
14.根据权利要求4所述的正电子发射计算机断层摄影装置，还具备图像重构部，所述图像重构部使用所述被检体的体动随时间的变化，从所述同时计数信息的时间序列列表重构所述被检体的图像。
15.一种图像处理装置，具备：
同时计数信息存储部，存储同时计数信息的时间序列列表，所述同时计数信息是大致同时计数了从被检体发射的一对湮没放射线而得到的计数信息的组；
体动检测部，基于所述同时计数信息的时间序列列表，检测所述被检体的体动随时间的变化；和
图像重构部，使用所述被检体的体动随时间的变化，从所述同时计数信息的时间序列列表重构所述被检体的图像。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011-262756 | 2011-11-30 | ||
| JP2011262756 | 2011-11-30 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103126701A true CN103126701A (zh) | 2013-06-05 |
| CN103126701B CN103126701B (zh) | 2016-04-27 |
Family
ID=48465950
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210496705.8A Active CN103126701B (zh) | 2011-11-30 | 2012-11-28 | 正电子发射计算机断层摄影装置和图像处理装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US9239393B2 (zh) |
| JP (1) | JP6253880B2 (zh) |
| CN (1) | CN103126701B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108634974A (zh) * | 2018-04-03 | 2018-10-12 | 沈阳东软医疗系统有限公司 | 一种门控信号的确定方法和装置 |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6545565B2 (ja) * | 2014-08-25 | 2019-07-17 | キヤノンメディカルシステムズ株式会社 | 検出器、核医学イメージング装置、ｐｅｔ−ｃｔ装置及びｐｅｔ−ｍｒｉ装置 |
| JP2023032314A (ja) * | 2021-08-26 | 2023-03-09 | キヤノンメディカルシステムズ株式会社 | 核医学診断装置、データ処理方法及びプログラム |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004111946A1 (en) * | 2003-06-18 | 2004-12-23 | Philips Intellectual Property & Standards Gmbh | Motion compensated reconstruction technique |
| CN101610719A (zh) * | 2007-02-07 | 2009-12-23 | 皇家飞利浦电子股份有限公司 | 治疗计划中的运动估计 |
| US20100166274A1 (en) * | 2006-02-28 | 2010-07-01 | Philips Intellectual Property & Standards Gmbh | Local motion compensation based on list mode data |
| WO2011070465A2 (en) * | 2009-12-10 | 2011-06-16 | Koninklijke Philips Electronics, N.V. | Method and apparatus for using time of flight information to detect and correct for motion in imaging scans |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5251128A (en) * | 1990-11-19 | 1993-10-05 | General Electric Company | Motion artifact reduction in projection imaging |
| US7117026B2 (en) * | 2002-06-12 | 2006-10-03 | Koninklijke Philips Electronics N.V. | Physiological model based non-rigid image registration |
| JP4467987B2 (ja) | 2004-01-05 | 2010-05-26 | 株式会社東芝 | 核医学診断装置 |
| JP2007003264A (ja) | 2005-06-22 | 2007-01-11 | Shimadzu Corp | Ｐｅｔ装置 |
| US20070197927A1 (en) * | 2005-12-22 | 2007-08-23 | Physical Logic Ag | Method for detecting cardiovascular problems using micro or nano vibrations |
| JP2009156856A (ja) | 2007-12-03 | 2009-07-16 | Shimadzu Corp | 核医学診断ｘ線ｃｔ装置 |
-
2012
- 2012-11-28 CN CN201210496705.8A patent/CN103126701B/zh active Active
- 2012-11-29 US US13/688,531 patent/US9239393B2/en active Active
- 2012-11-30 JP JP2012262782A patent/JP6253880B2/ja active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004111946A1 (en) * | 2003-06-18 | 2004-12-23 | Philips Intellectual Property & Standards Gmbh | Motion compensated reconstruction technique |
| US20100166274A1 (en) * | 2006-02-28 | 2010-07-01 | Philips Intellectual Property & Standards Gmbh | Local motion compensation based on list mode data |
| CN101610719A (zh) * | 2007-02-07 | 2009-12-23 | 皇家飞利浦电子股份有限公司 | 治疗计划中的运动估计 |
| WO2011070465A2 (en) * | 2009-12-10 | 2011-06-16 | Koninklijke Philips Electronics, N.V. | Method and apparatus for using time of flight information to detect and correct for motion in imaging scans |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108634974A (zh) * | 2018-04-03 | 2018-10-12 | 沈阳东软医疗系统有限公司 | 一种门控信号的确定方法和装置 |
| CN108634974B (zh) * | 2018-04-03 | 2022-03-04 | 东软医疗系统股份有限公司 | 一种门控信号的确定方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US9239393B2 (en) | 2016-01-19 |
| JP6253880B2 (ja) | 2017-12-27 |
| CN103126701B (zh) | 2016-04-27 |
| US20130134314A1 (en) | 2013-05-30 |
| JP2013137303A (ja) | 2013-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106539591B (zh) | Pet飞行时间状态质量检测方法和pet扫描装置 | |
| EP3508885A1 (en) | Multi-gamma ray photon simultaneous medicine emission time-coincidence nuclear medicine imaging system and method | |
| US9014330B2 (en) | Radiation diagnostic apparatus and control method | |
| CN102439483B (zh) | 核医学成像装置以及核医学成像方法 | |
| US9844351B2 (en) | Positron CT apparatus and a timing correction method | |
| EP2748641A2 (en) | Data-driven optimization of event acceptance/rejection logic | |
| JP2005140783A (ja) | 検出器モジュール | |
| US20160073976A1 (en) | Nuclear medicine diagnostic apparatus and nuclear medicine image generating method | |
| CN102727237B (zh) | 医用图像诊断装置以及控制方法 | |
| JP5595738B2 (ja) | 核医学診断装置 | |
| US10436915B2 (en) | Medical image diagnosis apparatus | |
| JP6545565B2 (ja) | 検出器、核医学イメージング装置、ｐｅｔ−ｃｔ装置及びｐｅｔ−ｍｒｉ装置 | |
| CN103126701A (zh) | 正电子发射计算机断层摄影装置和图像处理装置 | |
| JP6502105B2 (ja) | 核医学画像撮影装置及び制御方法 | |
| JP6878131B2 (ja) | 医用画像診断装置 | |
| US20120228511A1 (en) | Positron emission computed tomography apparatus | |
| JP5900740B2 (ja) | 放射線断層撮像装置用のデータ収集器およびそれを備えた放射線断層撮像装置、放射線断層撮像用データ収集方法 | |
| JP6825706B2 (ja) | 放射線検出装置およびそれを備えた核医学診断装置 | |
| CN102439485A (zh) | 核医学成像装置及核医学成像系统 | |
| JP7221623B2 (ja) | 医用画像診断装置 | |
| Mendes et al. | Evaluation of monolithic detector blocks for high-sensitivity PET imaging of the human brain | |
| Moradi | Front-end data processing of new positron emission tomography demonstrator | |
| JP2010190669A (ja) | 核医学診断装置及び画像処理装置 |
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: 20160705 Address after: Tochigi County, Japan Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |