发明名称 |
一种超快速分段式单次激发水脂分离方法 |
摘要 |
本发明公开了一种超快速分段式单次激发水脂分离方法,通过超快速化学位移成像序列、信号分析模型、回波幅度衰减补偿方法和相位快速解缠及精确校正技术用于在磁共振成像系统上实现水脂分离增强扫描和薄层水脂分离成像;可大幅度提高扫描效率的同时通过相位校正消除图像伪影;二维版首次使得化学位移成像可以在增强扫描中实现临床应用,也可用于婴儿和多动症患者扫描。 |
申请公布号 |
CN106324537A |
申请公布日期 |
2017.01.11 |
申请号 |
CN201610658204.3 |
申请日期 |
2016.08.12 |
申请人 |
江亲瑜 |
发明人 |
江亲瑜 |
分类号 |
G01R33/561(2006.01)I |
主分类号 |
G01R33/561(2006.01)I |
代理机构 |
北京方安思达知识产权代理有限公司 11472 |
代理人 |
徐淑东;王宇杨 |
主权项 |
一种超快速分段式单次激发水脂分离方法,其特征在于,所述超快速分段式单次激发水脂分离方法包括以下步骤:步骤S100:根据预扫描需要,通过控制台在工控机序列发生器上加载表观横向弛豫时间<img file="FDA0001076437790000011.GIF" wi="118" he="70" />测试序列;步骤S200:通过控制台在序列发生器上加载二维化学位移成像序列模块;在二维梯度回波序列基础上利用2N个正负极性交替的频率编码梯度形成由N组水脂同相位回波和反相位回波交替构成的梯度回波串;N个同相回波再分别进行梯度相位编码,与之相应地k空间划分为N个区域,第一个回波填充中心区域,下一个回波依次填充紧邻区域,并采用部分傅立叶方式重建出同相图;N个反相回波串分别进行梯度相位编码,与之相应地k空间划分为N个区域,第一个回波填充中心区域,下一个回波依次填充紧邻区域,并采用部分傅立叶方式重建出反相图;或在三维梯度回波序列基础上利用2N个正负极性交替的频率编码梯度形成由N个水脂同相位回波和N个反相位回波交替构成的梯度回波串;N个同相回波分别进行梯度相位编码,与之相应地k空间划分为N个区域,第一个回波填充中心区域,下一个回波依次填充紧邻区域,然后进行三维离散逆傅立叶变换产生一系列薄层同相图;N个反相回波分别进行梯度相位编码,与之相应地k空间划分为N个区域,第一个回波填充中心区域,下一个回波依次填充紧邻区域,然后进行三维离散逆傅立叶变换产生一系列薄层同相图;步骤S300:利用在空间频率域和图像域对梯度回波串建立的通用理论模型,校正射频脉冲激发到第一个回波形成的TE期间磁场不均匀性对回波相位和幅度的影响:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>2</mn></mrow><mi>t</mi></msubsup><mo>=</mo><munderover><mo>∫</mo><mrow><mo>-</mo><mi>∞</mi></mrow><mrow><mo>+</mo><mi>∞</mi></mrow></munderover><mo>{</mo><mrow><mo>(</mo><msub><mi>M</mi><mi>w</mi></msub><mo>(</mo><mi>r</mi><mo>)</mo><mo>+</mo><msub><mi>M</mi><mi>f</mi></msub><mo>(</mo><mi>r</mi><mo>)</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mn>2</mn><mi>π</mi><mi>Δ</mi><mi>f</mi><mo>·</mo><mo>[</mo><mi>T</mi><mi>E</mi><mo>+</mo><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mi>Δ</mi><mi>τ</mi><mo>]</mo></mrow></msup><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>Δ</mi><mi>τ</mi><mo>/</mo><msubsup><mi>T</mi><mn>2</mn><mo>*</mo></msubsup></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mrow><mo>[</mo><mrow><msub><mi>φ</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>φ</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></msup><mo>}</mo><msup><mi>d</mi><mn>3</mn></msup><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001076437790000012.GIF" wi="1765" he="135" /></maths><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>I</mi><mrow><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mi>t</mi></msubsup><mo>=</mo><munderover><mo>∫</mo><mrow><mo>-</mo><mi>∞</mi></mrow><mrow><mo>+</mo><mi>∞</mi></mrow></munderover><mo>{</mo><mrow><mo>(</mo><msub><mi>M</mi><mi>w</mi></msub><mo>(</mo><mi>r</mi><mo>)</mo><mo>+</mo><msub><mi>M</mi><mi>f</mi></msub><mo>(</mo><mi>r</mi><mo>)</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mn>2</mn><mi>π</mi><mi>Δ</mi><mi>f</mi><mo>·</mo><mo>[</mo><mi>T</mi><mi>E</mi><mo>+</mo><mrow><mo>(</mo><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>Δ</mi><mi>τ</mi><mo>]</mo></mrow></msup><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mn>2</mn><mi>N</mi><mo>·</mo><mi>Δ</mi><mi>τ</mi><mo>/</mo><msubsup><mi>T</mi><mn>2</mn><mo>*</mo></msubsup></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mrow><mo>[</mo><mrow><msub><mi>φ</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mi>N</mi><mo>·</mo><mi>φ</mi><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></msup><mo>}</mo><msup><mi>d</mi><mn>3</mn></msup><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001076437790000013.GIF" wi="1790" he="143" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>S</mi><mn>0</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>w</mi></msub><mo>-</mo><msub><mi>S</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>·</mo><msup><mi>e</mi><mrow><mi>i</mi><mover><mi>φ</mi><mo>‾</mo></mover></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001076437790000014.GIF" wi="1149" he="79" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>S</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>w</mi></msub><mo>+</mo><msub><mi>S</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>·</mo><mi>A</mi><mo>·</mo><msup><mi>e</mi><mrow><mi>i</mi><mrow><mo>(</mo><mover><mi>φ</mi><mo>‾</mo></mover><mo>+</mo><mi>φ</mi><mo>)</mo></mrow></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001076437790000015.GIF" wi="1157" he="93" /></maths>其中,t指当前的相位编码循环步,N在1到2L范围取整数,φ<sub>0</sub>(r)是质子磁化矢量的初始相位,φ(r)是主要由磁场不均匀效应产生的相位误差,<img file="FDA0001076437790000021.GIF" wi="54" he="70" />项用于对回波幅度进行优化补偿;S<sub>w</sub>和S<sub>f</sub>分别表示图像像素内水和脂肪成份,它们与M<sub>w</sub>和M<sub>f</sub>成正比,<img file="FDA0001076437790000022.GIF" wi="35" he="68" />包括反相图的初始相位以及在TE和(2N‑2)·Δτ期间场不均匀性、涡流场和麦克斯韦场产生的相位误差,<img file="FDA0001076437790000023.GIF" wi="37" he="70" />主要由第一个反相回波(即k<sub>y</sub>=0线)在TE期间的相位误差决定;M<sub>w</sub>(r)和M<sub>f</sub>(r)分别表示成像区域内水和脂肪的质子磁化矢量,在三维空间编码条件下可描述为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>M</mi><mi>w</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>M</mi><mi>w</mi><mn>0</mn></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>γ</mi><msubsup><mo>∫</mo><mn>0</mn><mi>t</mi></msubsup><mo>[</mo><mi>x</mi><mo>·</mo><msub><mi>G</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>+</mo><mi>y</mi><mo>·</mo><msub><mi>G</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>+</mo><mi>z</mi><mo>·</mo><msub><mi>G</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>]</mo><mi>d</mi><mi>τ</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001076437790000024.GIF" wi="1189" he="100" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>M</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>M</mi><mi>f</mi><mn>0</mn></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>γ</mi><msubsup><mo>∫</mo><mn>0</mn><mi>t</mi></msubsup><mo>[</mo><mi>x</mi><mo>·</mo><msub><mi>G</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>+</mo><mi>y</mi><mo>·</mo><msub><mi>G</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>+</mo><mi>z</mi><mo>·</mo><msub><mi>G</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow><mo>]</mo><mi>d</mi><mi>τ</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001076437790000025.GIF" wi="1238" he="103" /></maths>步骤S400:水像S<sub>w</sub>和脂肪像S<sub>f</sub>按照下式计算产生:<maths num="0007"><math><![CDATA[<mrow><msub><mi>S</mi><mi>w</mi></msub><mo>=</mo><mfrac><mrow><mi>κ</mi><mo>·</mo><mo>|</mo><msubsup><mi>S</mi><mn>0</mn><mo>`</mo></msubsup><mo>|</mo><mo>+</mo><mo>|</mo><msubsup><mi>S</mi><mn>1</mn><mo>`</mo></msubsup><mo>|</mo><mo>/</mo><mi>A</mi></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0001076437790000026.GIF" wi="390" he="142" /></maths><maths num="0008"><math><![CDATA[<mrow><msub><mi>S</mi><mi>f</mi></msub><mo>=</mo><mfrac><mrow><mo>|</mo><msubsup><mi>S</mi><mn>1</mn><mo>`</mo></msubsup><mo>|</mo><mo>/</mo><mi>A</mi><mo>-</mo><mi>κ</mi><mo>·</mo><mo>|</mo><msubsup><mi>S</mi><mn>0</mn><mo>`</mo></msubsup><mo>|</mo></mrow><mn>2</mn></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0001076437790000027.GIF" wi="434" he="142" /></maths> |
地址 |
116000 辽宁省大连市沙河口区中长东五街10号 |