CN107638179B - 磁共振翻转角的计算方法、***和计算机可读存储介质 - Google Patents

磁共振翻转角的计算方法、***和计算机可读存储介质 Download PDF

Info

Publication number
CN107638179B
CN107638179B CN201710752718.XA CN201710752718A CN107638179B CN 107638179 B CN107638179 B CN 107638179B CN 201710752718 A CN201710752718 A CN 201710752718A CN 107638179 B CN107638179 B CN 107638179B
Authority
CN
China
Prior art keywords
flip angle
angle configuration
criterion
echo
radio frequency
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
CN201710752718.XA
Other languages
English (en)
Other versions
CN107638179A (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
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201710752718.XA priority Critical patent/CN107638179B/zh
Priority to US15/692,148 priority patent/US10739430B2/en
Priority to EP17204380.4A priority patent/EP3451009B1/en
Priority to PCT/CN2017/114107 priority patent/WO2019041607A1/en
Publication of CN107638179A publication Critical patent/CN107638179A/zh
Application granted granted Critical
Publication of CN107638179B publication Critical patent/CN107638179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • 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/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/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • 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/288—Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room
    • 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/4808—Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]

Landscapes

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

本申请公开了一种磁共振翻转角的计算方法、***和计算机可读存储介质。该方法包括：获取回聚射频脉冲的初始翻转角配置；判断所述初始翻转角配置是否满足第一准则；基于第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；判断所述第一翻转角配置是否满足第二准则，所述第二准则与所述第一翻转角配置对应的回波链的第二参数有关；基于第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及基于所述第二翻转角配置，获取核磁共振信号。本发明中公开的方法和***可以减小回聚射频脉冲的特殊吸收率，同时保持MR图像信噪比和对比度基本不变。

Description

