正电子PET动态心肌线粒体显像的心肌血流定量分析方法及用途技术

本发明专利技术是关于一种心脏核医学领域，利用正电子计算机断层显像仪PET的动态心肌线粒体显像，通过图像处理、图像校正与心肌血流定量计算进行心肌血流定量分析的技术方法，以及该技术方法的用途，具体实施方式包括采用：1)一种正电子PET动态心肌线粒体显像的图像处理技术方法，应用于产生未校正靶心图与未校正左心血池时间活度曲线；2)一种正电子PET动态心肌线粒体显像的图像校正技术方法，应用于校正心肌溢出干扰、右心血池溢出干扰及心肌摄取残留；3)一种正电子PET动态线粒体显像的心肌血流定量计算技术方法，应用于心肌血流定量及提高拟合运算效能；4)一种正电子PET动态心肌线粒体显像的心肌血流定量分析技术方法的用途，通过心肌血流级数图的建立，应用于更直观评估整体心肌及血管分区的静息血流量，负荷血流量及血流储备的高低。本发明专利技术的技术方法能够解决以正电子PET动态心肌线粒体显像进行心肌血流定量分析的技术与应用问题。

【技术实现步骤摘要】
正电子PET动态心肌线粒体显像的心肌血流定量分析方法及用途
本专利技术是关于一种心脏核医学领域，利用正电子计算机断层显像仪(positronemissiontomography，正电子PET)的动态心肌线粒体显像，通过图像处理、图像校正与心肌血流定量计算进行心肌血流定量分析的技术方法，以及该技术方法的用途。
技术介绍
在心脏核医学领域，以核素标记的正电子PET显像药物(特别是使用18F核素)进行正电子PET静态心肌线粒体显像的方法逐渐受到国内外的重视，这类心肌显像方法的特征为心肌滞留时间延长，心肌成像轮廓清晰及显像时间点较为弹性，在相对较短的时间内能够同时完成定性的静息与负荷的正电子PET静态心肌线粒体显像，其显像特征符合实际应用要求。但这类正电子PET线粒体显像方法同时也存在一定的缺陷，例如正电子PET显像药物停滞于心肌的过程为通过心肌血流将药物运送到达心肌细胞附近，再通过心肌细胞膜的穿透进入心肌细胞内与线粒体产生作用，而进入心肌细胞内可能还同时发生了从心肌洗脱的现象，这种随时间实时变化的复杂动态过程，当进行正电子PET动态心肌线粒体显像时，由于成像时所使用的PET或PET/CT设备的图像空间分辨率有限，在正电子PET药物未到达心肌细胞前，早期左心血池及右心血池的放射线滞留能够对心肌部位产生放射线溢出干扰，而晚期滞留于左心血池的放射线清除趋近本底后，心肌摄取的放射线开始反向对左心产生放射线溢出干扰，从而发生一种随着时间改变的动态耦合溢出干扰现象，进而导致以正电子PET动态线粒体显像测量左心血池及心肌的放射性浓度，能达到要求的准确性仍然存在很大的挑战难点，即使目前正电子PET动态显像能提供量化的图像格式(Bq/ml或KBq/ml)应用于测量左心血池及心肌的放射线浓度，但该测量仍然存在过大的误差，加上心肌本身对PET显像药物摄取的动态过程过于复杂，导致目前仍然无法通过正电子PET动态心肌线粒体显像，有效实现心肌血流定量测定，从而阻碍了正电子PET动态线粒体心肌显像方法的实际应用发展。鉴于此，我们发展一种解决以正电子PET动态心肌线粒体显像进行心肌血流定量分析的技术方法及其应用。
技术实现思路
本专利技术的目的是提供一种能够通过正电子PET动态心肌线粒体显像进行心肌血流定量分析的技术方法及该技术方法的用途。本专利技术所采用的技术方法之一为一种正电子PET动态心肌线粒体显像的图像处理技术方法，应用于产生未校正靶心图与未校正左心血池时间活度曲线(TimeActivityCurve，TAC)，该技术方法包括采用：(1)将正电子PET动态心肌线粒体图像在一段时间内进行叠加，以产生正电子PET静态心肌线粒体图像；(2)将正电子PET静态心肌线粒体图像以3D平面切割，以区分心肌部位及非心肌部位；(3)将PET静态心肌线粒体图像由直角坐标系转换至球体坐标系，以产生正电子PET静态心肌线粒体图像在球体坐标系的心肌部位图像；(4)以射线采样方式(raytracing)，从经球体坐标系转换的正电子PET静态心肌线粒体图像中取得径向的心肌最大放射线浓度点，以产生3D心肌中线；(5)以球体坐标系、直角坐标系及2D靶心图座标系中的心尖部位为中心参考点，对球体座标径向的心肌最大放射线浓度点进行连续空间向量转换，以取得径向的心肌最大放射线浓度点在2D靶心图坐标系的几何关联性，并用于产生靶心图转换几何矩阵；(6)通过该靶心图几何转换矩阵应用于正电子PET动态心肌线粒体图像的每一帧图像，以产生心肌部位的未校正动态靶心图；(7)通过一采样长方盒自动置于PET静态心肌线粒体图像的左心房与左心室之间，并将该位置应用至正电子PET动态心肌线粒体图像的每一帧图像，以获得未校正左心血池时间活度曲线；本图像处理技术方法同时适用于处理静息及负荷正电子PET动态心肌线粒体图像。本专利技术所采用的技术方法之二为一种正电子PET动态心肌线粒体显像的图像校正技术方法，应用于校正心肌溢出干扰、右心血池溢出干扰及心肌摄取残留，该技术方法包括采用：(1)将左心血池测量的放射线浓度视为其真实放射线浓度及心肌放射线溢出干扰的组合，计算在左心血池受随时间改变的心肌放射线溢出干扰程度，并从未校正左心血池时间活度曲线中减去该随时间改变的心肌放射线溢出干扰程度，从而获得心肌溢出校正左心血池时间活度曲线；(2)以未校正左心血池时间活度曲线中最大放射线浓度作为正电子PET药物进入左心的时间分界点，通过感兴趣区测量未校正动态靶心图中的左心室间隔壁放射线浓度及左心室侧壁放射线浓度，以两者比值校正右心室血池放射线对于左心室间隔的溢出干扰，而获得右心室溢出校正动态靶心图；(3)针对先静息后负荷的正电子PET动态心肌线粒体显像，通过两显像时间起点时间的时间差、正电子PET药物的核素物理半衰期及指数模型，推算产生心肌残留程度靶心图，并将负荷右心室溢出校正动态靶心图扣除该心肌残留程度靶心图，从而产生负荷完整校正动态靶心图；(4)针对先负荷后静息的正电子PET动态心肌线粒体显像，通过两者显像时间起点时间的时间差、正电子PET药物的核素物理半衰期及指数模型，推算产生心肌残留程度靶心图，并将静息右心室溢出校正动态靶心图扣除该心肌残留程度靶心图，从而产生静息完整校正动态靶心图；(5)以一对静息右心室溢出校正动态靶心图及负荷完整校正动态靶心图，或一对负荷右心室溢出校正动态靶心图及静息完整校正动态靶心图，用于产生一组静息完整校正心肌时间活度曲线及一组负荷完整校正心肌时间活度曲线，并进一步产生一条静息完整校正整体心肌时间活度曲线及一条负荷完整校正整体心肌时间活度曲线。本专利技术所采用的技术方法之三为一种正电子PET动态线粒体显像的心肌血流定量计算技术方法，应用于心肌血流定量及提高拟合运算效能，包括采用：(1)双组织三腔室模型及一组合系数C，并将心肌各时间点所测量的心肌放射线浓度视为真实心肌摄取及左心血池对心肌溢出干扰的组合，以产生一项具备(K1，k2，k3，C)动力学参数的动力学组成方程，其中C同时代表左心血池对心肌干扰成度的系数，(1-C)代表真实心肌摄取对心肌放射线浓度测量值的贡献系数，以K1代表心肌血流数值；(2)先采用一组固定初始值，将一组静息完整校正整体心肌时间活度曲线及静息心肌溢出校正左心血池时间活度曲线，与动力学组成方程进行拟合运算，以产生静息整体心肌的动力学参数作为变动初始值，再将一组静息完整校正心肌时间活度曲线及静息心肌溢出校正左心血池时间活度曲线，与动力学组成方程进行拟合运算，从而计算出静息心肌部位各点的动力学参数，再通过心率及血压乘机校正静息心肌血流值，并以靶心图显示，从而产生静息血流靶心图；(3)先采用一组固定初始值，将一组负荷完整校正整体心肌时间活度曲线及静息心肌溢出校正左心血池时间活度曲线，与动力学组成方程进行拟合运算，以产生负荷整体心肌的动力学参数作为变动初始值，再将一组负荷完整校正心肌时间活度曲线及负荷心肌溢出校正左心血池时间活度曲线，与动力学组成方程进行拟合运算，从而计算出负荷心肌部位各点的动力学参数，并以靶心图显示，从而产生负荷血流靶心图，将负荷血流靶心图与静息血流靶心图本文档来自技高网...

