基于二维傅里叶变化法的MRI成像原理的Matlab仿真（2）

二、MRI仿真数据的获取

        基于（1）中的分析，我们知道K空间和图像空间是一对傅里叶变换的关系。因此，仿真数据获取可以通过对原始图像做2-D FT来实现。而为了验证仿真数据获取的正确性，可以对仿真数据做2-D IFT，得到重建后的图像，比较原始图像和重建图像，从而验证重建算法的正确性。

        在Matlab中，Shepp-Logan模型和人脑结构相似，是医学影像学中最先用于CT图像重建仿真领域的经典模型[1]。我们以Shepp-Logan模型为基础，仿真MRI成像的实现，模拟得到K空间的原始数据。模型图像的像素点可以选择为128*128、256*256或512*512，像素点决定了图像的理论分辨率。每个椭圆对应的位置和灰度值不同，分别模拟人的大脑的不同结构。为了更完整地模拟人的头部结构，体现头皮和脑脊液，我们在原有的10个椭圆的基础上，增加了两个椭圆（标号为1*和2*），并对各椭圆的灰度值进行了改进，以提高图像的对比度。创建的256*256的模型图像见图1.

图1 改进的Shepp-Logan模型图像

        用2-D FT重建后的图像如图2所示。