一种获取造影剂迟豫时间的方法技术

本发明专利技术涉及一种获取造影剂迟豫时间的方法。自动化识别图像中造影剂样品，精确定位样品中心点计算样品试管半径大小，并以试管中心点根据试管半径大小自动选取ROI，计算ROI内像素信号平均值，实验标准偏差。最后对自动化提取的信号做非线性拟合，准确快速高效地计算出造影剂的自旋-晶格弛豫时间T1和自旋-自旋弛豫时间T2值。在图像中识别多个样品选取样品ROI过程为计算机自动化过程，不需要人工过程，在定位精确性、ROI形状一致性、结果准确性远高于传统方法；该方法从ROI选取、ROI信号提取到非线性拟合整个流程为全自动化过程，在时间效率上远远高于传统方法。

【技术实现步骤摘要】
一种获取造影剂迟豫时间的方法
本专利技术涉及一种磁共振成像领域，特别涉及一种自动化提取磁共振造影剂成像中信号并用以拟合计算造影剂迟豫时间的方法。
技术介绍
磁共振中图像对比度反映两组织间信号强度差，对比度增强是改变这信号强度的过程。造影剂，又称对比剂（Contrastagents）是一种药剂，能改变组织的物理特性，从而影响物理参数，进而增强图像对比度。磁共振造影剂通过改变组织的局部磁环境而改变图像对比度。在磁共振成像中，质子所产生的磁共振信号及其自旋-晶格迟豫时间（纵向迟豫时间）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（RegionofInterest）还是之后在T1或T2的map上选取ROI都是人工手动过程，对于多个样品成像（如一次成像试管盒中装有几十支甚至上百支试管样品），人工手动过程非常耗时而且对于ROI中心点定位准确性以及形状选取的一致性很低；第二个缺点为上述第一种方法提取ROI信号后还要利用Origin软件做非线性拟合，其过程分两步进行，需将第一步获取的各个数据人工输入到Origin中，其过程非常耗时。而上述第二种方法虽然直接做非线性拟合，但做T1或T2的map由于是逐个对像素进行非线性拟合计算，一方面效率很低，另一方面磁共振噪声对单个像素信号影响很大，对于先做拟合后取ROI计算平均值的方法准确性不高。
技术实现思路
本专利技术是针对现在图像处理方法中获取造影剂迟豫时间耗时、效率低、并且受主观影响的问题，提出了一种获取造影剂迟豫时间的方法，自动化识别图像中造影剂样品，精确定位样品中心点计算样品试管半径大小，并以试管中心点根据试管半径大小自动选取ROI，计算ROI内像素信号平均值，实验标准偏差。最后对自动化提取的信号做非线性拟合，准确快速高效地计算出造影剂的T1和T2值。本专利技术的技术方案为：一种获取造影剂迟豫时间的方法，具体包括如下步骤：1）读取磁共振样品图像数据：对于图像以样品信号与背景信号对比度最好的图像为基准，其中，对于测量自旋-晶格迟豫时间T1的图像数据，以恢复时间TR或反转时间TI最大的图像为基准；对于测量自旋-自旋迟豫时间T2的图像数据，以单调递增的回波时间TE最小的图像为基准；2）自动化识别样品并定位中心选取ROI提取信号：对步骤1）选取的图像做简单二值化处理，将样品与背景区分；然后利用圆检测霍夫变换的方法，将图像数据变换为圆心累计点的参数图像数据，并将参数图像数据中的峰值点所对应的图像数据点确定为各个造影剂试管的中心；再结合二值化图像，准确计算出各个试管的半径；对每幅图像以试管中心点为ROI中心，相对于试管半径大小以固定比例自动化选取ROI；计算ROI内像素信号的平均值以及实验标准偏差；最后对自动识别后每个样品管进行数字编号；3）非线性拟合计算造影剂T1和T2值：利用步骤2）中提取得到的每幅图像各个造影剂样品的ROI的信号值，以及饱和恢复法对应的每幅图像的成像参数值恢复时间TR、反转恢复法对应的每幅图像的成像参数值反转时间TI做非线性拟合可得到T1值；或多回波序列法对应的每幅图像的成像参数值TE做非线性曲线拟合可得到T2值，对于饱和恢复法测T1的拟合公式为：M=M0*[1-exp(-TR/T1)]；对于反转恢复法测T1的拟合公式为：M=M0*[1-2exp(-TI/T1)]；对于多回波序列测T2的拟合公式为：M=M0*exp(-TE/T2)；公式中M表示提取得到的ROI内的造影剂图像信号值，M0表示M演化过程中的最大值，将不同TR时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0、T1进行拟合计算；对于反转恢复法来说，将不同TI时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0、T1进行拟合计算；对于多回波序列法来说，将不同TE时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0、T2进行拟合计算；非线性拟合的数学方法采用Levenberg-Marquardt迭代方法，最终测出各造影剂的T1或T2值。本专利技术的有益效果在于：本专利技术获取造影剂迟豫时间的方法，在图像中识别多个样品选取样品ROI过程为计算机自动化过程，不需要人工过程，在定位精确性、ROI形状一致性、结果准确性远高于传统方法；该方法从ROI选取、ROI信号提取到非线性拟合整个流程为全自动化过程，在时间效率上远远高于传统方法。附图说明图1为本专利技术实施例中反转恢复SE序列TI参数为5000ms时的造影剂样品磁共振成像图；图2为本专利技术实施例中对图1进行圆检测霍夫变换后的圆心累计点参数域图像图；图3为本专利技术实施例中对各个自动检测出的造影剂样品的编号图。具体实施方式本专利技术的具体方案如下：1、自动化识别样品并定位中心选取ROI提取信号读取磁共振样品图像数据，对于图像以样品信号与背景信号对比度最好的图像为基准。其中，对于测量自旋-晶格迟豫时间T1的图像数据，以恢复时间TR或反转时间TI最大的图像为基准；对于测量自旋-自旋迟豫时间T2的图像数据，以单调递增的回波时间TE最小的图像为本文档来自技高网...

