基于纵向弛豫时间的低场磁共振无创测温技术的研究.pdf

《基于纵向弛豫时间的低场磁共振无创测温技术的研究.pdf》由会员分享，可在线阅读，更多相关《基于纵向弛豫时间的低场磁共振无创测温技术的研究.pdf（5页珍藏版）》请在三一文库上搜索。

1、第2 7 卷第6 期 2 0 0 8 年1 2 月 北京生物医学工程 B e u i n gB i o m e d i c a lE n g i n e e r i n g V 0 1 2 7N o 6 D e e e m b e r2 0 0 8 基于纵向弛豫时问的低场磁共振无创 测温技术的研究 王莹h 2 宋涛1王慧贤1 胡丽丽 摘要磁共振成像( M R I ) 的一些成像参数对温度变化敏感，通过测量成像参数平衡磁化强度肘。和纵向弛 豫时间r 1 ，可以实现无创的温度测量该方法尤其适合在低场M R I 中使用。本文在介绍磁共振无创测温技术基本 原理的基础上，着重分析了低场磁共振无创测温技术

2、中r l 值的测量。肘。与温度的关系及测温时最佳T R 值的选取 等问题，提出了增加修正系数等解决方案，并对方案进行检验，结合仿真计算，验证了解决方案的可行性。 关键词无创测温；纵向弛豫时间；平衡磁化强度；磁共振成像 中图分类号R 3 1 8 0 4文献标识码A文章编号 1 0 0 2 3 2 0 8 ( 2 0 0 8 ) 0 6 0 5 8 7 0 4 N o n - i n v a s i v et h e r m o m e t r yu s i n gl o wf i e l dM R Ib a s e do np r o t o nr e l a x a t i o nt i m

3、eW A N GY m 9 1 一，S O N GT a 0 1 ，W A N G H u t x t a n l ，H UL i l t l 1I n s t i t u t eo fE l e c t r i c a lE n g m e e r i n g ，C h m e s eA c a d e m yo fS c i e n c e s ，B e t t i n g 1 0 0 0 8 0 ；2G r a d u a t e U n i v e r s i t yo fC h m e s eA c a d e m yo fS c i e n c e s ，B e t j t n g

4、 1 0 0 0 4 9 【A b s t r a c t 】 M R Ic a nb eu s e df o rn o n l n v a s l v et e m p e r a t u r et h e r m o m e t r y T h ee q u i l i b n u mm a g n e t i z a t i o nM o a n dt h e s p i n l a t t i c er e l a x a t i o nt i m eT I a r ed e p e n d e n to nt e m p e r a t u r ea n dc a nb e c a

5、l c u l a t e df r o mt h ev a l u e sf o rt h er e p e t i t i o nt i m e T Ra n dt h es i g n a li n t e n s i t i e s T h i sm e t h o di se s p e c i a l l ye f f e c t i v ei nt h el o wf i e l dMR I B a s e do nt h eI n t r o d u c t m no ft h e t h e o r yo fn o n - m v a s i v et h e r m o m

6、 e t r yu s i n gMR I ，t h ep a p e ra n a l y s e st h em e a s u r e m e n to fT I ，t h er e l a t i o n s h i po fM oa n d t e m p e r a t u r e ，a n dt h es e l e c t i o no ft h eb e s tT R ，t h e nt h es o l u t i o n so ft h ep r o b l e m sa r ep r o p o s e d T h ee x p e r i m e n tr e s u

7、 l ta n d s i m u l a t i o nr e s u l ta r ep r e s e n t e dt op r o v et h ef e a s i b i l i t yo ft h ep r o p o s e ds o l u t i o n s 【K e yw o r d s 】n o n m v a s i v et h e r m o m e t r y ；s p m l a t t i cr e l a x a t i o nt i m e ；e q u i l i b r i u mm a g n e t iz a t i o n ；m a g n