【技术保护点】
1.一种正电子PET动态心肌线粒体显像的图像处理技术方法，应用于产生未校正靶心图与未校正左心血池时间活度曲线，该技术方法采用以下步骤：/n(1)图像叠加步骤：将正电子PET动态心肌线粒体图像在一段时间内进行叠加以产生正电子PET静态心肌线粒体图像，叠加时将图像中每个像素的空间及时间坐标先进行标记，并根据叠加的时间范围，将各像素的空间坐标对各时间点进行叠加，以产生正电子PET静态心肌线粒体图像；/n(2)图像分割步骤：将正电子PET静态心肌线粒体图像，在代表影像空间的直角坐标系，以3D平面方程切割区分心肌部位及非心肌部位；/n(3)图像坐标系转换步骤：心肌部位图像以心肌的质量中心为原点，将正电子PET静态心肌线粒体图像由直角坐标系转换至球体坐标系，在坐标转换同时以各像素为中心，通过临近像素的空间座标与像素值，对各中心像素的像素值进行内差计算，以保证球体坐标系中各像素值的连续性；/n(4)心肌中线产生步骤：采用经球体坐标系转换的PET静态心肌线粒体图像，以平行于半径的径向射线，对该径向射线所通过每一个角度坐标对应的像素值进行对比，以取得径向的最大放射线浓度点，并将所有径向的最大放射线浓度点连成一个3D心肌中线；/n(5)靶心图几何转换矩阵产生步骤：在球体座标中，将所有径向的最大放射线浓度点离原点的距离，先排成一个2D极方图，并根据心尖部位在球体座标的位置、直角坐标系的位置及2D靶心图的位置，以心尖部位为中心参考点，对所有径向的最大放射线浓度点进行连续空间向量转换，以计算所有径向的最大放射线浓度点在靶心图中与心尖部位的距离及方位角，从而取得心肌所有径向的最大放射线浓度点在2D靶心图坐标系的几何关联性，并应用该几何关联性于产生靶心图几何转换矩阵；/n(6)动态靶心图产生步骤：通过步骤(3)的图像坐标系转换步骤，将正电子PET动态心肌线粒体图像进行每一帧图像的坐标系转换，再藉由步骤(5)所产生的靶心图几何转换矩阵，将正电子PET动态心肌线粒体图像的每一帧图像转换为靶心图，从而产生了未校正动态靶心图；/n(7)左心血池时间活度曲线产生步骤：在直角坐标中，藉由3D平面及心尖部位的空间信息，以一采样长方盒自动置于正电子PET静态心肌线粒体图像中的左心房与左心室之间，并将该采样长方盒应用至全部正电子PET动态心肌线粒体图像的每一帧图像，以获得未校正左心血池时间活度曲线。/n...

