发明名称 一种磁共振测温的方法和装置
摘要 本发明公开了一种磁共振测温方法,包括:建立以脂肪为参照物的磁共振温度模型;在热疗过程中使用多回波梯度回波序列进行图像采集;利用磁共振兼容的热电偶或热敏光纤在首次热疗或预扫描阶段进行实际温度测量用以标定;根据所述温度模型,利用扩展Prony算法对数据进行拟合,获取水和油的信号幅度、初始相位、横向弛豫率和扫描仪解调后的共振频率;根据水和脂肪质子的共振频率,获取油水化学位移;根据所述油水化学位移和实测温度值,获取温度系数α和系数β;根据所述油水化学位移、温度系数α和系数β,获取绝对温度。本发明还公开了一种磁共振测温的装置。本发明能够快速获取高分辨率的绝对温度图像,并大幅减轻运动、形变、场漂、脂肪等因素对测温准确度的影响。
申请公布号 CN101507603B 申请公布日期 2012.11.14
申请号 CN200910001229.6 申请日期 2009.01.04
申请人 清华大学;西门子迈迪特(深圳)磁共振有限公司 发明人 潘昕怿;李诚;应葵;白逸轩
分类号 A61B5/01(2006.01)I;A61B5/055(2006.01)I;A61N5/00(2006.01)I;A61N7/02(2006.01)I 主分类号 A61B5/01(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 张国良
主权项 1.一种磁共振测温的装置,其特征在于,所述装置包括:温度模型建立单元,用于根据公式<maths num="0001"><![CDATA[<math><mrow><mi>s</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>water</mi><mo>,</mo><mi>fat</mi></mrow></munder><msub><mi>&rho;</mi><mi>i</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&phi;</mi><mi>i</mi></msub></mrow></msup><msup><mi>e</mi><mrow><mrow><mo>(</mo><mo>-</mo><msubsup><mi>R</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow><mo>*</mo></msubsup><mo>+</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>t</mi><mi>n</mi></msub></mrow></msup><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></math>]]></maths>建立以脂肪为参照物的磁共振温度模型,其中s(t<sub>n</sub>)是采得的磁共振信号,ρ<sub>i</sub>是水和油的质子密度,φ<sub>i</sub>是初始相位,<img file="FDA00001703473000012.GIF" wi="52" he="54" />是横向弛豫率,f<sub>i</sub>是扫描仪解调后的共振频率,w(n)是复高斯白噪声,t<sub>n</sub>为回波时间,n=0,1,2,...N-1,N为单次测温所需图像数;图像采集单元,用于在热疗过程中使用多回波梯度回波序列进行图像采集;实际温度测量单元,用于利用磁共振兼容的热电偶或热敏光纤在首次热疗或预扫描阶段进行实际温度测量;共振频率获取单元,用于根据所述温度模型,利用公式<maths num="0002"><![CDATA[<math><mrow><mi>s</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>water</mi><mo>,</mo><mi>fat</mi></mrow></munder><msub><mi>&rho;</mi><mi>i</mi></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>T</mi><msub><mi>E</mi><mn>0</mn></msub><msubsup><mi>R</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow><mo>*</mo></msubsup><mi></mi></mrow></msup><msup><mi>e</mi><mrow><mi>j</mi><mrow><mo>(</mo><msub><mi>&phi;</mi><mi>i</mi></msub><mo>+</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>i</mi></msub><msub><mi>TE</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><msup><mi>e</mi><mrow><mrow><mo>(</mo><mo>-</mo><msubsup><mi>R</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow><mo>*</mo></msubsup><mo>+</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>&Delta;TEn</mi></mrow></msup><mo>,</mo></mrow></math>]]></maths>在回波时间间隔相同时利用扩展Prony算法对数据进行拟合,获取水和油的质子密度ρ<sub>i</sub>、初始相位φ<sub>i</sub>、横向弛豫率<img file="FDA00001703473000014.GIF" wi="61" he="61" />和扫描仪解调后的共振频率f<sub>i</sub>,其中TE<sub>0</sub>是多回波梯度回波序列的首个回波时间,ΔTE是回波间隔,S(n)是采得的磁共振信号,n=0,1,2,...N-1,N为单次测温所需图像数;油水化学位移获取单元,用于根据水和脂肪质子的共振频率f<sub>water</sub>和f<sub>fat</sub>,利用公式<img file="FDA00001703473000015.GIF" wi="476" he="135" />获取油水化学位移<img file="FDA00001703473000016.GIF" wi="180" he="59" />其中<img file="FDA00001703473000017.GIF" wi="48" he="42" />是氢质子旋磁比,B<sub>0</sub>是主磁场场强;系数获取单元,用于根据所述油水化学位移<img file="FDA00001703473000021.GIF" wi="149" he="59" />和首次热疗或预扫描阶段得到的实测温度值T<sub>s</sub>,利用公式<maths num="0003"><![CDATA[<math><mrow><msub><mi>&delta;</mi><mrow><msub><mi>H</mi><mn>2</mn></msub><mi>O</mi><mo>-</mo><msub><mi>CH</mi><mn>2</mn></msub></mrow></msub><mo>=</mo><mi>&alpha;</mi><msub><mi>T</mi><mi>s</mi></msub><mo>+</mo><mi>&beta;</mi><mo>,</mo></mrow></math>]]></maths>获取温度系数α和系数β;绝对温度获取单元,用于根据所述油水化学位移<img file="FDA00001703473000023.GIF" wi="184" he="60" />温度系数α和系数β,利用公式<maths num="0004"><![CDATA[<math><mrow><mi>T</mi><mo>=</mo><mfrac><mrow><msub><mi>&delta;</mi><mrow><msub><mi>H</mi><mn>2</mn></msub><mi>O</mi><mo>-</mo><msub><mi>CH</mi><mn>2</mn></msub></mrow></msub><mo>-</mo><mi>&beta;</mi></mrow><mi>&alpha;</mi></mfrac></mrow></math>]]></maths>获取整个热疗过程中病灶区域的绝对温度T。
地址 100084 北京市海淀区清华园北京100084-82信箱