8、e t i cr e s o n a n c e i m a g i n g 磁共振无创测温技术是一种新兴的温度测量技 术。鉴于低场开放式磁共振成像( M R I ) 系统价格上 的优势和已经大量应用到介入式治疗中的现状，无 创测温在低场M R I 中将具有更加广阔的应用前景。 目前，磁共振无创测温技术主要有三种基本方法川 测分子扩散系数D ，测质子共振频率的漂移( P R F ， p r o t o nr e s o n a n c ef r e q u e n c y ) 以及测纵向弛豫时间 r 1 。其中测D 的方法已经不常用。测P R F 虽然在 高场M R I 中能够达到较高的测量精

9、度( 1o C ) ，但 是这种方法的温度敏感度与场强成正比，在低场 M R I 中应用有很多困难“。测n 的方法，温度敏感 度在低场中比高场中更高，并且在脂肪中测r 1 的方 法温度准确性要好于测P R F 。“。采用测r 1 的方法 时，利用平衡磁化强度( e q u i l i b r i u mm a g n e t i z a t i o n ) 肘。 作者单位：l中国科学院电丁研究所( 北京1 0 0 0 8 0 ) 2 中国科学院研究生院( 北京1 0 0 0 4 9 ) 作者简介：王莹( 1 9 8 2 一) ，男，硕士研究生 与温度的关系可以增加信号S 的温度敏感度，鉴于 本

10、文中应用的是0 2 3 T 低场磁共振成像系统，因此 采用肘。一r l 的方法。 本文针对低场磁共振无创测温技术中r l 值的 测量，肘。与温度的关系及测温时最佳T R 值的选取 等问题，提出了增加修正系数等解决方案。仿真计 算结果和实验结果的比较分析表明，这些方案是可 行的。 l 磁共振无创测温的原理 1 1 M o 和n 与温度的关系 设样品初始绝对温度值为z ，温度变化后绝对 温度值为巩，单位为K ；温度变化值为丁，单位 为。根据居里定律，在所有组织中，肘。与组织的 绝对温度r 成反比。卜4 。，即M 。( T ) = C T ，C 为与磁 场有关的常数，由此可得肼。与r 的关系列如下：

11、 M o b = M o 。( 1 + T c ) ( 1 ) 万方数据 5 8 8 北京生物医学工程 第2 7 卷 式中，M 。和M 。分别代表温度变化前后样品的 平衡磁化强度。由于非理想条件的影响，需要对理 论公式进行修正，定义k 。为磁化强度修正系数，修 正后公式如下： M 。= M 。，( 1 + 后。t ) ( 2 ) 温度和丁l 之间的关系可由下式来描述“： T 1O Ce x p ( 一E 。k T ) ( 3 ) 式中，E 。是纵向弛豫过程的活化能；k 是玻尔兹 曼常数；T 是绝对温度。温度变化A T 较小时，温度 和r 1 可近似成如下关系“： T 1h = P l r +

12、p 2 丁+ T 1 。 ( 4 ) 式中，温度敏感校正系数P 。和P ：为常数，由组 织类型决定；T 1 。和r 1 。分别表示温度变化前后的 n 值。 当组织温度为r ，采用自旋回波( S E ，s p i n e c h o ) 序列，回波时间T E 较小时，信号强度s 主要反 映的是r 1 的变化，可将横向弛豫时间死忽略，s 可 以表示为： s ( T ) = M 。( r ) 1 一e x p 一T R T 1 ( r ) ( 5 ) 式中，s ( r ) 是温度为r 时的信号强度；T R 是重 复时间；T 1 ( z ) 是温度为? 时的r 1 值；M 。( r ) 是温 度为丁时

13、的平衡磁化强度。 1 2 S E 序列信号数学表达式的修正 目前已有的文献刊发现，主磁场的不均匀性， 射频脉冲的不理想性和噪声的影响等都会造成r l 值测量的误差。在低场中这些非理想情况的影响更 加显著。 由于非理想条件的影响会造成r l 值测量的不 准确，因此有必要对式( 5 ) 进行修正。B e r t s c h3 。采 用了增加修正参数的方法，将其应用到I R 序列中以 补偿非理想的1 8 0 。脉冲造成的影响。本文将这种 方法应用到S E 序列中，用修正系数b 来修正数据受 非理想条件影响所带来的误差。修正后的信号公式 可以表示为： S ( T ) = M o ( 丁) l b (

14、T ) e x p 一T R T 1 ( T ) ( 6 ) 实验中发现，b 与温度有如下关系： b 。= k b T + 6 。( 7 ) 式中，修正系数温度敏感校正常数后。由采用的 成像序列决定；6 。和b 。分别为温度变化前后的修正 系数。 综合起来，温度变化后的信号值s 。与T 的关 系可以表示为： S h = M o h 1 一( k b 7 + b 。) e x p ( 一T R T Ih ) = M ( 1 + | 。A T T c ) 1 一( | | b r + b 。) ( 8 ) e x p 一T R ( p 严+ p 2 7 + T 1 。) 式中，S 。为温度变化7

15、1 后的信号值。该式建 立了5 。和r 之间的非线性关系，温度变化值7 1 可 通过该式求解得出。当r 较小时，非线性简化为 线性关系，可以表示为： T = 一k S h + S o ( 9 ) 2 磁共振无创测温实验 实验采用的样品箱如图1 所示，圆柱形的样品 箱由有机玻璃制成，直径为1 2 0 m m ，外罩保温材料 橡塑海绵，样品箱内由试管架承载着四个5 0 m l 离心 管。离心管固定在试管架上，不与样品箱底部接触， 以保证样品受热均匀。离心管内装满2 5 m m o l L 的C u S O 。溶液，边沿用硫化硅橡胶密封，其中一个 样品管内置热敏电阻温度传感器测温。样品箱通过 两根硅

16、橡胶管与水浴箱相连，外加蠕动泵使蒸馏水 在水浴箱与样品箱之间进行循环，以实现温度的梯 度变化。 图1 实验样品箱 F i g1E x p e r i m e n t a la p p a r a t u s 将样品箱放入由中国科学院电工研究所研制的 0 2 3 T 磁共振成像系统中( O P E NS C A N2 0 0 0 。l E E C A S ，C h i n a ) ，通过进水管和出水管实现水浴箱与 样品箱之间的水循环，水浴箱温度分别选取2 0 、2 5 、 3 0 、3 5 、4 0 、4 5 、5 0 ，当样品箱的温度分别达到上述 各个温度时，关闭蠕动泵停止水循环。蠕动泵循环

17、时间大于1 5 r a i n ，样品受热均匀。由于保温效果良 万方数据 第6 期 基于纵向弛豫时间的低场磁共振无创测温技术的研究 5 8 9 好，样品箱内的温度稳定值( 从抗电磁干扰的温度 频率模式热敏电阻温度传感器读出) 可以在十多分 钟内保持不变，满足该稳定温度下的成像要求。 在每一个稳定温度下对样品管所在的样品箱断 面进行磁共振成像，采用s E 序列，根据样品T l 的 理论值( 约为2 0 0 3 0 0 m s 刘) ，分别选取重复时间 T R = 1 3 5 、2 0 0 、2 5 0 、3 3 0 、4 5 0 、6 5 0 、1 2 0 0 m s ，T E = 2 0 m

18、s ，层厚为1 5 m m ，重复次数N S = 2 ，F O V 为 2 4 0 r a m 2 4 0 r a m ，数据采集矩阵大小为1 2 8 1 2 8 。 除T R 外，其他参数在整个实验期间内保持不变。 3 实验结果与讨论 3 1 T 1 与温度的关系 本文采用的S E 序列由于1 8 0 。射频脉冲的作 用，在自旋回波中心，所有自旋核的相位都得到了重 聚，能够将温度变化引起的相位变化1 列全部消除 掉。 实际测量样品箱的稳定温度为2 0 3 、2 5 9 、 3 0 7 、3 5 4 、3 9 7 、4 5 5 、5 0 4 ，在每一个稳定温度 下得到了7 幅T R 值不同的S

