CN108983133B - 磁共振血管壁成像方法、装置、设备及存储介质 - Google Patents

磁共振血管壁成像方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN108983133B
CN108983133B CN201810813453.4A CN201810813453A CN108983133B CN 108983133 B CN108983133 B CN 108983133B CN 201810813453 A CN201810813453 A CN 201810813453A CN 108983133 B CN108983133 B CN 108983133B
Authority
CN
China
Prior art keywords
echo
flip angle
echoes
pulse
phase
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
Application number
CN201810813453.4A
Other languages
English (en)
Other versions
CN108983133A (zh
Inventor
郑海荣
刘新
张娜
张磊
贾琳
贾文霄
邹超
梁栋
万丽雯
赵世华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201810813453.4A priority Critical patent/CN108983133B/zh
Publication of CN108983133A publication Critical patent/CN108983133A/zh
Priority to PCT/CN2019/088327 priority patent/WO2020019851A1/zh
Application granted granted Critical
Publication of CN108983133B publication Critical patent/CN108983133B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/5635Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Condensed Matter Physics & Semiconductors (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

本发明适用磁共振成像技术领域,提供了一种磁共振血管壁成像方法、装置、设备及存储介质,该方法包括:将回波链划分为由多个回波组成的期相,根据每个期相内第一个回波对应的第一回聚脉冲翻转角度、最后一个回波对应的第二回聚脉冲翻转角度以及翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链,根据可变翻转角链计算回波链对应的回波信号演化,当回波信号演化满足约束条件时,根据回波信号演化重构磁共振血管壁图像,否则,进行参数调整并继续计算每个期相内所有回波对应的回聚脉冲翻转角度,以获得最优的可变翻转角链,从而提高了在3T磁场中磁共振血管壁成像的信噪比。

Description

磁共振血管壁成像方法、装置、设备及存储介质
技术领域
本发明属于磁共振成像技术领域,尤其涉及一种磁共振血管壁成像方法、装置、设备及存储介质。
背景技术
心脑血管疾病已成为威胁人类健康的首要疾病,而动脉粥样硬化斑块是导致心脑血管疾病发生的一个重要风险因素,因此准确地评估斑块的稳定性、鉴别易损斑块对于预防脑卒事件具有重要意义,基于影像学手段的血管斑块监测,对于心脑血管疾病的预测、分期和预后评估有着非常重要的意义。目前常用于头颈部血管成像的方法有多普勒超声、磁共振血管造影术(Magnetic Resonance Angiography,MRA)、CT血管造影术(ComputedTomography Angiography,CTA)、数字减影血管造影术(Digital SubtractionAngiography,DSA)等,而以上这些技术只局限于对血管管腔通畅程度的显示,在血管壁成分特征信息的评估方面存在很大局限性。磁共振血管壁成像(Magnetic Resonance VesselWall Imaging)是目前唯一能够无创、全景显示头颈联合血管壁和斑块的技术。该技术基于磁共振原理,通过抑制血管内流动血液的信号对动脉粥样硬化斑块的形态和成分(例如,斑块内出血(Intra-plaque Hemorrhage,IPH)、钙化(Calcification,CA)、脂质核(LipidRich Necrotic Core,LRNC)、纤维帽(Fibrous cap,FC)等)进行评估,进而达到对血管斑块定量分析的目的,对于评估斑块的稳定性和风险预测有积极的作用。因此,磁共振血管壁成像作为传统血管成像的补充和优化,可以全面评价血管管壁,帮助鉴别管腔狭窄原因并进一步评价易损斑块的组成成分,具有较好的发展前景。
随着磁共振技术的发展,越来越多的三维成像序列被开发和应用于临床。可变翻转角快速自旋回波成像(Sampling Perfection with Application-optimized Contrastsby using different flip angle Evolutions,SPACE)序列具备各向同性扫描、分辨率高的特点,开始被应用于动脉粥样硬化斑块成像。
SPACE可变翻转角的设计目的是使目标信号按照“既定回波信号演化(PrescribedSignal Evolution)”的要求进行衰减,根据目标信号的演化和拓展相位图算法(ExtendedPhase Graph Algorithm,EPG Algorithm),逆解得到上述目标信号所需的翻转角链,最后根据解出的翻转角链,利用布洛赫方程(Bloch Equation),便可得到不同生物体组织的实际信号演化。而上述可变翻转角所依据的“既定回波信号演化”是根据1.5T磁场下大脑灰质的弛豫特性设计的,这对于3T磁场下的血管壁成像不是最佳选择。
发明内容
本发明的目的在于提供一种磁共振血管壁成像方法、装置、设备及存储介质,旨在解决由于现有技术无法提供一种有效的磁共振血管壁成像方法,导致在3T磁场下磁共振血管壁成像不清晰、且血管壁斑块识别度低的问题。
一方面,本发明提供了一种磁共振血管壁成像方法,所述方法包括下述步骤:
当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成;
根据所述每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算所述每个期相内所有回波对应的回聚脉冲翻转角度,以构成所述回波链对应的可变翻转角链;
根据所述可变翻转角链计算所述回波链对应的回波信号演化;
判断所述回波信号演化是否满足预设的约束条件,是则,根据所述回波信号演化重构磁共振血管壁图像,否则,调整所述第一回聚脉冲翻转角度、所述第二回聚脉冲翻转角度以及所述翻转角控制参数,跳转至计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤。
另一方面,本发明提供了一种磁共振血管壁成像装置,所述装置包括:
回波链划分单元,用于当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成;
翻转角度计算单元,用于根据所述每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算所述每个期相内所有回波对应的回聚脉冲翻转角度,以构成所述回波链对应的可变翻转角链;
信号演化计算单元,用于根据所述可变翻转角链计算所述回波链对应的回波信号演化;以及
信号演化判断单元,用于判断所述回波信号演化是否满足预设的约束条件,是则,根据所述回波信号演化重构磁共振血管壁图像,否则,调整所述第一回聚脉冲翻转角度、所述第二回聚脉冲翻转角度以及所述翻转角控制参数,触发所述翻转角度计算单元执行计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤。
另一方面,本发明还提供了一种医疗设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上述磁共振血管壁成像方法所述的步骤。
另一方面,本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如上述磁共振血管壁成像方法所述的步骤。
本发明将回波链划分为由多个回波组成的期相,根据每个期相内第一个回波对应的第一回聚脉冲翻转角度、最后一个回波对应的第二回聚脉冲翻转角度以及翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链,根据可变翻转角链计算回波链对应的回波信号演化,判断回波信号演化是否满足预设的约束条件,是则,根据回波信号演化重构磁共振血管壁图像,否则,进行参数调整并继续计算每个期相内所有回波对应的回聚脉冲翻转角度,以获得最优的可变翻转角链,从而提高了在3T磁场中磁共振血管壁成像的信噪比,进而提高了血管壁病变特别是斑块的识别度和治疗效果。
附图说明
图1是本发明实施例一提供的磁共振血管壁成像方法的实现流程图;
图2是本发明实施例二提供的磁共振血管壁成像装置的结构示意图;
图3是本发明实施例二提供的磁共振血管壁成像装置的优选结构示意图;以及
图4是本发明实施例三提供的医疗设备的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以下结合具体实施例对本发明的具体实现进行详细描述:
实施例一:
图1示出了本发明实施例一提供的磁共振血管壁成像方法的实现流程,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在步骤S101中,当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成。
本发明实施例适用于医学图像处理平台、系统或医学设备,例如个人计算机、服务器等。当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为一个或者多个期相,每个期相内包含多个回波,这些回波是由多个连续的回聚脉冲产生的。
在将预先采集到的回波链划分为预设数量个期相之前,优选地,发射预设的回聚射频脉冲以产生回波,对产生的回波进行采集,获得对应的回波链,从而提高后续可变翻转角链的优化效果。作为示例地,在一个重复时间内,首先发射90°的射频脉冲,然后发射多个180°回聚脉冲,以产生多个回波,再对产生的回波进行采集,由采集的回波形成一个回波链。
在步骤S102中,根据每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链。
在本发明实施例中,在任意一个期相内,根据第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚翻转角度以及预设翻转角控制参数,计算该期相内所有回波对应的回聚脉冲翻转角度,由所有期相内所有回波对应的回聚脉冲翻转角度构成整个回波链对应的可变翻转角链。其中,每个期相内第一个回波对应的第一回聚脉冲翻转角度(或者最后一个回波对应的第二回聚脉冲翻转角度)可以预先设置为相同的角度,也可以设置为不同的角度。翻转角控制参数用于控制每个期相内回聚脉冲翻转角度在期相端点附近的变化速度,通常地,从低值的端点到高值的端点,回聚脉冲翻转角度由慢到快发生变化,不同期相内的翻转角控制参数可以是相同的,也可以是不同的。
优选地,翻转角控制参数为大于1的实数,从而提高对回聚脉冲翻转角度变化速度的控制精确度以及回聚脉冲翻转角度变化的适应度。
在计算每个期相内所有回波对应的回聚脉冲翻转角度时,优选地,当第一回聚脉冲翻转角度大于第二回聚脉冲翻转角度时,根据公式αn=αN-1+(α0N-1)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,从而降低了T2衰减带来的模糊效应程度,并提高了在3T磁场中血液信号的抑制程度,进而提高后续在3T磁场中磁共振血管壁成像的信噪比。其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000061
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1}。
在计算每个期相内所有回波对应的回聚脉冲翻转角度时,又一优选地,当第一回聚脉冲翻转角度不超过第二回聚脉冲翻转角度时,根据公式αn=α0+(αN-10)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,从而降低了T2衰减带来的模糊效应程度,并提高了在3T磁场中血液信号的抑制程度,进而提高后续在3T磁场中磁共振血管壁成像的信噪比。其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000062
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1}。
在步骤S103中,根据可变翻转角链计算回波链对应的回波信号演化。
在本发明实施例中,根据得到的可变翻转角链,再利用布洛赫方程(BlochEquation)计算得到回波链对应的回波信号演化。
在步骤S104中,判断回波信号演化是否满足预设的约束条件。
在步骤S105中,根据回波信号演化重构磁共振血管壁图像。
在步骤S106中,调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数。
在本发明实施例中,判断回波信号演化是否满足预设的约束条件(例如,根据回波信号演化分析官腔、管壁及周围组织的对比度是否满足预设的对比度阈值,或者根据回波信号演化分析射频能量吸收率(Specific Absorption Rate,SAR)是否满足预设的吸收率阈值),是则,执行步骤S105根据该回波信号演化重构磁共振血管壁图像,否则,执行步骤S106调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数,并跳转至步骤S102继续执行计算每个期相内所有回波对应的回聚脉冲翻转角度,以获取最优的可变翻转角链。
在调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数时,优选地,在预设的翻转角度范围内对第一回聚脉冲翻转角度和第二回聚脉冲翻转角度进行调整,同时在1至预设的最大临界值范围内对翻转角控制参数进行调整,具体地,在0°至180°范围内对第一回聚脉冲翻转角度和第二回聚脉冲翻转角度进行随机调整,从而提高可变翻转角链的优化速度。
在本发明实施例中,将回波链划分为由多个回波组成的期相,根据每个期相内第一个回波对应的第一回聚脉冲翻转角度、最后一个回波对应的第二回聚脉冲翻转角度以及翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链,根据可变翻转角链计算回波链对应的回波信号演化,当回波信号演化满足约束条件时,根据回波信号演化重构磁共振血管壁图像,否则,进行参数调整,并继续计算每个期相内所有回波对应的回聚脉冲翻转角度,以获得最优的可变翻转角链,从而提高了在3T磁场中磁共振血管壁成像的信噪比,进而提高了血管壁病变特别是斑块的识别度和治疗效果。
实施例二:
图2示出了本发明实施例二提供的磁共振血管壁成像装置的结构,为了便于说明,仅示出了与本发明实施例相关的部分,其中包括:
回波链划分单元21,用于当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成。
本发明实施例适用于医学图像处理平台、系统或设备,例如个人计算机、服务器等。当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为一个或者多个期相,每个期相内包含多个回波,这些回波是由多个连续的回聚脉冲产生的。
在将预先采集到的回波链划分为预设数量个期相之前,优选地,发射预设的回聚射频脉冲以产生回波,对产生的回波进行采集,获得对应的回波链,从而提高后续可变翻转角链的优化效果。作为示例地,在一个重复时间内,首先发射90°的射频脉冲,然后发射多个180°回聚脉冲,以产生多个回波,再对产生的回波进行采集,由采集的回波形成一个回波链。
翻转角度计算单元22,用于根据每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链。
在本发明实施例中,在任意一个期相内,根据第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚翻转角度以及预设翻转角控制参数,计算该期相内所有回波对应的回聚脉冲翻转角度,由所有期相内所有回波对应的回聚脉冲翻转角度构成整个回波链对应的可变翻转角链。其中,每个期相内第一个回波对应的第一回聚脉冲翻转角度(或者最后一个回波对应的第二回聚脉冲翻转角度)可以预先设置为相同的角度,也可以设置为不同的角度。翻转角控制参数用于控制每个期相内回聚脉冲翻转角度在期相端点附近的变化速度,通常地,从低值的端点到高值的端点,回聚脉冲翻转角度由慢到快发生变化,不同期相内的翻转角控制参数可以是相同的,也可以是不同的。
优选地,翻转角控制参数为大于1的实数,从而提高对回聚脉冲翻转角度变化速度的控制精确度以及回聚脉冲翻转角度变化的适应度。
在计算每个期相内所有回波对应的回聚脉冲翻转角度时,优选地,当第一回聚脉冲翻转角度大于第二回聚脉冲翻转角度时,根据公式αn=αN-1+(α0N-1)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,从而降低了T2衰减带来的模糊效应程度,并提高了在3T磁场中血液信号的抑制程度,进而提高后续在3T磁场中磁共振血管壁成像的信噪比。其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000081
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1}。
在计算每个期相内所有回波对应的回聚脉冲翻转角度时,又一优选地,当第一回聚脉冲翻转角度不超过第二回聚脉冲翻转角度时,根据公式αn=α0+(αN-10)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,从而降低了T2衰减带来的模糊效应程度,并提高了在3T磁场中血液信号的抑制程度,进而提高后续在3T磁场中磁共振血管壁成像的信噪比。其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000091
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1}。
信号演化计算单元23,用于根据可变翻转角链计算回波链对应的回波信号演化。
在本发明实施例中,根据得到的可变翻转角链,再利用布洛赫方程(BlochEquation)计算得到回波链对应的回波信号演化。
信号演化判断单元24,用于判断回波信号演化是否满足预设的约束条件,是则,根据回波信号演化重构磁共振血管壁图像,否则,调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数,触发翻转角度计算单元22执行通过预设的回聚脉冲翻转角度计算公式计算每个期相内所有回波对应的回聚脉冲翻转角度的步骤。
在本发明实施例中,判断回波信号演化是否满足预设的约束条件(例如,根据回波信号演化分析官腔、管壁及周围组织的对比度是否满足预设的对比度阈值,或者根据回波信号演化分析射频能量吸收率(Specific Absorption Rate,SAR)是否满足预设的吸收率阈值),是则,根据该回波信号演化重构磁共振血管壁图像,否则,调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数,触发翻转角度计算单元22执行计算每个期相内所有回波对应的回聚脉冲翻转角度的步骤。
在调整第一回聚脉冲翻转角度、第二回聚脉冲翻转角度以及翻转角控制参数时,优选地,在预设的翻转角度范围内对第一回聚脉冲翻转角度和第二回聚脉冲翻转角度进行调整,同时在1至预设的最大临界值范围内对翻转角控制参数进行调整,具体地,在0°至180°范围内对第一回聚脉冲翻转角度和第二回聚脉冲翻转角度进行随机调整,从而提高可变翻转角链的优化速度。
因此,如图3所示,优选地,翻转角度计算单元22包括:
第一角度计算单元221,用于当第一回聚脉冲翻转角度大于第二回聚脉冲翻转角度时,根据公式αn=αN-1+(α0N-1)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000101
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1};以及
第二角度计算单元222,用于当第一回聚脉冲翻转角度不超过第二回聚脉冲翻转角度时,根据公式αn=α0+(αN-10)×C(n)计算每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure BDA0001739759810000102
N为每个期限内回波的个数,p为翻转角控制参数,n∈{0,1,2,...,N-1}。
在本发明实施例中,磁共振血管壁成像装置的各单元可由相应的硬件或软件单元实现,各单元可以为独立的软、硬件单元,也可以集成为一个软、硬件单元,在此不用以限制本发明。
实施例三:
图4示出了本发明实施例三提供的医疗设备的结构,为了便于说明,仅示出了与本发明实施例相关的部分。
本发明实施例的医疗设备4包括处理器40、存储器41以及存储在存储器41中并可在处理器40上运行的计算机程序42。该处理器40执行计算机程序42时实现上述磁共振血管壁成像方法实施例中的步骤,例如图1所示的步骤S101至S106。或者,处理器40执行计算机程序42时实现上述各装置实施例中各单元的功能,例如图2所示单元21至24的功能。
在本发明实施例中,将回波链划分为由多个回波组成的期相,根据每个期相内第一个回波对应的第一回聚脉冲翻转角度、最后一个回波对应的第二回聚脉冲翻转角度以及翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链,根据可变翻转角链计算回波链对应的回波信号演化,判断回波信号演化是否满足预设的约束条件,是则,根据回波信号演化重构磁共振血管壁图像,否则,进行参数调整,并继续计算每个期相内所有回波对应的回聚脉冲翻转角度,以获得最优的可变翻转角链,从而提高了在3T磁场中磁共振血管壁成像的信噪比,进而提高了血管壁病变特别是斑块的识别度和治疗效果。
本发明实施例的医疗设备可以为个人计算机、服务器。该医疗设备4中处理器40执行计算机程序42时实现磁共振血管壁成像方法时实现的步骤可参考前述方法实施例的描述,在此不再赘述。
实施例四:
在本发明实施例中,提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现上述磁共振血管壁成像方法实施例中的步骤,例如,图1所示的步骤S101至S106。或者,该计算机程序被处理器执行时实现上述各装置实施例中各单元的功能,例如图2所示单元21至24的功能。
在本发明实施例中,将回波链划分为由多个回波组成的期相,根据每个期相内第一个回波对应的第一回聚脉冲翻转角度、最后一个回波对应的第二回聚脉冲翻转角度以及翻转角控制参数,计算每个期相内所有回波对应的回聚脉冲翻转角度,以构成回波链对应的可变翻转角链,根据可变翻转角链计算回波链对应的回波信号演化,判断回波信号演化是否满足预设的约束条件,是则,根据回波信号演化重构磁共振血管壁图像,否则,进行参数调整,并继续计算每个期相内所有回波对应的回聚脉冲翻转角度,以获得最优的可变翻转角链,从而提高了在3T磁场中磁共振血管壁成像的信噪比,进而提高了血管壁病变特别是斑块的识别度和治疗效果。
本发明实施例的计算机可读存储介质可以包括能够携带计算机程序代码的任何实体或装置、记录介质,例如,ROM/RAM、磁盘、光盘、闪存等存储器。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种磁共振血管壁成像方法,其特征在于,所述方法包括下述步骤:
当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成;
根据所述每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算所述每个期相内所有回波对应的回聚脉冲翻转角度,以构成所述回波链对应的可变翻转角链;
根据所述可变翻转角链计算所述回波链对应的回波信号演化;
判断所述回波信号演化是否满足预设的约束条件,是则,根据所述回波信号演化重构磁共振血管壁图像,否则,调整所述第一回聚脉冲翻转角度、所述第二回聚脉冲翻转角度以及所述翻转角控制参数,跳转至计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤;
其中,所述判断所述回波信号演化是否满足预设的约束条件具体为根据回波信号演化分析管腔、管壁及周围组织的对比度是否满足预设的对比度阈值,或者根据回波信号演化分析射频能量吸收率是否满足预设的吸收率阈值;
所述调整具体为在0°至180°范围内对所述第一回聚脉冲翻转角度和所述第二回聚脉冲翻转角度进行随机调整,在1至预设的最大临界值范围内对翻转角控制参数进行调整。
2.如权利要求1所述的方法,其特征在于,计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤,包括:
当所述第一回聚脉冲翻转角度大于所述第二回聚脉冲翻转角度时,根据公式αn=αN-1+(α0N-1)×C(n)计算所述每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure FDA0002423794740000011
N为所述每个期限内回波的个数,p为所述翻转角控制参数,n∈{0,1,2,...,N-1}。
3.如权利要求1所述的方法,其特征在于,计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤,还包括:
当所述第一回聚脉冲翻转角度不超过所述第二回聚脉冲翻转角度时,根据公式αn=α0+(αN-10)×C(n)计算所述每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure FDA0002423794740000021
N为所述每个期限内回波的个数,p为所述翻转角控制参数,n∈{0,1,2,...,N-1}。
4.如权利要求1所述的方法,其特征在于,所述翻转角控制参数为大于1的实数。
5.一种磁共振血管壁成像装置,其特征在于,所述装置包括:
回波链划分单元,用于当接收到磁共振血管壁成像请求时,将预先采集到的回波链划分为预设数量个期相,每个期相由多个回聚脉冲产生的回波组成;
翻转角度计算单元,用于根据所述每个期相内第一个回波对应的预设第一回聚脉冲翻转角度、最后一个回波对应的预设第二回聚脉冲翻转角度以及预设的翻转角控制参数,计算所述每个期相内所有回波对应的回聚脉冲翻转角度,以构成所述回波链对应的可变翻转角链;
信号演化计算单元,用于根据所述可变翻转角链计算所述回波链对应的回波信号演化;以及
信号演化判断单元,用于判断所述回波信号演化是否满足预设的约束条件,是则,根据所述回波信号演化重构磁共振血管壁图像,否则,调整所述第一回聚脉冲翻转角度、所述第二回聚脉冲翻转角度以及所述翻转角控制参数,触发所述翻转角度计算单元执行计算所述每个期相内所有回波对应的回聚脉冲翻转角度的步骤;
其中,所述判断所述回波信号演化是否满足预设的约束条件具体为根据回波信号演化分析管腔、管壁及周围组织的对比度是否满足预设的对比度阈值,或者根据回波信号演化分析射频能量吸收率是否满足预设的吸收率阈值;
所述调整具体为在0°至180°范围内对所述第一回聚脉冲翻转角度和所述第二回聚脉冲翻转角度进行随机调整,在1至预设的最大临界值范围内对翻转角控制参数进行调整。
6.如权利要求5所述的装置,其特征在于,所述翻转角度计算单元包括:
第一角度计算单元,用于当所述第一回聚脉冲翻转角度大于所述第二回聚脉冲翻转角度时,根据公式αn=αN-1+(α0N-1)×C(n)计算所述每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure FDA0002423794740000031
N为所述每个期限内回波的个数,p为所述翻转角控制参数,n∈{0,1,2,...,N-1}。
7.如权利要求5所述的装置,其特征在于,所述翻转角度计算单元还包括:
第二角度计算单元,用于当所述第一回聚脉冲翻转角度不超过所述第二回聚脉冲翻转角度时,根据公式αn=α0+(αN-10)×C(n)计算所述每个期相内所有回波对应的回聚脉冲翻转角度,其中,αn为第n+1个回波对应的回聚脉冲翻转角度,
Figure FDA0002423794740000032
N为所述每个期限内回波的个数,p为所述翻转角控制参数,n∈{0,1,2,...,N-1}。
8.如权利要求5所述的装置,其特征在于,所述翻转角控制参数为大于1的实数。
9.一种医疗设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至4任一项所述方法的步骤。
CN201810813453.4A 2018-07-23 2018-07-23 磁共振血管壁成像方法、装置、设备及存储介质 Active CN108983133B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810813453.4A CN108983133B (zh) 2018-07-23 2018-07-23 磁共振血管壁成像方法、装置、设备及存储介质
PCT/CN2019/088327 WO2020019851A1 (zh) 2018-07-23 2019-05-24 磁共振血管壁成像方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810813453.4A CN108983133B (zh) 2018-07-23 2018-07-23 磁共振血管壁成像方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN108983133A CN108983133A (zh) 2018-12-11
CN108983133B true CN108983133B (zh) 2020-06-12

