专利

申请号

CN201410057539.0

专利链接

http://www.soopat.com/Patent/201410057539

专利名称

基于单扫描正交时空编码磁共振成像的图像畸变校正方法

英文名称

申请日

2014-02-20

主分类号

G01R33/565

分类号

G01R33/565

申请人

厦门大学

公告号

CN103885017B

当前法律状态

授权

专利有效性

有效

法律状态

法律状态公告日：20140625;法律状态：公开;描述信息：公开;法律状态公告日：20140716;法律状态：实质审查的生效;描述信息：实质审查的生效IPC(主分类):G01R 33/565;  申请日:20140220;法律状态公告日：20170222;法律状态：授权;描述信息：授权;

专利类型

发明授权

公告日

2017-02-22

专利代理机构

厦门市首创君合专利事务所有限公司 35204

代理人

张松亭

地址

361000 福建省厦门市思明南路４２２号

优先权

国省代码

中国

摘要

本发明公开一种基于单扫描正交时空编码磁共振成像的图像畸变校正方法。该方法通过90度和180度的线性扫频脉冲（chirp脉冲）联合正交分布的空间编码梯度，在激发阶段使空间内的质子自旋获得一个和空间位置相关的二次相位；在采样阶段，通过快速切换的正负梯度回波链进行数据采样，超快速地获得具有T2*加权的空间域磁共振成像数据；最后通过高分辨重建算法和图像畸变校正方法对畸变的磁共振图像进行重建校正，从而获得高分辨率的高品质磁共振图像。本发明所提出的基于正交时空编码磁共振成像的图像畸变校正方法，不仅不需要额外的参考扫描，根据自身的数据就可以进行畸变校正，并且还具有克服局部不均匀场的能力，即使在极强的不均匀场下仍能正常工作。

主权项

基于单扫描正交时空编码磁共振成像的图像畸变校正方法，分为以下步骤进行操作：(1)首先对成像物体进行感兴趣区域定位，然后进行调谐、自动匀场、功率和频率校正；(2)用测量线性扫频脉冲功率的脉冲序列依次测量并记录90度和180度线性扫频脉冲(chirp脉冲)的功率；(3)在核磁共振成像仪上，导入事先编译好的单扫描正交时空编码序列；并通过事先定义的脉冲变量调入所需的90度和180度线性扫频脉冲(chirp脉冲)，将步骤(2)中所测量的90度和180度线性扫频脉冲(chirp脉冲)的功率值赋值到所对应的功率变量中；所述单扫描正交时空编码序列的结构依次为：高扫频率的90度线性扫频脉冲、TE₁/2、低扫频率的180度线性扫频脉冲、TE₁/2、TE₂/2、180度的sinc脉冲、TE₂/2、采样回波链；所述高扫频率的90度线性扫频脉冲结合90度的空间编码梯度G₉₀对低带宽维，即y方向进行空间编码，重聚梯度G_er紧接着作用在所述空间编码梯度G₉₀之后,其梯度面积为所述空间编码梯度G₉₀面积的一半，正负与所述空间编码梯度G₉₀相反；所述低扫频率的180度线性扫频脉冲结合180度的空间编码梯度G₁₈₀对高带宽维，即x方向进行空间编码；所述180度的sinc脉冲和层选梯度G_ss进行层选；所述低扫频率的180度线性扫频脉冲前后分别施加了回波延时TE₁/2，所述180度的sinc脉冲前后分别施加了回波延时TE₂/2,所述低扫频率的180度线性扫频脉冲与180度的sinc脉冲前后都有x,y,z三个方向的破坏梯度作用；所述采样回波链是由分别作用在x,y方向的梯度链组成；x方向的梯度链是由一系列正负切换的梯度构成，且每个梯度的面积是所述空间编码梯度G₁₈₀面积的二倍；y方向的梯度链是由一系列大小相等的“blips”梯度构成，且所述“blips”梯度的总面积和等于所述空间编码梯度G₉₀的面积，正负和所述空间编码梯度G₉₀一致；在所述采样回波链之前，x和y方向分别施加了重聚梯度G_ror和G_ar，所述G_ror的面积是x方向第一个梯度面积的一半，方向与之相反；所述G_ar的面积是所有所述“blips”梯度的总面积的一半，方向与所述“blips”梯度相反；(4)执行步骤(3)设置好的单扫描正交时空编码序列，进行数据采样，数据采样完成后，执行下一步骤，否则继续采样；从而得到正交时空编码磁共振数据；(5)当数据全部采样完成后，首先对所述数据进行高分辨重建，得到高分辨率的正交时空编码磁共振数据；并对所述高分辨率的正交时空编码磁共振数据进行相位处理，得到畸变校正所需的高分辨率的正交时空编码磁共振相位图；所述高分辨率重建的步骤是将获得的二维单扫描正交时空编码数据先进行相位平滑处理，然后进行二维的插值，最后逐维进行高分辨率重建；(6)对步骤(5)得到的高分辨率的正交时空编码磁共振相位图进行解缠绕处理，得到场分布的场图；(7)对步骤(6)得到的场图进行平面拟合处理，得到各阶不均匀场磁化率梯度值；所述平面拟合处理的公式为：φ_inh(x,y,tt)＝Ω₀·(tt+TE)+γg_inx(tt+TE)·x+γg_iny(tt+TE)·y+γg_inx2(tt+TE)·x²+γg_inxy(tt+TE)·xy+γg_iny2(tt+TE)·y²其中G_inh＝[g₀,g_inx,g_iny,g_inx2,g_inxy,g_iny2]^T为等效的不均匀场磁化率梯度矢量，tt为从所述180度的sinc脉冲中心到采样时间的时间距离，x”为x方向的空间位置，“y”为y方向的空间位置，TE＝(T₉₀-TE₂)/2,φ_inh(x,y,tt)为解缠绕后的二维相位图；其中T₉₀为90度线性扫频脉冲的时间；(8)将步骤(7)得到的各阶不均匀场磁化率梯度值，代入步骤(5)的高分辨重建中，得到经过畸变校正的高分辨率的正交时空编码磁共振数据，然后重复步骤(5)和(6)得到场图，依据校正是否成功的标准，做以下判断处理：如果不满足校正成功标准则循环进行(5-8)步骤，直至校正成功；如果满足校正成功的标准，则得到能高分辨率的、能正确反映空间信息的高分辨率的正交时空编码磁共振相位图。

发明人

蔡淑惠; 李敬; 蔡聪波; 陈林; 陈忠

inpadoc同族

厦大机构

更新时间

20170305

英文摘要