【技术特征摘要】
1.一种正电子PET动态心肌线粒体显像的图像处理技术方法，应用于产生未校正靶心图与未校正左心血池时间活度曲线，该技术方法采用以下步骤：
(1)图像叠加步骤：将正电子PET动态心肌线粒体图像在一段时间内进行叠加以产生正电子PET静态心肌线粒体图像，叠加时将图像中每个像素的空间及时间坐标先进行标记，并根据叠加的时间范围，将各像素的空间坐标对各时间点进行叠加，以产生正电子PET静态心肌线粒体图像；
(2)图像分割步骤：将正电子PET静态心肌线粒体图像，在代表影像空间的直角坐标系，以3D平面方程切割区分心肌部位及非心肌部位；
(3)图像坐标系转换步骤：心肌部位图像以心肌的质量中心为原点，将正电子PET静态心肌线粒体图像由直角坐标系转换至球体坐标系，在坐标转换同时以各像素为中心，通过临近像素的空间座标与像素值，对各中心像素的像素值进行内差计算，以保证球体坐标系中各像素值的连续性；
(4)心肌中线产生步骤：采用经球体坐标系转换的PET静态心肌线粒体图像，以平行于半径的径向射线，对该径向射线所通过每一个角度坐标对应的像素值进行对比，以取得径向的最大放射线浓度点，并将所有径向的最大放射线浓度点连成一个3D心肌中线；
(5)靶心图几何转换矩阵产生步骤：在球体座标中，将所有径向的最大放射线浓度点离原点的距离，先排成一个2D极方图，并根据心尖部位在球体座标的位置、直角坐标系的位置及2D靶心图的位置，以心尖部位为中心参考点，对所有径向的最大放射线浓度点进行连续空间向量转换，以计算所有径向的最大放射线浓度点在靶心图中与心尖部位的距离及方位角，从而取得心肌所有径向的最大放射线浓度点在2D靶心图坐标系的几何关联性，并应用该几何关联性于产生靶心图几何转换矩阵；
(6)动态靶心图产生步骤：通过步骤(3)的图像坐标系转换步骤，将正电子PET动态心肌线粒体图像进行每一帧图像的坐标系转换，再藉由步骤(5)所产生的靶心图几何转换矩阵，将正电子PET动态心肌线粒体图像的每一帧图像转换为靶心图，从而产生了未校正动态靶心图；
(7)左心血池时间活度曲线产生步骤：在直角坐标中，藉由3D平面及心尖部位的空间信息，以一采样长方盒自动置于正电子PET静态心肌线粒体图像中的左心房与左心室之间，并将该采样长方盒应用至全部正电子PET动态心肌线粒体图像的每一帧图像，以获得未校正左心血池时间活度曲线。


2.如权利要求项1中该步骤(1)至步骤(8)可同时适用于静息正电子PET动态心肌线粒体显像及负荷正电子PET动态心肌线粒体显像的图像处理，用于产生静息未校正动态靶心图、静息未校正左心血池时间活度曲线、负荷未校正动态靶心图、负荷未校正左心血池时间活度曲线。


3.一种正电子PET动态心肌线粒体显像的图像校正技术方法，应用于校正心肌溢出干扰、右心血池溢出干扰及心肌摄取残留，该技术方法采用以下步骤：
(1)心肌溢出校正步骤：在左心房及左心室的腔中，以一采样长方盒测量放射线浓度，并将视为(真实左心血池放射线浓度+心肌放射线溢出干扰)的组合，再根据该采样长方盒像素的空间坐标及心肌部位中各像素的空间坐标，以距离平方反比公式计算心肌部位像素与长方盒像素的几何效应矩阵，同时将心肌像素的放射线浓度视为独立点源，针对正电子PET动态心肌线粒体图像中每一帧图像，计算各自的心肌部位像素对长方盒内单一像素的放射线溢出干扰程度，并对各像素进行加总与平均，以取得随时间改变的心肌溢出干扰程度，并从未校正左心血池时间活度曲线中减去该随时间改变的心肌溢出干扰程度，从而获得心肌溢出校正左心血池时间活度曲线；
(2)右心血池溢出校正步骤：以未校正左心血池时间活度曲线中最大放射线浓度作为正电子PET药物进入左心的时间分界点，以2个感兴趣区分别测量该时间分界点之前未校正动态靶心图的每一帧图像中的左心室间隔壁的放射线浓度及左心室侧壁的放射线浓度，并计算各时间点两者间的比值，通过该比值校正右心室血池放射线对于左心室间隔的溢出干扰，而获得右心室溢出校正动态靶心图；
(3)静息心肌摄取残留校正：针对静息正电子PET动态心肌线粒体显像在前，负荷正电子PET动态心肌线粒体显像在后的情况，先以经步骤(2)校正的静息右心室溢出校正动态靶心图中接近或等于最晚期时间点作为静息心肌摄取靶心图，通过静息正电子PET动态心肌线粒体显像的起点时间及负荷正电子PET动态心肌线粒体显像的起点时间的时间差，以正电子PET药物的核素物理半衰期及指数模型，推算静息心肌摄取在负荷正电子PET动态心肌线粒体显像时间起点的心肌残留系数，并将系数与静息心肌摄取靶心图相乘，从而成产生心肌残留程度靶心图，将负荷右心室溢出校正动态靶心图扣除该心肌残留程度靶心图，从而产生负荷完整校正动态靶心图；
(4)负荷心肌摄取残留校正：针对是负荷正电子PET动态心肌线粒体显像在前，静息正电子PET动态心肌线粒体显像在后的情况，先以经步骤(2)校正的负荷右心室溢出校正动态靶心图中接近最晚期时间点作为负荷心肌摄取靶心图，通过负荷正电子PET动态心肌线粒体显像的起点时间及静息PET动态心肌线粒体显像的起点时间的时间差，以正电子PET药物的核素物理半衰期及指数模型，推算负荷心肌摄取在静息正电子PET动态心肌线粒体显像时间起点的心肌残留程度系数，并将系数与负荷心肌摄取靶心图相乘，从产生心肌残留程度靶心图，...

【专利技术属性】
技术研发人员：汪蕾，方纬，许百灵，
申请(专利权)人：方纬，许百灵，刘丽，
类型：发明
国别省市：北京;11