专利名称：一种用于磁共振成像基于相位一致性的图像去噪方法
技术领域：
本发明属于磁共振成像领域，具体涉及ー种用于磁共振成像基于相位一致性的图像去噪方法。
背景技术：
信噪比是磁共振图像中最重要的參数之一，然而在磁共振成像中由于电子环境以及电子接收链路不可避免的引入一定的背景噪声，另ー方面供磁共振成像的人体核磁信号是ー种比较微弱的电磁信号，所以磁共振信号很容易受到噪声的干扰，提高图像信噪比ー直人们不断追求的目标。为了提高磁共振成像的信噪比，一方面要増加接收到的信号強度，通过优化成像序列，提高主磁场以及射频场的均匀度，应用灵敏度更高的多通道接收线圈等方法都可以 改善接收到的信号強度；另ー方面在电子电路的硬件设备上尽可能的减小引入的背景噪声。然而在经过上述方法的优化之后，磁共振图像的信噪比往往还存在进ー步提高的空间，特别是ー些低信号区域容易受到噪声的影响，一些小的信号被本底噪声掩盖，不能被有效检出，降低了图像的质量。在序列优化和硬件条件都达到一定水平后，磁共振成像的信噪比还可以通过图像后处理去噪的方法进ー步得到提高。传统的图像去噪算法有高斯滤波，小波滤波等方法，然而这些方法在滤波的同时会损失ー些图像细节，而且在信号与噪声处于同等量级的情况下，这些传统方法很难把有效信号从噪声中甄别出来。更重要的是，这些方法都是一般的图像处理算法，并没有利用磁共振图像中特殊的信息，比如传统方法一般只处理幅度信号而没有用到磁共振图像中的相位信息；传统方法都是处理単一的一幅图像而没有用到磁共振成像的多通道图像信息。
发明内容
本发明目的是提供ー种利用磁共振图像中的相位信息以及多通道图像信息、成像效果好的用于磁共振成像基于相位一致性的图像去噪方法。本发明的技术方案是ー种用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，利用图像相位的标准差与信噪比关系曲线计算出像素点信噪比并计算得到模值放大修正系数，根据模值放大修正系数对多通道幅度图像进行模值修正。优选的，具体包括以下步骤(I)通过均值为零的相位标准差与信噪比的关系曲线计算出每个像素点的信噪比；(2)通过相位图像的幅值计算每个像素点的信噪比上限；(3)通过幅度图像的幅值结合各个通道的本底纯高斯噪声的标准差，计算每个像素点的信噪比的上限和下限；(4)根据上述步骤(I)、(2)、(3)的结果计算出每个像素点的信噪比值；
(5)根据信号模值均值放大系数随信噪比的变化曲线得到每个像素点对应的模值放大修正系数，依此系数对所有像素点的模值进行修正；(6)将所有通道图像合成，得到最終的去噪图像。进ー步的，所述步骤(3)具体为关闭射频发射，采集各个通道的纯噪声，计算出各个通道的噪声标准差，根据得到各个通道的噪声标准差结合多通道幅度图像计算出每个像素点的信噪比的上限和下限。其中各个通道的噪声标准差通过采集各个通道的纯噪声数据统计分析得到。采集纯噪声时只打开采集窗，不发送射频脉冲和梯度。进ー步的，所述步骤(I)具体为利用多通道合成的复数图像的相位剔除所有通道图像中由于成像序列导致的相位累积，再通过傅里叶级数拟合各个通道线圈灵敏度的相位分布并从图像相位中去除，得到相位真值为零的多通道相位图，再通过窗ロ滑动加权平均的方法计算该图像每个像素点的相位标准差，再通过均值为零的相位标准差与信噪比的关系曲线计算出每个像素点的信噪比。 进ー步的，所述步骤(6)图像合成的方法为平方和开根号的方法或自适应多通道合成的方法。本发明的优点是充分利用磁共振成像多通道的相位信息，通过相位的一致性分析计算像素点的信噪比，依次对低信噪比的像素加以较大的幅度值修正。本发明可以有效地将本底噪声抑制到最低限度，与此同时保留了有效信号，将原本被噪声掩盖的小信号恢复出来，不但不会损失图像细节，而且提取出更多有用信号，大大提高了图像的质量。
下面结合附图及实施例对本发明作进ー步描述图I为本发明的去噪方法流程图。图2是去噪前的图像。图3是用本发明方法去噪后的图像。
具体实施例方式实施例本发明的用于磁共振成像基于相位一致性的图像去噪方法，如图I所示，包括以下步骤(I)磁共振序列采集多通道复数图像得到多通道幅度图像；(2)关闭射频发射，采集各个通道的纯噪声，计算出各个通道的噪声标准差，根据得到各个通道的噪声标准差结合多通道幅度图像计算出每个像素点的信噪比的上限和下限；(3)利用多通道合成的复数图像的相位剔除所有通道图像中由于成像序列导致的相位累积，再通过傅里叶级数拟合各个通道线圈灵敏度的相位分布并从图像相位中去除，得到相位真值为零的多通道磁共振图像；(4)用窗ロ滑动加权平均的方法计算该图像每个像素点的相位标准差，再通过均值为零的相位标准差与信噪比的关系曲线计算出每个像素点的信噪比；(5)通过相位图像的幅值计算每个像素点的信噪比上限；(6)综合上述步骤(3)、(4)、(5)的結果，确定每个像素点的信噪比计算值；
(7)利用信号模值均值放大系数随信噪比的变化曲线得到每个像素点对应的模值放大修正系数，依此系数对所有像素点的模值进行修正，抑制噪声，保留信号；(8)将所有通道图像合成，得到最終的去噪图像。本发明的去噪方法通过对于磁共振多通道图像相位信息一致性的分析，计算像素点的信噪比，对信噪比低的像素点加上较大的幅度修正，在有效抑制噪声的同时保留信号，特别是原本已被噪声掩盖的小信号，显著提高图像质量。具体算法如下(I)任意第j通道的复数图像ふ=Pe^Ci，其中p为真实图像，p为序列累积的相位，Cj为第j通道的复数灵敏度分布。合成图像S =。首先利用多通道合成的复数图像的相位剔除所有通道图像中由于成像序列导致的相位累积,ろ矿=再通过傅里叶级数拟合各个通道线圈灵敏度的相位分布，即 吸/< 并从图像相位中去除，得到相位真值为零的多通道磁共振图像相位舞がイも)= OTgZe(Af)-OTgte(CyK)11 (2)用窗ロ滑动加权平均的方法计算该图像每个像素点的相位标准差，再通
过均值为零的相位标准差与信噪比的关系曲线计算出每个像素点的信噪比。相位标准差与信噪比的关系曲线通过以下公式计算
ap=Jjf —)—f e—制—^tan2 f\d0d(^-)
9 Ui{sm+cos0j UJ(3)通过相位图像的幅值计算每个像素点的信噪比上限max I SNR}=——-,-—,
v sin(a g/e(ろ))(4)通过幅度图像的幅值结合各个通道的本底纯高斯噪声的标准差O，计算每个像素点的信噪比的上限和下限
r n|,VJ + 3cr ,、，でmaxl|,S'..|-3<j.o)max、SNR、m = -，以及 min {SNR、m =-;~--(5)综合上述步骤(2)、(3)、(4)的結果，确定每个像素点的信噪比计算值
min {SNR) SNR (cr,) < min ISNR!び;)min {SNR} < SNR (c-) < max {SNR ；,
max {SNR] SNR (び)> max 丨 SNK；其中min {SNR} = min {SNR}m, max {SNR} = min (max {SNR}m, max (SNRjp)(6)利用信号模值均值放大系数随信噪比的变化曲线得到每个像素点对应的模值放大修正系数，依此系数对所有像素点的模值进行修正，抑制噪声，保留信号。信号模值均值放大系数随信噪比的变化曲线通过以下公式计算
权利要求
1.ー种用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，利用图像相位的标准差与信噪比关系曲线计算出像素点信噪比并计算得到模值放大修正系数，根据模值放大修正系数对多通道幅度图像进行模值修正。
2.根据权利要求I所述的用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，具体包括以下步骤 (1)通过均值为零的相位标准差与信噪比的关系曲线计算出每个像素点的信噪比； (2)通过相位图像的幅值计算每个像素点的信噪比上限； (3)通过幅度图像的幅值结合各个通道的本底纯高斯噪声的标准差，计算每个像素点的信噪比的上限和下限； (4)根据上述步骤(I)、(2)、(3)的结果计算出每个像素点的信噪比值；(5)根据信号模值均值放大系数随信噪比的变化曲线得到每个像素点对应的模值放大修正系数，依此系数对所有像素点的模值进行修正； (6)将所有通道图像合成，得到最終的去噪图像。
3.根据权利要求2所述的用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，所述步骤(3)具体为关闭射频发射，采集各个通道的纯噪声，计算出各个通道的噪声标准差，根据得到各个通道的噪声标准差结合多通道幅度图像计算出每个像素点的信噪比的上限和下限。
4.根据权利要求2所述的用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，所述步骤(I)具体为利用多通道合成的复数图像的相位剔除所有通道图像中由于成像序列导致的相位累积，再通过傅里叶级数拟合各个通道线圈灵敏度的相位分布并从图像相位中去除，得到相位真值为零的多通道相位图，再通过窗ロ滑动加权平均的方法计算该图像每个像素点的相位标准差，再通过均值为零的相位标准差与信噪比的关系曲线计算出姆个像素点的信噪比。
5.根据权利要求2所述的用于磁共振成像基于相位一致性的图像去噪方法，其特征在 于，所述步骤(6)图像合成的方法为平方和开根号的方法或自适应多通道合成的方法。
全文摘要
本发明公开了一种用于磁共振成像基于相位一致性的图像去噪方法，其特征在于，利用图像相位的标准差与信噪比关系曲线计算出像素点信噪比并计算得到模值放大修正系数，根据模值放大修正系数对多通道幅度图像进行模值修正。本发明的去噪方法通过对于磁共振多通道图像相位信息一致性的分析，计算像素点的信噪比，对信噪比低的像素点加上较大的幅度修正，在有效抑制噪声的同时保留信号，特别是原本已被噪声掩盖的小信号，显著提高图像质量。
文档编号G06T5/00GK102800057SQ20121022966
公开日2012年11月28日 申请日期2012年9月18日 优先权日2012年9月18日
发明者姜忠德, 唐昕, 李鹏宇, 陈铭明 申请人:苏州安科医疗系统有限公司