CN112630713A - 磁共振成像方法、系统、电子装置和存储介质 - Google Patents

磁共振成像方法、系统、电子装置和存储介质 Download PDF

Info

Publication number
CN112630713A
CN112630713A CN202011630562.6A CN202011630562A CN112630713A CN 112630713 A CN112630713 A CN 112630713A CN 202011630562 A CN202011630562 A CN 202011630562A CN 112630713 A CN112630713 A CN 112630713A
Authority
CN
China
Prior art keywords
magnetic resonance
magnetic field
information
field distribution
distribution information
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
Application number
CN202011630562.6A
Other languages
English (en)
Other versions
CN112630713B (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 CN202011630562.6A priority Critical patent/CN112630713B/zh
Publication of CN112630713A publication Critical patent/CN112630713A/zh
Application granted granted Critical
Publication of CN112630713B publication Critical patent/CN112630713B/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/543Control 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
    • 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
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities

Landscapes

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

本申请涉及一种磁共振成像方法、系统、电子装置和存储介质,其中,该磁共振成像方法包括:通过局部线圈获取受检者扫描部位的磁共振信号;通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息;根据梯度磁场分布信息,得到磁共振信号的空间编码信息;根据空间编码信息,对磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。通过本申请,解决了相关技术中梯度磁场发生不可预料的变化所导致的磁共振图像重建的准确度低的问题,实现了提高磁共振成像准确度的技术效果。

Description

磁共振成像方法、系统、电子装置和存储介质
技术领域
本申请涉及磁共振成像技术领域,特别是涉及一种磁共振成像方法、系统、电子装置和存储介质。
背景技术
磁共振共像(Magnetic Resonance Imaging,简称为MRI)是利用原子核在强磁场内发生共振产生的信号经图像重建的一种成像技术,是一种核物理现象。它是利用射频脉冲对置于磁场中含有自旋不为零的原子核进行激励,射频脉冲停止后,原子核进行弛豫,在其弛豫过程中用感应线圈采集信号,按一定的数学方法重建形成数学图像。
MRI成像技术不同于其他成像技术,它提供的信息量远远大于医学影像学中的其他许多成像技术。因此,对疾病的诊断具有很大的明显优越性。可以直接作出横断面、矢状面、冠状面和各种斜面的体层图像,不会产生电子计算机断层扫描(Computed Tomography,简称为CT)检测中的伪影。
在磁共振中无论是静磁场还是梯度磁场都非100%受控的,在空间和时间上都会由于各种机制产生不可准确预料和估算的变化和变形,这使得完全依赖精准磁场信息进行时空编码的磁共振图像重建变得不准确,而造成各种图像伪影和变形。因此,如何避免出现由于梯度磁场变化所导致的磁共振图像出现伪影的情况成为了亟待解决的问题。
相关技术中的磁共振成像系统往往通过在磁共振成像系统中的各类线圈上施加电流,并根据电流在线圈中的分布预测磁场分布来确定磁共振系统中的磁场信息,然而,这类技术方案测得的磁场信息往往与实际磁场信息偏差较大,这会导致磁共振图像出现图像伪影和变形,降低磁共振图像重建的准确率。
目前针对相关技术中梯度磁场发生不可预料的变化所导致的磁共振图像重建的准确度低的问题,尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种磁共振成像方法、系统、电子装置和存储介质,以至少解决相关技术中梯度磁场发生不可预料的变化所导致的磁共振图像重建的准确度低的问题。
第一方面,本申请实施例提供了一种磁共振成像方法,包括:通过局部线圈获取受检者扫描部位的磁共振信号;通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息;根据所述梯度磁场分布信息,得到磁共振信号的空间编码信息;根据所述空间编码信息,对所述磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。
在其中一些实施例中,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息包括:通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息;对所述实时磁场信息进行球谐函数展开,得到实时磁场分布信息;从所述实时磁场分布信息中提取所述实时梯度磁场分布信息。
在其中一些实施例中,通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息包括:通过预设磁场将设置于所述钻石NV色心传感器中的钻石NV色心的总磁共振峰劈裂得到四对磁共振峰,其中,所述四对磁共振峰分别对应所述钻石NV色心的四种主轴方向;根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对所述磁共振峰对应的主轴方向下的磁场矢量信息;根据四种所述主轴方向下的磁场矢量信息,计算得到所述受检者扫描部位的实时磁场信息。
在其中一些实施例中,根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对所述磁共振峰对应的主轴方向下的磁场矢量信息包括:根据磁共振峰劈裂宽度与磁场大小的对应关系,分别计算得到所述钻石NV色心在四种主轴方向下的色心磁场矢量信息;根据所述色心磁场矢量信息,计算得到所述预设磁场在四种主轴方向下的磁场矢量信息。
在其中一些实施例中,所述磁场分布信息包括静磁场分布信息以及梯度磁场分布信息;从所述实时磁场分布信息中提取所述实时梯度磁场分布信息包括:根据预设静磁场的磁场信息从所述磁场分布信息中剔除所述静磁场分布信息,得到所述实时梯度磁场分布信息。
在其中一些实施例中,所述空间编码信息包括:所述磁共振信号的层面信息以及所述磁共振信号的位置信息;根据所述梯度磁场分布信息,得到所述磁共振信号的空间编码信息包括:在所述空间编码信息包括所述层面信息的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行层面选择,得到所述磁共振信号的层面信息;在所述空间编码信息包括所述位置信息的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行位置编码,得到所述磁共振信号的位置信息。
在其中一些实施例中,所述磁共振信号的位置信息包括:所述磁共振信号的行位置信息以及所述磁共振信号的列位置信息;所述位置编码包括相位编码以及频率编码;根据所述梯度磁场分布信息,对所述磁共振信号进行位置编码,得到所述磁共振信号的位置信息包括:在所述位置编码包括相位编码的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行相位编码,得到所述磁共振信号的行位置信息;在所述位置编码包括频率编码的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行频率编码,得到所述磁共振信号的列位置信息。
第二方面,本申请实施例提供了一种磁共振成像系统,所述磁共振成像系统包括:磁共振扫描仪;以及被配置为在受检者位于磁共振扫描仪中时操作所述磁共振扫描仪,通过从受检者感兴趣区域中采集磁共振信号来执行诊断扫描的处理器;以及设置于受检者感兴趣区域位置的局部线圈;以及存储有计算机程序的存储器;其中,所述局部线圈上设置有预设数量的钻石NV色心传感器,所述处理器还被配置为运行所述计算机程序以执行如上述第一方面所述的磁共振成像方法。
第三方面,本申请实施例提供了一种电子装置,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上述第一方面所述的磁共振成像方法。
第四方面,本申请实施例提供了一种存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述第一方面所述的磁共振成像方法。
相比于相关技术,本申请实施例提供的一种磁共振成像方法、系统、电子装置和存储介质,解决了相关技术中梯度磁场发生不可预料的变化所导致的磁共振图像重建的准确度低的问题,实现了提高磁共振成像准确度的技术效果。
本申请的一个或多个实施例的细节在以下附图和描述中提出,以使本申请的其他特征、目的和优点更加简明易懂。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是根据本申请实施例的磁共振成像系统的结构示意图;
图2是根据本申请实施例的磁共振成像方法的流程图;
图3是根据本申请实施例的第一光探测磁共振图谱示意图;
图4是根据本申请实施例的第二光探测磁共振图谱示意图;
图5是根据本申请实施例的电子装置的硬件结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行描述和说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。基于本申请提供的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。此外,还可以理解的是,虽然这种开发过程中所作出的努力可能是复杂并且冗长的,然而对于与本申请公开的内容相关的本领域的普通技术人员而言,在本申请揭露的技术内容的基础上进行的一些设计,制造或者生产等变更只是常规的技术手段,不应当理解为本申请公开的内容不充分。
在本申请中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域普通技术人员显式地和隐式地理解的是,本申请所描述的实施例在不冲突的情况下,可以与其它实施例相结合。
除非另作定义,本申请所涉及的技术术语或者科学术语应当为本申请所属技术领域内具有一般技能的人士所理解的通常意义。本申请所涉及的“一”、“一个”、“一种”、“该”等类似词语并不表示数量限制,可表示单数或复数。本申请所涉及的术语“包括”、“包含”、“具有”以及它们任何变形,意图在于覆盖不排他的包含;例如包含了一系列步骤或模块(单元)的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元,而是可以还包括没有列出的步骤或单元,或可以还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。本申请所涉及的“连接”、“相连”、“耦接”等类似的词语并非限定于物理的或者机械的连接,而是可以包括电气的连接,不管是直接的还是间接的。本申请所涉及的“多个”是指大于或者等于两个。“和/或”描述关联对象的关联关系,表示可以存在三种关系,例如,“A和/或B”可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。本申请所涉及的术语“第一”、“第二”、“第三”等仅仅是区别类似的对象,不代表针对对象的特定排序。
本申请所涉及的系统和方法不仅可用于非侵入成像,其所涉及的处理系统可以包括磁共振成像系统(MR系统)以及正电子发射计算机断层显像-磁共振多模态混合系统(PET-MR系统)等。本申请所涉及的方法、装置、系统或者存储介质既可以与上述的处理系统集成在一起,也可以是相对独立的。
下面将以磁共振成像系统为例对本申请实施例进行说明。
本申请实施例提供了一种磁共振成像系统。图1是根据本申请实施例的磁共振成像系统的结构示意图,如图1所示,该磁共振成像系统包括:扫描仪和计算机,其中计算机包括存储器125、处理器122及存储在存储器125上并可在处理器122上运行的计算机程序。其中,处理器122被配置为运行计算机程序以执行本申请实施例的磁共振成像方法。
扫描仪具有孔腔,其通常包括磁共振机架,机架内有主磁体101,主磁体101可以是由超导线圈构成,用来产生静磁场,在一些情况下也可以采用永磁体。主磁体101可以用来产生0.2特斯拉、0.5特斯拉、1.0特斯拉、1.5特斯拉、3.0特斯拉或者更高强度的静磁场。在磁共振成像时,成像对象150会由患者床106进行承载,随着床板的移动,将成像对象150移入静磁场磁场分布较为均匀的区域105内。通常对于磁共振成像系统,如图1所示,空间坐标系(即磁共振成像系统的坐标系)的z方向设置为与磁共振成像系统机架的轴向相同,通常将患者的身长方向与z方向保持一致进行成像,磁共振成像系统的水平平面设置为xz平面,x方向与z方向垂直,y方向与x和z方向均垂直。
在磁共振成像中,脉冲控制单元111控制射频脉冲产生单元116产生射频脉冲,射频脉冲由放大器放大后,经过开关控制单元117,最终由体线圈103或者局部线圈104发出,对成像对象150进行射频激发。成像对象150根据射频激发,会由共振产生相应的射频信号。在接收成像对象150根据激发产生的射频信号时,可以由体线圈103或者局部线圈104进行接收,射频接收链路可以有很多条,射频信号发送到射频接收单元118后,进一步发送到图像重建单元121进行图像重建,形成磁共振图像。
局部线圈104上设置有预设数量的钻石NV(氮空位,Nitrogen Vacancy,简称为NV)色心传感器107,其中,钻石NV色心传感器107可以设置在局部线圈104的内表面或者外表面,局部线圈104上可以设置四个或者四个以上的钻石NV色心传感器107,通过多个钻石NV色心传感器107共同采集成像对象150感兴趣区域的磁场信息,有利于剔除噪声,提高磁共振图像重建的准确度,钻石NV色心传感器107的数据交互可以通过光纤传输进行,实现对钻石NV色心传感器107所在位置的磁场的实时监控。
磁共振扫描仪还包括梯度线圈102,梯度线圈可以用来在磁共振成像时对射频信号进行空间编码。脉冲控制单元111控制梯度信号产生单元112产生梯度信号,梯度信号通常会分为三个相互正交方向的信号:x方向、y方向和z方向,不同方向的梯度信号经过梯度放大器(113、114、115)放大后,由梯度线圈102发出,在区域105内产生梯度磁场。
脉冲控制单元111、图像重建单元121与处理器122、显示单元123、输入/输出设备124、存储器125、通信端口126之间可以通过通信总线127进行数据传输,从而实现对磁共振成像过程的控制。
其中,处理器122可以由一个或多个处理器组成,可以包括中央处理器(CPU),或者特定集成电路(Application Specific Integrated Circuit,简称为ASIC),或者可以被配置成实施本申请实施例的一个或多个集成电路。
其中,显示单元123可以是提供给用户用来显示图像的显示器。
其中,输入/输出设备124可以是键盘、鼠标、控制盒等相关设备,支持输入/输出相应数据流。
其中,存储器125可以包括用于数据或指令的大容量存储器。举例来说而非限制,存储器125可包括硬盘驱动器(Hard Disk Drive,简称为HDD)、软盘驱动器、闪存、光盘、磁光盘、磁带或通用串行总线(Universal Serial Bus,简称为USB)驱动器或者两个或更多个以上这些的组合。在合适的情况下,存储器125可包括可移除或不可移除(或固定)的介质。在合适的情况下,存储器125可在数据处理装置的内部或外部。在特定实施例中,存储器125是非易失性固态存储器。在特定实施例中,存储器125包括只读存储器(ROM)。在合适的情况下,该ROM可以是掩模编程的ROM、可编程ROM(PROM)、可擦除PROM(EPROM)、电可擦除PROM(EEPROM)、电可改写ROM(EAROM)或闪存或者两个或更多个以上这些的组合。存储器125可以用来存储需要处理和/或通信使用的各种数据文件,以及处理器122所执行的可能的程序指令。当处理器122执行存储器125中的存储的指定程序时,该处理器122可执行被本申请提出的磁共振成像方法。
其中,通信端口126可以实现与其他部件例如:外接设备、图像采集设备、数据库、外部存储以及图像处理工作站等之间进行数据通信。
其中,通信总线127包括硬件、软件或两者,将磁共振成像系统的部件彼此耦接在一起。举例来说而非限制,总线可包括加速图形端口(AGP)或其他图形总线、增强工业标准架构(EISA)总线、前端总线(FSB)、超传输(HT)互连、工业标准架构(ISA)总线、无限带宽互连、低引脚数(LPC)总线、存储器总线、微信道架构(MCA)总线、外围组件互连(PCI)总线、PCI-Express(PCI-X)总线、串行高级技术附件(SATA)总线、视频电子标准协会局部(VLB)总线或其他合适的总线或者两个或更多个以上这些的组合。在合适的情况下,通信总线127可包括一个或多个总线。尽管本申请实施例描述和示出了特定的总线,但本申请考虑任何合适的总线或互连。
在其中一些实施例中,处理器122被配置为通过局部线圈获取受检者扫描部位的磁共振信号;通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息;根据梯度磁场分布信息,得到磁共振信号的空间编码信息;根据空间编码信息,对磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。
在其中一些实施例中,处理器122被配置为通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息;对实时磁场信息进行球谐函数展开,得到实时磁场分布信息;从实时磁场分布信息中提取实时梯度磁场分布信息。
在其中一些实施例中,处理器122被配置为通过预设磁场将设置于钻石NV色心传感器中的钻石NV色心的总磁共振峰劈裂得到四对磁共振峰,其中,四对磁共振峰分别对应钻石NV色心的四种主轴方向;根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对磁共振峰对应的主轴方向下的磁场矢量信息;根据四种主轴方向下的磁场矢量信息,计算得到受检者扫描部位的实时磁场信息。
在其中一些实施例中,处理器122被配置为根据磁共振峰劈裂宽度与磁场大小的对应关系,分别计算得到钻石NV色心在四种主轴方向下的色心磁场矢量信息;根据色心磁场矢量信息,计算得到预设磁场在四种主轴方向下的磁场矢量信息。
在其中一些实施例中,处理器122被配置为根据预设静磁场的磁场信息从磁场分布信息中剔除静磁场分布信息,得到实时梯度磁场分布信息。
在其中一些实施例中,空间编码信息包括:磁共振信号的层面信息以及磁共振信号的位置信息;处理器122被配置为在空间编码信息包括层面信息的情况下,根据梯度磁场分布信息,对磁共振信号进行层面选择,得到磁共振信号的层面信息;在空间编码信息包括位置信息的情况下,根据梯度磁场分布信息,对磁共振信号进行位置编码,得到磁共振信号的位置信息。
在其中一些实施例中,磁共振信号的位置信息包括:磁共振信号的行位置信息以及磁共振信号的列位置信息;位置编码包括相位编码以及频率编码;处理器122被配置为在位置编码包括相位编码的情况下,根据梯度磁场分布信息,对磁共振信号进行相位编码,得到磁共振信号的行位置信息;在位置编码包括频率编码的情况下,根据梯度磁场分布信息,对磁共振信号进行频率编码,得到磁共振信号的列位置信息。
本实施例提供了一种磁共振成像方法,图2是根据本申请实施例的磁共振成像方法的流程图,如图2所示,该流程包括如下步骤:
步骤S201,通过局部线圈获取受检者扫描部位的磁共振信号。
步骤S202,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息。
步骤S203,根据梯度磁场分布信息,得到磁共振信号的空间编码信息。
步骤S204,根据空间编码信息,对磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。
在本实施例中,局部线圈用于作为射频接收线圈采集受检者扫描部位的磁共振信号,局部线圈包括但不限于以下至少之一:头部线圈、脊柱线圈、腹部线圈、下肢线圈,其中,局部线圈上还可以设置有预设数量的钻石NV色心传感器,钻石NV色心传感器可以设置在局部线圈的内表面或者外表面,钻石NV色心传感器的数据交互可以通过光纤传输进行,实现对钻石NV色心传感器所在位置的磁场的实时监控;局部线圈上可以设置四个或者四个以上的钻石NV色心传感器,通过多个钻石NV色心传感器共同采集受检者扫描部位的梯度磁场分布信息,有利于剔除噪声,提高磁共振图像重建的准确度。
目前,基于钻石NV色心的磁测量技术相对于传统的磁测量技术,例如霍尔效应传感器、磁力显微镜等,具有诸多优势:工作温度范围宽,空间分辨率和灵敏度高,同时对样品本身磁场没有任何侵扰。
另一方面,作为一个传感器,钻石NV色心由于其原子尺度的大小,能够靠近待测样品到纳米量级,结合其高灵敏度的优势就能够实现单电子自旋、单核自旋的测量。
因此,借鉴传统的磁共振技术,钻石NV色心能够实现单分子核磁共振和单分子顺磁共振,基于钻石NV色心的磁测量技术往往被应用于磁学和磁学材料研究中。
然而,钻石NV色心传感器还并未应用于磁共振成像系统中,相关技术中的磁共振成像系统往往通过在磁共振成像系统中的各类线圈上施加电流,并根据电流在线圈中的分布预测磁场分布来确定磁共振系统中的磁场信息,这类技术方案测得的磁场信息往往与实际磁场信息偏差较大,这会导致磁共振图像出现图像伪影和变形,降低磁共振图像重建的准确率,相比传统的磁共振成像技术,基于钻石NV色心传感器的磁共振成像方法不仅具有纳米量级的空间分辨率和微特斯拉量级的磁灵敏度,同时对受检者没有空间或环境上的要求,测量的频率带宽也可以满足要求,可以在磁共振环境下实时测得磁场信息。
在其中一些实施例中,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息包括如下步骤:
步骤1,通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息。
步骤2,对实时磁场信息进行球谐函数展开,得到实时磁场分布信息。
步骤3,从实时磁场分布信息中提取实时梯度磁场分布信息。
在本实施例中,磁场分布信息可以包括静磁场分布信息以及梯度磁场分布信息,其中,静磁场可以由磁共振系统的机架中设置的主磁体产生,主磁体可以为超导线圈或者永磁体,梯度磁场可以由梯度线圈产生,可以根据预设静磁场的磁场信息,从磁场分布信息中剔除静磁场分布信息,进而得到梯度磁场分布信息。
在其中一些实施例中,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息包括如下步骤:
步骤1,通过预设磁场将设置于钻石NV色心传感器中的钻石NV色心的总磁共振峰劈裂得到四对磁共振峰,其中,四对磁共振峰分别对应钻石NV色心的四种主轴方向。
步骤2,根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对磁共振峰对应的主轴方向下的磁场矢量信息。
步骤3,根据四种主轴方向下的磁场矢量信息,计算得到受检者扫描部位的实时磁场信息。
在本实施例中,四种主轴方向可以分别对应于钻石晶格结构的四个晶轴方向,即(111),(1-11),(-111)和(11-1)四个晶轴方向。
图3是根据本申请实施例的第一光探测磁共振图谱示意图,图中Y轴为荧光强度(Fluorescence intensity),X轴为探头微波频率(Probe microwave frequency),单位为(MHz),如图3所示,可以通过预设磁场将设置于钻石NV色心传感器中的钻石NV色心的总磁共振峰劈裂得到a、b、c、d四对磁共振峰,从八个磁共振峰的中心分开每一对左右两个磁共振峰,图中中心点为2.87GHz,每一对磁共振峰对应钻石NV色心轴在碳晶格中的四种主轴方向。
图4是根据本申请实施例的第二光探测磁共振图谱示意图,图中Y轴为荧光强度(Fluorescence intensity),X轴为探头微波频率(Probe microwave frequency),单位为(MHz),如图4所示,可以根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对磁共振峰对应的主轴方向下的磁场矢量信息,例如:以对应关系f±=2γBz为例,其中,f±=f+-f-,即该对磁共振峰的劈裂宽度,f+为该对磁共振峰中右边磁共振峰的探头微波频率,f-为该对磁共振峰中左边磁共振峰的探头微波频率,γ为钻石NV色心的旋磁比,Bz即预设磁场在Z轴方向下的磁场矢量信息,可以通过该对应关系,得到预设磁场在Z轴方向下的磁场矢量信息,图中△v为中心线宽,C为实测得到的磁共振峰的荧光强度。
在其中一些实施例中,根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对磁共振峰对应的主轴方向下的磁场矢量信息包括:根据磁共振峰劈裂宽度与磁场大小的对应关系,分别计算得到钻石NV色心在四种主轴方向下的色心磁场矢量信息;根据色心磁场矢量信息,计算得到预设磁场在四种主轴方向下的磁场矢量信息。
在本实施例中,得到钻石NV色心在碳晶格中的四种主轴方向下的色心磁场矢量信息后,对钻石NV色心在碳晶格中的四种主轴方向下的色心磁场矢量信息进行计算,可以综合得到预设磁场在四种主轴方向下的磁场矢量信息,进一步得到受检者扫描部位的磁场信息。
在其中一些实施例中,空间编码信息包括:磁共振信号的层面信息以及磁共振信号的位置信息;根据梯度磁场分布信息,得到磁共振信号的空间编码信息包括如下步骤:
步骤1,在空间编码信息包括层面信息的情况下,根据梯度磁场分布信息,对磁共振信号进行层面选择,得到磁共振信号的层面信息。
步骤2,在空间编码信息包括位置信息的情况下,根据梯度磁场分布信息,对磁共振信号进行位置编码,得到磁共振信号的位置信息。
在本实施例中,可以通过层面选择梯度线圈发射层面选择梯度脉冲产生层面选择梯度场以选择扫描层面(沿Z轴的横断层面,其中Z轴方向与静磁场方向一致)。
磁共振信号的位置信息包括:磁共振信号的行位置信息以及磁共振信号的列位置信息;位置编码包括相位编码以及频率编码;根据梯度磁场分布信息,对磁共振信号进行位置编码,得到磁共振信号的位置信息包括:在位置编码包括相位编码的情况下,根据梯度磁场分布信息,对磁共振信号进行相位编码,得到磁共振信号的行位置信息;在位置编码包括频率编码的情况下,根据梯度磁场分布信息,对磁共振信号进行频率编码,得到磁共振信号的列位置信息。
在本实施例中,可以通过相位编码梯度线圈发射相位编码梯度脉冲产生相位编码梯度场,以及通过频率编码梯度线圈发射频率编码梯度脉冲产生频率编码梯度场,获取在已选择的扫描层面上的位置信息(扫描层面上分别沿Y轴和X轴的位置信息),相位编码梯度场用于获取扫描层面上各像素所在行的位置,频率编码梯度场用于获取扫描层面上各像素所在列的位置。
通过上述步骤S201至S204,通过局部线圈获取受检者扫描部位的磁共振信号,并通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息,根据梯度磁场分布信息,得到磁共振信号的空间编码信息,最后根据空间编码信息,对磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像,钻石NV色心传感器的尺寸非常小,可以达到1mm3,同时精度以及动态范围较高,测量的频率带宽也可以满足要求,可以在磁共振环境下实时测得磁场信息,避免出现相关技术中测得的磁场信息往往与实际磁场信息偏差较大所导致的磁共振图像出现图像伪影和变形的情况。通过本申请,解决了相关技术中梯度磁场发生不可预料的变化所导致的磁共振图像重建的准确度低的问题,实现了提高磁共振成像准确度的技术效果。
本实施例还提供了一种电子装置,图5是根据本申请实施例的电子装置的硬件结构示意图,如图5所示,该装置包括存储器504和处理器502,该存储器504中存储有计算机程序,该处理器502被设置为运行计算机程序以执行上述任一项方法实施例中的步骤。
具体地,上述处理器502可以包括中央处理器(CPU),或者特定集成电路(Application Specific Integrated Circuit,简称为ASIC),或者可以被配置成实施本申请实施例的一个或多个集成电路。
其中,存储器504可以包括用于数据或指令的大容量存储器504。举例来说而非限制,存储器504可包括硬盘驱动器(Hard Disk Drive,简称为HDD)、软盘驱动器、固态驱动器(Solid State Drive,简称为SSD)、闪存、光盘、磁光盘、磁带或通用串行总线(UniversalSerial Bus,简称为USB)驱动器或者两个或更多个以上这些的组合。在合适的情况下,存储器504可包括可移除或不可移除(或固定)的介质。在合适的情况下,存储器504可在数据处理装置的内部或外部。在特定实施例中,存储器504是非易失性(Non-Volatile)存储器。在特定实施例中,存储器504包括只读存储器(Read-Only Memory,简称为ROM)和随机存取存储器(Random Access Memory,简称为RAM)。在合适的情况下,该ROM可以是掩模编程的ROM、可编程ROM(Programmable Read-Only Memory,简称为PROM)、可擦除PROM(ErasableProgrammable Read-Only Memory,简称为EPROM)、电可擦除PROM(Electrically ErasableProgrammable Read-Only Memory,简称为EEPROM)、电可改写ROM(ElectricallyAlterable Read-Only Memory,简称为EAROM)或闪存(FLASH)或者两个或更多个以上这些的组合。在合适的情况下,该RAM可以是静态随机存取存储器(Static Random-AccessMemory,简称为SRAM)或动态随机存取存储器(Dynamic Random Access Memory,简称为DRAM),其中,DRAM可以是快速页模式动态随机存取存储器504(Fast Page Mode DynamicRandom Access Memory,简称为FPMDRAM)、扩展数据输出动态随机存取存储器(ExtendedDate Out Dynamic Random Access Memory,简称为EDODRAM)、同步动态随机存取内存(Synchronous Dynamic Random-Access Memory,简称SDRAM)等。
存储器504可以用来存储或者缓存需要处理和/或通信使用的各种数据文件,以及处理器502所执行的可能的计算机程序指令。
处理器502通过读取并执行存储器504中存储的计算机程序指令,以实现上述实施例中的任意一种磁共振成像方法。
可选地,上述电子装置还可以包括传输设备506以及输入输出设备508,其中,该传输设备506和上述处理器502连接,该输入输出设备508和上述处理器502连接。
可选地,在本实施例中,上述处理器502可以被设置为通过计算机程序执行以下步骤:
S1,通过局部线圈获取受检者扫描部位的磁共振信号。
S2,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息。
S3,根据梯度磁场分布信息,得到磁共振信号的空间编码信息。
S4,根据空间编码信息,对磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。
需要说明的是,本实施例中的具体示例可以参考上述实施例及可选实施方式中所描述的示例,本实施例在此不再赘述。
另外,结合上述实施例中的磁共振成像方法,本申请实施例可提供一种存储介质来实现。该存储介质上存储有计算机程序;该计算机程序被处理器执行时实现上述实施例中的任意一种磁共振成像方法。
本领域的技术人员应该明白,以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种磁共振成像方法,其特征在于,包括:
通过局部线圈获取受检者扫描部位的磁共振信号;
通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息;
根据所述梯度磁场分布信息,得到所述磁共振信号的空间编码信息;
根据所述空间编码信息,对所述磁共振信号进行傅里叶变换,得到受检者扫描部位的磁共振图像。
2.根据权利要求1所述的磁共振成像方法,其特征在于,通过钻石NV色心传感器获取受检者扫描部位的实时梯度磁场分布信息包括:
通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息;
对所述实时磁场信息进行球谐函数展开,得到实时磁场分布信息;
从所述实时磁场分布信息中提取所述实时梯度磁场分布信息。
3.根据权利要求2所述的磁共振成像方法,其特征在于,通过钻石NV色心传感器获取受检者扫描部位的实时磁场信息包括:
通过预设磁场将设置于所述钻石NV色心传感器中的钻石NV色心的总磁共振峰劈裂得到四对磁共振峰,其中,所述四对磁共振峰分别对应所述钻石NV色心的四种主轴方向;
根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对所述磁共振峰对应的主轴方向下的磁场矢量信息;
根据四种所述主轴方向下的磁场矢量信息,计算得到所述受检者扫描部位的实时磁场信息。
4.根据权利要求3所述的磁共振成像方法,其特征在于,根据磁共振峰劈裂宽度与磁场大小的对应关系,分别确定在与每对所述磁共振峰对应的主轴方向下的磁场矢量信息包括:
根据磁共振峰劈裂宽度与磁场大小的对应关系,分别计算得到所述钻石NV色心在四种主轴方向下的色心磁场矢量信息;
根据所述色心磁场矢量信息,计算得到所述预设磁场在四种主轴方向下的磁场矢量信息。
5.根据权利要求2所述的磁共振成像方法,其特征在于,所述磁场分布信息包括静磁场分布信息以及梯度磁场分布信息;从所述实时磁场分布信息中提取所述实时梯度磁场分布信息包括:
根据预设静磁场的磁场信息从所述磁场分布信息中剔除所述静磁场分布信息,得到所述实时梯度磁场分布信息。
6.根据权利要求1所述的磁共振成像方法,其特征在于,所述空间编码信息包括:所述磁共振信号的层面信息以及所述磁共振信号的位置信息;根据所述梯度磁场分布信息,得到所述磁共振信号的空间编码信息包括:
在所述空间编码信息包括所述层面信息的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行层面选择,得到所述磁共振信号的层面信息;
在所述空间编码信息包括所述位置信息的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行位置编码,得到所述磁共振信号的位置信息。
7.根据权利要求6所述的磁共振成像方法,其特征在于,所述磁共振信号的位置信息包括:所述磁共振信号的行位置信息以及所述磁共振信号的列位置信息;所述位置编码包括相位编码以及频率编码;根据所述梯度磁场分布信息,对所述磁共振信号进行位置编码,得到所述磁共振信号的位置信息包括:
在所述位置编码包括相位编码的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行相位编码,得到所述磁共振信号的行位置信息;
在所述位置编码包括频率编码的情况下,根据所述梯度磁场分布信息,对所述磁共振信号进行频率编码,得到所述磁共振信号的列位置信息。
8.一种磁共振成像系统,其特征在于,所述磁共振成像系统包括:磁共振扫描仪;以及被配置为在受检者位于磁共振扫描仪中时操作所述磁共振扫描仪,通过从受检者感兴趣区域中采集磁共振信号来执行诊断扫描的处理器;以及设置于受检者感兴趣区域位置的局部线圈;以及存储有计算机程序的存储器;其中,所述局部线圈上设置有预设数量的钻石NV色心传感器,所述处理器还被配置为运行所述计算机程序以执行权利要求1至7中任一项所述的磁共振成像方法。
9.一种电子装置,包括存储器和处理器,其特征在于,所述存储器中存储有计算机程序,所述处理器被设置为运行所述计算机程序以执行权利要求1至7中任一项所述的磁共振成像方法。
10.一种存储介质,其特征在于,所述存储介质中存储有计算机程序,其中,所述计算机程序被设置为运行时执行权利要求1至7中任一项所述的磁共振成像方法。
CN202011630562.6A 2020-12-30 2020-12-30 磁共振成像方法、系统、电子装置和存储介质 Active CN112630713B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011630562.6A CN112630713B (zh) 2020-12-30 2020-12-30 磁共振成像方法、系统、电子装置和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011630562.6A CN112630713B (zh) 2020-12-30 2020-12-30 磁共振成像方法、系统、电子装置和存储介质

Publications (2)

Publication Number Publication Date
CN112630713A true CN112630713A (zh) 2021-04-09
CN112630713B CN112630713B (zh) 2022-11-15

Family

ID=75290320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011630562.6A Active CN112630713B (zh) 2020-12-30 2020-12-30 磁共振成像方法、系统、电子装置和存储介质

Country Status (1)

Country Link
CN (1) CN112630713B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113466279A (zh) * 2021-06-30 2021-10-01 中国科学技术大学 宽场超分辨自旋磁成像装置及方法
CN114236439A (zh) * 2021-11-09 2022-03-25 深圳市联影高端医疗装备创新研究院 线圈定位方法、装置、磁共振设备及存储介质
WO2023028786A1 (zh) * 2021-08-30 2023-03-09 深圳先进技术研究院 一种设计非线性梯度线圈的方法以及非线性空间编码方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW201525502A (zh) * 2013-12-31 2015-07-01 Univ Nat Taiwan 一種產生核磁共振影像傾倒角空間分佈的方法及裝置
WO2016161278A1 (en) * 2015-04-01 2016-10-06 The General Hospital Corporation System and method for imaging nanodiamonds as dynamic nuclear polarization agent
CN107356820A (zh) * 2017-06-07 2017-11-17 南京邮电大学 一种基于脉冲光探测磁共振的电磁场近场成像系统及方法
CN107374719A (zh) * 2017-07-24 2017-11-24 南京邮电大学 一种热消融针及应用方法
CN107624162A (zh) * 2015-05-12 2018-01-23 皇家飞利浦有限公司 具有场探头的磁共振检查系统
CN109917315A (zh) * 2019-04-30 2019-06-21 上海联影医疗科技有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质
CN111166321A (zh) * 2019-12-31 2020-05-19 上海联影医疗科技有限公司 心电导联线和心电信号检测装置
CN111650543A (zh) * 2020-06-16 2020-09-11 宿迁学院 一种基于金刚石nv色心的微波近场矢量测量方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW201525502A (zh) * 2013-12-31 2015-07-01 Univ Nat Taiwan 一種產生核磁共振影像傾倒角空間分佈的方法及裝置
WO2016161278A1 (en) * 2015-04-01 2016-10-06 The General Hospital Corporation System and method for imaging nanodiamonds as dynamic nuclear polarization agent
US20180113185A1 (en) * 2015-04-01 2018-04-26 The General Hospital Corporation System and method for imaging nanodiamonds as dynamic nuclear polarization agent
CN107624162A (zh) * 2015-05-12 2018-01-23 皇家飞利浦有限公司 具有场探头的磁共振检查系统
CN107356820A (zh) * 2017-06-07 2017-11-17 南京邮电大学 一种基于脉冲光探测磁共振的电磁场近场成像系统及方法
CN107374719A (zh) * 2017-07-24 2017-11-24 南京邮电大学 一种热消融针及应用方法
CN109917315A (zh) * 2019-04-30 2019-06-21 上海联影医疗科技有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质
CN111166321A (zh) * 2019-12-31 2020-05-19 上海联影医疗科技有限公司 心电导联线和心电信号检测装置
CN111650543A (zh) * 2020-06-16 2020-09-11 宿迁学院 一种基于金刚石nv色心的微波近场矢量测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王红理,等: "《综合与近代物理实验》", 31 August 2015 *
邹德江: "《最新医院大型医疗设备操作规范与诊断技术实务全书》", 31 May 2006 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113466279A (zh) * 2021-06-30 2021-10-01 中国科学技术大学 宽场超分辨自旋磁成像装置及方法
CN113466279B (zh) * 2021-06-30 2024-03-29 中国科学技术大学 宽场超分辨自旋磁成像装置及方法
WO2023028786A1 (zh) * 2021-08-30 2023-03-09 深圳先进技术研究院 一种设计非线性梯度线圈的方法以及非线性空间编码方法
CN114236439A (zh) * 2021-11-09 2022-03-25 深圳市联影高端医疗装备创新研究院 线圈定位方法、装置、磁共振设备及存储介质

Also Published As

Publication number Publication date
CN112630713B (zh) 2022-11-15

Similar Documents

Publication Publication Date Title
CN112630713B (zh) 磁共振成像方法、系统、电子装置和存储介质
US9488709B2 (en) Restriction of the imaging region for MRI in an inhomogeneous magnetic field
EP2689261B1 (en) Compressed sensing mr image reconstruction using constraint from prior acquisition
EP2730935B1 (en) Motional error correction in functional magnetic resonance imaging
RU2616984C2 (ru) Магнитно-резонансная (mr) томография электрических свойств
US10234523B2 (en) MRI with dixon-type water/fat separation with estimation of the main magnetic field variations
JP5848713B2 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
EP2610632A1 (en) MRI with Dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
US10203387B2 (en) MR imaging with enhanced susceptibility contrast
WO2014053927A1 (en) Metal resistant mr imaging reference scan
US20140368195A1 (en) Mr imaging with suppresion of flow artifacts
CN106574954A (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
CN112798995B (zh) 应用于磁共振成像的运动监测方法以及磁共振成像系统
US9274197B2 (en) Magnetic resonance imaging data sampling methods and systems
US9720063B2 (en) Diffusion spectrum imaging systems and methods
US10126395B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US20140121492A1 (en) Mri with separation of different chemical species using a spectral model
US20220308148A1 (en) Mr electric properties tomography without contrast agent
CN112617797B (zh) 应用于磁共振成像的生理信号检测方法以及电子装置
JP4447104B2 (ja) 磁気共鳴装置
US9618595B2 (en) Magnetic resonance imaging system, data processing apparatus, and method for generating magnetic resonance image
US9255979B2 (en) Measuring diffusional anisotropy of ODF lobes having maxima peaks and minima troughs with diffusion weighted MRI
US8285066B2 (en) System and method for generating high resolution images
CN112790863B (zh) 手术器械定位方法、磁共振成像系统以及电子装置
CN112244812B (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