【技术保护点】
一种获取造影剂迟豫时间的方法，其特征在于，具体包括如下步骤：1）读取磁共振样品图像数据：对于图像以样品信号与背景信号对比度最好的图像为基准，其中，对于测量自旋‑晶格弛豫时间T1的图像数据，以恢复时间TR或反转时间TI最大的图像为基准；对于测量自旋‑自旋弛豫时间T2的图像数据，以单调递增的回波时间TE最小的图像为基准；2）自动化识别样品并定位中心选取ROI提取信号：对步骤1）选取的图像做简单二值化处理，将样品与背景区分；然后利用圆检测霍夫变换的方法，将图像数据变换为圆心累计点的参数图像数据，并将参数图像数据中的峰值点所对应的图像数据点确定为各个造影剂试管的中心；再结合二值化图像，准确计算出各个试管的半径；对每幅图像以试管中心点为ROI中心，相对于试管半径大小以固定比例自动化选取ROI；计算ROI内像素信号的平均值以及实验标准偏差；最后对自动识别后每个样品管进行数字编号；3）非线性拟合计算造影剂T1和T2值：利用步骤2）中提取得到的每幅图像各个造影剂样品的ROI的信号值，以及饱和恢复法对应的每幅图像的成像参数值TR、反转恢复法对应的每幅图像的成像参数值TI做非线性拟合可得到T1值；或多回波序列法对应的每幅图像的成像参数值TE做非线性曲线拟合可得到T2值，对于饱和恢复法测T1的拟合公式为：M = M0* [1‑exp(‑TR/T1)]；        对于反转恢复法测T1的拟合公式为：M = M0* [1‑2exp(‑TI/T1)] ；     对于多回波序列测T2的拟合公式为：M = M0*exp(‑TE/T2) ；      公式中M表示提取得到的ROI内的造影剂图像信号值，M0表示M演化过程中的最大值，对于饱和恢复法来说，将不同TR时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0 、T1进行拟合计算；对于反转恢复法来说，将不同TI时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0 、T1进行拟合计算；对于多回波序列法来说，将不同TE时间点以及提取得到的对应的M图像信号值分别代入对应拟合公式中，对M0 、T2进行拟合计算；非线性拟合的数学方法采用Levenberg‑Marquardt迭代方法，最终测出各造影剂的T1或T2值。...

【技术特征摘要】
1.一种获取造影剂迟豫时间的方法，其特征在于，具体包括如下步骤：1）读取磁共振样品图像数据：对于图像以样品信号与背景信号对比度最好的图像为基准，其中，对于测量自旋-晶格迟豫时间T1的图像数据，以恢复时间TR或反转时间TI最大的图像为基准；对于测量自旋-自旋迟豫时间T2的图像数据，以单调递增的回波时间TE最小的图像为基准；2）自动化识别样品并定位中心选取ROI提取信号：对步骤1）选取的图像做简单二值化处理，将样品与背景区分；然后利用圆检测霍夫变换的方法，将图像数据变换为圆心累计点的参数图像数据，并将参数图像数据中的峰值点所对应的图像数据点确定为各个造影剂试管的中心；再结合二值化图像，准确计算出各个试管的半径；对每幅图像以试管中心点为ROI中心，相对于试管半径大小以固定比例自动化选取ROI；计算ROI内像素信号的平均值以及实验标准偏差；最后对自动识别后每个样品管进行数字编号；3）非线性拟合计算造影剂T1和T2值：利用步骤2）中提取得到的每幅图像各个造影剂样品的ROI的信号值，以及饱和恢复法对应的每幅图像的成像参数值...

【专利技术属性】
技术研发人员：王丽嘉，裴孟超，李建奇，聂生东，王远军，
申请(专利权)人：上海理工大学，
类型：发明
国别省市：上海;31