[发明专利]一种获取造影剂迟豫时间的方法

说明书

技术领域

本发明涉及一种磁共振成像领域，特别涉及一种自动化提取磁共振造影剂成像中信号并用以拟合计算造影剂迟豫时间的方法。

背景技术

磁共振中图像对比度反映两组织间信号强度差，对比度增强是改变这信号强度的过程。造影剂，又称对比剂（Contrast agents）是一种药剂，能改变组织的物理特性，从而影响物理参数，进而增强图像对比度。磁共振造影剂通过改变组织的局部磁环境而改变图像对比度。 在磁共振成像中，质子所产生的磁共振信号及其自旋-晶格迟豫时间（纵向迟豫时间）T1和自旋-自旋迟豫时间（横向迟豫时间）T2决定着不同组织在磁共振图像上的对比，磁共振造影剂与质子相互作用来影响T1和T2迟豫时间。

磁共振造影剂根据其磁特性分为顺磁性、超顺磁性和铁磁性三类。 顺磁性对比剂由顺磁性金属元素组成，如Gd、Mn。对比剂浓度低时，主要使T2缩短并使信号增强，浓度高时，则组织T1缩短超过T2效应，使MR信号降低，常用其T1效应作为T1加权像中的阳性对比剂；铁磁性及超顺磁性对比剂由氧化铁组成，为不同大小微晶金属粒子，二者均影响局部磁场均匀性且产生磁化率效应，使质子失相位加速，T2迟豫时间缩短。

因此对于造影剂合成制备以及应用中，T1和T2是造影剂的重要属性，也是衡量造影剂作用与影响的重要指标，利用磁共振技术对造影剂迟豫时间的准确测量成为关键。磁共振成像相比于磁共振波谱技术，对试剂样品的测量有时间快，效率高等优势，目前已成为造影剂迟豫时间测量的重要手段。

磁共振成像技术中对样品试剂的T1迟豫时间通常可采用饱和恢复法或反转恢复法进行测量。饱和恢复法利用SE序列，通过改变成像参数中的恢复时间（TR）值进行成像得到不同TR时间的图像；反转恢复法通过施加一个180度反转预脉冲后紧接一个SE序列构成，把180度翻转预脉冲中点到90度脉冲中点的时间间隔定义为反转时间（TI），通过改变成像参数TI值进行成像得到不同TI时间的图像；对样品的T2迟豫时间的测量可使用多回波SE序列，通过设定单调递增的回波时间（TE）对每个回波信号进行采集成像，从而得到不同TE时间的图像。

对上述三种成像方法得到的不同成像参数时间的图像中的样品选取ROI，抽取图像信号，再通过做非线性拟合方法可以测得T1和T2值。现有的一些图像处理软件如MRIcro是人工手动选取样品的ROI，得到各个时间点ROI内像素信号的平均值，然后再用Origin软件做非线性拟合得到样品T1或T2值；还有一些软件如ImageJ的测量迟豫时间的方法是对图像中的样品信号所有像素做非线性拟合，最终得到一个T1和T2的map图，最后在map图上手动画ROI，以ROI内像素平均值作为样品T1或T2值。

通过上述图像后处理方法可以计算出造影剂的T1和T2时间，但是上述方法的第一个缺点是：不管是在之前原图像选取ROI（Region of Interest）还是之后在T1或T2的map上选取ROI都是人工手动过程，对于多个样品成像（如一次成像试管盒中装有几十支甚至上百支试管样品），人工手动过程非常耗时而且对于ROI中心点定位准确性以及形状选取的一致性很低；第二个缺点为上述第一种方法提取ROI信号后还要利用Origin软件做非线性拟合，其过程分两步进行，需将第一步获取的各个数据人工输入到Origin中，其过程非常耗时。而上述第二种方法虽然直接做非线性拟合，但做T1或T2的map由于是逐个对像素进行非线性拟合计算，一方面效率很低，另一方面磁共振噪声对单个像素信号影响很大，对于先做拟合后取ROI计算平均值的方法准确性不高。

发明内容

本发明是针对现在图像处理方法中获取造影剂迟豫时间耗时、效率低、并且受主观影响的问题，提出了一种获取造影剂迟豫时间的方法，自动化识别图像中造影剂样品，精确定位样品中心点计算样品试管半径大小，并以试管中心点根据试管半径大小自动选取ROI，计算ROI内像素信号平均值，实验标准偏差。最后对自动化提取的信号做非线性拟合，准确快速高效地计算出造影剂的T1和T2值。

本发明的技术方案为：一种获取造影剂迟豫时间的方法，具体包括如下步骤：

1）读取磁共振样品图像数据：对于图像以样品信号与背景信号对比度最好的图像为基准，其中，对于测量自旋-晶格迟豫时间T1的图像数据，以恢复时间TR或反转时间TI最大的图像为基准；对于测量自旋-自旋迟豫时间T2的图像数据，以单调递增的回波时间TE最小的图像为基准；