Family

ID=64550208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810813453.4A Active CN108983133B (zh) 2018-07-23 2018-07-23 磁共振血管壁成像方法、装置、设备及存储介质

Country Status (2)

Country Link
CN (1) CN108983133B (zh)
WO (1) WO2020019851A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108983133B (zh) * 2018-07-23 2020-06-12 中国科学院深圳先进技术研究院 磁共振血管壁成像方法、装置、设备及存储介质
CN111381204B (zh) * 2018-12-28 2022-07-12 西门子(深圳)磁共振有限公司 基于二维快速自旋回波的磁共振成像方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102217935A (zh) * 2010-04-13 2011-10-19 西门子公司 用于产生血管造影的磁共振图像的方法
JP2013078562A (ja) * 2011-09-22 2013-05-02 Toshiba Corp 磁気共鳴イメージング装置
EP2741096A1 (en) * 2012-12-04 2014-06-11 Samsung Electronics Co., Ltd Fast B1 mapping based on interleaved three flip angle excitation
EP2844143A1 (en) * 2012-05-04 2015-03-11 Emory University Methods, systems and computer readable storage media storing instructions for imaging and determining information associated with regions of the brain
CN105232046A (zh) * 2015-11-26 2016-01-13 厦门大学 一种基于重叠回波的单扫描定量磁共振t2成像方法
CN107638179A (zh) * 2017-08-28 2018-01-30 上海联影医疗科技有限公司 磁共振翻转角的计算方法、系统和计算机可读存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8676296B2 (en) * 2008-09-15 2014-03-18 Univerity of Pittsburgh—Of The Commonwealth System of Higher Education Echo-specific K-space reordering approach to compatible dual-echo arteriovenography
US9625548B2 (en) * 2012-06-19 2017-04-18 The Johns Hopkins University System and method for magnetic resonance imaging of intracranial vessel walls
CN105596003B (zh) * 2015-12-24 2019-08-02 深圳先进技术研究院 一种下肢深静脉血栓磁共振成像方法和装置
WO2018085969A1 (zh) * 2016-11-08 2018-05-17 中国科学院深圳先进技术研究院 磁共振血管壁成像方法和设备
CN108983133B (zh) * 2018-07-23 2020-06-12 中国科学院深圳先进技术研究院 磁共振血管壁成像方法、装置、设备及存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102217935A (zh) * 2010-04-13 2011-10-19 西门子公司 用于产生血管造影的磁共振图像的方法
JP2013078562A (ja) * 2011-09-22 2013-05-02 Toshiba Corp 磁気共鳴イメージング装置
EP2844143A1 (en) * 2012-05-04 2015-03-11 Emory University Methods, systems and computer readable storage media storing instructions for imaging and determining information associated with regions of the brain
EP2741096A1 (en) * 2012-12-04 2014-06-11 Samsung Electronics Co., Ltd Fast B1 mapping based on interleaved three flip angle excitation
CN105232046A (zh) * 2015-11-26 2016-01-13 厦门大学 一种基于重叠回波的单扫描定量磁共振t2成像方法
CN107638179A (zh) * 2017-08-28 2018-01-30 上海联影医疗科技有限公司 磁共振翻转角的计算方法、系统和计算机可读存储介质

Also Published As

Publication number Publication date
WO2020019851A1 (zh) 2020-01-30
CN108983133A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
US11635514B2 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
Romero et al. Cardiac imaging for assessment of left atrial appendage stasis and thrombosis
JP5750452B2 (ja) 時間分解磁気共鳴血管造影と潅流画像のシステムと方法
JP2981675B2 (ja) 磁気共鳴結像
US9335410B2 (en) System and method for multiple spotlight synthetic radar imaging using random beam steering
EP1269212B1 (en) Magnetic resonance angiography with automated vessel segmentation
US20110208061A1 (en) Ultrasonic lesion identification using temporal parametric contrast images
US20140147013A1 (en) Direct echo particle image velocimetry flow vector mapping on ultrasound dicom images
US20160135775A1 (en) System And Method For Time-Resolved, Three-Dimensional Angiography With Physiological Information
US10134127B2 (en) Method for post-processing flow-sensitive phase contrast magnetic resonance images
EP2291121A1 (en) System and method for characterizing tissue based upon homomorphic deconvolution of backscattered ultrasound
JP2022510217A (ja) 磁気共鳴温度測定における測定アーチファクトを補正するためのシステムおよび方法
JP2018512238A (ja) 傾斜磁場パルスmriにより得られた信号パターンから観察対象組織のシグネチャ指標を判定するmri法
CN108983133B (zh) 磁共振血管壁成像方法、装置、设备及存储介质
Richy et al. Blood velocity estimation using compressive sensing
Stanziola et al. Ultrasound imaging with microbubbles [life sciences]
Taghavi et al. Ultrasound super-resolution imaging with a hierarchical Kalman tracker
Zachiu et al. An adaptive non-local-means filter for real-time MR-thermometry
Park et al. Pulse inversion chirp coded tissue harmonic imaging (PI-CTHI) of zebrafish heart using high frame rate ultrasound biomicroscopy
US20150196223A1 (en) System and method for neutral contrast magnetic resonance imaging of calcifications
Yiu et al. High frame rate doppler ultrasound bandwidth imaging for flow instability mapping
Baltes et al. Determination of peak velocity in stenotic areas: echocardiography versus kt SENSE accelerated MR Fourier velocity encoding
US20120268128A1 (en) Method for Principal Frequency Magnetic Resonance Elastography Inversion
Pan et al. Improving axial resolution based on the deconvolution recovery method combined with adaptive weighting techniques for ultrasound imaging
JP7357007B2 (ja) 熱アブレーションのレベルを推定するための装置及び方法

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
TR01 Transfer of patent right

Effective date of registration: 20230223

Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 518000 No. 1068, Xue Yuan Avenue, Shenzhen University Town, Nanshan District, Shenzhen, Guangdong.

Patentee before: SHENZHEN INSTITUTES OF ADVANCED TECHNOLOGY CHINESE ACADEMY OF SCIENCES

TR01 Transfer of patent right