CN113171076B - 基于外插磁化转移信号的cest数据拟合方法、装置及介质 - Google Patents

基于外插磁化转移信号的cest数据拟合方法、装置及介质 Download PDF

Info

Publication number
CN113171076B
CN113171076B CN202110475243.0A CN202110475243A CN113171076B CN 113171076 B CN113171076 B CN 113171076B CN 202110475243 A CN202110475243 A CN 202110475243A CN 113171076 B CN113171076 B CN 113171076B
Authority
CN
China
Prior art keywords
fitting
cest
round
value
magnetization
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
CN202110475243.0A
Other languages
English (en)
Other versions
CN113171076A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110475243.0A priority Critical patent/CN113171076B/zh
Publication of CN113171076A publication Critical patent/CN113171076A/zh
Priority to JP2023566736A priority patent/JP2024522017A/ja
Priority to PCT/CN2021/140311 priority patent/WO2022227639A1/zh
Application granted granted Critical
Publication of CN113171076B publication Critical patent/CN113171076B/zh
Priority to US18/495,715 priority patent/US20240053423A1/en
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/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5605—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by transferring coherence or polarization from a spin species to another, e.g. creating magnetization transfer contrast [MTC], polarization transfer using nuclear Overhauser enhancement [NOE]
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Algebra (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Public Health (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开了一种基于外插磁化转移信号的CEST数据拟合方法、装置及介质，属于磁共振成像领域。本发明提供的基于外插磁化转移信号的CEST数据拟合方法，将MT效应的线型融入到Bloch‑McConnell方程中，通过对该方程代表的模型的拟合，可以得出MT效应，将拟合得出的曲线与实验采集到的曲线作差，就可以去除干扰效应。本发明的基于Bloch‑McConnell的拟合方法，可以考虑实验中实际使用的射频脉冲波形，并且在系统非稳态下也可以使用。

Description

基于外插磁化转移信号的CEST数据拟合方法、装置及介质
技术领域
本发明属于磁共振成像技术领域，具体涉及一种基于外插磁化转移信号的CEST数据拟合方法、装置及介质。
背景技术
化学交换饱和转移(Chemical Exchange Saturation Transfer,CEST)成像是一种磁共振分子成像技术，可以检测活体里游离的蛋白质和多肽。CEST成像已经被应用于多个临床领域，比如脑肿瘤分级、监测放化疗的早期治疗反应、预测肿瘤标志物、检测缺血性中风等。但是CEST信号容易受多种效应的影响，国内外研究者提出了多种计算CEST信号的方法以去除这些干扰效应，这些方法可以大致分成两类，一是使用参考信号与标记信号作差得到最终的信号，二是模型拟合得到最终的信号。非对称磁化转移(Magnetizationtransfer ration asymmetry,MTRasym)假设(1)参考信号中不包含CEST效应；(2)标记信号包含CEST效应；(3)其他效应在参考信号和标记信号中相等。然后用参考信号减去实验信号，就消除了其他效应只保留CEST效应。但是其他干扰效应在参考信号和标记信号中并不都相等，比如大分子磁化转移(Magnetization transfer,MT)效应，所以MTRaysm方法得到的信号包含了其他效应。多池洛伦兹拟合假设每个效应的曲线是一条洛伦兹曲线，这些洛伦兹曲线叠加起来形成了实验采集的曲线，通过拟合实验采集到的曲线，可以得出各个效应的值，但是其他效应的曲线并不都是洛伦兹曲线，比如MT，它的曲线是超洛伦兹。外插磁化转移信号(Extrapolated Semisolid Magnetization Transfer Reference,EMR)拟合出的曲线代表除了CEST效应之外的其他效应，将EMR拟合的曲线与实验得到的曲线作差，就可以得到CEST效应。EMR假设实验使用的饱和脉冲是方波而且系统已经到达稳态，但是这两个条件并不总是满足的。
发明内容
现有技术中的模型拟合方法对实验条件作了近似，但是这种近似在许多实验条件下并不成立。因此，本发明的目的在于解决现有技术中存在的问题，并在去除CEST成像的干扰效应方面，提供一种鲁棒的、适合于多种实验条件的基于外插磁化转移信号的CEST数值拟合(Numerical Fit of Extrapolated Semisolid Magnetization Transfer ReferenceSignal,简称NEMR)方法。
本发明所采用的具体技术方案如下：
第一方面，本发明提供了一种基于外插磁化转移信号的CEST数据拟合方法，该方法针对CEST图像中的每一个体素，依次执行S1～S3，得到每一个体素的消除背景干扰信号的CEST效应值；其中S1～S3为：
S1、对获得的原始Z谱进行主磁场频率B0偏移校正，得到校正后Z谱；
S2、取校正后Z谱中仅含有MT效应的部分数据点作为拟合样本数据，对预先设置有待定参数初值和上下限的两池模型进行拟合，得到待定参数的拟合值；
所述两池模型表示为
Figure BDA0003047176270000021
式中：exp表示以自然常数e为底的指数函数；Δt表示时刻t的增量；
Figure BDA0003047176270000022
分别为t+Δt时刻、t时刻的磁化矢量
Figure BDA0003047176270000023
其中
Figure BDA0003047176270000024
矩阵M＝[Mxa,Mya,Mza,Mzb]T，Mxa、Mya、Mza、Mzb分别为水集团的x方向磁化强度分量、水集团的y方向磁化强度分量、水集团的z方向磁化强度分量、MT集团的z方向磁化强度分量；系数矩阵
Figure BDA0003047176270000025
其中矩阵
Figure BDA0003047176270000026
向量C＝[0,0,R1aM0a,R1bM0b]T，Δωa表示外界施加的射频脉冲的频率和水的频率差，ω1表示所施加的射频脉冲强度，kab表示从水集团交换到MT集团的速率，kba表示从MT集团交换到水集团的速率，Rrfb为MT集团对射频脉冲的吸收速率，Δωb表示外界施加的射频脉冲的频率和MT的频率差，R1a＝1/T1a和R2a＝1/T2a分别为水集团的纵向弛豫速率和横向弛豫速率，R1b＝1/T1b和R2b＝1/T2b分别为MT集团的纵向弛豫速率和横向弛豫速率，M0a和M0b分别为水集团和MT集团的稳态磁化强度；
所述两池模型中的待定参数为T1a,T2a,T2b,M0b,kba；
S3、基于代入所述拟合值的两池模型得到包含MT效应的拟合曲线，将所述拟合曲线与所述校正后Z谱上感兴趣频率处的差值，作为消去背景干扰信号的CEST效应值。
作为第一方面的优选，所述S2中，两池模型中待定参数的拟合具有两轮，第一轮拟合时对5个待定参数在其预先设定的上下限范围内均利用所述拟合样本数据进行参数优化，获得第一轮拟合值；第二轮拟合时，以第一轮拟合值为中心，选择相对于第一轮拟合上下限范围更窄的上下限，重新使用和第一轮拟合中相同的拟合样本数据对所述两池模型进行第二轮精度更高的拟合，得到待定参数的最终拟合值。
再进一步的优选是，在进行第二轮拟合时，在所述两池模型的5个待定参数中仅选择部分参数调整其上下限，剩余参数的上下限与第一轮拟合保持一致；
更进一步的优选是，在进行第二轮拟合时，需改变上下限的参数组合为T1a,T2b,M0b或T1a,T2b或T2b,M0b。
作为第一方面中两轮拟合方案的优选，所述S2中，仅含有MT效应的数据点为频率位于80～6ppm范围内的数据点，优选为频率位于80～20ppm范围内的数据点。
作为第一方面中两轮拟合方案的优选，所述S2中，进行第一轮拟合时，待定参数T1a,T2a,T2b,M0b,kba预设的上限是[1.5，0.25，12e-6，0.16，40]，预设的下限是[0.6，0.065，8e-6，0.001，20]。
作为第一方面中两轮拟合方案的优选，所述S2中，在进行第二轮拟合时，对于需改变上下限的参数，其上限值改为该参数的第一轮拟合值的120％～150％，下限值改为该参数的第一轮拟合值的80％～50％；进一步的，上限值优选改为该参数的第一轮拟合值的120％，下限值改为该参数的第一轮拟合值的80％。
作为第一方面中两轮拟合方案的优选，所述S2中，MT集团对射频脉冲的吸收速率
Figure BDA0003047176270000041
其中g(·)是一个线形函数，包括洛伦兹型、高斯型或超洛伦兹型。
进一步的优选是，所述S3中，所述感兴趣频率为3.5ppm或-3.5ppm，3.5ppm的感兴趣频率对应的CEST效应值为APT#值，-3.5ppm的感兴趣频率对应的CEST效应值为NOE#。
第二方面，本发明提供了一种计算机可读存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如第一方面中任一项所述的基于外插磁化转移信号的CEST数据拟合方法。
第三方面，本发明提供了一种基于外插磁化转移信号数值拟合的CEST数据拟合装置，其包括处理器和存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如第一方面中任一项所述的基于外插磁化转移信号的CEST数据拟合方法。
第四方面，本发明提供了一种去除CEST成像干扰效应的磁共振成像装置，其特征在于，包括磁共振扫描器以及控制单元，所述磁共振扫描器用于通过磁共振CEST成像获取CEST图像；所述控制单元能获取所述CEST图像且控制单元中存储有计算机程序，当所述计算机程序被执行时，用于实现如第一方面中任一项所述的基于外插磁化转移信号的CEST数据拟合方法，输出每一个体素消除背景干扰信号的CEST效应值。
本发明相对于现有技术而言，具有以下有益效果：
本发明提出的拟合方法，本质上是对BM方程的分步求解，该方法有两大优点，首先，由于NEMR在时间维度上进行迭代，可以将实际使用的射频脉冲的波形纳入到求解过程中，而不是假定波形为矩形波。其次，NEMR方法采用BM方程仿真来描述CEST效应，而不是经过简化的解析表达式，使得在系统非稳态时也适用，BM方程是描述CEST效应最基本的方程，在此基础上进行仿真和迭代可以准确地还原生理参数值。
附图说明
图1为本发明的拟合方法流程图；
图2为实施例中实验采集的一条曲线和它的拟合曲线图；
图3为实施例中一位缺血性中风患者的T2加权图、弥散加权图、APT#图和NOE#图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步阐述和说明。
本发明提供了一种去除CEST成像干扰效应的方法，称为外插磁化转移信号的数值拟合(Numerical Fit of Extrapolated Semisolid Magnetization Transfer ReferenceSignal,NEMR)。该基于外插磁化转移信号的CEST数据拟合方法将MT效应的线型融入到Bloch-McConnell(BM)方程中，通过对该方程代表的模型的拟合，可以得出MT效应，将拟合得出的曲线与实验采集到的曲线作差，就可以去除干扰效应。
本发明的核心是通过对BM方程的合理近似从而将MT的线型融入到BM方程，形成一个用于进行后续拟合的两池模型。下面先对该两池模型的原理进行详细描述。
磁共振的成像过程可由如下的Bloch方程来描述，
dMx/dt＝-MxR2-(ωw-ω)My
dMy/dt＝(ωw-ω)Mx-MyR2-ω1Mz
dMz/dt＝ω1My-MzR1+M0R1
其中Mx、My和Mz分别为x、y和z三个方向的磁化强度分量，ωw和ω分别为水的磁共振频率和外界施加的射频脉冲的频率，ω1为所施加的射频脉冲的强度，R1＝1/T1和R2＝1/T2分别为纵向和横向弛豫速率，M0表示稳态磁化强度。该方程只描述了水的磁化强度在射频脉冲影响下随时间改变的情况。对于具有两个交换集团(水集团a；MT集团b)的系统，可由以下方程来描述(需注意的是下列公式中下标a代表水集团的相关参数，而下标b代表MT集团的相关参数)：
Figure BDA0003047176270000061
Figure BDA0003047176270000062
Figure BDA0003047176270000063
Figure BDA0003047176270000064
Figure BDA0003047176270000065
Figure BDA0003047176270000066
其中Δωa＝ω-ωa表示外界施加的射频脉冲的频率和水的频率差；kij(i＝a或b；j＝a或b)则表示了i和j两个集团之间的交换速率(例如kab表示从水集团a交换到MT集团b的速率)。
由于MT集团的T2在微秒量级，Mxb,Myb会快速达到稳态，也即
0＝-ΔωbMyb-R2bMxb
0＝ΔωbMxb-R2bMyb-ω1Mzb
Figure BDA0003047176270000067
可以解得
Figure BDA0003047176270000068
如果定义
Figure BDA0003047176270000069
为MT集团对射频脉冲的吸收速率，Myb可以进一步表示为ω1Myb＝-RrfbMzb，那么关于Mzb的微分方程可以写为
Figure BDA00030471762700000610
上述MT集团对射频脉冲的吸收速率Rrfb还可以用一个线型函数来表示，
Figure BDA00030471762700000611
线形函数g(·)可以使用洛伦兹型、高斯型、超洛伦兹型等多种形式。其中：
洛伦兹型为
Figure BDA00030471762700000612
高斯型为：
Figure BDA0003047176270000071
超洛伦兹型为：
Figure BDA0003047176270000072
经过上述简化之后，原有的六个方程变成四个：
Figure BDA0003047176270000073
Figure BDA0003047176270000074
Figure BDA0003047176270000075
Figure BDA0003047176270000076
这四个微分方程可以表示为：
Figure BDA0003047176270000077
其中，M＝[Mxa,Mya,Mza,Mzb]T，
Figure BDA0003047176270000078
C＝[0,0,R1aM0a,R1bM0b]T
进一步可以将其转化为齐次微分方程：
Figure BDA0003047176270000079
其中，
Figure BDA00030471762700000710
上述齐次微分方程的解可以表示为
Figure BDA00030471762700000711
其中
Figure BDA00030471762700000712
为t+Δt时刻的磁化矢量。
上述齐次微分方程的解即为本发明的两池模型，该模型有5个未知数，分别是T1a,T2a,T2b,M0b,kba，其余参数都是已知量。通过该模型的拟合，可以得出包含MT效应的拟合曲线，将此拟合曲线和经过主磁场频率B0偏移校正后的Z谱曲线作差，就可以得到更为干净的CEST效应。
需说明的是，本发明可以用于获取不同的CEST效应值，但由于不同的CEST效应值对应于曲线上的不同频率，因此为了便于表示，本发明中将Z谱上能够反映某种CEST效应的频率称为感兴趣频率。例如，假如需要计算的CEST效应值为APT#(Amide proton transfersignal from NEMR analysis，使用外插磁化转移信号分析得到的酰胺质子转移信号)值，那么其对应的感兴趣频率为3.5ppm；假如需要计算的CEST效应值为NOE#(Numerical Fitof Extrapolated Semisolid Magnetization Transfer Reference Signal，使用外插磁化转移信号分析得到的核奥氏效应信号)，那么其对应的感兴趣频率为-3.5ppm。当然，假如Z谱上具有其他频率也可以反映某种CEST效应，那么也可以作为感兴趣频率，通过计算拟合曲线与校正后Z谱上感兴趣频率处的差值，即可得到该感兴趣频率对应的消去背景干扰信号的CEST效应值。本发明在实现过程中，可以单次计算一种CEST效应值，亦可同时计算多种CEST效应值，对此不作限定。
基于上述原理，本发明提供了两种基于外插磁化转移信号的CEST数据拟合方法。如图1中的(a)和(b)所示，两种做法的区别在于对于两池模型的参数拟合方法不同，第一种方式为单轮拟合，第二种方式为两轮拟合，其余的做法均相同。下面对两种做法分别通过实施例进行展示说明。
实施例1
在本实施例中，如图1中(a)所示展示了采用上述第一种拟合方式的NEMR拟合方法实现流程。该方法需要先获取CEST图像，然后针对CEST图像中的每一个体素，依次执行S1～S3，得到每一个体素的APT#值和NOE#值，从而组成APT#图和NOE#图。下面对于S1～S3的具体实现过程进行展开描述。
S1、对获得的原始Z谱进行主磁场频率B0偏移校正，得到校正后Z谱。主磁场频率B0偏移校正属于现有技术，对此不再赘述。
S2、取校正后Z谱中仅含有MT效应的部分数据点作为拟合样本数据，对预先设置有待定参数初值和上下限的两池模型进行拟合，得到待定参数的拟合值。
如前所述，本发明中的两池模型表示为
Figure BDA0003047176270000081
式中：exp表示以自然常数e为底的指数函数；Δt表示时刻t的增量；
Figure BDA0003047176270000091
分别为t+Δt时刻、t时刻的磁化矢量
Figure BDA0003047176270000092
其中
Figure BDA0003047176270000093
矩阵M＝[Mxa,Mya,Mza,Mzb]T，Mxa、Mya、Mza、Mzb分别为水集团的x方向磁化强度分量、水集团的y方向磁化强度分量、水集团的z方向磁化强度分量、MT集团的z方向磁化强度分量；系数矩阵
Figure BDA0003047176270000094
其中矩阵
Figure BDA0003047176270000095
向量C＝[0,0,R1aM0a,R1bM0b]T，Δωa表示外界施加的射频脉冲的频率和水的频率差，ω1表示所施加的射频脉冲强度，kab表示从水集团交换到MT集团的速率，kba表示从MT集团交换到水集团的速率，参数
Figure BDA0003047176270000096
其中g(·)为超洛伦兹型，Δωb表示外界施加的射频脉冲的频率和MT的频率差，R1a＝1/T1a和R2a＝1/T2a分别为水集团的纵向弛豫速率和横向弛豫速率，R1b＝1/T1b和R2b＝1/T2b分别为MT集团的纵向弛豫速率和横向弛豫速率，M0a和M0b分别为水集团和MT集团的稳态磁化强度。
上述两池模型中，需要对5个待定参数T1a,T2a,T2b,M0b,kba进行后续的参数拟合。仅含有MT效应的部分数据点可以根据实际的原始Z谱进行分析确定，本实施例中优选为频率位于80～20ppm范围内的数据点。
对于第一种方式而言，可以直接利用拟合样本数据对5个待定参数T1a,T2a,T2b,M0b,kba进行拟合，在5个待定参数各自的上下限范围内，获取使整个两池模型方程拟合程度最高的5个待定参数的拟合值。具体的拟合方法可以利用MATLAB、SPSS等软件或者其他现有技术实现，拟合程度可以用最小均方误差来表示。
在进行拟合时，待定参数T1a,T2a,T2b,M0b,kba预设的上下限以及初值可以根据实际情况调整，只要符合各参数的正常取值范围即可。在本实施例中，数T1a,T2a,T2b,M0b,kba预设的上限设置为[1.5，0.25，12e-6，0.16，40]，预设的下限设置为[0.6，0.065，8e-6，0.001，20]，初值均设为各自的下限。
S3、基于代入5个拟合值的两池模型，即可得到包含MT效应的拟合曲线，双池模型中的频率ω即为曲线的横坐标，双池模型中的Mza即为曲线的纵坐标。用该拟合曲线减去主磁场B0频率校正后的Z谱，即可得到一条差值曲线，在这条差值曲线上取横坐标为3.5ppm和-3.5ppm对应的纵坐标值，分别记为APT#值和NOE#值。当然，从简化角度来看，差值曲线不属于必须要获取的，只要能够计算拟合曲线上3.5ppm对应的纵坐标值与主磁场B0频率校正后的Z谱上3.5ppm对应的纵坐标值间的差值，即可作为当前体素的APT#值，同时计算拟合曲线上-3.5ppm对应的纵坐标值与主磁场B0频率校正后的Z谱上-3.5ppm对应的纵坐标值间的差值，即可作为当前体素的NOE#值。
整张CEST图像中每一个体素均获得各自对应的APT#值和NOE#值后，即可形成一张APT#图和一张NOE#图，APT#图中的每个位置的值为该位置体素对应的APT#值，NOE#图中的每个位置的值为该位置体素对应的NOE#值。
为了进一步展示上述方法的技术效果，将上述NEMR拟合方法在一例缺血性中风病人的大脑CEST图像上进行了测试。在本实施例中，采集图像所使用的CEST成像序列包含三个模块：(1)CEST饱和模块，该模块包含10个高斯波形饱和脉冲，每个脉冲持续时间约为100ms，幅度为1μT；(2)频谱预饱和反转恢复压脂模块；(3)快速自旋回波采集模块。
如图1(a)所示，首先对采集到的Z谱数据(如图2中的实验曲线所示，频率范围80～-6ppm)进行主磁场频率B0偏移校正，然后进行上述NEMR拟合。第一轮拟合使用频率为80～20ppm的数据作为拟合样本数据(拟合用到的点如图2中所示)，未知参数T1a,T2a,T2b,M0b,kba上限是[1.5，0.25，12e-6,0.16,40],下限是[0.6，0.065，8e-6,0.001,20]，初值等于下限。最终拟合曲线如图2所示。
本实施例的实验结果如附图3所示，可以看出，一步拟合的NEMR方法拟合得出的APT#图和NOE#图除了能够清楚地显示出其他模态上可以看到的病灶(T2加权图像上的病灶，白色实线箭头)，还能够显示出其他模态上不可见的病灶(NOE#图上的病灶，白色虚线箭头)。
实施例2
在本实施例中，如图1中(b)所示展示了采用上述第二种拟合方式的NEMR方法实现流程。同样的，该方法需要先获取CEST图像，然后针对CEST图像中的每一个体素，依次执行S1～S3，得到每一个体素的APT#值和NOE#值，从而组成APT#图和NOE#图。
本实施例的S1和S3与实施例1做法相同，区别仅在于步骤S2中，两池模型中待定参数的拟合具有两轮：
第一轮拟合与实施例1相同，先对5个待定参数在其预先设定的上下限范围内均利用前述拟合样本数据进行参数优化，获得第一轮拟合值。5个待定参数的初始值也设置为各自的下限。
但是多参数拟合受变量上下限影响大，有些参数会等于给定的边界值，因此本实施例需要通过设置第二步拟合来减轻这一问题。第二轮拟合时，可以以第一轮拟合值为中心，选择相对于第一轮拟合上下限范围更窄的上下限，重新使用和第一轮拟合中相同的拟合样本数据对上述两池模型进行第二轮精度更高的拟合，得到待定参数的最终拟合值。在具体实现过程中，两池模型的5个待定参数不需要全部进行上下限调整，可以在两池模型的5个待定参数中选择部分参数进一步扩大其上下限值，剩余参数的上下限依然与第一轮拟合保持一致。本实施例在第二轮拟合时，需改变上下限值的参数组合为T1a,T2b,M0b。
本实施例中，通过将部分参数(如T1a,T2b,M0b)限制在第一轮结果的附近，可以使得这些参数在一个较小的范围内变动，从而减轻上下限的不当设置对最终结果的影响。
另外，对于需要改变其上下限值的第一轮拟合值，其上下限值的改变幅度不宜过大但也不宜过小。本实施例中，对于需要改变其上下限值的三个参数T1a,T2b,M0b，第二轮拟合所用的参数上限值优选改为该参数第一轮拟合值的120％，第二轮拟合所用的参数下限值优选改为该参数第一轮拟合值的80％，从而在第一轮拟合值附近对这几个参数进行进一步优化。
为了进一步展示上述方法的技术效果，将上述NEMR拟合方法在一例缺血性中风病人的大脑CEST图像上进行了测试。在本实施例中，采集图像所使用的CEST成像序列与实施例1相同。
如图1(b)所示，首先对采集到的数据进行主磁场频率B0偏移校正，然后进行NEMR拟合。第一轮拟合使用频率为80～20ppm的数据作为拟合样本数据，未知的待拟合参数T1a,T2a,T2b,M0b,kba上限是[1.5，0.25，12e-6,0.16,40]，下限是[0.6，0.065，8e-6,0.001,20]，初值等于下限。第二轮拟合时，使用和第一轮拟合中一样的拟合样本数据，将T1a,T2b,M0b的上限设为第一轮得到的T1a,T2b,M0b的120％，下限设为第一轮得到的T1a,T2b,M0b的80％，其余参数的上下限与第一轮拟合保持一致，所有参数的初值设定为该参数值的下限。
本实施例的实验结果如附图3所示，可以看出，两步拟合的NEMR方法拟合得出的APT#图和NOE#图除了能够清楚地显示出其他模态上可以看到的病灶(T2加权图像上的病灶，白色实线箭头)，还能够显示出其他模态上不可见的病灶(NOE#图上的病灶，白色虚线箭头)，且其显示效果比实施例1更优。
需要注意的是，上述两个实施例仅仅是本发明的两个实例，但并非仅限于此。例如，本发明的上述S2中的拟合样本数据的选择频率范围，可以选用其他范围，如80～20ppm，80～6ppm等。同样的，在进行参数拟合时，对于参数的上下限限制，可以根据实际情况调整，上限值可以改为该参数在第一轮得到的第一轮拟合值的120％～150％，下限值可以改为该参数的第一轮拟合值的80％～50％。再另外，其中第二轮拟合时，需改变上下限值的参数组合可以是多样的，例如除了T1a,T2b,M0b之外，还可以选择T1a,T2b或T2b,M0b等等。在进行拟合时，参数的初值设为上下限范围内的合理数值，也可以通过在上下限范围内随机或者均匀设置几个不同的待选值，所得拟合误差最小的待选值即可选定为最终的初值。
另外，在本发明的另一优选实施例中，上述基于外插磁化转移信号的CEST数据拟合方法，可以以计算机程序的形式存储与计算机可读存储介质。当计算机程序被处理器调用并执行时，可以按照上述方法实现基于外插磁化转移信号的CEST数据拟合方法。
同样的，在本发明的另一优选实施例中，还提供了一种基于外插磁化转移信号数值拟合的CEST数据拟合装置，其包括处理器和存储介质，存储介质上存储有计算机程序，当计算机程序被处理器执行时，实现上述基于外插磁化转移信号的CEST数据拟合方法。
计算机可读存储介质一般以存储器硬件形式提供，存储器可以包括随机存取存储器(Random Access Memory，RAM)，也可以包括非易失性存储器(Non-Volatile Memory，NVM)，例如至少一个磁盘存储器。
上述处理程序的处理器可以是通用处理器，包括中央处理器(CentralProcessing Unit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(Digital Signal Processing，DSP)、专用集成电路(Application SpecificIntegrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
当然，随着云服务器的广泛应用，上述软件程序也可以搭载于云平台上，提供相应的服务，因此计算机可读存储介质并不限于本地硬件的形式。
在另一优选实施例中，可以将上述基于外插磁化转移信号的CEST数据拟合方法，以程序的形式集成于磁共振成像装置的控制单元中。磁共振成像装置应当包括常规的磁共振扫描器以及控制单元，磁共振扫描器可采用现有技术实现，属于成熟商用产品，不再赘述。控制单元中除了存储有上述计算机程序之外，还应当具有实现CEST成像所必要的成像序列以及其他软件程序。磁共振扫描器用于通过磁共振CEST成像获取CEST图像，而控制单元能获取磁共振扫描器得到的CEST图像，且控制单元中存储有前述的计算机程序，当计算机程序被执行时，用于实现上述基于外插磁化转移信号的CEST数据拟合方法，输出每一个体素消除背景干扰信号的CEST效应值。
本领域的技术人员应当知道，本发明中所涉及的各模块、功能可以通过电路、其他硬件或者可执行的程序代码来完成，只要能够实现相应功能即可。若采用代码，则代码可存储于存储装置中，并有计算装置中的相应元件执行。本发明的实现便不限制于任何特定的硬件和软件结合。本发明中的各硬件型号均可采用市售产品，可根据实际用户需求进行选择。当然，磁共振CEST成像序列及装置中，也需要配合必要的其他硬件或软件，此处不再赘述。
以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (13)

1.一种基于外插磁化转移信号的CEST数据拟合方法，其特征在于，针对CEST图像中的每一个体素，依次执行S1～S3，得到每一个体素的消除背景干扰信号的CEST效应值；其中S1～S3为：
S1、对获得的原始Z谱进行主磁场频率B0偏移校正，得到校正后Z谱；
S2、取校正后Z谱中仅含有MT效应的部分数据点作为拟合样本数据，对预先设置有待定参数初值和上下限的两池模型进行拟合，得到待定参数的拟合值；
所述两池模型表示为
Figure FDA0003376476950000011
式中：exp表示以自然常数e为底的指数函数；Δt表示时刻t的增量；
Figure FDA0003376476950000012
Figure FDA0003376476950000013
分别为t+Δt时刻、t时刻的磁化矢量
Figure FDA0003376476950000014
其中
Figure FDA0003376476950000015
矩阵M＝[Mxa,Mya,Mza,Mzb]T，Mxa、Mya、Mza、Mzb分别为水集团的x方向磁化强度分量、水集团的y方向磁化强度分量、水集团的z方向磁化强度分量、MT集团的z方向磁化强度分量；系数矩阵
Figure FDA0003376476950000016
其中矩阵
Figure FDA0003376476950000017
向量C＝[0,0,R1aM0a,R1bM0b]T，Δωa表示外界施加的射频脉冲的频率和水的频率差，ω1表示所施加的射频脉冲强度，kab表示从水集团交换到MT集团的速率，kba表示从MT集团交换到水集团的速率，Rrfb为MT集团对射频脉冲的吸收速率，R1a＝1/T1a和R2a＝1/T2a分别为水集团的纵向弛豫速率和横向弛豫速率，R1b＝1/T1b和R2b＝1/T2b分别为MT集团的纵向弛豫速率和横向弛豫速率，M0a和M0b分别为水集团和MT集团的稳态磁化强度；
所述两池模型中的待定参数为T1a,T2a,T2b,M0b,kba；
S3、基于代入所述拟合值的两池模型得到包含MT效应的拟合曲线，将所述拟合曲线与所述校正后Z谱上感兴趣频率处的差值，作为消去背景干扰信号的CEST效应值。
2.如权利要求1所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，两池模型中待定参数的拟合具有两轮，第一轮拟合时对5个待定参数在其预先设定的上下限范围内均利用所述拟合样本数据进行参数优化，获得第一轮拟合值；第二轮拟合时，以第一轮拟合值为中心，选择相对于第一轮拟合上下限范围更窄的上下限，重新使用和第一轮拟合中相同的拟合样本数据对所述两池模型进行第二轮精度更高的拟合，得到待定参数的最终拟合值。
3.如权利要求2所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，在进行所述第二轮拟合时，在所述两池模型的5个待定参数中仅选择部分参数调整其上下限，剩余参数的上下限与第一轮拟合保持一致。
4.如权利要求2所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，在进行所述第二轮拟合时，需改变上下限的参数组合为T1a,T2b,M0b或T1a,T2b或T2b,M0b。
5.如权利要求1所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，仅含有MT效应的数据点为频率位于6～80ppm范围内的数据点。
6.如权利要求5所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，仅含有MT效应的数据点为频率位于20～80ppm范围内的数据点。
7.如权利要求2所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，进行第一轮拟合时，待定参数T1a,T2a,T2b,M0b,kba预设的上限是[1.5，0.25，12e-6，0.16，40]，预设的下限是[0.6，0.065，8e-6，0.001，20]。
8.如权利要求2所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，在进行第二轮拟合时，对于需改变上下限的参数，其上限值改为该参数的第一轮拟合值的120％～150％，下限值改为该参数的第一轮拟合值的50％～80％。
9.如权利要求2所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S2中，MT集团对射频脉冲的吸收速率
Figure FDA0003376476950000031
其中g(·)是一个线形函数，包括洛伦兹型、高斯型或超洛伦兹型，Δωb表示外界施加的射频脉冲的频率和MT的频率差。
10.如权利要求9所述的基于外插磁化转移信号的CEST数据拟合方法，其特征在于，所述S3中，所述感兴趣频率为3.5ppm或-3.5ppm，3.5ppm的感兴趣频率对应的CEST效应值为APT#值，-3.5ppm的感兴趣频率对应的CEST效应值为NOE#。
11.一种计算机可读存储介质，其特征在于，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如权利要求1～10任一项所述的基于外插磁化转移信号的CEST数据拟合方法。
12.一种基于外插磁化转移信号数值拟合的CEST数据拟合装置，其特征在于，包括处理器和存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现如权利要求1～10任一项所述的基于外插磁化转移信号的CEST数据拟合方法。
13.一种去除CEST成像干扰效应的磁共振成像装置，其特征在于，包括磁共振扫描器以及控制单元，所述磁共振扫描器用于通过磁共振CEST成像获取CEST图像；所述控制单元能获取所述CEST图像且控制单元中存储有计算机程序，当所述计算机程序被执行时，用于实现如权利要求1～10任一项所述的基于外插磁化转移信号的CEST数据拟合方法，输出每一个体素消除背景干扰信号的CEST效应值。
CN202110475243.0A 2021-04-29 2021-04-29 基于外插磁化转移信号的cest数据拟合方法、装置及介质 Active CN113171076B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202110475243.0A CN113171076B (zh) 2021-04-29 2021-04-29 基于外插磁化转移信号的cest数据拟合方法、装置及介质
JP2023566736A JP2024522017A (ja) 2021-04-29 2021-12-22 外挿磁化移動信号に基づくｃｅｓｔデータフィッティング方法、装置及び媒体
PCT/CN2021/140311 WO2022227639A1 (zh) 2021-04-29 2021-12-22 基于外插磁化转移信号的cest数据拟合方法、装置及介质
US18/495,715 US20240053423A1 (en) 2021-04-29 2023-10-26 Cest data fitting method and apparatus based on extrapolated semisolid magnetization transfer reference signal, and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110475243.0A CN113171076B (zh) 2021-04-29 2021-04-29 基于外插磁化转移信号的cest数据拟合方法、装置及介质

Publications (2)

Publication Number Publication Date
CN113171076A CN113171076A (zh) 2021-07-27
CN113171076B true CN113171076B (zh) 2022-01-04

Family

ID=76925346

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110475243.0A Active CN113171076B (zh) 2021-04-29 2021-04-29 基于外插磁化转移信号的cest数据拟合方法、装置及介质

Country Status (4)

Country Link
US (1) US20240053423A1 (zh)
JP (1) JP2024522017A (zh)
CN (1) CN113171076B (zh)
WO (1) WO2022227639A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113171076B (zh) * 2021-04-29 2022-01-04 浙江大学 基于外插磁化转移信号的cest数据拟合方法、装置及介质
CN115349847B (zh) * 2022-10-19 2023-01-31 之江实验室 基于分离式定量apt成像的乳腺肿瘤辨别系统
WO2024108545A1 (zh) * 2022-11-25 2024-05-30 中国科学院深圳先进技术研究院 一种用于cest信号精准定量的b1场校正方法及装置
CN116879338A (zh) * 2023-06-12 2023-10-13 汕头大学医学院第二附属医院 一种vdmp-cest结合非线性拟合检测gaba的方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104204839A (zh) * 2012-04-03 2014-12-10 皇家飞利浦有限公司 使用apt对比增强和多回波时间采样的mr成像
CN104997511A (zh) * 2015-06-01 2015-10-28 中国科学院深圳先进技术研究院 用于磁共振化学交换饱和转移成像的cestr测量方法和系统
JP2017086667A (ja) * 2015-11-13 2017-05-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
CN110706199A (zh) * 2019-08-29 2020-01-17 西北大学 基于伪佛克脱线型的z谱拟合方法及cest mri量化方法
CN111413655A (zh) * 2020-03-31 2020-07-14 浙江大学 磁共振cest成像频率漂移校正方法、装置、介质及成像设备
CN111566497A (zh) * 2018-01-11 2020-08-21 皇家飞利浦有限公司 用于化学交换饱和转移mri的基于磁化转移的度量

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8536866B2 (en) * 2007-09-26 2013-09-17 The Johns Hopkins University Frequency referencing for chemical exchange saturation transfer (CEST) MRI
US9709511B2 (en) * 2011-11-28 2017-07-18 New York University Apparatus, system, method and computer-readable medium for isolating chemical exchange saturation transfer contrast from magnetization transfer asymmetry under two-frequency RF irradiation
JP6348449B2 (ja) * 2015-05-14 2018-06-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
EP3304060A4 (en) * 2015-05-29 2019-01-23 The Regents of the University of California PH-WEIGHTED MRI USING QUICK IMAGERY OF AMINE EXCHANGE SATURATION TRANSFER (CEST) OF AMINES
US10775462B2 (en) * 2017-07-05 2020-09-15 The General Hospital Corporation System and method for direct saturation-corrected chemical exchange saturation transfer (DISC-CEST)
US11474174B2 (en) * 2019-08-30 2022-10-18 Regents Of The University Of Minnesota Systems and methods for utilizing periodic irradiation to enhance sensitivity of magnetic resonance detection
CN113171076B (zh) * 2021-04-29 2022-01-04 浙江大学 基于外插磁化转移信号的cest数据拟合方法、装置及介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104204839A (zh) * 2012-04-03 2014-12-10 皇家飞利浦有限公司 使用apt对比增强和多回波时间采样的mr成像
CN104997511A (zh) * 2015-06-01 2015-10-28 中国科学院深圳先进技术研究院 用于磁共振化学交换饱和转移成像的cestr测量方法和系统
JP2017086667A (ja) * 2015-11-13 2017-05-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
CN111566497A (zh) * 2018-01-11 2020-08-21 皇家飞利浦有限公司 用于化学交换饱和转移mri的基于磁化转移的度量
CN110706199A (zh) * 2019-08-29 2020-01-17 西北大学 基于伪佛克脱线型的z谱拟合方法及cest mri量化方法
CN111413655A (zh) * 2020-03-31 2020-07-14 浙江大学 磁共振cest成像频率漂移校正方法、装置、介质及成像设备

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Influences of experimental parameters on chemical exchange saturation transfer (CEST) metrics of brain tumors using animal models at 4.7T;Heo, Hye-Young et al;《Magnetic Resonance In Medicine》;20200101;第81卷(第1期);全文 *
Quantifying amide proton exchange rate and concentration in chemical exchange saturation transfer imaging of the human brain;Heo, Hye-Young et al;《NeuroImage》;20200401;第189卷;全文 *
Simplified and scalable numerical solution for describing multi-pool chemical exchange saturation transfer (CEST) MRI contrast;Phillip Zhe Sun;《Journal of magnetic resonance》;20100510;第205卷(第2期);全文 *
化学交换饱和转移磁共振造影方法与应用研究;袁亚平;《中国博士学位论文全文数据库 基础科学辑》;20210215(第2期);全文 *
基于CEST理论模型的优化采样方案研究;梁立;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20160315(第3期);全文 *
脑化学交换饱和转移成像技术与临床应用;杨永贵;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20200215(第2期);全文 *

Also Published As

Publication number Publication date
JP2024522017A (ja) 2024-06-07
WO2022227639A1 (zh) 2022-11-03
CN113171076A (zh) 2021-07-27
US20240053423A1 (en) 2024-02-15

Similar Documents

Publication Publication Date Title
CN113171076B (zh) 基于外插磁化转移信号的cest数据拟合方法、装置及介质
Langkammer et al. Quantitative susceptibility mapping: report from the 2016 reconstruction challenge
Buonincontri et al. MR fingerprinting with simultaneous B1 estimation
Schweser et al. An illustrated comparison of processing methods for phase MRI and QSM: removal of background field contributions from sources outside the region of interest
De Beer et al. Analysis of NMR data using time domain fitting procedures
US10143389B2 (en) Distinguishing diseased tissue from healthy tissue based on tissue component fractions using magnetic resonance fingerprinting (MRF)
US10663545B2 (en) Method and apparatus for low-artifact magnetic resonance fingerprinting scan
Manjón et al. Diffusion weighted image denoising using overcomplete local PCA
US10698055B2 (en) Method, neural network, and magnetic resonance apparatus for assigning magnetic resonance fingerprints
Rodgers et al. Coil combination for receive array spectroscopy: Are data‐driven methods superior to methods using computed field maps?
Robison et al. Correction of B0 eddy current effects in spiral MRI
US10145917B2 (en) Multi-component voxel separation using magnetic resonance fingerprinting with compartment exchange
CN110573897B (zh) 亚体素分辨率磁共振指纹成像
CN110807492B (zh) 一种磁共振多参数同时定量成像方法及系统
Dayan et al. Tractography of the optic radiation: a repeatability and reproducibility study
US11137467B2 (en) System and method for improved diffusion-weighted imaging
Rosenblatt et al. Revisiting multi-subject random effects in fMRI: Advocating prevalence estimation
Positano et al. Fast generation of T2⁎ maps in the entire range of clinical interest: Application to thalassemia major patients
Borkowski et al. Assessment of the systematic errors caused by diffusion gradient inhomogeneity in DTI‐computer simulations
Yu et al. Delay-period activity in frontal, parietal, and occipital cortex tracks noise and biases in visual working memory
Frank et al. Dynamic multiscale modes of resting state brain activity detected by entropy field decomposition
Wang et al. Stochastic optimization of three‐dimensional non‐Cartesian sampling trajectory
Kuczera et al. Optimized bias and signal inference in diffusion‐weighted image analysis (OBSIDIAN)
Zhang et al. Multi‐task convolutional neural network‐based design of radio frequency pulse and the accompanying gradients for magnetic resonance imaging
Shi et al. An improved multi-objective optimization-based CICA method with data-driver temporal reference for group fMRI data analysis

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