19、 E 序列磁共振图像。 图2 表示T R = 4 5 0 m s 时样品在两个不同温度下的 磁共振图像。其中( a ) 为样品稳定温度是2 0 3 时 的图像，( b ) 为样品稳定温度是5 0 4 时的图像， ( a ) 中样品的信号强度值比( b ) 略强。 针对每一个稳定温度下7 幅T R 值不同的图 像，选取感兴趣区R O I = 1 5 m m 1 5 r a m ，计算数据信 号强度平均值S ，采用原始公式( 5 ) 和采用修正公式 ( 6 ) 对7 组S 和T R 值进行拟合，得到肘。和丁l 值。 图3 是稳定温度为2 0 3 时肘。和r l 的计算 结果，虚线和实线分别表示拟合

20、时采用原始公式 ( 5 ) 和采用修正公式( 6 ) 的结果对比。可以看出，采 用修正公式，其决定系数( r 2 ) 比不采用时大( 在每一 个稳定温度下都得到了这样的结果) ，因此，修正公 式在一定程度上去除了非理想情况带来的影响。 选定r ，为2 9 3 4 5 K ，丁的变化为5 6 、1 0 4 、 1 5 1 、1 9 4 、2 5 2 、3 0 1o C 。根据前面通过式( 6 ) 拟 合得到的不同稳定温度下样品的r l 值，可以得到 r 1 随3 T 的变化曲线( 图4 ) 。实验结果符合式( 4 ) ， 温度敏感校正系数P ，= 0 0 3 1 3 ，P 2 = 3 1 9 8

21、 3 。 根据实验结果，修正系数b 确实也受温度变化 的影响，与温度近似成线性关系，对不同稳定温度下 的修正系数6 进行计算，可以确定式( 7 ) 中k 。= 图2 不同温度下的样品图像 F i g2I m a m si nd i f f e r e n tt e m p e r a t u r e 图3T 1 拟合的对比曲线 F i g3 C o n t r a s t o fT 1b yd i f f e r e n te q u a t m n 一0 0 0 1 7 7 ，b 。= 1 0 9 4 7 。 3 2 M 。与温度的关系 用帆一r l 方法进行温度测量，必须注意到非理 想条件

22、带来的r l 值测量的不准确性，实验结果表 万方数据 5 9 0 - 北京生物医学工程第2 7 卷 图4T I 与A T 的关系图 F J g4R e l a h o n s h l po f7 1 a n dA T 明采用不同的拟合公式得到的丁1 值有一定的差 别，M 。和A T 的关系与理论上也有一定的差别。这 个影响可能是将砣忽略或主磁场的不均匀性造成 的，可以通过添加修正系数后。的方法将这些影响都 放入肘。和A T 的关系中考虑。 本实验中，L 为2 9 3 4 5 K ，r 变化范围为0 3 0 。根据前面通过式( 6 ) 拟合得到的不同稳定温 度下样品的肘。，发现膨。随温度的变化情

23、况与理论 公式( 1 ) 有一定差异，按照式( 2 ) 进行拟合，得到磁 化强度修正系数l | 。= 1 0 6 。可以计算出温度每增 加1 ，M o h M 。约减少0 0 0 3 3 ，即T 。= 2 9 3 4 5 K 时， 肘。h M 。，的温度变化率约为一0 3 3 o C 。 3 3 最佳T R 值的选取 在磁共振无创温度测量中，采用s E 序列时，T R 是成像前唯一需要选择的参数。该参数的选取直接 影响到温度测量的准确性。 式( 8 ) 建立了信号s 。与A T 之间的非线性模 型。该模型中，T R 选择越大时，指数项e x p ( 一T R r l 。) 对n 的变化越不敏感

24、，即S 。对温度的变化也 不敏感，S 。与A T 之间的非线性关系不显著；T R 选 择过小时，S 。与r 之间的非线性关系显著。如何 选取适当的T R 值，使5 。既具有较高的温度敏感 度，又与r 之间具有较好的线性关系，是进行磁共 振无创温度测量时必须解决的问题。 根据式( 8 ) ，对于该信号S 。与r 之间的非线 性模型而言，T R 取1 0 0 1 5 0 0 之间的任意常数时， 在每一个确定的T R 下都可以得到r 在0 3 0 之间变化时对应的信号值s 。，用式( 9 ) 对每一组取 定T R 后的S 。和丁进行拟合，得到J | i ( 1 几可以表 征s 。的温度敏感度) 。1

25、 k 。值随T R ( 范围1 0 0 1 5 0 0 ) 的变化情况如图5 中曲线所表示。7 个数据 点分别表示T R = 1 3 5 、2 0 0 、2 5 0 、3 3 0 、4 5 0 、6 5 0 、 1 2 0 0 m s 时，用式( 9 ) 直接对实验采得的信号S 。与 丁进行拟合得到的l 后值。可以看出拟合实验数 据点的1 k 。值与拟合非线性模型的1 k 。值符合得 很好，因此可以用1 k 。来选择测温时最佳T R 值。 从图中可知，T R 的取值范围在3 3 0 4 5 0 m s 时，比例 系数1 k 。最大，此时信号J s 。的温度敏感度最高。 图5线性模型的斜率l k

26、 与T R 关系曲线 F l g5R e l a t i o n s h i po f1 k 。a n dT Ri nh n e a rm o d e l 用式( 9 ) 对式( 8 ) 取定T R 后的S 。和A T 进行拟 合，决定系数r 2 可以表征式( 9 ) 对式( 8 ) 线性拟合 的好坏。r 2 与T R ( 范围1 0 0 1 5 0 0 ) 的关系如图6 所示，T R 5 0 0 时，r 2 随着T R 的增加而增加。T R = 3 3 0 m s 时，r 2 = 0 9 2 4 1 9 ，T R = 4 5 0 m s 时，r 2 = 0 9 8 1 5 4 ，表明线性关系

27、很好。当选定T R = 4 5 0 m s 时，有S h = 一2 3 9 5 4 A T + 2 9 3 6 6 4 3 。 1 0 9 5 09 08 5 0 8 、O7 5 0 7 06 5 0 6 05 5 0 5 图6 线性模型可决系数r 2 与T R 关系曲线 F i g6R e l a t i o n s h i po fr 2a n dT Ri nh n e a rm o d e l 综上，根据得到的l k 。、r 2 与T R 的关系，T R 值 的取值范围在3 3 0 4 5 0 m s 之间时，S 。既具有较高 ( 下转第5 9 5 页) 万方数据 第6 期种基于完全匹配

28、和分词匹配的混合分词匹配算法 5 9 5 3 4 5 6 7 8 9 1 0 1 1 R a n t a n e nA ，M i e h k a m e nT R o u s uJ ，e ta 1 P J a n n m go p t i m a l m e a s u r e m e n t so f l s o t o p o m e r d l s t r t h u t l o n sf o re s t i m a t i o no f m e t a b o l i cf l u x e s O x f o r dU m vP r e s s 2 0 0 6 ：l1 9 8 一1

29、2 0 6 G a a s t e r l a n dT ，S e l k o vE R e c o n s t r u c t , o no fm e t a b o l i cn e t w o r k s u s i n gi n c o m p l e t ei n f o r m a t i o n P r o cI n tC o n “ I n t e l lS y s tM o lB l o l ， 1 9 9 5 3 ：1 2 7 3 5 P a l s s o nB O ，S y s t e m sB i o l o g y P r o p e m e so fR e c o

30、 n s t r u c t e d N e t w o r k s N e wY o r k N Y U S A ：C a m b r i d g eU n l v e m f l yP r e s s ， 2 0 0 6 B o r o d i n aI ，K r a b b e n P ，N i e l s e nJ G e n o m e s c a l ea n a l y s i so f S t r e p t o m y c e se o e h c o l o rA 3 ( 2 ) m e t a b o l i s m G e n o m eR e s ，2 0 0 5 ，

