一种图像重建方法、装置及设备与流程

文档序号:14949733发布日期:2018-07-17 22:16阅读:173来源:国知局

本申请涉及数据处理领域,具体涉及一种图像重建方法、装置及设备。



背景技术:

磁共振设备的特点是扫描时间长,而无论是在快速成像还是在提高成像质量方面,都需要加快磁共振设备的成像速度。

sense为业界较为流行的利用线圈通道接收器的空间分布特性来减少扫描时间的一种并行成像技术。具体的,在进行扫描时,通过减小相位编码步数,加快磁共振设备的成像速度,扫描出一个卷叠的图像,进而利用线圈通道接收器的空间分布特性对卷叠的图像进行去卷叠处理,重建出一个完整的图像。

其中,利用线圈通道接收器的空间分布特性对折叠的图像进行去卷叠处理,决定了重建出的图像质量,而线圈通道接收器的空间分布特性通常用敏感度谱数据表征,所以,如何得到准确的敏感度谱数据是图像重建的关键影响因素。



技术实现要素:

有鉴于此,本申请提供了一种图像重建方法、装置及设备,能够得到较准确的敏感度谱数据,进而得到质量较好的重建图像。

第一方面,本发明提供了一种图像重建方法,所述方法应用于磁共振设备,所述磁共振设备具有阵列线圈和大体线圈,利用所述磁共振设备对目标对象进行扫描,所述方法包括:

获取所述阵列线圈采集到的预设第一视野的第一线圈数据,并基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;

获取所述阵列线圈采集到的预设第二视野的第二线圈数据,以及所述大体线圈采集到的所述预设第二视野的第三线圈数据,并基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据;

根据所述第二敏感度谱数据,计算敏感度谱校正数据;

利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到校正后的敏感度谱数据;

利用所述校正后的敏感度谱数据进行图像重建。

可选的,所述基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据,包括:

计算所述第一线圈数据的绝对值的平方和,作为第一中间值;

将所述第一线圈数据与所述第一中间值的商,作为第一敏感度谱数据。

可选的,所述基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据,包括:

将所述第二线圈数据与所述第三线圈数据的商,作为第二敏感度谱数据。

可选的,所述根据所述第二敏感度谱数据,计算敏感度谱校正数据,包括:

计算所述第二敏感度谱数据的绝对值的平方和,作为第二中间值;

对所述第二中间值进行开平方处理,得到敏感度谱校正数据。

可选的,所述利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到较正后的敏感度谱数据,包括:

将所述第一敏感度谱数据与所述敏感度谱校正数据进行点乘处理,得到较正后的敏感度谱数据。

可选的,所述利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到较正后的敏感度谱数据之前,还包括:

对所述敏感度谱校正数据进行预处理,所述预处理包括多项式拟合、低通滤波和\或阈值处理。

第二方面,本发明提供了一种图像重建装置,所述装置应用于磁共振设备,所述磁共振设备具有阵列线圈和大体线圈,利用所述磁共振设备对目标对象进行扫描,所述装置包括:

第一敏感度谱数据获取单元,用于获取所述阵列线圈采集到的预设第一视野的第一线圈数据,并基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;

第二敏感度谱数据获取单元,用于获取所述阵列线圈采集到的预设第二视野的第二线圈数据,以及所述大体线圈采集到的所述预设第二视野的第三线圈数据,并基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据;

敏感度谱校正数据获取单元,用于根据所述第二敏感度谱数据,计算敏感度谱校正数据;

敏感度谱数据校正单元,用于利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到校正后的敏感度谱数据;

图像重建单元,用于利用所述校正后的敏感度谱数据进行图像重建。

可选的,所述第一敏感度谱数据获取单元包括:

第一敏感度谱数据计算子单元,用于计算所述第一线圈数据的绝对值的平方和,作为第一中间值;并将所述第一线圈数据与所述第一中间值的商,作为第一敏感度谱数据。

可选的,所述第二敏感度谱数据获取单元包括:

第二敏感度谱数据计算子单元,用于将所述第二线圈数据与所述第三线圈数据的商,作为第二敏感度谱数据。

可选的,所述敏感度谱校正数据获取单元包括:

敏感度谱校正数据计算子单元,用于计算所述第二敏感度谱数据的绝对值的平方和,作为第二中间值;并对所述第二中间值进行开平方处理,得到敏感度谱校正数据。

可选的,所述敏感度谱数据校正单元包括:

敏感度谱数据计算子单元,用于将所述第一敏感度谱数据与所述敏感度谱校正数据进行点乘处理,得到较正后的敏感度谱数据。

可选的,所述图像重建设备还包括:

敏感度谱校正数据预处理单元,用于对所述敏感度谱校正数据进行预处理,所述预处理包括多项式拟合、低通滤波和\或阈值处理。

第三方面,本发明提供了一种图像重建设备,所述设备包括:存储器和处理器,

所述存储器用于存储程序代码,并将所述程序代码传输给所述处理器;

所述处理器用于根据所述程序代码中的指令,执行以下步骤:

获取所述阵列线圈采集到的预设第一视野的第一线圈数据,并基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;

获取所述阵列线圈采集到的预设第二视野的第二线圈数据,以及所述大体线圈采集到的所述预设第二视野的第三线圈数据,并基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据;

根据所述第二敏感度谱数据,计算敏感度谱校正数据;

利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到校正后的敏感度谱数据;

利用所述校正后的敏感度谱数据进行图像重建。

本申请提供的图像重建方法获取阵列线圈采集到的第一线圈数据,根据该第一线圈数据计算得到第一敏感度谱数据,由于第一线圈数据为仅利用阵列线圈采集的数据,因此,若直接利用该第一敏感度谱数据对图像进行重建,得到的重建图像不存在伪影,但是基于阵列线圈采集到的数据均匀性较差的原因,该重建图像均匀性较差,中心区域较暗。为了得到显示效果更好的重建图像,本申请提供的图像重建方法进一步获取阵列线圈采集到的预设第二视野的第二线圈数据,和大体线圈采集到的预设第二视野的第三线圈数据,根据第二线圈数据和第三线圈数据计算得到第二敏感度谱数据,将该第二敏感度谱数据与第一敏感度谱数据进行结合,获得敏感度谱数据,由于第二敏感度谱数据是基于阵列线圈和大体线圈采集的数据计算得到的,大体线圈采集到的数据能够改善阵列线圈采集到的数据的均匀性,所以,结合第二敏感度谱数据能够改善重建图像的均匀性,因此,利用该敏感度谱数据进行图像重建,可以得到均匀性较好,且不存在伪影的重建图像。

附图说明

为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。

图1为采用sense技术进行图像扫描得到的扫描结果的示意图;

图2为采用并行成像技术对sense技术扫描结果进行图像重建的示意图;

图3为采用sense技术进行图像扫描,各线圈通道接收器接收的图像数据示意图;

图4为采用并行成像技术,进行图像重建的过程的示意图;

图5为本申请实施例提供的一种图像重建方法的流程图;

图6为本申请实施例提供的一种在具体应用场景中的图像重建方法的流程图;

图7为采用本申请实施例提供的图像重建方法得到的重建图像示意图;

图8为本申请实施例提供的一种图像重建装置的结构图;

图9为本申请实施例提供的一种图像重建设备的硬件结构示意图。

具体实施方式

下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。

由于磁共振设备具有扫描时间长的特点,当利用磁共振设备进行扫描时,会出现成像速度慢、成像质量低等问题,因此,为了解决因磁共振设备扫描时间长而出现的问题,业界内现普遍采用sense技术减少磁共振设备的扫描时间。

sense技术利用线圈通道接收器的空间分布特性,减少磁共振设备的扫描时间。具体的,在进行扫描时,通过减少相位编码步数,加快磁共振设备的成像速度,但是减少相位编码步数,会导致扫描得到的图像为一个卷叠的图像。因此,还需要根据线圈通道接收器的空间分布特性,对卷叠的图像进行去卷叠处理,以重建出一个完整的不存在卷叠的图像。

为了便于理解,现结合附图对sense技术进行介绍:

如图1所示,针对一幅磁共振图像,在不改变扫描分辨率的前提下,减少相位编码步数,磁共振图像会发生卷叠。如图2所示,利用图像空间的并行成像技术对有卷叠的磁共振图像进行重建,得到一幅完整的不存在卷叠的图像。

其中,利用图像空间的并行成像技术对有卷叠的磁共振图像进行重建,需要根据每个线圈通道接收器的空间分布特性,对有卷叠的磁共振图像进行去卷叠处理,以重建图像。

具体的,每个线圈通道接收器在空间中所处的位置均不相同,因此会导致每个线圈通道接收器所接收到的图像数据的灰度值不同,在各个线圈通道接收器接收到的图像数据中,越靠近线圈通道接收器,图像数据的灰度值越亮,相反,越远离线圈通道接收器,图像数据的灰度值越暗。由于每个线圈通道接收器在空间中所处的位置不同,因此,对于同一扫描对象,不同线圈通道接收器采集到的图像数据的灰度值不同。相应地,对于有卷叠的图像数据,可以根据每个线圈通道接收器接收到的卷叠图像数据灰度值的不同,求解出一幅完整的不存在卷叠的图像。

如图3所示,图3左侧为线圈通道接收器和被扫描人体的空间位置示意图,线圈通道接收器1位于被扫描人体的前方,线圈通道接收器2位于被扫描人体的后方。图3右侧分别为线圈通道接收器1接收到的图像数据和线圈通道接收器2接收到的图像数据,通过观察可以发现,由于线圈通道接收器1和线圈通道接收器2所处的位置不同,每个线圈通道接收器所接收到的图像数据的灰度值也不同,在线圈通道接收器1接收到的图像数据中,左侧的灰度较亮,右侧的灰度较暗,在线圈通道接收器2接收到的图像数据中,左侧的灰度较暗,右侧的灰度较亮。

如图4所示,减少相位编码步数后,图像会发生卷叠,但是由于每个线圈通道接收器在空间中所处的位置不同,每个线圈通道接收器接收到的卷叠图像数据的灰度值也不完全一样,可以基于各个线圈通道接收器接收到的卷叠图像数据的灰度值的不同,相应地建立对应于各个卷叠图像数据的敏感度谱,并利用空间分布的敏感度谱对图像进行去卷叠处理,以重建图像。

具体的,去卷叠处理重建图像的方法如下:

1)计算敏感度谱s。

2)计算卷叠点对应敏感度谱,计算公式如下:

sγ,ρ=sγ(rρ)公式(5)

其中,γ代表线圈通道,ρ代表卷叠点,rρ代表卷叠点位置。

3)针对每个卷叠点位置rρ,具有如下公式:

a=sγ,ρu公式(6)

其中,a为卷叠图像数据,如4图中线圈通道接收器1或线圈通道接收器2接收到的卷叠图像;u为重建图像数据。

采用正则化的最小二乘法计算重建图像数据u的公式如下:

其中,ψ-1为噪声矩阵,r为正则化因子,h用于表示共轭转置,

为对sγ,ρ进行共轭转置处理。

为了获得显示效果更好的磁共振扫描图像,本申请提供的图像重建方法采用了一种新的敏感度谱计算方法,基于该敏感度谱计算方法,对卷叠图像进行去卷叠处理,重建扫描图像。该图像重建方法应用于磁共振设备,该磁共振设备具有阵列线圈和大体线圈,利用该磁共振设备对目标对象进行扫描。

具体的,获取阵列线圈采集到的预设第一视野的第一线圈数据,根据该第一线圈数据计算第一敏感度谱数据;获取阵列线圈采集到的预设第二视野的第二线圈数据,以及大体线圈采集到的预设第二视野的第三线圈数据,根据第二线圈数据和第三线圈数据,计算第二敏感度谱数据;根据该第二敏感度谱数据,计算敏感度谱校正数据;再利用敏感度谱校正数据,校正第一敏感度谱数据得到敏感度谱数据,进而,利用该敏感度谱数据进行图像重建。

本申请提供的图像重建方法获取阵列线圈采集到的第一线圈数据,根据该第一线圈数据计算得到第一敏感度谱数据,由于第一线圈数据为仅利用阵列线圈采集的数据,因此,若直接利用该第一敏感度谱数据对图像进行重建,得到的重建图像虽然不存在伪影,但是基于阵列线圈采集到的数据均匀性较差的原因,该重建图像均匀性较差,中心区域较暗。为了得到显示效果更好的重建图像,本申请提供的图像重建方法进一步获取阵列线圈采集到的预设第二视野的第二线圈数据,和大体线圈采集到的预设第二视野的第三线圈数据,根据第二线圈数据和第三线圈数据计算得到第二敏感度谱数据,将该第二敏感度谱数据与第一敏感度谱数据进行结合,获得敏感度谱数据,由于第二敏感度谱数据是基于阵列线圈和大体线圈采集的数据计算得到的,大体线圈采集到的数据能够改善阵列线圈采集到的数据的均匀性,所以,结合第二敏感度谱数据能够改善重建图像的均匀性,因此,利用该敏感度谱数据进行图像重建,可以得到均匀性较好,且不存在伪影的重建图像。

方法实施例一

参见图5,为本实施例提供的一种图像重建方法的流程图,该图像重建方法应用于磁共振设备,该磁共振设备具有阵列线圈和大体线圈,利用该磁共振设备对目标对象进行扫描时,该图像重建方法包括:

步骤501:获取阵列线圈采集到的预设第一视野的第一线圈数据,并基于第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据。

实际应用中,在利用磁共振设备对待扫描部位进行扫描时,首先确定预设第一视野,并在利用该磁共振设备中的阵列线圈采集待扫描部位在该预设第一视野的第一线圈数据后,获取该第一线圈数据,用于敏感度谱数据的计算。

获取阵列线圈采集到的预设第一视野的第一线圈数据,将基于该第一线圈数据计算得到的敏感度谱数据,作为第一敏感度谱数据,具体的,如基于第一线圈数据计算第一敏感度谱数据的公式(1)所示,计算第一线圈数据的绝对值的平方和,作为第一中间值;将第一线圈数据与第一中间值的商,作为第一敏感度谱数据:

其中,slγ表示第一敏感度谱数据,cγ表示第一线圈数据。

需要说明的是,阵列线圈中包含有γ个线圈通道接收器,第一线圈数据cγ代表阵列线圈中γ个线圈通道接收器采集得到的线圈数据,因此,cγ为对应于阵列线圈中γ个线圈通道接收器所采集的线圈数据的γ个矩阵。相应的,经公式(1)计算得到的第一敏感度谱数据slγ,代表阵列线圈中γ个线圈通道接收器的敏感度谱数据,因此,slγ为对应于阵列线圈中γ个线圈通道接收器敏感度谱数据的γ个矩阵。

此外,还可以采用其他方法,根据第一线圈数据计算第一敏感度谱数据,在此不对第一敏感度谱数据的获取方法进行任何限定。

步骤502:获取阵列线圈采集到的预设第二视野的第二线圈数据,以及大体线圈采集到的预设第二视野的第三线圈数据,并基于第二线圈数据和第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据。

在实际应用中,利用磁共振设备对待扫描部位进行扫描时,预先确定预设第二视野,利用磁共振设备中的阵列线圈对处于预设第二视野中的待扫描部位进行扫描,以便采集得到第二线圈数据,然后由阵列线圈切换至大体线圈,再利用磁共振设备中的大体线圈对处于预设第二视野中的待扫描部位进行扫描,以便采集得到第三线圈数据。本发明实施例基于采集到的第二线圈数据和第三线圈数据,计算第二敏感度谱数据。

获取阵列线圈采集到的预设第二视野的第二线圈数据,和大体线圈采集到的预设第二视野的第三线圈数据,根据第二线圈数据和第三线圈数据,计算第二敏感度谱数据时,具体可以将第二线圈数据与第三线圈数据的商,作为第二敏感度谱数据,公式如下:

其中,sgγ表示第二敏感度谱数据,c′γ表示第二线圈数据,q表示第三线圈数据。

需要说明的是,阵列线圈中包含有γ个线圈通道接收器,第二线圈数据c′γ代表阵列线圈中γ个线圈通道接收器扫描处于预置大视野中的待扫描部位时,采集到的γ个线圈数据,因此,c′γ为对应于阵列线圈中所有线圈通道接收器所采集的线圈数据的γ个矩阵。而大体线圈中仅包含一个线圈通道接收器,因此,第三线圈数据q代表大体线圈中的线圈通道接收器扫描处于预置大视野中的待扫描部位时,采集到的线圈数据,因此,第三线圈数据q为对应于该大体线圈中的线圈通道接收器采集到的线圈数据的矩阵。

相应地,经公式(2)计算得到的第二敏感度谱数据sgγ,仍代表阵列线圈中γ个线圈通道接收器的敏感度谱数据,因此,sgγ为对应于阵列线圈中所有线圈通道接收器敏感度谱数据的γ个矩阵。

