CN107456225A - A kind of cerebral blood flow (CBF) partial volume effect bearing calibration - Google Patents

A kind of cerebral blood flow (CBF) partial volume effect bearing calibration Download PDF

Info

Publication number
CN107456225A
CN107456225A CN201710639032.XA CN201710639032A CN107456225A CN 107456225 A CN107456225 A CN 107456225A CN 201710639032 A CN201710639032 A CN 201710639032A CN 107456225 A CN107456225 A CN 107456225A
Authority
CN
China
Prior art keywords
cbf
brain
value
white
blood flow
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
CN201710639032.XA
Other languages
Chinese (zh)
Other versions
CN107456225B (en
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.)
Peking University Third Hospital
Original Assignee
Peking University Third Hospital
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 Peking University Third Hospital filed Critical Peking University Third Hospital
Priority to CN201710639032.XA priority Critical patent/CN107456225B/en
Publication of CN107456225A publication Critical patent/CN107456225A/en
Application granted granted Critical
Publication of CN107456225B publication Critical patent/CN107456225B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • A61B5/0263Measuring blood flow using NMR
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

The present invention relates to a kind of cerebral blood flow (CBF) partial volume effect bearing calibration, the purpose of the present invention is exactly to carry out partial volume correction using the CBF values obtained in ASL technologies, and accurate algorithm is provided for clear and definite AD diagnosis.Using the CBF values obtained in ASL and T1 values as inputting, the CBF values of full brain white matter of brain, ectocinerea and cerebrospinal fluid are drawn by PVE algorithms.The CBF values carried out according to this method after PVE corrections are divided into two groups of white matter of brain, ectocinerea pixel values.Compared with full brain CBF values, result of calculation is more accurate, is advantageous to the diagnosis and analysis of clinic.

Description

