CN104299239B - 一种基于阿尔法散度的动态pet图像因子处理方法 - Google Patents
一种基于阿尔法散度的动态pet图像因子处理方法 Download PDFInfo
- Publication number
- CN104299239B CN104299239B CN201410570293.7A CN201410570293A CN104299239B CN 104299239 B CN104299239 B CN 104299239B CN 201410570293 A CN201410570293 A CN 201410570293A CN 104299239 B CN104299239 B CN 104299239B
- Authority
- CN
- China
- Prior art keywords
- factor
- image
- dynamic pet
- divergence
- dynamic
- 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.)
- Active
Links
- 238000011282 treatment Methods 0.000 title abstract description 3
- 238000000556 factor analysis Methods 0.000 claims abstract description 25
- 238000012879 PET imaging Methods 0.000 claims abstract description 11
- 238000012636 positron electron tomography Methods 0.000 claims abstract 18
- 238000000034 method Methods 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 16
- 238000003672 processing method Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 4
- 238000000354 decomposition reaction Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000002600 positron emission tomography Methods 0.000 description 35
- 239000008280 blood Substances 0.000 description 11
- 210000004369 blood Anatomy 0.000 description 11
- 238000012905 input function Methods 0.000 description 10
- 230000007547 defect Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000010191 image analysis Methods 0.000 description 3
- 102000002262 Thromboplastin Human genes 0.000 description 2
- 108010000499 Thromboplastin Proteins 0.000 description 2
- 238000010241 blood sampling Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000003703 image analysis method Methods 0.000 description 2
- 210000005240 left ventricle Anatomy 0.000 description 2
- 239000000700 radioactive tracer Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009509 drug development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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
- G06T2207/10104—Positron emission tomography [PET]
-
- 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
- G06T2207/30048—Heart; Cardiac
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
一种基于 散度的动态PET图像因子处理方法，包括如下步骤，（1）利用PET成像设备进行动态扫描并进行图像重建，得到动态PET图像；（2）建立因子分析模型；（3）定义动态PET图像与因子分析模型的散度；（4）最小化散度得到初级因子图像和对应的初级因子；（5）最小化初级因子图像间的重叠程度，得到最终因子图像和对应的最终因子。本发明利用了动态PET图像与因子分析模型的散度测度，可以根据动态PET图像噪声不同分布特性，进行因子分析，得到相应的因子图像及因子。由于对分解得到的因子图像加入了唯一性约束使得得到的因子图像之间重叠度最小，可以有效提高因子图像的准确性。
Description
技术领域
本发明涉及一种医学影像的图像分析方法，特别涉及一种基于阿尔法（）散度的动态PET图像因子处理方法。
背景技术
正电子发射计算机断层成像（Positron Emission Tomography, PET）是医学影像领域最先进的临床检查技术，可以在分子水平对疾病进行诊断。动态PET成像，不仅能显示放射性示踪剂浓度的空间分布，还揭示了示踪剂随时间代谢的动态过程。通过应用动力学模型，即可获得具有丰富意义的动力学参数，从而为疾病诊断、治疗以及药物开发提供指导。
在动力学模型的应用中，血输入函数的估计是非常关键的。传统方法通过连续的血液采样估计血输入函数，然而长时间的血液采样经常会导致病人不适。在动态PET成像技术中，亟待无创性、高效准确的方法来获得血输入函数。
目前广泛采用的估计血输入函数方法为基于感兴趣区域法，此方法通过有经验人士在动态PET图像上勾画出感兴趣区域（如左心室），从而得到相应的时间活度曲线作为血输入函数。感兴趣区域法简单易行，但该方法存在一定的缺陷，如不精确性，易受部分容积效应影响等。1982年，Di Paola等首次提出在动态序列图像中利用因子分析模型提取组织的时间活度曲线。该方法假设各组织的时间活度曲线为因子，在对动态序列图像进行主成分分析后非负约束下旋转投影轴，进而获得因子图像和因子。1999年，Attias等提出基于独立成分分析的方法提取组织的时间活度曲线。2000年，Sitek等在因子分析模型的基础上引入先验信息，并将其应用于动态SPECT图像。实验结果表明该方法获得的组织时间活度曲线与真实值较为一致。随后，Sitek 与Fakhri等提出用最小二乘法对因子分析模型进行求解，并提出唯一性约束最小化各因子图像间的重叠程度。该方法很好地解决了因子分析模型中非唯一性解的问题。2007年，Yi Su等提出最大似然因子分析法（MLFA），并基于唯一性约束将其应用于心脏动态PET图像数据上。
然而，感兴趣区域法极度依赖于医生手工勾画的感兴趣区域，导致感兴趣区域的准确性受医生个人经验及部分容积效应影响。导致动态PET成像中血输入函数估计的准确度有限。
因此，针对现有技术不足，提供一种基于散度的动态PET图像因子处理方法以克服现有技术不足甚为必要。
发明内容
本发明的目的在于避免现有技术的不足之处而提供一种基于阿尔法散度的动态PET图像因子处理方法，可以准确得到各组织因子图像及对应因子，为动态PET成像中血输入函数估计提供技术支持。
本发明的上述目的通过如下技术手段实现。
提供一种基于阿尔法散度的动态PET图像因子处理方法，包括如下步骤，
（1）利用PET成像设备进行动态扫描并进行图像重建，得到动态PET图像；
（2）建立因子分析模型；
（3）定义动态PET图像与因子分析模型的散度；
（4）最小化散度得到初级因子图像和对应的初级因子；
（5）最小化初级因子图像间的重叠程度，得到最终因子图像和对应的最终因子。
上述步骤（2）中建立的因子分析模型为线性模型，具体是：
假设动态PET图像I上每个像素点的像素值可由各个因子F的线性组合表示，每个因子的系数由因子图像L确定，动态PET图像I可表示为：
；
其中因子为组织的时间活度曲线，动态PET图像的大小为，其中表示每一帧PET图像的像素点数目，表示动态PET图像的帧数，因子矩阵的大小为，因子图像的大小为，其中为所需提取的因子数目。
上述步骤（3）中定义的动态PET图像与因子分析模型的散度为：
；
其中 ，表示在当前因子图像和当前因子下所得到估计值；
为任意值，且当=1时，； 当=-1时，。
上述步骤（4）具体是：非负约束下最小化散度，令以及，计算得到的因子、因子图像具体为：
；
；
其中，表示迭代步数，根据以上迭代公式依次对因子、因子图像进行更新，以更新后的结果作为初级因子图像L’和对应的初级因子F’。
上述步骤（5）中最小化初级因子图像间的重叠程度，具体是：
各初级因子图像间的重叠程度可用以下式子表示：
；
用梯度下降法最小化 ，求得最终因子图像L 0 和对应的最终因子F 0 。
上述步骤（1）具体是采用滤波反投影方法进行图像重建，得到动态PET图像的。
本发明的基于阿尔法散度的动态PET图像因子处理方法，包括如下步骤，（1）利用PET成像设备进行动态扫描并进行图像重建，得到动态PET图像；（2）建立因子分析模型；（3）定义动态PET图像与因子分析模型的散度；（4）最小化散度得到初级因子图像和对应的初级因子；（5）最小化初级因子图像间的重叠程度，得到最终因子图像和对应的最终因子。本发明利用了动态PET图像与因子分析模型的散度测度，可以根据动态PET图像噪声不同分布特性，进行因子分析，得到相应的因子图像及因子。由于对分解得到的因子图像加入了唯一性约束使得得到的因子图像之间重叠度最小，可以有效提高因子图像的准确性从而能够准确得到各组织因子图像及对应的因子。
附图说明
利用附图对本发明作进一步的说明，但附图中的内容不构成对本发明的任何限制。
图1是本发明一种基于阿尔法散度的动态PET图像因子处理方法的流程图。
图2（a）是本发明实施例2使用的仿真的体模图像; 图2（b ）是本发明实施例2使用的仿真的体模图像对应的时间活度曲线。
图3（a）是时得到的最终因子图像; 图3（b ）是时得到的时间活度曲线。
图4（a）是时得到的最终因子图像; 图4（b）是时得到的时间活度曲线。
图5（a）是时得到的最终因子图像; 图5（b）是时得到的时间活度曲线。
图6（a）是时得到的最终因子图像; 图6（b）是时得到的时间活度曲线。
图7是在不同取值情况下血输入函数的均方误差表。
具体实施方式
结合以下实施例对本发明作进一步描述。
实施例1。
一种基于散度的动态PET图像因子处理方法，如图1所示，包括如下步骤。
（1）利用PET成像设备进行动态扫描并进行图像重建，得到动态PET图像。具体可以采用滤波反投影方法进行图像重建，得到动态PET图像。
（2）建立因子分析模型。
步骤（2）中建立的因子分析模型为线性模型，具体是：
假设动态PET图像I上每个像素点的像素值可由各个因子F的线性组合表示，每个因子的系数由因子图像L确定，动态PET图像I可表示为：
；
其中因子为组织的时间活度曲线，动态PET图像的大小为，其中表示每一帧PET图像的像素点数目，表示动态PET图像的帧数，因子矩阵的大小为，因子图像的大小为，其中为所需提取的因子数目。
（3）定义动态PET图像与因子分析模型的散度。
步骤（3）中定义的动态PET图像与因子分析模型的散度为：
；
其中，表示在当前因子图像L和当前因子F下所得到估计值；
为任意值，且当=1时，；当=-1时，；当=2时，。
（4）最小化散度得到初级因子图像和对应的初级因子。
步骤（4）具体是：非负约束下最小化散度，令以及，计算得到的因子、因子图像具体为：
；
；
其中，表示迭代步数，根据以上迭代公式依次对因子、因子图像进行更新，以更新后的结果作为初级因子图像L’和对应的初级因子F’。
（5）最小化初级因子图像间的重叠程度，得到最终因子图像和对应的最终因子。
步骤（5）中最小化初级因子图像间的重叠程度，具体是：
各初级因子图像间的重叠程度可用以下式子表示：；
；
用梯度下降法最小化 ，求得最终因子图像L 0 和对应的最终因子F 0 。
本发明的基于散度的动态PET图像因子处理方法，利用了动态PET图像与因子分析模型的散度测度，可以根据动态PET图像噪声不同分布特性，进行因子分析，得到相应的因子图像及因子。由于对分解得到的因子图像加入了唯一性约束使得得到的因子图像之间重叠度最小，可以有效提高因子图像的准确性。
实施例2。
为了验证本发明技术方案的效果，采用仿真体模进行验证。
采用图2所示的体模进行，该体模包含两个组织成分，第一个组织成分是位于图像中心的圆圈，该组织成分由血液组织（左心室）组成。第二个组织成分是围绕圆圈的圆环，该组织成分由心肌组织组成。
对体模进行2 min的动态扫描，每帧间隔12秒，共10帧。重建得到的体模动态序列图像数据为，图像像素点大小为。在仿真的动态PET图像上加入高斯噪声模拟现实中的情况，并用FWHM=12的高斯滤波器进行滤波。
图3（a）、（b）至图6（a）、（b）分别是为1、-1、2和1.5时得到的最终因子图像及相应的最终因子。从图3（a）、（b）至图6（a）、（b）的结果可以看出，本发明的方法所得到的最终因子图像准确性高，所得到对应的最终因子也具有准确性高的特点。
图7是在不同取值情况下血输入函数的均方差表。从中可以看出，血输入函数的均方差在不同取值情况下均较小。
最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非对本发明保护范围的限制，尽管参照较佳实施例对本发明作了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。
Claims (4)
1.一种基于阿尔法散度的动态PET图像因子处理方法，其特征在于：包括如下步骤，
(1)利用PET成像设备进行动态扫描并图像重建，得到动态PET图像；
(2)建立因子分析模型；
(3)定义动态PET图像与因子分析模型的α散度；
(4)最小化α散度得到初级因子图像和对应的初级因子；
(5)最小化初级因子图像间的重叠程度，得到最终因子图像和对应的最终因子；
所述步骤(2)中建立的因子分析模型为线性模型，具体是：
假设动态PET图像I上每个像素点的像素值可由各个因子F的线性组合表示，每个因子的系数由因子图像L确定，动态PET图像I可表示为：
I＝LF；
其中因子F为组织的时间活度曲线，动态PET图像I的大小为N×T，其中N表示每一帧PET图像的像素点数目，T表示动态PET图像的帧数，因子矩阵F的大小为K×T，因子图像L的大小为N×K，其中K为所需提取的因子数目；
所述步骤(3)中定义的动态PET图像与因子分析模型的α散度Dα为：
其中i＝1...N，f＝1...T，表示在当前因子图像L和当前因子F下所得到估计值；
α为任意值，且当α＝1时，当α＝-1时，
2.根据权利要求1所述的基于阿尔法散度的动态PET图像因子处理方法，其特征在于：
所述步骤(4)具体是：非负约束下最小化α散度，令以及 计算得到的因子Fkf、因子图像Lik具体为：
其中，m表示迭代步数，根据以上迭代公式依次对因子Fkf、因子图像Lik进行更新，以更新后的结果作为初级因子图像L’和对应的初级因子F’。
3.根据权利要求2所述的基于阿尔法散度的动态PET图像因子处理方法，其特征在于：
所述步骤(5)中最小化初级因子图像间的重叠程度，具体是：
各初级因子图像间的重叠程度可用以下式子表示：
用梯度下降法最小化funi(L’)，求得最终因子图像L0和对应的最终因子F0。
4.根据权利要求3所述的基于阿尔法散度的动态PET图像因子处理方法，其特征在于：所述步骤(1)具体是采用滤波反投影方法进行图像重建，得到动态PET图像的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410570293.7A CN104299239B (zh) | 2014-10-23 | 2014-10-23 | 一种基于阿尔法散度的动态pet图像因子处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410570293.7A CN104299239B (zh) | 2014-10-23 | 2014-10-23 | 一种基于阿尔法散度的动态pet图像因子处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104299239A CN104299239A (zh) | 2015-01-21 |
| CN104299239B true CN104299239B (zh) | 2017-08-08 |
Family
ID=52318961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410570293.7A Active CN104299239B (zh) | 2014-10-23 | 2014-10-23 | 一种基于阿尔法散度的动态pet图像因子处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104299239B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108765318A (zh) * | 2018-05-15 | 2018-11-06 | 南方医科大学 | 一种基于动力学聚类的动态pet图像因子处理方法 |
| CN112052882B (zh) * | 2020-08-14 | 2023-08-22 | 北京师范大学 | 磁共振脑结构影像的分类模型构建、分类与可视化方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102314698A (zh) * | 2011-08-10 | 2012-01-11 | 南方医科大学 | 基于阿尔法散度约束的全变分最小化剂量ct重建方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7653264B2 (en) * | 2005-03-04 | 2010-01-26 | The Regents Of The University Of Michigan | Method of determining alignment of images in high dimensional feature space |
| WO2013116709A1 (en) * | 2012-02-01 | 2013-08-08 | The Research Foundation of States University of New York | Computerized image reconstruction method and apparatus |
-
2014
- 2014-10-23 CN CN201410570293.7A patent/CN104299239B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102314698A (zh) * | 2011-08-10 | 2012-01-11 | 南方医科大学 | 基于阿尔法散度约束的全变分最小化剂量ct重建方法 |
Non-Patent Citations (4)
| Title |
|---|
| -divergence for PET imaging.《Computerized Medical Imaging and Graphics》.2011,第35卷第294-301页. * |
| Information Divergence Constrained Total Variation Minimization for Positron Emission Tomography Image Reconstruction;Lingling Tian, et al.;《2011 IEEE Nuclear Science Symposium Conference Record》;20111231;第2587-2592页 * |
| Yueyang Teng, et al..Iterative reconstruction algorithms with ˛ * |
| 基于全变分α散度最小化的PET优质重建;田玲玲 等;《电子学报》;20120630;第40卷(第6期);第1263-1268页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104299239A (zh) | 2015-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109009200B (zh) | 用于正电子发射断层成像图像重建的系统和方法 | |
| JP5959699B2 (ja) | パノラマエックス線写真を利用した診断支援システム、及びパノラマエックス線写真を利用した診断支援プログラム | |
| CN103136773B (zh) | 一种稀疏角度x射线ct成像方法 | |
| CN112435307B (zh) | 一种深度神经网络辅助的四维锥束ct影像重建方法 | |
| JP7324195B2 (ja) | ディープイメージングを使用したポジトロンエミッショントモグラフィシステム設計の最適化 | |
| TW201444534A (zh) | 動態單光子放射電腦斷層攝影或是單光子放射電腦斷層攝影／電腦斷層攝影成像之心肌血流量化系統及方法 | |
| KR20140055896A (ko) | 소수의 저선량 ｃｔ 영상을 이용하여 ｐｅｔ 영상을 움직임 보상 및 감쇠 보정하는 방법 | |
| US8655040B2 (en) | Integrated image registration and motion estimation for medical imaging applications | |
| WO2014172927A1 (zh) | 一种基于rpca的pet图像动态重建方法及系统 | |
| CN103186882A (zh) | Pet系统中图像的衰减校正方法及装置 | |
| EP2147409A1 (en) | Model-based spect heart orientation estimation | |
| CN108596861A (zh) | 一种基于深度残差学习的ct图像金属伪影去除方法 | |
| CN109961419A (zh) | 对pet活度分布图像进行衰减校正的校正信息获取方法 | |
| Yoder | Basic PET data analysis techniques | |
| CN104299239B (zh) | 一种基于阿尔法散度的动态pet图像因子处理方法 | |
| Aladl et al. | Automated image registration of gated cardiac single‐photon emission computed tomography and magnetic resonance imaging | |
| Izadi et al. | Enhanced direct joint attenuation and scatter correction of whole-body PET images via context-aware deep networks | |
| WO2014122236A1 (en) | Synchronous physiological measurements for cardiac acquisitions | |
| CN105279777B (zh) | 基于改进序贯滤波的静态pet图像重建方法 | |
| Jiao et al. | Spatial-temporal pharmacokinetic model based registration of 4d brain pet data | |
| Barna et al. | A direct method for reconstructing dynamic SPECT images | |
| Zhang et al. | Self-supervised Vessel Segmentation from X-ray Images using Digitally Reconstructed Radiographs | |
| Yang et al. | Assessment of Spatial Resolution in PET Images Denoised with Deep Convolutional Neural Networks | |
| US20230045406A1 (en) | System and method for hybrid imaging | |
| Xie et al. | Noise-aware Dynamic Image Denoising and Positron Range Correction for Rubidium-82 Cardiac PET Imaging via Self-supervision |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |