CN106680750A - 磁共振匀场图像获取方法、匀场方法及磁共振系统 - Google Patents
磁共振匀场图像获取方法、匀场方法及磁共振系统 Download PDFInfo
- Publication number
- CN106680750A CN106680750A CN201611250997.1A CN201611250997A CN106680750A CN 106680750 A CN106680750 A CN 106680750A CN 201611250997 A CN201611250997 A CN 201611250997A CN 106680750 A CN106680750 A CN 106680750A
- Authority
- CN
- China
- Prior art keywords
- shimming
- length
- magnetic resonance
- minimum
- image
- 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.)
- Granted
Links
Classifications
-
- 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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
-
- 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
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种磁共振匀场图像获取方法、磁共振系统匀场方法及磁共振系统,磁共振匀场图像获取方法包括:获取成像对象的预扫描图像;基于待匀场区域确定候选方向集合;根据所述预扫描图像计算所述成像对象在各所述候选方向的第一长度;计算待匀场区域在各所述候选方向的第二长度;基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度,确定匀场扫描的相位编码方向和相位编码方向的视野;由所述匀场扫描获取包含所述待匀场区域的匀场图像。本发明的技术方案可以在较短的时间内采集较高分辨率的B0场,采集到的图像无卷褶,大大提高了B0匀场的准确性和稳定性。
Description
技术领域
本发明涉及医学图像领域,特别涉及磁共振匀场调节方法及系统。
背景技术
在磁共振成像系统中,主磁场(B0场)的均匀性对图像质量有非常重要的影响。例如B0的均匀性直接影响压脂图像的压脂效果,B0场不均匀会导致压脂图像中脂肪亮起、肌肉组织发黑或者压脂不均匀等现象的产生。
BO场的不均匀一方面是因为磁体本身的B0具有一定空间分布,随着偏离中心程度的增加,B0不均匀度随之增加;另一方面源于被扫描患者进入磁体后自身带来的磁化率分布改变了B0场的分布。
为提高临床扫描区域内的B0场的均匀程度,就需要针对扫描区域进行B0匀场。B0匀场一般先通过序列采集需要匀场区域内的B0场图,然后通过匀场算法对B0场图做球谐函数的拟合展开,最后将拟合所得的各阶系数转化为各方向的梯度偏置和二阶线圈的电流,反馈到系统中,从而达到匀场的效果。
发明人在实现本发明的过程中发现,在上述匀场过程中,第一步采集场图尤其重要,因为这是匀场算法赖以分析计算的数据。一方面,采集的B0场图需要有较高的分辨率,这决定了匀场结果的准确性和稳定性;另一方面,采集B0的过程属于扫描之前的校准项,不应花费太长时间而影响用户体验。然而,高分辨率和快速采集B0场图本身是一对矛盾体。
发明内容
本发明要解决的问题是提供一种磁共振匀场图像获取方法、磁共振系统匀场方法及磁共振系统,缓解匀场图像采集速度慢并且分辨率不高的问题。
为解决上述问题,本发明提供了一种磁共振匀场图像获取方法,包括:获取成像对象的预扫描图像;基于待匀场区域确定候选方向集合,所述候选方向集合包括至少两个候选方向;根据所述预扫描图像计算所述成像对象在各所述候选方向的第一长度;计算待匀场区域在各所述候选方向的第二长度;基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度,将各所述最小无卷褶长度中的最小值所对应的候选方向作为所述匀场扫描的相位编码方向,将所述相位编码方向对应的最小无卷褶长度作为所述匀场扫描相位编码方向的视野;由所述匀场扫描获取包含所述待匀场区域的匀场图像。
优选的,所述基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度包括:计算所述第一长度和第二长度的重叠区域;基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度。
优选的,所述基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度包括:以所述第一长度和第二长度之和减去所述重叠区域作为最小无卷褶长度。
优选的,还包括:将所述待匀场区域在非所述相位编码方向的第二长度作为匀场扫描在非相位编码方向的视野。
优选的,所述成像对象为水模样品。
为解决上述问题,本发明还提供了一种磁共振系统匀场方法,包括:根据磁共振匀场图像获取方法获取包含待匀场区域的匀场图像;基于所述匀场图像计算匀场参数,根据所述匀场参数进行磁共振系统匀场。
优选的,所述基于所述匀场图像计算匀场参数包括:基于所述匀场图像进行球函数展开,计算球函数的各阶系数;基于所述各阶系数计算匀场参数。
优选的,所述匀场参数是关于梯度偏置和二阶电流的匀场参数。
为解决上述问题,本发明还提供了一种磁共振系统,包括:存储器,用于存储可由处理器执行的指令;处理器,用于执行指令,所述指令包括:获取样品的预扫描图像;基于待匀场区域确定候选方向集合,所述候选方向集合包括至少两个候选方向;根据所述预扫描图像计算所述样品在各所述候选方向的第一长度;计算待匀场区域在各所述候选方向的第二长度;基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度,将各所述最小无卷褶长度中的最小值所对应的候选方向作为所述匀场扫描的相位编码方向,将所述相位编码方向对应的最小无卷褶长度作为所述匀场扫描相位编码方向的视野;由所述匀场扫描获取包含所述待匀场区域的匀场图像。
优选的,所述基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度包括:计算所述第一长度和第二长度的重叠区域;基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度。
与现有技术相比,本发明的技术方案提供一种磁共振匀场图像获取方法、磁共振系统匀场方法及磁共振系统,可以在较短的时间内采集较高分辨率的B0场,采集到的图像无卷褶,大大提高了B0匀场的准确性和稳定性。进一步的,特别是提升了小关节压脂的压脂图像效果及稳定性。
附图说明
图1是本发明一些实施例的磁共振系统的结构示意图;
图2是本发明一些实施例的磁共振匀场图像获取方法的流程图;
图3是本发明一些实施例的预扫描图像的示意图;
图4是本发明一些实施例的计算候选方向对应长度的示意图;
图5是本发明一些实施例的匀场扫描视野的示意图;
图6是本发明一些实施例的匀场扫描脉冲序列的示意图;
图7是本发明一些实施例的根据采集到的k空间数据重建磁共振图像的示意图;
图8是是本发明一些实施例的磁共振系统匀场方法的示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
如本申请和权利要求书中所示,除非上下文明确提示例外情形,“一”、“一个”、“一种”和/或“该”等词并非特指单数,也可包括复数。一般说来,术语“包括”与“包含”仅提示包括已明确标识的步骤和元素,而这些步骤和元素不构成一个排它性的罗列,方法或者设备也可能包含其他的步骤或元素。
本申请中使用了流程图用来说明根据本申请的实施例的系统所执行的操作。应当理解的是,前面或下面操作不一定按照顺序来精确地执行。相反,可以按照倒序或同时处理各种步骤。同时,也可以将其他操作添加到这些过程中,或从这些过程移除某一步或数步操作。
图1是本发明一些实施例的磁共振系统的结构示意图,如图1所示,磁共振系统100通常包括磁共振机架,机架内有主磁体101,主磁体101可以是由超导线圈构成,用来产生主磁场,在一些情况下也可以采用永磁体。主磁体101可以用来产生0.2特斯拉、0.5特斯拉、1.0特斯拉、1.5特斯拉、3.0特斯拉或者更高的主磁场强度。在磁共振成像时,成像对象150会由患者床106进行承载,随着床板的移动,将成像对象150移入主磁场磁场分布较为均匀的区域105内。通常对于磁共振系统,如图1所示,空间坐标系(即设备的坐标系)的z方向设置为与磁共振系统机架的轴向相同,通常将患者的身长方向与z方向保持一致进行成像,磁共振系统的水平平面设置为xz平面,x方向与z方向垂直,y方向与x和z方向均垂直。
在磁共振成像,脉冲控制单元111控制射频脉冲产生单元116产生射频脉冲,射频脉冲由放大器放大后,经过开关控制单元117,最终由体线圈103或者局部线圈104发出,对成像对象150进行射频激发。成像对象150根据射频激发,会由共振产生相应的射频信号。在接收成像对象150根据激发产生的射频信号时,可以是由体线圈103或者局部线圈104进行接收,射频接收链路可以有很多条,射频信号发送到射频接收单元118后,进一步发送到图像重建单元121进行图像重建,形成磁共振图像。
磁共振系统100还包括梯度线圈102,梯度线圈可以用来在磁共振成像时对射频信号进行空间编码。脉冲控制单元111控制梯度信号产生单元112产生梯度信号,梯度信号通常会分为三个相互正交方向的信号:x方向、y方向和z方向,不同方向的梯度信号经过梯度放大器(113、114、115)放大后,由梯度线圈102发出,在区域105内产生梯度磁场。
脉冲控制单元111、图像重建单元121与处理器122、显示单元123、输入/输出设备124、存储单元125、通信端口126之间可以通过通信总线125进行数据传输,从而实现对磁共振成像过程的控制。其中,处理器122可以由一个或多个处理器组成。显示单元123可以是提供给用户用来显示图像的显示器。输入/输出设备124可以是键盘、鼠标、控制盒等相关设备,支持输入/输出相应数据流。存储单元125可以是只读存储器(ROM)、随机存取存储器(RAM)、硬盘等,存储单元125可以用来存储需要处理和/或通信使用的各种数据文件,以及处理器122所执行的可能的程序指令。通信端口105可以实现与其他部件例如:外接设备、图像采集设备、数据库、外部存储以及图像处理工作站等之间进行数据通信。
图2是本发明一些实施例的磁共振匀场图像获取方法流程图,参见图2,步骤S201获取成像对象的预扫描图像。成像对象可以是样品(水模)、患者、动物,或者是样品、患者、动物的一部分,例如成像目标仅是患者头部、身体或者四肢。预扫描是对成像对象进行的概略扫描,通常用来反映成像对象的整体情况,例如对成像对象扫描得到冠状位图像或者矢状位图像,预扫描的扫描时间非常短,通常只有几十毫秒,并且预扫描图像的分辨率通常较低。
在对磁场强度进行匀场扫描时,通常会采用一些扫描样本进行扫描,扫描样本可以是水模(phantom)。结合图3,水模300可以是截面为椭圆形的圆柱形状,预扫描水模300得到的预扫描图像可以是如图像301所示,或者图像302所示。在预扫描图像301中成像对象(水模)部分311在图中用灰色底色的网格线进行示意性表示,在预扫描图像302中成像对象(水模)部分312在图中用灰色底色的网格线进行示意性表示。
步骤S202,在匀场扫描时,通常匀场扫描会根据临床扫描来确定匀场扫描的方位,,匀场扫描是对待匀场区域的磁共振扫描,通常需要包含待扫描区域。匀场扫描是临床扫描之前的磁场均匀性校准项目,匀场扫描可以是和临床扫描有相同的方位,例如临床扫描的方位包括方向a、方向b和方向c,则匀场扫描的方位也可以包括方向a、方向b和方向c。对于匀场扫描来说,方向a、方向b、方向c可以作为一个候选方向集合,在其中选取各方向分别作为匀场扫描的相位编码方向(PE,Phase Encoding)、选层方向(SS,Slice Selection)和读出方向(RO,Readout),但具体如何选择并没有固定规则。发明人发现通常RO方向存在过采样,SS为选层方向,通常不会发生卷褶,因此重点需要在候选方向集合中确认哪一个方向作为PE方向。
步骤S203,根据预扫描图像计算成像对象在各候选方向的第一长度。结合图4,由于上一步骤确定的候选方向集合包含三个候选方向:方向a、方向b和方向c,可以计算得到成像对象311对于各候选方向的长度分别为:Lal、Lb1和Lc1(Lc1图4中未示出),这里记为第一长度。
步骤S204,对于待匀场区域,通常是认为磁场均匀性不够,后续需要根据匀场扫描以便开展匀场的区域。待匀场区域可以是完整的B0场区域,例如图1中的区域105,在局部匀场时,待匀场区域也可以只是B0场区域的一小块区域。当匀场扫描的视野确定后,对应的方向a、方向b和方向c也相应确定,因此可以计算待匀场区域对于各候选方向的长度,这里记为第二长度。如图4所示,待匀场区域400在各候选方向可以分别为La2、Lb2和Lc2(Lc2图4中未示出)。
步骤S205,根据之前计算得到的第一长度和第二长度,计算各候选方向的最小无卷褶长度。具体可以根据公式[1]计算对应方向a、方向b、方向c的最小无卷褶长度,分别记为Lma、Lmb和Lmc。
公式[1]:
其中,如图4所示,Δa、Δb、Δc分别为对应a、b、c三个方向上成像对象对应的第一长度(La1、Lb1、Lc1)与匀场扫描视野的第二长度(La2、Lb2、Lc2)的重叠区域。
之后比较各候选方向的最小无卷褶长度:Lma、Lmb、Lmc,将其中的最小值所对应的方向确定为后续匀场扫描的相位编码方向(PE,Phase Encoding),而其余两个方向则分别确定为选层方向(SS,Slice Selection)和读出方向(RO,Readout)。同时,将选取为相位编码方向的候选方向所对应的最小无卷褶长度设定为在匀场扫描时,对应相位编码方向的视野。例如若根据第一长度和第二长度计算得到的最小无卷褶长度Lma相比于Lmb和Lmc最小,则设定方向a为匀场扫描的PE方向,同时设定PE方向的视野FOVPE=Lma,而其余的方向b和方向c可以设定为匀场扫描的SS方向、RO方向。对于方向b和方向c,具体可以是方向b作为RO方向,方向c作为SS方向,也可以是方向b作为SS方向,方向c作为RO方向。
匀场扫描在SS、RO方向的视野可以设为待匀场区域在对应方向的长度,例如若选择方向a为相位编码方向,方向b为选层方向,方向c为读出方向,则FOVSS=Lb2,FOVRO=Lc2。在一些实施例中,匀场扫描在SS、RO方向的视野还可以比待匀场区域在对应方向的长度稍大一些。
在一些实施例中,在计算各候选方向的最小无卷褶长度时,可以不考虑第一长度和第二长度之间的重叠区域,从而将公式[1]进行简化为公式[2],公式[2]如下。
公式[2]
在一些实施例中,若待匀场区域包含在成像对象的区域内,则公式[2]可以进一步简化为Lma=La1,Lmb=Lb1,Lmc=Lc1。
在一些实施例中,对于匀场扫描各方向设定视野时,为进一步确保无卷褶产生,可以适当放大视野,例如可以在原有基础上乘以一定的系数因子,系数因子可以是1.1、1.15、1.2或者1.25。
步骤S206,由所述匀场扫描获取包含所述待匀场区域的匀场图像。匀场图像是对待匀场区域进行磁共振扫描得到的图像。具体的,可以根据在磁场区域105内设定好的匀场扫描的视野,例如根据之前步骤计算得到方向a为相位编码方向,方向b为选层方向,方向c为读出方向,各方向对应的视野FOVPE=Lma,FOVSS=Lb2,FOVRO=Lb2,另根据其他一些扫描参数,由脉冲控制单元111控制对应的脉冲序列(包含梯度脉冲和射频脉冲),如图6所示,控制对应的射频信号(RF)、不同方向(SS方向、PE方向、RO方向)的梯度信号,实现对待匀场区域信号的激发与采集(在ADC区间采集),可以采集到对应dTE间隔时间的磁共振信号,如图7所示,之后将采集到的磁共振信号填充到K空间,将K空间进行傅里叶变换重建得到图像,由于在ADC区间采集得到两个磁共振信号,因此得到对应的两幅图像,将两幅图像进行相位做差再除以时间间隔dTE就可以得到匀场区域的B0场匀场图像。。
本发明还提供了一种磁共振系统匀场方法,如图8所示,包括步骤S801,由磁共振匀场图像获取方法获取磁共振匀场图像,具体获取方法可以参考之前实施例的描述。
步骤S802,根据磁共振匀场图像计算匀场参数,并基于匀场参数开展匀场。具体计算过程可以是和现有过程相同,根据磁共振匀场图像进行球函数的拟合展开,求得球函数的各阶系数,例如一阶项系数、二阶项系数等。其中,一阶项系数对应三个垂直方向的梯度,二阶项系数对应磁共振系数的二阶匀场线圈。磁共振系统可以将一阶项系数对应设置为梯度偏置,二阶项系数对应设置为对应二阶匀场线圈的偏置电流,从而实现匀场。
本发明的技术方案提供了一种磁共振成像匀场调节方法及装置,可以在较短的时间内采集较高分辨率的B0场,大大提高了B0匀场的准确性和稳定性。进一步的,特别是提升了小关节压脂的压脂图像效果及稳定性。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
同时,本申请使用了特定词语来描述本申请的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的某一特征、结构或特点。因此,应强调并注意的是,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一替代性实施例”并不一定是指同一实施例。此外,本申请的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外,本领域技术人员可以理解,本申请的各方面可以通过若干具有可专利性的种类或情况进行说明和描述,包括任何新的和有用的工序、机器、产品或物质的组合,或对他们的任何新的和有用的改进。相应地,本申请的各个方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“子模块”、“引擎”、“单元”、“子单元”、“组件”或“系统”。此外,本申请的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品,该产品包括计算机可读程序编码。
计算机可读信号介质可能包含一个内含有计算机程序编码的传播数据信号,例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式,包括电磁形式、光形式等等、或合适的组合形式。计算机可读信号介质可以是除计算机可读存储介质之外的任何计算机可读介质,该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机可读信号介质上的程序编码可以通过任何合适的介质进行传播,包括无线电、电缆、光纤电缆、射频信号、或类似介质、或任何上述介质的组合。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写,包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等,常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP,动态编程语言如Python、Ruby和Groovy,或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下,远程计算机可以通过任何网络形式与用户计算机连接,比如局域网(LAN)或广域网(WAN),或连接至外部计算机(例如通过因特网),或在云计算环境中,或作为服务使用如软件即服务(SaaS)。
此外,除非权利要求中明确说明,本申请所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本申请流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例,但应当理解的是,该类细节仅起到说明的目的,附加的权利要求并不仅限于披露的实施例,相反,权利要求旨在覆盖所有符合本申请实施例实质和范围的修正和等价组合。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在现有的服务器或移动设备上安装所描述的系统。
同理,应当注意的是,为了简化本申请披露的表述,从而帮助对一个或多个发明实施例的理解,前文对本申请实施例的描述中,有时会将多种特征归并至一个实施例、附图或对其的描述中。但是,这种披露方法并不意味着本申请对象所需要的特征比权利要求中提及的特征多。实际上,实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字,应当理解的是,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,说明书和权利要求中使用的数值参数均为近似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施例中用于确认其范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。
针对本申请引用的每个专利、专利申请、专利申请公开物和其他材料,如文章、书籍、说明书、出版物、文档等,特此将其全部内容并入本申请作为参考。与本申请内容不一致或产生冲突的申请历史文件除外,对本申请权利要求最广范围有限制的文件(当前或之后附加于本申请中的)也除外。需要说明的是,如果本申请附属材料中的描述、定义、和/或术语的使用与本申请所述内容有不一致或冲突的地方,以本申请的描述、定义和/或术语的使用为准。
最后,应当理解的是,本申请中所述实施例仅用以说明本申请实施例的原则。其他的变形也可能属于本申请的范围。因此,作为示例而非限制,本申请实施例的替代配置可视为与本申请的教导一致。相应地,本申请的实施例不仅限于本申请明确介绍和描述的实施例。
Claims (10)
1.一种磁共振匀场图像获取方法,其特征在于,包括:
获取成像对象的预扫描图像;基于待匀场区域确定候选方向集合,所述候选方向集合包括至少两个候选方向;
根据所述预扫描图像计算所述成像对象在各所述候选方向的第一长度;
计算待匀场区域在各所述候选方向的第二长度;
基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度,将各所述最小无卷褶长度中的最小值所对应的候选方向作为匀场扫描的相位编码方向,将所述相位编码方向对应的最小无卷褶长度作为所述匀场扫描相位编码方向的视野;
由所述匀场扫描获取包含所述待匀场区域的匀场图像。
2.根据权利要求1所述的磁共振匀场图像获取方法,其特征在于,所述基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度包括:
计算所述第一长度和第二长度的重叠区域;
基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度。
3.根据权利要求2所述的磁共振匀场图像获取方法,其特征在于,所述基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度包括:以所述第一长度和第二长度之和减去所述重叠区域作为最小无卷褶长度。
4.根据权利要求1所述的磁共振匀场图像获取方法,其特征在于,还包括:将所述待匀场区域在非所述相位编码方向的第二长度作为匀场扫描在非相位编码方向的视野。
5.根据权利要求1所述的磁共振匀场图像获取方法,其特征在于,所述成像对象为水模样品。
6.一种磁共振系统匀场方法,其特征在于,包括:
根据权利要求1-5任一项所述的磁共振匀场图像获取方法获取包含待匀场区域的匀场图像;
基于所述匀场图像计算匀场参数,根据所述匀场参数进行磁共振系统匀场。
7.根据权利要求6所述的磁共振系统匀场方法,其特征在于,所述基于所述匀场图像计算匀场参数包括:
基于所述匀场图像进行球函数展开,计算球函数的各阶系数;
基于所述各阶系数计算匀场参数。
8.根据权利要求7所述的磁共振系统匀场方法,其特征在于,所述匀场参数是关于梯度偏置和二阶电流的匀场参数。
9.一种磁共振系统,其特征在于,包括:
存储器,用于存储可由处理器执行的指令;
处理器,用于执行指令,所述指令包括:
获取样品的预扫描图像;
基于待匀场区域确定候选方向集合,所述候选方向集合包括至少两个候选方向;
根据所述预扫描图像计算所述样品在各所述候选方向的第一长度;
计算待匀场区域在各所述候选方向的第二长度;
基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度,将各所述最小无卷褶长度中的最小值所对应的候选方向作为所述匀场扫描的相位编码方向,将所述相位编码方向对应的最小无卷褶长度作为所述匀场扫描相位编码方向的视野;
由所述匀场扫描获取包含所述待匀场区域的匀场图像。
10.根据权利要求9所述的磁共振系统,其特征在于,所述基于所述第一长度和第二长度计算各所述候选方向的最小无卷褶长度包括:
计算所述第一长度和第二长度的重叠区域;
基于所述第一长度、第二长度和所述重叠区域计算最小无卷褶长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611250997.1A CN106680750B (zh) | 2016-12-29 | 2016-12-29 | 磁共振匀场图像获取方法、匀场方法及磁共振系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611250997.1A CN106680750B (zh) | 2016-12-29 | 2016-12-29 | 磁共振匀场图像获取方法、匀场方法及磁共振系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106680750A true CN106680750A (zh) | 2017-05-17 |
CN106680750B CN106680750B (zh) | 2019-07-19 |
Family
ID=58872688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611250997.1A Active CN106680750B (zh) | 2016-12-29 | 2016-12-29 | 磁共振匀场图像获取方法、匀场方法及磁共振系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106680750B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107358588A (zh) * | 2017-07-12 | 2017-11-17 | 上海联影医疗科技有限公司 | 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备 |
CN109709503A (zh) * | 2019-02-13 | 2019-05-03 | 上海联影医疗科技有限公司 | 一种磁共振系统控制方法、磁共振系统、计算机设备 |
CN109839608A (zh) * | 2019-03-15 | 2019-06-04 | 上海联影医疗科技有限公司 | 磁共振场图确定方法、装置、电子设备及存储介质 |
CN110074786A (zh) * | 2019-04-30 | 2019-08-02 | 上海东软医疗科技有限公司 | 核磁共振匀场方法、装置、计算设备及核磁共振成像系统 |
CN110415213A (zh) * | 2019-06-24 | 2019-11-05 | 上海联影医疗科技有限公司 | 磁场均匀性检测方法、装置、计算机设备和存储介质 |
CN113466765A (zh) * | 2020-03-31 | 2021-10-01 | 通用电气精准医疗有限责任公司 | 磁共振扫描方法及系统、计算机可读存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0585973A1 (en) * | 1986-10-29 | 1994-03-09 | Hitachi Medical Corporation | Method for correcting position deviation due to static magnetic field change in an NMR imaging device |
CN1435155A (zh) * | 2002-01-28 | 2003-08-13 | 西门子公司 | 通过自动匹配测量场的磁共振成像方法 |
US20060122485A1 (en) * | 2004-11-02 | 2006-06-08 | Oliver Heid | Optimized method for prevention of foldover artifacts in magnetic resonance tomography |
CN1977181A (zh) * | 2004-06-28 | 2007-06-06 | 皇家飞利浦电子股份有限公司 | 平行磁共振成像 |
US20070276221A1 (en) * | 2004-03-12 | 2007-11-29 | Koninklijke Philips Electronics N.V. | Prescan for optimization of mri scan parameters |
CN102768347A (zh) * | 2012-08-09 | 2012-11-07 | 厦门大学 | 核磁共振波谱仪上缩减相位编码数快速三维梯度匀场方法 |
-
2016
- 2016-12-29 CN CN201611250997.1A patent/CN106680750B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0585973A1 (en) * | 1986-10-29 | 1994-03-09 | Hitachi Medical Corporation | Method for correcting position deviation due to static magnetic field change in an NMR imaging device |
CN1435155A (zh) * | 2002-01-28 | 2003-08-13 | 西门子公司 | 通过自动匹配测量场的磁共振成像方法 |
US20070276221A1 (en) * | 2004-03-12 | 2007-11-29 | Koninklijke Philips Electronics N.V. | Prescan for optimization of mri scan parameters |
CN1977181A (zh) * | 2004-06-28 | 2007-06-06 | 皇家飞利浦电子股份有限公司 | 平行磁共振成像 |
US20060122485A1 (en) * | 2004-11-02 | 2006-06-08 | Oliver Heid | Optimized method for prevention of foldover artifacts in magnetic resonance tomography |
CN102768347A (zh) * | 2012-08-09 | 2012-11-07 | 厦门大学 | 核磁共振波谱仪上缩减相位编码数快速三维梯度匀场方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107358588A (zh) * | 2017-07-12 | 2017-11-17 | 上海联影医疗科技有限公司 | 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备 |
CN107358588B (zh) * | 2017-07-12 | 2020-07-17 | 上海联影医疗科技有限公司 | 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备 |
CN109709503A (zh) * | 2019-02-13 | 2019-05-03 | 上海联影医疗科技有限公司 | 一种磁共振系统控制方法、磁共振系统、计算机设备 |
CN109709503B (zh) * | 2019-02-13 | 2021-03-23 | 上海联影医疗科技股份有限公司 | 一种磁共振系统控制方法、磁共振系统、计算机设备 |
US11313927B2 (en) | 2019-02-13 | 2022-04-26 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for magnetic resonance imaging |
US11977137B2 (en) | 2019-02-13 | 2024-05-07 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for magnetic resonance imaging |
CN109839608A (zh) * | 2019-03-15 | 2019-06-04 | 上海联影医疗科技有限公司 | 磁共振场图确定方法、装置、电子设备及存储介质 |
CN109839608B (zh) * | 2019-03-15 | 2021-04-30 | 上海联影医疗科技股份有限公司 | 磁共振场图确定方法、装置、电子设备及存储介质 |
CN110074786A (zh) * | 2019-04-30 | 2019-08-02 | 上海东软医疗科技有限公司 | 核磁共振匀场方法、装置、计算设备及核磁共振成像系统 |
CN110074786B (zh) * | 2019-04-30 | 2022-12-06 | 上海东软医疗科技有限公司 | 核磁共振匀场方法、装置、计算设备及核磁共振成像系统 |
CN110415213A (zh) * | 2019-06-24 | 2019-11-05 | 上海联影医疗科技有限公司 | 磁场均匀性检测方法、装置、计算机设备和存储介质 |
CN113466765A (zh) * | 2020-03-31 | 2021-10-01 | 通用电气精准医疗有限责任公司 | 磁共振扫描方法及系统、计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106680750B (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106680750A (zh) | 磁共振匀场图像获取方法、匀场方法及磁共振系统 | |
CN104583799B (zh) | 基于快速图像采集的运动跟踪 | |
JP5558737B2 (ja) | 磁気共鳴イメージング装置 | |
JP6270709B2 (ja) | 磁気共鳴生データの再構成方法および装置 | |
CN109658468B (zh) | 磁共振参数成像方法、装置、设备及存储介质 | |
CN109310361A (zh) | 用于快速磁共振成像的系统和方法 | |
US11327132B2 (en) | System and method for magnetic resonance imaging acceleration | |
CN105074491A (zh) | 具有使用压缩感知的图像重建的动态mri | |
US10466330B2 (en) | Magnetic resonance diffusion weighted imaging method and apparatus | |
CN107607895A (zh) | 磁共振成像方法和系统 | |
JP6417406B2 (ja) | 強調磁化率コントラストによるmrイメージング | |
CN107367703A (zh) | 磁共振扫描方法、系统、装置及存储介质 | |
JP5678163B2 (ja) | 磁気共鳴イメージング装置 | |
CN106680748A (zh) | 磁共振图像获取方法及磁共振系统 | |
US11415656B2 (en) | Medical information processing apparatus, magnetic resonance imaging apparatus, and medical information processing method | |
JP6552804B2 (ja) | 磁気共鳴イメージング装置 | |
CN111012349A (zh) | 磁共振成像中信号表征确定的系统与方法 | |
EP4071494A1 (en) | Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset | |
CN110226100A (zh) | 用于磁共振成像的系统和方法 | |
CN107656225A (zh) | 磁共振频率校准方法、磁共振成像方法及系统 | |
Koppers et al. | Delimit pytorch-an extension for deep learning in diffusion imaging | |
JP2021512678A (ja) | 脂肪/水分離を使用したmri | |
US20210223343A1 (en) | Methods and systems for fingerprinting magnetic resonance imaging and machine learning | |
CN113050009B (zh) | 三维磁共振快速参数成像方法和装置 | |
CN206671522U (zh) | 磁共振局部线圈底座、局部线圈组件以及脚踝线圈组件 |
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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
CP01 | Change in the name or title of a patent holder |