A kind of cerebral blood flow (CBF) partial volume effect bearing calibration
Technical field
The present invention relates to data processing field, and in particular to a kind of cerebral blood flow (CBF) partial volume effect bearing calibration.
Background technology
Cerebral blood flow (CBF)(Cerebral Blood Flow, CBF)It is that a certain amount of brain tissue blood vessel knot is flowed through in the unit interval The CBF of structure.The blood conveying capacity for the tissue being defined as in time per unit to unit mass, unit ml/100g/ min.CBF change is the early sign that brain is activated, and is the important symbol for evaluating brain health status.It is diagnosis and controlled Treat the main of the cerebrovascular disease such as Alzheimer's disease, cerebral infarction, cerebral hemorrhage, aneurysm and congenital artery and cerebral venous malformation Foundation.At present, the mode of brain blood flow is measured mainly by PET, SPECT, MRI technique etc..With nmr imaging technique not It is disconnected to improve, nuclear-magnetism artery rotation flag sequence(Artery Spin Label, ASL)Technology is as a kind of completely non-invasive new blood Infusate flow technology has been increasingly used for CBF measurement.Alzheimer ' is being carried out using the CBF values of magnetic resonance ASL technical limit spacings Silent disease(AD)During research, compared with recognizing normal aging people, because different degrees of encephalatrophy be present especially with Medial Temporal Lobe in patient AD Atrophy is notable, and the degree difference of encephalatrophy will influence whether the registration of image, because the interlamellar spacing of most of ASL images is more For 3 ~ 4mm, the mixing voxel of a large amount of grey matters and white matter be present.If these voxels are calculated using the brain blood flow of grey matter or white matter Method, then it will necessarily cause error, it is therefore desirable to which partial volume correction is carried out to these regions(Partial Volume Correction, PVE), so as to obtain accurate CBF values.
The content of the invention
The purpose of the present invention is exactly to carry out partial volume correction using the CBF values obtained in ASL technologies, for examining for clear and definite AD It is disconnected that accurate algorithm is provided.Using the CBF values obtained in ASL and T1 values as input, show that full brain brain is white by PVE algorithms The CBF values of matter, ectocinerea and cerebrospinal fluid.Because CBF value pixels are low, on the position of full ectocinerea white matter cerebrospinal fluid boundary, one White matter of brain, the Multiple components such as ectocinerea and cerebrospinal fluid may be contained in individual pixel value, and cerebral blood flow (CBF) is in each composition(In vain Matter, grey matter, cerebrospinal fluid)Value differ greatly, therefore carry out partial volume correction after, it is more accurate to the CBF values of frontier district.It is main Want formula as follows:
CBFpve= CBF/(PGM*ΔMGM+ PWM*ΔMWM+PCSF*ΔMCSF).
In above-mentioned formula, CBF be correction before cerebral blood perfusion value, PGM, PWMAnd PCSFRepresent ectocinerea, white matter of brain and brain The probability graph of spinal fluid.ΔMGM, ΔMWM and ΔMCSFThe arterial spin labeling for dividing table to represent ectocinerea, white matter and cerebrospinal fluid fills Note weighted signal, wherein Δ MCSFIt is worth for 0, because the free marking signal of the artery of cerebrospinal fluid is almost 0.The CBF of adultgrayValue is about For WMCBF3 times of value.Therefore above-mentioned equation can be reduced to:
CBFpve= CBF/ (PGM*3+ PWM).
PVE corrections are carried out to ASL data with the formula after this simplification.
The present invention is achieved through the following technical solutions:
A kind of cerebral blood flow (CBF) partial volume effect bearing calibration, it is characterised in that comprise the following steps:
Step 1)The brain blood flow value and ectocinerea number of full brain are obtained by the testing result of nuclear-magnetism arterial spin labeling sequence Value and white matter of brain numerical value;
Step 2)The ectocinerea value in region of the ectocinerea numerical value less than or equal to 0.25 is set to 0, obtained ectocinerea value is not 0 region is GM;The white matter of brain value in region of the white matter of brain numerical value less than or equal to 0.25 is set to 0, obtained white matter of brain value is not Region for 0 is WM;
Step 3)For each numerical point in the GM of region, it is corrected as follows
CBFgray=CBF/(GMgray+0.33*GMwhite
Wherein CBFgrayFor the correction result of the brain blood flow value in the ectocinerea of the point, CBF is the brain blood flow value of the point, GMgrayFor the ectocinerea numerical value of the point, GMwhiteFor the white matter of brain numerical value of the point;
For each numerical point in the WM of region, it is corrected as follows
CBFwhite=CBF/(WMwhite+3*WMgray
Wherein CBFwhiteFor the correction result of the brain blood flow value in the white matter of brain of the point, CBF is the brain blood flow value of the point, WMwhiteFor the white matter of brain numerical value of the point, WMgrayFor the ectocinerea numerical value of the point;
Step 4)According to CBFgrayNumerical value make ectocinerea brain blood flow spirogram, according to CBFwhiteNumerical value to make brain white Matter brain blood flow spirogram.
Compared with prior art, the present invention has advantages below:
The CBF values carried out according to this method after PVE corrections are divided into two groups of white matter of brain, ectocinerea pixel values.With full brain CBF value phases Than result of calculation is more accurate, is advantageous to the diagnosis and analysis of clinic.
Brief description of the drawings
Fig. 1 schemes for full brain CBF;
Fig. 2 is the grey matter CBF figures after correction.
Embodiment
Embodiment 1
This method can use matlab programming realizations, and matlab programs are as follows:
function FG_PVE_correction_for_perfusiondata(Imgs,Grays,Whites,smooth_ kernel)
A PVE function is defined, function variable has four,(Full brain, white matter, grey matter, smooth nuclear parameter)
sub=[1,2,3,4,5,6,7,….,n] (Picture numbers are read one by one, and n is the data number into program)
for k=1:x(X indicates entry into the data amount check of program)
n=sub(k) (N represents the numbering of data)
S=num2str(n);(Numerical variable is changed to by character variable)
Imgs=strcat (' local path ', ' means', S, ' .img') obtain full brain pixel value
Grays=strcat (' local path ' ', ' c1s', S, ' .img') obtain ectocinerea pixel value
Whites=strcat (' local path ', ' c2s', S, ' .img') obtain white matter of brain pixel value
smooth_kernel=[2 2 2];It is 2 to define smooth check figure
for i=1:Size (Grays, 1) reads brain data
Vmat=spm_vol(Imgs(i,:));Full brain pixel value is assigned to Vmat variables
V=spm_read_vols(Vmat);Vmat variables are read in V variables
GM=spm_read_vols(spm_vol(Grays (i,:)));Grey matter data are assigned to GM
WM=spm_read_vols(spm_vol(Whites(i,:)));White matter data are assigned to WM
GM(find(GM<=0.25))=0;Grey matter is set to 0 by grey matter numerical value less than or equal to 0.25
WM(find(WM<=0.25))=0;White matter is set to 0 by white matter numerical value less than or equal to 0.25
V21=zeros(size(V));Define V21 variable spaces size and be equal to V(Brain space size), assignment 0
V22=zeros(size(V));Define V22 variable spaces size and be equal to V(Brain space size), assignment 0
V21(find(GM))=V(find(GM))./(GM(find(GM))+0.4*WM(find(GM)));
V21 (grey matter)=grey matter/(Grey matter+0.33* white matters)
V22(find(WM))=V(find(WM))./(WM(find(WM))+2.5*GM(find(WM)));V22 (white matter)= White matter/(White matter+3* grey matters)
Vmat.fname=deblank(PVE_gray_imgs(i,:));
spm_write_vol(Vmat,V21);
V21 numerical value is assigned to PVE_gray_imgs
Vmat.fname=deblank(PVE_white_imgs(i,:));
spm_write_vol(Vmat,V22);
V22 numerical value is assigned to PVE_ white _ imgs
end
fprintf ('\n-----PVE correction is done!......\n')
Gui interface prompting PVE calibrations have been finished
end
Fig. 1 and Fig. 2 indicates the calibration result of this method, and it is white that the CBF values carried out according to this method after PVE corrections are divided into brain Two groups of matter, ectocinerea cerebral blood flow (CBF) images.Fig. 2 is the grey matter brain blood flow spirogram after correction, compared with full brain CBF values, calculates knot Fruit is more accurate, is advantageous to the diagnosis and analysis of clinic.

Claims (1)

1. a kind of cerebral blood flow (CBF) partial volume effect bearing calibration, it is characterised in that comprise the following steps:
Step 1)The brain blood flow value and ectocinerea number of full brain are obtained by the testing result of nuclear-magnetism arterial spin labeling sequence Value and white matter of brain numerical value;
Step 2)The ectocinerea value in region of the ectocinerea numerical value less than or equal to 0.25 is set to 0, obtained ectocinerea value is not 0 region is GM;The white matter of brain value in region of the white matter of brain numerical value less than or equal to 0.25 is set to 0, obtained white matter of brain value is not Region for 0 is WM;
Step 3)For each numerical point in the GM of region, it is corrected as follows
CBFgray=CBF/(GMgray+0.33*GMwhite
Wherein CBFgrayFor the correction result of the brain blood flow value in the ectocinerea of the point, CBF is the brain blood flow value of the point, GMgrayFor the ectocinerea numerical value of the point, GMwhiteFor the white matter of brain numerical value of the point;
For each numerical point in the WM of region, it is corrected as follows
CBFwhite=CBF/(WMwhite+3*WMgray
Wherein CBFwhiteFor the correction result of the brain blood flow value in the white matter of brain of the point, CBF is the brain blood flow value of the point, WMwhiteFor the white matter of brain numerical value of the point, WMgrayFor the ectocinerea numerical value of the point;
Step 4)According to CBFgrayNumerical value make ectocinerea brain blood flow spirogram, according to CBFwhiteNumerical value to make brain white Matter brain blood flow spirogram.
CN201710639032.XA 2017-07-31 2017-07-31 Partial volume effect correction method for cerebral blood flow Expired - Fee Related CN107456225B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710639032.XA CN107456225B (en) 2017-07-31 2017-07-31 Partial volume effect correction method for cerebral blood flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710639032.XA CN107456225B (en) 2017-07-31 2017-07-31 Partial volume effect correction method for cerebral blood flow

Publications (2)

Publication Number Publication Date
CN107456225A true CN107456225A (en) 2017-12-12
CN107456225B CN107456225B (en) 2020-06-09

Family

ID=60547074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710639032.XA Expired - Fee Related CN107456225B (en) 2017-07-31 2017-07-31 Partial volume effect correction method for cerebral blood flow

Country Status (1)

Country Link
CN (1) CN107456225B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116029983A (en) * 2022-12-05 2023-04-28 中国科学院心理研究所 Image quality evaluation method and system for magnetic resonance arterial spin labeling perfusion imaging

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1062649A (en) * 1991-12-04 1992-07-15 中国人民解放军海军医学研究所 Brain figure imaging system and cerebral blood flow landform drawing generating method
US20100014735A1 (en) * 2008-07-21 2010-01-21 Siemens Medical Solutions Usa, Inc. Enhanced Contrast MR System Accommodating Vessel Dynamic Fluid Flow
WO2015161156A1 (en) * 2014-04-18 2015-10-22 The General Hospital Corporaton System and method for correcting intrinsic heterogeneity in magnetic resonance imaging
CN105654485A (en) * 2015-12-30 2016-06-08 中国人民解放军第四军医大学 Magnetic resonance arterial spin labeling sequence partial volume correction method using spatio-temporal information
US20160239966A1 (en) * 2008-05-28 2016-08-18 The Trustees Of Columbia University In The City Of New York Voxel-based methods for assessing subjects using molecular magnetic resonance imaging
CN106127792A (en) * 2016-07-22 2016-11-16 杭州师范大学 Magnetic resonance arterial spin labeling brain perfusion imaging data artefact figure minimizing technology

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1062649A (en) * 1991-12-04 1992-07-15 中国人民解放军海军医学研究所 Brain figure imaging system and cerebral blood flow landform drawing generating method
US20160239966A1 (en) * 2008-05-28 2016-08-18 The Trustees Of Columbia University In The City Of New York Voxel-based methods for assessing subjects using molecular magnetic resonance imaging
US20100014735A1 (en) * 2008-07-21 2010-01-21 Siemens Medical Solutions Usa, Inc. Enhanced Contrast MR System Accommodating Vessel Dynamic Fluid Flow
WO2015161156A1 (en) * 2014-04-18 2015-10-22 The General Hospital Corporaton System and method for correcting intrinsic heterogeneity in magnetic resonance imaging
CN105654485A (en) * 2015-12-30 2016-06-08 中国人民解放军第四军医大学 Magnetic resonance arterial spin labeling sequence partial volume correction method using spatio-temporal information
CN106127792A (en) * 2016-07-22 2016-11-16 杭州师范大学 Magnetic resonance arterial spin labeling brain perfusion imaging data artefact figure minimizing technology

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HENRI J.M.M. MUTSAERTS ET AL: "Gray matter contamination in arterial spin labeling white matter perfusion measurements in patients with dementia", 《NEUROIMAGE: CLINICAL》 *
JESSICA E SIMON ET AL: "A novel method to derive separate gray and white matter cerebral blood flow measures from MR imaging of acute ischemic stroke patients", 《JOURNAL OF CEREBRAL BLOOD FLOW & METABOLISM》 *
WEN-CHAU WU ET AL: "Measurement of Cerebral White Matter Perfusion Using Pseudocontinuous Arterial Spin Labeling 3T Magnetic Resonance Imaging – an Experimental and Theoretical Investigation of Feasibility", 《PLOS ONE》 *
XIAOYUN LIANG ET AL: "Improved Partial Volume Correction for Single Inversion Time Arterial Spin Labeling Data", 《MAGNETIC RESONANCE IN MEDICINE》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116029983A (en) * 2022-12-05 2023-04-28 中国科学院心理研究所 Image quality evaluation method and system for magnetic resonance arterial spin labeling perfusion imaging

Also Published As

Publication number Publication date
CN107456225B (en) 2020-06-09

Similar Documents

Publication Publication Date Title
De Guio et al. Reproducibility and variability of quantitative magnetic resonance imaging markers in cerebral small vessel disease
Park et al. Reliable and robust method for abdominal muscle mass quantification using CT/MRI: an explorative study in healthy subjects
US20240055126A1 (en) Medical Image Segmentation Using Deep Learning Models Trained with Random Dropout and/or Standardized Inputs
Jahanshad et al. Disrupted brain networks in the aging HIV+ population
Marshall et al. Altered cerebral hemodyamics and cortical thinning in asymptomatic carotid artery stenosis
Gold et al. Water exchange rate across the blood‐brain barrier is associated with CSF amyloid‐β 42 in healthy older adults
Shi et al. Diffusion tensor imaging–based characterization of brain neurodevelopment in primates
US20130289387A1 (en) Image processing apparatus and magnetic resonance imaging apparatus
Rapacchi et al. Reducing view‐sharing using compressed sensing in time‐resolved contrast‐enhanced magnetic resonance angiography
JP2015208679A (en) Magnetic resonance imaging apparatus, method, and program
Pai et al. BRAHMA: Population specific T1, T2, and FLAIR weighted brain templates and their impact in structural and functional imaging studies
US9971009B2 (en) Magnetic resonance imaging (MRI) with auto-detection and adaptive encodings for offset frequency scanning
CN107456225A (en) A kind of cerebral blood flow (CBF) partial volume effect bearing calibration
Haggenmüller et al. Pictorial review on imaging findings in cerebral CTP in patients with acute stroke and its mimics: a primer for general radiologists
US20240053421A1 (en) Systems and methods for magnetic resonance imaging
Ayaz et al. Brain MR image simulation for deep learning based medical image analysis networks
CN106137198A (en) A kind of MR imaging method and device
CN109472803B (en) Intracranial artery blood vessel segmentation method and system
Madsen et al. Mapping dynamic changes in ventricular volume onto baseline cortical surfaces in normal aging, mci, and Alzheimer’s disease
Zeinali et al. Brain volume estimation enhancement by morphological image processing tools
Suh et al. Development and validation of a deep learning-based automatic segmentation model for assessing intracranial volume: comparison with NeuroQuant, FreeSurfer, and SynthSeg
Ngamsombat et al. Evaluation of ultrafast Wave-CAIPI 3D FLAIR in the visualization and volumetric estimation of cerebral white matter lesions
Dietrich et al. Integrated intravoxel incoherent motion tensor and diffusion tensor brain MRI in a single fast acquisition
CN110739078B (en) Aneurysm rupture risk assessment method and system
KR20150073519A (en) Method of classifying cerebral infarction region

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200609

CF01 Termination of patent right due to non-payment of annual fee