一种利用MRI动态图像进行多采样的灰度波谱分析方法技术

本发明专利技术公开了一种利用MRI动态图像进行多采样的灰度波谱分析方法，包括：在MRI视频上截取一个心脏的多个周期运动图像；利用一种结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法，获取图像周期灰度值变化数据；将图像灰度周期值的波谱数据同对应的患者的心动周期心电图重叠，选取心电图R波峰顶作为等容收缩期开始，标记出之后图像变化率最小区域；利用核回归方法或反向传播神经网络和时频阶谱方法建立心脏等容收缩期信息，从而判定心脏相对静止期。与既往在心脏MRI电影图像回放中人工观察右冠状动脉最小运动时相识别心脏相对静止期相比，本发明专利技术极大地降低了人工偏倚，使心脏相对不应期的判定更加简捷。使心脏相对不应期的判定更加简捷。使心脏相对不应期的判定更加简捷。

【技术实现步骤摘要】
一种利用MRI动态图像进行多采样的灰度波谱分析方法


[0001]本专利技术涉及数字图像分析处理
，尤其涉及一种利用MRI动态图像进行多采样的灰度波谱分析方法。

技术介绍

[0002]现今心血管疾病仍是全球居民死亡的主要原因，其中冠状动脉粥样硬化性疾病(CAD)发生率较为显著。已知X射线冠状动脉造影和冠状动脉CT血管造影(CCTA)是评估和诊断CAD的首选成像方式(“金标准”)。遗憾的是，以上方法存有电离辐射和使用碘化造影剂的风险。冠状动脉磁共振血管成像(MRCA)作为一种安全有效的无创冠状动脉成像技术已逐步应用于临床，无电离辐射、可不注射对比剂成像的优势使其在冠状动脉疾病诊断中具有独特价值，临床需求逐步增强。在儿童、肾功能不全、孕妇、对比剂过敏等特殊患者中的临床应用价值更为突出。
[0003]在心脏磁共振成像过程中，心脏随心动周期不断做有规律的搏动，而冠状动脉紧贴心脏表面运动，常规成像必然出现运动伪影。因此，为了获得真实反映冠状动脉的图像，必须选择心脏运动幅度最小的时相(心脏相对静止期)。既往是通过心脏电影回放，然后人工观察右冠状动脉中段运动最小的时相，从而判断心脏相对静止期，但这种人工观察法存在一定的人工偏倚。

技术实现思路

