CN110058185B - 一种磁共振旋转成像方法及*** - Google Patents
一种磁共振旋转成像方法及*** Download PDFInfo
- Publication number
- CN110058185B CN110058185B CN201910262523.6A CN201910262523A CN110058185B CN 110058185 B CN110058185 B CN 110058185B CN 201910262523 A CN201910262523 A CN 201910262523A CN 110058185 B CN110058185 B CN 110058185B
- Authority
- CN
- China
- Prior art keywords
- echo
- gradient
- nth
- encoding
- rotation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 110
- 239000013598 vector Substances 0.000 claims abstract description 93
- 238000012937 correction Methods 0.000 claims abstract description 67
- 238000002592 echocardiography Methods 0.000 claims abstract description 23
- 238000012790 confirmation Methods 0.000 claims abstract description 18
- 238000010606 normalization Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 8
- 230000004807 localization Effects 0.000 claims description 20
- 238000000034 method Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 4
- 230000005415 magnetization Effects 0.000 description 7
- 238000012986 modification Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- High Energy & Nuclear Physics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种磁共振旋转成像方法及***。该方法包括如下步骤:获取任意相邻两个回波之间的旋转角;将旋转角反向旋转,以计算梯度场的旋转矩阵;根据旋转矩阵,确认回波与主磁场之间的相对空间定位信息,并对其进行归一化处理,以得到回波的修正向量;根据空间定位确认步骤得到的修正向量以及回波的梯度脉冲,计算回波的实际物理梯度;根据回波的翻转角以及梯度值修正步骤中得到的实际物理梯度,填充K空间。本发明通过翻转角和修正向量计算的实际物理梯度填充K空间,以得到完整的K空间数据,以此主磁场和梯度场的旋转效果被消除,避免了消除了各个相位之间的差异,以避免旋转造成的伪影的出现,减轻或者消除了伪影对扫描图像的影响。
Description
技术领域
本发明涉及磁共振技术领域,具体而言,涉及一种磁共振旋转成像方法及***。
背景技术
磁共振产生的信号由原子核的自旋运动产生。一般情况下,用于人体磁共振成像的是氢质子,每个质子自旋都会产生一个小磁场,人体中有着数以亿计的这种小磁场,正常情况下,由于这种小磁场的分布是杂乱随机的,每个质子自旋产生的磁场相互抵消,因此人体没有磁性。
当这种带小磁场的质子进入一个大的外加磁场时,这种杂乱的质子分布便趋于平行磁场。由于与磁场方向同向的质子数量大于与磁场方向相反的质子,人体组织便产生了一个与主磁场方向一致的宏观纵向磁化量。
从宏观角度来说,磁共振便是通过射频脉冲激发,使纵向磁化量发生偏转。射频脉冲的能量越大,纵向磁化量的偏转越大。纵向的磁化量偏转后产生的横向磁化量切割接收线圈,产生信号接收填充数据K空间,生成磁共振影像。
在射频脉冲的作用下产生的信号来自整个成像空间,没有把成像空间区分,无法形成有效的磁共振影像。为了实现核磁共振成像,必须把收集到的信号进行空间定位,通常采用的空间梯度编码实现这一目的。MRI的梯度磁场有3种:选层梯度场Gz、频率编码梯度场Gx、相位编码梯度场Gy。
正如上述所言,产生磁共振信号的条件为纵向磁化量偏移主磁场方向,一幅完整的磁共振图像,需要N次相位编码对K空间进行填充。N次重复扫描过程中,主磁场和梯度场的稳定是保证信号准确无伪影的前提。而在一些特定的应用场景(如磁共振引导的手术导航或放疗等情况),会导致病患在扫描的过程中,相对主磁场发生相对旋转(主磁场旋转或病患旋转),这会对纵向磁化的方向产生额外的影响,即旋转的磁场会导致相位编码的相位存在差异,对成像带来伪影。
发明内容
鉴于此,本发明提出了一种磁共振旋转成像方法及***,旨在解决现有患者和主磁场之间发生相对旋转时导致相位编码的相位存在差异对成像带来伪影的问题。
一方面,本发明提出了一种磁共振旋转成像方法,该方法包括如下步骤:旋转角获取步骤,获取任意相邻两个回波之间的旋转角;旋转矩阵计算步骤,将所述旋转角获取步骤获取的旋转角反向旋转,以计算梯度场的旋转矩阵;空间定位确认步骤,根据所述旋转矩阵计算步骤得到的旋转矩阵,确认回波与主磁场之间的相对空间定位信息,并对其进行归一化处理,以得到回波的修正向量;梯度值修正步骤,根据所述空间定位确认步骤得到的修正向量以及回波的梯度脉冲,计算回波的实际物理梯度;填充步骤,根据回波的翻转角以及所述梯度值修正步骤中得到的实际物理梯度,填充K空间。
进一步地,上述磁共振旋转成像方法,在所述旋转矩阵计算步骤中,如果患者或主磁场方向绕***的z轴旋转,第n个回波的旋转矩阵Rzn利用如下公式计算:
如果患者或主磁场方向绕***的x轴旋转,第n个回波的旋转矩阵Rxn利用如下公式计算:
如果患者或主磁场方向绕***的y轴旋转,第n个回波的旋转矩阵Ryn利用如下公式计算:
其中,αn-1为第n-1个回波和第n个回波之间的旋转角,n=2,3,4…N,N为磁共振信号中回波的个数。
进一步地,上述磁共振旋转成像方法,在所述空间定位确认步骤中,所述空间定位信息包括:选层矩阵、频率编码矩阵和相位编码矩阵;利用如下公式计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz);第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz);第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz)。
进一步地,上述磁共振旋转成像方法,在所述空间定位确认步骤中,所述修正向量利用如下公式计算:
其中,Tzn为第n个回波的选层修正向量;Txn为第n个回波的频率编码修正向量;Tyn为第n个回波的相位编码修正向量。
进一步地,上述磁共振旋转成像方法,在所述梯度值修正步骤中,所述回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
进一步地,上述磁共振旋转成像方法,在所述修正翻转角步骤中,第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;
其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
另一方面,本发明还提出了一种磁共振旋转成像***,该***包括:旋转角获取模块,用以获取任意相邻两个回波之间的旋转角;旋转矩阵计算模块,其与所述旋转角获取模块电连接,用以接收所述旋转角获取模块获取的旋转角并将所述旋转角反向旋转,以计算梯度场的旋转矩阵;空间定位确认模块,其与所述旋转矩阵计算模块电连接,用以接收所述旋转矩阵计算模块计算的旋转矩阵,根据所述旋转矩阵确认与主磁场的相对空间定位信息,并对其进行归一化处理,以得到回波的修正向量;梯度值修正模块,其与所述空间定位确认模块电连接,用以接收所述空间定位确认模块得到的修正向量,并根据所述修正向量及回波的梯度脉冲,计算回波的实际物理梯度;填充模块,其与所述梯度值修正模块电连接,用以接收所述梯度值修正模块计算的实际物理梯度,并根据回波的翻转角以及所述实际物理梯度,填充K空间。
进一步地,上述磁共振旋转成像***,在所述旋转矩阵计算模块中,如果患者或主磁场方向绕***的z轴旋转,第n个回波的旋转矩阵Rzn利用如下公式计算:
如果患者或主磁场方向绕***的x轴旋转,第n个回波的旋转矩阵Rxn利用如下公式计算:
如果患者或主磁场方向绕***的y轴旋转,第n个回波的旋转矩阵Ryi利用如下公式计算:
其中,αn-1为第n-1个回波和第n个回波之间的旋转角,i=2,3,4…N,N为磁共振信号中回波的个数。
进一步地,上述磁共振旋转成像***,在所述空间定位确认模块中,所述空间定位信息包括:选层矩阵、频率编码矩阵和相位编码矩阵;
利用如下公式计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz);第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz);第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz)。
进一步地,上述磁共振旋转成像***,在所述空间定位确认模块中,所述修正向量利用如下公式计算:
其中,Tzn为第n个回波的选层修正向量;Txn为第n个回波的频率编码修正向量;Tyn为第n个回波的相位编码修正向量。
进一步地,上述磁共振旋转成像***,在所述梯度值修正步骤中,所述回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
进一步地,上述磁共振旋转成像***,该***还包括:修正翻转角模块,其与所述旋转角获取模块电连接,用以接收所述旋转角获取模块获取的旋转角,并在射频线圈与***的中心磁场之间发生偏转时,根据所述旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿,并将补偿后的翻转角信息发送给所述填充模块。
进一步地,上述磁共振旋转成像***,在所述修正翻转角模块中,第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;
其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
本发明提供的磁共振旋转成像方法及***,通过计算各个回波与主磁场之间的相对空间定位信息,以便计算各个回波的修正向量,进而以此对回波的梯度脉冲进行修正,从而通过翻转角和修正向量计算的实际物理梯度填充K空间,以得到完整的K空间数据,以此主磁场和梯度场的旋转效果被消除,得到的K空间数据具有一致的空间定位信息,进而对其进行傅里叶变换以得到磁共振旋转图像,从而避免了消除了各个相位之间的差异,以避免旋转造成的伪影的出现,减轻或者消除了伪影对扫描图像的影响,并且提高了图像的信噪比和对比度,以保证扫描结果的准确,提高扫描图像的质量。
进一步地,为避免旋转造成导致翻转角的不足,翻转角不准确将导致信噪比和对比度的差异,优选地,在射频线圈与***的中心磁场之间发生偏转时,根据所述旋转角获取步骤获取的旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿,以此填充K空间,进而进一步避免旋转影响成像的准确性,以进一步提高扫描结果的准确,提高扫描图像的质量。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本发明实施例提供的主磁场相对于患者进行旋转的示意图;
图2为本发明实施例提供的患者相对于主磁场进行旋转的示意图;
图3为本发明实施例提供的旋转K空间的示意图;
图4为本发明实施例提供的磁共振旋转成像方法的流程示意图;
图5为本发明实施例提供的磁共振旋转成像***的结构框图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
参见图1和图2,其示出了本实施例中患者和主磁场之间发生相对旋转的示意图。常规磁共振扫描中,如图1a和图2c所示,图中主磁场B0方向1和患者2之间的相对位置固定不动。而在旋转扫描过程中,如图1b和图2d所示,主磁场B0方向1和患者2绕***中心点旋转。
参见图3,其示出了本实施例中旋转K空间的示意图。如图所示,一幅完整磁共振图像通常需要N次相位编码梯度填充K空间,以得到完整的K空间数据。磁共振成像的每一个回波都需要一个激发脉冲和梯度,后续的1,2,3,4一直到N行结束的回波都是激发脉冲产生信号,梯度产生空间定位,并接收数据填充K空间,以便得到完整的K空间数据。
方法实施例:
参见图4,其示出了本发明实施例提供的磁共振旋转成像方法的流程示意图。如图所示,该方法包括如下步骤:
旋转角获取步骤S1,获取任意相邻两个回波之间的旋转角。
具体地,可通过编码器或光栅尺等方式得到任意相邻两个回波之间的旋转角,如图2所示,第n-1个回波和第n个回波之间的旋转角记为αn-1,则其中,n=2,3,4…N,即n为2到N之间的正整数,N为磁共振信号中回波的个数。当然,旋转角亦可通过其他方式得到,本实施例中对其不做任何限定。
修正翻转角步骤S2,在射频线圈与***的中心磁场之间发生偏转时,根据旋转角获取步骤获取的旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿。
具体地,翻转角即为偏转角,其通过射频线圈产生的射频脉冲产生。旋转的磁场导致翻转角的不足,翻转角不准确将导致信噪比和对比度的差异,故为避免信噪比和对比度的差异,在射频线圈与***的中心磁场之间发生偏转时,因为***的中心磁场发生偏转,射频线圈不偏转,所以两者之间会逐步产生一个角度差,为避免其影响成像的准确性,需对该角度差进行修正。
采集中第一次激发的翻转角θ即第一个回波的翻转角θ作为基准不作校正,其他回波在前一个回波基础上进行校正,以完成翻转角不足的补偿。其中,各个回波的翻转角的修正可根据如下公式进行计算:
第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
在射频线圈跟随***同步旋转时,即射频线圈与***的中心磁场之间发生不发生相对偏转,则不需要修正这个角度,即各个回波的修正角均为零,也就是说各个回波的翻转角均为θ。其中,第一个回波的翻转角θ可由用户根据实际情况对其进行设定,本实施例中对其不做任何限定。
旋转矩阵计算步骤S3,将旋转角获取步骤S1获取的旋转角反向旋转,以计算梯度场的旋转矩阵。
具体地,以右手坐标系为例,将旋转角反方向旋转-αn-1(,n=2,3,4…N),可以得到梯度场的旋转矩阵。其中,
如果患者或主磁场方向绕***的z轴旋转,第n个回波的旋转矩阵Rzn利用如下公式计算:
如果患者或主磁场方向绕***的x轴旋转,第n个回波的旋转矩阵Rxn利用如下公式计算:
如果患者或主磁场方向绕***的y轴旋转,第n个回波的旋转矩阵Ryn利用如下公式计算:
其中,αn-1为第n-1个回波和第n个回波之间的旋转角,n=2,3,4…N,N为磁共振信号中回波的个数。
通过如上公式计算各个回波的旋转矩阵,即第二个回波至第N个回波的旋转矩阵,进而计算确认磁共振信号的空间定位信息,从而对回波的实际物理梯度进行计算,以获取完整的K空间数据。
空间定位确认步骤S4,根据旋转矩阵计算步骤得到的旋转矩阵,确认回波与主磁场之间的相对空间定位信息,并对其进行归一化处理,以得到回波的修正向量。
具体地,空间定位信息息包括:选层矩阵、频率编码矩阵和相位编码矩阵。首先,空间定位信息可依据旋转矩阵计算步骤S3计算的旋转矩阵利用如下公式进行计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个,也就是说,如果患者或主磁场方向绕***的z轴旋转,第n个回波的旋转矩阵Rn=Rzn,如果患者或主磁场方向绕***的x轴旋转,第n个回波的旋转矩阵Rn=Rxn,如果患者或主磁场方向绕***的y轴旋转,第n个回波的旋转矩阵Rn=Ryn;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz),由选层梯度Gz在***物理坐标系的三个方位投影得出;第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz),由频率编码梯度场Gx在***物理坐标系的三个方位投影得出;第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz),由相位编码梯度场Gy在***物理坐标系的三个方位投影得出。
然后,对上述计算的空间定位信息息即选层矩阵、频率编码矩阵和相位编码矩阵进行归一化处理,以得到各个回波的修正向量,进而对各个回波的梯度脉冲进行修正。
其中,修正向量可利用如下公式计算:
其中,Tzn为第n个回波的选层修正向量;Txn为第n个回波的频率编码修正向量;Tyn为第n个回波的相位编码修正向量。
通过如上公式计算各个回波的三个修正向量即选层修正向量、频率编码修正向量和相位编码修正向量,即第二个回波至第N个回波的选层修正向量、频率编码修正向量和相位编码修正向量,进而据此进行实际物力梯度的计算。
梯度值修正步骤S5,根据空间定位确认步骤得到的修正向量以及回波的梯度脉冲,计算回波的实际物理梯度。
具体地,根据空间定位确认步骤得到的修正向量以及回波的梯度脉冲,将选层、频率、相位三个方向的梯度脉冲分别乘以其对应的修正向量,以便得到选层、频率和相位的实际物理梯度坐标。其中,回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
通过如上公式计算各个回波的三个实际梯度值即实际选层梯度、实际频率编码梯度和实际相位编码梯度,即第二个回波至第N个回波的实际选层梯度、实际频率编码梯度和实际相位编码梯度,进而据此进行K空间的填充。
填充步骤S6,根据回波的翻转角以及梯度值修正步骤中得到的实际物理梯度,填充K空间。
具体地,根据修正翻转角步骤S2中回波的修正后的翻转角以及梯度值修正步骤S5中得到的实际物理梯度,填充K空间,以得到完整的K空间数据,以此主磁场和梯度场的旋转效果被消除,得到的K空间数据具有一致的空间定位信息,进而对其进行傅里叶变换以得到磁共振旋转图像,从而避免了消除了各个相位之间的差异,以避免旋转造成的伪影的出现,减轻或者消除了伪影对扫描图像的影响,并且提高了图像的信噪比和对比度,以保证扫描结果的准确,提高扫描图像的质量。
在本实施例中,修正翻转角步骤S2与旋转矩阵计算步骤S3至梯度值修正步骤S5之间没有先后顺序。
综上,本实施例提供的磁共振旋转成像方法,通过计算各个回波与主磁场之间的相对空间定位信息,以便计算各个回波的修正向量,进而以此对回波的梯度脉冲进行修正,从而通过翻转角和修正向量计算的实际物理梯度填充K空间,以得到完整的K空间数据,以此主磁场和梯度场的旋转效果被消除,得到的K空间数据具有一致的空间定位信息,进而对其进行傅里叶变换以得到磁共振旋转图像,从而避免了消除了各个相位之间的差异,以避免旋转造成的伪影的出现,减轻或者消除了伪影对扫描图像的影响,并且提高了图像的信噪比和对比度,以保证扫描结果的准确,提高扫描图像的质量。
进一步地,为避免旋转造成导致翻转角的不足,翻转角不准确将导致信噪比和对比度的差异,优选地,在射频线圈与***的中心磁场之间发生偏转时,根据旋转角获取步骤获取的旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿,以此填充K空间,进而进一步避免旋转影响成像的准确性,以进一步提高扫描结果的准确,提高扫描图像的质量。
***实施例:
参见图5,其为本发明实施方式提供的磁共振旋转成像***的结构框图。如图5所示,本发明的实施方式提供的磁共振旋转成像***,包括:旋转角获取模块100、修正翻转角模块200、旋转矩阵计算模块300、空间定位确认模块400、梯度值修正模块500和填充模块600;其中,
旋转角获取模块100,用以获取任意相邻两个回波之间的旋转角。
修正翻转角模块200,其与旋转角获取模块100电连接,用以接收旋转角获取模块100获取的旋转角,并在射频线圈与***的中心磁场之间发生偏转时,根据旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿,并将补偿后的翻转角信息发送给填充模块600。
旋转矩阵计算模块300,其与旋转角获取模块100电连接,用以接收旋转角获取模块100获取的旋转角并将旋转角反向旋转,以计算梯度场的旋转矩阵。
空间定位确认模块400,其与旋转矩阵计算模块300电连接,用以接收旋转矩阵计算模块300计算的旋转矩阵,根据旋转矩阵确认与主磁场的相对空间定位信息,并对其进行归一化处理,以得到回波的修正向量。
梯度值修正模块500,其与空间定位确认模块400电连接,用以接收空间定位确认模块400得到的修正向量,并根据修正向量及回波的梯度脉冲,计算回波的实际物理梯度。
填充模块600,其与修正翻转角模块200和梯度值修正模块500电连接,用以接收修正翻转角模块200修正的翻转角和梯度值修正模块500计算的实际物理梯度,并根据回波的翻转角以及实际物理梯度,填充K空间。
优选地,在旋转矩阵计算模块中300,如果患者或主磁场方向绕***的z轴旋转,第n个回波的旋转矩阵Rzn利用如下公式计算:
如果患者或主磁场方向绕***的x轴旋转,第n个回波的旋转矩阵Rxn利用如下公式计算:
如果患者或主磁场方向绕***的y轴旋转,第n个回波的旋转矩阵Ryi利用如下公式计算:
其中,αn-1为第n-1个回波和第n个回波之间的旋转角,i=2,3,4…N,N为磁共振信号中回波的个数。
优选地,在空间定位确认模块400中,空间定位信息包括:选层矩阵、频率编码矩阵和相位编码矩阵;
利用如下公式计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz);第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz);第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz)。
进一步优选地,在空间定位确认模块400中,修正向量利用如下公式计算:
其中,Tzn为第n个回波的选层修正向量;Txn为第n个回波的频率编码修正向量;Tyn为第n个回波的相位编码修正向量。
优选地,在梯度值修正步骤中,回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
优选地,在修正翻转角模块200中,第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;
其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
本发明的实施例的磁共振旋转成像***与本发明的另一个实施例的磁共振旋转成像方法相对应,在此不再赘述。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (14)
1.一种磁共振旋转成像方法,其特征在于,包括如下步骤:
旋转角获取步骤,获取任意相邻两个回波之间的旋转角;
旋转矩阵计算步骤,将所述旋转角获取步骤获取的旋转角反向旋转,以计算梯度场的旋转矩阵;
空间定位确认步骤,根据所述旋转矩阵计算步骤得到的旋转矩阵,确认回波与主磁场之间的相对空间定位信息,并对回波与主磁场的相对空间定位信息进行归一化处理,以得到回波的修正向量;
梯度值修正步骤,根据所述空间定位确认步骤得到的修正向量以及回波的梯度脉冲,计算回波的实际物理梯度;
填充步骤,根据回波的翻转角以及所述梯度值修正步骤中得到的实际物理梯度,填充K空间。
3.根据权利要求2所述的磁共振旋转成像方法,其特征在于,
在所述空间定位确认步骤中,所述空间定位信息包括:选层矩阵、频率编码矩阵和相位编码矩阵;
利用如下公式计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz);第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz);第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz)。
5.根据权利要求1至4任一项所述的磁共振旋转成像方法,其特征在于,
在所述梯度值修正步骤中,所述回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
6.根据权利要求1至4任一项所述的磁共振旋转成像方法,其特征在于,在所述旋转角获取步骤和所述填充步骤之间还包括如下步骤:
修正翻转角步骤,在射频线圈与***的中心磁场之间发生偏转时,根据所述旋转角获取步骤获取的旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿。
7.根据权利要求6所述的磁共振旋转成像方法,其特征在于,在所述修正翻转角步骤中,
第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;
其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
8.一种磁共振旋转成像***,其特征在于,包括:
旋转角获取模块,用以获取任意相邻两个回波之间的旋转角;
旋转矩阵计算模块,其与所述旋转角获取模块电连接,用以接收所述旋转角获取模块获取的旋转角并将所述旋转角反向旋转,以计算梯度场的旋转矩阵;
空间定位确认模块,其与所述旋转矩阵计算模块电连接,用以接收所述旋转矩阵计算模块计算的旋转矩阵,根据所述旋转矩阵确认回波与主磁场的相对空间定位信息,并对回波与主磁场的相对空间定位信息进行归一化处理,以得到回波的修正向量;
梯度值修正模块,其与所述空间定位确认模块电连接,用以接收所述空间定位确认模块得到的修正向量,并根据所述修正向量及回波的梯度脉冲,计算回波的实际物理梯度;
填充模块,其与所述梯度值修正模块电连接,用以接收所述梯度值修正模块计算的实际物理梯度,并根据回波的翻转角以及所述实际物理梯度,填充K空间。
10.根据权利要求9所述的磁共振旋转成像***,其特征在于,
在所述空间定位确认模块中,所述空间定位信息包括:选层矩阵、频率编码矩阵和相位编码矩阵;
利用如下公式计算:
其中,Mzn为第n个回波的选层矩阵;Mxn为第n个回波的频率编码矩阵;Myn为第n个回波的相位编码矩阵;(Zxn、Zyn、Zzn)为第n个回波的选层编码梯度场的方向向量;(Xxn、Xyn、Xzn)为第n个回波的频率编码梯度场的方向向量;(Yxn、Yyn、Yzn)为第n个回波的相位编码梯度场的方向向量;Rn为第n个回波的旋转矩阵,其对应取值为Rzn、Rxn和Ryn中的一个;第一个回波的选层梯度场的方向向量为(Zx、Zy、Zz);第一个回波的频率编码梯度场的方向向量为(Xx、Xy、Xz);第一个回波的相位编码梯度场的方向向量为(Yx、Yy、Yz)。
12.根据权利要求8至11任一项所述的磁共振旋转成像***,其特征在于,
在所述梯度值修正步骤中,所述回波的实际物理梯度利用如下公式计算:
Gzn=Gzn-1×Tzn,
Gxn=Gxn-1×Txn,
Gyn=Gyn-1×Tyn,
其中,Gzn为第n个回波的实际选层梯度;Gxn为第n个回波的实际频率编码梯度;Gyn为第n个回波的实际相位编码梯度;第一个回波的实际选层梯度场为Gz×(Zx、Zy、Zz);第一个回波的实际频率编码梯度为Gx×(Xx、Xy、Xz);第一个回波的实际相位编码梯度为Gy×(Yx、Yy、Yz);Gz为选层梯度场;Gx为频率编码梯度场;Gy为相位编码梯度场。
13.根据权利要求8至11任一项所述的磁共振旋转成像***,其特征在于,还包括:
修正翻转角模块,其与所述旋转角获取模块电连接,用以接收所述旋转角获取模块获取的旋转角,并在射频线圈与***的中心磁场之间发生偏转时,根据所述旋转角对各个回波的翻转角进行修正,以进行翻转角的补偿,并将补偿后的翻转角信息发送给所述填充模块。
14.根据权利要求13所述的磁共振旋转成像***,其特征在于,在所述修正翻转角模块中,
第n-1个回波的翻转角θn-1作为基准,第n个回波的翻转角校正为θn=θn-1-αn-1;
其中,第一个回波的翻转角为θ;αn-1为第n-1个回波和第n个回波之间的翻转角;n=2,3,4…N;N为磁共振信号中回波的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910262523.6A CN110058185B (zh) | 2019-04-02 | 2019-04-02 | 一种磁共振旋转成像方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910262523.6A CN110058185B (zh) | 2019-04-02 | 2019-04-02 | 一种磁共振旋转成像方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110058185A CN110058185A (zh) | 2019-07-26 |
CN110058185B true CN110058185B (zh) | 2021-01-15 |
Family
ID=67318152
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910262523.6A Active CN110058185B (zh) | 2019-04-02 | 2019-04-02 | 一种磁共振旋转成像方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110058185B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116152367B (zh) * | 2022-11-02 | 2023-10-20 | 佛山瑞加图医疗科技有限公司 | 一种基于弥散加权的自旋回波磁共振成像方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101063712A (zh) * | 2005-12-20 | 2007-10-31 | 西门子公司 | 用于在磁共振数据中确定横向松弛时间t2*的方法 |
CN101357063A (zh) * | 2008-08-29 | 2009-02-04 | 华东师范大学 | 一种磁共振快速自旋回波成像方法 |
CN103584864A (zh) * | 2012-08-15 | 2014-02-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种磁共振成像方法和装置 |
CN103777162A (zh) * | 2012-10-19 | 2014-05-07 | 上海联影医疗科技有限公司 | 磁共振成像k空间运动伪影矫正并行采集重建方法 |
WO2014163370A1 (ko) * | 2013-04-01 | 2014-10-09 | 고려대학교 산학협력단 | 하이브리드 자기 공명 영상 처리 장치 및 방법 |
CN105044636A (zh) * | 2014-04-16 | 2015-11-11 | 西门子公司 | 用于补偿磁场的不均匀性的方法、***和磁共振设备 |
CN109477877A (zh) * | 2016-07-14 | 2019-03-15 | 上海联影医疗科技有限公司 | 磁共振成像***和方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8064982B2 (en) * | 2006-11-21 | 2011-11-22 | Battelle Memorial Institute | Methods for magnetic resonance analysis using magic angle technique |
-
2019
- 2019-04-02 CN CN201910262523.6A patent/CN110058185B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101063712A (zh) * | 2005-12-20 | 2007-10-31 | 西门子公司 | 用于在磁共振数据中确定横向松弛时间t2*的方法 |
CN101357063A (zh) * | 2008-08-29 | 2009-02-04 | 华东师范大学 | 一种磁共振快速自旋回波成像方法 |
CN103584864A (zh) * | 2012-08-15 | 2014-02-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种磁共振成像方法和装置 |
CN103777162A (zh) * | 2012-10-19 | 2014-05-07 | 上海联影医疗科技有限公司 | 磁共振成像k空间运动伪影矫正并行采集重建方法 |
WO2014163370A1 (ko) * | 2013-04-01 | 2014-10-09 | 고려대학교 산학협력단 | 하이브리드 자기 공명 영상 처리 장치 및 방법 |
CN105044636A (zh) * | 2014-04-16 | 2015-11-11 | 西门子公司 | 用于补偿磁场的不均匀性的方法、***和磁共振设备 |
CN109477877A (zh) * | 2016-07-14 | 2019-03-15 | 上海联影医疗科技有限公司 | 磁共振成像***和方法 |
Non-Patent Citations (2)
Title |
---|
MRI of the coronary arteries: flip angle train optimization for 3D sequences;VolkerHietschold;《Zeitschrift für Medizinische Physik》;20020331;第12卷(第3期);第177-181页 * |
功能磁共振图像小波域分析结果的修正;支联合;《中国医学影像技术》;20060131;第22卷(第1期);第31-34页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110058185A (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6842001B2 (en) | MRI systems with parallel receivers for phase correction | |
US7323873B2 (en) | Phase error measuring method, MR imaging method, and MRI system | |
EP1585991B1 (en) | Magnetic resonance method and device | |
US8934691B2 (en) | System for motion compensated MR parallel imaging | |
US9223001B2 (en) | MR imaging using navigators | |
US8890530B2 (en) | Magnetic resonance imaging method, magnetic resonance imaging apparatus, and control device of magnetic resonance imaging apparatus | |
WO2010035569A1 (ja) | 磁気共鳴イメージング装置 | |
EP0965855A2 (en) | Correction of errors in magnetic resonance images due to magnetic gradient non-uniformity | |
US10048343B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
US10705172B2 (en) | Magnetic resonance apparatus and method for dynamic adjustment thereof with multiple adjustment parameters | |
US6611143B2 (en) | MRI apparatus generating odd and/or even echo images with sensitivity distribution of coils | |
US11194003B2 (en) | MRI apparatus and MRI method | |
US20120025825A1 (en) | Method and magnetic resonance apparatus for dynamic phase correction in a multi channel rf transmission | |
US20110115487A1 (en) | Method and magnetic resonance system for imaging particles | |
CN110058185B (zh) | 一种磁共振旋转成像方法及*** | |
CN109983358B (zh) | Propeller mr成像 | |
CN109716155B (zh) | 具有迪克逊型水/脂肪分离的mr成像 | |
US10746832B2 (en) | DTI with correction of motion-induced diffusion gradient inconsistency | |
CN110730912A (zh) | 并行多切片mr成像 | |
CN111693916B (zh) | 用于控制磁共振成像***的方法和装置 | |
CN111164444B (zh) | 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像 | |
CN114902062A (zh) | 具有内部运动检测的3d mr成像 | |
Park et al. | Three-Dimensional Variable Slab-Selective Projection Acquisition Imaging | |
EP3480617A1 (en) | Diffusion weighted turbo spin echo mr imaging with motion compensation | |
US11892532B2 (en) | Whole-blade acquisition and phase correction in magnetic resonance imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |