学习MRI成像的螺旋桨填充

论文《Motion Correction With PROPELLER MRI: Application to Head Motion and Free-Breathing Cardiac Imaging》的个人理解，存在错误请指正：
一、MRI成像过程

1. 频率、相位的编码
2. 填充K空间
3. 傅里叶反变换

二、工作内容
　　作者提出了一种K空间的填充方式，利用多条围绕K空间原点旋转的同心矩形条中收集数据，即K空间的螺旋桨式填充，减少了因运动产生的伪影。伪影产生的原因主要是运动产生的，具体可分为：
1、I型 ：数据采样周期和随后的射频激励之间的静止周期内组织位移；
2、II型：激励脉冲和随后的数据采样周期之间磁场梯度运动而产生的自旋相位。

三、何为K空间
　　K空间即二维空间的频率矩阵，将一组连续空间频率的波进行二维傅里叶变换，从图像角度可以理解成图像的二维傅里叶谱。在图像的二维傅里叶谱中，中心部分代表了图像的内容即存在对比度的部分，而边缘部分则是图像的细节部分即没有对比度的部分。下图为一幅二维脑部图像的K空间采样过程与反变换的过程，从该过程中可以看出：

• 当采样圆的直径的越来越大时，右侧图像也逐渐清晰
• 当圆的直径大到一定程度时，即达到采样边缘时，图像的变化肉眼几乎不可见

注：动图源于https://zhuanlan.zhihu.com/p/139047870

四、数据采样
　　常见的填充方式有条状的水平或垂直方向上的步进式填充，螺旋式填充以及螺旋桨式填充。上述动图中的只对中心部分采样一次，由于K空间中心部分对应于图像的内容，那么每旋转采样一次则会对中心采样那么对多次采样，对多次采样的结果做平均则会得到锐化的效果，减少了内容部分的伪影。此外，K空间边缘没有对比度，肉眼很难观察。螺旋桨式的采样打散了K空间，难以形成完整的伪影。

五、一些相关矫正
1、相位矫正
K空间的中心位移导致每个条带的线性相位变化，去除每个条带在图像空间中所有低频空间变化的相位，同时去除K空间位移，以及去除在图像变化缓慢且可能随条带变化的运动的相关相位。

2、旋转矫正
　　去除每个条带之间物体的整体旋转，图像的空间域旋转与傅里叶谱图旋转的角度一致。图像旋转时，旋转中心位移较小，边缘的旋转位移较大。将每个条带在极坐标上进行网格化并加权平均，然后每一个都旋转一个设定范围的角度，然后为每个条带测量相关性。
3、平移矫正
　　将每个条带得到的低分辨率的图像进行配准求均值。当低分辨率的卷积取最大时相关性最大，对卷积进行傅里叶变换，找到峰值大小，并在x和y方向上对峰值的大小进行三点抛物线拟合。