MRI与其它影像学根本区别在于使用者能够控制采集数据及图像重建的方式与方法,即通过改变一些软件控制、脉冲时间、数据采集顺序、辅助磁场的强度与变化率,就可以改变对比度、分辨力、采集速度、视野(FOV)、伪影效应等等[1],这些控制的核心是K-空间(K-space)。

K-空间与付理叶变换
  K-空间是一个抽象空间(三维空间)或平面(二维空间),MR成像数据根据不同的空间频率编排在特定的K-空间位置,最后被变换成图像。由于K-空间以空间频率为单位(Hz/cm),空间频率K又是由空间互垂的3个分量Kx、Ky、Kz来描述,Kx、Ky、Kz正好对应一个三维频率空间,所以将该抽象空间称为K-空间。所谓空间频率是指在一定方向上的单位空间(距离)波动的周期数,它不仅具有大小而且具有方向,是一个矢量。如一列波在Kx方向频率为3Hz/cm,在Ky方向上频率为4Hz/cm,则该波在Kx-Ky二维空间的空间频率为(42+32)1/2=5Hz/cm,其相位(即波动传播方向)由Kx与Ky值决定,为arctg4/3。空间频率用来描述某些以波的形式在空间传播的能量,与时间频率(Hz/s)不同。
  K-空间与付理叶变换有密切关系,所以认识K-空间首先应认识付理叶变换。付理叶变换是将强度-时间的对应关系,转变为强度-频率的对应关系。某种含有多种频率的信号,如复色光、射频、MR信号等,强度随时间的变化关系可转变为强度-频率的相关关系(图1),转换后的每条线位置对应于其频率,高度对应相应频率信号的幅度,这种将混合频率信号的强度与频率的对应关系在一维数组上的表达方法称一维付理叶转换[2]。

图1 一维付理叶变换。时间-强度关系与频率-强度关系变换示意图。

  有些信号的频率,具有矢量性,如空间频率,一维付理叶变换不能完全表达其物理意义,必须将信号与频率的对应关系在二维空间上才能完全表达出来,这种方法称二维付理叶变换。如将一组连续空间频率的波进行二维付理叶变换,结果是产生一个二维空间频率矩阵,即二维K-空间,二维K-空间中每一点对应一个具有一定空间频率的信号,信号强度SI(Kx,Ky)决定于该点的亮度(图2)。

图2 二维付理叶变换对应产生二维K-空间示意图。

K-空间内每一点对应Kx、Ky及信号强度,越近K-空间中心部,频率越低,反之则越高。

  更加复杂的信号必须在三维空间才能表达其信号,这种方法称三维付理叶变换。三维付理叶变换产生一个三维空间频率矩阵,即三维K-空间,矩阵中每一点数据决定于MR信号强度。
  K-空间是一个无法观察到的空间,但我们可以通过适当的方法认识它,了解它在MR中的作用,并充分利用对K-空间的控制技术优化MR图像质量。

K-空间与图像的关系
  为了进一步认识K-空间,首先认识一下简单的光学成像(图3),在第一个透镜前一倍焦距处的物体发射的光线经第一透镜折射后,在其后一倍焦距平面处放一纸即可观察到规律分布的所谓干涉图形,再于第一透镜后二倍焦距处放置第二透镜,则其后一倍焦距处产生图像。这个过程中,第一透镜的作用既是所谓的二维付理叶变换,它将物体发射的连续频率的光波折射后在一倍焦距处形成一个以空间频率顺序分布的图形(干涉图形),其中心点光波频率最低,并以此中心由近及远的顺序依次排列低频率至高频率光波。由此可见,该平面即为一个二维K-空间,也称付理叶平面。

图3 透镜成像原理示意图。1f为一倍焦距,第一透镜后一倍焦距光线分布平面即为K-空间。

  第二透镜将K-空间发射的光线处理后形成图像。它将K-空间中每一点所代表的光波按其空间频率和相位的不同进行还原,并相互叠加在一点产生图像。这个过程称为付理叶逆变换或反向付理叶变换。
  在以上透镜成像中,图像的特征如大小、空间分辨力、对比度等等在很大程度上取决于付理叶变换的过程,如图像的分辨力与透镜的大小成正比,即透镜越大,图像的分辨力越高。
  MR成像与透镜成像类似,也分2个过程:首先,射频激发样体在二维梯度场的作用下,每次产生一组具有Ky空间频率相同,而Kx空间不同的MR信号,并按其Kx空间频率顺序依次写入一行K-空间线,改变相位编码梯度磁场强度,Ky值随之改变,每改变一个Ky值填写一行K-空间线,直至所有Ky对应的K-空间线完全填写(图4),便产生了一个由MR信号组成的二维K-空间或二维付理叶平面[3]。二维付理叶平面与光学成像的付理叶平面相似,其内的每一点都对应于一列具体的波,因此将K-空间的数据进行逆向付理叶变换,即将所有的波相互叠加产生MR图像,这个过程与透镜产生图像的过程完全相同。只是MRI中K-空间数据的填写可人为控制。

图4 SE序列中K-空间的填写顺序示意图。每一个TR填写一行K-空间线(Ky)。

  无论是MR图像还是光学影像,图像上每一点与K-空间内每一点不是一一对应关系,图像上每一点的信号来源于K-空间所有点;K-空间内每一点参与图像上所有点信号的形成。
  当我们使用照像机照像时,缩小照像机的光圈后,原来可以成像的物体,缩小光圈后仍可成像,但其亮度下降,如果在镜头正中央贴上一片不透光物,则图像上无法显示对比,只能显示边缘影像。由此可见,镜头的边缘区域(高频率区域)信号主要影响图像的分辨力,而中心区域(低频率区域)则决定图像的对比,且镜头越大则图像越清晰。
  同样,MR成像中K-空间的中央区域(低频率区域)数据决定MR图像的对比,而边缘区域(高频率区域)则影响图像的空间分辨力,且K-空间越大(包括的空间频率范围越大)则图像越清晰。因此我们可以借用光学成像原理来理解抽象的K-空间对MR图像质量的控制作用,并通过对K-空间的控制,优化MR图像质量。

K-空间不同填写方式对图像质量的影响
  1.自旋回波序列中K-空间对图像的影响
  在SE序列中,一个90°激发脉冲后填写一行K-空间线,K-空间线的长度与读出梯度和其持续时间的积成正比,K-空间线的位置由相位编码梯度决定,所以填写一个K-空间所需的时间为N×TR,其中N为相位编码行数,TR为脉冲重复时间。由于每填写一条K-空间线需要一个TR时间,减少相位编码行数(K-空间线)即可相应减少成像时间。如相位编码行从512减至256,则成像时间节约50%,如减到128,则成像时间降到25%。
  K-空间线的减少对图像质量所产生的影响,取决于K-空间线填写的两种不同方式。如果保持相位编码梯度磁场的增幅不变,则相位编码线即K-空间线的距离保持不变,当K-空间线从512减至128时,相位编码方向的K-空间减小75%,所有K-空间线都集中K-空间中心部,由于K-空间的大小控制图像的分辨力,所以图像的空间分辨力下降。但是,K-空间只是在相位编码方向减少,而频率编码方向未变,所以空间分辨力也只有一维的下降,因为图像大小并没有变化。
  如果保持K-空间的大小不变,增加K-空间线间距,则由于K-空间大小决定图像分辨力,所以图像分辨力不变,但图像的大小(FOV)发生了改变,即相位编码方向FOV减至25%,信噪比也下降了50%。这是由于K-空间的不同填写方式所致的两种完全不同的结果,FOV与K-空间的关系如下:
  FOV∝1/D,其中D为K-空间线间距。
  2.重叠影
  在光学成像中,付理叶平面的数据是连续的频率分布,而MR图像的K-空间数据却是不连续的节段有限频率,这就相当于在成像透镜表面放置一个单频散射光栅(图5),由于光栅的散射效应,使物体影像两端产生多个“复制影像”,技术上称之为高序影像(higher orders)。高序影像之间的距离与光栅间距呈反比,即栅间距越大则高序影像距离越近,反之亦然。

图5 光栅散射效应产生高序影像(伪影)。

  MR图像的K-空间数据,由于不是连续频率,所以也存在“光栅效应”。在MR成像中也经常见到“高序”影像,即我们平时所见的“回卷”伪影。K-空间线的间距与高序影像间呈反比,所以减小K-空间线间距可使高序影像与主图像间距增加,从而使其消失在视野之外。其方法有:增加K-空间线数或增加FOV。
  K-空间与影像之间还存在如下关系:
  Kmax∝1/d    Nk=Ni
  Kmax表示最大空间频率,d则表示像素行间距,即最大空间频率越高则像素行间距越小,即空间分辨力越高;Nk表示K-空间线数,Ni表示像素行数,K-空间线数越多,则空间分辨力越高。
  3.半付理叶变换及半回波序列中K-空间对图像的影响
  在前面的讨论中已知,可以通过减少K-空间线以减少采集时间,但减少K-空间大小会使空间分辨力下降;而增加K-空间线的距离则导致FOV减少。如果仔细认识一下K-空间的数据特点,就会发现,K-空间的数据不是孤立的,而是一个对称的频率空间,例如,如果知道了K-空间中K(3,4),则一定还有(-3,-4)、(3,-4)及(-3,4)三个点存在,也就是说,如果知道了一个相限K-空间(1/4K-空间)的数据,理论上就可以完全推算出整个K-空间的数据,这样就不必去耗费大量时间填写K-空间,又可使K-空间大小保持不变,既提高了成像速度又保持了图像的空间分辨力,但由于数据采集量减少,所以信噪比下降。所以一般采用填写一半K-空间,而另一半则根据对称法则推算出来。由于利用对称法则,所以K-空间的填写采用上下1/2或左右1/2,这两种K-空间填写方式在MR中方法不同,前一种为部分相位编码采集,每一次采集一个完整回波信号,一共填写一半K-空间线,称为半付理叶变换采集,或半激发次数(1/2NEX)采集;而后一种为部分频率编码采集,每次仅采集半个回波或3/4个回波,所以又称部分回波或3/4激发次数(3/4NEX)采集[4]。这两种K-空间数据填写,一般都应超过一半K-空间,由于K-空间中心部数据对图像的贡献至关重要。
  半付理叶变换减少K-空间线的填写而提高速度,部分回波技术则是通过缩短回波时间,从而缩短重复时间TR,以提高成像速度。

K-空间填写对图像性质的影响
  在SE序列中由于TE相同,所以无论K-空间线填写顺序如何,中心K-空间线对应的回波时间都不变,但在快速多回波SE序列(TSE)中,每个回波的TE都不相同,如果K-空间填写的顺序不同,则产生的图像性质也不同。如中心K-空间线对应的回波时间TE为20ms,则图像为质子加权图像,若TE为100ms,则图像为T2加权像。一般序列中,将相位编码梯度为0时对应回波写在中心K-空间线,若TSE序列有128个回波,第一回波时间为20ms,按照相位编码由低至高的顺序,相位编码梯度为0对应的回波时间TE=64×20=1280ms(有效回波时间,Teff),则图像为极度T2加权,可能只有边缘增强效应信号。因此,在K-空间填写时,应根据需要设定填写顺序,如需要质子加权像则将第一回波填入中心K-空间线;如需要T2加权像则将回波时间为100或120ms左右的回波时间填入中心K-空间线。由于K-空间线的位置是由相位编码梯度决定,所以可通过控制相位编码梯度,随意安排K-空间线的写入顺序,以满足我们对图像特定性质的要求[5]。
  综上所述,MR成像与其它影像学的主要区别在于,MR操作者对K-空间数据的控制力。通过对K-空间数据的控制,达到控制图像对比度、分辨力、FOV、成像速度。通过改变K-空间数据的写入方式,可以控制伪影,换取对比度、分辨力。作为MRI工作者应对抽象的K-空间概念有一个较具体的认识,并进一步理解K-空间对MR图像的影响作用及其原因。

作者单位:430030 武汉,同济医科大学附属协和医院放射科MR室

参考文献
 1 Mezrich R.A Perspective on K-Space.Radiology,1995,195∶297-315.
 2 Edelman RR,Kleefield J,Wentz KU,et al.Basic Principles of Magnetic Rasonance Imaging.in:Edelman RR and Hesselink JR.Clinical Magnetic Resonance Imaging.Philadelphia W.B.Saunders Company,1990.23.
 3 Martin JF,Edelman RR.Fast Imaging.in:Edelman RR and Hesselink JR.Clinical Magnetic Resonance Imaging.Philadelphia W.B.Saunders Company,1990.214-215.
 4 Twieg DB.The K-trajectory Formulation of NMR Imaging process with applications in analysis and synthesis of Imaging methods.Med Phys,1983,10∶620-621.
 5 Petersson JS,Christoffersson JO,Golman K.MRI Simulation using the K-space Formalism.Magn Rason Imaging,1993,11∶557-558.

转载自http://www.qnr.cn/med/data/yxyx/yxyx/201003/380429.html

K-空间 在图像变换及重建相关推荐

  1. matlab mri的k空间,理解磁共振K空间,自己动手还原和处理K空间数据

    本文内容.图片与涉及的源码均为作者原创,未经许可不得转载.版权声明或联系作者请移步 "关于" 继续计算机与医学影像的跨界之旅,本文所述内容是理解磁共振成像的K空间,并利用计算机还原 ...

  2. 了不起的K空间-cover懋式百科全书

    前言:该文转自"懋式百科全书",作者李懋,飞利浦磁共振首席培训工程师,工作出差培训,全国满天飞,坐飞机如打的士的优质精英商务男士. 这次我们要来聊聊一个看起来高大上,又悬而又悬的名 ...

  3. Albumentations 中的空间级图像变换

    GitHub 官方网站 Featurize Albumentations 中的空间级图像变换 Albumentations 在计算机视觉的图像增强中使用频率很高,无论是分类人物.语义分割人物.目标检测 ...

  4. 半导体器件物理【5】固体量子 —— 能带与k空间

    前言 没时间了,进度拖了两周 ! 好想游泳! 目录 前言 量子数 泡利不相容原理 核外电子排布 能级分裂成能带 允带禁带价带导带 固体中的导带价带 Kronig-Penney模型 单个原子势函数 一维 ...

  5. <Principles of fMRI 1>课程笔记7 K空间

    这一模块的主要作用是增强学生对K空间(K space)与图像空间(image space)之间关系的感性理解.由于K空间的数据分布实际上是图像空间中数据的二维傅立叶变换结果,因此如果你对傅立叶变换比较 ...

  6. 采样频率和带宽的关系_图像视野(FOV)与K空间的关系

    K空间对于整个磁共振成像来说算是最难的部分了,在前文啃一啃硬骨头:K空间我们说过了K空间是一个二维频率空间,在我们的思维中可能有这么一种想法:K空间的频率编码轴对应图像二维视野的编码轴,相位编码轴对应 ...

  7. 基于改进的k最近邻算法的单体型重建问题An Improved KNN Algorithm for Haplotype Reconstruction Problem

    基于改进的k最近邻算法的单体型重建问题 An Improved KNN Algorithm for Haplotype Reconstruction Problem DOI: 10.12677/csa ...

  8. K空间的理解 倒空间

    魏格纳塞茨原胞就是第一布里渊区 一,晶格的Fourier变换 一个物理问题既可以在正(坐标)空间描写,也可以在倒(动量)空间描写 坐标表象,动量表象         为什么选择不同的表象?为什么是动量 ...

  9. oracle 10G 表空间迁移 索引需要重建

    昨天进行了表空间迁移的维护,维护后出现一个问题:MYTEST用户下的表表空间迁移后索引存在失效. 针对该问题和其他DBA进行了一些沟通,大家理解并不统一.后通过实际测试发现,表空间迁移后没有数据的表索 ...

最新文章

  1. (原+转)linux安装rtl 8812au驱动
  2. 修改代码的艺术----- 2.2 高层测试 2.3 测试覆盖
  3. 工厂方法模式适用场景
  4. cocos3.2中如何创建一个场景
  5. 网络编程常用接口的内核实现----sys_listen()
  6. 3. GD32F103C8T6 GPIO的位带操作
  7. 使用 C-JDBC 给 Mysql 集群
  8. 东大OJ-1544: GG的战争法则
  9. 剑指offer32-III.从上到下打印二叉树(蛇形打印)
  10. 解决虚拟机中桥接模式下ip地址的设置
  11. MySQL索引之全文索引(FULLTEXT)
  12. USB(十)2022-03-03
  13. js中利用prompt和parseFloat来实现用户体温华氏和摄氏的提取(18)
  14. 世界上最全的防醉酒+解酒法(为了关心的人,一定要看看哦)
  15. PTA 6-10 二分查找
  16. imagemagick安装问题
  17. 2023安徽农业大学计算机考研信息汇总
  18. 大学生新生报到管理系统
  19. GhostNet代码解析
  20. Pausing and Resuming an Activity 暂停和恢复活动

热门文章

  1. Android 生成自己的签名key(releasekey platform shared media networkstack verify等)
  2. 恶搞 GitHub 贡献图
  3. 原生Ajax五个基本步骤(面试题)
  4. 职场小白新建SSM项目
  5. 原生JS写三级联动 --省--市--区
  6. android更换开机动画,修改安卓开机动画(除了部分系统 如MIUI等)
  7. php html5 cms,建站教程|CMS教程|PHP教程|html5教程 - 站长图库
  8. Please disable your ad blocker to download
  9. 灰色预测模型--两秒直接上手
  10. Pandas数据视图 — groupby 和 pivot