此外,还可以采用其他方法,根据第二线圈数据和第三线圈数据,计算第二敏感度谱数据,在此不对第二敏感度谱数据的获取方法做任何限定。

需要说明的是,步骤501和步骤502为两个并列的执行步骤,步骤501和步骤502的执行顺序不分先后,可以先执行步骤501,后执行步骤502,也可以先执行步骤502,后执行步骤501。

步骤503:根据第二敏感度谱数据,计算敏感度谱校正数据。

由于在步骤501中获取到的第一敏感度谱数据,为仅基于阵列线圈采集到的第一线圈数据计算得到的敏感度谱数据,利用该第一敏感度谱数据对图像进行重建,得到的图像虽然不存在伪影,但是图像的均匀性较差。而在步骤502中,将阵列线圈扫描得到的第二线圈数据和大体线圈扫描得到的第三线圈数据进行结合,经过计算得到了第二敏感度谱数据,由于大体线圈采集到的数据能够改善阵列线圈采集到的数据的均匀性,因此,将第二敏感度谱数据与第一敏感度谱数据进行结合,能够得到显示效果更好的重建图像。

为了将步骤502中获取到的第二敏感度谱数据,与步骤501中获取到的第一敏感度谱数据结合,现需要根据第二敏感度谱数据,计算敏感度谱校正数据。

具体的,计算敏感度谱校正数据时,可以先计算第二敏感度谱数据的绝对值的平方和,作为第二中间值;对第二中间值进行开平方处理,得到敏感度谱校正数据,如下公式:

其中,t表示敏感度谱校正数据。

具体的,计算各个采集通道γ的第二敏感度谱数据sgγ矩阵的模的平方,将各个第二敏感度谱数据sgγ矩阵的模的平方进行求和计算后,再开平方,得到用于校正第一敏感度谱数据slγ的敏感度谱校正数据t。

需要说明的是,根据公式(3)计算得到的敏感度谱校正数据t为一个用于校正的矩阵。

当然,还可以采用其他计算方法,根据第二敏感度谱数据计算敏感度谱校正数据,在此对计算敏感度谱校正数据的方法不做任何限定。

步骤504:利用敏感度谱校正数据,对第一敏感度谱数据进行校正,得到较正后的敏感度谱数据。

利用在步骤503中得到的敏感度谱校正数据t,对在步骤501中获取到的第一敏感度谱数据slγ进行校正,得到校正后的敏感度谱数据。

具体的,可以将第一敏感度谱数据与敏感度谱校正数据进行点乘处理,得到较正后的敏感度谱数据,如下公式:

sγ=slγ·t公式(4)

其中,sγ表示校正后的敏感度谱数据,slγ表示第一敏感度谱数据。

需要说明的是,公式(4)中,slγ与t是进行点乘计算,即矩阵的对应位置的元素相乘,由于第一敏感度谱数据slγ为对应于阵列线圈中所有线圈通道接收器的γ个矩阵,相应地,校正后的敏感度谱数据sγ也为对应于阵列线圈中所有线圈通道接收器的γ个矩阵。

由于敏感度谱校正数据t是根据第二敏感度谱数据sgγ计算获得的,因此,利用在步骤503中获得的敏感度谱校正数据t,对第一敏感度谱数据slγ进行校正,即实现了第二敏感度谱数据sgγ与第一敏感度谱数据slγ的结合,相应地,利用结合后得到的敏感度谱数据sγ对图像进行重建,得到的重建图像不存在伪影,且均匀性较好。

当然,也可以采用其他计算方法,对第一敏感度谱数据进行校正,在此不对第一敏感度谱数据的校正方法做任何限定。

步骤505:利用校正后的敏感度谱数据进行图像重建。

得到敏感度谱数据sγ后,可利用该敏感度谱数据sγ进行图像重建。

具体的,可以先计算与每个采集通道所采集的数据图像上的卷叠点对应的敏感度谱,一般采用如下公式:

sy,ρ=sγ(rρ)公式(5)

其中,γ代表采集通道,ρ代表采集图像上的卷叠点,rρ代表卷叠点在图像上所处的位置。

对于每个卷叠点所处的位置,均满足如下公式:

a=sγ,ρu公式(6)

其中,a代表卷叠图像数据,u代表需要重建图像数据,该重建图像数据中不存在卷叠。

采用正则化的最小二乘法计算重建图像数据,采用如下公式:

其中,ψ-1为噪声矩阵,r为正则化因子,h用于表示共轭转置,为对sγ,ρ进行共轭转置处理。

需要说明的是,在步骤504执行之前,即在利用敏感度谱校正数据,对第一敏感度谱数据进行校正,得到校正后的敏感度谱数据之前,还可以先对敏感度谱校正数据进行预处理,具体的,可以对敏感度谱校正数据进行多项式拟合、低通滤波和\或阈值处理,以保证得到的敏感度谱校正数据表现平滑。

本实施例提供的图像重建方法获取阵列线圈采集到的第一线圈数据,根据该第一线圈数据计算得到第一敏感度谱数据,由于第一线圈数据为仅利用阵列线圈采集的数据,因此,若直接利用该第一敏感度谱数据对图像进行重建,得到的重建图像不存在伪影,但是基于阵列线圈采集到的数据均匀性较差的原因,该重建图像均匀性较差,中心区域较暗。为了得到显示效果更好的重建图像,本实施例提供的图像重建方法进一步获取阵列线圈采集到的预设第二视野的第二线圈数据,和大体线圈采集到的预设第二视野的第三线圈数据,根据第二线圈数据和第三线圈数据计算得到第二敏感度谱数据,将该第二敏感度谱数据与第一敏感度谱数据进行结合,获得敏感度谱数据,由于第二敏感度谱数据是基于阵列线圈和大体线圈采集的数据计算得到的,大体线圈采集到的数据能够改善阵列线圈采集到的数据的均匀性,所以,结合第二敏感度谱数据能够改善重建图像的均匀性,因此,利用该敏感度谱数据进行图像重建,可以得到均匀性较好,且不存在伪影的重建图像。

下面结合实际应用,对上述方法实施例一中的图像重建方法进行介绍:

方法实施例二

参见图6,为本实施例提供的图像重建方法的流程图。

步骤601:利用磁共振设备对待扫描部位进行扫描后,阵列线圈采集到预设第一视野内的第一线圈数据。

在实际应用中,利用磁共振设备对待扫描部位进行扫描时,可以将待扫描部位所处的范围作为预设第一视野,利用磁共振设备对处于预设第一视野内的待扫描部位进行扫描,利用阵列线圈采集到预设第一视野内的第一线圈数据。

步骤602:在不改变扫描分辨率的前提下,减少相位编码步数,利用磁共振设备对待扫描部位进行再次扫描后,阵列线圈采集到图像卷叠数据。

由于磁共振设备的扫描时间过长,导致磁共振设备的成像质量较低,因此,为了提高磁共振设备的成像质量,实际应用中,在不改变扫描分辨率的前提下,通过减少相位编码步数,提高磁共振设备的扫描速度,进而提高磁共振设备的成像质量。

具体的,为了提高最终重建后图像的质量,需要尽量保证步骤602和步骤601中阵列线圈采集到的数据所属部位的一致性,例如,假设步骤601中阵列线圈采集到的是人体头部中左脑的数据,则步骤602中需要尽量保证阵列线圈采集到的也是左脑的数据,因此需要尽量在步骤602执行结束后立即执行步骤601,以避免因时间拖长造成人体头部不可避免的移动导致的数据所属部位不一致的问题,例如,避免上述步骤602中阵列线圈采集到的不是左脑的数据,而是因人体头部的移动而采集到的右脑的数据。

因此,实际应用中,在步骤s601之后,磁共振设备尽量立即减少相位编码步数,对预设第一视野中的待扫描部位再次进行扫描,阵列线圈采集得到图像卷叠数据。由于减少相位编码步数后得到的图像数据中存在有卷叠,因此,还需要对该图像卷叠数据进行去卷叠处理,得到显示效果较好的重建图像数据。

步骤603:利用磁共振设备对待扫描部位进行再次扫描后,阵列线圈采集到预设第二视野内的第二线圈数据,以及大体线圈采集到第二视野内的第三线圈数据。

在步骤602之后,磁共振设备将视野设置为第二视野,并再次对待扫描部位进行正常扫描,阵列线圈采集在第二视野内的第二线圈数据后,为了保证第二线圈数据和第三线圈数据的一致性,避免伪影问题,磁共振设备需要立即从阵列线圈切换至大体线圈,并利用大体线圈采集第二视野内的第三线圈数据。

步骤604:基于第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;以及基于第二线圈数据和第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据。

步骤605:根据第二敏感度谱数据,计算敏感度谱校正数据,并利用敏感度谱校正数据对第一敏感度谱数据进行校正,得到较正后的敏感度谱数据。

步骤606:利用较正后的敏感度谱数据,对图像卷叠数据进行重建。

由于第一敏感度谱数据为,仅根据阵列线圈扫描得到第一线圈数据计算得到的敏感度谱数据,因此,直接利用该第一敏感度谱数据对上述图像卷叠数据进行去卷叠处理,得到的重建图像虽然不存在伪影,但是该重建图像的均匀性较差。因此,为了获得显示效果更好的重建图像,还需要对第一敏感度谱数据进行校正。

具体的,利用该第二线圈数据和第三线圈数据进行计算,得到第二敏感度谱数据。由于大体线圈采集的数据具有较好的均匀性,因此,根据该第二敏感度谱数据,计算敏感度谱校正数据,利用该敏感度谱校正数据对第一敏感度谱数据进行校正,利用校正后的敏感度谱数据对图像进行重建,能够得到不存在伪影,且均匀性较好的重建图像。

值得注意的是,上述方法实施例二只是本发明中图像重建方法的一种具体实现方式,不作为本发明具体实现方式的限定,实际应用中还存在其他实现方式,在此不再介绍。

需要说明的是,上述方法实施例一、二中,预设第一视野和预设第二视野的大小关系不做限定,具体的,预设第二视野可以大于预设第一视野,也可以等于预设第一视野,还可以小于预设第一视野。

一种情况下,当预设第二视野大于预设第一视野时,基于阵列线圈在预设第二视野采集到的第二线圈数据,以及大体线圈在预设第二视野采集到的第三线圈数据,计算得到的第二敏感度谱数据,不仅可以用于阵列线圈在第一视野内采集到的图像卷叠数据的重建,还可以用于阵列线圈采集到的其他部位(在第二视野内)的图像卷叠数据的重建,节省对后续其他部位的图像卷叠数据的重建步骤,提高效率。

例如,要利用该磁共振设备扫描人体头部中的左脑,将左脑具体所处的视野范围设置为预设第一视野,将人体头部所处的视野范围设置为预设第二视野,分别利用阵列线圈和大体线圈采集人体头部所处视野范围内的第二线圈数据和第三线圈数据,基于该第二线圈数据和第三线圈数据计算得到第二敏感度谱数据,在后续扫描该人体头部中的其他部位(如右脑)时,可以直接使用该第二敏感度谱数据,用于对其他部位的图像卷叠数据的重建,无需再次对人体头部所处视野范围内的第二线圈数据和第三线圈数据进行采集。

另一种情况下,当预设第二视野等于预设第一视野时,步骤602之后不需要对视野进行调整,直接利用磁共振设备对待扫描部位进行再次扫描,利用阵列线圈采集预设第一视野(即第二视野)得到第二线圈数据,以及利用大体线圈采集该预设第一视野得到第三线圈数据,并计算得到的第二敏感度谱数据。

再一种情况下,当预设第二视野小于预设第一视野时,基于步骤603中阵列线圈和大体线圈采集到的第二线圈数据和第三线圈数据计算得到的第二敏感度谱数据,只是与步骤601中计算得到第一敏感度谱数据的一部分相对应,为了获得能够与第一敏感度谱数据全部对应的第二敏感度谱数据,本发明实施例可以基于敏感度谱数据的平滑特性以及待扫描部位图像本身的特性,利用外插处理的方法,对第二敏感度谱数据进行补偿计算,以得到较为准确的能够与第一敏感度谱数据全部对应的第二敏感度谱数据。

综上所述,利用方法实施例一的图像重建方法在实际应用中进行图像重建,能够获得显示效果较好的重建图像,重建得到的图像不存在伪影,且均匀性较好。

如图7所示,为采用本实施例提供的图像重建方法,对图像进行重建得到的显示效果,通过观察可以发现,采用本实施例提供的图像重建方法重建得到的图像中不存在伪影,且图像的均匀性较好。

装置实施例

参见图8,为本实施例提供的一种图像重建装置的结构图,该装置应用于磁共振设备,所述磁共振设备具有阵列线圈和大体线圈,利用所述磁共振设备对目标对象进行扫描,所述装置包括:

第一敏感度谱数据获取单元801,用于获取所述阵列线圈采集到的预设第一视野的第一线圈数据,并基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;

第二敏感度谱数据获取单元802,用于获取所述阵列线圈采集到的预设第二视野的第二线圈数据,以及所述大体线圈采集到的所述预设第二视野的第三线圈数据,并基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据;

敏感度谱校正数据获取单元803,用于根据所述第二敏感度谱数据,计算敏感度谱校正数据;

敏感度谱数据校正单元804,用于利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到较正后的敏感度谱数据;

图像重建单元805,用于利用所述校正后的敏感度谱数据进行图像重建。

可选的,所述第一敏感度谱数据获取单元包括:

第一敏感度谱数据计算子单元,用于计算所述第一线圈数据的绝对值的平方和,作为第一中间值;并将所述第一线圈数据与所述第一中间值的商,作为第一敏感度谱数据;

可以利用公式(1)计算敏感度谱数据,作为第一敏感度谱数据;

其中,slγ表示所述第一敏感度谱数据,cγ表示所述第一线圈数据。

可选的,所述第二敏感度谱数据获取单元包括:

第二敏感度谱数据计算子单元,用于将所述第二线圈数据与所述第三线圈数据的商,作为第二敏感度谱数据;

可以利用公式(2)计算敏感度谱数据,作为第二敏感度谱数据;

其中,sgγ表示所述第二敏感度谱数据,c′γ表示所述第二线圈数据,q表示所述第三线圈数据。

可选的,所述敏感度谱校正数据获取单元包括:

敏感度谱校正数据计算子单元,用于计算所述第二敏感度谱数据的绝对值的平方和,作为第二中间值;并对所述第二中间值进行开平方处理,得到敏感度谱校正数据;

可以利用公式(3)计算敏感度谱校正数据;

其中,t表示所述敏感度谱校正数据。

可选的,所述敏感度谱数据校正单元包括:

敏感度谱数据计算子单元,用于将所述第一敏感度谱数据与所述敏感度谱校正数据进行点乘处理,得到较正后的敏感度谱数据;

可以利用公式(4)计算较正后的敏感度谱数据;

sγ=slγ·t公式(4)

其中,sγ表示较正后的敏感度谱数据,slγ表示所述第一敏感度谱数据。

可选的,所述图像重建设备还包括:

敏感度谱校正数据预处理单元,用于对所述敏感度谱校正数据进行预处理,所述预处理包括多项式拟合、低通滤波和\或阈值处理。

本实施例图8所示的装置是与方法实施例所述的方法对应的装置,具体实现方法类似,参考方法实施例中的描述,这里不再赘述。

相应的,本发明实施例还提供一种图像重建设备,参见图9所示,可以包括:

处理器901、存储器902、输入装置903和输出装置904。处理器801的数量可以一个或多个,图9中以一个处理器为例。在本发明的一些实施例中,处理器901、存储器902、输入装置903和输出装置904可通过总线或其它方式连接,其中,图9中以通过总线连接为例。

存储器902可用于存储软件程序以及模块,处理器901通过运行存储在存储器902的软件程序以及模块,从而执行光射野一致性检测设备的各种功能应用以及数据处理。存储器902可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序等。此外,存储器902可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。输入装置903可用于接收输入的数字或字符信息,以及产生与光射野一致性检测设备的用户设置以及功能控制有关的信号输入。

具体在本实施例中,处理器901会按照如下的指令,将一个或一个以上的应用程序的进程对应的可执行文件加载到存储器902中,并由处理器901来运行存储在存储器902中的应用程序,从而实现各种功能:

获取所述阵列线圈采集到的预设第一fov的第一线圈数据,并基于所述第一线圈数据计算敏感度谱数据,作为第一敏感度谱数据;

获取所述阵列线圈采集到预设第二fov的第二线圈数据,以及所述大体线圈采集到的所述预设第二fov的第三线圈数据,并基于所述第二线圈数据和所述第三线圈数据计算敏感度谱数据,作为第二敏感度谱数据;

根据所述第二敏感度谱数据,计算敏感度谱校正数据;

利用所述敏感度谱校正数据,对所述第一敏感度谱数据进行校正,得到校正后的敏感度谱数据;

利用所述校正后的敏感度谱数据进行图像重建。

对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。

需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

以上对本申请实施例所提供的一种图像重建方法、装置及设备进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1