CN116740066B - 基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 - Google Patents
基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 Download PDFInfo
- Publication number
- CN116740066B CN116740066B CN202311014862.5A CN202311014862A CN116740066B CN 116740066 B CN116740066 B CN 116740066B CN 202311014862 A CN202311014862 A CN 202311014862A CN 116740066 B CN116740066 B CN 116740066B
- Authority
- CN
- China
- Prior art keywords
- oef
- cbf
- brain
- image
- cmro2
- 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
Links
- 210000004556 brain Anatomy 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000008557 oxygen metabolism Effects 0.000 title claims abstract description 25
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 title claims abstract description 22
- 239000001301 oxygen Substances 0.000 title claims abstract description 22
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 33
- 230000000284 resting effect Effects 0.000 claims abstract description 9
- 239000003550 marker Substances 0.000 claims abstract description 7
- 210000004884 grey matter Anatomy 0.000 claims abstract description 5
- 230000000694 effects Effects 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 9
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 8
- 239000008280 blood Substances 0.000 claims description 8
- 210000004369 blood Anatomy 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 8
- 210000001519 tissue Anatomy 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 230000001419 dependent effect Effects 0.000 claims description 7
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 claims description 6
- 241000257303 Hymenoptera Species 0.000 claims description 6
- 108010002255 deoxyhemoglobin Proteins 0.000 claims description 6
- 230000010412 perfusion Effects 0.000 claims description 6
- 238000010561 standard procedure Methods 0.000 claims description 6
- 238000000844 transformation Methods 0.000 claims description 6
- 210000001175 cerebrospinal fluid Anatomy 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 108010054147 Hemoglobins Proteins 0.000 claims description 4
- 102000001554 Hemoglobins Human genes 0.000 claims description 4
- 101100045546 Arabidopsis thaliana TFCE gene Proteins 0.000 claims description 3
- 238000000692 Student's t-test Methods 0.000 claims description 3
- 210000004727 amygdala Anatomy 0.000 claims description 3
- 210000004227 basal ganglia Anatomy 0.000 claims description 3
- 230000002146 bilateral effect Effects 0.000 claims description 3
- 210000003710 cerebral cortex Anatomy 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 210000001320 hippocampus Anatomy 0.000 claims description 3
- 238000006213 oxygenation reaction Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000013179 statistical model Methods 0.000 claims description 3
- 238000012353 t test Methods 0.000 claims description 3
- 210000004885 white matter Anatomy 0.000 claims description 3
- 102000006386 Myelin Proteins Human genes 0.000 claims description 2
- 108010083674 Myelin Proteins Proteins 0.000 claims description 2
- 210000004204 blood vessel Anatomy 0.000 claims description 2
- 230000005821 brain abnormality Effects 0.000 claims description 2
- 230000001054 cortical effect Effects 0.000 claims description 2
- 230000005292 diamagnetic effect Effects 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 210000001652 frontal lobe Anatomy 0.000 claims description 2
- 238000005534 hematocrit Methods 0.000 claims description 2
- 210000005012 myelin Anatomy 0.000 claims description 2
- 230000008520 organization Effects 0.000 claims description 2
- 230000005298 paramagnetic effect Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims description 2
- 210000003478 temporal lobe Anatomy 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 claims 2
- 239000004484 Briquette Substances 0.000 claims 2
- 238000012360 testing method Methods 0.000 claims 2
- 230000008859 change Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000009499 grossing Methods 0.000 abstract 2
- 230000003727 cerebral blood flow Effects 0.000 description 47
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000002600 positron emission tomography Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000019771 cognition Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007658 neurological function Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012307 MRI technique Methods 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000008455 cerebrovascular function Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000005291 magnetic effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000626 neurodegenerative effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001558 permutation test Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008844 regulatory mechanism Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 230000004218 vascular function Effects 0.000 description 1
Classifications
-
- 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
- A61B5/004—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 adapted for image acquisition of a particular organ or body part
- A61B5/0042—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 adapted for image acquisition of a particular organ or body part for the brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/37—Determination of transform parameters for the alignment of images, i.e. image registration using transform domain methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Quality & Reliability (AREA)
- Neurology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
基于QSM评估脑氧摄取分数及脑氧代谢率的方法及装置，通过获取受试者的三维T1加权快速梯度回波3D‑FSPGR序列，静息态三维伪连续动脉自旋标记3D‑PCASL序列，三维梯度多回波3D‑mGRE序列。CBF从3D‑PCASL中获取，OEF从3D‑mGRE中重建获得。CBF及OEF图像刚性对齐到各自T1空间，根据Fick原理计算CMRO2；CBF、OEF和CMRO2图像归一化到MNI标准空间，平滑处理用于全脑体素分析，未平滑处理参数图像用于区域分析，能够探索OEF、CBF和CMRO2的全脑改变模式，找到病例组和对照组之间深灰质区域的OEF、CBF和CMRO2的差异，研究氧代谢与临床特征之间的关系。
Description
技术领域
本发明涉及医学图像处理的技术领域，尤其涉及一种基于QSM评估脑氧摄取分数及脑氧代谢率的方法，以及基于QSM评估脑氧摄取分数及脑氧代谢率的装置。
背景技术
大脑氧利用、血液供应和能量消耗整体状态的几个生理变量包括氧摄取分数（OEF）、脑血流（CBF）和脑氧代谢率（CMRO2）。这些参数提供了关于脑组织代谢水平和血液动力学的综合生理信息，直接反映了脑内神经和血管功能以及脑平衡血管反应性和神经功能的能力。因此，这些参数是神经退行性基本中的潜在感兴趣的生物标志物。大多数研究显示患者皮质区域CBF减少。然而， OEF值和CMRO2的变化在不同研究中各异，可以显示出增加、减少或无显著变化的混合模式。这些混合结果暗示了复杂的脑血管功能和氧代谢调节机制。因此，在全脑整体和局部区域水平上综合评估上述参数是一个极有吸引力的研究目标。
传统上，使用正电子发射断层扫描（PET）来测量OEF、CBF和CMRO2，被广泛认为是定量成像OEF和CMRO2的金标准。然而，PET OEF和CMRO2成像需要复杂的过程，并且放射性同位素示踪剂15O具有短半衰期（2分钟），对患者不可避免地会有辐射暴露，这在很大程度上限制了15O PET的临床应用。随着磁共振成像技术的进步，提出了各种定量技术，可以基于MR幅度和相位信号与血氧化相关性来估算OEF。一些代表性的方法包括双校准fMRI、定量血氧水平依赖性成像（qBOLD）、基于T2的方法和定量磁化率映像（QSM）方法。我们提出了一种基于QSM的方法，名为QSM+qBOLD（QQ）模型，使用来自多回波梯度回波（mGRE）序列的相位和幅度信号的综合信息。通过时间演化和组织成分算法（QQ-CCTV）的聚类分析，该方法是一种准确评估OEF的便捷和可行方法，具有改善的信噪比。由于QQ方法已被证明与双气体校准BOLD和PET测量相媲美，基于QQ的OEF映像是临床实践中测量OEF的重要工具。
发明内容
为克服现有技术的缺陷，本发明要解决的技术问题是提供了一种基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其能够探索AD中OEF、CBF和CMRO2的整体改变模式，找到病例组和对照组之间深灰质区域的OEF、CBF和CMRO2值的差异，并进一步研究了氧代谢与临床特征之间的关系。
本发明的技术方案是：这种基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其包括以下步骤：
（1）获取受试者的三维T1加权快速梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度回波多回波3D-mGRE序列；
（2）CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN算法重新采样到T1加权图像的空间分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖的相位信号和qBOLD模型中mGRE数据的幅度信号；
（3）所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析。
本发明通过获取受试者的三维T1加权快速梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度多回波3D-mGRE序列，CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间，所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析，从而能够探索AD中OEF、CBF和CMRO2的整体改变模式，找到病例组和对照组之间深灰质区域的OEF、CBF和CMRO2值的差异，并进一步研究了氧代谢与综合认知之间的关系。
还提供了一种基于QSM评估脑氧摄取分数及脑氧代谢率的装置，其包括：
数据获取模块，其配置来获取受试者的三维T1加权快速破裂梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度多回波3D-mGRE序列；
归一化模块，其配置来CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN算法重新采样到T1加权图像的空间分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖的相位信号和qBOLD模型中mGRE数据的幅度信号；
分析模块，其配置来所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析。
附图说明
图1示出了根据本发明的基于QSM评估脑氧摄取分数及脑氧代谢率的方法的流程图。
具体实施方式
如图1所示，这种基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其包括以下步骤：
（1）获取受试者的三维T1加权快速梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度回波多回波3D-mGRE序列；
（2）CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN例程重新采样到T1加权图像的分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖相位信号和qBOLD模型中mGRE数据的幅度信号；
（3）所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析。
本发明通过获取受试者的三维T1加权快速破裂梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度回波多回波3D-mGRE序列，CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间，所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析，从而能够探索AD中OEF、CBF和CMRO2的整体改变模式，找到AD和对照组之间深灰质区域的OEF、CBF和CMRO2值的差异，并进一步研究了氧代谢与综合认知之间的关系。
优选地，所述步骤（1）中，扫描参数如下：3D-FSPGR：回波时间TE= 2.9毫秒；重复时间TR= 6.7毫秒；翻转角度= 12°；带宽= 31.25 Hz/像素；切片厚度= 1.0毫米；视野FOV=256毫米× 256毫米；体素尺寸= 1× 1× 1 mm³；3D-mGRE：回波时间TE：TE1st/ΔTE/TE8th = 3.19毫秒/2.37毫秒/19.77毫秒；TR = 22.9毫秒；带宽= 62.5 Hz/像素；切片厚度= 1.0毫米；FOV = 256毫米× 256毫米；体素尺寸= 1× 1× 1 mm³；3D-PCASL：TR = 4,817毫秒；TE = 14.6毫秒；PLD = 1,525毫秒；读出中的螺旋线有12个臂，包含1024个采样点；切片厚度= 4毫米；FOV = 240毫米× 240毫米；体素尺寸= 1.875× 1.875× 4 mm³；轴向T2WI、T2-FLAIR和DWI的常规MR序列也包括在内以检测脑部异常。
优选地，所述步骤（2）中，首先通过对mGRE相位数据进行非线性拟合来估计总场图像，然后使用FMRIB软件库版本5.0.9中的脑提取工具根据平均幅度图像生成二值脑掩膜；使用Laplacian边界值LBV方法和脑掩膜计算局部场图像以去除背景频率，并且通过LBV和Spherical mean value方法部分修正背景B0非均匀性；最后，应用MEDI方法和自动统一的脑脊液CSF零参考MEDI+0算法估计QSM图像。
优选地，所述步骤（2）中，qBOLD模型区分由柱状磁敏源和扩散磁敏源引起的独特信号衰减模式，包括R2效应横向弛豫率、气体-组织界面效应以及由脱氧血液体积v、非血液组织磁敏度x_nb和静脉氧合度Yv引起的FBOLD效应；而QSM模型考虑组织磁敏度的空间变化，并部分将强抗磁性的柱状髓鞘束与顺磁性的脱氧血液分离，其中包括x_nb、血浆磁敏度和由v和Yv确定的血红蛋白磁敏度；通过结合qBOLD模型和QSM模型这两个模型，使用mGRE幅度和QSM进行迭代求解来估计OEF，v、x_nb、R2效应和初始信号强度；使用体素传播函数方法校正mGRE数据中的B0非均匀性；为了提高QQ模型基于QEF对噪声的鲁棒性和OEF的准确性，使用时态聚类、组织成分和总变差算法。
优选地，所述步骤（2）中，通过使用ANTS，将经射频偏差校正N4ITK的第一回波幅度图像与其T1图像进行仿射配准；OEF图像通过ANTS中的antsApplyTransforms算法沿上述变换被进一步映射到其对应的解剖图像中的本地空间；根据Fick原理，根据以下公式计算CMRO2图像：
CMRO2 = CBF × OEF × [H]a
[H]a = 7.377μmol/ml，是在红细胞比容Hct= 0.357的动脉血管中氧化血红蛋白的摩尔浓度。
优选地，所述步骤（2）中，个体的T1结构图像通过ANTs中的可变形b值共配准算法，被映射到MNI标准空间，这些变形用于将原生空间的CBF、OEF和CMRO2图像转换到MNI标准空间，使用BNA246图谱，提取246个脑区中CBF、OEF和CMRO2的组平均值，并使用BrainNetViewer工具箱将其投影到一个脑表面模板上。
优选地，所述步骤（3）中，T1解剖图像使用Statistical Parametric Mapping版本12软件分割为灰质和白质的概率图，这些概率图通过上述变形进行空间归一化到MNI空间；最终的二进制脑掩模在组织概率上使用0.6的阈值生成，排除非实质性体素，最终的二进制脑掩模用于限制归一化的OEF、CBF和CMRO2图像的实质边界，并更准确地消除掉脑脊液区域；为了提高统计效率并减少配准不准确的影响，归一化的CBF、OEF和CMRO2图像使用标准差为3 mm的3D高斯核进行平滑处理，并随后使用最终的脑掩模进行平滑补偿；
在FSL中建立CBF、OEF和CMRO2组间比较的一般线性模型，以整个脑体素水平为基础；年龄和性别都被去平均处理，并作为这些一般线性统计模型中的协变量，统计区域限制在所有受试者的最终的平均脑掩模中，在FSL中进行10,000次随机非参数置换检验和团块增强TFCE，采用FWE算法进行多重比较校正，显著性统计的阈值设置为α = 0.05；基于AAL图谱，使用SPM的xjView工具箱来确定和报告显著集群的解剖位置。
优选地，所述步骤（3）中，
从MNI空间的BNA246图谱中确定包括双侧杏仁核、海马和基底节等二十个深灰质区域，用于基于感兴趣区域ROI的分析；为减少部分容积效应的影响，使用FSL将灰质ROI进行边缘修正，并与最后的脑掩模进行点乘，脑实质的密度阈值设置为0.5；然后，使用Matlab从所有受试者的数据中自动分割、提取每个ROI的未平滑的平均OEF、CBF和CMRO2值；
采用t检验探索病例组和对照组在每个ROI的CBF、OEF和CMRO2值之间的差异；然后，使用偏相关分析评估MMSE分数与显示显著差异的ROI中的变量之间的相关性，调整年龄；采用Benjamini-Hochberg假发现率（FDR）方法进行多重比较校正，α值为0.05；统计分析和图形创建使用GraphPad Prism 8软件进行。
优选地，所述步骤（3）中，在病例组和对照组中，CBF、OEF和CMRO2在整个大脑皮层的分布模式在视觉上相似；在病例组和对照组中，CBF和CMRO2较高的脑区位于额叶和颞叶，而OEF图谱相对平坦；AD组整体皮层CBF和CMRO2值明显低于HC组，然而，OEF值在病例组和对照组之间并没有显示出高度差异。
本领域普通技术人员可以理解，实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，所述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，包括上述实施例方法的各步骤，而所述的存储介质可以是：ROM/RAM、磁碟、光盘、存储卡等。因此，与本发明的方法相对应的，本发明还同时包括一种基于QSM评估脑氧摄取分数及脑氧代谢率的装置，该装置通常以与方法各步骤相对应的功能模块的形式表示。该装置包括：
数据获取模块，其配置来获取受试者的三维T1加权快速梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度多回波3D-mGRE序列；
归一化模块，其配置来CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN例程重新采样到T1加权图像的分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖相位信号和qBOLD模型中mGRE数据的幅度信号；
分析模块，其配置来所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属本发明技术方案的保护范围。
Claims (8)
1.基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：其包括以下步骤：
（1）获取受试者的三维T1加权快速梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度回波多回波3D-mGRE序列；
（2）CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN算法重新采样到T1加权图像的空间分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖的相位信号和qBOLD模型中mGRE数据的幅度信号；
（3）所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析；
所述步骤（3）中，T1解剖图像使用Statistical Parametric Mapping版本12软件分割为灰质和白质的概率图，将这些概率图空间归一化到MNI空间；最终的二进制脑掩模在组织概率上使用0.6的阈值生成，排除非实质性体素，最终的二进制脑掩模用于限制归一化的OEF、CBF和CMRO2图像的实质边界，并更准确地消除掉脑脊液区域；为了增强统计效率并减少配准不准确的影响，归一化的CBF、OEF和CMRO2图像使用标准差为3 mm的3D高斯核进行平滑处理，并随后使用最终的脑掩模进行平滑补偿；
在FSL中建立CBF、OEF和CMRO2组间比较的一般线性模型，以整个脑体素水平为基础；年龄和性别都去平均化处理，并作为一般线性统计模型中的协变量，统计区域限制在所有受试者的最终平均脑掩模内，在FSL中进行10,000次的随机非参数置换检验和团块增强TFCE算法，采用FWE方法进行多重比较校正，显著性统计的阈值设置为α = 0.05；基于AAL图谱，使用SPM的xjView工具箱来确定和报告显著集群的解剖位置；
从MNI空间的BNA246图谱中确定包括双侧杏仁核、海马和基底节的二十个深灰质区域，用于基于感兴趣区域ROI的分析；为减少部分容积效应的影响，使用FSL将灰质ROI进行边缘修正，并与最后的脑掩模进行点乘，脑实质的密度阈值设置为0.5；然后，使用Matlab软件从所有受试者的数据中自动分割、提取每个ROI的未平滑的平均OEF、CBF和CMRO2值；
采用t检验探索病例组和对照组在每个ROI的CBF、OEF和CMRO2值之间的差异；然后，使用偏相关分析评估MMSE分数与显示显著差异的ROI中的变量之间的相关性，调整年龄；采用Benjamini-Hochberg假发现率FDR方法进行多重比较校正，α值为0.05；统计分析和图形创建使用GraphPad Prism 8软件进行。
2.根据权利要求1所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（1）中，扫描参数如下：3D-FSPGR：回波时间TE= 2.9毫秒；重复时间TR= 6.7毫秒；翻转角度= 12°；带宽= 31.25 Hz/像素；切片厚度= 1.0毫米；视野FOV= 256毫米× 256毫米；体素尺寸= 1× 1× 1 mm³；3D-mGRE：回波时间TE：TE1st/ΔTE/ TE8th = 3.19毫秒/2.37毫秒/19.77毫秒；TR = 22.9毫秒；带宽= 62.5 Hz/像素；切片厚度= 1.0毫米；FOV =256毫米× 256毫米；体素尺寸= 1× 1× 1 mm³；3D-PCASL：TR = 4,817毫秒；TE = 14.6毫秒；PLD = 1,525毫秒；读出中的螺旋线有12个臂，包含1024个采样点；切片厚度= 4毫米；FOV = 240毫米× 240毫米；体素尺寸= 1.875× 1.875× 4 mm³；轴向T2WI、T2-FLAIR和DWI的常规MR序列也包括在内以检测脑部异常。
3.根据权利要求2所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（2）中，首先通过对mGRE相位数据进行非线性拟合来估计总场图像，然后使用FMRIB软件库版本5.0.9中的脑提取工具根据平均幅度图像生成二值脑掩膜；使用Laplacian边界值LBV方法和脑掩膜计算局部场图像以去除背景频率，并且通过LBV和Spherical mean value方法部分修正背景B0非均匀性；最后，应用MEDI方法和自动统一的脑脊液CSF零参考MEDI+0算法估计QSM图像。
4.根据权利要求3所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（2）中，qBOLD模型区分由柱状磁敏源和扩散磁敏源引起的独特信号衰减模式，包括R2效应横向弛豫率、气体-组织界面效应以及由脱氧血液体积v、非血液组织磁敏度x_nb和静脉氧合度Yv引起的FBOLD效应；而QSM模型考虑组织磁敏度的空间变化，并部分将强抗磁性的柱状髓鞘束与顺磁性的脱氧血液分离，其中包括x_nb、血浆磁敏度和由v和Yv确定的血红蛋白磁敏度；通过结合qBOLD模型和QSM模型这两个模型，使用mGRE幅度和QSM进行迭代求解来估计OEF，v、x_nb、R2效应和初始信号强度；使用体素传播函数方法校正mGRE数据中的B0非均匀性；为了提高QQ模型基于QEF对噪声的鲁棒性和OEF的准确性，使用时态聚类、组织成分和总变差算法。
5.根据权利要求4所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（2）中，通过使用ANTS，将经射频偏差校正N4ITK的第一回波幅度图像与其T1图像进行仿射配准；OEF图像通过ANTS中的antsApplyTransforms算法沿上述变换被进一步映射到其对应的解剖图像中的本地空间；根据Fick原理，根据以下公式计算CMRO2图像：
CMRO2 = CBF × OEF × [H]a
[H]a = 7.377μmol/ml，是在红细胞比容Hct= 0.357的动脉血管中氧化血红蛋白的摩尔浓度。
6.根据权利要求5所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（2）中，个体的T1结构图像通过ANTs中的可变形b值共配准算法，被映射到MNI标准空间，这些变形用于将原生空间的CBF、OEF和CMRO2图像转换到MNI标准空间，使用BNA246图谱，提取246个脑区中CBF、OEF和CMRO2的组平均值，并使用BrainNet Viewer工具箱将其投影到一个脑表面模板上。
7.根据权利要求6所述的基于QSM评估脑氧摄取分数及脑氧代谢率的方法，其特征在于：所述步骤（3）中，在病例组和健康对照组中，CBF、OEF和CMRO2在整个大脑皮层的分布模式在视觉上相似；在病例和对照组中，CBF和CMRO2较高的脑区位于额叶和颞叶，而OEF图相对平坦；病例组整体皮层CBF和CMRO2值明显低于对照组，然而，OEF值在病例组和对照组之间并没有显示出高度差异。
8.基于QSM评估脑氧摄取分数及脑氧代谢率的装置，其特征在于：其包括：
数据获取模块，其配置来获取受试者的三维T1加权快速破裂梯度回波3D-FSPGR序列，用于CBF计算的静息态三维伪连续动脉自旋标记3D-PCASL序列，以及用于QSM和OEF分析的三维梯度回波多回波3D-mGRE序列；
归一化模块，其配置来CBF和OEF图像经过刚性对齐到各自的T1空间，然后根据Fick原理计算出CMRO2；CBF、OEF和CMRO2图像最终被归一化到MNI空间；
CBF图像是使用GE的perfusion FuncTool从PCASL数据中重建的，通过从控制图像中减去标记图像，计算出ASL差异图像，然后根据ASL差异图像和质子密度加权的参考图像进一步计算得到CBF图像，CBF图像随后使用刚性和仿射变换与其对应的T1解剖图像在本地空间进行配准，并通过Advanced Normalization Tools中的antsRegistrationSyN例程重新采样到T1加权图像的分辨率；
QSM图像是通过MEDI工具箱中的标准流程从3D-mGRE数据生成的；
OEF图像是基于QQ模型估计的，该模型结合了QSM中静脉脱氧血红蛋白依赖的相位信号和qBOLD模型中mGRE数据的幅度信号；
分析模块，其配置来所有参数图像被平滑处理用于全脑体素分析，未经平滑处理的参数图像用于区域分析；
所述分析模块中，T1解剖图像使用Statistical Parametric Mapping版本12软件分割为灰质和白质的概率图，将这些概率图空间归一化到MNI空间；最终的二进制脑掩模在组织概率上使用0.6的阈值生成，排除非实质性体素，最终的二进制脑掩模用于限制归一化的OEF、CBF和CMRO2图像的实质边界，并更准确地消除掉脑脊液区域；为了增强统计效率并减少配准不准确的影响，归一化的CBF、OEF和CMRO2图像使用标准差为3 mm的3D高斯核进行平滑处理，并随后使用最终的脑掩模进行平滑补偿；
在FSL中建立CBF、OEF和CMRO2组间比较的一般线性模型，以整个脑体素水平为基础；年龄和性别都去平均化处理，并作为一般线性统计模型中的协变量，统计区域限制在所有受试者的最终平均脑掩模内，在FSL中进行10,000次的随机非参数置换检验和团块增强TFCE算法，采用FWE方法进行多重比较校正，显著性统计的阈值设置为α = 0.05；基于AAL图谱，使用SPM的xjView工具箱来确定和报告显著集群的解剖位置；
从MNI空间的BNA246图谱中确定包括双侧杏仁核、海马和基底节的二十个深灰质区域，用于基于感兴趣区域ROI的分析；为减少部分容积效应的影响，使用FSL将灰质ROI进行边缘修正，并与最后的脑掩模进行点乘，脑实质的密度阈值设置为0.5；然后，使用Matlab软件从所有受试者的数据中自动分割、提取每个ROI的未平滑的平均OEF、CBF和CMRO2值；
采用t检验探索病例组和对照组在每个ROI的CBF、OEF和CMRO2值之间的差异；然后，使用偏相关分析评估MMSE分数与显示显著差异的ROI中的变量之间的相关性，调整年龄；采用Benjamini-Hochberg假发现率FDR方法进行多重比较校正，α值为0.05；统计分析和图形创建使用GraphPad Prism 8软件进行。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311014862.5A CN116740066B (zh) | 2023-08-14 | 2023-08-14 | 基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311014862.5A CN116740066B (zh) | 2023-08-14 | 2023-08-14 | 基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116740066A CN116740066A (zh) | 2023-09-12 |
| CN116740066B true CN116740066B (zh) | 2023-10-27 |
Family
ID=87902947
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311014862.5A Active CN116740066B (zh) | 2023-08-14 | 2023-08-14 | 基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116740066B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103957778A (zh) * | 2012-11-27 | 2014-07-30 | 株式会社东芝 | 阶段判定支援系统 |
| WO2016094461A1 (en) * | 2014-12-09 | 2016-06-16 | The University Of North Carolina At Chapel Hill | Methods, systems, and computer readable media for estimating the index of cerebral oxygen metabolism (i-com) using magnetic resonance (mr) imaging |
| CN107613915A (zh) * | 2015-01-15 | 2018-01-19 | 海博亚医疗有限责任公司 | 根据热生成估算器官的局部代谢率及估算用于灌注组织体积的局部血流量的装置 |
| CN108693491A (zh) * | 2017-04-07 | 2018-10-23 | 康奈尔大学 | 稳健的定量磁化率成像系统和方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012130249A1 (en) * | 2011-03-31 | 2012-10-04 | Aarhus Universitet | System for assessing tissue substance extraction |
| US20210338092A1 (en) * | 2019-04-01 | 2021-11-04 | The Regents Of The University Of California | Portable device for quantitative measurement of tissue autoregulation and neurovascular coupling using eeg, metabolism, and blood flow diagnostics |
-
2023
- 2023-08-14 CN CN202311014862.5A patent/CN116740066B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103957778A (zh) * | 2012-11-27 | 2014-07-30 | 株式会社东芝 | 阶段判定支援系统 |
| WO2016094461A1 (en) * | 2014-12-09 | 2016-06-16 | The University Of North Carolina At Chapel Hill | Methods, systems, and computer readable media for estimating the index of cerebral oxygen metabolism (i-com) using magnetic resonance (mr) imaging |
| CN107613915A (zh) * | 2015-01-15 | 2018-01-19 | 海博亚医疗有限责任公司 | 根据热生成估算器官的局部代谢率及估算用于灌注组织体积的局部血流量的装置 |
| CN108693491A (zh) * | 2017-04-07 | 2018-10-23 | 康奈尔大学 | 稳健的定量磁化率成像系统和方法 |
Non-Patent Citations (2)
| Title |
|---|
| Cerebral oxygen extraction fraction (OEF): Comparison of challenge-free gradient echo QSMþqBOLD (QQ) with 15O PET in healthy adults;Junghun Cho 等;《JCBFM》;全文 * |
| 定量磁化率成像技术对抑郁症患者脑脊液铁沉积的纵向研究;谢杰 等;《放射学实践》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116740066A (zh) | 2023-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zaretskaya et al. | Advantages of cortical surface reconstruction using submillimeter 7 T MEMPRAGE | |
| Gurney-Champion et al. | Quantitative imaging for radiotherapy purposes | |
| CN102525466B (zh) | 图像处理装置和磁共振成像装置 | |
| Henriksen et al. | Estimation of intersubject variability of cerebral blood flow measurements using MRI and positron emission tomography | |
| US9339239B2 (en) | Methods and devices for optimization of magnetic resonance imaging protocols | |
| US7343193B2 (en) | Background suppression method for time-resolved magnetic resonance angiography | |
| US7583992B2 (en) | Magnetic resonance angiography with automated vessel segmentation | |
| US8099149B2 (en) | MRI method for quantification of cerebral perfusion | |
| Bissell et al. | 4D Flow cardiovascular magnetic resonance consensus statement: 2023 update | |
| US20100030062A1 (en) | System and method to analyze blood parameters using magnetic resonance imaging | |
| WO2018083952A1 (ja) | 磁気共鳴イメージング装置および酸素摂取率算出方法 | |
| US8824766B2 (en) | Systems and methods for automated magnetic resonance imaging | |
| US10684339B2 (en) | Dual flip angle multi-echo ultra-short echo time (DUFA-MUTE) magnetic resonance imaging (MRI) systems and methods | |
| Lévy et al. | Intravoxel incoherent motion at 7 tesla to quantify human spinal cord perfusion: limitations and promises | |
| US10846861B2 (en) | Image processor, image processing method, program for image processing, and magnetic resonance imaging apparatus | |
| US20230320611A1 (en) | System and method of accurate quantitative mapping of biophysical parameters from mri data | |
| CN116740066B (zh) | 基于qsm评估脑氧摄取分数及脑氧代谢率的方法及装置 | |
| Bydder | The Mackenzie Davidson Memorial Lecture: Detection of small changes to the brain with serial magnetic resonance imaging | |
| Pandey et al. | Multiresolution imaging using golden angle stack‐of‐stars and compressed sensing for dynamic MR urography | |
| Nomura et al. | Registration method between phase-contrast magnetic resonance angiography and time-of-flight magnetic resonance angiography—a preliminary study | |
| Tatewaki et al. | Regional reliability of quantitative signal targeting with alternating radiofrequency (STAR) labeling of arterial regions (QUASAR) | |
| Oh et al. | Unpaired deep learning for pharmacokinetic parameter estimation from dynamic contrast-enhanced MRI without AIF measurements | |
| Ahmed | Advancements to the reference region model for dynamic contrast enhanced MRI | |
| Yu | Imaging markers of cognitive deficits in RRMS using a multimodal MR approach | |
| Kossen et al. | Redefining Hemodynamic Imaging in Stroke: Perfusion Parameter Map Generation from TOF-MRA using Artificial Intelligence |
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 |