[0004]本专利技术实施例所要解决的技术问题在于，针对通过心脏电影回放，然后人工观察右冠状动脉中段运动最小的时相，从而判断心脏相对静止期，存在一定的人工偏倚的问题，提出了一种利用MRI动态图像进行多采样的灰度波谱分析方法。
[0005]为了解决上述技术问题，本专利技术实施例提供了一种利用MRI动态图像进行多采样的灰度波谱分析方法，该利用MRI动态图像进行多采样的灰度波谱分析方法包括：
[0006]步骤一：利用患者的MRI视频以及在该MRI视频对应时间内的患者心电图，在MRI视频上截取一个心脏的多个周期运动图像；
[0007]步骤二、利用一种结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法，获取图像周期灰度值变化数据；
[0008]步骤三、将图像灰度周期值的波谱数据同对应的患者的心动周期心电图重叠，选取心电图R波峰顶作为等容收缩期开始，标记出之后图像变化率最小区域；
[0009]步骤四、利用核回归方法或反向传播神经网络和时频阶谱方法建立心脏等容收缩期信息，从多个周期灰度值变化数据中找出等容收缩期中点，从而判定心脏相对静止期。
[0010]其中，步骤一中利用患者的MRI视频，在视频上截取一个心脏的多个周期运动图像的具体方法为：在心脏磁共振电影成像视频中的多个心动周期内，连续采集心脏多个期相的亮血系列图像；选取心脏四腔心层面电影图像进行下一步分析。
[0011]其中，步骤二中一种结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方
法为：于心脏电影图像上手工选取多个采样点，且每个采样点均局限在心脏范围内，并且有肉眼可见的灰度变化；利用结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法对各个随机采样点进行信号采集；获得各个采样点对应的灰度波谱。
[0012]其中，步骤二中采样的方法为：判断图像中存在周期性变化的区域；选取周期性区域像素点周围3*3的像素点进行像素采集，并且计算其平均值。
[0013]其中，获取图像周期灰度值变化数据的具体方法为：
[0014]S101：手动在图像上选取3
‑
10个肉眼能够看到明显的周期变化的区域；
[0015]S102：记录该区域的灰度值变化，并且将灰度值变化和时间t之间的关系形成一个函数；
[0016]S103：利用短时傅立叶变换对图像进行分析；
[0017]S104：将获得的灰度值变化周期分为N(3<＝N<＝10)个时间窗,每个时间窗用傅立叶变换进行处理，其表示方法为：
[0018][0019]上式中,ω[n]是窗函数，x[n]是待变换的信号；
[0020]S105：由上面的函数可以计算获得灰度随时间分段变化的振幅和频率：
[0021][0022]也可以表示为三角函数的形式：
[0023][0024]即复杂周期性变换信号可以表示为一系列正弦或者余弦函数的和。
[0025]其中，步骤四中从多个周期灰度值变化数据中自动判定心脏相对静止期的方法为：利用核回归方法/反向传播神经网络和时频阶谱方法建立心脏等容收缩期信息；在心脏MRI电影视频中自动判定等容特征区间；联合心电图信息，在心脏MRI电影视频中自动识别心脏相对静止期。
[0026]其中，自动判定心脏相对静止期的具体方法为：在公式(3)的基础上，利用每一段图像灰度变化的加权平均值进行估计：
[0027][0028]其中K
h
是一个带宽为h的核，分母是一个总和为1的加权像，由于此时函数可以被看成是一系列正弦和余弦的和的形式，因此上面的公式可以变化为：
[0029][0030]其中由以上可以获得灰度随时间分段变化的振幅和频率；
[0031]一般情况下，等容收缩期大致等于心动周期的1/16，因此，公式5中选择
[0032][0033]式中：p为心率，将一个周期中不同的时间t代入公式(5)，选取最大ω1所在的时间t为等容收缩期中点。
[0034]其中，等容收缩期中点的具体判断方法为：计算每个点的ω值，如果多个采样点ω值存在相同的时间t位置，且其时间不超过0.01秒，那么可以认定该ω可以标定等容收缩期的中点。
[0035]其中，根据等容收缩期起点和中点信息，在MRI图像中标记出整个等容收缩期。
[0036]实施本专利技术实施例，具有如下有益效果：
[0037]应用本利用MRI动态图像进行多采样的灰度波谱分析方法，通过对心脏MRI运动周期图像的灰度波谱分析，可自动识别心脏相对静止期；与既往在心脏MRI电影图像回放中人工观察右冠状动脉最小运动时相识别心脏相对静止期相比，本专利技术极大地降低了人工偏倚，使心脏相对不应期的判定更加简捷。
附图说明
[0038]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0039]图1是本专利技术提供的利用MRI动态图像进行多采样的灰度波谱分析方法的流程示意图。
具体实施方式
[0040]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0041]本专利技术提供了一种利用MRI动态图像进行多采样的灰度波谱分析方法，请参见图1，图1是本专利技术提供的利用MRI动态图像进行多采样的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用MRI动态图像进行多采样的灰度波谱分析方法，包括：步骤一：利用患者的MRI视频以及在该MRI视频对应时间内的患者心电图，在MRI视频上截取一个心脏的多个周期运动图像；步骤二、利用一种结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法，获取图像周期灰度值变化数据；步骤三、将图像灰度周期值的波谱数据同对应的患者的心动周期心电图重叠，选取心电图R波峰顶作为等容收缩期开始，标记出之后图像变化率最小区域；步骤四、利用核回归方法或反向传播神经网络和时频阶谱方法建立心脏等容收缩期信息，从多个周期灰度值变化数据中找出等容收缩期中点，从而判定心脏相对静止期。2.根据权利要求1所述的利用MRI动态图像进行多采样的灰度波谱分析方法，其特征在于，所述步骤一中利用患者的MRI视频，在视频上截取一个心脏的多个周期运动图像的具体方法为：在心脏磁共振电影成像视频中的多个心动周期内，连续采集心脏多个期相的亮血系列图像；选取心脏四腔心层面电影图像进行下一步分析。3.根据权利要求1所述的利用MRI动态图像进行多采样的灰度波谱分析方法，其特征在于，所述步骤二中一种结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法为：于心脏电影图像上手工选取多个采样点，且每个采样点均局限在心脏范围内，并且有肉眼可见的灰度变化；利用结合短时傅里叶变换和速度频率阶次法的时间频率阶次谱方法对各个随机采样点进行信号采集；获得各个采样点对应的灰度波谱。4.根据权利要求3所述的利用MRI动态图像进行多采样的灰度波谱分析方法，其特征在于，所述步骤二中采样的方法为：判断图像中存在周期性变化的区域；选取周期性区域像素点周围3*3的像素点进行像素采集，并且计算其平均值。5.根据权利要求4所述的利用MRI动态图像进行多采样的灰度波谱分析方法，其特征在于，所述获取图像周期灰度值变化数据的具体方法为：S101：手动在图像上选取3
‑
10个肉眼能够看到明显的周期变化的区域；S102：记录该区域的灰度值变化，并且将灰度值变化和时间t之间的...

【专利技术属性】
技术研发人员：文凌仪，彭谨，胡雷，伏川，周忠琴，阿折十干莫，郭应坤，
申请(专利权)人：四川大学华西第二医院，
类型：发明
国别省市：