MRI图像噪声的最为主要的来源是检测对象和电路元件,产生的噪声主要为热噪声,有时候还存在生理学噪声。磁共振设备采集到的原始数据是K空间数据,其包含实部和虚部这两部分的信号,且两者之间的相位差为90°。实部和虚部分别带有均值为0,方差相同且独立的加性高斯白噪声。对K空间数据进行傅里叶变换,得到图像域数据。傅里叶变换不改变噪声形式,因此我们得到的图像域数据是复数信号,其表达式为:
其中和
分别为实部和虚部信号,
和
分别服从均值为0,方差为
的高斯分布。通常情况下我们观测的图像是经过模运算得到的幅值图像,即
从该式子可以看出来,最后得到的复制图像已经不再是简单的加性高斯噪声了,而是与信号相关的莱斯噪声(Rician)。在低信号区域,Rician噪声近似于瑞利分布,在高信号区域,Rician噪声近似于高斯分布。随着多线圈成像技术的快速发展,MR成像速度和质量得到了提高,但是由此得到的图像噪声不再服从Rician分布了,而是更为复杂的非偏心卡方分布。其产生的图像表示为
其中C为采集通道数。
只有在对K空间数据进行全采样,并且假设各个线圈得到的数据是不相关的前提下,MRI图像才服从Rician或是非中心卡方分布。
关于如何添加莱斯噪声,请浏览下面网址:
这里面的ricernd函数下载地址: