CN100423691C - Method for analyzing functional MRI data by integration of time domain and space domain information - Google Patents
Method for analyzing functional MRI data by integration of time domain and space domain information Download PDFInfo
- Publication number
- CN100423691C CN100423691C CNB2005100869592A CN200510086959A CN100423691C CN 100423691 C CN100423691 C CN 100423691C CN B2005100869592 A CNB2005100869592 A CN B2005100869592A CN 200510086959 A CN200510086959 A CN 200510086959A CN 100423691 C CN100423691 C CN 100423691C
- Authority
- CN
- China
- Prior art keywords
- voxel
- data
- resting state
- functional mri
- time
- 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.)
- Expired - Fee Related
Links
- 238000002599 functional magnetic resonance imaging Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000010354 integration Effects 0.000 title 1
- 230000000284 resting effect Effects 0.000 claims description 33
- 230000000694 effects Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 230000003925 brain function Effects 0.000 claims 4
- 230000003340 mental effect Effects 0.000 abstract 1
- 210000004556 brain Anatomy 0.000 description 21
- 238000011160 research Methods 0.000 description 9
- 230000002146 bilateral effect Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000007619 statistical method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000006213 oxygenation reaction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 210000001103 thalamus Anatomy 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000133 brain stem Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 210000001652 frontal lobe Anatomy 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 210000001259 mesencephalon Anatomy 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001936 parietal effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000002360 prefrontal effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 230000000392 somatic effect Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a functional magnetic resonance imaging technique, which is characterized by the following: predisposing data; metering single local consistence through Kendall coefficient; obtaining fluctuation strength in the effective frequent band through standard variance of time sequence; multiplying Kendall coefficient by standard variance as mental level factor to display the active level.
Description
Technical field
The present invention relates to the functional mri technical field, a kind of functional MRI data time domain and spatial-domain information are combined detected the method for cerebration.
Background technology
Based on (the Blood Oxygenation Level Dependent of oxygenation level in the blood, BOLD) functional MRI (Functional Magnetic ResonanceImaging, FMRI) technology has been widely used among the research of live body brain, such as the neural activity mechanism of research live body brain in perception, cognition and affective activity process.In traditional functional mri research, obtain brain action message with respect to the baseline state under this task in the similarities and differences of BOLD signal under the particular task and BOLD response model by comparing the live body brain.
Even people such as nineteen ninety-five Biswal discover under resting state, low-frequency oscillation (the Low FrequencyFluctuation of BOLD signal in the human brain, LFFs 0.01-0.08Hz) has concordance (Biswal et al., 1995) highly between the different motion cortex of human brain.People such as Biswal think that the low-frequency component of BOLD signal has reflected the self-organization of brain.Cerebration research under the resting state has very important meaning: though the reaction under various particular tasks has had more understanding for brain at present, also know few for people's cerebration of resting state; And in the data acquisition of clinical research, need not patient and cooperate the complicated task of finishing, only under resting state, count the minute data collection, promptly may find the situation of change of the relative normal control of the specific brain zone function of patient.
Owing to there is not particular task to stimulate, the analysis of resting state BOLD signal is lacked necessary reference model to obtain cerebration information.Situation analysis provides good enlightenment to the resting state cerebration but task status activates two features in brain district: one, on spatial domain, each voxel of movable brain district show high locally coherence (Regional Homogeneity, ReHo).The brain that people such as Zang Yufeng once were applied to this characteristic under the motion task activates detection (Zang et al., 2004).Its two, on time domain, each voxel signal of movable brain district shows strong undulatory property in effective band.Such as one 60 seconds to be under the tile designs task in cycle, activating the voxel time series should have very strong fluctuation on this frequency content of 1/60Hz.In other words, we can extract each voxel BOLD signal 1/60Hz composition, select those voxels that big variance is arranged as possible activation voxel under this radio-frequency component.Based on this characteristic, people such as Fransson find out under this state remarkable activity (Fransson, 2005) such as brain districts such as posterior cingutate, veutro prefrontal lobes according to resting state BOLD signal low-frequency component.
Any remarkable active voxel need possess above two characteristics simultaneously.60 seconds being that the tile designs task in cycle is an example, the brain district of apparent altitude locally coherence only under this task, this concordance may be owing to causing with the irrelevant factor of task.Equally, do not possess the height locally coherence and only be the BOLD signal at the voxel of the strong undulatory property of 1/60Hz performance, its strong undulatory property may since random noise cause.To sum up, take all factors into consideration the locally coherence information of resting state BOLD signal space territory performance and the fluctuation information in the time domain effective band (0.01-0.08Hz), can effectively detect the cerebration situation under this state.Up to the present, we do not retrieve the locally coherence information of comprehensive voxel and the report that the interior fluctuation information of BOLD signal effective band is analyzed the cerebration level as yet.
List of references:
Biswal?BB，Yetkin?FZ，Haughton?VM，Hyde?JS.Functional?connectivity?in?themotor?cortex?of?resting?human?brain?using?echo-planar?MRI.Magn?ResonMed?1995；34：537-541.
Zang?YF，Jiang?TZ，Lv?YL，He?Y，Tian?LX.Regional?homogeneity?approach?tofMRI?data?analysis.Neuroimage?2004；22：394-400.
Fransson?P.Spontaneous?low-frequency?BOLD?signal?fluctuations-an?fMRIinvestigation?of?the?resting-state?default?mode?of?brain?function?hypothesis.Hum?Brain?Mapp?2005；26：15-29.
Summary of the invention
Core of the present invention is, utilize computer equipment, the method that adopts time domain information (effectively in-band signal cymomotive force) to combine with spatial-domain information (locally coherence) is carried out quantitative analysis to each voxel active level of functional MRI data, to reach the purpose that detects the cerebration level.
Described to based on resting state hypencephalon functional MRI data, the locally coherence information of each voxel is combined with the cymomotive force of signal low-frequency component (0.01-0.08Hz), to detect the cerebration level under this state.
After the functional MRI original time series carried out pretreatment, for each voxel, measure the locally coherence of its performance with Kendall, measure the undulatory property of its low frequency signal (0.01-0.08Hz) with standard variance, with the two product as activity index (the Resting-State Activity Index of this voxel under resting state, RSAI), thus the cerebration under this state is carried out quantitative assessment.To full brain data one by one voxel calculate its RSAI, can obtain full cerebration level.This method realizes simple, calculates fast, and hardware requirement is low, can be widely used in the functional MRI basic research and clinical disease research of resting state.
Description of drawings
Fig. 1 is the method flow diagram of comprehensive time domain of the present invention and spatial information (si) analytic function MR data.
Fig. 2 is an adolescents with normal RSAI cartogram.
Fig. 3 is patient ADHD and normal control RSAI statistical discrepancy figure.
The specific embodiment
Its implementation procedure can be divided into following 4 steps, as shown in Figure 1:
Step 1, resting state functional MRI data are gathered.Being captured on the magnetic resonance scanner that possesses plane echo-wave imaging (EPI) sequence of functional MRI data finished.The concrete parameter of imaging does not have specific (special) requirements, and the sampling time is preferably in 6-8 minute, and the repetition time, (Repetition Time TR) generally selected for use below 2 seconds or 2 seconds, and spatial resolution is generally below 5 millimeters.Only need tested peace and quiet close order in the data acquisition and have a rest, keep head still as far as possible, and try not to collect one's thoughts and ponder a problem;
Step 2, data pretreatment.Generally need carry out the pretreatment of following steps to the data that collect, comprise that acquisition time correction, a dynamic(al) correction of different time points image data, the space criteriaization between aspect arrives regulation mould plate, space resampling, space smoothing etc.After these basic processes are finished, time series is carried out bandpass filtering (0.01-0.08Hz), to obtain the signal low-frequency component.With the composition filtering below the 0.01Hz in the signal mainly is to remove some wait the extremely low frequency composition such as drift influence;
Step 3, voxel calculates the standard deviation (to measure the undulatory property of its time domain effective ingredient) of its Kendall (to measure the locally coherence of its spatial domain) and its BOLD signal low-frequency component one by one.Wherein, it is as follows to be used for measuring local conforming Kendall computing formula:
Here KCC represents Kendall, and n is the number of this voxel time series time point, and K is a voxel number in this voxel neighborhood, generally selects K=7, and 19,27, R
iBe the order sum of K voxel at time point i, R is R
iAverage at n time point;
Step 4, resting state cerebration index (Resting State Activity Index, RSAI) calculate, the standard deviation of the Kendall of voxel and its BOLD signal low-frequency component multiplies each other, product promptly is the activity index RSAI of this voxel at resting state, so far, the level of activation that obtains under this voxel resting state is measured.
Any that should be noted that is when the cerebration index RSAI to the burst data resting state carries out statistical analysis, because this index is disobeyed normal distribution or distributed such as T distribution, χ with closely-related other of normal state substep
2Distribute etc., suggestion is used robust comparatively and is widely used in distributing the nonparametric statistical method of unknown data.
Provide two experimental results that adopt the present invention to handle the resting state functional MRI data below.
Fig. 2 is the cerebration statistical picture that obtains according to 10 adolescents with normal resting states of methods analyst of the present invention functional MRI data.The statistical method that adopts is the Wilcoxon-Mann-Whitney rank test.As can be seen from the figure this group teenager has remarkable cerebration at bilateral posterior cingutate/cuneus, bilateral fusiform gyrus, inboard time of bilateral veutro anterior cingutate/frontal lobe, bilateral superior temporal gyrus, bilateral thalamus and left side inferior parietal lobule under resting state.
Fig. 3 is according to method of the present invention moving obstacle (ADHD) teenager and 10 statistical discrepancy images that adolescents with normal resting state cerebration obtains of 8 attention deficits/how relatively.The statistical method that adopts also is the Wilcoxon-Mann-Whitney rank test.As can be seen from the figure midbrain part, left side superior temporal gyrus (auxiliary auditory area) and the left side postcentral gyrus (basic somatic sensory area) at bilateral cuneus (basic visual area, accessorial visual district), bilateral thalamus, left side brain stem shows the cerebration stronger than normal control under the ADHD teenager resting state.
Patient ADHD and normal control RSAI statistical discrepancy figure.On behalf of patient ADHD, black region more normally contrast stronger leave cerebration.Compare with patient ADHD, normal control does not have the leave activity significantly to be better than patient's ADHD brain district.The threshold value of statistical test is p＜0.05, and minimum agglomerate volume is greater than 1188mm
3
Carry out functional MRI research under resting state following major advantage is arranged: (1) data acquisition is simple, and tested needs rest to get final product, and this is for some disease research significant (patient need not to cooperate and finishes complicated experimental duties); (2) multiformity of relative task status experimental design, the simplicity of resting state self can increase different experiments result's comparability; (3) compare with the particular task design, resting state more approaches the naturalness of human brain.Comprehensive time domain of the present invention and spatial information (si) provide a kind of new effective tool for analyzing the resting state functional MRI data.
Claims (4)
1. method that the brain function MR data is handled, it is characterized in that, utilize computer equipment, the method that adopts time domain information to combine with spatial-domain information is carried out quantitative analysis to each voxel active level of functional MRI data, reaching the purpose that detects the cerebration level, to based on resting state hypencephalon functional MRI data, the locally coherence information of each voxel is combined with the cymomotive force of signal low-frequency component, to detect the cerebration level under this state, its step is as follows:
Step 1, resting state functional MRI data are gathered;
Step 2, data pretreatment need be carried out the pretreatment of following steps to the data that collect, and comprise that acquisition time correction, a dynamic(al) correction of different time points image data, the space criteriaization between aspect resamples to regulation mould plate, space, space smoothing;
Step 3, voxel calculates the standard deviation of its Kendall and its BOLD signal low-frequency component one by one, and wherein, it is as follows to be used for measuring local conforming Kendall computing formula:
Here KCC represents Kendall, and n is the number of this voxel time series time point, and K is a voxel number in this voxel neighborhood, selects K=7, and 19,27, R
iBe the order sum of K voxel at time point i, R is R
iAverage at n time point;
Step 4, resting state cerebration Index for Calculation, the standard deviation of the Kendall of voxel and its BOLD signal low-frequency component multiplies each other, product promptly is the activity index RSAI of this voxel at resting state, so far, obtains the level of activation tolerance under this voxel resting state.
2. the method that the brain function MR data is handled according to claim 1 is characterized in that step 1, the collection of resting state functional MRI data; Sampling time, the repetition time selected for use below 2 seconds or 2 seconds at 6-8 minute, and spatial resolution is below 5 millimeters.
3. the method that the brain function MR data is handled according to claim 1 is characterized in that, step 2, data pretreatment are carried out bandpass filtering to time series, to obtain the signal low-frequency component.
4. the method that the brain function MR data is handled according to claim 3 is characterized in that, time series is carried out bandpass filtering, to obtain 0.01-0.08Hz low frequency composition signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005100869592A CN100423691C (en) | 2005-11-23 | 2005-11-23 | Method for analyzing functional MRI data by integration of time domain and space domain information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005100869592A CN100423691C (en) | 2005-11-23 | 2005-11-23 | Method for analyzing functional MRI data by integration of time domain and space domain information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1969746A CN1969746A (en) | 2007-05-30 |
| CN100423691C true CN100423691C (en) | 2008-10-08 |
Family
ID=38110913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2005100869592A Expired - Fee Related CN100423691C (en) | 2005-11-23 | 2005-11-23 | Method for analyzing functional MRI data by integration of time domain and space domain information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100423691C (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103345749B (en) * | 2013-06-27 | 2016-04-13 | 中国科学院自动化研究所 | A kind of brain network function connectivity lateralization detection method based on modality fusion |
| CN103886591A (en) * | 2014-03-06 | 2014-06-25 | 西安电子科技大学 | Brain nuclei Granger causal analysis method based on RYGB surgery weight losing |
| CN104434109B (en) * | 2014-12-19 | 2017-02-22 | 大连海事大学 | Functional nuclear magnetic resonance time sequence matching method |
| CN105708462A (en) * | 2016-01-14 | 2016-06-29 | 内蒙古医科大学附属医院 | Data processing method based on rfMRI (resting-state functional magnetic resonance imaging) of idiopathic epilepsy |
| CN112528743B (en) * | 2020-11-09 | 2024-05-10 | 江苏海洋大学 | Spectrum difference mapping frame method for decoding brain activities based on functional magnetic resonance technology |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003339670A (en) * | 2003-07-01 | 2003-12-02 | Toshiba Corp | Magnetic resonance imaging system |
| CN1628608A (en) * | 2003-12-15 | 2005-06-22 | 中国科学院自动化研究所 | Functional magnetic resonance data processing method utilizing partial uniformity method |
-
2005
- 2005-11-23 CN CNB2005100869592A patent/CN100423691C/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003339670A (en) * | 2003-07-01 | 2003-12-02 | Toshiba Corp | Magnetic resonance imaging system |
| CN1628608A (en) * | 2003-12-15 | 2005-06-22 | 中国科学院自动化研究所 | Functional magnetic resonance data processing method utilizing partial uniformity method |
Non-Patent Citations (2)
| Title |
|---|
| Functional Connectivity in the Motor Cortex of RestingHumanBrain Using Echo-Planar MRI. Bharat Biswal, F. Zerrin Yetkin, Victor M. Haughton, JamesS.Hyde.Magnetic Resonance in Medicine,Vol.vol.34 No.,no.4. 1995 * |
| Modulation of Functional Connectivity During theRestingStateand the Motor Task. Tianzi Jiang, Yong He, Yufeng Zang,and Xuchu Weng.Human Brain Mapping,Vol.22 . 2004 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1969746A (en) | 2007-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huotari et al. | Sampling rate effects on resting state fMRI metrics | |
| Filippi et al. | Correlations between structural CNS damage and functional MRI changes in primary progressive MS | |
| Steele et al. | Neuroimaging measures of error-processing: Extracting reliable signals from event-related potentials and functional magnetic resonance imaging | |
| Stevens et al. | Hemispheric differences in hemodynamics elicited by auditory oddball stimuli | |
| Proskovec et al. | Load modulates the alpha and beta oscillatory dynamics serving verbal working memory | |
| Trevarrow et al. | The developmental trajectory of sensorimotor cortical oscillations | |
| Logothetis | The neural basis of the blood–oxygen–level–dependent functional magnetic resonance imaging signal | |
| Basten et al. | Trait anxiety and the neural efficiency of manipulation in working memory | |
| Shaw et al. | Abnormal functional connectivity in posttraumatic stress disorder | |
| Spooner et al. | Aberrant oscillatory dynamics during somatosensory processing in HIV-infected adults | |
| Hale et al. | Comparison of functional connectivity in default mode and sensorimotor networks at 3 and 7T | |
| Liu et al. | Baseline BOLD correlation predicts individuals' stimulus-evoked BOLD responses | |
| DE102016202254B4 (en) | Model-free determination of image areas with anomalous diffusion based on diffusion-weighted magnetic resonance image data | |
| CN100423691C (en) | Method for analyzing functional MRI data by integration of time domain and space domain information | |
| Pope et al. | Force related activations in rhythmic sequence production | |
| Dorokhov et al. | Experimental model of study of consciousness at the awakening: fMRI, EEG and behavioral methods | |
| Shaw et al. | Evaluating subject specific preprocessing choices in multisubject fMRI data sets using data-driven performance metrics | |
| Aso et al. | Comparison of diffusion-weighted fMRI and BOLD fMRI responses in a verbal working memory task | |
| Spooner et al. | Methodological considerations for a better somatosensory gating paradigm: The impact of the inter-stimulus interval | |
| Kyathanahally et al. | A realistic framework for investigating decision making in the brain with high spatiotemporal resolution using simultaneous EEG/fMRI and joint ICA | |
| He et al. | Alteration of resting-state functional connectivity in the sensorimotor network in patients with thalamic infarction | |
| Machner et al. | Resting-state functional connectivity in the dorsal attention network relates to behavioral performance in spatial attention tasks and may show task-related adaptation | |
| Niu et al. | Comparing test-retest reliability of entropy methods: Complexity analysis of resting-state fMRI | |
| Liu et al. | An introduction to normalization and calibration methods in functional MRI | |
| Labounek et al. | Exploring task-related variability in fMRI data using fluctuations in power spectrum of simultaneously acquired EEG |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20081008 Termination date: 20171123 |