磁共振翻转角的计算方法、***和计算机可读存储介质
技术领域
本申请涉及核磁共振成像领域，特别涉及一种用于确定目标翻转角配置的方法及***。
背景技术
核磁共振成像(MRI)***通过利用高强度磁场和射频技术广泛用于诊断和医学治疗。射频技术可以通过利用激发射频脉冲和一定数量的回聚射频脉冲获取核磁共振信号，用于重建核磁共振图像。在一些实施例中，多个回聚射频脉冲可以具有相同或可变的翻转角。若使用大量高频(即大翻转角)的回聚射频脉冲(RF)，会导致回聚射频脉冲的特殊吸收率(SAR)增加。若使用可变翻转角的回聚射频脉冲(RF)，可能会导致核磁共振(MR)图像的信噪比和/ 或对比度降低。因此，本申请针对以上问题，提出用于确定目标翻转角配置的方法及***，在降低特殊吸收率(SAR)的同时，保持图像对比度和信噪比基本不变。
发明内容
针对上述使用可变翻转角的回聚射频脉冲可能导致MR图像中信噪比或对比度下降的问题，本发明在于提供一种确定核磁共振成像中翻转角的方法的***和方法，可以减小回聚射频脉冲的特殊吸收率，同时保持MR图像信噪比和对比度基本不变。
为达到上述发明目的，本发明提供的技术方案如下：
一种用于获取核磁共振信号的方法，所述方法在至少一个机器上执行，所述至少一个机器中的每一个机器具有至少一个处理器和一个存储器。所述方法包括：获取回聚射频脉冲的初始翻转角配置，所述回聚射频脉冲用于产生回波链；判断所述初始翻转角配置是否满足第一准则以获得第一判断结果，所述第一准则与所述初始翻转角配置对应的回波链的第一参数有关；基于所述第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；判断所述第一翻转角配置是否满足第二准则以获得第二判断结果，所述第二准则与所述第一翻转角配置对应的回波链的第二参数有关；基于所述第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及基于所述第二翻转角配置，获取核磁共振信号。
在本发明中，所述获取回聚射频脉冲的初始翻转角配置包括确定所述回聚射频脉冲的初始条件，所述初始条件包括第一翻转角和第二翻转角的初始值；以及基于所述第一翻转角和所述第二翻转角的初始值，确定所述初始翻转角配置。
在本发明中，进一步包括获取与至少部分的所述回波链对应的函数，所述函数与所述第一翻转角和所述第二翻转角中至少一个相关；以及基于所述函数确定至少一个回聚射频脉冲的翻转角，所述至少一个回聚射频脉冲的翻转角与所述回波链中的一个回波对应。
在本发明中，所述第一参数包括所述回波链中至少一个回波的强度值。
在本发明中，所述基于所示第一判断结果，确定第一翻转角配置，进一步包括确定所述初始翻转角配置不满足所述第一准则；以及调整所述初始翻转角配置以确定第一所述第一翻转角配置。
在本发明中，所述判断所述初始翻转角配置是否满足第一准则以获得第一判断结果包括基于所述初始翻转角配置，确定第一信号演化；以及通过确定所述第一信号演变不满足所述第一准则确定所述初始翻转角配置是不满足所述第一准则，其中，所述第一信号演变包括所述回波链中至少部分的回波的强度值。
在本发明中，所述第一信号演变与横向弛豫时间或纵向弛豫时间有关。
在本发明中，所述第一准则包括第一阈值，所述第一阈值与所述第一信号演变有关，以及所述确定所述初始翻转角配置不满足所述第一准则包括确定所述第一信号演变中最大强度值小于所述第一阈值。
在本发明中，所述第一参数包括具有所述初始翻转角配置的回聚射频脉冲的总能量。
在本发明中，所述判断所述初始翻转角配置是否满足第一准则以获得第一判断结果包括基于所述初始翻转角配置，确定所述回聚射频脉冲的总能量；以及基于所述第一准则评价所述回聚射频脉冲的总能量。
在本发明中，所述第一准则包括第二阈值，所述第二阈值对应所述回聚射频脉冲的总能量，以及确定所述初始翻转角配置不满足所述第一准则包括确定所述回聚射频脉冲的总能量等于或大于所述第二阈值。
在本发明中，所述基于所述第二判断结果，确定第二翻转角配置进一步包括确定所述第一翻转角配置不满足所述第二准则；以及调整所述第一翻转角配置以确定所述第二翻转角配置。
在本发明中，所述第二参数包括与所述第一翻转角配置对应的回波链的等效回波时间。
在本发明中，所述第二准则包括第三阈值，所述第三阈值与所述第一翻转角配置对应的回波链的等效回波时间对应。
在本发明中，所述判断所述第一翻转角配置是否满足第二准则以获得第二判断结果包括确定第二信号演变，所述第二信号演变包括与所述第一翻转角配置对应的回波链中至少部分的回波的信号强度；基于所述第二信号演变，确定所述第一翻转角配置对应的回波链的等效回波时间；以及通过确定所述第一翻转角配置对应的回波链的等效回波时间不满足所述第二准则确定所述第一翻转角配置不满足所述第二准则。
在本发明中，所述判断所述第一翻转角配置是否满足所述第二准则以获得第二判断结果包括确定所述第一翻转角配置对应的回波链的等效回波时间小于所述第三阈值。
在本发明中，进一步包括确定所述初始翻转角配置满足所述第一准则；以及确定所述初始翻转角配置为所述第一翻转角配置。
根据本发明的一些实施例，进一步包括确定所述第一翻转角配置满足所述第二准则；以及确定所述第一翻转角配置为所述第二翻转角配置。
一种用于获取核磁共振信号的***。该***包括：至少一个处理器以及可执行指令。所述指令被所述至少一个处理器执行时，导致所述***实现一种方法。所述方法包括：获取回聚射频脉冲的初始翻转角配置，所述回聚射频脉冲用于产生回波链；判断所述初始翻转角配置是否满足第一准则，所述第一准则与所述初始翻转角配置对应的回波链的第一参数有关；基于所述第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；判断所述第一翻转角配置是否满足第二准则，所述第二准则与所述第一翻转角配置对应的回波链的第二参数有关；基于所述第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及基于所述第二翻转角配置，获取核磁共振信号。
本申请的一部分附加特性可以在下面的描述中进行说明。通过对以下描述和相应附图的检查或者对实施例的生产或操作的了解，本申请的一部分附加特性对于本领域技术人员是显而易见的。本披露的特性可以通过对以下描述的具体实施例的各种方面的方法、手段和组合的实践或使用得以实现和达到。
与现有技术相比，本发明的有益效果表现如下：
可以减小回聚射频脉冲的特殊吸收率，同时保持MR图像信噪比和对比度基本不变。
附图说明
根据示例性实施例可以进一步描述本申请。参考附图可以详细描述所述示例性实施例。所述实施例并非限制性的示例性实施例，其中相同的附图标记代表附图的几个视图中相似的结构，并且其中：
图1是根据本发明的一些实施例所示的核磁共振成像(MRI)***的工作环境的一个示意图；
图2是根据本发明的一些实施例所示的一个核磁共振扫描装置的结构示意图；
图3是根据本发明的一些实施例所示的一个计算设备的结构示意图；
图4是根据本发明的一些实施例所示的一个移动设备的结构示意图；
图5是根据本发明的一些实施例所示的处理引擎的一个示意图；
图6是根据本发明的一些实施例所示的处理模块的一个示意图；
图7是根据本发明的一些实施例所示的处理核磁共振(MR)信号的一种示例性流程图；
图8是根据本发明的一些实施例所示的获取核磁共振(MR)信号的一种示例性流程图；
图9是根据本发明的一些实施例所示的获取初始翻转角配置的一种示例性流程图；
图10是根据本发明的一些实施例所示的确定目标翻转角配置的一种示例性流程图；
图11是根据本发明的一些实施例所示的基于第一准则评价初始翻转角配置的一种示例性流程图；
图12是根据本发明的一些实施例所示的基于第二准则评价第一翻转角配置的一种示例性流程图；
图13是根据本发明的一些实施例所示的回波链与翻转角关系曲线的一种示例性示意图；
图14是根据本发明的一些实施例所示的基于固定翻转角配置获取的回波链重建的MR图像；
图15是根据本发明的一些实施例所示的基于可变翻转角配置获取的回波链重建的MR图像；
具体实施方式
在下面的详细描述中，为了提供对相关申请的透彻理解，通过实施例阐述了实施例的具体细节。对本申请的实施例的各种修改对于本领域技术人员将是显而易见的，并且在不脱离本申请的精神和范围的情况下，本申请定义的一般原理可以应用于其他实施例和应用场景中。因此，本申请不限于所示的实施例，而是符合与权利要求一致的最广范围。
本申请使用的术语仅用于描述特定示例性实施例的目的，而不是限制性的。如本申请所使用的，除非上下文明确提示例外情形，单数形态的“一”，“一个”和“该”也可以包括复数。应当进一步理解的是，在本说明书中使用的术语“包括”和 /或“包含”时，指存在所述特征、整形常量、步骤、操作、元素和/或组件，但不排除存在或添加一个或多个其他特征、整形常量、步骤、操作、元件、组件和/ 或其组合。
应当理解的是，本申请使用的术语“***”、“引擎”、“单元”、“模块”和/或“区块”是一种以升序的形式来区分不同组件、元件、部件、部件或组件的级别方法。然而，如果其他表达方式达到相同的目的，则这些术语可能被其他表达方式所取代。
通常，如本申请所使用的“模块”、“单元”或“区块”是指以硬件或固件或软件指令的集合体现的逻辑。本申请描述的模块、单元或区块可以在软件和/或硬件上被执行，并且可以被存储在任何类型的非暂时性计算机可读介质或其他存储设备中。在一些实施例中，软件模块、单元、区块可以被编译并连接到可执行程序中。应当理解的是，软件模块可以从其他模块、单元、区块或其自身调用和/ 或可以响应于检测到的事件或中断而被调用。配置用于在计算设备上执行的软件模块/单元/区块(例如，如图1所示的处理引擎120)可以被提供在计算机可读介质上，诸如光盘、数字视频盘、闪存驱动器、磁盘或任何其他有形介质或作为数字下载(并且可以原始地以压缩或可安装的格式存储，在执行之前需要安装、解压缩或解密)。软件代码可以部分或全部存储在执行计算设备的存储设备上供计算设备执行。软件指令可以被嵌入到固件当中，例如EPROM。应当理解的是，硬件模块、单元或区块可以包括在连接的逻辑组件中，例如门和触发器和/或可以包括在诸如可编程门阵列或处理器之类的可编程单元中。本申请描述的模块、单元、区块或计算设备功能可以被实现为软件模块/单元/区块，但是可以用硬件或固件来表示。通常，这里描述的模块、单元、区块是指可以与其他模块、单元、区块组合或者分成子模块、子单元、子区块的逻辑模块、单元、区块，尽管它们的物理组织或存储。所述描述可以适用于***、引擎或其一部分。
应当理解的是，当单元、引擎、模块或区块被称为“在…上”、“连接到”或“耦合到”另一单元、引擎、模块或区块时，其可以直接在另一单元、引擎、模块或区块通信、或者可以存在单元、引擎、模块或区块，除非上下文明确提示例外情形。如本申请所使用的，术语“和/或”包括一个或多个相关所列项目的任何一种和所有的组合。
参照附图并考虑以下描述，本申请的这些和其他特征以及相关的结构元件以及制造的部件和经济的结合的操作和功能的方法可以变得更加明显，且都构成本申请的一部分。然而，应当明确地理解，附图仅仅是为了说明和描述的目的，并不意图限制本申请的范围。应当理解附图不是按比例的。
图1是根据本发明的一些实施例所示的核磁共振成像(MRI)***的工作环境的一个示意图。如图1所示，MRI***100可以包括MR扫描仪110、处理引擎120、存储设备130、一个或多个终端140以及网络150。在一些实施例中，MR扫描仪110、处理引擎120、存储设备130和/或终端140之间的连接和/或通信方式可以是有线连接，无线连接(例如，网络150)以及以上两种组合。例如，MR扫描仪110可以通过网络150连接处理引擎120。又例如，MR 扫描仪110可以直接连接处理引擎120。
MR扫描仪110可以通过扫描目标或部分目标生成或提供与MR信号相关的图像数据。在一下实施例中，MR扫描仪110可以包括，例如主磁体220、梯度线圈230、射频线圈240等。关于MR扫描仪110的详细描述可以参见图 2。在一些实施例中，根据磁场的强度，MR扫描仪110可以包括高场强MR扫描仪、中场强MR扫描仪、低场强MR扫描仪等。
所述目标可以包括生物的或非生物的。例如，所述目标可以包括患者、人造物体等。又例如，所述目标可以包括患者的特定部分、器官和/或组织。例如，所述目标可包括头、脑、颈、体、肩、臂、胸、心、胃、血管、软组织、膝、脚等其中一种或几种的组合。在一些实施例中，MR扫描仪110可接收与至少部分的所述目标相关的MR信号。
处理引擎120可以处理从MR扫描仪110、存储设备130和/或终端140获得的数据和/或信息。例如，处理引擎120可以通过处理由MR扫描仪110采集的图像数据(例如，编码的MR信号)以生成MR图像。又例如，处理引擎 120可以根据，例如扫描协议，确定一个或多个成像参数(例如，与回聚射频脉冲信号相关的目标翻转角度配置)。在一些实施例中，处理引擎120可以是单个服务器或服务器组。服务器组可以是集中式的，也可以是分布式的。在一些实施例中，处理引擎120可以是本地的或远程的。例如，处理引擎120可以通过网络150从MR扫描仪110、存储设备130和/或终端140访问信息和/或数据。又例如，处理引擎120可以直接连接MR扫描器110、终端140和/或存储设备130以访问信息和/或数据。在一些实施例中，处理引擎120可以在云平台上实现。例如，云平台可以包括私有云、公共云、混合云、社区云、分布式云、跨云、多云等其中一种或几种的组合。在一些实施例中，处理引擎120可以由图3所描述的具有一个或多个组件的计算装置300实现。
存储设备130可以存储数据、指令和/或任何其他信息。在一些实施例中，存储设备130可以存储从处理引擎120和/或终端140获得的数据。在一些实施例中，存储设备130可以存储处理引擎120可以执行或使用的数据和/或指令，以执行本申请中描述的示例性方法。在一些实施例中，存储设备130可包括大容量存储器、可移除存储器、易失性读写存储器、只读存储器(ROM)等其中一种或几种的组合。大容量存储可以包括磁盘、光盘、固态硬盘、移动存储等。可移除存储器可以包括闪存驱动器、软盘、光盘、存储卡、ZIP磁盘、磁带等。易失性读写存储器可以包括随机存取存储器(RAM)。RAM可以包括动态随机存储器(DRAM)、双数据率同步动态随机存取存储器(DDR SDRAM)、静态随机存取存储器(SRAM)、可控硅随机存取存储器(t-ram)、零电容随机存取存储器(Z-RAM)等。ROM可以包括掩模只读存储器 (MROM)、可编程的只读存储器(PROM)、可擦除可编程只读存储器 (EPROM)，电可擦除可编程只读存储器(EEPROM)、光盘只读存储器(CD- ROM)、数字多功能光盘的光盘等。在一些实施例中，存储装置130可以通过本申请中描述的云平台实现。
在一些实施例中，存储设备130可以连接网络150，以与MRI***100中的一个或多个组件(例如，处理引擎120、终端140等)之间实现通信。MRI ***100中的一个或多个组件可以通过网络150读取存储设备130中的数据或指令。在一些实施例中，存储设备130可以是处理引擎120的一部分。
终端140可与MR扫描仪110、处理引擎120和/或存储装置130之间相互连接和/或通信。例如，处理引擎120可以从终端140获取扫描协议。又例如，终端140可以从MR扫描仪110、处理引擎120和/或存储设备130获取图像数据。在一些实施例中，终端140可以包括移动设备140-1、平板电脑140-2、笔记本电脑140-3等其中一种或几种的组合。例如，移动设备140-1可以包括移动电话、个人数字助理(PDA)、游戏设备、导航设备、销售点(POS)设备、笔记本电脑、平板电、台式机等其中一种或几种的组合。在一些实施例中，终端140可以包括输入设备、输出设备等。输入设备可以包括键盘、触摸屏(例如，触觉或触觉反馈)、语音输入、眼睛跟踪输入、大脑的监测***等任何其他类似的输入机制。通过输入设备接收的输入信息可通过，例如总线等传送到处理引擎120，以便进一步处理。其他类型的输入设备可以包括光标控制设备，如鼠标、轨迹球、光标方向键等。输出设备可以包括显示器、音箱、打印机等其中一种或几种的组合。在一些实施例中，终端140可以是处理引擎120 的一部分。
网络150可以包括能够促进磁共振成像***100的信息和/或数据交换的任何合适的网络。在一些实施例中，MRI***100的一个或多个组件(例如， MR扫描仪110、处理引擎120、存储设备130、终端140等)可以通过网络 150与MRI***100的一个或多个组件之间交换信息和/或数据。例如，处理引擎120可以通过网络150从MR扫描仪110获取图像数据(例如，MR信号)。又例如，处理引擎120可以通过网络150从终端140获取用户指令。网络150 可以包括公共网络(如互联网)、私人网络(例如，局域网(LAN)、广域网 (WAN))等)、有线网络(如以太网)、无线网络(例如，802.11网络、无线 Wi-Fi网络等)、蜂窝网络(例如，长期演进(LTE)网络)、帧中继网络、虚拟专用网络(“VPN”)、卫星网络、电话网络、路由器、集线器、服务器计算机等其中一种或几种组合。例如，网络150可以包括有线网络、光纤网络、电信网络、局域网、无线局域网(WLAN)、城域网(MAN)，公用电话交换网(PSTN)、蓝牙TM网络，ZigBeeTM网络、近场通信(NFC)网络等其中一种或几种的组合。在一些实施例中，网络150可以包括一个或多个网络接入点。例如，网络150可以包括有线和/或无线网络接入点，例如基站和/或因特网交换点，通过所述接入点，MRI***100的一个或多个组件可以连接网络150以交换数据和/或信息。
需要注意的是，以上对于MRI***100的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对各个模块进行任意组合或者构成子***与其他模块连接，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，存储设备130可以配置在具有数据存储功能的云计算平台中，包括但不限于公用云、私有云、社区云和混合云等。又如，成像设备110、数据处理引擎120、存储设备130和终端140等中的两个或多个可以不通过网络150通信，而直接配置在一个设备中。诸如此类的变形，均在本申请的保护范围之内。
图2是根据本发明的一些实施例所示的一个MR扫描仪110的模块示意图。如图2所示，MR扫描器110可以包括主磁体220、梯度线圈230和射频线圈240和脉冲序列250。
主磁体220可以在扫描至少部分的目标过程中产生静态磁场。主磁体112 可以包括各种类型的，例如，例如永久磁铁、超导电磁铁、电阻电磁铁等。
梯度线圈230可向x方向、y方向和/或z方向上的主磁场提供磁场梯度。这里所述的x方向、y方向和z方向可以表示坐标系中的x轴、y轴和z轴。例如，x轴和z轴可以在水平面上，x轴和y轴可以在垂直平面内，z轴可以沿着主磁体112的轴线。在一些实施例中，梯度线圈230可以包括X方向线圈，用于向X方向的主磁场提供磁场梯度，Y方向线圈，用于向Y方向的主磁场提供磁场梯度，以及Z方向线圈，用于向X方向的主磁场提供磁场梯度。在一些实施例中，X方向线圈、Y方向线圈和/或Z向线圈可以包括各种形状或结构。例如，Z方向线圈可以是循环(Maxwell)线圈。又例如，X方向和Y方向的线圈可以是马鞍状(Golay)线圈。
射频线圈240可以向被检测的人体210发射射频脉冲信号和/或从被检测的人体210接收MR信号。在一些实施例中，射频线圈240可以包括发射线圈和接收线圈。发射线圈可以发出信号(例如，射频脉冲)，所述信号可以激发目标(例如，人体210)中的原子核发生共振。接收线圈可以接收从所述目标发射的MR信号。在一些实施例中，射频发射线圈和射频接收线圈可以集成到同一线圈中。在一些实施例中，射频线圈240可以包括各种类型，例如，QD正交线圈、相阵列线圈等。
在一些实施例中，射频线圈240可以周期性地产生一系列射频脉冲。例如，射频线圈240可以产生一个激发射频脉冲和一定数量的回聚射频脉冲。在一些实施例中，激发射频脉冲和回聚射频脉冲可以包括一个或多个参数，例如带宽(也被称为频率范围)、振幅或强度、施加射频脉冲的时间、施加射频脉冲的持续时间、回聚时间(也可以称为两个射频脉冲之间的时间间隔)、射频脉冲翻转角、射频脉冲数量等。例如，射频线圈240可产生一个翻转角度为 90°的激发射频脉冲以及多个翻转角度为180°的回聚射频脉冲。在一些实施例中，激发射频脉冲可以具有90°以外的翻转角，例如，从0°到180°的任何角度。回聚射频脉冲翻转角可以是除了180°的其它角度。此外，射频线圈240可以周期性地产生一系列射频脉冲。例如，射频线圈240可以产生一个翻转角度为90°的激发射频脉冲以及多个具有相同的翻转角或可变的翻转角的回聚射频脉冲，所述可变翻转角的范围是0°到180°。激励射频脉冲的翻转角也可以是可变的。
在一些实施例中，激发射频脉冲可用于相应的梯度线圈230生成磁场梯度，以及多个回聚射频脉冲可以用于生成一个或多个MR信号(例如，具有多个回波的回波链)。在一些实施例中，与MR信号相关的参数可以包括，例如，MR信号类型(自旋回波、快速自旋回波(FSE)、快速恢复自旋回波、单发快速自旋回波、梯度回波、稳态进动快速成像等)、回波数、等效回波时间 (TE)、重复时间(TR)、回波链长度(ETL)、回波链的期相数、回波信号强度等。例如，回波链的期相可以指的是回声链的一节、一部分或片段，可以用于表明在回波信号的变化趋势。期相的数目和/或每个期相的回波数可以取决于临床需要。在一些实施例中，回波链可以被分为多个期相(例如，三个期相)。对应于所述回波链的多个回聚射频脉冲也可以被分为与回波链期相对应的多个期相。在一些实施例中，在每个期相的回聚射频脉冲可以单调变化，例如，递增或递减。回波链长度(ETL)可以指回波序列中的回波数。回波链长度(ETL)可以是固定的，也可以是变化的。例如，对于相同的组织进行成像，ETL可以是不变的。对于不同的组织，ETL可以是变化的。此外，即使对于相同的组织，ETL也可能是可变的。重复时间(TR)可以指施加两个连续激发射频脉冲的时间间隔。回波时间(TE)是指激发射频脉冲的中点与自旋回波产生的中点之间的时间间隔。本申请中所使用的“中点”可以是指与脉冲，例如，激发射频脉冲，回聚射频脉冲等相对应的回波信号强度到达最大值时的时间点。在一些实施例中，在快速自旋回波中，用于填充K空间中心线的回波可以确定图像的对比度。等效回波时间(等效TE或TEeff)，也可以称为有效回波时间可以指的是激发射频脉冲的中点与用于填充K空间中心的回波(也被称为中心回波)的中点之间的时间间隔。
脉冲序列250可以包括至少部分的成像参数以及所述至少部分的成像参数在时序上的排列。在一些实施例中，成像参数可以包括射频线圈220发射的射频脉冲相关参数(例如，激励次数(NEX)、带宽等)、梯度线圈230产生的梯度场相关参数，MR信号相关参数(例如，回波时间(TE)、回波链长度 (ETL)、自旋回波类型、相数等)。在一些实施例中，成像参数可以包括图像的对比度和/或比、感兴趣的目标、层厚、成像的类型(例如，T1加权成像、 T2加权成像、质子密度加权成像等)、T1、T2、采集时间(TA)、反转时间等。本申请中所述的T1(即纵向弛豫时间)是指纵向磁化矢量达到(1-1/e)或其最大值的63％左右所需的时间。T2(即横向弛豫时间)是指横向磁化矢量降到1/e或其最大值的37％左右所需的时间。需要注意的是，对于不同的物体 (例如，组织)，即使在同一磁场下，它们的T1和/或T2通常也不相同。还需要注意的是，在同一磁场下，同一物体的同一组织的T1和T2可能是不同的。在一些实施例中，脉冲序列250可以包括自旋回波序列、梯度回波序列、扩散序列、反转恢复序列等其中一种或几种的组合。例如，自旋回波序列可以包括快速自旋回波序列(FSE、TSE)、快速采集弛豫增强序列(RARE)、半傅里叶采集单次激发快速自旋回波序列(HASTE)、Turbo梯度自旋回波序列 (TGSE)等其中一种或几种的组合。
在一些实施例中，脉冲序列250可以连接和/或与处理引擎120通信。例如，在MRI扫描过程之前，可根据临床需求或扫描协议，处理引擎120可以设计和/或确定至少一个部分的脉冲序列250(例如，与射频脉冲相关的参数、与梯度场相关的参数等)。在MRI扫描过程中，MR扫描仪110可以基于脉冲序列250扫描物体(例如，人体210)。例如，射频线圈240可以基于脉冲序列 250中的射频脉冲的特定参数发射射频脉冲，并且可以根据脉冲序列250(例如，回波时间)接收MR信号。
需要注意的是，以上对于MR扫描仪110的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对各个模块进行任意组合或者构成子***与其他模块连接，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，脉冲序列250可以集成于处理引擎120中。诸如此类的变形，均在本申请的保护范围之内。
图3是根据本发明的一些实施例所示的可以实现处理引擎120的示例性计算设备300的硬件和/或软件组件的示例性示意图。如图3所示，计算装置300 可包括处理器310、存储器320、输入/输出(I/O)330和通信端口340。
处理器310可以执行计算机指令(例如，程序代码)并可以根据申请中描述的技术执行处理引擎120的功能。所述计算机指令可以用于执行本申请中描述的特定功能，所述计算机指令可以包括例如程序、对象、组件、数据结构、程序、模块和功能。例如，处理器310可以处理从MR扫描仪110、存储设备 130、终端140、和/或MR***100的任何其它组件获取的图像数据。在一些实施例中，处理器310可以包括一个或多个硬件处理器，例如微控制器、微处理器、精简指令集计算机(reduced instruction set computer(RISC))、特定应用集成电路(application specific integrated circuit(ASIC))、应用程序特定的指令集处理器(application-specific instruction-set processor(ASIP))、中央处理单元 (centralprocessing unit(CPU))、图形处理单元(graphics processing unit (GPU))、物理处理单元(physics processing unit(PPU))、数字信号处理器 (digital signal processor(DSP))、现场可编程门阵列(field programmable gate array(FPGA))、先进的RISC机器(advanced RISC machine(ARM))、可编程逻辑器件(programmable logic device(PLD))、能够执行一个或多个功能的任何电路或处理器等其中一种或几种的组合。
仅用于说明，在计算设备300中仅描述一个处理器。然而，需要说明的是，计算装置300也可以包括多个处理器。由本申请中描述一个处理器执行的操作和/或方法也可以由多个处理器共同或分别执行。例如，如果本申请中描述的计算设备300的处理器执行操作A和操作B，应当理解的是，操作A和操作 B也可以由计算装置中的300中的两个或两个以上不同处理器共同或分别执行 (例如，第一处理器执行操作A和第二处理器执行操作B，或第一处理器和第二处理器共同执行操作A和B)。
存储器320可以存储从MR扫描仪110、存储设备130、终端140、和/或 MRI***100的任何其它组件获取的数据/信息。在一些实施例中，存储器320 可包括大容量存储器、可移除存储器、易失性读写存储器、只读存储器(ROM)等其中一种或几种的组合。大容量存储可以包括磁盘、光盘、固态硬盘、移动存储等。可移除存储器可以包括闪存驱动器、软盘、光盘、存储卡、 ZIP磁盘、磁带等。易失性读写存储器可以包括随机存取存储器(RAM)。 RAM可以包括动态随机存储器(DRAM)、双数据率同步动态随机存取存储器(DDR SDRAM)、静态随机存取存储器(SRAM)、可控硅随机存取存储器(t- ram)、零电容随机存取存储器(Z-RAM)等。ROM可以包括掩模只读存储器 (MROM)、可编程的只读存储器(PROM)、可擦除可编程只读存储器(EPROM)，电可擦除可编程只读存储器(EEPROM)、光盘只读存储器(CD- ROM)、数字多功能光盘的光盘等。在一些实施例中，存储器320可以存储一个或多个程序和/或指令，用于执行本申请中描述的示例性方法。例如，存储 320可以存储程序，所述程序可以用于处理引擎120确定目标翻转角配置。
输入/输出330可以输入和/或输出信号、数据、信息等。在一些实施例中，输入/输出330可以实现用户与处理引擎120之间的交互。在一些实施例中，输入/输出330可以包括输入设备和输出设备。输入设备可以包括键盘、鼠标、触摸屏、麦克风等其中一种或几种的组合。输出装置可以包括显示装置、扬声器、打印机、投影仪等其中一种或几种的组合。所述显示装置可以包括液晶显示器(LCD)、发光二极管(LED)显示器、平板显示器、弧形屏幕、电视装置、阴极射线管(CRT)、触摸屏等其中一种或几种的组合。
通信端口340可以连接网络(例如，网络120)，以便于数据通信。通信端口340可以在处理引擎120和MR扫描仪110、终端140和/或存储设备130之间建立连接。所述连接可以是有线连接、无线连接、任何能够实现数据传输和/ 或接收的连接等其中一种或几种的组合。所述有线连接可以包括例如电缆、光缆、电话线等其中一种或几种的组合。所述无线连接可以包括，例如，蓝牙TM链接、Wi-FiTM链接、WiMAXTM链路、无线局域网链接、ZigBee链接、移动网络链接(例如，3G、4G、5G、等)其中一种或几种的组合。在一些实施例中，通信端口340可以是和/或包括标准化通信端口，如RS232、RS485等。在一些实施例中，通信端口340可以是专门设计的通信端口。例如，通信端口 340可以根据医学(DICOM)协议中的数字成像和通信来设计。
图4是根据本发明的一些实施例所示的可以实现终端140的示例性移动设备400的硬件和/或软件组件的示例性示意图。如图4所示，移动设备400可包括通信平台410、显示器420、图形处理器(GPU)430、中央处理器(CPU) 440、输入/输出450、内存460和存储器490。在一些实施例中，移动设备 400也可以包括任何其它合适的组件，包括但不限于***总线或控制器(图中未显示)。在一些实施例中，移动操作***470(例如，iOSTMTM，Android，WindowsPhoneTM等)和一个或多个应用程序480可以从存储器490装载入内存 460，以便能够由中央处理器440执行。应用程序480可以包括浏览器或任何其他合适的移动应用程序，用于从处理引擎120接收和呈现与图像处理或其他信息相关的信息。信息流的用户交互可以通过输入/输出450实现，并且通过网络150提供给处理引擎120和/或MR***100的其他组件
为了实现本申请中描述的各种模块、单元及其功能，计算机硬件平台可以用作本申请中描述的一个或多个元素的硬件平台。具有用户界面元素的计算机可用于实现个人计算机(PC)或任何其他类型的工作站或终端设备。如果适当编程，计算机也可以充当服务器。
图5是根据本发明的一些实施例所示的处理引擎120的模块示意图。处理引擎120可以包括获取模块502、控制模块504、处理模块506和存储模块 508。如图3所示，处理引擎120的至少一部分可以通过图3所示的计算设备或图4所示的移动设备执行。
获取模块502可以获取数据。在一些实施例中，获取模块502可以从MR 扫描仪110、存储设备130和/或终端140获取数据。在一些实施例中，所述数据可以包括扫描协议、至少部分的成像参数、图像数据(例如，MR信号、MR 图像等)、指令等其中一种或几种的组合。所述指令可以由处理引擎120的处理器执行，以运行本申请中描述的示例性方法。在一些实施例中，所获取的数据可以被传送到处理模块506以进行进一步的处理或存储在存储模块508中。
控制模块504可以控制获取模块502、存储模块508和/或处理模块506的操作(例如，通过生成一个或多个控制参数)。例如，控制模块504可以控制处理模块506确定回聚射频脉冲的翻转角配置。又例如，控制模块504可以控制获取模块502获取图像数据(例如，MR信号)。再例如，控制模块504可以控制处理模块506通过处理获取的MR信号来生成MR图像。在一些实施例中，控制模块504可接收实时命令或检索由用户(例如，医生)提供的预定指令，以控制获取模块502和/或处理模块506的一个或多个操作。例如，控制模块504可以调整获取模块502和/或处理模块506，以根据所述实时指令和/或所述预定指令生成与所述MR信号相关联的图像数据。在一些实施例中，控制模块504可以与处理引擎120的一个或多个其它模块通信，以交换信息和/或数据。
处理模块506可以处理处理引擎120的各个模块提供的数据。在一些实施例中，处理模块506可以处理获取模块502获取的或从存储模块508和/或存储设备130中检索的MR信号以重建MR图像。在一些实施例中，处理模块506 可以确定和/或调整至少部分的成像参数，例如回聚射频脉冲的翻转角配置。
存储模块508可以存储信息。所述信息可以包括程序、软件、算法、数据、文本、数字、图像和其他一些信息。例如，所述信息可以包括扫描协议、成像参数、图像数据(例如MR信号、MR图像等)等其中一种或几种的组合。在一些实施例中，存储模块508可以存储由处理引擎120的处理器执行的程序和/或指令，以获取数据、确定成像参数、重建MR图像和/或显示任何中间结果或结果图像。
在一些实施例中，图5中所示的一个或多个模块可以通过如图1所示的 MRI***100中的至少部分组件实施。例如，获取模块502、控制模块504、处理模块506和/或存储模块508可以集成为控制台(图中未示出)。通过所述控制台，用户可以设置用于扫描物体的参数、控制成像过程的参数、控制图像重建的参数等。在一些实施例中，控制台可以通过处理引擎120和/或终端140 来实施。
图6是根据本发明的一些实施例所示的处理模块506的示例性示意图。如图所示，处理模块506可包括初始化单元602、评估单元604、调整单元606、计算单元608和存储单元610。处理模块506可以通过各种组件(例如，图3 所示的计算设备300的处理器310)实现。例如，处理模块506的至少部分单元可以在图3所示的计算设备或图4所示的移动设备上实现。
初始化单元602可以初始化本申请中其他地方所描述的一个或多个成像参数。例如，初始化单元602可以提供的初始翻转角配置，所述初始翻转角配置可以包括多个初始翻转角度值。在一些实施例中，初始化单元602可以初始化与初始翻转角度配置有关的一个或多个成像参数(例如，在一个回波链的期相中的起始翻转角和终点翻转角)。在一些实施例中，由初始化单元602初始化的成像参数可以被发送至处理模块506中的其他单元以进行进一步处理。例如，初始化单元602可以将初始翻转角配置发送到存储单元510以进行存储。
评估单元604可以在确定目标翻转角配置的过程中执行评估功能。例如，评估单元604可以确定翻转角配置(例如，由初始化单元602确定的初始翻转角配置)是否满足准则、标准、要求或阈值。在一些实施例中，评估单元604 可以将评估结果传送至处理模块506中的其他单元以进行进一步处理。例如，评估单元604可以将评估结果传送至调整单元606以调整翻转角配置。
调整单元606可调整至少部分的成像参数。例如，调整单元606可以调整至少部分的翻转角度配置(例如，由初始化单元602确定的初始翻转角配置或中间翻转角度配置)。在一些实施例中，调整单元606可以向处理模块506中的其他单元发送调整后的翻转角度配置，以便进行进一步处理。例如，调整单元606可以向存储单元510发送调整后的翻转角度配置以存储所述翻转角配置。
计算单元608可以计算从初始化单元602、评估单元604、调整单元606 和/或MRI***100中的其他模块或单元接收的不同类型的信息。例如，计算单元608可以通过利用不同的函数和/或基于初始信息或预设条件的算法计算初始翻转角配置中的翻转角度。所述函数可以包括布洛赫(Bloch)方程、多项式函数、线性函数、三角函数、反三角函数、指数函数、幂函数、对数函数等其中一种或几种的组合。所述算可以包括递归法、平分法、穷举搜索(或蛮力搜索)法、贪心算法、分治算法、动态规划法、迭代法、分支定界算法、回溯算法等其中一种或几种的组合。
存储单元610可以存储用于确定目标翻转角配置的信息。所述信息可以包括程序、软件、算法、数据、文本、数字和其他一些信息。例如，存储单元 610可存储由初始化单元602确定的初始翻转角配置、由计算单元608和/或调整单元608生成的中间翻转角配置和/或目标翻转角配置。在一些实施例中，存储单元610可以存储用于评估初始翻转角配置和/或中间翻转角配置的要求、准则、阈值或标准。
需要注意的是，以上对于处理模块506的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对各个模块进行任意组合或者构成子***与其他模块连接，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，评估单元620和调整单元640可以集成到一个单元中。又例如，计算单元608可以分别集成到初始化单元602、评估单元604和/或调整单元606中。诸如此类的变形，均在本申请的保护范围之内。
图7是根据本发明的一些实施例所示用于处理MR信号的流程700的一种示例性流程图。在一些实施例中，图7所示的用于处理MR信号的流程700中的一个或多个操作可以通过图1所示的MRI***100实现。例如，图7所示的过程可以以指令的形式存储于存储设备130中，并由处理引擎120执行调用和/ 或执行(例如，图3所示的计算设备300的处理器310、图4所示的移动设备 400的中央处理器440)。
在702中，可以获取MR信号。操作702可以由获取模块502执行。在一些实施例中，所述MR信号可以从MR扫描仪110、存储设备130、终端140 和/或外部数据源获取。例如，所述MR信号可以由MR扫描仪110通过向待检查物体施加激发射频脉冲和多个回聚射频脉冲产生。在一些实施例中，所述多个回聚射频脉冲可以包括目标翻转角配置，所述目标偏转角配置的确定可以参考本申请图8、10、11、和/或12中的详细描述。在一些实施例中，所述MR 信号可以包括一个或多个回波链，所述回波链可以包括多个回波。所述回波的类型可以包括自旋回波、快速自旋回波(FSE)，快速恢复自旋回波、单发快速自旋回波、梯度回波、稳态进动快速成像等，在一些实施例中，所述MR信号可以以数值的形式存储于存储设备130和/或存储模块508。在一些实施例中，所述MR信号可以表示为K空间的数据值(也被称为频率域)。
在704中，可以通过处理MR信号来生成图像。操作704可以由处理模块 506执行。在一些实施例中，可以基于图像重建技术来处理所述MR信号以生成所述图像。所述图像重建技术可以包括傅立叶变换(FT)、频率编码、相位编码、迭代重建、后投影重建等其中一种或几种的组合。例如，所述傅里叶变换可以包括快速傅里叶变换(FFT)、二维傅里叶变换、三维离散傅里叶变换、逆傅里叶变换、反向快速傅里叶变换(IFFT)等其中一种或几种的组合。
在706中，可以输出所述图像。操作706可以由处理模块506执行。在一些实施例中，所述图像可以输出至终端140以显示所述图像。在一些实施例中，所述图像可以输出至存储设备130和/或存储模块508以存储所述图像。
需要注意的是，以上对于处理MR信号的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，处理MR信号的流程可以进一步包括预处理所述MR信号。诸如此类的变形，均在本申请的保护范围之内。
图8是根据本申请的一些实施例所示用于获取MR信号的流程800的一种示例性流程图。图7所示的操作702可以根据图8所示的流程800执行。在一些实施例中，图8所示的用于获取MR信号的流程800中的一个或多个操作可以通过图1所示的MRI***100实现。例如，图8所示的流程800可以以指令的形式存储于存储设备130中，并由处理引擎120执行调用和/或执行(例如，图3所示的计算设备300的处理器310、图4所示的移动设备400的中央处理器 440)。
在802中，可以提供与回聚射频脉冲相关的初始翻转角配置。操作802可以由初始化单元602执行。在一些实施例中，初始翻转角度调度可以包括特定数目的初始翻转角，所述特定数目的初始翻转角可以包括相同的初始值或不同的初始值。在一些实施例中，所述初始翻转角配置可以分为一个或多个部分，所述一个或多个部分可以对应于回波链的一个或多个期相。这里所述的回波链的期相可以指回波链的趋势。例如，所述回波链可以分为三个期相，包括第一期相、第二期相和第三期相。在每一个期相中，回波链中的回波信号的变化趋势可以不同于其他两个期相。例如，第一期相中回波信号的趋势可以是稳定状态，表明第一期相中的回波信号没有显著变化。因此，所述期相可以表明回聚射频脉冲的翻转角被划分的组数或期相。例如，所述初始翻转角配置可以分为三个部分，所述三个部分分别与回波链的所述第一期相、所述第二期相以及所述第三期相对应。在每一个期相中，初始翻转角可以单调变化。例如，对应于所述第一期相的初始翻转角度可以是单调递减的，对应于所述第二期相的初始翻转角可以单调递增，以及对应第三期相的初始翻转角可以单调递减。在一些实施例中，可以通过确定每个期相的翻转角度来确定初始翻转角配置。在一些实施例中，不同期相的成像参数可以不同或相同的。例如，同一物体中相同组织产生的回波链，不同期相对应的T1或T2值可以是相同的。在一些实施例中，可以通过基于至少部分的成像参数，利用函数和/或算法确定所述初始翻转角配置。例如，所述函数可以包括布洛赫(Bloch)方程、多项式函数、线性函数、三角函数、反三角函数、指数函数、幂函数、对数函数等其中一种或几种的组合。所述算法可以包括递归法、平分法、穷举搜索(或蛮力搜索)、贪心算法、分治算法、动态规划法、迭代法、分支定界算法、回溯算法等其中一种或几种的组合。在一些实施例中，所述初始翻转角配置可以由用户通过终端 140设置，也可以是MRI***100的默认设置。例如，所述初始翻转角配置中的翻转角度可以具有合适的数值，所述合适的数值可以在特定数值范围内基于特定规则变化，例如在90°到180°范围内单调递增。
在804中，可以评估所述初始翻转角配置。操作804可以由评估单元604 执行。在一些实施例中，可以根据与所述初始翻转角配置相对应的信号演化来评估所述初始翻转角配置。所述信号演化可以包括具有多个回波信号的回波链，所述回波信号可以基于具有所述初始翻转角配置的多个回聚射频脉冲获得。在一些实施例中，所述评估可以基于所述信号演化中的一个或多个回波信号是否满足准则、标准、要求或阈值。在一些实施例中，所述评估可以基于所述信号演化的部分(例如一个期相或一片段)是否满足要求、标准、阈值或标准。在一些实施例中，所述评估可以基于所述信号演化与参考信号演化之间的比较。所述信号演化和所述参考信号演化可以分别对应一个回波链，所述回波链可以分为一个或多个期相。所述参考信号演化可以对应一个或多个回波链，所述回波链中的每个回波可以具有期望的信号强度。另外，所述参考信号演化可以包括限制一个或多个回波信号强度，例如，每个期相中起始回波的信号强度、末端回波的信号强度、回波链中第n个回波的信号强度等。又例如，所述参考信号演化可以指定期望的信号强度，而不指定回波链中相应的回波(对应于指定信号强度的回波)如何或何时产生。在一些实施例中，所述参考信号演化可以包括限制一个或多个回波链中的每一个期相的趋势。所述趋势可以是增加、减小、稳定状态等其中一种或几种的组合。需要注意的是，所述趋势的限制可以应用于回声链中的任何期相。
在一些实施例中，可以根据翻转角度和回波信号之间的关系来确定所述信号演化。例如，可以通过例如布洛赫(Bloch)方程、EPG算法等其中一种或几种的组合描述所述关系。在一些实施例中，可以根据被检查目标的组织的T1 和T2来确定所述信号演化。可选的，所述信号演化的计算可以不考虑组织的松弛时间。在一些实施例中，所述参考信号演化可以由用户选择或设置。在一些实施例中，所述参考信号演化可以由MRI***100根据用户提供的信息确定。所述信息可以包括要执行的实际成像、需要成像的目标、所述目标的T1 值、所述目标的T2值、所述目标的质子密度等其中一种或几种的组合。关于信号演化的确定方法，可以参考文献，例如，国际专利申请 PCT/CN2015/087818，申请日为2015/08/21，标题为SYSTEM AND METHOD FOR FLIP ANGLE DETERMINATION IN MAGNETIC RESONANCE IMAGING，文献中的内容以引用方式被包含于此。
在一些实施例中，可以基于一个或多个准则来评估所述初始翻转角配置。例如，所述初始翻转角配置的评估可以包括确定所述初始翻转角配置是否满足第一准则和/或第二准则。如果初始翻转角配置满足第一准则，则可以根据第二准则评估所述初始翻转角配置。如果所述初始翻转角配置不满足所述第一准则，则可以调整所述初始翻转角配置以获得满足所述第一准则的中间翻转角配置。然后，可以根据所述第二准则评估所述中间翻转角配置。在一些实施例中，所述初始翻转角度的评估可以包括基于所述第一准则和/或第二准则确定和 /或评估一个或多个中间翻转角配置。
在806中，可以基于所述初始翻转角配置的评估结果来确定目标翻转角配置。操作806可以由调整单元606执行。在一些实施例中，所述目标翻转角配置的确定可以根据一个或多个准则反复调整所述初始翻转角配置(例如，第一个标准，第二准则等)。例如，如果所述初始翻转角配置不满足第一准则，则可以调整所述初始翻转角配置以获得中间翻转角配置(例如，第一翻转角配置)。如果所述中间翻转角配置不满足第二准则，则可以调整所述中间翻转角配置以获得另一中间翻转角配置(例如，第二翻转角配置)。在一些实施例中，可以生成和调整多个中间翻转角配置，以得到目标翻转角配置，所述目标翻转配置可以同时满足第一准则和第二准则。
在808中，可以基于所述目标翻转角配置使得射频线圈240产生回聚射频脉冲。操作808可以由控制模块504执行。在一些实施例中，可以基于所述目标翻转角配置确定一个或多个与所述回聚射频脉冲的参数。所述回聚射频脉冲的参数可以包括回聚射频脉冲的强度、回聚间隔(即两个回聚射频脉冲之间的时间间隔)、回聚射频脉冲数目等。例如，回聚射频脉冲的翻转角度值越大，所述回聚射频脉冲的强度越大。在一些实施例中，控制模块504可以控制射频线圈240基于所述回聚射频脉冲的参数产生回聚射频脉冲。
在810中，可以基于所述回聚射频脉冲获取MR信号。操作810可以由获取模块502执行。在一些实施例中，所述MR信号可以从MR扫描仪110中获取，所述MR信号可以由MR扫描仪110通过具有所述目标翻转角配置的回聚射频脉冲扫描被检查物体产生。在一些实施例中，所述MR信号可以从存储装置130、存储模块508和/或其它外部存储器获取。在一些实施例中，所述MR 信号可以包括自旋回波、快速自旋回波、快速恢复自旋回波、单发自旋回波、梯度回波等。在一些实施例中，所述MR信号可以包括回波链，所述回波链可以包括多个回波信号。每个回波可以对应一个具有目标翻转角的回聚射频脉冲。在一些实施例中，可以基于所述目标翻转转角配置和/或所述回聚射频脉冲的参数确定所述MR信号的参数。例如，回波的信号强度可以与回聚射频脉冲的翻转角有关。又例如，回波时间可以与两个回聚射频脉冲之间的时间间隔 (即回聚时间)有关。
本实施例提供一种获取核磁共振信号的方法，本方法通过基于一个或多个准则确定回聚射频脉冲的可变翻转角配置，以获取核磁共振信号，所述准则与回聚射频脉冲的总能量、回波信号强度、等效回波时间等参数有关，可以确保降低回聚射频脉冲的特殊吸收率(SAR)的情况下保持信噪比以及图像对比度基本不变。
需要注意的是，以上对于获取MR信号的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，操作802和804 可以同时执行。在一些实施例中，期相的数目和长度可以是可变的或可调整的。诸如此类的变形，均在本申请的保护范围之内。
图9是根据本申请的一些实施例所示用于确定初始翻转角配置的流程900 的一种示例性流程图。图8所示的操作802可以根据图9所示的流程900执行。在一些实施例中，流程900中的一个或多个操作可以通过图1所示的MRI ***100实现。例如，图9所示的流程900可以以指令的形式存储于存储设备 130中，并由处理引擎120执行调用和/或执行(例如，图3所示的计算设备 300的处理器310、图4所示的移动设备400的中央处理器440)。
在902中，可以确定与回波链相关的一个或多个参数。操作902可以由初始化单元602执行。所述与回波链相关的参数可以包括在本申请的其他地方所描述的至少部分的成像参数。例如，所述成像参数可以包括图像的对比度和/或比、感兴趣目标、层厚、成像的类型(例如，T1加权成像、T2加权成像、质子密度加权成像等)，T1、T2、自旋回波类型(例如，自旋回波快速自旋回波、快速恢复自旋回波、单发自旋回波、梯度回波、与稳态进动快速成像等)、翻转角度值、采集时间(TA)、回波时间(TE)、重复时间(TR)、回波链长度(ETL)、相数、激励次数(NEX)、反转时间、带宽等其中一种或多种的组合。
在一些实施例，可以确定初始条件。所述初始条件可以包括设置所述初始翻转角配置中部分翻转角的初始值。在一些实施例中，所述初始条件可以包括起始翻转角度和/或末端翻转角度的范围，并且所述起始翻转角度和末端翻转角度的大小可以在一定范围内。本申请中所使用的起始翻转角指的是MR信号 (例如，参考信号演变、回波链等)的一个期相的起点对应的翻转角度；末端翻转角指的是MR信号(例如，参考信号演变、回波链等)的一个期相的结束点对应的翻转角。在一些实施例中，所述初始条件(例如，起始翻转角和/或末端翻转角)可由用户通过终端140确定。在一些实施例中，所述初始条件(例如，起始翻转角和/或末端翻转角)可由处理引擎120根据临床需要确定。例如，在颈椎成像或膝关节成像中，回波链对应的最大翻转角可以设置为低于 120°。又例如，在腹部成像中，最小翻转角度可以大于80°。又例如，在髋关节成像中，最大翻转角可设置为140°。在一些实施例中，在每个期相中，起始翻转角可以是最大翻转角度，并且结束翻转角可以是最小翻转角度。
在904中，可以通过处理所述参数确定函数。操作902可以由初始化单元 602执行。所述函数可用于描述初始翻转角配置或其部分(例如，初始翻转角配置的一个期相)。在一些实施例中，所述函数可以包括布洛赫方程、EPG算法、多项式、线性函数、三角函数、反三角函数、指数函数、幂函数、对数函数等其中一种或几种的组合。在一些实施例中，可以根据一个或多个参数，例如，起始翻转角、末端翻转角和/或特征参数来确定与初始翻转角配置对应的函数。本申请中描述的特征函数可以控制或影响翻转角度变化率。例如，回波链可分为一个或多个期相。每个期相的翻转角可以按照不同的指数函数变化。例如，假设N是回波链长度(ETL)，α0是一个期相的起始翻转角，αN-1是所述期相的末端翻转角，所述期相内的翻转角可以用下面的函数来描述：如果αN-1≥α0，可以通过如下描述的方程(1)计算所述期相内的其余翻转角度：
Figure BDA0001391490010000311
如果α_0>α_(N-1),可以通过如下描述的方程(2)计算所述期相内的其余翻转角度：
Figure BDA0001391490010000321
其中n＝0，1，…，n-1，p是特征参数。在一些实施例中，p可以控制或影响一个期相内的起始点和结束点周围的翻转角度的变化率。p可能是大于1 的实数。不同期相对应的P值可以相同或不同。
在一些实施例中，回波链可以分为一个或多个期相。在一个期相内的回聚射频脉冲的翻转角可以根据线性函数变化。例如，假设N是回波链长度 (ETL)，α0是一个期相内的起始翻转角，αN-1是所述期相内的末端翻转角，所述期相内的翻转角可以根据如下描述的方程(3)计算：
Figure BDA0001391490010000322
其中n＝0，1，…，N-1。
在一些实施例中，回波链可以分为一个或多个期相。在一个期相内的回聚射频脉冲的翻转角可以根据多项式函数变化。例如，假设N是回波链长度 (ETL)，α0是一个期相内的起始翻转角，αN-1是所述期相内的末端翻转角，所述期相内的翻转角可以根据如下描述的方程(4)计算：
Figure BDA0001391490010000323
其中n＝0,1,…,N-1.P＝[P1，P2，…，PK]是特征参数向量。在一些实施例中，k可以是小于10的整数。P1，P2，…，PK可以基于翻转角度单调增加或单调递减确定，并符合以下标准：
Figure BDA0001391490010000324
在906中，可以根据所述函数来确定初始角度调度中的一个或多个翻转角度。操作902可以由初始化单元602执行。在一些实施例中，可以根据所述函数(例如方程(1)、方程(2)、方程(3)、方程(4)和/或方程(5)来计算一个期相中的其余翻转角度。关于初始翻转角配置的确定方法，可以参考文献，例如，国际专利申请PCT/CN2015/087818，申请日为2015/08/21，标题为 SYSTEM AND METHOD FOR FLIP ANGLE DETERMINATION IN MAGNETICRESONANCE IMAGING，文献中的内容以引用方式被包含于此。
本实施例提供一种获取确定初始翻转角配置的方法，本方法基于回波信号演化与翻转角配置之间的关系确定所述初始翻转角配置。由于组织的T2衰减特性，使得回波链上的信号强度是依次衰减的，因此可以通过逐渐增加翻转角大小，以确保信号强度保持基本不变，增加图像信噪比。
需要注意的是，以上对于确定初始翻转角配置的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，操作904 可以不是必须的。在一些实施例中，所述函数可以包括其它的形式。诸如此类的变形，均在本申请的保护范围之内。
图10是根据本申请的一些实施例所示用于确定目标翻转角配置的流程 1000的一种示例性流程图。图8所示的操作804可以根据流程1000执行。在一些实施例中，流程1000中的一个或多个操作可以通过图1所示的MRI*** 100实现。例如，流程1000可以以指令的形式存储于存储设备130中，并由处理引擎120执行调用和/或执行(例如，图3所示的计算设备300的处理器 310、图4所示的移动设备400的中央处理器440)。
在1002中，可以确定回聚射频脉冲的初始翻转角配置。操作1002可以由初始化单元602执行。在一些实施例中，可以根据图8和9中的描述确定所述初始翻转角配置
在1004中，可以判断所述初始翻转角配置是否满足第一准则。操作1004 可以由评估单元604执行。如果所述初始翻转角度时间表满足第一标准，则流程1000可以进行操作1008。如果所述初始翻转角配置不满足所述第一准则，则流程1000可以执行操作1006。在一些实施例中，所述第一准则可以包括第一阈值。所述第一阈值可以与回波链有关，所述回波链可以基于所述初始翻转角配置下的多个回聚射频脉冲获得。例如，所述第一参数可以包括所述回波链中回波的信号强度。在一些实施例中，如果所述第一参数小于所述第一阈值，则所述初始翻转角配置不满足所述第一准则。如果所述第一参数等于或大于所述第一个阈值，则所述初始翻转配置不满足所述第一准则。在一些实施例中，如果所述第一参数和所述第一阈值之间的差值大于第四阈值(例如，常数)，则所述初始翻转角配置不满足所述第一准则。如果所述第一参数和所述第一阈值之间的差值小于或等于所述第四阈值(例如，常数)，则所述初始翻转角配置满足所述第一准则。
在一些实施例中，所述第一准则可以包括第二阈值。所述第二阈值与所述初始翻转角配置下的多个回聚射频脉冲有关。例如，所述第二参数可以包括所述初始翻转角配置下的多个回聚射频脉冲的总能量。在一些实施例中，如果所述第二参数大于所述第二阈值，则所述初始翻转角配置不满足所述第一准则。如果所述第二参数小于或等于所述第二阈值，则所述初始翻转角配置满足第一准则。在一些实施例中，如果所述第二参数和所述第二阈值之间的差值小于第五阈值(例如，常数)，则所述初始翻转角配置不满足所述第一准则。如果所述第二参数和所述第二阈值之间的差值等于或大于所述第五阈值(例如，常数)，则所述初始翻转角配置满足所述第一准则。在一些实施例中，如果所述第一参数不满足所述第一阈值或所述第二参数不满足所述第二阈值，则所述初始翻转角配置不满足所述第一准则。在一些实施例中，当所述第一参数满足所述第一阈值且所述第二参数满足所述第二阈值，则所述初始翻转角配置满足所述第一准则。在一些实施例中，所述第一阈值和所述第二阈值可以由用户根据临床需要(例如，组织的类型)设置。在一些实施例中，所述第一阈值和所述第二阈值可以由处理引擎120根据参考信号演化确定，如图8和/或图11中所述。本申请中描述的第四阈值和第五阈值可以是相对较小的值(例如，取值范围从0到1)。在一些实施例中，所述第四阈值和所述第五阈值可以由用户设置或根据MRI***100的默认设置而设定。
在一些实施例中，可以根据与所述初始翻转角配置对应的第一信号演化确定所述初始翻转角配置是否满足所述第一准则。所述第一信号演化可以包括具有多个回波的回波链，所述回波链可以基于所述初始翻转角配置下的多个回聚射频脉冲获得。在一些实施例中，可以基于所述第一信号演化中的一个或多个信号是否满足所述第一准则确定所述初始翻转角配置是否满足所述第一准则。
在1006中，可以调整所述初始翻转角配置。操作1006可以由调整单元 606执行。在一些实施例中，可以调整所述初始翻转角配置中的一个或多个翻转角度(例如，增加或减少)。在一些实施例中，可以基于算法，通过调整与所述初始角配置相对应的回波链的一个或多个参数(例如，起始翻转角、结束翻转角、特征参数等)来调整所述初始翻转角配置。典型的算法可以包括对分算法、穷举搜索(或蛮力搜索)算法、贪婪算法、分治算法、动态规划算法、迭代算法、分支定界算法、回溯算法等其中一种或几种的组合。在一些实施例中，可以根据操作1004中描述的，基于所述第一准则评估所述调整后的初始翻转角配置。
在1008中，可以确定第一个翻转角配置。操作1008可以由评估单元604 执行。所述第一个翻转角配置满足所述第一准则。例如，基于所述第一翻转角配置下的多个回聚射频脉冲获得的回波链相关的第一参数满足所述第一阈值且与所述回聚射频脉冲相关的所述第二参数满足所述第二阈值。在一些实施例中，当所述初始翻转角配置满足所述第一准则，则所述初始翻转角配置可以被指定为所述第一翻转角配置。
在1010中，可以确定所述第一翻转角配置是否满足所述第二准则。操作 1010可以由评估单元604执行。当所述第一翻转角配置满足所述第二准则，流程1000可以执行操作1018。当所述第一翻转角配置不满足所述第二准则，流程1000可以执行操作1012。在一些实施例中，所述第二准则可以包括第三阈值。所述第三阈值可以与回波链的第三参数相关。在一些实施例中，所述第三参数可以包括基于所述第一翻转角配置下的回聚射频脉冲产生的回波链的等效回波时间(TE)。在一些实施例中，当所述回波链的等效TE小于所述第三阈值，则所述第一翻转角配置不满足所述第二准则。当所述回波链的等效TE等于或大于所述第三阈值，所述第一翻转角配置满足所述第二准则。在一些实施例中，当所述回波链的等效TE与所述第三阈值之间的差值大于第六阈值(例如，常数)，则所述第一翻转角配置不满足所述第二准则。当所述第一回波链的等效TE与所述第三阈值之间的差值小于或等于第六阈值(例如，常数)，所述第一翻转角配置满足所述第二准则。在一些实施例中，可以根据参考信号演化确定所述第三阈值，如图8和/或11中所述。本申请中所描述的第六阈值可以是相对较小的值(例如，从0到1的范围内的常数)。在一些实施例中，所述第六阈值可以由用户设置，也可以根据MRI***100的默认设置来设置。
在一些实施例中，可以基于与所述第一翻转角配置相对应的第二信号演化确定所述第一翻转角配置是否满足第二准则。所述第二信号演化可以包括由一个或多个回波信号的回波链，所述回波信号可以基于所述第一翻转角配置下的回聚射频脉冲获得。在一些实施例中，可以基于所述第二信号演化的中心回波确定所述第一翻转角配置是否满足所述第二准则。例如，可以基于所述第二信号演化确定等效TE，基于所述等效TE确定所述第一翻转角配置是否满足所述第二准则。
在1012中，可以调整第一翻转角配置。操作1012可以由调整单元606执行。在一些实施例中，可以增加或减小所述第一翻转角配置中的一定数量的翻转角度以获得第二翻转角配置。可以根据所述第二准则确定需要调整的翻转角数量，如图12所述。在一些实施例中，可以通过调整与所述第一翻转角配置相对应的回波链的一个或多个参数(例如，开始翻转角、结束翻转角、特征参数)来调整所述第一翻转角配置。典型的算法包括递归法、平分法、穷举搜索 (或蛮力搜索)、贪心算法、分治算法、动态规划法、迭代算法、分支定界算法、回溯算法等其中一种或几种的组合。在一些实施例中，可以根据操作1010 中所述，基于所述第二准则评估调整后的第一翻转角配置。
在1014中，可以确定第二翻转角配置。操作1014可以由评估单元604执行。所述第二个翻转角配置满足所述第二准则。例如，基于所述第二翻转角配置的多个回聚射频脉冲而获得的回波链的等效TE满足所述第三阈值(例如，期望的等效TE)。
在1016中，所述第二翻转角配置可以被智能为初始翻转角配置。操作 1014可以由初始化单元602执行。在一些实施例中，可以根据操作10014- 1018，基于所述第一准则和/或所述第二准则评估所述第二翻转角配置。
在1018中，可以确定目标翻转角配置。操作1018可以由调整单元606执行。在一些实施例中，操作1004-1016可以被反复执行直到确定目标翻转角配置。所述目标翻转角配置可以同时满足所述第一准则和所述第二准则。在一些实施例中，当所述第一翻转角配置满足所述第二准则，所述第一翻转角配置可以被确定为所述目标翻转角配置。在一些实施例中，当所述第二翻转角配置满足所述第一准则，所述第二翻转角配置可以被确定为所述目标翻转角配置。
本实施例提供一种确定目标翻转角配置的方法，通过第一准则及第二准则调整所述目标翻转角配置中的翻转角。所述第一准则是基于回聚射频脉冲的总能量以及回波信号强度调整初始翻转角配置，在保证回聚射频脉冲总能量符合一定条件的情况下保证回聚射频脉冲的强度符合一定条件，以确保降低回聚射频脉冲的特殊吸收率(SAR)的情况下信噪比保持基本不变。所述第二准则是基于等效回波时间确定目标翻转角，以确保图像对比度保持基本不变。
需要注意的是，以上对于确定目标翻转角配置的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，操作1002 不是必须的。又例如，操作1006和1008可以同时执行。
图11是根据本申请的一些实施例所示用于评估初始翻转角配置的流程1100的一种示例性流程图。在一些实施例中，流程1100中的一个或多个操作可以通过图1所示的MRI***100实现。例如，流程1000可以以指令的形式存储于存储设备130中，并由处理引擎120执行调用和/或执行(例如，图3所示的计算设备300的处理器310、图4所示的移动设备400的中央处理器 440)。图10中描述的操作1004可以根据流程1100执行。
在1102中，可以确定第一参数，所述第一参数与初始翻转角配置相对应的回波链相关。操作1102可以由计算单元608执行。在一些实施例中，可以根据图8、9和/或10中所述来确定所述初始翻转角配置。所述回波链可以基于多个具有所述初始翻转角配置的回聚射频脉冲获得。在一些实施例中，所述具有初始角配置的多个回聚射频脉冲可以通过其他参数定义，如相邻两个回聚射频脉冲之间的回聚间隔、回聚射频脉冲的序列号、多个回聚射频脉冲的总能量等。
在一些实施例中，所述第一参数与所述初始翻转角配置对应的回波链中的回波信号强度有关。在一些实施例中，所述回波信号强度可以根据回波和回聚脉冲的翻转角度之间的关系确定。例如，所述回波和回聚脉冲的翻转角度之间的关系可以通过布洛赫方程、EPG算法等来描述。在一些实施例中，所述回波信号强度可以基于被检测的组织的纵向弛豫时间T1和横向弛豫时间T2来确定。例如，可以根据如下所述的方程(6)确定回波链中的回波信号强度：
S(TE)＝S0*fα*ft(TE,T1,T2)， (6)
其中S(TE)表示回波信号强度，S0是常数，TE表示回波时间，T1表示纵向弛豫时间T2表示横向弛豫时间以及fα表示信号强度因子。信号强度因子fα与回聚射频脉冲的翻转角度有关。例如，回聚射频脉冲的翻转角度越大，信号强度因子fα越大。此外，当回聚射频脉冲的翻转角小于180°时，信号强度的因素fα小于1。
在一些实施例中，所述第一参数可以包括特定信号强度。例如，所述特定信号强度可以是最大信号强度。所述最大信号强度可以基于与所述初始翻转角配置对应的第一信号演化确定。所述第一信号演化可以包括基于所述初始翻转角配置得到的回波链中所有回波的信号强度。在一些实施例中，可以基于，例如方程(6)确定所述第一信号演化。又例如，所述特定信号强度可以是所述回波链的中心回波的信号强度。在一些实施例中，可以根据与所述初始翻转角配置相对应的回波链的初始等效TE(也称中心回波的回波时间)确定所述回波链的中心回波。所述初始等效TE可以由用户通过终端140确定，也可以根据MRI***100的默认设置确定。可以基于初始等效TE确定中心回聚射频脉冲。基于所述初始翻转角配置中所述中心回聚射频脉冲的翻转角确定所述中心回聚射频脉冲对应的信号强度因子fα。然后可以根据方程(6)以及所述中心回聚射频脉冲对应的信号强度因子fα和所述初始等效TE确定所述中心回波的信号强度。又例如，所述特定信号强度可以是在所述初始翻转角配置中最大翻转角对应的回波信号强度。然后根据如上所述的方程(6)确定所述特定信号强度。
在1104中，可以确定与所述回波链有关的第二参数。操作1104可以由计算单元608执行。在一些实施例中，所述第二参数可以包括所述初始翻转角配置下的回聚射频脉冲的总能量。在一些实施例中，所述回聚射频脉冲的总能量可以基于本申请中披露的一个或多个回聚射频脉冲相关参数确定，例如回聚射频脉冲的总数量、回聚射频脉冲的翻转角度、回聚射频脉冲的频率等。例如，具有较大翻转角度的回聚射频脉冲的能量较大。又例如，回聚射频脉冲的总数量越大，所述回聚射频脉冲的总能量越大。
在1106中，可以确定所述第一参数和/或第二参数是否满足第一准则。操作1106可以由评估单元604执行。如果所述第一参数和所述第二参数满足所述第一准则，则流程1100可以执行操作1110。如果所述第一参数和/或第二参数不满足所述第一准则，则过程1100可以执行操作1108。在一些实施例中，所述第一准则可以包括与所述第一参数相对应的第一阈值和/或与所述第二参数相对应的第二阈值。所述第一参数和/或第二参数是否满足所述第一准则可以包括所述第一参数是否满足所述第一阈值和/或所述第二参数是否满足所述第二阈值。在一些实施例中，所述第一阈值和第二阈值可以根据MRI***100的默认设置来设置，或者由用户通过终端140设置。在一些实施例中，所述第一阈值可以基于参考信号演化(例如，期望的信号演化)来确定。在一些实施例中，所述参考信号演化可以是基于具有固定翻转角配置的多个回聚射频脉冲确定。所述具有固定翻转角配置的多个回聚射频脉冲可以具有相同的翻转角，例如，所述初始翻转角配置中的最大翻转角。所述参考信号演化可以包括期望信号强度。在一些实施例中，所述第一阈值可以是所述参考信号演化中的最大信号强度。在一些实施例中，所述第一个阈值可以是基于所述固定翻转角配置产生的回波链的中心回波的信号强度。
在一些实施例中，如果所述第一参数等于或大于所述第一阈值且所述第二参数小于所述第二阈值，则所述初始翻转角配置满足所述第一准则。例如，当所述初始翻转角配置对应的回波链的中心回波信号强度大于或等于所述固定翻转角配置对应的回波链的中心回波信号强度且所述初始翻转角配置下的回聚射频脉冲的总能量小于所述第二阈值，则所述第一准则被满足。在一些实施例中，如果所述第一个参数和所述第一阈值之间的差值小于第四阈值(例如，一个常数)且所述第二参数和所述第二阈值之间的差值小于第五阈值(例如，一个常数)，则所述第一准则被满足。所述第四阈值和第五阈值可以根据MRI***100的默认设置来设置，或者可以由用户通过终端140设置。例如，所述第一准则包括所述初始翻转角配置对应的回波链的中心回波信号强度接近所述固定翻转角配置对应的回波链的中心回波信号强度且所述初始翻转角配置下的回聚射频脉冲的总能量接近所述第二阈值(例如，期望值)。本申请中所使用的第一个值“接近”第二值可以表明所述第一值和第二值之间的偏差小于40％，或30％，或20％，或15％，或10％，或8％，或5％，或3％等。又例如，所述第一准则可以是所述第一信号演化中最大信号强度接近于所述固定翻转角配置对应的回波链的中心回波信号强度且所述初始翻转角配置下的回聚射频脉冲的总能量接近所述第二阈值(例如，期望值)。
在1108中，可以调整所述初始翻转角配置。操作1108可以由调整单元 606执行。在一些实施例中，可以调整至少部分的初始翻转角配置。例如，所述初始翻转角配置中对应于中心回波的翻转角可以增加或减少，从而使得所述初始翻转角配置下的回波链的中心回波信号强度可以接近、等于或大于所述固定翻转角配置下的回波链的中心回波信号强度和/或所述初始翻转角配置下的回聚射频脉冲的总能量可以小于所述第二阈值。
在1110中，可以确定第一翻转角配置。操作1110可以由调整单元608执行。所述第一个翻转角配置满足所述第一准则。在一些实施例中，在1110中确定的调整后的初始翻转角配置可以被指定为第一个翻转角配置。在一些实施例中，当所述初始翻转角配置满足所述第一准则，则所述初始翻转角配置可以被指定为所述第一翻转角配置。
本实施例提供一种调整可变翻转角配置的方法，通过第一参数及第二参数调整可变翻转角配置中的翻转角。所述第一参数包括回波信号强度有关，所述第二参数包括回聚射频脉冲的总能量，本实施例可以在保证回聚射频脉冲总能量符合一定条件的情况下保证回聚射频脉冲的强度符合一定条件，以确保降低回聚射频脉冲的特殊吸收率(SAR)的情况下信噪比保持基本不变。
需要注意的是，以上对于评估初始翻转角配置的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，操作1108 以及1110可以同时执行。又例如，操作1102和1104可以同时执行。
图12是根据本申请的一些实施例所示用于评估第一翻转角配置的流程 1200的一种示例性流程图。在一些实施例中，流程1200中的一个或多个操作可以通过图1所示的MRI***100实现。例如，流程1200可以以指令的形式存储于存储设备130中，并由处理引擎120执行调用和/或执行(例如，图3所示的计算设备300的处理器310、图4所示的移动设备400的中央处理器 440)。图10中描述的操作1008可以根据流程1200执行。
在1202中，可以确定与第一翻转角配置相对应的第二信号演化。操作 1202可以由计算单元608执行。在一些实施例中，所述第一翻转角配置可以根据图11中的流程1100确定。在一些实施例中，所述第二信号演化可以多个回波信号强度，所述回波对应于所述第一翻转角配置下的多个回聚射频脉冲。在一些实施例中，所述第二信号演化可以基于本申请中描述的回波和翻转角之间的关系确定。例如，所述第二信号演化可以根据图11所示的公式(6)来确定。
在1204中，可以基于所述第二信号演化来确定与所述第一翻转角配置有关的等效TE。操作1202可以由计算单元608执行。所述等效TE指的是从激发射频脉冲中点到第一中心回波对应的回聚射频脉冲中点的时间间隔。所述第一中心回波是基于所述第一翻转角配置下的回聚射频脉冲产生的回波链的中心回波。在一些实施例中，所述第二信号演化中具有最大信号强度的回波可以被指定为第一中心回波。所述第一个翻转角配置相关的等效TE可以基于所述第一中心回波对应的回聚射频脉冲确定。例如，可以根据如下所述的公式(7)来确定所述等效TE：
Figure BDA0001391490010000441
其中，i表示所示第一中心回波对应的回聚射频脉冲的序列号，ESP表示回聚间隔，也可以表示两个相邻回聚射频脉冲之间的时间间隔。在一些实施例中，所述等效TE可以是所述第一中心回波对应的回聚射频脉冲之前的所有回聚间隔的总和。
在1206中，可以确定与所述第一个翻转角配置有关的等效TE是否满足第二准则。操作1206可以由评估单元606执行。在一些实施例中，如果与所述第一翻转角配置有关的等效TE满足所述第二准则，则流程1200可以执行操作 1212。如果与所述第一翻转角配置有关的等效TE不满足所述第二准则，则流程1200可以执行操作1208。在一些实施例中，所述第二准则可以包括与所述等效TE有关的第三阈值。在一些实施例中，所述第二准则可以是所述等效TE 等于或大于所述第三阈值。在一些实施例中，所述第二准则可以是所述等效 TE和所述第三阈值之间的差值可以小于第六阈值(例如，常数)。
在一些实施例中，所述第三阈值可以包括期望等效TE。在一些实施例中，用户可以根据临床需求(例如，被检查对象的类型)，通过终端140设置期望等效TE。在一些实施例中，所述期望等效TE可以基于固定翻转角配置的回聚射频脉冲确定。例如，所述固定翻转角配置的回聚射频脉冲的翻转角可以与所述第一中心回波对应的回聚射频脉冲的翻转角相同。然后，可以基于所述固定翻转角配置的回聚射频脉冲确定回波信号强度。具有最大信号强度的回波可以被确定为与所述固定翻转角配置相对应的第二中心回波。可以基于所述第二中心回波确定与所述第二中心回波相对应的回聚射频脉冲的序列号。然后，可以根据，例如方程(7)，确定所述期望等效TE。
在1208中，可以基于所述第二准则和与所述第一翻转角配置有关的等效 TE确定需要调整的回聚射频脉冲的翻转角数量。操作1208可以由计算单元 608执行。在一些实施例中，所述第二准则可以包括期望等效TE，例如，在 1206中确定的期望等效TE。然后可以基于所述期望等效TE与所述第一翻转角配置有关的等效TE确定需要调整的回聚射频脉冲的翻转角数量。例如，所述需要调整的回聚射频脉冲的翻转角数量可以根据如下所述的方程确定(8)：
Figure BDA0001391490010000451
其中j表示需要调整的回聚射频脉冲的翻转角数量，ESP表示回聚间隔，也指的是一个两个相邻回聚射频脉冲之间的时间间隔。
在1210中，可以基于需要调整的回聚射频脉冲的翻转角数量调整所述第一翻转角配置。操作1210可以由调整单元606执行。在一些实施例中，可以调整所述第一翻转角配置，使得与所述第一翻转角配置有关的等效TE可以满足所述第二准则，例如，所述第一翻转角配置有关的等效TE可以接近或等于期望等效TE。在一些实施例中，所述期望等效TE可以对应于期望中心回波 (例如，第二中心回波)。在一些实施例中，可以调整所述第一翻转角配置，使得第一回波中心的序列号可以等于期望中心回波的序列号(例如，第二中心回波)。例如，第一中心回波可以是第i回波，所述第i回波对应第i回聚射频脉冲(RFi)。期望中心回波可以是第(i+j)回波，所述第(i+j)回波对应第 (i+j)回聚射频脉冲(RFi+j)。然后调整第i回聚射频脉冲(RFi)与第(i+ j)回聚射频脉冲(RFi+j)之间的回聚射频脉冲的翻转角度，可以使得RFi+j满足期望中心回波(例如，第二中心回波)。此外，第i回聚射频脉冲(RFi)与第 (i+j)回聚射频脉冲(RFi+j)之间的回聚射频脉冲的翻转角度可以单调递增，从而使得RFi+j的翻转角度最大。
在一些实施例中，RFi和RFi+j可以处于不同的期相。例如，所述第一翻转角配置可分为三个期相，包括第一期相、第二期相和第三期相。所述第一期相中的翻转角度可以单调递减，所述第二期相中的翻转角度可以单调递增，以及所述第三期相中的翻转角度可以单调递减。RFi可以设置于所述第二期相中，并且是所述第二期相中的最后的回聚射频脉冲。RFi+j可以设置在第三期相中。 RFi和RFi+j之间的回聚射频脉冲的翻转角度小于RFi的翻转角并单调递减。则可以通过调整RFi和RFi+j之间的回聚射频脉冲的翻转角度使得RFi和RFi+j之间的回聚射频脉冲的翻转角度单调递增。从而使得RFi+j的翻转角最大。调整后的 RFi和RFi+j之间的回聚射频脉冲位于第二期相。所述第二期相中的回聚射频脉冲数量会增加，所述第三期相中的回聚射频脉冲数量会减少，但是总的回聚射频脉冲数量不变。
在1212，可以确定第二个翻转角配置。操作1212可以由评估单元604执行。所述第二个翻转角配置可以满足所述第二准则。例如，所述第二翻转角配置对应的回波链的等效TE等于或接近所述第三阈值(例如，期望等效TE)，所述回波链包括多个回波信号，所述回波是基于在所述第二翻转角配置下的多个回聚射频脉冲获得的。在一些实施例中，调整后的第一翻转角配置可以被指定为所述第二翻转角配置。在一些实施例中，当所述第一翻转角配置对应的等效TE满足所述第二准则，则所述第一翻转角配置可以被指定为所述第二翻转角配置。
本实施例提供一种调整可变翻转角配置的方法，通过回波信号的等效回波时间调整可变翻转角配置的翻转角，以确保图像对比度保持基本不变。
需要注意的是，以上对于评估第一翻转角配置的描述，仅为描述方便，并不能把本申请限制在所举实施例范围之内。可以理解，对于本领域的技术人员来说，在了解该***的原理后，可能在不背离这一原理的情况下，对实施上述方法和***的应用领域进行形式和细节上的各种修正和改变。例如，操作1210 以及1212可以同时执行。又例如，操作1202可以不需要。
具体实施例
以下具体实施例只是为了方便说明，并不会限制本申请的范围。
具体实施例一：翻转角曲线示意图
图13是根据本申请的一些实施例所示回波链对应的翻转角变化曲线的一种示意图。如图13所示，横坐标表示回波链长度(ETL)，也可以称为回波序列号，纵坐标表示回聚射频脉冲的翻转角度。实线代表固定翻转角度配置，虚线代表可变翻转角度配置。所述固定翻转角配置包括相同的145°翻转角度。在一些实施例中，所述固定翻转角配置中的第一翻转角度(也称为起始翻转角度)可以与其他固定翻转角度不同。例如，所述第一翻转角可以等于90°和固定翻转角度(例如145°)之和的一半。所述固定翻转角配置下的回波链的等效 TE可以为从激励脉冲的中点到第八个回聚射频脉冲的中点(即第八个回波) 之间的时间间隔。所述可变翻转角配置下的回波链的等效TE可以为从激励脉冲的中点到第十个回聚射频脉冲的中点(即第十个回波)之间的时间间隔。可变翻转角配置下的回聚射频脉冲的总能量相对于固定翻转角配置下的回聚射频脉冲的总能量降低36％。
具体实施例二：基于不同回波链重建的MR图像
图14和15是根据本申请的一些实施例所示基于不同翻转角配置情况下获取的回波链重建的MR图像。图14中所示的图像是基于，如图13所示的固定翻转角配置下的回聚射频脉冲获得的MR信号重建的。图15中所示的图像是基于，如图13所示的可变翻转角配置下的回聚射频脉冲获得的MR信号重建的。图14和15所示图像的信噪比和对度相似，但是可变翻转角配置下的回聚射频脉冲的总能量相对于固定翻转角配置下的回聚射频脉冲的总能量降低 36％。
以上概述了血管提取方法的不同方面和/或通过程序实现其他步骤的方法。技术中的程序部分可以被认为是以可执行的代码和/或相关数据的形式而存在的“产品”或“制品”，是通过计算机可读的介质所参与或实现的。有形的、永久的储存介质包括任何计算机、处理器、或类似设备或相关的模块所用到的内存或存储器。例如各种半导体存储器、磁带驱动器、磁盘驱动器或者类似任何时间能够为软件提供存储功能的设备。
所有软件或其中的一部分有时可能会通过网络进行通信，如互联网或其他通信网络。此类通信能够将软件从一个计算机设备或处理器加载到另一个。例如：从按需服务***的一个管理服务器或主机计算机加载至一个计算机环境的硬件平台，或其他实现***的计算机环境，或与提供按需服务所需要的信息相关的类似功能的***。因此，另一种能够传递软件元素的介质也可以被用作局部设备之间的物理连接，例如光波、电波、电磁波等，通过电缆、光缆或者空气实现传播。用来载波的物理介质如电缆、无线连接或光缆等类似设备，也可以被认为是承载软件的介质。在这里的用法除非限制了有形的“储存”介质，其他表示计算机或机器“可读介质”的术语都表示在处理器执行任何指令的过程中参与的介质。
因此，一个计算机可读的介质可能有多种形式，包括但不限于，有形的存储介质，载波介质或物理传输介质。稳定的储存介质包括：光盘或磁盘，以及其他计算机或类似设备中使用的，能够实现图中所描述的***组件的存储***。不稳定的存储介质包括动态内存，例如计算机平台的主内存。有形的传输介质包括同轴电缆、铜电缆以及光纤，包括计算机***内部形成总线的线路。载波传输介质可以传递电信号、电磁信号，声波信号或光波信号，这些信号可以由无线电频率或红外数据通信的方法所产生的。通常的计算机可读介质包括硬盘、软盘、磁带、任何其他磁性介质；CD-ROM、DVD、DVD-ROM、任何其他光学介质；穿孔卡、任何其他包含小孔模式的物理存储介质；RAM、 PROM、EPROM、FLASH-EPROM，任何其他存储器片或磁带；传输数据或指令的载波、电缆或传输载波的连接装置、任何其他可以利用计算机读取的程序代码和/或数据。这些计算机可读介质的形式中，会有很多种出现在处理器在执行指令、传递一个或更多结果的过程之中。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写，包括面向对象编程语言如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 (13)

1.一种用于获取核磁共振信号的方法，所述方法包括：
获取回聚射频脉冲的初始翻转角配置，所述回聚射频脉冲用于产生回波链；
判断所述初始翻转角配置是否满足第一准则以获得第一判断结果，所述第一准则与所述初始翻转角配置对应的回波链的第一参数和第二参数有关，所述第一参数包括所述回波链中至少一个回波的强度值，所述第二参数包括所述初始翻转角配置的回聚射频脉冲的总能量；
基于所述第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；
判断所述第一翻转角配置是否满足第二准则以获得第二判断结果，所述第二准则与所述回波链的第三参数有关，所述第三参数包括与所述第一翻转角配置对应的回波链的等效回波时间；
基于所述第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及
基于所述第二翻转角配置，获取核磁共振信号。
2.权利要求1所述的方法，其特征在于，所述基于所述第一判断结果，确定第一翻转角配置，进一步包括：
确定所述初始翻转角配置不满足所述第一准则；以及
调整所述初始翻转角配置以确定所述第一翻转角配置。
3.权利要求2所述的方法，其特征在于，所述判断所述初始翻转角配置是否满足第一准则以获得第一判断结果，包括：
基于所述初始翻转角配置，确定第一信号演变；以及
通过确定所述第一信号演变不满足所述第一准则确定所述初始翻转角配置不满足所述第一准则，其中，所述第一信号演变包括所述回波链中至少部分的回波的强度值。
4.权利要求3所述的方法，其特征在于，所述第一信号演变与横向弛豫时间或纵向弛豫时间有关。
5.权利要求3所述的方法，其特征在于，所述第一准则包括第一阈值，所述第一阈值与所述第一信号演变有关，以及所述确定所述初始翻转角配置不满足所述第一准则包括确定所述第一信号演变中最大强度值小于所述第一阈值。
6.权利要求1所述的方法，其特征在于，所述判断所述初始翻转角配置是否满足第一准则以获得第一判断结果包括：
基于所述初始翻转角配置，确定所述回聚射频脉冲的总能量；以及
基于所述第一准则评价所述回聚射频脉冲的总能量。
7.权利要求6所述的方法，其特征在于，所述第一准则包括第二阈值，所述第二阈值对应所述回聚射频脉冲的总能量，以及确定所述初始翻转角配置不满足所述第一准则包括确定所述回聚射频脉冲的总能量等于或大于所述第二阈值。
8.权利要求1所述的方法，其特征在于，所述基于所述第二判断结果，确定第二翻转角配置进一步包括：
确定所述第一翻转角配置不满足所述第二准则；以及
调整所述第一翻转角配置以确定所述第二翻转角配置。
9.权利要求1所述的方法，其特征在于，所述第二准则包括第三阈值，所述第三阈值与所述第一翻转角配置对应的回波链的等效回波时间对应。
10.权利要求9所述的方法，其特征在于，所述确定所述第一翻转角配置不满足第二准则包括：
确定第二信号演变，所述第二信号演变包括与所述第一翻转角配置对应的回波链中至少部分的回波的信号强度；
基于所述第二信号演变，确定所述第一翻转角配置对应的回波链的等效回波时间；以及
通过确定所述第一翻转角配置对应的回波链的等效回波时间不满足所述第二准则确定所述第一翻转角配置不满足所述第二准则。
11.权利要求1所述的方法，其特征在于，所述基于所述第二翻转角配置获取核磁共振信号进一步包括：
基于所述第一准则，评价所述第二翻转角配置；以及
基于所述第二翻转角配置的评价结果，调整所述第二翻转角配置以确定目标翻转角配置，所述目标翻转角配置满足所述第一准则以及所述第二准则。
12.一种用于获取核磁共振信号的***，所述***包括：
至少一个处理器；以及可执行指令，所述可执行指令可以由所述至少一个处理器执行，导致所述***实现一种方法，包括：
获取回聚射频脉冲的初始翻转角配置，所述回聚射频脉冲用于产生回波链；
判断所述初始翻转角配置是否满足第一准则以获得第一判断结果，所述第一准则与所述初始翻转角配置对应的回波链的第一参数和第二参数有关，所述第一参数包括所述回波链中至少一个回波的强度值，所述第二参数包括具有所述初始翻转角配置的回聚射频脉冲的总能量；
基于所述第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；
判断所述第一翻转角配置是否满足第二准则以获得第二判断结果，所述第二准则与所述回波链的第三参数有关，所述第三参数包括与所述第一翻转角配置对应的回波链的等效回波时间；
基于所述第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及
基于所述第二翻转角配置，获取核磁共振信号。
13.一种非暂时性的计算机可读介质，包括可执行指令，所述可执行指令被至少一个处理器执行时，导致所述至少一个处理器实现一种方法，包括：
获取回聚射频脉冲的初始翻转角配置，所述回聚射频脉冲用于产生回波链；
判断所述初始翻转角配置是否满足第一准则以获得第一判断结果，所述第一准则与所述初始翻转角配置对应的回波链的第一参数和第二参数有关，所述第一参数包括所述回波链中至少一个回波的强度值，所述第二参数包括具有所述初始翻转角配置的回聚射频脉冲的总能量；
基于所述第一判断结果，确定第一翻转角配置，所述第一翻转角配置满足所述第一准则；
判断所述第一翻转角配置是否满足第二准则以获得第二判断结果，所述第二准则与所述回波链的第三参数有关，所述第三参数包括与所述第一翻转角配置对应的回波链的等效回波时间；
基于所述第二判断结果，确定第二翻转角配置，所述第二翻转角配置满足所述第二准则；以及
基于所述第二翻转角配置，获取核磁共振信号。
CN201710752718.XA 2017-08-28 2017-08-28 磁共振翻转角的计算方法、***和计算机可读存储介质 Active CN107638179B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201710752718.XA CN107638179B (zh) 2017-08-28 2017-08-28 磁共振翻转角的计算方法、***和计算机可读存储介质
US15/692,148 US10739430B2 (en) 2017-08-28 2017-08-31 System and method for determining flip angles in magnetic resonance imaging
EP17204380.4A EP3451009B1 (en) 2017-08-28 2017-11-29 System and method for determining flip angles of refocusing radio frequency pulses in magnetic resonance imaging
PCT/CN2017/114107 WO2019041607A1 (en) 2017-08-28 2017-11-30 SYSTEM AND METHOD FOR DETERMINING ROCK ANGLES IN MAGNETIC RESONANCE IMAGING

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710752718.XA CN107638179B (zh) 2017-08-28 2017-08-28 磁共振翻转角的计算方法、***和计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN107638179A CN107638179A (zh) 2018-01-30
CN107638179B true CN107638179B (zh) 2020-09-04

Family

ID=61110827

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710752718.XA Active CN107638179B (zh) 2017-08-28 2017-08-28 磁共振翻转角的计算方法、***和计算机可读存储介质

Country Status (3)

Country Link
US (1) US10739430B2 (zh)
CN (1) CN107638179B (zh)
WO (1) WO2019041607A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107638179B (zh) * 2017-08-28 2020-09-04 上海联影医疗科技有限公司 磁共振翻转角的计算方法、***和计算机可读存储介质
CN108983133B (zh) * 2018-07-23 2020-06-12 中国科学院深圳先进技术研究院 磁共振血管壁成像方法、装置、设备及存储介质
CN111381204B (zh) * 2018-12-28 2022-07-12 西门子(深圳)磁共振有限公司 基于二维快速自旋回波的磁共振成像方法和装置
US11796618B2 (en) * 2019-07-12 2023-10-24 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for magnetic resonance imaging
CN112731235B (zh) * 2020-12-09 2024-04-19 中国科学院深圳先进技术研究院 一种磁共振化学交换饱和转移成像方法以及相关设备

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6492809B1 (en) * 1998-12-04 2002-12-10 Schlumberger Technology Corporation Preconditioning spins near a nuclear magnetic resonance region
US6272369B1 (en) * 1999-01-22 2001-08-07 Ge Medical Systems Global Technology Company Llc Method for optimizing fat suppression using the chemical shift selective MR imaging technique
US6809518B2 (en) * 1999-06-29 2004-10-26 Gilles Beaudoin System and method for converting adiabatic RF pulses into pseudo adiabatic RF pulses
JP3576069B2 (ja) * 2000-04-27 2004-10-13 株式会社東芝 Ｍｒｉ装置
USRE47178E1 (en) * 2000-12-21 2018-12-25 University Of Virginia Patent Foundation Method and apparatus for spin-echo-train MR imaging using prescribed signal evolutions
JP5002099B2 (ja) * 2001-08-31 2012-08-15 株式会社東芝 磁気共鳴イメージング装置
US6956374B2 (en) * 2003-07-02 2005-10-18 General Electric Company Method and apparatus to reduce RF power in high field MR imaging incorporating multi-phase RF pulse flip angles
US7227356B1 (en) * 2004-09-03 2007-06-05 General Electric Company Method and apparatus of multi-echo MR data acquisition with non-discrete flip angle train
DE102007021719B4 (de) * 2007-05-09 2010-09-30 Siemens Ag Turbospinechobildgebungssequenz mit langen Echozügen und optimiertem T1-Kontrast
US7737690B2 (en) * 2007-05-18 2010-06-15 General Electric Company System and method for amplitude reduction in RF pulse design
US8228060B2 (en) * 2007-06-25 2012-07-24 General Electric Company Method and apparatus for generating a flip angle schedule for a spin echo train pulse sequence
DE102008032155B4 (de) * 2008-07-08 2010-05-12 Siemens Aktiengesellschaft Verfahren zur Erzeugung einer Anregung in einem Untersuchungsobjekt und Magnetresonanzanlage
DE102008046022A1 (de) * 2008-09-05 2010-03-25 Siemens Aktiengesellschaft Phasensensitive Flüssigkeitsunterdrückung
US8212562B2 (en) * 2010-01-28 2012-07-03 Kabushiki Kaisha Toshiba Method and apparatus for designing and/or implementing variable flip angle MRI spin echo train
DE102011005649B4 (de) 2011-03-16 2013-07-04 Siemens Aktiengesellschaft Bestimmung des tatsächlichen Kippwinkels und Justierung der Transmitterspannung bei der MR-Bildgebung eines kontinuierlich verfahrenen Untersuchungsobjekts
CN102680928B (zh) * 2011-03-18 2014-09-03 西门子(深圳)磁共振有限公司 磁共振图像处理方法和磁共振图像处理装置
WO2012169350A1 (ja) * 2011-06-09 2012-12-13 株式会社 日立メディコ 磁気共鳴イメージング装置及びフリップ角決定方法
JP5978431B2 (ja) * 2011-06-09 2016-08-24 株式会社日立製作所 磁気共鳴イメージング装置およびフリップ角決定方法
EP2762910A1 (en) 2013-01-30 2014-08-06 Samsung Electronics Co., Ltd Susceptibility-weighted magnetic resonance imaging
DE102013226246B4 (de) * 2013-12-17 2015-06-25 Siemens Aktiengesellschaft Verfahren zur schichtselektiven Magnetresonanz-Bildgebung und Magnetresonanz-Anlage
KR101649275B1 (ko) * 2014-10-21 2016-08-18 삼성전자주식회사 자기 공명 영상 생성 방법 및 그에 따른 의료 영상 장치
WO2017031638A1 (en) * 2015-08-21 2017-03-02 Shanghai United Imaging Healthcare Co., Ltd. System and method for flip angle determination in magnetic resonance imaging
CN107638179B (zh) * 2017-08-28 2020-09-04 上海联影医疗科技有限公司 磁共振翻转角的计算方法、***和计算机可读存储介质

Also Published As

Publication number Publication date
CN107638179A (zh) 2018-01-30
US20190064295A1 (en) 2019-02-28
US10739430B2 (en) 2020-08-11
WO2019041607A1 (en) 2019-03-07

Similar Documents

Publication Publication Date Title
CN109658471B (zh) 一种医学图像重建方法和***
CN107638179B (zh) 磁共振翻转角的计算方法、***和计算机可读存储介质
US7649354B2 (en) Method and apparatus for acquiring magnetic resonance imaging data
CN109310361B (zh) 用于快速磁共振成像的***和方法
US10324153B2 (en) System and method for flip angle determination in magnetic resonance imaging
US20180231631A1 (en) System and method for magnetic resonance imaging
CN112370040B (zh) 磁共振成像方法、装置、存储介质及电子设备
CN112763957B (zh) 磁共振成像***和方法
CN111012349B (zh) 磁共振成像中信号表征确定的***与方法
CN114076911B (zh) 磁共振成像的方法及***
WO2018137191A1 (en) System and method for magnetic resonance imaging
US11022667B2 (en) System and method for image reconstruction
US11067654B2 (en) Systems and methods for determining field map
CN108896940B (zh) 用于补偿mri的梯度脉冲的***和方法
US10634749B2 (en) System and method for amplitude reduction in RF pulse design
US11940516B2 (en) Systems and methods for magnetic resonance imaging
EP3451009B1 (en) System and method for determining flip angles of refocusing radio frequency pulses in magnetic resonance imaging
JP7230149B2 (ja) 磁気共鳴イメージングで得た画像の処理方法、画像処理プロブラム、及び、計算機
US11408952B2 (en) Systems and methods for pulse sequence configuration in magnetic resonance imaging
CN106943143A (zh) 磁共振成像方法和装置
US12032045B2 (en) Systems and methods for pulse sequence configuration in magnetic resonance imaging
EP3451007B1 (en) System and method for amplitude reduction in rf pulse design
CN116205999A (zh) 磁共振t1映射***和方法

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, north of the city of Jiading District Road No. 2258

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.