31、 1 5 ( 6 ) ：8 2 0 9 D u a r t eN C ，B e c k e rS A ，J a m s h l d lN ，e ta 1 G l o h a lr e c o n s t r u c t i o n o ft h eh u m a nm e t a b o l i cn e t w o r kb a s e do ng e n o m l ca n db l b h o m w d a t a ，P r o cN a t lA c a dS o lU S A 。2 0 0 7 1 0 4 ( 6 ) ：1 7 7 7 8 2 F e i s tA M ，S c

32、h o h e nJ C ，P a l s s o n B O ， e ta 1 M o d e l i n g m e t h a n o g e u e s l sw i t h ag e n o m e s c a l em e t a b o l i cr e c o n s t r u c t i o no f M e t h a n o s a r c i n ah a r k e n M o lS y s tB i o l ，2 0 0 6 ，2 ：2 0 0 6 4 F o r s t e rJ 。F a m l hI F uP ，e la 1 G e n o m e s c

33、a l er e c o n s t r u c t i o no f t h eS a c c h a r o m y c e sc o r e v l s l a em e t a b o l i cn e t w o r k G e n o m eR e s ， 2 0 0 3 ，1 3 ( 2 ) ：2 4 4 5 3 B e c k e tS A P a l s s o nB O G e n o m e s c a l er e c o n s t r u c t i o no ft h e m e t a b o l i cn e t w o r ki nS t a p h y l

34、 o c o c c u sa u r e u sN 3 1 5 ：a nl n l t l a ld r a f t t ot h et w o d i m e n s i o n a la n n o t a t i o n B M CM i c r o b l o i ，2 0 0 5 ，5 ( 1 ) ： 8 F e t s tA M ，H e n r yC S ，R e e d 儿，e ta 1 Ag e n o m e - s c a l em e t a b o l i c 1 2 1 3 1 4 1 5 1 6 【1 7 j 1 8 1 9 【2 0 r e c o n s t

35、 r u c t i o nf o rE s e h e n c h i ac o l 】K 一】2M G l 6 5 5t h a ta c c o n n t sf o r I2 6 0O R F sa n d t h e r m o d y n a m i cm f o r r n a t l o n M o lS y s tB x o l ， 2 0 0 7 3 ：1 2 1 S t o e s s e rG ，T u hM A ，L o p e zR ，e ta 1 T h eE M B LN u c l e o t l d e S e q u e n c eD a t a b a

36、s e O x f o r dU m vP r e s s ：1 8 2 4 B a l r o c hA T h eE N Z Y M Ed a t a b a s ei n2 0 0 0 N u c l e i cA c i d s R e s e a r c h 2 0 0 0 2 8 ( 1 ) ：3 0 4 3 0 5 K a n e h m aM ，T h eK E G Gd a t a b a s e N o v a r h sF o u n dS y m p 2 0 0 2 ， 2 4 7 ：9 1 一1 0 1 K a r pP D ，R i l e yM ，S a l e

37、rM ，e ta 1 T h eE c o C y ca n dM e t a C y e d a t a b a s e s N u c l e , cA c i d sR e s e a r c h 2 0 0 0 ，2 8 ( 1 ) ：5 6 5 9 P a l e yS M ，K a r pP D E v a l u a t i o no fc o m p u t a t i o n a lm e t a b o l i c - p a t h w a yp r e d i c t i o n sf o r H e h c o b a c t e rp y l o n B i o m

38、 f o r m a t i c s ， 2 0 0 2 ，1 8 ( 5 ) ：7 1 5 7 2 4 B l l o f s k yH S ，B u r k sC ，F l c k e t tJ W ，e ta l ，T h eG e n B a n kg e n e t m s e q u e n c ed a t a b a n k O x f o r d ：U n i vP r e s s ，1 9 8 6 S c h o m b u r gI ，C h a n gA ，E b e l m gC ，e ta lB R E N D A ，t h ee n z y m e d a t a

