CN109001660A - 电影成像方法及磁共振成像*** - Google Patents
电影成像方法及磁共振成像*** Download PDFInfo
- Publication number
- CN109001660A CN109001660A CN201810601448.7A CN201810601448A CN109001660A CN 109001660 A CN109001660 A CN 109001660A CN 201810601448 A CN201810601448 A CN 201810601448A CN 109001660 A CN109001660 A CN 109001660A
- Authority
- CN
- China
- Prior art keywords
- phase
- data
- magnetic resonance
- specified
- resonance scanning
- 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
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/023—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (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
【技术领域】
本方案涉及磁共振成像技术领域,尤其涉及一种电影成像方法及磁共振成像***。
【背景技术】
径向MRI(Magnetic Resonance Imaging,磁共振成像)采集技术相比于传统的笛卡尔采集方式具有更好的对于运动的不敏感度性等优点。当相邻两个径向采样的角度以黄金角(111.25°)增加时,即为黄金角径向采样。黄金角径向采样具有非常特殊的数学性质,即在连续的采集当中截取任意一段数据时,其在k空间的轨迹都会在角度方向形成相对均匀的分布。由于黄金角径向采样具有的这个优点,可以允许在重建时选取任意的时间分辨率,而非像传统的技术需要在序列开始时即选择好期相数等。
然而,径向采集技术用于心脏信号触发的、回顾性的心脏电影成像时,连续采集的径向数据会被根据心电信号分配到不同的期相中。此时,同一期相的数据来自于连续的采集数据中的某些分离片段的集合,而非一段连续的片段,这使得k空间的角度方向形成极端不均匀的分布,造成径向数据在角度方向的“聚集”。这种“聚集”导致重建图像中出现严重的放射状伪影,因此降低了图像质量。
【发明内容】
有鉴于此,本方案实施例提供了一种电影成像方法及磁共振成像***,用以解决现有技术中心脏电影成像时放射状伪影严重、图像质量较低的问题。
第一方面,本发明实施例提供一种电影成像方法,所述方法包括:
径向采集心脏的磁共振扫描数据,并记录所述磁共振扫描数据对应的心跳信号;
基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中;
对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;根据所述待重建数据进行图像重建,得到该期相的重建图像。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:
根据所述心跳信号,确定所述磁共振扫描数据对应的各个心跳周期;
将所述各个心跳周期中的每个心跳周期平均分成指定期相数目个期相;
对于每一个期相,从所述磁共振扫描数据中提取采集时间处于该期相的数据,作为该期相对应的第一分配数据。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:
根据所述心跳信号,确定所述磁共振扫描数据对应的各个心跳周期;
将所述各个心跳周期中的每个心跳周期归一化到标准周期上;
将所述磁共振扫描数据的采集时间归一化到所述标准周期上,得到归一化采集时间;
将所述标准周期平均分成指定期相数目个期相;
对于每一个期相,从所述磁共振扫描数据中提取归一化采集时间处于该期相的数据,作为该期相对应的第一分配数据。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔,包括:
在该期相k空间中查找大于角度间隔阈值的角度间隔,作为指定角度间隔;
从该期相的相邻期相对应的第二分配数据中,查找与所述指定角度间隔匹配的数据,作为填充数据;
用所述填充数据填充所述指定角度间隔。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,从该期相的相邻期相对应的第二分配数据中,查找与所述指定角度间隔匹配的数据,作为填充数据,包括:
获取所述第二分配数据对应的数据线,作为预选数据线;
从所述预选数据线中,查找角度处于所述指定角度间隔内的数据线,作为匹配数据线,所述匹配数据线上的数据,即为与所述指定角度间隔匹配的数据;
获取所述匹配数据线上的数据,作为填充数据。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述指定条件为:该期相k空间中的所有角度间隔均小于角度间隔阈值;或者,
所述指定条件为:在该期相的相邻期相对应的第二分配数据中,不存在与所述指定角度间隔匹配的数据。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,所述角度间隔阈值等于该期相k空间中的所有角度间隔的平均角度间隔与指定倍数的乘积。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,在满足指定条件时,获取该期相填充后的数据,作为待重建数据之前,所述方法还包括:
判断是否满足所述指定条件。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中之前,所述方法还包括:
根据所述心跳信号,查找处于心跳周期范围外的心跳周期,作为待删除心跳周期;
从所述磁共振扫描数据中,删除所述待删除心跳周期对应的数据,得到合格磁共振扫描数据;
基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:基于指定期相数目和所述心跳信号,将所述合格磁共振扫描数据分配到对应的期相中。
如上所述的方面和任一可能的实现方式,进一步提供一种实现方式,径向采集心脏的磁共振扫描数据,包括:黄金角径向采集心脏的磁共振扫描数据。
第二方面,本发明实施例提供一种磁共振成像***,所述***包括处理器以及存储器;所述存储器用于存储指令,所述指令被所述处理器执行时,导致所述***实现如第一方面任一项所述的方法。
本发明实施例具有以下有益效果:
本发明实施例,通过径向采集心脏的磁共振扫描数据,并记录磁共振扫描数据对应的心跳信号,基于指定期相数目和心跳信号,将磁共振扫描数据分配到对应的期相中,对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;根据待重建数据进行图像重建,得到该期相的重建图像,利用期相的相邻期相的数据对该期相k空间中的较大角度间隔进行填充,使得该期相k空间的角度方向分布趋于均匀,从而减弱了重建图像中的放射状伪影,因此提高了图像质量。
【附图说明】
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的电影成像方法的第一流程示例图。
图2为本发明实施例提供的心跳信号的示例图。
图3为本发明实施例提供的电影成像方法的第二流程示例图。
图4为根据本发明实施例提供的电影成像方法得到的心脏短轴图像与根据现有技术得到的心脏短轴图像的对比示例图。
图5为根据本发明实施例提供的电影成像方法得到的剔除心律不齐数据时的心脏图像与根据现有技术得到的剔除心律不齐数据时的心脏图像的对比示例图。
【具体实施方式】
为了更好的理解本发明的技术方案,下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。
应当理解,本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。
取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地,取决于语境,短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。
实施例一
本发明实施例提供了一种电影成像方法,该电影成像方法可以应用于心脏的磁共振扫描回顾性电影成像过程中。
图1为本发明实施例提供的电影成像方法的第一流程示例图。如图1所示,本实施例中,电影成像方法可以包括如下步骤:
S101,径向采集心脏的磁共振扫描数据,并记录磁共振扫描数据对应的心跳信号。
S102,基于指定期相数目和心跳信号,将磁共振扫描数据分配到对应的期相中。
S103,对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;根据待重建数据进行图像重建,得到该期相的重建图像。
根据上述方法可获得每一个期相心脏对应的磁共振图像,按照期相的顺序可将心脏图像动态显示。
通过步骤S103,利用期相的相邻期相的数据对该期相k空间中的较大角度间隔进行填充后,使得各个期相k空间的角度方向分布趋于均匀,从而减弱了重建图像中的放射状伪影,因此能够提高重建图像的图像质量。
其中,磁共振扫描数据是通过磁共振扫描的数据线记录到的。每一条数据线都具有一个时间戳(timestamp),该时间戳所显示的时间点为数据线上数据的采集时间。
其中,心跳信号(也即ECG(Electrocardiograph,心电图)信号)一方面可以提供磁共振扫描期间各个心跳周期的时间信息,另一方面可以作为磁共振扫描数据所对应的采集时间的参照信息,以便将磁共振扫描数据按照采集时间划分到各个心跳周期中,并进一步将同一个心跳周期对应的磁共振扫描数据分配到该心跳周期的各个期相中。
可选地,心跳周期/心动周期(cardiac cycle)指从一次心跳的起始到下一次心跳的起始,心血管***所经历的过程。心动周期各时相心室内压、心室容积、血流与瓣膜活动的变化如以心室的舒缩活动为中心,整个心动周期可按8个时相进行活动:等容收缩期、快速射血期、减慢射血期、舒张前期、等容舒张期、快速充盈期、减慢充盈期以及心房收缩期。
图2为本发明实施例提供的心跳信号的示例图。请参见图2,心跳信号可通过心电监控设备(electrocardiogram,ECG)获取,心跳信号的波形是周期性的,波形中的最高峰是R波,相邻两个R波之间的波形对应一个心跳周期/心动周期。
在实际采集的数据中,各个心跳周期的时间长度有可能是不一致的。比如,第一个心跳周期长度为873ms(毫秒,下同),第二个心跳周期长度为771ms,……。
在通过心跳信号获得各个心跳周期的时间长度后,还可以根据指定的正常心跳周期的时间长度范围,将心律不齐等异常心跳周期的磁共振扫描数据剔除。例如,时间长度小于正常心跳周期的时间长度范围的最小值的心跳周期,或者时间长度大于正常心跳周期的时间长度范围的最大值的心跳周期等。
其中,指定期相数目可以由用户选择或输入。
在一个示例性的实现过程中,基于指定期相数目和心跳信号,将磁共振扫描数据分配到对应的期相中,可以包括:根据心跳信号,确定磁共振扫描数据对应的各个心跳周期;将各个心跳周期中的每个心跳周期平均分成指定期相数目的期相;对于每一个期相,从磁共振扫描数据中提取采集时间处于该期相的数据,作为该期相对应的第一分配数据。
举例说明。假设根据心跳信号确定的第一个心跳周期长度为873ms,第二个心跳周期长度为771ms。心电信号R波到达后的第一条数据线的时间戳为3ms,第二条数据线的时间戳为6ms,第三条数据线的时间戳为9ms,…,第二百九十一条数据线的时间戳为873ms。此时新的R波到达,此后的第一条数据线的时间戳为3ms,第二条数据线的时间戳为6ms,第三条数据线的时间戳为9ms,…,第二百五十七条条数据线的时间戳为771ms。此时新的R波到达,此后的第一条数据线timestamp为3ms,……。
假设用户选择的期相数目为20,那么第一个心跳周期中每个期相的长度为873ms/20=43.65ms,这样,第一个心跳周期中第一个期相对应时间段0ms~43.65ms,第二个期相对应时间段43.65ms~87.3ms,……,第二十个期相对应829.35ms~873ms。第一个心跳周期中第一条数据线的时间戳为3ms,3ms处于0ms~43.65ms这个时间段,因此第一个心跳周期中第一条数据线的数据属于第一个心跳周期的第一期相。第一个心跳周期中第二条数据线的时间戳为6ms,6ms也处于0ms~43.65ms这个时间段,因此第一个心跳周期中第二条数据线的数据属于第一个心跳周期的第一期相,……,以此类推,可以将全部磁共振扫描数据分配到各个心跳周期的各个期相中,然后将各个心跳周期的同一期相的数据提取出来,即可得到每一期相对应的第一分配数据。
在一个示例性的实现过程中,基于指定期相数目和心跳信号,将磁共振扫描数据分配到对应的期相中,可以包括:根据心跳信号,确定磁共振扫描数据对应的各个心跳周期;将各个心跳周期中的每个心跳周期归一化到标准周期上;将磁共振扫描数据的采集时间归一化到标准周期上,得到归一化采集时间;将标准周期平均分成指定期相数目的期相;对于每一个期相,从磁共振扫描数据中提取归一化采集时间处于该期相的数据,作为该期相对应的第一分配数据。
举例说明。仍然假设根据心跳信号确定的第一个心跳周期长度为873ms,第一个心跳周期中,心电信号R波到达后的第一条数据线的时间戳为3ms,第二条数据线的时间戳为6ms,第三条数据线的时间戳为9ms,…,第二百九十一条数据线的时间戳为873ms。将第一个心跳周期归一化到标准周期上,假设标准周期为1000ms。
首先,计算第一个心跳周期的归一化系数=1000ms/873ms,该归一化系数与每一个数据线的timestamp的乘积,即为该数据线上数据的归一化采集时间。例如,第一个心跳周期中第一条数据线上数据的归一化采集时间为3*1000/873=3.43ms,第二条数据线上数据的归一化采集时间为6*1000/873=6.87ms,…,第二百九十一条数据线上数据的归一化采集时间为873*1000/873=1000ms。
假设用户选择的指定期相数目为20,那么每个期相的长度为1000ms/20=50ms。也就是说第一个心跳周期第一期相对应0ms~50ms,第二期相对应50ms~100ms,…,第二十期相对应950ms~1000ms。
对于第一个心跳周期归一化之后的数据来说,第一条数据线上数据的归一化采集时间为3.43ms,处于0ms~50ms这个时间段,因此属于第一期相。第二条数据线上数据的归一化采集时间为6.87ms,也处于0ms~50ms这个时间段,因此属于第一期相,……,以此类推,第二百九十一条数据线上数据的归一化采集时间线为1000ms,处于950ms~1000ms这个时间段,属于第二十期相。
对于其他心跳周期的数据,也按照第一个心跳周期的数据的处理方式,分配到各自心跳周期的相应期相中。然后将各个心跳周期的同一期相的数据提取出来,即可得到每一期相对应的第一分配数据。
在一个示例性的实现过程中,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔,可以包括:在该期相k空间中查找大于角度间隔阈值的角度间隔,作为指定角度间隔;从该期相的相邻期相对应的第二分配数据中,查找与指定角度间隔匹配的数据,作为填充数据;用填充数据填充指定角度间隔。
其中,用填充数据填充指定角度间隔,指将填充数据复制到指定角度间隔所在的期相中。
其中,角度间隔阈值可以等于该期相k空间中的所有角度间隔的平均角度间隔与指定倍数的乘积。
举例说明。假设当前期相共有300条数据线,那么平均的角度间隔为360度/300=1.2度。假设设置的阈值的系数为1.5,那么角度间隔阈值=1.2度*1.5=1.8度。那么,该期相k空间中所有超过1.8度的角度间隔均为指定角度间隔。
在一个示例性的实现过程中,从该期相的相邻期相对应的第二分配数据中,查找与指定角度间隔匹配的数据,作为填充数据,可以包括:获取第二分配数据对应的数据线,作为预选数据线;从预选数据线中,查找角度处于指定角度间隔内的数据线,作为匹配数据线,匹配数据线上的数据,即为与指定角度间隔匹配的数据;获取匹配数据线上的数据,作为填充数据。
在一个示例性的实现过程中,指定条件为:该期相k空间中的所有角度间隔均小于角度间隔阈值;或者,指定条件为:在该期相的相邻期相对应的第二分配数据中,不存在与指定角度间隔匹配的数据。
其中,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔的过程可以是循环进行的一个过程,直到满足指定条件时停止。也即,在利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的较大角度间隔后,改变角度间隔阈值(减小角度间隔阈值),重新检测该期相中是否还有较大角度间隔,如果有,继续利用该期相的相邻期相中的磁共振扫描数据来填充该期相k空间中的较大角度间隔,直到满足指定条件时停止。
举例说明。假设当前期相共有300条数据线,平均角度间隔为360度/300=1.2度,角度间隔阈值=1.2度*1.5=1.8度。那么,该期相k空间中所有超过1.8度的角度间隔均为指定角度间隔。假设当前期相中有22个角度超过1.8度的角度间隔,那么需要去相邻的两个期相里面遍历全部的数据线,看看数据线的角度是否能够落在这22个角度间隔当中。假设其中20个角度间隔能够被‘填上’,2个不能,则将对应的20条数据线的数据复制到当前期相当中。对于不能被‘填上’的空隙,则不做任何操作。
填充之后,当前期相有300+20=320条数据线,平均角度间隔为360度/320=1.125度。新的角度间隔阈值为1.125度*1.5=1.6875度。重复上述步骤,寻找当前期相的k空间角度分布中所有超过1.6875度的角度间隔,用相邻的两个期相里角度落在超过1.6875度的角度间隔中的数据线的数据填充相应的超过1.6875度的角度间隔,……,填充之后,重新计算角度间隔阈值和指定角度间隔的值,重复上述步骤,直到所有的角度间隔都小于角度间隔阈值时停止,或者所有的角度间隔都不可能再被‘填上’(即该期相的相邻期相对应的第二分配数据中,不存在与指定角度间隔匹配的数据)时停止。
在一个示例性的实现过程中,在满足指定条件时,获取该期相填充后的数据,作为待重建数据之前,电影成像方法还包括:判断是否满足指定条件。
在一个示例性的实现过程中,径向采集心脏的磁共振扫描数据,可以包括:黄金角径向采集心脏的磁共振扫描数据。在黄金角径向采集时,支持任意选择时间分辨率重建。
图3为本发明实施例提供的电影成像方法的第二流程示例图。如图3所示,本实施例中,电影成像方法可以包括如下步骤:
S301,径向采集心脏的磁共振扫描数据,并记录磁共振扫描数据对应的心跳信号。
S302,根据心跳信号,查找处于心跳周期范围外的心跳周期,作为待删除心跳周期。
S303,从磁共振扫描数据中,删除待删除心跳周期对应的数据,得到合格磁共振扫描数据。
S304,基于指定期相数目和心跳信号,将合格磁共振扫描数据分配到对应的期相中。
S305,对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;根据待重建数据进行图像重建,得到该期相的重建图像。将上述每一个期相的重建图像按照对应的期相顺序显示,即可获得心脏电影图像。
图3所示实施例的重建数据中,剔除了心律不齐时的数据(即待删除心跳周期对应的数据)进一步减弱了放射状伪影,从而进一步提高了图像质量。
图4为根据本发明实施例提供的电影成像方法得到的心脏短轴图像与根据现有技术得到的心脏短轴图像的对比示例图。图4中,a和c是舒张末期的心脏图像,b和d是收缩末期的心脏图像,a和b是根据现有技术方法得到的,c和d是根据本发明实施例提供的电影成像方法得到的。由图4可见,由于k空间非常不均匀的分布,a和b中有明显的放射状伪影,而c和d中的放射状伪影则减轻很多,c和d的图像质量明显优于a和b的图像质量。
图5为根据本发明实施例提供的电影成像方法得到的剔除心律不齐数据时的重建心脏图像与根据现有技术得到的剔除心律不齐数据时的重建心脏图像的对比示例图。图5中,第一行为根据现有技术得到的剔除心律不齐数据时的重建心脏图像,从左往右数据丢弃程度依次增加;第二行为根据本发明实施例提供的电影成像方法得到的剔除心律不齐数据时的重建心脏图像,从左往右数据丢弃程度依次增加。第一行和第二行中上下对应图像的数据丢弃程度相同。由图5可见,根据现有技术得到的图像随着数据丢弃程度的增加而显示出明显的伪影、图像质量严重受损,而根据本发明实施例提供的电影成像方法得到的图像则始终较好、可以被用于诊断。第二行最后一个图显示即使丢弃10个心跳周期时也有较好的质量,而这几乎意味着丢弃了59%的数据。
本发明实施例提供的电影成像方法,通过径向采集心脏的磁共振扫描数据,并记录磁共振扫描数据对应的心跳信号,基于指定期相数目和心跳信号,将磁共振扫描数据分配到对应的期相中,对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;根据待重建数据进行图像重建,得到该期相的重建图像,利用期相的相邻期相的数据对该期相k空间中的较大角度间隔进行填充,使得该期相k空间的角度方向分布趋于均匀,从而减弱了重建图像中的放射状伪影,因此提高了图像质量。
实施例二
本发明实施例提供一种磁共振成像***,该***包括处理器以及存储器;存储器用于存储指令,指令被处理器执行时,导致***实现实施例一中任一种电影成像方法。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***,装置和模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本发明所提供的几个实施例中,应该理解到,所揭露的***,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示以有另外的划分方式,例如,多个模块或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或模块的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能模块可以集成在一个处理单元中,也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元,可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中,包括若干指令用以使得一台计算机装置(可以是个人计算机,服务器,或者网络装置等)或处理器(Processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (10)
1.一种电影成像方法,其特征在于,所述方法包括:
径向采集心脏的磁共振扫描数据,并记录所述磁共振扫描数据对应的心跳信号;
基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中;
对于每一个期相,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔;在满足指定条件时,获取该期相填充后的数据,作为待重建数据;
根据所述待重建数据进行图像重建,得到该期相的重建图像。
2.根据权利要求1所述的方法,其特征在于,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:
根据所述心跳信号,确定所述磁共振扫描数据对应的各个心跳周期;
将所述各个心跳周期中的每个心跳周期平均分成指定期相数目个期相;
对于每一个期相,从所述磁共振扫描数据中提取采集时间处于该期相的数据,作为该期相对应的第一分配数据。
3.根据权利要求1所述的方法,其特征在于,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:
根据所述心跳信号,确定所述磁共振扫描数据对应的各个心跳周期;
将所述各个心跳周期中的每个心跳周期归一化到标准周期上;
将所述磁共振扫描数据的采集时间归一化到所述标准周期上,得到归一化采集时间;
将所述标准周期平均分成指定期相数目个期相;
对于每一个期相,从所述磁共振扫描数据中提取归一化采集时间处于该期相的数据,作为该期相对应的第一分配数据。
4.根据权利要求1所述的方法,其特征在于,利用该期相的相邻期相中的指定磁共振扫描数据,填充该期相k空间中的指定角度间隔,包括:
在该期相k空间中查找大于角度间隔阈值的角度间隔,作为指定角度间隔;
从该期相的相邻期相对应的第二分配数据中,查找与所述指定角度间隔匹配的数据,作为填充数据;
用所述填充数据填充所述指定角度间隔。
5.根据权利要求4所述的方法,其特征在于,从该期相的相邻期相对应的第二分配数据中,查找与所述指定角度间隔匹配的数据,作为填充数据,包括:
获取所述第二分配数据对应的数据线,作为预选数据线;
从所述预选数据线中,查找角度处于所述指定角度间隔内的数据线,作为匹配数据线,所述匹配数据线上的数据,即为与所述指定角度间隔匹配的数据;
获取所述匹配数据线上的数据,作为填充数据。
6.根据权利要求4所述的方法,其特征在于,
所述指定条件为:该期相k空间中的所有角度间隔均小于角度间隔阈值;或者,
所述指定条件为:在该期相的相邻期相对应的第二分配数据中,不存在与所述指定角度间隔匹配的数据。
7.根据权利要求1所述的方法,其特征在于,在满足指定条件时,获取该期相填充后的数据,作为待重建数据之前,所述方法还包括:
判断是否满足所述指定条件。
8.根据权利要求1所述的方法,其特征在于,基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中之前,所述方法还包括:
根据所述心跳信号,查找处于心跳周期范围外的心跳周期,作为待删除心跳周期;
从所述磁共振扫描数据中,删除所述待删除心跳周期对应的数据,得到合格磁共振扫描数据;
基于指定期相数目和所述心跳信号,将所述磁共振扫描数据分配到对应的期相中,包括:基于指定期相数目和所述心跳信号,将所述合格磁共振扫描数据分配到对应的期相中。
9.根据权利要求1所述的方法,其特征在于,径向采集心脏的磁共振扫描数据,包括:黄金角径向采集心脏的磁共振扫描数据。
10.一种磁共振成像***,其特征在于,所述***包括处理器以及存储器;所述存储器用于存储指令,所述指令被所述处理器执行时,导致所述***实现如权利要求1~9任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810601448.7A CN109001660B (zh) | 2018-06-12 | 2018-06-12 | 电影成像方法及磁共振成像*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810601448.7A CN109001660B (zh) | 2018-06-12 | 2018-06-12 | 电影成像方法及磁共振成像*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109001660A true CN109001660A (zh) | 2018-12-14 |
CN109001660B CN109001660B (zh) | 2020-07-28 |
Family
ID=64601257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810601448.7A Active CN109001660B (zh) | 2018-06-12 | 2018-06-12 | 电影成像方法及磁共振成像*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109001660B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109633503A (zh) * | 2018-12-28 | 2019-04-16 | 上海联影医疗科技有限公司 | 磁共振图像重建方法、装置、计算机设备和存储介质 |
CN112494029A (zh) * | 2019-11-29 | 2021-03-16 | 上海联影智能医疗科技有限公司 | 实时mr电影数据重建方法和*** |
CN116228802A (zh) * | 2023-05-05 | 2023-06-06 | 济南科汛智能科技有限公司 | 一种心脏mri辅助成像控制方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1245047A (zh) * | 1998-08-14 | 2000-02-23 | 通用电气公司 | 用于三维磁共振成象的分段k-空间方法 |
JP4347788B2 (ja) * | 2004-12-01 | 2009-10-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
CN101843487A (zh) * | 2009-03-23 | 2010-09-29 | 西门子公司 | 借助并行采集技术建立图像的方法、磁共振设备和程序 |
CN101843482A (zh) * | 2009-01-28 | 2010-09-29 | 西门子公司 | 确定心电图信号中r形脉冲的方法、测量装置和磁共振仪 |
CN101855564A (zh) * | 2007-11-09 | 2010-10-06 | 皇家飞利浦电子股份有限公司 | Mr-pet周期运动门控和校正 |
CN101874733A (zh) * | 2009-04-30 | 2010-11-03 | 西门子公司 | 用于触发的测量的方法以及磁共振断层造影设备 |
CN102215749A (zh) * | 2008-11-18 | 2011-10-12 | 株式会社日立医疗器械 | 磁共振成像装置及磁共振成像方法 |
CN103006218A (zh) * | 2012-12-22 | 2013-04-03 | 中国科学院深圳先进技术研究院 | 心脏磁共振实时电影成像的数据处理方法及*** |
CN103099618A (zh) * | 2011-11-15 | 2013-05-15 | 西门子公司 | 用于生成磁共振血管造影图像的方法和相应的磁共振设备 |
CN104181487A (zh) * | 2013-11-19 | 2014-12-03 | 上海联影医疗科技有限公司 | K空间重建方法及装置 |
CN105793722A (zh) * | 2013-12-02 | 2016-07-20 | 皇家飞利浦有限公司 | 用于稳态mr序列的实时自适应生理同步和门控 |
CN206209898U (zh) * | 2016-12-06 | 2017-05-31 | 中国科学院深圳先进技术研究院 | 三维心脏图像重建*** |
CN106910182A (zh) * | 2015-12-22 | 2017-06-30 | 上海联影医疗科技有限公司 | 心脏心功能磁共振图像中舒张末期图像中血池的分割方法 |
CN107576924A (zh) * | 2017-08-07 | 2018-01-12 | 上海东软医疗科技有限公司 | 一种磁共振动态成像方法和装置 |
-
2018
- 2018-06-12 CN CN201810601448.7A patent/CN109001660B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1245047A (zh) * | 1998-08-14 | 2000-02-23 | 通用电气公司 | 用于三维磁共振成象的分段k-空间方法 |
JP4347788B2 (ja) * | 2004-12-01 | 2009-10-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
CN101855564A (zh) * | 2007-11-09 | 2010-10-06 | 皇家飞利浦电子股份有限公司 | Mr-pet周期运动门控和校正 |
CN102215749A (zh) * | 2008-11-18 | 2011-10-12 | 株式会社日立医疗器械 | 磁共振成像装置及磁共振成像方法 |
CN101843482A (zh) * | 2009-01-28 | 2010-09-29 | 西门子公司 | 确定心电图信号中r形脉冲的方法、测量装置和磁共振仪 |
CN101843487A (zh) * | 2009-03-23 | 2010-09-29 | 西门子公司 | 借助并行采集技术建立图像的方法、磁共振设备和程序 |
CN101874733A (zh) * | 2009-04-30 | 2010-11-03 | 西门子公司 | 用于触发的测量的方法以及磁共振断层造影设备 |
CN103099618A (zh) * | 2011-11-15 | 2013-05-15 | 西门子公司 | 用于生成磁共振血管造影图像的方法和相应的磁共振设备 |
CN103006218A (zh) * | 2012-12-22 | 2013-04-03 | 中国科学院深圳先进技术研究院 | 心脏磁共振实时电影成像的数据处理方法及*** |
CN104181487A (zh) * | 2013-11-19 | 2014-12-03 | 上海联影医疗科技有限公司 | K空间重建方法及装置 |
CN105793722A (zh) * | 2013-12-02 | 2016-07-20 | 皇家飞利浦有限公司 | 用于稳态mr序列的实时自适应生理同步和门控 |
CN106910182A (zh) * | 2015-12-22 | 2017-06-30 | 上海联影医疗科技有限公司 | 心脏心功能磁共振图像中舒张末期图像中血池的分割方法 |
CN206209898U (zh) * | 2016-12-06 | 2017-05-31 | 中国科学院深圳先进技术研究院 | 三维心脏图像重建*** |
CN107576924A (zh) * | 2017-08-07 | 2018-01-12 | 上海东软医疗科技有限公司 | 一种磁共振动态成像方法和装置 |
Non-Patent Citations (1)
Title |
---|
李硕: "自门控心脏磁共振成像技术的发展与应用", 《中国医学物理学杂志》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109633503A (zh) * | 2018-12-28 | 2019-04-16 | 上海联影医疗科技有限公司 | 磁共振图像重建方法、装置、计算机设备和存储介质 |
CN112494029A (zh) * | 2019-11-29 | 2021-03-16 | 上海联影智能医疗科技有限公司 | 实时mr电影数据重建方法和*** |
CN116228802A (zh) * | 2023-05-05 | 2023-06-06 | 济南科汛智能科技有限公司 | 一种心脏mri辅助成像控制方法 |
CN116228802B (zh) * | 2023-05-05 | 2023-07-04 | 济南科汛智能科技有限公司 | 一种心脏mri辅助成像控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109001660B (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9918651B2 (en) | Electrocardiogram data analysis method for rapid diagnosis | |
JP5060719B2 (ja) | 生体臓器の画像処理装置、方法及びコンピュータプログラム | |
US9737229B1 (en) | Noninvasive electrocardiographic method for estimating mammalian cardiac chamber size and mechanical function | |
CN109001660A (zh) | 电影成像方法及磁共振成像*** | |
US8099158B2 (en) | System and method for selecting end of diastole and end of systole frames | |
JP2007516815A5 (zh) | ||
JP2007529272A (ja) | 心拍同期ct画像取得 | |
US10545208B2 (en) | System, method and computer-accessible medium for rapid real-time cardiac magnetic resonance imaging utilizing synchronized cardio-respiratory sparsity | |
US11269036B2 (en) | System and method for phase unwrapping for automatic cine DENSE strain analysis using phase predictions and region growing | |
US20220361799A1 (en) | Circuitless heart cycle determination | |
CN107978001B (zh) | 一种重建心脏ct图像的方法和装置 | |
CN106163396A (zh) | 用于呈现跳动心脏的运动模型的***和方法 | |
Suksaranjit et al. | Incidental LV LGE on CMR imaging in atrial fibrillation predicts recurrence after ablation therapy | |
CN113269845A (zh) | 图像重建方法、装置、存储介质及电子设备 | |
JP5060720B2 (ja) | 診断支援装置、方法、及びコンピュータプログラム | |
US20160345926A1 (en) | Acquisition of projection data for motion-corrected computed tomography images | |
CN110769748B (zh) | 伪迹宽容的脉搏率变异性测量 | |
CN103153197A (zh) | 超声波诊断装置、图像处理装置以及图像处理方法 | |
Sörnmo et al. | Characterization of f Waves | |
Chava et al. | Segmented radial cardiac MRI during arrhythmia using retrospective electrocardiogram and respiratory gating | |
WO2020221632A1 (en) | Electrocardiogram signal segmentation | |
JP6030002B2 (ja) | 画像処理プログラム、画像処理装置、および画像処理方法 | |
US8126902B2 (en) | Operating method for a computer having sequences of data sets changing over time and corresponding devices | |
Torres et al. | Flow-function relation in patients with chronic coronary artery disease and reduced regional function: a positron emission tomographic and two-dimensional echocardiographic study with coronary vasodilator stress | |
CN110897626A (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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |