CN102764135A - 心肌spect灌注显像的心脏轴向移动校正方法 - Google Patents
心肌spect灌注显像的心脏轴向移动校正方法 Download PDFInfo
- Publication number
- CN102764135A CN102764135A CN2012101390513A CN201210139051A CN102764135A CN 102764135 A CN102764135 A CN 102764135A CN 2012101390513 A CN2012101390513 A CN 2012101390513A CN 201210139051 A CN201210139051 A CN 201210139051A CN 102764135 A CN102764135 A CN 102764135A
- Authority
- CN
- China
- Prior art keywords
- projection
- myocardium
- axial
- perfusion imaging
- image
- 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 28
- 230000010412 perfusion Effects 0.000 title claims abstract description 18
- 238000003384 imaging method Methods 0.000 title claims abstract description 17
- 238000002603 single-photon emission computed tomography Methods 0.000 title claims abstract description 16
- 210000004165 myocardium Anatomy 0.000 title claims 10
- 230000033001 locomotion Effects 0.000 title abstract description 14
- 238000006073 displacement reaction Methods 0.000 claims abstract description 43
- 238000012937 correction Methods 0.000 claims abstract description 23
- 230000000747 cardiac effect Effects 0.000 claims abstract description 11
- 238000005314 correlation function Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims 3
- 230000002107 myocardial effect Effects 0.000 abstract description 9
- 238000013480 data collection Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000009278 visceral effect Effects 0.000 description 1
- 230000002747 voluntary effect Effects 0.000 description 1
Images
Landscapes
- Nuclear Medicine (AREA)
Abstract
一种心肌SPECT灌注显像的心脏轴向移动校正方法，包括以下步骤：(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；(3)使用经初步轴向位移校正的投影数据，重建三维图像；(4)将重建图像根据实际投影的采集条件进行再投影；(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；(7)使用经上一步轴向位移校正的投影数据，重建三维图像；(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。
Description
技术领域
本发明涉及核医学影像仪器领域，特别适用于使用单光子计算机断层装置进行心肌灌注显像时对心脏轴向移动的校正方法。
背景技术
由于SPECT心肌灌注显像采集的时间较长，采集过程中受检者身体移动、呼吸、内脏的运动等自主或不自主的因素都会造成心脏的移动，使重建结果中出现移动伪影，影响诊断结果。据统计，心肌灌注断层显像采集过程中，约25％的受检者会发生不同程度的体位移动。心脏的移动可能导致心肌血流灌注影像出现缺损而导致假阳性，造成误诊。
发明内容
本发明所要解决的技术问题是提供一种心肌SPECT灌注显像的心脏轴向移动校正方法，以最大限度地降低运动伪影对正确诊断的影响。
为解决上述技术问题本发明的技术方案是：一种心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于，包括以下步骤：
(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；
(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；
(3)使用经初步轴向位移校正的投影数据，重建三维图像；
(4)将重建图像根据实际投影的采集条件进行再投影；
(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；
(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；
(7)使用经上一步轴向位移校正的投影数据，重建三维图像；
(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。
作为改进，所述步骤(1)具体为：首先将各个角度的投影图像沿横向求和，然后对这个一维序列进行互相关分析；一维序列数据Pn与Pn-1的互相关函数可表示为： 其中：d代表偏移，M代表每组一维数据中数据的总个数；-k≤d≤k，k为最大可能的位移，一般设为10个像素；当j+d＜1或j+d＞M时，令Pn-1(j+d)＝0；互相关函数CCn(d)表示Pn与偏移了d的Pn-1的相互匹配程度，CCn(d)越大，匹配得越好；因此对CCn(d)进行抛物线拟合，求出抛物线顶点对应的d，即可计算出位移。
作为改进，所述步骤(2)具体为：通过线性插值的方法，实现非像素大小整数倍的图像的轴向位移，插值后对图像的总计数进行归一，即每个像素的计数乘以原始投影的总计数后，除以插值后图像的总计数。
作为改进，所述步骤(3)具体为：图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正。
作为改进，所述步骤(4)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。
作为改进，所述步骤(5)具体为：对每个角度投影的计算值与测量值根据互信息进行图像配准，计算出各个角度投影沿轴向位移的距离。
作为改进，所述步骤(6)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。
作为改进，所述步骤(7)具体为：图像重建采用滤波反投影法，使用But terworth滤波器重建图像，不加衰减校正。
作为改进，所述步骤(8)具体为：重复第(4)-第(5)步，比较两次校正值，当对应校正值最大差别小于0.5mm时，校正完成。
本发明与现有技术相比所带来的有益效果是：
本发明的优点是最大限度地降低了再投影的图像中位移的影响，使基于再投影的校正为有根据的校正。以前的投影再投影方法，直接对有位移的投影进行重建，重建结果是基于错误投影的，以此重建结果的再投影做为校正的标准显然是没有根据的。本发明的优点是拓宽了再投影方法的应用范围，提高了运动校正的精度。
附图说明
图1为实施例1从第8帧开始加偏移，检测位移与设定位移的对比；
图2为实施例1从第12帧开始加偏移，检测位移与设定位移的对比。
具体实施方式
下面结合说明书附图对本发明作进一步说明。
实施例1
对单次移动的校正
本实施例用来模拟采集过程中从某帧开始以后的投影均有同样的位移的情况。
投影的采集：
显像仪器SPECT，在心肌模型中注射99Tcm1mCi，注射30min后进行心肌灌注显像，探头旋转从右前斜45°到左后斜45°，帧间角度间隔6°，每个探头16帧，共采集32帧。选用低能高分辨准直器，采集矩阵为64×64，放大倍数为1.46，每帧采集20s。
数据采集过程中，手工移动模型，通过固定在检查床上的尺子来确定实际移动的大小。
数据处理：
(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；
首先将各个角度的投影图像沿横向求和，然后对这个一维序列进行互相关分析；
一维序列数据Pn与Pn-1的互相关函数可表示为：
其中：d代表偏移，M代表每组一维数据中数据的总个数；-k≤d≤k，k为最大可能的位移，一般设为10个像素；当j+d＜1或j+d＞M时，令Pn-1(j+d)＝0；
互相关函数CCn(d)表示Pn与偏移了d的Pn-1的相互匹配程度，CCn(d)越大，匹配得越好；
因此对CCn(d)进行抛物线拟合，求出抛物线顶点对应的d，即可计算出位移；
(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；
通过线性插值的方法，实现非像素大小整数倍的图像的轴向位移，插值后对图像的总计数进行归一，即每个像素的计数乘以原始投影的总计数后，除以插值后图像的总计数；
(3)使用经初步轴向位移校正的投影数据，重建三维图像；
图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正；
(4)将重建图像根据实际投影的采集条件进行再投影；
通过拉登变换计算实际采集条件的、上一步重建图像的投影；
(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；
对每个角度投影的计算值与测量值根据互信息进行图像配准，计算出各个角度投影沿轴向位移的距离；
(6)根据上一步的轴向位移值，对各个投影进行校正；
通过拉登变换计算实际采集条件的、上一步重建图像的投影；
(7)使用经上一步轴向位移校正的投影数据，重建三维图像；
图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正；
(8)重复第(4)-第(5)步，比较两次校正值，当对应校正值最大差别小于0.5mm时，校正完成。否则，重复第(6)步-第(8)步。
如图1、2所示，综合所有模拟结果，本方法对轴向移动的校正的误差小于3mm。
实施例2
与实施例1所不同的是，步骤(1)、(2)，实施过程中采用手工移动的方法，对投影的位移进行校正。
实施例3
与实施例1所不同的是，步骤(8)，采用其他标准停止迭代的过程，如比较本次校正的投影与再投影的差异，如各个像素的均方差，利用制定的差异标准停止迭代。
Claims (9)
1.一种心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于，包括以下步骤：
(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；
(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；
(3)使用经初步轴向位移校正的投影数据，重建三维图像；
(4)将重建图像根据实际投影的采集条件进行再投影；
(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；
(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；
(7)使用经上一步轴向位移校正的投影数据，重建三维图像；
(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。
2.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(1)具体为：首先将各个角度的投影图像沿横向求和，然后对这个一维序列进行互相关分析；
一维序列数据Pn与Pn-1的互相关函数可表示为：
其中：d代表偏移，M代表每组一维数据中数据的总个数；-k≤d≤k，k为最大可能的位移，一般设为10个像素；当j+d＜1或j+d＞M时，令Pn-1(j+d)＝0；互相关函数CCn(d)表示Pn与偏移了d的Pn-1的相互匹配程度，CCn(d)越大，匹配得越好；因此对CCn(d)进行抛物线拟合，求出抛物线顶点对应的d，即可计算出位移。
3.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(2)具体为：通过线性插值的方法，实现非像素大小整数倍的图像的轴向位移，插值后对图像的总计数进行归一，即每个像素的计数乘以原始投影的总计数后，除以插值后图像的总计数。
4.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(3)具体为：图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正。
5.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(4)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。
6.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(5)具体为：对每个角度投影的计算值与测量值根据互信息进行图像配准，计算出各个角度投影沿轴向位移的距离。
7.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(6)具体为：通过拉登变换计算实际采集条件的、上一步重建图像的投影。
8.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(7)具体为：图像重建采用滤波反投影法，使用Butterworth滤波器重建图像，不加衰减校正。
9.根据权利要求1所述的心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于：所述步骤(8)具体为：重复第(4)-第(5)步，比较两次校正值，当对应校正值最大差别小于0.5mm时，校正完成。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210139051.3A CN102764135B (zh) | 2012-04-23 | 2012-04-23 | 心肌spect灌注显像的心脏轴向移动校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210139051.3A CN102764135B (zh) | 2012-04-23 | 2012-04-23 | 心肌spect灌注显像的心脏轴向移动校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102764135A true CN102764135A (zh) | 2012-11-07 |
| CN102764135B CN102764135B (zh) | 2017-02-22 |
Family
ID=47091879
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210139051.3A Expired - Fee Related CN102764135B (zh) | 2012-04-23 | 2012-04-23 | 心肌spect灌注显像的心脏轴向移动校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102764135B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105997125A (zh) * | 2016-06-15 | 2016-10-12 | 刘丽 | 多针孔单光子spect心肌血流绝对定量方法与用途 |
| TWI611795B (zh) * | 2016-07-18 | 2018-01-21 | 許百靈 | 多針孔單光子ｓｐｅｃｔ心肌血流絕對定量方法與用途 |
| CN112785588A (zh) * | 2021-02-05 | 2021-05-11 | 南京钺曦医疗科技有限公司 | 一种ct与mr脑灌注数据的运动幅度自动估计方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060000983A1 (en) * | 2004-06-18 | 2006-01-05 | Martin Charron | Fast dynamic imaging protocol using a multi-head single photon emission computed tomography system |
| US7332722B1 (en) * | 2006-02-21 | 2008-02-19 | Jefferson Science Associates, Llc | Simultaneous multi-headed imager geometry calibration method |
| US20080087829A1 (en) * | 2006-04-11 | 2008-04-17 | John Hoppin | Single-photon emission computed tomography (SPECT) using helical scanning with multiplexing multi-pinhole apertures |
| CN101449980A (zh) * | 2008-12-30 | 2009-06-10 | 北京海思威科技有限公司 | 单光子断层成像方法和单光子断层成像系统 |
| CN101512602A (zh) * | 2006-08-29 | 2009-08-19 | 皇家飞利浦电子股份有限公司 | 胸部ct成像中心脏运动伪影的减少 |
-
2012
- 2012-04-23 CN CN201210139051.3A patent/CN102764135B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060000983A1 (en) * | 2004-06-18 | 2006-01-05 | Martin Charron | Fast dynamic imaging protocol using a multi-head single photon emission computed tomography system |
| US7332722B1 (en) * | 2006-02-21 | 2008-02-19 | Jefferson Science Associates, Llc | Simultaneous multi-headed imager geometry calibration method |
| US20080087829A1 (en) * | 2006-04-11 | 2008-04-17 | John Hoppin | Single-photon emission computed tomography (SPECT) using helical scanning with multiplexing multi-pinhole apertures |
| CN101512602A (zh) * | 2006-08-29 | 2009-08-19 | 皇家飞利浦电子股份有限公司 | 胸部ct成像中心脏运动伪影的减少 |
| CN101449980A (zh) * | 2008-12-30 | 2009-06-10 | 北京海思威科技有限公司 | 单光子断层成像方法和单光子断层成像系统 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105997125A (zh) * | 2016-06-15 | 2016-10-12 | 刘丽 | 多针孔单光子spect心肌血流绝对定量方法与用途 |
| CN105997125B (zh) * | 2016-06-15 | 2021-09-17 | 刘丽 | 多针孔单光子spect心肌血流绝对定量方法与用途 |
| TWI611795B (zh) * | 2016-07-18 | 2018-01-21 | 許百靈 | 多針孔單光子ｓｐｅｃｔ心肌血流絕對定量方法與用途 |
| CN112785588A (zh) * | 2021-02-05 | 2021-05-11 | 南京钺曦医疗科技有限公司 | 一种ct与mr脑灌注数据的运动幅度自动估计方法 |
| CN112785588B (zh) * | 2021-02-05 | 2021-09-14 | 南京钺曦医疗科技有限公司 | 一种ct与mr脑灌注数据的运动幅度自动估计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102764135B (zh) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103136773B (zh) | 一种稀疏角度x射线ct成像方法 | |
| CN103679642B (zh) | 一种ct图像金属伪影校正方法、装置及ct设备 | |
| CN102737392B (zh) | 一种低剂量x线ct图像的非局部正则化先验重建方法 | |
| CN104574416B (zh) | 一种低剂量能谱ct图像去噪方法 | |
| CN104408758A (zh) | 一种低剂量能谱ct图像处理方法 | |
| Kyme et al. | Practical aspects of a data-driven motion correction approach for brain SPECT | |
| CN105147312A (zh) | Pet图像获取方法及系统 | |
| CN111260748B (zh) | 基于神经网络的数字合成x射线层析成像方法 | |
| CN103106676B (zh) | 一种基于低剂量投影数据滤波的x射线ct图像重建方法 | |
| Chan et al. | End-expiration respiratory gating for a high-resolution stationary cardiac SPECT system | |
| CN102831627A (zh) | 一种基于gpu多核并行处理的pet图像重建方法 | |
| CN102156974B (zh) | 解剖信息约束下基于h∞滤波的动态pet浓度重建方法 | |
| Sun et al. | A motion correction approach for oral and maxillofacial cone-beam CT imaging | |
| CN103186882A (zh) | Pet系统中图像的衰减校正方法及装置 | |
| McQuaid et al. | A novel method for incorporating respiratory-matched attenuation correction in the motion correction of cardiac PET–CT studies | |
| Mukherjee et al. | Improved frame‐based estimation of head motion in PET brain imaging | |
| CN110211193B (zh) | 三维ct层间图像插值修复与超分辨处理方法及装置 | |
| CN109512448A (zh) | 一种ct探测器位置校准方法及其系统 | |
| CN103325141A (zh) | 基于非等中心c形臂2d投影图像的3d模型构建方法 | |
| Friot et al. | Iterative tomographic reconstruction with TV prior for low-dose CBCT dental imaging | |
| Ebner et al. | Super‐resolution for upper abdominal MRI: Acquisition and post‐processing protocol optimization using brain MRI control data and expert reader validation | |
| CN102764135B (zh) | 心肌spect灌注显像的心脏轴向移动校正方法 | |
| JP7238134B2 (ja) | Ｐｅｔ画像化中の自動的な動き補正 | |
| Cui et al. | Prior knowledge-guided triple-domain transformer-GAN for direct PET reconstruction from low-count sinograms | |
| Thies et al. | A gradient-based approach to fast and accurate head motion compensation in cone-beam 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 | ||
| C53 | Correction of patent of invention or patent application | ||
| CB02 | Change of applicant information |
Address after: Lin Xi Lu Dragon Dragon Street Industrial Park in Longgang District of Shenzhen City, Guangdong Province, No. three building 518000 Applicant after: SHENZHEN BASDA MEDICAL APPARATUS CO.,LTD. Address before: Lin Xi Lu Dragon Dragon Street Industrial Park in Longgang District of Shenzhen City, Guangdong Province, No. three building 518000 Applicant before: SHENZHEN BASDA MEDICAL APPARATUS Co.,Ltd. |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: SHENZHEN BASDA MEDICAL APPARATUS CO., LTD. TO: SHENZHEN BASDA MEDICAL APPARATUS CO., LTD. |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: 518000 workshop 1a1901, building a, beisida medical equipment building, 28 Nantong Avenue, Baolong community, Baolong street, Longgang District, Shenzhen City, Guangdong Province Patentee after: SHENZHEN BASDA MEDICAL APPARATUS Co.,Ltd. Address before: Lin Xi Lu Dragon Dragon Street Industrial Park in Longgang District of Shenzhen City, Guangdong Province, No. three building 518000 Patentee before: SHENZHEN BASDA MEDICAL APPARATUS Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170222 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |