CN109633503A - 磁共振图像重建方法、装置、计算机设备和存储介质 - Google Patents
磁共振图像重建方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN109633503A CN109633503A CN201811623082.XA CN201811623082A CN109633503A CN 109633503 A CN109633503 A CN 109633503A CN 201811623082 A CN201811623082 A CN 201811623082A CN 109633503 A CN109633503 A CN 109633503A
- Authority
- CN
- China
- Prior art keywords
- interphase
- phase
- data line
- interphase phase
- space data
- 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5619—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请涉及一种磁共振图像重建方法、装置、计算机设备和存储介质。方法包括:获取磁共振信号,根据磁共振信号得到K空间数据;将K空间数据划分到一个或多个时间期相内,计算K空间数据各时间期相的时间分辨率以及各时间期相的数据线数目;根据各时间期相的时间分辨率以及各时间期相的数据线数目,调整K空间数据的时间期相;并进行图像重建,得到医学图像。通过各时间期相的时间分辨率对K空间数据的时间期相进行调整,能够针对造影剂进入器官的时机适当的调整时间期相。通过各时间期相的数据线数目,能够尽可能的减小由于呼吸门控技术导致的连续采样条件被破坏的情况。
Description
技术领域
本申请涉及磁共振图像重建技术领域,特别是涉及一种磁共振图像重建方法、装置、计算机设备和存储介质。
背景技术
MRI动态灌注成像是一种重要的非介入医学检查手段。经过一段时间的连续扫描(例如:连续扫描40秒),产生多个期相的图(例如:四幅图,每10秒钟一幅图)。此项检查展示造影剂随血液进入器官的过程,为疾病诊断提供不可替代的信息。其中,期相是指在时间上将数据分在多幅图像中,每幅图像的时间就是期相。
目前的传统技术,都是将全部扫描时间等分为多个期相,并重建得到图像。但在利用MRI动态灌注成像进行检查时,造影剂灌注过程有两个特点:1)在造影剂没有进入器官之前和造影剂完全通过器官之后的图像对比度随时间变化缓慢,造影剂逐渐通过器官的过程图像对比度随时间变化迅速。因此传统技术无法针对需要适当的调整时间期相。2)由于MRI动态灌注成像需要进行连续扫描,在连续扫描的过程中,如果连续时间过长,待检测患者无法进行屏气。因此呼吸运动将会引入图像模糊和图像伪影,降低图像质量。为了消除呼吸运动的影响,需要使用呼吸门控技术对采集的数据进行分类。但使用呼吸门控技术导致连续采样条件被破坏。
发明内容
基于此,有必要针对上述技术问题,提供一种能够适当的调整时间分辨率并且减小连续采样条件被破坏的磁共振图像重建方法、装置、计算机设备和存储介质。
一种磁共振图像重建方法,所述方法包括:获取磁共振信号,并根据所述磁共振信号得到K空间数据;将所述K空间数据划分到一个或多个时间期相内,并计算各时间期相的时间分辨率以及各时间期相的数据线数目;根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;根据调整后的K空间数据的时间期相以及各时间期相的数据线进行图像重建,得到医学图像。
在其中一个实施例中,所述磁共振信号为采用黄金角径向采集模式获得。
在其中一个实施例中,所述计算各时间期相的时间分辨率包括:获取各时间期相第一条数据线的采集时间以及最后一条数据线的采集时间;将所述第一条数据线的采集时间以及最后一条数据线的采集时间做差,得到各时间期相的时间分辨率。
在其中一个实施例中,所述根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相包括:根据各时间期相的所述时间分辨率,设置第一约束条件;根据各时间期相的所述数据线数目,设置第二约束条件;根据所述第一约束条件以及第二约束条件,确定每个时间期相包含的K空间数据。
在其中一个实施例中,所述第一约束条件为:所述时间分辨率小于等于预设时间分辨率阈值。
在其中一个实施例中,所述第二约束条件为:所述数据线数目大于等于预设数据线数目阈值。
在其中一个实施例中,所述方法还包括在将所述K空间数据划分到一个或多个时间期相内之前对所述K空间数据进行筛选,所述筛选包括:
根据所述K空间数据,确定运动曲线;
根据所述运动曲线筛选出满足重建条件的K空间数据。
在其中一个实施例中,所述方法还包括:根据所述K空间数据,计算K空间数据线对应的亮度值;根据K空间数据线的亮度值,设置第三约束条件;根据所述第一约束条件、第二约束条件以及第三约束条件,确定每个时间期相包含的K空间数据。
一种磁共振图像重建装置,所述装置包括:获取模块,用于获取磁共振信号,并根据所述磁共振信号得到K空间数据;计算模块,用于将所述K空间数据划分到一个或多个时间期相内,并计算各时间期相的时间分辨率以及各时间期相的数据线数目;调整模块,用于根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;重建模块,用于根据调整后的K空间数据的时间期相以及个时间期相的数据线进行重建,得到医学图像。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述任一种所述方法的步骤。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述任一种所述的方法的步骤。
上述磁共振图像重建方法、装置、计算机设备和存储介质,通过获取磁共振信号,得到K空间数据;将所述K空间数据划分到一个或多个时间期相内,计算各时间期相的时间分辨率以及各时间期相的数据线数目,根据各时间期相的时间分辨率以及各时间期相的数据线数目,调整K空间数据的时间期相,最后根据调整后的K空间数据的时间期相以及各时间期相的数据线进行重建,得到医学图像。通过各时间期相的时间分辨率对K空间数据的时间期相进行调整,能够针对造影剂进入器官的时机适当的调整时间期相。通过各时间期相的数据线数目,能够尽可能的减小由于呼吸门控技术导致的连续采样条件被破坏的情况。
附图说明
图1为一个实施例中磁共振图像重建方法的流程示意图;
图2为一个实施例中磁共振图像重建装置的结构框图;
图3为一个实施例中计算机设备的内部结构图;
图4为采用现有技术重建得到的不同时间期相的磁共振图像;
图5为采用如图1所述的方法重建得到的不同时间期相的磁共振图像。
附图标记:获取模块100、计算模块200、调整模块300、重建模块400。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
核磁共振(magnetic resonance,MR)检查,是将人体置于特殊的磁场中,用无线电射频脉冲激发人体内氢原子核,引起氢原子核共振,并吸收能量。在停止射频脉冲后,氢原子核按特定频率发出射电信号,并将吸收的能量释放出来,被体外的接受器收录,经电子计算机处理获得图像,这就叫做核磁共振成像。核磁共振检查彻底摆脱了电离辐射对人体的损害,又有参数多,信息量大,可多方位成像,以及对软组织有高分辨力等突出的特点,被广泛用于临床疾病的诊断,对有些病变成为必不可少的检查方法。
磁共振灌注成像是指用来反映组织的微血管灌注分布及血流灌注情况的磁共振检查技术。目前依据其成像原理可大致分为二种类型,即对比剂首过灌注成像以及动脉血质子自旋标记技术。其中,对比剂首过灌注成像属于动态增强磁共振成像范畴。动态增强磁共振成像技术注入顺磁性造影剂如二乙三胺五醋酸钆(Gd-DTPA)等。顺磁性的造影剂一进入毛细血管床便在毛细血管内外建立起多个小的局部磁场,即形成一定的磁敏感性差别,类似于在毛细管与组织间建立了无数小梯度磁场,这样不仅使组织质子所经历的磁场均匀性降低,而且导致质子相位相干的损失,即加速了质子的失相位过程,从而使组织的T1、T2时间均缩短。这时使用T2敏感序列进行测量,即可观察到组织信号的显著减小,即所谓的负性增强;如果用对T1时间敏感的序列检查,则表现为组织的正性增强。
传统的MRI动态灌注成像通常采用笛卡尔坐标采集模式。其中,笛卡尔坐标采集模式是指在固定读出编码梯度幅值、方向和整数倍相位编码梯度面积的情况下,在时间上均匀地采集数据,所得数据在K空间上对应的坐标为整数值。而在本申请的实施例中使用径向采集模式,其相对于传统的笛卡尔坐标采集模式可以重建多种相位分辨率的动态图。在连续采集的条件下,每个期相的时间分辨率正比于包括的K空间线数目。例如,每条线的脉冲重复时间TR是5ms,一个期相包括N=200条K空间线,时间分辨率就是N*TR=1000ms。一个期相的时间越长,说明时间分辨率越低;一个期相的时间越短,说明时间分辨率越高。也就是说,如果减少一个期相包含的K空间线的数目,时间分辨率就会升高;如果增加一个期相包含的K空间线的数目,时间分辨率就会降低。而径向采集模式使得N等于任意值的时候径向采集线分布都比较均匀,有利于重建。
提高时间分辨率需要减少每个期相内包含的K空间线的数目。但是减少K空间线的数目会降低重建之后的图像质量。当一个期相内包含的k空间线的数目下降到一个阈值以下,重建之后的图像的已经无法满足诊断需要。使用呼吸门控技术之后,一个期相的时间分辨率和其对应的k空间线的数目之间已经没有简单的正比关系。所以,仅仅考虑造影剂的影响提高时间分辨率可能造成某些期相包含过少的k空间线,从而造成重建后图像的质量难以满足临床诊断需要。
在一个实施例中,如图1所示,提供了一种磁共振图像重建方法,包括以下步骤:
步骤S102,获取磁共振信号,并根据磁共振信号得到K空间数据。
具体地,磁共振信号可以是通过磁共振设备利用基于黄金角(golden angle)径向采集模式实时采集的磁共振信号,也可以是存储介质中存储的利用黄金角(golden angle)径向采集模式采集的磁共振信号。对采集到的磁共振信号进行相位编码,得到多条数据线,再将多条数据线填充至K空间,得到K空间数据。其中,K空间也就是傅里叶变换空间,是带有空间定位编码信息的磁共振信号的填充空间。每一幅磁共振图像都有其相应的K空间数据。对K空间数据进行傅里叶转换,就能对原始磁共振信号中的空间定位编码信息进行解码,得到磁共振的图像数据,即把不同信号强度的磁共振信息分配到相应的空间位置上,即可重建出磁共振图像。
步骤S104,将参与重建的所述K空间数据划分到一个或多个时间期相内,并计算各时间期相的时间分辨率以及各时间期相的数据线数目。
具体地,时间期相为生成一幅磁共振图像所需要的时间。一个时间期相的时间越长,说明时间分辨率越低;一个时间期相的时间越短,说明时间分辨率越高。在通过K空间数据计算各时间期相的时间分辨率时,首先获取每个时间期相第一条数据线的采集时间,再获取每个时间期相最后一条数据线的采集时间,将第一条数据线的采集时间和最后一条数据线的采集时间做差,得到每个时间期相的时间分辨率。也即,利用最后一条数据线的采集时间减去第一条数据线的采集时间,得到时间期相的时间分辨率。各时间期相的数据线数目为各时间期相内参与重建的K空间数据线数目。
在其中一个实施例中,所述方法还包括在将所述K空间数据划分到一个或多个时间期相内之前对所述K空间数据进行筛选,所述筛选包括:根据所述K空间数据,计算运动曲线,根据所述运动曲线确定对应K空间数据是否参与重建。
具体地,根据所述K空间数据,确定运动曲线;根据所述运动曲线筛选出满足重建条件的K空间数据。运动曲线为通过呼吸门控技术得到的曲线。运动曲线表示随时间的变化人体组织的运动情况。根据运动曲线确定对应K空间数据是否参与重建,也就是将人体组织运动幅度大的时间段所对应的K空间数据不参与重建。从而得到参与重建的K空间数据。
步骤S106,根据各时间期相的时间分辨率以及各时间期相的数据线数目,调整K空间数据的时间期相。
首先根据所有时间期相的时间分辨率设置第一约束条件,以及根据所有时间期相的数据线数目设置第二约束条件,再根据第一约束条件和第二约束条件,确定每个时间期相包含的K空间数据线数。其中,第一约束条件为:时间分辨率小于等于预设时间分辨率阈值;第二约束条件为:数据线数目大于等于预设数据线数目阈值。例如:每条数据线的耗时为0.2s,将所有时间期相中最高的时间分辨率作为时间分辨率阈值,当最高时间分辨率阈值要求每幅图像的成像时间不得长于10s,也就是说,每个时间期相的数据线数目小于等于50条数据线。根据图像质量的要求根据每个时间期相的数据线数目,设置数据线数目阈值,当数据线数目阈值设置为20条数据线时,也就是说,每个时间期相的数据线数目大于等于20条数据线。结合两个条件,也就是时间期相需要满足,其数据线数目大于等于20条数据线,小于等于50条数据线。确定每个时间期相包含的K空间数据线数,从而达到调整K空间数据的时间期相的目的。
在其中一个实施例中,还可以加入亮度值的约束条件。根据K空间数据,计算K空间数据线对应的亮度值,也就是计算K空间数据各时间期相的亮度值与前一时间期相亮度值的差值绝对值,再计算得到的差值绝对值与前一时间期相亮度值的比值,根据个时间期相的比值,设置第三约束条件,根据第一约束条件、第二约束条件以及第三约束条件,确定每个时间期相包含的K空间数据线数。其中,第三约束条件还可以是每条K空间数据线的亮度值的函数。根据每条K空间数据线的亮度值的函数确定造影剂进入器官的时间,然后再确定所需要的时间分辨率。
步骤S108,根据调整后的K空间数据的时间期相以及各时间期相的数据线进行图像重建,得到医学图像。
上述磁共振图像重建方法通过综合考虑对比度随时间的变化率以及使用呼吸门控技术后的采样率,确定自适应的时间分辨率以及每个时间期相包括的采样线数目,进行自适应可变时间分辨率的重建。也就是自动在造影剂进入组织时适当增加时间分辨率,同时在呼吸运动强烈的时候适当降低时间分辨率,可以在保证图像质量的情况下提高时间分辨率。
上述磁共振图像重建方法,通过获取磁共振信号,得到K空间数据,根据K空间数据,计算各时间期相的时间分辨率以及各时间期相的数据线数目,根据各时间期相的时间分辨率以及各时间期相的数据线数目,调整K空间数据的时间期相,最后根据调整后的K空间数据的时间期相以及各时间期相的数据线进行重建,得到医学图像。通过各时间期相的时间分辨率对K空间数据的时间期相进行调整,能够针对造影剂进入器官的时机适当的调整时间期相。通过各时间期相的数据线数目,能够尽可能的减小由于呼吸门控技术导致的连续采样条件被破坏的情况。
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中,如图2所示,提供了一种磁共振图像重建装置,包括:获取模块100、计算模块200、调整模块300和重建模块400,其中:
获取模块100,用于获取磁共振信号,并根据所述磁共振信号得到K空间数据。
计算模块200,用于将所述K空间数据划分到一个或多个时间期相内,并计算K空间数据各时间期相的时间分辨率以及各时间期相的数据线数目。
调整模块300,用于根据各时间期相的时间分辨率以及各时间期相的数据线数目,调整K空间数据的时间期相。
重建模块400,用于根据调整后的K空间数据的时间期相以及个时间期相的数据线进行重建,得到医学图像。
获取模块100包括:磁共振信号获取单元、编码单元以及填充单元。
磁共振信号获取单元,用于获取磁共振信号。
编码单元,用于对磁共振信号进行相位编码,得到多条数据线。
填充单元,用于将多条数据线填充至K空间,得到K空间数据。
计算模块200包括:采集时间获取单元以及计算单元。
采集时间获取单元,用于获取各时间期相第一条数据线的采集时间以及最后一条数据线的采集时间。
计算单元,用于将第一条数据线的采集时间以及最后一条数据线的采集时间做差,得到各时间期相的时间分辨率。
调整模块300包括:约束条件设置单元以及调整单元。
约束条件设置单元,用于根据各时间期相的时间分辨率,设置第一约束条件;根据各时间期相的数据线数目,设置第二约束条件。
调整单元,用于根据第一约束条件以及第二约束条件,确定每个时间期相包含的K空间数据。
约束条件设置单元,还用于根据所述K空间数据,计算K空间数据线对应的亮度值;根据K空间数据线的亮度值,设置第三约束条件。
调整单元,还用于根据所述第一约束条件、第二约束条件以及第三约束条件,确定每个时间期相包含的K空间数据。
计算模块200,还用于根据所述K空间数据,确定运动曲线;根据所述运动曲线筛选出满足重建条件的K空间数据。
关于磁共振图像重建装置的具体限定可以参见上文中对于磁共振图像重建方法的限定,在此不再赘述。上述磁共振图像重建装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端,其内部结构图可以如图3所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种磁共振图像重建方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图3中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
获取磁共振信号,并根据所述磁共振信号得到K空间数据;将所述K空间数据划分到一个或多个时间期相内,并计算K空间数据各时间期相的时间分辨率以及各时间期相的数据线数目;根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;根据调整后的K空间数据的时间期相以及各时间期相的数据线进行图像重建,得到医学图像。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
获取磁共振信号;对磁共振信号进行相位编码,得到多条数据线;将多条数据线填充至K空间,得到K空间数据。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
获取各时间期相第一条数据线的采集时间以及最后一条数据线的采集时间;将第一条数据线的采集时间以及最后一条数据线的采集时间做差,得到各时间期相的时间分辨率。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
根据各时间期相的时间分辨率,设置第一约束条件;根据各时间期相的数据线数目,设置第二约束条件;根据第一约束条件以及第二约束条件,确定每个时间期相包含的K空间数据。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
根据所述K空间数据,确定运动曲线;根据所述运动曲线筛选出满足重建条件的K空间数据。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
根据K空间数据,计算K空间数据各时间期相的亮度值与前一时间期相亮度值的差值绝对值;根据差值绝对值以及前一时间期相的亮度值,计算差值绝对值与前一时间期相亮度值的比值;根据各时间期相的比值,设置第三约束条件;其中,第三约束条件为:比值小于等于预设比值阈值;根据第三约束条件,确定每个时间期相包含的K空间数据。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
获取磁共振信号,并根据所述磁共振信号得到K空间数据;将所述K空间数据划分到一个或多个时间期相内,并计算K空间数据各时间期相的时间分辨率以及各时间期相的数据线数目;根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;根据调整后的K空间数据的时间期相以及各时间期相的数据线进行图像重建,得到医学图像。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
获取磁共振信号;对磁共振信号进行相位编码,得到多条数据线;将多条数据线填充至K空间,得到K空间数据。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
获取各时间期相第一条数据线的采集时间以及最后一条数据线的采集时间;将第一条数据线的采集时间以及最后一条数据线的采集时间做差,得到各时间期相的时间分辨率。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
根据各时间期相的时间分辨率,设置第一约束条件;根据各时间期相的数据线数目,设置第二约束条件;根据第一约束条件以及第二约束条件,确定每个时间期相包含的K空间数据。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
根据所述K空间数据,确定运动曲线;根据所述运动曲线筛选出满足重建条件的K空间数据。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
根据K空间数据,计算K空间数据各时间期相的亮度值与前一时间期相亮度值的差值绝对值;根据差值绝对值以及前一时间期相的亮度值,计算差值绝对值与前一时间期相亮度值的比值;根据各时间期相的比值,设置第三约束条件;其中,第三约束条件为:比值小于等于预设比值阈值;根据第三约束条件,确定每个时间期相包含的K空间数据。
图4为采用现有技术重建得到的不同时间期相的磁共振图像,扫描部位为腹部,多幅图像之间采用等时间间隔分辨率,从左至右依次为第一时间期相至第八时间期相采集得到的磁共振信号重建得到的磁共振图像。
图5为采用如图1所述的方法重建得到的不同时间期相的磁共振图像,扫描部位为腹部,多幅图像之间采用自适应时间间隔分辨率,从左至右依次为第一时间期相至第八时间期相采集得到的磁共振信号重建得到的磁共振图像。将图5与图4对比可知,采用本申请自适应重建方法,对比不同时间期相可清晰得到造影剂的渐变过程,对造影剂进入肝部的过程有更加精细的表现,同时兼顾了图像本身的质量。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种磁共振图像重建方法,其特征在于,所述方法包括:
获取磁共振信号,并根据所述磁共振信号得到K空间数据;
将所述K空间数据划分到一个或多个时间期相内,并计算各时间期相的时间分辨率以及各时间期相的数据线数目;
根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;
根据调整后的K空间数据的时间期相以及各时间期相的数据线进行图像重建,得到医学图像。
2.根据权利要求1所述的方法,其特征在于,所述磁共振信号为采用黄金角径向采集模式获得。
3.根据权利要求2所述的方法,其特征在于,所述计算各时间期相的时间分辨率包括:
获取各时间期相第一条数据线的采集时间以及最后一条数据线的采集时间;
将所述第一条数据线的采集时间以及最后一条数据线的采集时间做差,得到各时间期相的时间分辨率。
4.根据权利要求3所述的方法,其特征在于,所述根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相包括:
根据各时间期相的所述时间分辨率,设置第一约束条件;
根据各时间期相的所述数据线数目,设置第二约束条件;
根据所述第一约束条件以及第二约束条件,确定每个时间期相包含的K空间数据。
5.根据权利要求4所述的方法,其特征在于,
所述第一约束条件为:所述时间分辨率小于等于预设时间分辨率阈值;
所述第二约束条件为:所述数据线数目大于等于预设数据线数目阈值。
6.根据权利要求1所述的方法,其特征在于,所述方法还包括在将所述K空间数据划分到一个或多个时间期相内之前对所述K空间数据进行筛选,所述筛选包括:
根据所述K空间数据,确定运动曲线;
根据所述运动曲线筛选出满足重建条件的K空间数据。
7.根据权利要求4所述的方法,其特征在于,所述方法还包括:
根据所述K空间数据,计算K空间数据线对应的亮度值;
根据K空间数据线的亮度值,设置第三约束条件;
根据所述第一约束条件、第二约束条件以及第三约束条件,确定每个时间期相包含的K空间数据。
8.一种磁共振图像重建装置,其特征在于,所述装置包括:
获取模块,用于获取磁共振信号,并根据所述磁共振信号得到K空间数据;
计算模块,用于将所述K空间数据划分到一个或多个时间期相内,并计算各时间期相的时间分辨率以及各时间期相的数据线数目;
调整模块,用于根据各时间期相的所述时间分辨率以及各时间期相的所述数据线数目,调整K空间数据的时间期相;
重建模块,用于根据调整后的K空间数据的时间期相以及个时间期相的数据线进行重建,得到医学图像。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811623082.XA CN109633503B (zh) | 2018-12-28 | 2018-12-28 | 磁共振图像重建方法、装置、计算机设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811623082.XA CN109633503B (zh) | 2018-12-28 | 2018-12-28 | 磁共振图像重建方法、装置、计算机设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633503A true CN109633503A (zh) | 2019-04-16 |
CN109633503B CN109633503B (zh) | 2021-01-26 |
Family
ID=66078808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811623082.XA Active CN109633503B (zh) | 2018-12-28 | 2018-12-28 | 磁共振图像重建方法、装置、计算机设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633503B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146835A (zh) * | 2019-05-22 | 2019-08-20 | 山东颐邦齐鲁医生集团管理有限公司 | 一种基于并行成像的自导航磁共振图像重建方法及装置 |
CN110176294A (zh) * | 2019-05-31 | 2019-08-27 | 数坤(北京)网络科技有限公司 | 一种血管cta影像数据的调度方法、装置及可读存储介质 |
CN111009020A (zh) * | 2019-12-02 | 2020-04-14 | 上海联影医疗科技有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
CN113100741A (zh) * | 2020-01-13 | 2021-07-13 | 上海联影医疗科技股份有限公司 | 一种磁共振扫描方法、设备及存储介质 |
CN115267632A (zh) * | 2022-05-12 | 2022-11-01 | 上海东软医疗科技有限公司 | 磁共振灌注成像方法、装置、计算机设备及可读存储介质 |
WO2023040918A1 (en) * | 2021-09-14 | 2023-03-23 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for vascular analysis |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540115A (zh) * | 2011-12-12 | 2012-07-04 | 中国科学院深圳先进技术研究院 | 磁共振成像方法及系统 |
CN103099618A (zh) * | 2011-11-15 | 2013-05-15 | 西门子公司 | 用于生成磁共振血管造影图像的方法和相应的磁共振设备 |
CN103384836A (zh) * | 2010-12-22 | 2013-11-06 | 皇家飞利浦电子股份有限公司 | 包含grappa算子的用于任意k空间轨迹的快速两步并行重建 |
CN103901379A (zh) * | 2012-12-28 | 2014-07-02 | 上海联影医疗科技有限公司 | 磁共振扫描成像方法及系统 |
CN106597333A (zh) * | 2016-12-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 一种磁共振并行成像方法及磁共振成像系统 |
US20170248671A1 (en) * | 2012-04-04 | 2017-08-31 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus and method |
CN109001660A (zh) * | 2018-06-12 | 2018-12-14 | 上海联影医疗科技有限公司 | 电影成像方法及磁共振成像系统 |
-
2018
- 2018-12-28 CN CN201811623082.XA patent/CN109633503B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103384836A (zh) * | 2010-12-22 | 2013-11-06 | 皇家飞利浦电子股份有限公司 | 包含grappa算子的用于任意k空间轨迹的快速两步并行重建 |
CN103099618A (zh) * | 2011-11-15 | 2013-05-15 | 西门子公司 | 用于生成磁共振血管造影图像的方法和相应的磁共振设备 |
CN102540115A (zh) * | 2011-12-12 | 2012-07-04 | 中国科学院深圳先进技术研究院 | 磁共振成像方法及系统 |
US20170248671A1 (en) * | 2012-04-04 | 2017-08-31 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus and method |
CN103901379A (zh) * | 2012-12-28 | 2014-07-02 | 上海联影医疗科技有限公司 | 磁共振扫描成像方法及系统 |
CN106597333A (zh) * | 2016-12-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 一种磁共振并行成像方法及磁共振成像系统 |
CN109001660A (zh) * | 2018-06-12 | 2018-12-14 | 上海联影医疗科技有限公司 | 电影成像方法及磁共振成像系统 |
Non-Patent Citations (2)
Title |
---|
STEPHEN R. YUTZY等: "Improvements in Multislice Parallel Imaging Using Radial CAIPIRINHA", 《MAGNETIC RESONANCE IN MEDICINE》 * |
黄小倩等: "基于同时多层激发和并行成像的心脏磁共振电影成像", 《波普学杂志》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146835A (zh) * | 2019-05-22 | 2019-08-20 | 山东颐邦齐鲁医生集团管理有限公司 | 一种基于并行成像的自导航磁共振图像重建方法及装置 |
CN110146835B (zh) * | 2019-05-22 | 2021-09-07 | 山东颐邦齐鲁医生集团管理有限公司 | 一种基于并行成像的自导航磁共振图像重建方法及装置 |
CN110176294A (zh) * | 2019-05-31 | 2019-08-27 | 数坤(北京)网络科技有限公司 | 一种血管cta影像数据的调度方法、装置及可读存储介质 |
CN111009020A (zh) * | 2019-12-02 | 2020-04-14 | 上海联影医疗科技有限公司 | 图像重建方法、装置、计算机设备和存储介质 |
CN113100741A (zh) * | 2020-01-13 | 2021-07-13 | 上海联影医疗科技股份有限公司 | 一种磁共振扫描方法、设备及存储介质 |
CN113100741B (zh) * | 2020-01-13 | 2024-02-27 | 上海联影医疗科技股份有限公司 | 一种磁共振扫描方法、设备及存储介质 |
WO2023040918A1 (en) * | 2021-09-14 | 2023-03-23 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for vascular analysis |
CN115267632A (zh) * | 2022-05-12 | 2022-11-01 | 上海东软医疗科技有限公司 | 磁共振灌注成像方法、装置、计算机设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109633503B (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633503A (zh) | 磁共振图像重建方法、装置、计算机设备和存储介质 | |
Lin et al. | Artificial intelligence for MR image reconstruction: an overview for clinicians | |
US10561337B2 (en) | Rapid 3D dynamic arterial spin labeling with a sparse model-based image reconstruction | |
Schloegl et al. | Infimal convolution of total generalized variation functionals for dynamic MRI | |
Do et al. | Reconstruction of multicontrast MR images through deep learning | |
Haji‐Valizadeh et al. | Validation of highly accelerated real‐time cardiac cine MRI with radial k‐space sampling and compressed sensing in patients at 1.5 T and 3T | |
Roy et al. | Accelerated MRI of the fetal heart using compressed sensing and metric optimized gating | |
Ramzi et al. | NC-PDNet: A density-compensated unrolled network for 2D and 3D non-Cartesian MRI reconstruction | |
CN109917315B (zh) | 磁共振成像扫描方法、装置、计算机设备和存储介质 | |
CN109658471A (zh) | 一种医学图像重建方法和系统 | |
Montesinos et al. | Application of the compressed sensing technique to self‐gated cardiac cine sequences in small animals | |
US11047942B2 (en) | Magnetic resonance imaging method and system with optimal variable flip angles | |
Zhao et al. | Rapid 3D dynamic arterial spin labeling with a sparse model-based image reconstruction | |
US11835612B2 (en) | System and method for motion correction of magnetic resonance image | |
Naresh et al. | Accelerated, first‐pass cardiac perfusion pulse sequence with radial k‐space sampling, compressed sensing, and k‐space weighted image contrast reconstruction tailored for visual analysis and quantification of myocardial blood flow | |
Bell et al. | Simultaneous MRI of lung structure and perfusion in a single breathhold | |
Chang et al. | Deep learning-based rigid motion correction for magnetic resonance imaging: a survey | |
EP4066205A1 (en) | Systems and methods for medical imaging | |
Lee et al. | Rapid three‐dimensional steady‐state chemical exchange saturation transfer magnetic resonance imaging | |
Wang et al. | High‐fidelity direct contrast synthesis from magnetic resonance fingerprinting | |
Ljungberg et al. | Motion corrected silent ZTE neuroimaging | |
Keerthivasan et al. | Abdominal T2‐weighted imaging and T2 mapping using a variable flip angle radial turbo spin‐echo technique | |
Wilferth et al. | Motion‐corrected 23Na MRI of the human brain using interleaved 1H 3D navigator images | |
Eichhorn et al. | Deep Learning-Based Detection of Motion-Affected k-Space Lines for T2*-Weighted MRI | |
Schneider et al. | Non‐rigid motion‐compensated 3D whole‐heart T2 mapping in a hybrid 3T PET‐MR system |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |