技术领域
[0001] 本
发明涉及磁共振技术领域，尤其涉及一种三维多参数加权磁共振成像的同步采集和校准方法。
背景技术
[0002] 临床上横向弛豫时间(T2)加权压脂图像有利于提高人体解剖细节的
对比度和
分辨率，从而提高
疾病的确诊率，常常与T2加权图像配合使用。但是，反转恢复法压脂会同时降低其它人体组织和病灶的图像
信噪比，而
频率选择性压脂易于受到
磁场不均匀性和射频场不均匀性的干扰，正因为如此，基于
水脂化学位移的
相位编码水脂分离成像技术作为一种高级成像技术是临床上实现T2加权压脂成像的最佳选择，具有信噪比高和解剖结构细节显示清晰等显著优点，在临床上可用于人体各部位的疾病诊断并提高确诊率，尤其为临床上具有独特诊断价值的磁共振水成像提供了一种更好的技术实现方式。不仅如此，水脂分离成像技术可与纵向弛豫时间(T1)加权和质子
密度(PD)加权方式结合实现T1或PD加权压脂成像，同样具有信噪比高和解剖结构细节显示清晰等优点，并为脂肪变性
鉴别和脂肪定量分析提供了有效的临床诊断方法。
[0003] 近年来基于化学位移成像序列的二维T2加权压脂成像技术在美国GE和德国Siemens高场成像系统上已初步获得临床应用，该技术是以快速自旋回波序列(FSE)为
基础对每个回波一次进行三个化学位移相位编码(-π，0，+π)构成三点Dixon回波组，每次分段激发的回波串的长度是FSE的三倍，并在图像重建时通过线性相位校正方式消除磁场不均匀性和双极性梯度极性切换引起的相位误差，可以获得
质量合格的图像，还通过非Dixon相位编码方式进行不对称采集来避免相位缠绕引起的水脂分离错误。然而，在
硬件性能不足够优异的成像系统上早期的
数据采集和处理方式缺乏普遍适用性，双极性梯度的使用可以导致三点Dixon回波组之间的化学位移相位与预设值偏差较大，回波幅度
波形的不对称性效应更为显著，场漂效应引起选层误差并导致图像模糊，
涡流场的长时间指数项还可以引起高阶相位误差，这些误差不能通过180°重聚脉冲和线性相位校正等方式消除；尤其是，在低场强条件下由于FSE的回波时间和化学位移编码回波的间隔时间较长且磁场均匀性不足够高，相位缠绕通常不可避免，T2*弛豫效应也不可忽略，从而导致水脂分离不完全；另外，早期的二维成像方式在梯度强度不高时难以获得选层方向的高分辨图像，而三维成像方式容易受到磁体和射频系统性能不稳定的影响。
发明内容
[0004] 为了解决上述难题，在此本发明基于单次激发化学位移成像序列提出一种三维T2加权图像和水脂分离图像同步采集方法，并扩展到T1加权和PD加权图像，一次扫描中最大限度获得多幅图像，即，同
相图、反相图、脂肪像、压脂水像、常规T2加权(或T1加权/PD加权)图像以及T2*加权图像，明显节省了临床扫描时间并增加了临床扫描方案的可选择性；并且，本发明进一步配备有预扫描自动校准功能和实时自动校正功能，可采用对称或不对称回波采集模式，对MRI系统的硬件性能的依赖性更小。
[0005] 本发明提供了一种三维多参数加权磁共振成像的同步采集和校准方法，在常规三维快速自旋回波基础上，射频脉冲波形选用最小相位SLR脉冲，
采样期施加三个极性交替的频率编码梯度，其中梯度G1的积分面积是梯度G2的一半，二者极性相反，且梯度G2、G3和G4的积分面积相等，G3和G2的极性相反，ΔG2和ΔG4为补偿梯度，每次射频激发采集多组回波，每组回波包括两个梯度回波和一个自旋回波，回波
顶点之间的时间间隔Δτ设置为1/Δf/n，其中Δf为水脂共振频率差，n一般选2、3或4，然后，将上述成像序列编译后加载到MRI成像系统的谱仪上，并按照工作流程进行数据采集和处理。三维/二维快速多回波水脂分离序列及其
信号调试方法、预扫描方法、扫描方法、数据预处理方法和图像重建方法可实现一次扫描获得水脂分离图像和T2加权(T1加权或PD加权)图像。
[0006] 优选的，三维/二维快速多回波水脂分离序列，以三维/二维快速自旋回波为基础，射频脉冲选用具有均匀激发特征的最小相位SLR脉冲，其时带积优先为8，脉宽优选为2ms，带内和带外纹波系数均不大于0.5％，并在信号采集期施加三个极性交替的读梯度，其信号调试方法是在没有相位编码条件下通过实时调试模式采集信号并通过调试频率编码方向的预备梯度(G1)和补偿梯度(ΔG2和ΔG4)校准回波串中每组三个回波的幅度和相位直到回波峰的时间间隔为1/Δf/n(n一般为2、3或4)。
[0007] 优选的，按下述方式采集参考数据并计算相位误差：
[0008] (1)在所有相位编码梯度关闭条件下，设置相位编码步数DIM2为回波链长度ETL，ETL在T1加权扫描时一般在3至5范围取值，在T2加权扫描时一般在8至24范围取值，采集ETL个回波组，每组包括A、B和C三个回波，其回波峰相位分别为-π、0和+π，并对每组的回波B取复数共轭和
时间反演，再沿频率编码方向对每组回波进行离散逆
傅立叶变换，得到一系列混合空间矩阵元PiA、PiB和PiC，这里i表示分段激发的各回波编号，取值范围为1至ETL；
[0009] (2)在所有相位编码梯度关闭且频率编码梯度极性反转条件下，设置DIM2＝ETL，ETL在T1加权扫描时一般在1至5范围取值，在T2加权扫描时一般在8至24范围取值，重新采集ETL个回波组，并对每组的回波B取复数共轭和时间反演，再沿频率编码方向对每组回波进行离散逆傅立叶变换，得到一系列混合空间矩阵元 和 这里i在1至ETL范围依次取值；
[0010] (3)计算PiA和 之间的
相位差、PiB和 之间的相位差以及PiC和 之间的相位差和 这里i在1至ETL范围依次取值；
[0011] (4)如果α、β和γ有相位缠绕则采用多项式拟合法或区域增长法进行相位解缠；
[0012] 优选的，在场漂效应或具有长时间常数的涡流效应显著情况下，按下述方式采集参考数据并计算相位误差：
[0013] (1)在相位编码梯度关闭条件下按照上述方式采集整个k空间矩阵并重排为DIM1×DIM2×DIM3×DIM4矩阵，其中DIM1为频率编码步数，DIM2是相位编码步数，DIM3是扫描层数，临床上可按通常方式设置，例如，分别为256、192和16，DIM4是化学位移相位编码数，这里DIM4设置为3，再沿选层方向进行一维离散逆傅里叶变换，获得一系列二维k空间矩阵，对其中回波B对应的k空间线进行复数共轭和时间反演，按DIM4＝3分离采集矩阵，再沿频率编码方向进行一维离散傅里叶变换，得到任一层面的复矩阵EA、EB和EC。
[0014] (2)EA的矩阵元沿相位编码方向分为DIM2/ETL组，通过四象限反正切函数计算两两相邻组之间的相位差 这里i取值范围为1至ETL，g取值范围为1至DIM2/ETL-1，这样得到相位矩阵 再将相位矩阵的各列拟合到其中i取值范围为1至ETL。
[0015] (3)同样，对回波B和回波C分别按类似方式得到相位差 和 并分别拟合到和 其中i取值范围为1至ETL。
[0016] 优选的，三维/二维快速多回波序列在三维成像模式下采用多叠交替射频激发方式，基于ky＝0对应的k空间线实时测试中心频率νs及其偏移Δν并按照 自动校准厚
块位置 对应的梯度Gs；根据回波串长度(ETL)、回波时间(TE)和重复时间(TR)选择T1加权或PD加权或T2加权方式；并在正常相位编码条件下采用分段激发方式采集DIM2·X％/ETL个回波组，每组包括A、B和C三个回波对应于回波时间TE、TE-Δt和TE+Δt，X％优选为55％；在
正交接收或多通道接收模式下，回波信号合成方式是是通道i接收的层面j磁共振信号，ai和ΔΦi分别是通道i的灵敏度权重因子和
相移，其标定方式如下：
[0017] (1)对通道i采集的ky＝0对应的矩阵行进行一维傅立叶变换，取模并计算最大值Imax；
[0018] (2)计算矩阵元素I>5％·Imax的横坐标范围[p1,p2]并求绝对值积分面积Ai；
[0019] (3)基于四象限反正切函数计算[p1,p2]范围的相位φi，基于Itoh
算法展开相位φi并计算其平均值，
[0020] (4)对每个通道重复上述过程，计算相移ΔΦi＝<Φi>-<Φi-1>，在计算灵敏度权重因子ai＝Ai/∑Ai。
[0021] 优选的，按下述方式进行相位校正和幅度校正：
[0022] (1)在正常相位编码条件下采集的整个k空间矩阵(DIM1×DIM2×DIM3×DIM4)沿选层方向进行一维离散逆傅里叶变换，获得各扫描层面的二维k空间矩阵，任一层面的k空间矩阵包括DIM2条k空间线，每条k空间线包含DIM1×DIM4个复数点，可按DIM4分离为和
[0023] (2)对KB取复数共轭和时间反演，然后将KA、KB和KC沿频率编码方向进行一维离散逆傅立叶变换后得到一系列混合空间的复数矩阵 和 再沿相位编码方向分为g＝DIM2/ETL组，每组分别乘以 和 其中i从1至ETL依次取值。
[0024] (3)按上述方式校正相位误差后将 和 按照相位编码梯度表重新排序，得到复矩阵 和 然后将 和 沿相位编码方向后进行一维离散逆傅里叶变换，由此获得任一层面的同相图和反相图 和 或者，在半傅立叶采集情况下，将 和 沿频率编码方向后进行一维离散傅里叶变换后按照Cuppen或POCS算法进行半傅立叶重建。
[0025] (4)分别对PiA、PiB和PiC取模，得到|PiA|、|PiB|和|PiC|，并分别对 和 取模，得到 和 这里i从1至ETL取值，计算幅度校正因子矩阵和 这里i从1至ETL
取值，再将 和 沿相位编码方向分为g＝DIM2/ETL组，每组ETL行矩阵元素分别除以ai，bi和ci，其中i从1至ETL取值。
[0026] (5)按上述方式校正幅度不对称性后将 和 按照相位编码梯度表重新排序，得到复矩阵 和 再将 和 沿相位编码方向进行一维离散逆傅里叶变换，由此获得任一层面的同相图和反相图 和 其中 等价于常规T2(或T1或PD)加权图像；或者，在半傅立叶采集情况下，将 和 沿频率编码方向后进行一维离散傅里叶变换后按照Cuppen或POCS算法进行半傅立叶重建。
[0027] 优选的，在场漂效应或具有长时间常数的涡流效应显著情况下，按下述方式进行相位校正：
[0028] (1)在正常相位编码条件下采集的整个k空间矩阵(DIM1×DIM2×DIM3×DIM4)沿选层方向进行一维离散逆傅里叶变换，获得一系列二维k空间矩阵，每个矩阵按照DIM4＝3分离为三个DIM1×DIM2矩阵并沿频率编码方向进行一维离散傅里叶变换得到复矩阵FA、FB和FC；
[0029] (2)将FA、FB和FC沿相位编码方向分为DIM2/ETL组，每组乘以 和进行残余相位误差校正，这里i取值范围为1至ETL；
[0030] (3)按相位编码梯度表重排k空间线，并沿频率编码方向进行一维离散逆傅里叶变换到图像域，对每个层面的k空间矩阵重复上述
数据处理过程，得到图像域复数矩阵和 这里j从1至DIM3取值。
[0031] 优选的，对于Δτ＝1/Δf/2情况，任一层面j的水像和脂肪像分别计算为[0032]
[0033] 同时，基于 获得T2*加权图像，并基于 获得磁场均匀度分布图。
[0034] 优选的，在低场MRI系统上设置Δτ＝1/Δf/3，每组回波的化学位移编码相位分别设置为-2π/3，0，2π/3，三个不同的回波时间满足条件t2-t1＝1/3/Δf和t3-t2＝1/3/Δf，在场强很低(例如0.2T)情况下，Δτ设置可为1/Δf/4，回波峰的相位分别设置为-π/2，0，+π/2，基于下式
[0035]
[0036] 采用
迭代最小二乘法拟合法获得水脂分离完全的水像和脂肪像。
[0037] 有益效果：本发明三维多参数加权同步扫描和校准方法在一次扫描中可以最大限度获得多幅图像，即，同相图、反相图、脂肪像、压脂水像、常规T2加权(或T1加权/PD加权)图像以及T2*加权图像，明显节省了临床扫描时间并增加了临床扫描方案的可选择性，对MRI系统的硬件性能的依赖性更小。
附图说明
[0038] 图1三维快速多回波水脂分离序列。其中，TE为回波时间，相邻回波的峰值之间的时间间隔Δτ设置为1/Δf/2，射频激发脉冲和重聚脉冲选用sinc脉冲或最小相位SLR脉冲，虚框内的部分重复ETL-1次，Gs1和Gs2是选层梯度，其两侧是散相梯度，Gp1和ΔGp1是选层方向相位编码梯度及其递增量，Gp2和ΔGp2是二维平面相位编码梯度及其递增量，虚线框内回波串对应的相位编码方式与常规三维FSE分段激发情况类似，G1是预备读梯度，G2、G3和G4是二维平面频率编码梯度，ΔG2和ΔG4是施加在读梯度方向的补偿梯度，用于校正回波中心位置，-ΔG2和-ΔG4用于对补偿梯度作用下磁化矢量的相位变化进行恢复。
[0039] 图2二维快速多回波水脂分离序列。其中，TE为回波时间，相邻回波的峰值之间的时间间隔Δτ设置为1/Δf/2，射频激发脉冲和重聚脉冲选用sinc脉冲或最小相位SLR脉冲，虚框内的部分重复ETL-1次，Gs1和Gs2是选层梯度，其两侧是散相梯度，Gp1和ΔGp1是选层方向相位编码梯度及其递增量，Gp2和ΔGp2是二维平面相位编码梯度及其递增量，虚线框内回波串对应的相位编码方式与常规三维FSE分段激发情况类似，G1是预备读梯度，G2、G3和G4是二维平面频率编码梯度，ΔG2和ΔG4是施加在读梯度方向的补偿梯度，用于校正回波中心位置，-ΔG2和-ΔG4用于对补偿梯度作用下磁化矢量的相位变化进行恢复。
[0041] 图4信号调试模块工作流程图。
[0042] 图5具有均匀激发特征的射频激发脉冲波形图。
[0043] 图6具有均匀激发特征的射频激发轮廓图。其中，波形特征参数：最小相位SLR，时带积为8，脉宽为2ms，带内和带外纹波系数均为0.5％，激发带宽可通过调节选层梯度幅度改变。
[0044] 图7三维多叠交替激发模式。
[0045] 图8预扫描模块工作流程图。
[0046] 图9扫描模块工作流程图。
[0047] 图10多通道参数标定方案。
具体实施方式
[0048] 为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚，下面结合附图和
实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部内容。
[0049] 医学磁共振成像仪主要由磁体、谱仪、控制台、
梯度线圈、射频线圈、射频功放和梯度功放等硬件单元构成，本发明涉及的一种多参数加权成像方法包括三维/二维快速多回波水脂序列、信号调试、预扫描、扫描、数据预处理、图像重建和水脂分离模块，安装在控制台主机内按照图3所示的扫描流程工作。
[0050] 实施例1：
[0051] 按通常方式在谱仪的序列发生器上加载常规梯度预加重序列并调试涡流场补偿参数进行梯度波形校正，然后在序列发生器上加载图1(或图2)所示的三维/二维快速多回波水脂序列控制各硬件单元实现同相和反相质子信号的激发、
空间编码和采集，这里，射频脉冲波形采用90°脉冲和180°脉冲选用具有均匀激发特征的最小相位SLR脉冲波形，例如图5和图6所示的时带积优先为8，脉宽优选为2ms，带内和带外纹波系数均为0.5％的最小相位SLR波形，Gs是厚块选择梯度，Gp1是选层方向的相位编码梯度，Gp2是二维平面相位编码梯度，G1是预备读梯度，G2、G3和G4是二维平面频率编码梯度，ΔG2和ΔG4是补偿梯度，虚线框表示可重复执行的序列部分。设置第一个频率编码梯度的面积为预备读梯度面积的两倍，设置每组频率编码梯度(G2、G3和G4)的积分面积初始值为G1梯度的两倍，其中G2和G4的极性与G1的极性相反。参照常规快速自旋回波序列的分段激发方式设置回波链长度ETL和相应的相位编码梯度值，在每组频率编码梯度(G2、G3和G4)施加期间采集一组回波(包括梯度回波A、自旋回波B和梯度回波C)，然后保存序列文件和相位编码梯度表。在序列参数表中，设置采集矩阵大小为DIM1×DIM2×DIM3×DIM4＝256×192×16×3，其中DIM1、DIM2、DIM3和DIM4分别表示成像序列中的频率编码步数、相位编码编码步数、扫描层数和化学位移相位编码数，采用多叠交替三维采集模式，选层方向扫描
视野切分为多个(例如9个)激发厚块并按照图7所示的编号次序进行交替激发和半傅立叶采集，每个厚块的厚度为16mm，近邻厚块之间的
叠加范围为2mm，相当于边缘2层重叠。其它参数设置如下：
[0052] 对于T1加权和水脂分离同步扫描，设置ETL为1～5范围的一个恰当值，设置回波时间TE＝10ms，设置序列重复时间TR＝450ms，对于T2加权和水脂分离同步扫描，设置ETL为8～24范围的一个恰当值，TE＝12ms，设置TR为2000～4000ms范围的一个恰当值，对于PD加权和水脂分离同步扫描，设置ETL为1～5范围的一个恰当值，TE＝12ms，设置TR为2000～4000ms范围的一个恰当值。另外，设置视野FOV＝250mm，块厚THK＝32mm，累加次数NEX＝1，按常规方式设置其它参数并调试90°脉冲射频功率。该优化方案在临床环境下不仅适合高场扫描也适合低场扫描，例如，对于T2加权和水脂分离同步扫描，扫描时间TACQ＝TR·DIM2/ETL·DIM3/60＝2.5·108/16·16/60＝4.5(min)。然后，保存序列文件和参数表文件并按下述方式对每个厚块进行信号采集和处理：
[0053] 首先，按照图4所示工作流程在相位编码梯度关闭条件下自动执行信号校准模块100，设置G2＝0，G4＝0，ΔG2＝0，ΔG4＝0，G3＝2·G1，微调G1使得回波B的顶点在采样窗中心位置；设置G2＝2·G1，G4＝2·G1，反转G3的极性，用图谱处理常用的标尺标记回波幅值顶点，测量顶点之间的时标间隔Δτ，分别调试补偿梯度ΔG2和ΔG4直到回波A的顶点与回波B的顶点的时标间隔Δτ＝1/Δf/2，以及回波B的顶点与回波C的顶点的时标间隔Δτ＝1/Δf/2，这里水和脂肪质子的化学位移差Δf在序列内设定，例如1.5T场强下210Hz，3.0T场强下
420Hz。
[0054] 其次，按照图8所示的工作流程执行预扫描模块110，测试场漂并限制场漂引起的中心频率
波动范围不超过扫描层厚对应的
频率范围；然后，按照图9所示工作流程执行扫描模块120，在序列循环之间基于ky＝0对应的k空间线测试中心频率νs及其偏移Δν并按照自动校准厚块位置 对应的梯度Gs。
[0055] 其三，在采样期间双极性频率编码梯度切换引起的涡流场或麦克斯韦场不一致，难以通过180°脉冲和梯度校准彻底消除，需要通过数据预处理模块130进行低阶和高阶相位校正，其实现方式说明如下：
[0056] (1)在所有相位编码梯度关闭条件下，设置DIM2＝ETL，采集ETL个回波组，每组包括A、B和C三个回波，其回波峰相位分别为-π、0和+π，并对每组的回波B取复数共轭和时间反A演，再沿频率编码方向对每组回波进行离散逆傅立叶变换，得到一系列混合空间矩阵元Pi、PiB和PiC，这里i表示分段激发的各回波编号，取值范围为1至ETL；
[0057] (2)在所有相位编码梯度关闭且频率编码梯度极性反转条件下，设置DIM2＝ETL，重新采集ETL个回波组，并对每组的回波B取复数共轭和时间反演，再沿频率编码方向对每组回波进行离散逆傅立叶变换，得到一系列混合空间矩阵元 和 这里i在1至ETL范围依次取值；
[0058] (3)计算PiA和 之间的相位差、PiB和 之间的相位差以及PiC和 之间的相位差 和 这里i在1至ETL范围依次取值，如果α、β和γ有相位缠绕则采用多项式拟合法或区域增长法进行相位解缠；
[0059] (4)在正常相位编码条件下，采集整个k空间矩阵并沿选层方向进行一维离散逆傅里叶变换，获得各扫描层面的二维k空间矩阵，任一层面的k空间矩阵包括DIM2条k空间线，每条k空间线包含DIM1×DIM4个复数点，可按DIM4分离为和
[0060] (5)对KB取复数共轭和时间反演，然后将KA、KB和KC沿频率编码方向进行一维离散逆傅立叶变换后得到一系列混合空间的复数矩阵 和 再沿相位编码方向分为g＝DIM2/ETL组，每组分别乘以 和 其中i从1至ETL依次取值。
[0061] (6)按上述方式校正相位误差后将 和 按照相位编码梯度表排序，得到复矩阵 和
[0062] (7)然后将 和 沿相位编码方向后进行一维离散逆傅里叶变换，由此获得任一层面的同相图和反相图 和 或者，在半傅立叶采集情况下，将和 沿频率编码方向后进行一维离散傅里叶变换后按照Cuppen或POCS算法进行半傅立叶重建。
[0063] 其三，在双极性梯度频率编码条件下，接收通道
滤波器幅度响应的不对称性会引起回波幅度不对称以至于水脂分离不完全，所以在上述相位误差校正基础上 和需要进行幅度校正如下：
[0064] (1)分别对PiA、PiB和PiC取模，得到|PiA|、|PiB|和|PiC|，这里i从1至ETL取值；
[0065] (2)分别对 和 取模，得到 和 这里i从1至ETL取值；
[0066] (3)计算幅度校正因子矩阵 和这里i从1至ETL取值；
[0067] (4)将 和 沿相位编码方向分为g＝DIM2/ETL组，每行矩阵的ETL个元素分别除以ai，bi和ci，其中i从1至ETL取值；
[0068] (5)按上述方式校正幅度不对称性后将 和 按照相位编码梯度表排序，得到复矩阵 和
[0069] (6)将 和 沿相位编码方向进行一维离散逆傅里叶变换，由此获得任一层面的同相图和反相图 和 其中 等价于常规T2(或T1或PD)加权图像；或者，在半傅立叶采集情况下，将 和 沿频率编码方向后进行一维离散傅里叶变换后按照Cuppen或POCS算法进行半傅立叶重建。
[0070] 其四，在场漂效应或具有长时间常数的涡流效应显著情况下，序列重复周期之间可能存在额外的相位误差，为此，按下述步骤进行相位误差消除：
[0071] (1)在相位编码梯度关闭条件下按照上述方式采集整个k空间矩阵并重排为DIM1×DIM2×DIM3×DIM4矩阵，再沿选层方向进行一维离散逆傅里叶变换，获得一系列二维k空间矩阵，对其中回波B对应的k空间线进行复数共轭和时间反演，按化学位移相位编码数DIM4＝3分离采集矩阵，再沿频率编码方向进行一维离散傅里叶变换，得到任一层面的复矩阵EA、EB和EC；
[0072] (2)EA的矩阵元沿相位编码方向分为DIM2/ETL组，通过四象限反正切函数计算各组与第一组的相位差 这里i取值范围为1至ETL，g取值范围为1至DIM2/ETL-1，这样得到相位矩阵
[0073]
[0074] (3)将相位矩阵的各列分别拟合到 其中i取值范围为1至ETL；
[0075] (4)同样，对回波B和回波C分别按类似方式得到相位差 和 并分别拟合到和 其中i取值范围为1至ETL；
[0076] (5)在正常相位编码条件下按照上述方式采集整个k空间矩阵(DIM1×DIM2×DIM3×DIM4)并沿选层方向进行一维离散逆傅里叶变换，获得一系列二维k空间矩阵，每个按照DIM4＝3分离为三个DIM1×DIM2矩阵并沿频率编码方向进行一维离散傅里叶变换得到复矩阵FA、FB和FC；
[0077] (6)将FA、FB和FC沿相位编码方向分为DIM2/ETL组，每组乘以 和进行残余相位误差校正，这里i取值范围为1至ETL；
[0078] (7)按相位编码梯度表重排k空间线，并沿频率编码方向进行一维离散逆傅里叶变换到图像域，对每个层面的k空间矩阵重复上述数据处理过程，得到图像域复数矩阵和 这里j从1至DIM3取值。
[0079] 在多通道接收情况下，按下述方式进行信号合成：
[0080] (1)对通道i采集的ky＝0对应的矩阵行进行一维傅立叶变换，取模并计算最大值Imax；
[0081] (2)计算矩阵元素I>5％·Imax的横坐标范围[p1,p2]并求绝对值积分面积Ai；
[0082] (3)基于四象限反正切函数计算[p1,p2]范围的相位φi，基于Itoh算法展开相位φi并计算其平均值，
[0083] (4)对每个通道重复上述过程，计算相移ΔΦi＝<Φi>-<Φi-1>，在计算灵敏度权重因子ai＝Ai/∑Ai。
[0084] 然后，对于Δτ设置为1/Δf/2的情况，按下述方式对 进行场不均匀性校正并进行水脂分离：
[0085]
[0086]
[0087]
[0088] 基于式(2)计算初始相位如下：
[0089]
[0090] 基于常用的区域增长法进行相位解缠，再对同相和反相图像消除φ0如下：
[0091]
[0092]
[0093]
[0094]
[0095]
[0096] 这里定义 用于确定Sw-Sf的符号。
[0097] 然后，按下述方式进行水脂分离，分别得到水像Sw和脂肪像Sf：
[0098]
[0099]
[0100] 其中
[0101] 最后，按照前述厚块交替激发次序选择近邻厚块，将其中叠加区域的图像
像素进行加和平均，并采用中值滤波或非局域均值滤波进行图像平滑去噪，得到完整的三维水像和脂肪像。
[0102] 上述成像方法最大程度上降低了硬件
缺陷或物理效应的影响，可方便地同时获得水脂分离图像和常规T2加权图像，在DIM2＝ETL且ETL<4情况下可获得T1加权或PD加权图像，也可同时基于下式
[0103]
[0104] 获得T2*加权图像，并基于下式
[0105]
[0106] 获得场图，单位为赫兹。在某些特殊情况下个别像素的水脂质子含量接近相等，这可导致式(13)的计算产生异常点，这些异常点可通过近邻插值或双线性插值方式弥补。
[0107] 实施例2：
[0108] 对于Δτ设置为1/Δf/3的情况，在上述扫描方案中将每组回波的化学位移编码相位分别设置为-2π/3，0，2π/3并采集数据，按类似方式校正相位和幅度误差并进行图像重建后进行数据处理如下：
[0109] 对于任一层面j所对应的图像的任一像素，其对应的磁共振信号可按下式描述：
[0110]
[0111] 这里，Sw和Sf分别表示人体组织中水和脂肪的质子磁化矢量的初始值，后者分别构成水像和脂肪像的复数矩阵，下标w和f分别表示水和脂肪，ν表示场不均匀性或涡流效应引起的
频率偏差，tn(n＝1，2或3)表示三个不同的回波时间，满足条件t2-t1＝1/3/Δf和t3-t2＝1/3/Δf。在正交接收或多通道接收模式下， 是通道i接收的层面j磁共振信号，ai和ΔΦi分别是通道i的灵敏度权重因子和相移，通过图10所示的标定方案测定。
[0112] 已知频率偏差ν是一个可测量，式(14)在消除其效应后改写为
[0113]
[0114] 上式中上标R和I分别表示复数的
实部和
虚部。这里采用线性最小二乘法拟合求解上述线性方程组得到水脂分离图像和场图。
[0115] 上述不对称回波采集和数据处理方式的额外优势是可避免水和脂肪含量相当的极个别像素出现异常的计算误差，但需要另行采集常规T2加权(或T1加权或PD加权)图像。在场强很低(例如0.2T)情况下，Δτ设置可为1/Δf/4，回波峰的相位分别设置为-π/2，0，+π/2，每次分段激发并采集三幅图像，然后按照上述相位和幅度校正方式和最小二乘法拟合迭代方式得到水脂分离图像和场图。
[0116] 最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其对前述各实施例所记载的技术方案进行
修改，或者对其中部分或者全部技术特征进行等同替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
