CN101658426B - 基于三维正电子发射探测器灵敏度特性的呼吸门控技术 - Google Patents
基于三维正电子发射探测器灵敏度特性的呼吸门控技术 Download PDFInfo
- Publication number
- CN101658426B CN101658426B CN2009100949595A CN200910094959A CN101658426B CN 101658426 B CN101658426 B CN 101658426B CN 2009100949595 A CN2009100949595 A CN 2009100949595A CN 200910094959 A CN200910094959 A CN 200910094959A CN 101658426 B CN101658426 B CN 101658426B
- Authority
- CN
- China
- Prior art keywords
- frame
- pet
- sensitivity
- characteristic
- positron emission
- 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
- 230000035945 sensitivity Effects 0.000 title claims abstract description 27
- 230000000241 respiratory effect Effects 0.000 title claims abstract description 26
- 238000005516 engineering process Methods 0.000 title claims abstract description 23
- 230000033001 locomotion Effects 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000003384 imaging method Methods 0.000 claims abstract description 13
- 238000012879 PET imaging Methods 0.000 claims abstract description 12
- 238000012636 positron electron tomography Methods 0.000 claims abstract 9
- 239000000284 extract Substances 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 7
- 230000002195 synergetic effect Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 abstract description 9
- 210000000056 organ Anatomy 0.000 abstract description 6
- 238000009206 nuclear medicine Methods 0.000 abstract description 3
- 230000011218 segmentation Effects 0.000 abstract 1
- 238000002600 positron emission tomography Methods 0.000 description 14
- 210000000481 breast Anatomy 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 5
- 206010028980 Neoplasm Diseases 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000003187 abdominal effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- 230000004323 axial length Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
Images
Landscapes
- Nuclear Medicine (AREA)
Abstract
本发明涉及一种基于三维正电子发射(PET)探测器灵敏度特性的呼吸门控技术，它属于核医学成像领域。本发明提出的门控方法能够有效地补偿呼吸导致的图像伪影从而提高医生的正确诊断率。在三维PET成像中，利用扫描探测器几何灵敏度的变化特性，将获取的动态数据进行帧分割，其帧内光子数可以反映器官运动所在的运动相位位置，提取每一个循环中相同的相位进行叠加和图像重建就可以改善在PET成像中由于呼吸而导致的图像伪影。与现有门控技术相比，本发明方法不涉及硬件设备，不需要患者及临床做额外的扫描前准备，在实际运用中将更简单、更有效和更可靠。
Description
技术领域
本发明涉及一种基于三维正电子发射PET(Positron Emission Tomography)探测器灵敏度特性的呼吸门控技术。属核医学成像技术领域。
背景技术
PET成像主要利用符合线路探测技术和正电子核素及其标记化合物进行成像，其图像可以从分子水平反映人体正常或疾病时的代谢状态。由于变化在疾病早期阶段就可以出现，因此正电子成像能够进行疾病的早期诊断并为治疗决策提供客观依据。这一技术的应用在肺癌患者早期的诊断和治疗计划方面已发挥着重要的作用。
然而，由于PET成像的特殊性，它扫描采集数据的时间相对于CT(ComputerTomography)成像，一般要很长时间。患者不可能像CT扫描成像那样，在短时间内能够屏住呼吸而获得胸腹部清晰的成像。患者的胸腹部PET成像是在扫描时间段内平均呼吸运动叠加情况下的成像，它的成像是模糊的(有伪影)。另外在图像融合时，CT和PET图像上病灶位置的不匹配以及衰减校正后的活性偏差也是产生图像伪影的原因之一。图像伪影给医生对胸腹肿瘤的诊断和治疗带来困难，这就需要有性能稳定可靠，空间分辨率高，简单速度快的方法来补偿由于呼吸运动而造成的PET/CT成像模糊。
目前有许多中技术可用于改进呼吸运动的成像，其中一种是呼吸门控的技术。这种技术是将呼吸运动视为类似的呼吸循环组成，每一个呼吸循环可划分为不同的相位，用一定的测量方法记录在一定相位的呼吸容量和位置，然后提取在相同相位上的所有获取的扫描数据进行叠加成像，就能得到改善的图像质量。
门控技术可以分为外部和内部两类方法。外部门控方法使用一个呼吸运动监测系统来跟踪呼吸运动的参数，然后将其运用到重建算法中进行运动的补偿。这样的呼吸运动跟踪系统已经被应用于PET/CT成像的研究中。例如胸腹传感带运动检测系统，肺呼吸容量检测仪和信号激光位置跟踪系统。但是它具有对硬件设备的额外需求以及测量运动参数的校准、归一化和医疗设备的接口等缺点。另外一种内部门控的方法是利用数据驱动的方法，在已获取的数据中利用一些特定的算法对数据进行分析处理以便估算出呼吸运动的参数，再提取相应的数据进行图像重建。例如中心质点算法就是典型的数据驱动方法。然而，它具有算法负责，计算量大等缺点。
发明内容
本发明的目的在于克服上述现有技术在PET扫描成像过程中，由于患者的呼吸会导致图像的模糊伪影，从而影响医生对病情的诊断和治疗的不足，发明一种在三维PET成像中，利用扫描探测器几何灵敏度的变化特性，将获取的动态数据进行帧分割，其帧内光子数可以反映器官运动所在的运动相位位置，提取每一个循环中相同的相位进行叠加和图像重建，就可改善在PET成像中由于呼吸而导致的图像伪影的基于三维正电子发射(PET)探测器灵敏度特性的呼吸门控技术，它能有效地补偿呼吸导致的图像伪影，提高医生的正确诊断率。
本发明基于三维正电子发射探测器灵敏度特性的呼吸门控技术是这样实现的：在三维PET成像过程中，首先利用动态模式获取成像的数据；然后根据不同的时间分辨率将数据划分为一定的数据帧，每一帧看作一个呼吸循环中的一个运动相位；结合扫描器空间灵敏度的特性就可分析观察出每帧中光子数的周期性变化情况；依据所有帧内光子数的变化，用门控方法进行图像校正，最后获得校正的图像。
动态扫描数据帧内的光子数会反映出呼吸运动沿扫描探测器轴向变化的位置。
对获取的动态模式数据划分为相同时间间隔的数据帧，其每一帧看作一个呼吸循环中的一个运动相位。
对分割的帧光子进行噪声滤波以凸显呼吸循环的特征。
在已经滤波的帧中，寻找那些光子数接近的帧，然后把它们进行叠加。
对叠加的帧进行图像重建。
具体步骤是：
第一步，用PET的动态数据模式获取两个呼吸循环的扫描成像数据；
第二步，将这两个呼吸循环划分成为16个帧，每个循环包含8个帧；
第三步，计算每个帧内的光子数，并画出帧光子曲线图；
第四步，对帧光子曲线图进行噪声滤波以凸显呼吸循环的特征；
第五步，在两个呼吸循环中寻找光子数接近的帧，然后进行提取；
第六步，将提取的帧进行叠加变成一个帧；
第七步，对已叠加的帧进行图像重建，最后获得校正的图像。
本发明基于三维正电子发射(PET)探测器灵敏度特性的呼吸门控技术是在三维PET成像过程中，首先使用动态模式获取成像的数据。然后根据不同的时间分辨率将数据划分为一定的数据帧，每一帧可以看作一个呼吸循环中的一个运动相位。结合扫描器空间灵敏度的特性就可以分析观察出每帧中光子数的周期性变化情况，根据所有帧内光子数的变化就可以使用门控方法进行图像校正，如图3所示。
工作原理：在三维PET成像中，对于获得高质量的图像，扫描探测器几何灵敏度是一个关键因素。它可以被定义为在一定的人体模型中，单位的放射性浓度(每秒光子数/贝克/毫升)可被检测出的光子数。
在三维PET成像中，几何灵敏度的变化主要沿轴线方向，如图3示例所示，当接近轴线中心时，灵敏度值达到最大；而靠近轴线边沿时，灵敏度值接近最小。这就意味着由于灵敏度沿轴向的变化，探测器检测到的那些从人体胸腹部器官发射的光子数目，将会随着灵敏度的变化而变化。当器官由于运动到轴向中心点时，探测到的光子数目最多；而当器官运动到轴向扫描器边沿时，探测到的光子数最少。由此可见，探测到的光子数可以反映器官运动所在的运动相位位置。呼吸是一个反复循环的运动，探测到的光子数也会呈现出周期性的变化，这样就可以提取每一个循环中相似的相位进行门控。
本发明基于三维正电子发射(PET)探测器灵敏度特性的呼吸门控技术与现有门控技术相比，具有不涉及硬件设备，不需要患者及临床做额外的扫描前准备，在实际运用中将更简单、更有效和更可靠，它能有效地补偿呼吸导致的图像伪影，提高医生的正确诊断率。
附图说明
图1为本发明门控方法的具体步骤流程示意图。
图2.为本发明PET扫描探测器灵敏度曲线。
图3为本发明门控方法，在两个呼吸的循环过程中计算每个帧光子的数目并画出其柱形图，它显示帧内光子的变化率是接近肺容量的变化曲线的。
图4为本发明运动振幅为5厘米、半径2厘米，长5厘米的圆柱体模用于仿真肿瘤运动变化的人体模门控校正示意图。
图5为本发明一个小肿瘤位于肺和肝之间的矢状面图，做呼吸运动振幅为2厘米的门控校正示意图。
在图2中扫描器轴向长度＝180毫米，探测器环直径＝800毫米。可见灵敏度在边缘视野的值是零(在-90毫米或90毫米)，而灵敏度的值在中心视野位置达到最大。
在图3中每个柱状代表一个帧(相位)，柱的高度代表帧的光子数，具有相同纹理的柱状即为帧光子数目相同，根据光子数就可以提取相同的帧相位进行门控。图中显示了8个已被门控的帧。
具体实施方式：
结合图1和图3给出一个本发明在核医学三维正电子发射成像中，利用扫描探测器灵敏度特性进行呼吸门控校正，实现方法的具体步骤：
第一步，用PET的动态数据模式获取两个呼吸循环的扫描成像数据；
第二步，将这两个呼吸循环划分成如图2所示的16个帧，每个循环包含8个帧；
第三步，计算每个帧内的光子数，并画出帧光子曲线图；
第四步，对帧光子曲线图进行噪声滤波以凸显呼吸循环的特征；
第五步，在两个呼吸循环中寻找光子数接近的帧，然后进行提取，例如提取第4帧和12帧；
第六步，将提取的帧进行叠加变成一个帧；例如第4帧加第12帧；
第七步，对已叠加的帧(4+12)进行图像重建得到最后校正的图像。
用实验来评价验证该门控方法。实验分别采用了运动振幅为5厘米和2厘米的人体模来仿真呼吸运动。然后将其对应收集到的动态数据划分成相应的帧来验证所提出的门控方法。其图像重建结果如图4和图5所示。比较重建图像可以看出，图4(a)和图5(a)是静态图像，也是校正质量的比较标准。可以看出运动的体模在轴向方向(如虚线标记)边沿产生了很严重的模糊如图4(b)和图5(b)，在经过门控校正后的图像和模糊运动图像相比显示，在边沿区域得到了较明显的恢复校正，如4图(c)和图5(c)所示。这样就有利于对呼吸运动目标的大小、形状和位置进行有效的判断。
Claims (7)
1.一种基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征在于在三维PET成像过程中，首先利用动态模式获取成像的数据；然后根据不同的时间分辨率将数据划分为数据帧，每一帧看作一个呼吸循环中的一个运动相位；结合扫描器空间灵敏度的特性就可分析观察出每帧中光子数的周期性变化情况；依据所有帧内光子数的变化，用门控方法进行图像校正，最后获得校正的图像。
2.根据权利要求1所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征在于动态扫描数据帧内的光子数会反映出呼吸运动沿扫描探测器轴向变化的位置。
3.根据权利要求1所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征是对获取的动态模式数据划分为相同时间间隔的数据帧，其每一帧看作一个呼吸循环中的一个运动相位。
4.根据权利要求1所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征是对每一帧光子进行噪声滤波以凸显呼吸循环的特征。
5.根据权利要求4所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征是在已经进行噪声滤波的帧中，寻找那些光子数接近的帧，然后把它们进行叠加。
6.根据权利要求5所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征是对叠加的帧进行图像重建。
7.根据权利要求1所述的基于三维正电子发射探测器灵敏度特性的呼吸门控技术，其特征在于具体步骤是：
第一步，用PET的动态数据模式获取两个呼吸循环的扫描成像数据；
第二步，将这两个呼吸循环划分成为16个帧，每个循环包含8个帧；
第三步，计算每个帧内的光子数，并画出帧光子曲线图；
第四步，对帧光子曲线图进行噪声滤波以凸显呼吸循环的特征；
第五步，在两个呼吸循环中寻找光子数接近的帧，然后进行提取；
第六步，将提取的帧进行叠加变成一个帧；
第七步，对已叠加的帧进行图像重建，最后获得校正的图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100949595A CN101658426B (zh) | 2009-09-11 | 2009-09-11 | 基于三维正电子发射探测器灵敏度特性的呼吸门控技术 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100949595A CN101658426B (zh) | 2009-09-11 | 2009-09-11 | 基于三维正电子发射探测器灵敏度特性的呼吸门控技术 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101658426A CN101658426A (zh) | 2010-03-03 |
| CN101658426B true CN101658426B (zh) | 2011-08-31 |
Family
ID=41786839
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100949595A Expired - Fee Related CN101658426B (zh) | 2009-09-11 | 2009-09-11 | 基于三维正电子发射探测器灵敏度特性的呼吸门控技术 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101658426B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102151142B (zh) * | 2011-04-14 | 2012-08-15 | 华中科技大学 | 一种正电子发射断层成像中的运动门控方法及系统 |
| EP2724319B1 (en) | 2011-06-21 | 2020-03-18 | Koninklijke Philips N.V. | Respiratory motion determination apparatus |
| WO2014054379A1 (ja) * | 2012-10-04 | 2014-04-10 | コニカミノルタ株式会社 | 画像処理装置及びプログラム |
| CN104000618B (zh) * | 2014-05-13 | 2016-03-30 | 昆明理工大学 | 一种用环真光子数门控方法进行呼吸运动门控校正方法 |
| CN107095691B (zh) * | 2017-07-11 | 2020-07-28 | 上海联影医疗科技有限公司 | 一种pet成像方法及系统 |
| CN107468267B (zh) * | 2017-08-11 | 2020-12-04 | 上海联影医疗科技股份有限公司 | 数据处理方法及医学成像设备 |
| US10690782B2 (en) * | 2017-09-30 | 2020-06-23 | Uih America, Inc. | Systems and methods for positron emission tomography image reconstruction |
| CN110192886B (zh) * | 2019-05-29 | 2022-12-23 | 明峰医疗系统股份有限公司 | 一种动态调用呼吸门控的重建方法 |
| CN111643082A (zh) * | 2020-06-22 | 2020-09-11 | 苏州高新区人民医院 | 一种乳腺mri成像呼吸伪影抑制方法 |
| CN113971659B (zh) * | 2021-09-14 | 2022-08-26 | 杭州微引科技有限公司 | 一种经皮肺以及腹部穿刺的呼吸门控系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1415275A (zh) * | 2002-11-22 | 2003-05-07 | 赵耀德 | 数字化呼吸门控ct导航手术系统 |
| CN101305921A (zh) * | 2007-05-17 | 2008-11-19 | 上海西门子医疗器械有限公司 | 一种实现ct呼吸门控的方法及装置 |
| US7519413B1 (en) * | 2003-02-20 | 2009-04-14 | S.A. Instruments, Inc. | Apparatus and method for measuring motion in a strong magnetic field |
-
2009
- 2009-09-11 CN CN2009100949595A patent/CN101658426B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1415275A (zh) * | 2002-11-22 | 2003-05-07 | 赵耀德 | 数字化呼吸门控ct导航手术系统 |
| US7519413B1 (en) * | 2003-02-20 | 2009-04-14 | S.A. Instruments, Inc. | Apparatus and method for measuring motion in a strong magnetic field |
| CN101305921A (zh) * | 2007-05-17 | 2008-11-19 | 上海西门子医疗器械有限公司 | 一种实现ct呼吸门控的方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 樊卫.PET/CT呼吸门控技术及其应用介绍.《中华核医学杂志》.2007,第27卷(第2期),126-127. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101658426A (zh) | 2010-03-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101658426B (zh) | 基于三维正电子发射探测器灵敏度特性的呼吸门控技术 | |
| CN101702232B (zh) | 正电子发射成像中呼吸校正技术 | |
| KR101428005B1 (ko) | 소수의 저선량 ｃｔ 영상을 이용하여 ｐｅｔ 영상을 움직임 보상 및 감쇠 보정하는 방법 | |
| WO2013140749A1 (ja) | 画像解析装置および方法並びにプログラム | |
| US8532357B2 (en) | Method and apparatus for reducing image artifacts | |
| CN103608845A (zh) | 呼吸运动确定装置 | |
| WO2006137294A1 (ja) | Ｘ線診断支援装置、プログラム及び記録媒体 | |
| JP2010501856A (ja) | 動きアーティファクト画像についてのアーティファクト補正 | |
| US11823384B2 (en) | CT image generation method for attenuation correction of pet images | |
| CN107348969B (zh) | 一种pet数据处理方法、系统及pet成像设备 | |
| CN102247156B (zh) | 图像处理装置以及x射线ct装置 | |
| Visvikis et al. | Respiratory motion in positron emission tomography for oncology applications: Problems and solutions | |
| US20230130015A1 (en) | Methods and systems for computed tomography | |
| JP2006102353A (ja) | 関節動作解析装置、関節動作解析方法及び関節動作解析プログラム | |
| CN114387364A (zh) | 用于pet图像重建的线性衰减系数获取方法及重建方法 | |
| Daouk et al. | Respiratory-gated positron emission tomography and breath-hold computed tomography coupling to reduce the influence of respiratory motion: methodology and feasibility | |
| US20100030572A1 (en) | Temporal registration of medical data | |
| Tanaka et al. | Computerized methods for determining respiratory phase on dynamic chest radiographs obtained by a dynamic flat-panel detector (FPD) system | |
| CN104000618A (zh) | 一种用环真光子数门控方法进行呼吸运动门控校正技术 | |
| CN101077303B (zh) | 测定身体部位的正电子发射测量信息的方法及所属装置 | |
| Malczewski | PET image reconstruction using compressed sensing | |
| EP3707677B1 (en) | Attenuation correction of pet data of moving object | |
| CN112991478A (zh) | 基于深度学习的影像分析多时相异常特征区域参数的方法 | |
| Nye et al. | Current state of hybrid imaging: attenuation correction and fusion | |
| JP2010169542A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110831 Termination date: 20140911 |
|
| EXPY | Termination of patent right or utility model |