CN102764135A - 心肌spect灌注显像的心脏轴向移动校正方法 - Google Patents

心肌spect灌注显像的心脏轴向移动校正方法 Download PDF

Info

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
axial
myocardium
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
Application number
CN2012101390513A
Other languages
English (en)
Other versions
CN102764135B (zh
Inventor
曹文田
陈文波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHENZHEN BASDA MEDICAL APPARATUS CO Ltd
Original Assignee
SHENZHEN BASDA MEDICAL APPARATUS CO Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by SHENZHEN BASDA MEDICAL APPARATUS CO Ltd filed Critical SHENZHEN BASDA MEDICAL APPARATUS CO Ltd
Priority to CN201210139051.3A priority Critical patent/CN102764135B/zh
Publication of CN102764135A publication Critical patent/CN102764135A/zh
Application granted granted Critical
Publication of CN102764135B publication Critical patent/CN102764135B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Nuclear Medicine (AREA)

Abstract

一种心肌SPECT灌注显像的心脏轴向移动校正方法，包括以下步骤：(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；(3)使用经初步轴向位移校正的投影数据，重建三维图像；(4)将重建图像根据实际投影的采集条件进行再投影；(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；(7)使用经上一步轴向位移校正的投影数据，重建三维图像；(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。

Description

心肌SPECT灌注显像的心脏轴向移动校正方法
技术领域
本发明涉及核医学影像仪器领域，特别适用于使用单光子计算机断层装置进行心肌灌注显像时对心脏轴向移动的校正方法。
背景技术
由于SPECT心肌灌注显像采集的时间较长，采集过程中受检者身体移动、呼吸、内脏的运动等自主或不自主的因素都会造成心脏的移动，使重建结果中出现移动伪影，影响诊断结果。据统计，心肌灌注断层显像采集过程中，约25％的受检者会发生不同程度的体位移动。心脏的移动可能导致心肌血流灌注影像出现缺损而导致假阳性，造成误诊。
发明内容
本发明所要解决的技术问题是提供一种心肌SPECT灌注显像的心脏轴向移动校正方法，以最大限度地降低运动伪影对正确诊断的影响。
为解决上述技术问题本发明的技术方案是：一种心肌SPECT灌注显像的心脏轴向移动校正方法，其特征在于，包括以下步骤：
(1)投影数据采集完成后，初步估计各个投影轴向位移的大小；
(2)根据各个投影的轴向位移初步估计值，对各个投影进行初步校正；
(3)使用经初步轴向位移校正的投影数据，重建三维图像；
(4)将重建图像根据实际投影的采集条件进行再投影；
(5)对原始投影和对应的再投影进行分析，计算出实际投影的轴向位移；
(6)根据上一步的实际投影的轴向位移值，对各个投影进行校正；
(7)使用经上一步轴向位移校正的投影数据，重建三维图像；
(8)重复第(4)-第(5)步，当对应校正值最大差别小于要求值时，校正完成；否则，重复第(6)步-第(8)步。
作为改进，所述步骤(1)具体为：首先将各个角度的投影图像沿横向求和，然后对这个一维序列进行互相关分析；一维序列数据Pn与Pn-1的互相关函数可表示为： CC n ( d ) = Σ j = 1 M P n ( j ) * P n - 1 ( j + d ) , 其中：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的互相关函数可表示为：
CC n ( d ) = Σ j = 1 M P n ( j ) * P n - 1 ( j + d )
其中：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的互相关函数可表示为： CC n ( d ) = Σ j = 1 M P n ( j ) * P n - 1 ( j + d )
其中：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时，校正完成。
CN201210139051.3A 2012-04-23 2012-04-23 心肌spect灌注显像的心脏轴向移动校正方法 Active CN102764135B (zh)

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 Active CN102764135B (zh) 2012-04-23 2012-04-23 心肌spect灌注显像的心脏轴向移动校正方法

Country Status (1)

Country Link
CN (1) CN102764135B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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成像中心脏运动伪影的减少

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
AU2009260422B2 (en) Assessing tumor response to therapy
CN105147312A (zh) Pet图像获取方法及系统
CN109961419B (zh) 对pet活度分布图像进行衰减校正的校正信息获取方法
Chan et al. End-expiration respiratory gating for a high-resolution stationary cardiac SPECT system
US20230121358A1 (en) Ct image generation method for attenuation correction of pet images
CN101658426A (zh) 基于三维正电子发射探测器灵敏度特性的呼吸门控技术
Mukherjee et al. Improved frame‐based estimation of head motion in PET brain imaging
CN102764135A (zh) 心肌spect灌注显像的心脏轴向移动校正方法
CN110060234B (zh) 基于spect断层图像计算肾脏gfr值的方法和装置
CN104000618B (zh) 一种用环真光子数门控方法进行呼吸运动门控校正方法
CN104361568A (zh) 基于配准的肺4d-ct图像呼气过程中间相位图像的重建方法
CN110458779A (zh) 对呼吸或心脏的pet图像进行衰减校正的校正信息获取方法
Khurshid et al. Automated cardiac motion compensation in PET/CT for accurate reconstruction of PET myocardial perfusion images
Pourmoghaddas et al. Respiratory phase alignment improves blood‐flow quantification in Rb82 PET myocardial perfusion imaging
CN110063739B (zh) Spect和pet的图像校正方法、装置以及电子设备
Zhang et al. Comparison of different attenuation correction methods for dual gating myocardial perfusion SPECT/CT
CN104361567A (zh) 基于配准的肺4d-ct图像吸气过程中间相位图像的重建方法
Dennerlein et al. Geometric jitter compensation in cone-beam CT through registration of directly and indirectly filtered projections
CN110428384B (zh) 对呼吸或心脏的pet图像进行衰减校正的校正信息获取方法
CN113436708A (zh) 一种基于深度学习算法的延迟ct图像生成方法
JP2003190147A (ja) Ｘ線ｃｔ装置及び画像処理装置
DE112020006198T5 (de) Vorrichtung, System, Verfahren und Computerprogramm zum Bereitstellen eines Kernbildes einer interessiernden Region eines Patienten
Dey et al. Multi-modal rigid and non-rigid registration for attenuation correction in cardiac SPECT/CT using emission scatter to CT conversion
Barnden et al. Detection and validation of the body edge in low count emission tomography images
Lamare et al. Incorporation of elastic transformations in list-mode based reconstruction for respiratory motion correction in PET

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 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