39、 b a s e ：u p d a t e sa n dm a j o rf l e wd e v e l o p m e n t s N u c l e i cA c i d s R e s e a r c h 2 0 0 4 K a r pP D O u z o u m sC A M o o r e K o c h l a e sC ，e ta 1 E x p a n s i o no f t h eB m C y cc o l l e c t i o no fp a t h w 8 y g e n o m ed a t a b a s e st o1 6 0 g e n o m e sN

40、 u c l e i cA c i d sR e s e a r c h 2 0 0 5 ，3 3 ( 1 9 ) ：6 0 8 3 6 0 8 9 K a r pP D P a l e yS ，R o m e r oP T h eP a t h w a yT o o l ss o f t w a r e B l m n f o r m a t i c s ，2 0 0 2 ，1 8 ( S u p p l1 ) ：$ 2 2 5 一$ 2 3 2 ( 2 0 0 8 0 l 0 4 收稿，2 0 0 8 一0 3 一1 9 惨回) ( 上接第5 9 0 页) 的温度敏感度又和r 有比较好的线

41、性关系。 在3 3 0 4 5 0 m s 范围内选择最佳的T R 值时，成 像时间为l 一2 m i n ，这对于实时的温度测量而言，时 间有些长，以后的研究中可以把本文的方法应用到 F S E 序列中以缩短成像时间。 4结论 在0 2 3 T 磁共振成像系统中，由于受到很多非 理想条件的限制，准确测量温度比较困难。本文采 用了修正S E 序列参数的方法，并对测温时最佳T R 值的选取问题做了研究，取得了比较令人满意的结 果。该方法应用于低场M R I 测温中是可行的。 参考文献 1 】G e r m a mD C h e v a l h e rP t h e r a p y M A G M

42、 A ，2 0 0 1 ， M Rm o m t o n n go ft o m o u rt h e r m a l 1 3 ：4 7 5 9 2 3 4 5 6 7 G e r m a mD ，V a h a a l aE M Rt e m p e r a t u r em e a s u r e m e n ti nh v e r t i s s u ea t02 3 Tw x t has t e a d y s t a t ef r e ep r e c e s s i o ns e q u e n c e M R M 2 0 0 2 4 7 ：9 4 7 9 4 9 B e r t

43、s c hF M a t t n e rJ N o n l n v a s l v et e m p e r a t u r em a p p i n gu s i n g M R I ：c o m p a r i s o no ft w om e t h o d sb a s e do nc h e m i c a ls h i f ta n dT 1 r e l a x a t m n M R I 1 9 9 8 1 6 ：3 9 3 4 0 4 G u h e k mD H ，G o r eJ C T e m p e r a t u r ed e p e n d e n c eo fn

44、u c l e a r m a g n e t i z a t i o na n dr e l a x a t i o n J M R 。2 0 0 5 ，1 7 2 ：1 3 3 1 4 1 N e l s o nT R ，T u n gS M T e m p e r a t u r ed e p e n d e n c eo fp r o t o n r e l a x a t m nt i m e si nV l l O M R I 1 9 8 7 5 ：1 8 9 1 9 9 G e r m a l nD ，c h e v a l i e rP M Rm o m t o n n go

45、fl a s e r i n d u c e dl e s l o n s o ft h el i v e ri n V l V Ol i tal o w f i e l do p e nm a g n e t ：t e r n p e r a t u r e m a p p i n ga n dl e s i o n s l z ep r e d i c t i o n J M R I ，2 0 0 1 ，1 3 ：4 2 4 9 K j a e rL T h o m s o nC E v a l u a t l o no fr e l a x a t i o n t i m em e a s u r e m e n t s b ym a g n e t i c r e s o n a n c ei m a g i n g A c t aR s d l o l o g l e a ，1 9 9 7 。2 8 ： 3 4 5 3 5 1 ( 2 0 0 8 一0 2 一1 9 收稿) 万方数据