一种心脏灌注磁共振图像的配准方法与流程

文档序号:12064745阅读:347来源:国知局
一种心脏灌注磁共振图像的配准方法与流程

本发明涉及医学图像领域,尤其是涉及心脏灌注磁共振图像的配准。



背景技术:

近年来,心血管疾病的发病率和死亡率正逐年增加。全世界每年几千万人死于心血管疾病。心肌缺血是引发心血管疾病发病和死亡的最重要原因。心脏灌注磁共振成像被认为是检查缺血性心脏病的首选无创型检查方法。根据心脏灌注磁共振图像的基本原理,放射科医生可以通过直接观察左心室心肌部分的信号强度随时间的变化情况进行心肌缺血的临床诊断。但是,在图像获取过程中,心脏运动、呼吸及病人位置移动在心脏灌注磁共振图像中产生的形变和运动伪影、快速时间采样、及相对较低的图像对比度,限制了直接进行视觉诊断的可靠性和效率。



技术实现要素:

本发明所要解决的技术问题是提供一种全自动化的心脏灌注磁共振图像的配准方法。

本发明为解决上述技术问题而采用的技术方案是:一种心脏灌注磁共振图像的配准方法,其包括以下步骤:

获取包含左心室心肌的若干层切片的磁共振图像INT,其中N表示在同一心 跳周期中切片所在层的序号,T表示不同心跳周期的序号,N、T均为大于或等于1的整数;

选定起始层切片在T个心跳周期内的磁共振图像I1T,分割磁共振图像I1T中的心肌内膜,确定第r个心跳周期所对应的磁共振图像为参考图像I1r以及感兴趣区域ROI1

分别选定第二层至第N层切片在T个心跳周期内的磁共振图像IMT,设定当前层切片在T个心跳周期内的磁共振图像为IMT,以前一层切片在T个心跳周期内的磁共振图像中选择出的参考图像I(M-1)r的心跳周期的序号r,感兴趣区域ROIM-1为参考依据,分割第M层切片磁共振图像IMT中的心肌内膜,获得第M层切片IMT磁共振图像中的参考图像IMr

以每层切片的磁共振图像中的参考图像IMr为基准,完成本层切片在T个心跳周期内的磁共振图像的配准。

优选的,所述起始层切片位于左心室中间层以下,如心底层。

优选的,所述分割起始层切片在不同心跳周期的磁共振图像I1T中的心肌内膜,包括以下步骤:

a)从心跳周期总数的1/4处开始到2/3处结束,以每间隔2个心跳周期所对应的磁共振图像I1T被选定为候选图像,并以候选图像中心确定初始感兴趣区域ROI1’,在初始感兴趣区域ROI1’中做模糊C均值聚类,聚类数目为2;

b)在每一个候选图像对应的初始感兴趣区域ROI1’的聚类二值图中选取圆度最大的区域,此区域被定义为候选图像的心肌内膜所界定的血池区域;

c)根据所有血池区域的质心点位置、平均灰度和圆度,选择出起始层切片的磁共振图像的血池区域;

d)对血池区域求凸包并平滑,其中心点和长轴可用于设置最终的感兴趣区 域ROI1的位置和大小。

优选的,所述分割第二层至第N层切片的在T个心跳周期的磁共振图像IiT中的心肌内膜,包括以下步骤:

a)从获取第M层切片的磁共振图像所需的心跳周期总数的1/4处开始到第M-1层切片磁共振图像中的参考图像对应的心跳周期序数加3处结束,每间隔2个心跳周期所对应的磁共振图像被选定为第M层切片的候选图像,第M层切片的候选图像分别以第M-1层切片所对应的参考图像的心肌内膜所界定的血池的中心和长轴加若干个像素作为第M层切片所对应的候选图像的心肌内膜所界定的血池的初始感兴趣区域ROIM’的中心和边长,在初始感兴趣区域中做模糊C均值聚类,聚类数目为2;

b)在每一个第M层切片的候选图像的聚类二值图中选取与第M-1层切片的血池重合度最大、且长轴小于第M-1层切片的参考图像的血池长轴的1.1倍的区域被定义为相对应的候选图像的血池区域;如果没找到相对应的候选图像的血池,聚类数目加1,继续做模糊C均值聚类,直到出现合适的候选图像的血池区域;

c)根据所有第M层切片的候选图像的血池区域的质心点位置、平均灰度和圆度选择出第M层切片的候选图像的血池区域,且该候选图像为第M层切片的参考图像IMr

d)对第M层切片的候选图像血池区域求凸包并平滑,其中心点和长轴可用于设置最终感兴趣区域的位置和大小,具体为以血池区域质心为中心点,血池长轴加上2*心肌厚度(20mm)为感兴趣区域的边长。

优选的,每层切片的磁共振图像的配准包括以本层中的参考图像IMr为基准,采用刚体配准或非刚体配准。

优选的,所述刚体配准包括:从第M层切片的参考图像IMr向两端各心跳周期的磁共振图像配准。

优选的,将第M层切片第t个心跳周期的磁共振图像IMt进行形变,使得形变图像与参考图像IMr及第t-1个(前一个)或t+1个(后一个)心跳周期磁共振图像配准后的图像IM(t-1)的相似度之和达到最大;

优选的,刚体配准分两步进行形变,针对第t个心跳周期的磁共振图像IMt,第一步为大尺度形变,第二步为小尺度形变。

优选的,包括判断第M层切片的磁共振图像IMt经过刚体配准后得到的磁共振图像IMt’是否需要非刚体配准。

优选的,经过刚体配准后获得的磁共振图像IMt’与参考图像INr的相似度变量s_ref,与其前一相位或后一相位的磁共振图像IM(t-1)、IM(t+1)的相似度变量s_pre,若两相似度变量s_ref、s_pre的一阶导的平均值>0.03,则该经过刚体配准后得到的磁共振图像IMt’还需要进行非刚体配准。

本发明对比现有技术有如下的有益效果:本发明专利可自动确定参考图像、感兴趣区域和左心室心肌相关的心内膜,结合采用刚体和非刚体配准方法实现四维心肌灌注磁共振图像在时序上的全局和局部配准。

【附图说明】

图1为本发明实施例中获取若干层切片在若干个心跳周期内的心脏磁共振图像的示意图;

图2为本发明实施例中心脏灌注磁共振图像分割过程所对应图像;

图3为本发明实施例中心脏灌注磁共振图像的配准方法流程图;

图4为本发明实施例中心脏灌注磁共振图像的配准前、配准后的图像。

【具体实施方式】

下面结合附图和实施例对本发明作进一步的描述。

一种心脏灌注磁共振图像的配准方法,其包括以下步骤:

获取包含左心室心肌的若干层切片的磁共振图像INT,其中N表示在同一心跳周期中切片所在层的序号(数),T表示不同心跳周期的序号(数),N、T均为大于或等于1的整数;对N层切片中的每一相应的层(M层)的切片的磁共振图像IMT进行以下处理:选定起始层(M=1)切片在T个心跳周期内的磁共振图像I1T,分割磁共振图像I1T中的心肌内膜,确定磁共振图像I1T中的第r个心跳周期所对应的磁共振图像为参考图像I1r以及感兴趣区域ROI1;其中,r为小于T的整数;

分别选定第二层至第N层切片在T个心跳周期内的磁共振图像IMT,设定第M层(当前层)切片在T个心跳周期内的磁共振图像为IMT,以第M-1层(其前一层)切片在T个心跳周期内的磁共振图像中选择出的参考图像I(M-1)r的心跳周期的序号r,感兴趣区域ROIM-1为参考依据,分割第M层切片磁共振图像IMT中的心肌内膜,获得第M层切片IMT磁共振图像中的参考图像IMr;其中,M为大于1且小于等于N的整数;

以每层切片的磁共振图像中的参考图像为基准,完成本层切片在T个心跳周期内的磁共振图像的配准。

其中,所述磁共振图像INT可通过以下方式获得,首先,向被扫描对象注射示踪剂(或其他药物),在对比剂到达心脏之前、之中、之后的时间,用磁共振成像设备获得心脏不同部位的切片(slice)的磁共振图像。本实施例中,沿心脏的长轴方向(大致为上下方向),获得6个不同层面的切片图像,即N 的取值为6,但是,可根据具体的需求,对N取不同的值,例如N的值也可以为4、8、10等。另外,需要在不同个连续的心跳周期T或心脏时相(phase)获得每一层面的切片磁共振图形,T的取值范围为40-60,假定每个心跳周期的时间为Δt,则完成磁共振扫描所需的总时间大约为Δt*T。在每个心跳周期内,需要完成一次N个层面的切片的磁共振图像,因此,一共需要有N*T幅磁共振图像。

其中,起始层切片位于左心室中间层以下,特别的所述起始层切片优选为心底层,此处由心肌内膜所界定(包围)的血池(blood pool)最大。

其中,分割起始层(定义M=1时的层为起始层)切片在不同心跳周期的磁共振图像I1T中的心肌内膜,包括以下步骤:

a)从心跳周期总数的1/4(T/4)处开始到2/3(2*T/3)处结束,以每间隔2个心跳周期所对应的磁共振图像I1T被选定为候选图像,并以候选图像中心确定初始感兴趣区域ROI1’(参图2a),在初始感兴趣区域ROI1’中做模糊C均值聚类;优选的,所述初始感兴趣区域ROI0的形状为正方形,所述初始感兴趣区域ROI0边长为111像素,设定聚类数目为2。

b)在每一个候选图像对应的初始感兴趣区域ROI1’的聚类二值图中选取圆度最大的区域,此区域被定义为候选图像的心肌内膜(参图2c)所界定的血池区域;圆度的定义为(周长*周长)/(4*PI*面积);

c)根据所有血池区域的质心点位置、平均灰度和圆度,选择出起始层切片的磁共振图像的血池区域(参图2b);其中,所述质心点位置越远,平均灰度越大且圆度越大的区域为起始层血池区域的概率越大,且以该血池区域所在的图像作为参考图像;

d)对血池区域求凸包并平滑,其中心点和长轴可用于设置最终的感兴趣区 域ROI1的位置和大小。

