本发明属于脂肪温度测量领域,尤其涉及一种脂肪温度的测量方法和装置。
背景技术:
高强度聚焦超声(highintensityfocusedultrasound,hifu)可利用超声波的穿透性和能量沉积性将体外低能量超声波汇聚于特定深度的靶组织,使焦域内温度迅速升高至60℃,引发蛋白质变性和组织凝固性坏死,而不损伤周围正常组织,hifu治疗已逐渐发展为一种重要的肿瘤治疗手段。hifu治疗过程中实时监控靶组织的温度是治疗成功的关键所在。
磁共振成像由于其多方位立体成像、特有的温度监控、多种组织对比等特点,已成为hifu治疗监控的重要手段。而脂肪组织的温度监控是hifu发展过程中必需要解决的问题,一方面,随着hifu治疗技术的不断发展,其作用部位也在不断扩展,如对富含脂肪的胸部肿瘤进行治疗;另一方面,非目标区域的组织灼伤仍是hifu治疗过程中的首要安全问题,其中以富含脂肪的皮下脂肪层灼伤最为常见。
在温度监控方面,基于质子共振频率转移(protonresonancefrequencyshift,prfs)原理的磁共振温度成像已广泛应用于多种组织治疗过程中的温度监控,该方法主要基于水分子质子共振频率与温度变化呈现较为理想的线性关系。但脂肪内质子的共振频率与水分子不同,且不随温度的改变而改变,使得prfs温度检测在应用于脂肪、及富含脂肪的组织时受限。除质子共振频率外,组织的固有特性t2也与组织温度相关,但目前常用的t2检测方法,利用快速自选回波序列(tse)检测t2从而实现脂肪温度监控,但tse序列检测得到的t2值受激励回波及组织t1影响,并不是脂肪组织的实际t2值。该tse序列连续使用多个180°射频脉冲,在高场应用时常常出现sar值超标问题,且tse序列一次只能扫描一个平面,扫描速度慢,无法实现对温度的实时监控。
技术实现要素:
本发明的目的在于提供一种脂肪温度的测量方法,以解决现有技术通过tse序列检测t2值不准确,并且经常出现sar超标,以及扫描速度慢,无法实现对温度的实时监测的问题。
第一方面,本发明实施例提供了一种脂肪温度的测量方法,所述方法包括:
基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。
结合第一方面,在第一方面的第一种可能实现方式中,所述基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2步骤包括:
基于间隔扫描三回波稳态itess序列扫描一次,得到三个不同阶次的稳态信号;
将三个信号值对应的信号值与不同的自旋-自旋弛豫时间t2值对应的解析解进行比较,均方差最小时得到的t2值即为脂肪组织的自旋-自旋弛豫时间t2值。
结合第一方面的第一种可能实现方式,在第一方面的第二种可能实现方式中,所述将三个信号值与不同的自旋-自旋弛豫时间t2值对应的解析解进行比较,均方差最小时得到的t2值即为脂肪组织的自旋-自旋弛豫时间t2值步骤具体为:
根据公式
根据公式t2=argmin{[s1-m0f1]2+[s0-m0f0]2+[s-1-m0f-1]2}计算得到脂肪组织的自旋-自旋弛豫时间t2值,其中:m0为平衡态时的纵向磁化矢量。
结合第一方面,在第一方面的第三种可能实现方式中,在所述基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2步骤之前,所述方法还包括:
将脂肪组织均匀的升温或降温至预定温度,通过间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
根据预设的多个温度与检测的多个脂肪组织的自旋-自旋弛豫时间t2之间的对应关系进行线性拟合,得到脂肪组织的自旋-自旋弛豫时间t2与脂肪温度的线性对应关系的相关参数。
结合第一方面的第三种可能实现方式,在第一方面的第四种可能实现方式中,所述将脂肪组织均匀的升温或降温至预定温度步骤具体为:
通过水浴法将所述脂肪组织升温或者降温至预定温度。
第二方面,本发明实施例提供了一种脂肪温度的测量方法,所述方法包括:
t2检测单元,用于基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
温度查找单元,用于根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。
结合第二方面,在第二方面的第一种可能实现方式中,所述t2检测单元包括:
扫描子单元,用于基于间隔扫描三回波稳态itess序列扫描一次,得到三个不同阶次的稳态信号;
比较子单元,用于将三个信号值与不同的自旋-自旋弛豫时间t2值对应的解析解进行比较,均方差最小时得到的t2值即为脂肪组织的自旋-自旋弛豫时间t2值。
结合第二方面的第一种可能实现方式,在第二方面的第二种可能实现方式中,所述比较子单元具体用于:
根据公式
p=1-e1cosα-e22(e1-cosα),q=e2(1-e1)(1+cosα);e1=exp(-tr/t1),e2=exp(-tr/t2),其中tr指重复时间,α指翻转角,,t1为自旋-晶格弛豫时间;
根据公式t2=argmin{[s1-m0f1]2+[s0-m0f0]2+[s-1-m0f-1]2}计算得到脂肪组织的自旋-自旋弛豫时间t2值,其中:m0为平衡态时的纵向磁化矢量。
结合第二方面,在第二方面的第三种可能实现方式中,所述装置还包括:
温度设定单元,用于将脂肪组织均匀的升温或降温至预定温度,通过间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
拟合单元,用于根据预设的多个温度与检测的多个脂肪组织的自旋-自旋弛豫时间t2之间的对应关系进行线性拟合,得到脂肪组织的自旋-自旋弛豫时间t2与脂肪温度的线性对应关系的相关参数。
结合第二方面的第三种可能实现方式,在第二方面的第四种可能实现方式中,所述温度设定单元具体用于:
通过水浴法将所述脂肪组织升温或者降温至预定温度。
在本发明中,基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2,根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。本方法基于itess序列t2检测的脂肪温度测量,不受单纯脂肪组织约束,可进行快速检测,序列所使用的射频脉冲翻转角小,不易出现sar值超标问题。
附图说明
图1是本发明实施例提供的脂肪温度的测量方法的实现流程图;
图2是本发明实施例提供的脂肪温度的测量装置的结构示意图;
图3是本发明实施例提供的间隔扫描三回波稳态itess序列的时序示意图;
图4是本发明实施例提供的自旋-自旋弛豫时间t2与脂肪温度的对应关系实验结果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例的目的在于提供一种脂肪温度的测量方法,以解决现有技术的磁共振温度监控的方法中,基于质子共振频率转移(protonresonancefrequencyshift,prfs)原理的磁共振温度成像已广泛应用于多种组织治疗过程中的温度监控,该方法主要基于水分子质子共振频率与温度变化呈现较为理想的线性关系。但脂肪内质子的共振频率与水分子不同,且不随温度的改变而改变,使得prfs温度检测在应用于脂肪、及富含脂肪的组织时受限。而现有技术中除质子共振频率外,目前通过检测自旋-自旋弛豫时间t2的方法,检测速度相对较慢,无法实现对温度的实时监测。下面结合附图,对本发明作进一步的说明。
图1示出了本发明第一实施例提供的脂肪温度的测量方法的实现流程,详述如下:
在步骤s101中,基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2。
具体的,所述基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2步骤包括:
基于间隔扫描三回波稳态itess序列扫描一次,得到三个不同阶次的稳态信号。
将三个信号值与不同的自旋-自旋弛豫时间t2值对应的解析解进行比较,均方差最小时得到的t2值即为脂肪组织的自旋-自旋弛豫时间t2值。
比如,基于间隔扫描三回波稳态itess序列扫描一次,得到三个不同阶次的稳态信号分别为s-1,s0,s1,那么三个不同阶次的稳态信号对应的解析解分别为:
根据公式t2=argmin{[s1-m0f1]2+[s0-m0f0]2+[s-1-m0f-1]2}计算得到脂肪组织的自旋-自旋弛豫时间t2值,其中:m0为平衡态时的纵向磁化矢量。
在步骤s102中,根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。
所述对应关系可以为脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应表,也可以为脂肪温度与脂肪组织的自旋-自旋弛豫时间t2对应的线性函数。在根据对应关系表或者线性函数查找脂肪组织的自旋-自旋弛豫时间t2对应的脂肪温度前,还需要建立脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系表,或者拟合脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的线性函数。具体可以为:
将脂肪组织均匀的升温或降温至预定温度,通过间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
根据预设的多个温度与检测的多个脂肪组织的自旋-自旋弛豫时间t2之间的对应关系进行线性拟合,得到脂肪组织的自旋-自旋弛豫时间t2与脂肪温度的线性对应关系的相关参数。
其中,通过间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2的具体过程,与步骤s101中获取脂肪组织的自旋-自旋弛豫时间t2的过程基本相同。
比如,通过多次测量后,得到脂肪组织的自旋-自旋弛豫时间t2与脂肪温度的线性对应关系可以表示为:temp=a*t2+b,其中temp为脂肪温度,t2为脂肪组织的自旋-自旋弛豫时间,a、b分别为线性函数的斜率和截距。
在高强度聚焦超声(highintensityfocusedultrasound,hifu)治疗过程中,可以利用本方法基于间隔扫描三回波稳态itess序列进行扫描,根据求解关系式temp=a*t2+b得到脂肪组织内的温度及其变化情况。
本方法基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2,根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。本方法基于itess序列t2检测的脂肪温度测量,不受单纯脂肪组织约束,可进行快速检测,序列所使用的射频脉冲翻转角小,不易出现sar值超标问题。
图2为本发明实施例提供的脂肪温度的测量装置的结构示意图。如图2所示,所述脂肪温度的测量装置,包括:
t2检测单元,用于基于间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
温度查找单元,用于根据预先标定的脂肪温度与脂肪组织的自旋-自旋弛豫时间t2的对应关系,查找检测的自旋-自旋弛豫时间t2所对应的脂肪温度。
优选的,所述t2检测单元包括:
扫描子单元,用于基于间隔扫描三回波稳态itess序列扫描一次,得到三个不同阶次的稳态信号;
比较子单元,用于将三个信号值对应的信号值与不同的自旋-自旋弛豫时间t2值对应的解析解进行比较,均方差最小时得到的t2值即为脂肪组织的自旋-自旋弛豫时间t2值。
优选的,所述比较子单元具体用于:
根据公式
p=1-e1cosα-e22(e1-cosα),q=e2(1-e1)(1+cosα);其中tr指重复时间,α指翻转角,t1为自旋-晶格弛豫时间
根据公式t2=argmin{[s1-m0f1]2+[s0-m0f0]2+[s-1-m0f-1]2}计算得到脂肪组织的自旋-自旋弛豫时间t2值,其中:其中:m0为平衡态时的纵向磁化矢量。e1=exp(-tr/t1),e2=exp(-tr/t2),p=1-e1cosα-e22(e1-cosα),q=e2(1-e1)(1+cosα);其中tr指重复时间,α指翻转角,t1为自旋-晶格弛豫时间;
根据公式t2=argmin{[s1-m0f1]2+[s0-m0f0]2+[s-1-m0f-1]2}计算得到脂肪组织的自旋-自旋弛豫时间t2值,其中:m0为平衡态时的纵向磁化矢量。
优选的,所述装置还包括:
温度设定单元,用于将脂肪组织均匀的升温或降温至预定温度,通过间隔扫描三回波稳态itess序列检测脂肪组织的自旋-自旋弛豫时间t2;
拟合单元,用于根据预设的多个温度与检测的多个脂肪组织的自旋-自旋弛豫时间t2之间的对应关系进行线性拟合,得到脂肪组织的自旋-自旋弛豫时间t2与脂肪温度的线性对应关系的相关参数。
优选的,所述温度设定单元具体用于:
通过水浴法将所述脂肪组织升温或者降温至预定温度。
图2所述脂肪温度的测量装置,与图1中所述的脂肪温度的测量方法对应,在此不作重复赘述。
为了验证本发明的可行性,在西门子3t磁共振成像系统(siemenstimtrio,erlangen,german)上进行实验,实验对象为脂肪仿体,先将仿体加热至约50°,使其在检测过程中自然冷却降温至大约室温20°,脂肪仿体中插有光纤温度计,对脂肪仿体的温度进行同步检测。磁共振检测序列为itess序列,序列时序图如图3所示,扫描的序列参数为:重复时间tr=10ms,回波时间te=5ms,翻转角fa=20°,层厚=3.0mm,采集层数=24,视野fov=192mm,矩阵大小matrixsize=192*192,带宽bandwidth=606hz/pixel。
图4为检测得到的脂肪仿体的自旋-自旋弛豫时间t2随温度改变情况,脂肪组织内温度变化与t2存在线性相关,拟合曲线斜率为1.274,即温度每升高1℃,检测到的t2值增大1.274ms。
在本发明所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:u盘、移动硬盘、只读存储器(rom,read-onlymemory)、随机存取存储器(ram,randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。