美文网首页
12.4 非均匀采样

12.4 非均匀采样

作者: 十年磨剑_莫回首 | 来源:发表于2021-05-16 18:33 被阅读0次

沿着k-space数据的任何一个方向,都存在这样的因素,涡流,模数转换误差,时间调换误差,这些都会导致数据的非均匀间隔。这种类型的误差会沿着读取方向,相位方向和分区编码线出现,这种非均匀采样会导致重建图像的混叠效应。

在k-space中间如果某些数据线的采样时非均匀的,我们通常会把这种非均匀采样获得的数据表示为多个均匀互相有位移的采样数据的叠加。为了探究混叠的存在性,傅里叶变换经常会沿着读取线或者相位编码线进行。

12.4.1 来自于交叉采样的混叠

我们对信号进行采样,采样间隔为\Delta k = 1/L, 分别对2m \Delta k的数据点和(2m+1)\Delta k的数据点分别采样,获得两组采样数据,一组是偶数位置数据集s_{e}(k),一组是奇数位置数据集s_{o}(k)。在每组采样数据中,其空间间隔为 2  \Delta k。这样的两组数据,拿出任意一组,去进行重建,图像都会有混叠现象,如果我们把这两组数据结合在一起,那么每个数据的空间间隔就是\Delta k,从而减少混叠的发生。

数据正确交叉

采样函数u(k)可以被重新写作两个采样函数的和:一个是奇采样函数u_{o}(k),另外一个是偶采样函数u_{e}(k)。有如下式子:

u(k)=u_{o}(k)+u_{e}(k)

对于奇彩样函数和偶采样函数的定义:

u_{e}(k)=\Delta k \sum_{m=-inf}^{inf} \delta (k-2m\Delta k)     \\         

u_{o}(k)=\Delta k \sum_{m=-inf}^{inf} \delta (k-(2m+1)\Delta k)

无限采样的数据的重建图像是实际的spin density \r  (x)和采样函数U(x)的卷积。我们也可以把这样的采样函数写成:

U(x)=U_{e}(x)+U_{o}(x)

然后我们对这样的采样函数进行进行反傅里叶变换,可以得到:

然后我们分别使用这两个采样函数对\r  (x)进行采样,可以得到如下两个结果:

上面得到的\hat  \r  _{o}(x)\hat \r  _{e}(x)其实是\hat \r _{inf}(x)的两个子成分。

由于奇数项在两个序列里面被消减了,最后得到的就是无穷采样数据的的重建结果:

我们可以假设在-L/2 \quad to \quad L/2 之前的区域,\r  (x)为0 ,\r _{inf}(x)=\r (x), 那么就不存在混叠的影响。但是在奇序列和偶序列中如果q取-1,0,1 这三个值的话,那么就机会存在混叠。我们可以这样定义\r  (x)=\r _{l}(x) + \r  _{r}(x) 我们这样定义:

这样定义的左右函数对于描述混叠特别有帮助。

如果我们把这两个重建结果加在一起,那么混叠的哪一步成分就会消除掉,而只保留\r  (x)

数据不正确交叉

数据不正确交叉,就是说数据之间的间隔不是统一的\Delta k
,我们还是以奇采样函数和偶采样函数为例。

如果偶采样的间隔不变,但是奇采样变为:

这里-1<\alpha<1, 当其取值为0的时候,就是一个完美的奇序列采样。这样的采样函数在空间域就会产生一个相移:

然后我们遵照以前定义的左右函数,可以定义由奇序列重建得到的图像:

偶序列重见得到图像依旧是

最后把奇偶序列得到的重建图像相加,得到:

这里可以看到,我们如果对\hat  \r   (x)取幅值,依旧无法去除相移带来的影响,也就无法去除混叠效应。

当然了 这是最简单的一种情况,还有其他更多的交叉采样,每个序列的采样间隔不同,都会带来无法消除的相移,这一点在图像最后重建时候,可以看出来。

12.4.2 在梯度上面的模数转换带来的混叠

对于相位编码方向的梯度间隔, \Delta G_{y} \equiv \Delta G, 它要足够小,这样才能获得一个比较大的FOV。 但是这样会带来一个约数取整的问题: 我们取一个12.5 \Delta G_{DAC,min}作为梯度间隔,

其中\Delta G_{DAC,min} 是MR 梯度放大器能够产生的最小的梯度间隔。,然后我们沿着相位编码方向的梯度场就是12.5,25,37.5,50,... 省略了梯度场基础单位。但是在数模转换器里面,由于只能传递整数编码,这个序列会变成13,25,38,50,...,这就会导致在k-space空间的采样间隔不统一,进而产生无法通过取幅值去除的相移干扰和幅值干扰,从而带来混叠的影响。

相关文章

  • 12.4 非均匀采样

    沿着k-space数据的任何一个方向,都存在这样的因素,涡流,模数转换误差,时间调换误差,这些都会导致数据的非均匀...

  • 圆盘均匀采样

    考虑在半径为1的圆盘上均匀的取随机点应该怎么做? Rejection Sampling 首先想象有一个2X2的单位...

  • PointRend实现细节

    如何选择采样点? 从均匀分布随机采样kN个点 重点采样BN个点 从均匀分布中采样(1-B)N个点 def get_...

  • 论文阅读笔记:Frequency Bias in Neural

    两层非均匀分布神经网络NTK的特征函数 我们开始研究从非均匀分布中采样的两层神经网络的NTK模型。首先关注1D目标...

  • 计算机实现非均匀采样的方式。

    拒绝采样https://blog.csdn.net/baimafujinji/article/details/51...

  • 点云采样

    原文链接 点云采样分类 点云采样的方法有很多种,常见的有均匀采样,几何采样,随机采样,格点采样等。下面介绍一些常见...

  • 【SQL】抽样

    随机采样 分层采样 hash 版 非hash 版

  • 不同形状内的均匀采样

    1. 两椭圆间采样[http://jsfiddle.net/SQERh/] 2. 椭圆绘制 3. 椭圆均匀采样 3...

  • 概率简要学习记录

    随机数问题 构造均匀的随机数发生器 要等概率才可以丢掉 不均匀的随机数产生器 采样问题 水库采样利用数组和随机数取...

  • 2019-03-28

    对于简单函数的分布,如均匀分布,高斯分布等,很多编程语言有直接的采样函数; 复杂的分布,往往没有直接的采样函数可供...

网友评论

      本文标题:12.4 非均匀采样

      本文链接:https://www.haomeiwen.com/subject/okcrjltx.html