进一步的,分割第二层至第N层切片在T个心跳周期内的磁共振图像IMT中的心肌内膜,包括以下步骤:

a)从获取第M层切片(当前层)的磁共振图像所需的心跳周期总数的1/4处(T/4)开始到第M-1(上一层)层切片磁共振图像中的参考图像对应的心跳周期总数(序数)加3处结束,每间隔2个心跳周期所对应的磁共振图像被选定为第M层切片的候选图像,第M层切片的候选图像分别以第M-1层切片所对应的参考图像的心肌内膜所界定的血池的中心和长轴加若干个像素作为第M层切片所对应的候选图像的心肌内膜所界定的血池的初始感兴趣区域ROIM’的中心和边长,在初始感兴趣区域中做模糊C均值聚类,目前设定聚类数目为2类;所述像素数为5个至20个,优选为10个。

b)在每一个第M层切片的候选图像的聚类二值图中选取与第M-1层切片的参考图像的血池重合度最大、且长轴小于第M-1层切片的参考图像的血池长轴的1.1倍的区域被定义为相对应的候选图像血池区域;如果没找到相对应的候选图像的血池,聚类数目加1,继续做模糊聚类,直到出现合适的候选图像血池区域;

c)根据所有第M层切片的候选图像的血池区域的质心点位置、平均灰度和圆度,即质心点位置越远,平均灰度越大且圆度越大的区域为血池区域的概率越大,选择出第M层切片的候选图像的血池区域,且该候选图像为第M层切片的参考图像IMr

d)对第M层切片的候选图像血池区域求凸包并平滑,其中心点和长轴可用于设置最终感兴趣区域的位置和大小,具体为以血池区域质心为中心点,血池长轴加上2*心肌厚度为感兴趣区域的边长,可选的,心肌厚度范围为6-20mm。

进一步的,对于每层切片的磁共振图像的配准包括以本层中的参考图像IMr为基准,采用刚体配准和非刚体配准,当刚体配准没有取得较好效果时,进一步采用非刚体配准。其中M为大于或等于1,小于或等于N的整数。

所述刚体配准包括:从第M层(当前层)切片的参考图像IMr向两端各心跳周期的磁共振图像配准;即以第t个心跳周期所对应的参考图像IMr分别向心跳周期序号数递减和心跳周期序号数递增所对应的同一层切片的磁共振图像配准。

具体的,将第M层切片在第t个心跳周期的磁共振图像IMt进行形变,使得形变图像IMt’与参考图像IMr及第t-1个(前一个)心跳周期配准后的磁共振图像IM(t-1)‘的相似度之和达到最大,此时t>r;或者使得形变图像IMt’与参考图像IMr及第t+1个(后一个)心跳周期配准后的磁共振图像IM(t+1)’的相似度之和达到最大,此时t<r。

其中,相似度准则为sobel梯度的幅值和角度信息,采用的优化方法为下降单纯型法。

其中,刚体配准分两步进行形变,针对第t个心跳周期的磁共振图像IMt,第一步为大尺度(具体为13像素)形变(平移+缩放),第二步为小尺度(具体为3像素)形变(平移+缩放+旋转)。

进一步的,包括判断第t个心跳周期的磁共振图像IMt经过刚体配准后得到的磁共振图像IMt’是否需要非刚体配准。

进一步的,经过刚体配准后获得的第t个心跳周期的磁共振图像IMt’与参考图像INr的相似度变量(s_ref),与其第t-1个(前一个)或第t+1个(后一个)心跳周期形变后的磁共振图像IN(t-1)’的相似度变量(s_pre),若两变量(s_ref、s_pre)的一阶导的平均值>0.03。

从第M层(当前层)切片的参考图像IMr向两端同一层切片在各心跳周期的磁共振图像配准。

利用demons非刚体配准方法对当前图像进行形变,与其相应的伪真实图像进行配准。相似度准则为sobel梯度的幅值和角度信息。

若当前图像为参考图像IMr前、后一个心跳周期的图像IM(r-1)、IM(r+1),则当前图像对应的伪真实图像PM(r-1)、PM(r+1)为参考图像INr,否则,当前磁共振图像IMt的伪真实图像PMt为第t-1个(前一个)心跳周期的磁共振图像的形变图像IM(t-1)’(*0.3)与第t-1个(前一个)心跳周期的伪真实图像PM(t-1)(*0.7)的加权平均,此时t>r。或者,当前磁共振图像IMt的伪真实图像为第t+1个(后一个)心跳周期的磁共振图像的形变图像IM(t+1)’(*0.3)与第t+1个(后一个)心跳周期的伪真实图像PM(t+1)(*0.7)的加权平均,此时t<r。

本发明专利提出方法步骤中涉及到的参数根据实际医学图像特点可以进行任意设置;

本发明专利中刚体配准和非刚体配准采用的相似性测度可以进行变化,如应用互信息。

本发明专利的主要优势在于自动确定参考图像、感兴趣区域和左心室心肌相关的心内膜,结合采用刚体和非刚体配准方法实现四维心肌灌注磁共振图像在时序上的全局和局部配准。

虽然本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。

当前第1页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1