Brain mapping and detection of functional patterns in fMRI using wavelet transform; application in detection of dyslexia
© Ji et al; licensee BioMed Central Ltd. 2009
Published: 3 November 2009
Functional Magnetic Resonance Imaging (fMRI) has been proven to be useful for studying brain functions. However, due to the existence of noise and distortion, mapping between the fMRI signal and the actual neural activity is difficult. Because of the difficulty, differential pattern analysis of fMRI brain images for healthy and diseased cases is regarded as an important research topic. From fMRI scans, increased blood ows can be identified as activated brain regions. Also, based on the multi-sliced images of the volume data, fMRI provides the functional information for detecting and analyzing different parts of the brain.
In this paper, the capability of a hierarchical method that performed an optimization algorithm based on modified maximum model (MCM) in our previous study is evaluated. The optimization algorithm is designed by adopting modified maximum correlation model (MCM) to detect active regions that contain significant responses. Specifically, in the study, the optimization algorithm is examined based on two groups of datasets, dyslexia and healthy subjects to verify the ability of the algorithm that enhances the quality of signal activities in the interested regions of the brain. After verifying the algorithm, discrete wavelet transform (DWT) is applied to identify the difference between healthy and dyslexia subjects.
We successfully showed that our optimization algorithm improves the fMRI signal activity for both healthy and dyslexia subjects. In addition, we found that DWT based features can identify the difference between healthy and dyslexia subjects.
The results of this study provide insights of associations of functional abnormalities in dyslexic subjects that may be helpful for neurobiological identification from healthy subject.
Experimental approaches towards the study of the brain functions often involve the detection of brain regions that respond differently to external stimuli. Thus, several researchers [1, 2] studied the brain activity when specific external stimuli is incited. An effective technology to study the human brain is functional magnetic resonance imaging (fMRI). It measures blood oxygenation-level-dependent (BOLD) signal changes caused by hemodynamic response of neural activity [3, 4]. However, due to the presence of different sources of noise and distortion, the mapping between the fMRI signal and the actual neural activity requires significant filtering and signal processing . For instance, unwanted changes in fMRI signal intensity can be caused by head motions as a main source of noise and distortion; this is commonly referred to motion artifact . Because of these difficulties, differential analysis of healthy and diseased cases using signal processing and pattern analysis of fMRI brain images remains as one of the main research challenges in the field .
Functional MRI is a well known technique for studying physiological conditions of the brains in human subjects. It is a non-invasive technology that tracks the changes in BOLD responses related with neural activities. From fMRI scans, the increased blood flow to the activated areas of brain is detected . Based on a series of multi-slice images of the brain, fMRI provides the functional information of detecting and analyzing different parts of the brain [9, 10]. Through the analysis of dynamic views of brain activities, the functions and high-level cognitive tasks can be studied. For instance, Ohnishi  shows that the classification of morphologic changes in the brain with normal aging and Alzheimer disease is different . Gallaghera  provides functional imaging studies on 'theory of mind' in verbal and non-verbal to identify distinct active region. Since fMRI signals are generated by the changes of BOLD in neural activity, the functional role of neural activity might be addressed by detecting the activated areas .
Dyslexia is a significant disease of reading acquisition . Habib  indicates that in spite of appropriate educational opportunities, five to ten percent of school students do not have reading ability simply because they suffer from dyslexia. Dyslexia, which is defined as a specific disability of reading and learning, was first identified about 100 years ago as a learning imparity, while at the early stage dyslexia is often taken as a vision problem [17, 18]. Nowadays, fMRI can determine that when people with dyslexic patients read, they do use different parts of their brain compared to people without dyslexia . Many previous studies [20–22] discuss investigations of the specific brain activities regions involved in dyslexia. According to the studies [1, 23], posterior part of the left superior temporal gyrus and the inferior parietal gyrus (BA (Broadmann area) 40) are the most likely to involve language processes. However, the most significant difficulty in dyslexia is the lack of identifying the patterns of dyslexia.
Therefore, a possible approach to identify dyslexia and healthy subjects is to adopt a signal processing technique to extract valuable information from the fMRI signal of interested regions. In our previous study , we proposed a signal processing technique using an optimization theory, called modified MCM, to filter, process, and classify the fMRI images. To analyze the fMRI time-series, pixel-regions in fMRI are used to explore the brain activities that significantly improves conventionally maximum correlation modeling (MCM) of the brain hemodynamic response across the healthy and dyslexia subjects. The intention of this paper thus as continuation of our previous work is that presenting the improvements of filtering technique developed in our previous study to compare the brain activities of the healthy and dyslexic subjects. In addition, the utility of discrete wavelet transform (DWT) to identify different patterns between healthy and dyslexia subjects is presented. Simple statistics based on mean and standard deviation quantitative analysis is used to provide knowledge of the signals' significance. However, standard deviation and mean may not provide an appropriate characterization of the rapid changes in a signal. For this reason, further measurement, Higuchi fractal dimension (FD) and discrete wavelet transform (DWT) features are used to describe the difference between before and after signal filtering applied. Before filtering indicates a raw time series of single activated pixel at regions of interest based on statistical parametric mapping (SPM) and filtering signal suggests that the results of linearly combined signal with neighbors' pixel of the activated single pixel using a practical optimization algorithm. FD analysis is a useful tool in the identification of complexity under different conditions. In other words, FD is useful for measuring self-similarity of the signals. For this study, Higuchi FD is used because it is easy and simple to implement . Wavelet transform analysis [26, 27] is a very promising signal processing technique to detect abnormal changes within the signals. In particular, DWT decomposes the signals at different scales and resolutions. Thus, it is suitable for analyzing non-stationary signals.
Our hypotheses for this study are as follows:
1. We hypothesize that signal filtering using a hierarchical optimization algorithm with MCM developed in our previous study helps to enhance the signal activity of interested regions.
2. We hypothesize that the features extracted from the signals using wavelet transformation may help to differentiate between healthy and dyslexia subjects.
Description of dataset
Each functional scan consists of 96 sequential images as an echo-planar sequence. Parameters for these scans are: TR = 2500 ms, TE = 50 ms, and flip angle is set to 90. Each image consists of 28 transaxial slices that are 5 mm thick (no space between slices) with a field of view of 24 cm (frequency) × 15 cm (phase) and an acquisition matrix of 64 × 40 resulting in an in-plane resolution of 3.75 mm. High-resolution anatomic images are collected using an inversion recovery 3D spoiled gradient echo sequence with the following parameters: matrix is 256 × 256, field of view is 24 cm, section thicknesses 3 mm with no gap between slices, number of sections is 60, and in-plane resolution is 0.94 mm. This sequence is used both in anatomic overlays of the functional data and in spatial normalization of the data sets with a standard atlas.
The detail description of the practical optimization is introduced in our previous study . In briefly describing the method, the practical optimization algorithm contains two part i.e. SPM imaging analysis and modified MCM. The identification of specifically informative regions is considered in terms of dyslexic and healthy subjects with external stimuli (i.e. showing words (Task) and not showing words (No task)). Figure 2 depicted this boxcar diagram processes. In order to detect the highly informative regions, SPM is used [6, 28, 29]. For the interesting region selection, the statistical threshold is set to p <0.05, uncorrelated to find the most activated regions. For this study, six interested pixels which is activated when extrenal stimuli is applied between two groups are selected and used for further analysis.
Then, a significantly improved and modified version of Maximum Correlation Modeling (MCM) method presented by Friman  to detect interesting brain region's activities is described. MCM method is based on time series modeling of pixel neighborhoods (3 × 3) for detecting active parts that contain a significant hemodynamic response. To use the method, first we consider the regions in fMRI images and the time series of the pixels in the region. In order to estimate the activity of center pixel in each region, the activities of the symmetric neighbors of the center point are linearly combined using a vector weight which is determined by maximum correlation between the signals (time series) and known hemodynamic wave pattern. In order to maximize the correlation coefficient between the parameters of the time series and the hemodynamic response, an optimization process was conducted. This optimization process finds the best set of parameters under given constraints that identifies the best match between the expected hemodynamic response and the combined time-series. A more detail description of this algorithm is explained in our previous paper .
Statistical evaluation for signal enhancement
We evaluate the performance of filtering of our proposed method introduced in previous study. In our previous study, we showed a different pattern of pixels activity and its correlation between healthy and dyslexic subjects. Now, the performance of signal improvement with the proposed method is evaluated. Before filtering (time series is presented using only a single pixel which is activated) and after filtering (time series is presented with linear combination of neighbor's pixels) signals are therefore compared with statistical analysis software (SAS).
where represents the normalization factor for the curve length and N is the total length of the signal. < L(k) > is defined as the length of the curve for the time series k and < L m (k) denotes the average value over k. Thus, if < L(k) > ∞k-D, then the curve has dimension D. In other words, FD identifies the slope of the best fit-line at the log-log plot for log(L(k) <versus log(k) .
Feature extraction with discrete wavelet transformation
In this section, first the results between optimal filtering with neighbors' pixels and raw time-series activated single pixel on fMRI data of five healthy and five dyslexic subjects are presented as described in Method section. Then, the results of DWT to differentiate between healthy and dyslexia with ANOVA are presented.
The comparison between optimal filtering and before filtering using FD and wavelet based features among within subjects (d i indicates the detail coefficient at level i).
Within Healthy Subjects
Within Dyslexia Subjects
Sum of d1
Variance of d1
Mean absolute deviation of d1
Sum of d2
Variance of d2
Mean absolute deviation of d2
Sum of d3
Mean absolute deviation of d3
Higuchi FD (win = 8)
Higuchi FD (win = 15)
In our previous study, we presented the brain activation using a hierarchical optimization algorithm. In particular, an optimization algorithm is applied to specific regions of the brain identified by SPM to find patterns that are significantly different across two groups of subjects.
In our study, fMRI analysis is performed based on the activated areas for the tasks (showing the words/not showing the words). It must be addressed that we do not emphasizes the issue of the statistical significance of activations that are detected in this paper. Our approach allows us to compare the improvements of filtering and provides the capability of identification between healthy and dyslexic subjects with wavelet transform method.
The main aim of this paper is to not only compare the performance between before and after filtering using practical optimization algorithm with fractal dimension and DWT features, but also evaluate the capability of wavelet based features to distinguish between two groups. DWT and FD can show the improvements of signal filtering by comparing with single pixel time-series of brain activation. However, fractal dimension does not show that they can distinguish the different conditions. This finding demonstrates that the practical optimization algorithm can help to improve the filtering to obtain the brain activation signals. In addition, we found that wavelet based features may be useful to distinguish the healthy and dyslexic subjects. This fact supports the opportunity of using wavelet analysis to distinct different condition. Thus, it may possible to use early decision strategy to identify the abnormal activation of the brain in fMRI data. The key point of applying the wavelet transform, which is a multi-resolution method is selecting of the appropriate wavelet bases functions. Since there is no ultimate way of choosing the wavelet, one of Daubechies families, db4, is used for this study. But, the question still remains, what is the best choice of wavelet to extract the valuable information from the signal. Also, the detail coefficients is contained data variation information at a particular subband. Thus, they can present a knowledge which is contained within the dataset. For this study, detail coefficients at level 2 and 3 are only significant to discover two subjects. The detail coefficient at level 1 which is associated with the most high frequency of data may not provide big difference between two subjects due to the pattern similarity.
A limitation of this study is the small size of the dataset used to identify the fMRI brain activity. However, despite the small size of the dataset, the method based on DWT features successfully differentiates the two groups, healthy and dyslexia subjects. Therefore, additional analysis with a large amount of dataset with many pixels should be performed to test between dyslexic and healthy cases.
Conclusion & future work
As a continuation of our previous work, the purpose of this work is not only to evaluate the improvement of filtering with the hierarchical optimization algorithm, but also apply a signal processing strategy to identify two subjects, healthy and dyslexia subjects.
In this paper, we present the utility of practical optimization algorithm based on FD and DWT features. In addition, we describe that wavelet based features may useful to differentiate two groups, healthy and dyslexic subjects, of fMRI time series. In this study, we perform wavelet transform as an appropriate non-stationary signal analysis method, which may be suitable for differentiating two different conditions. For future work, more datasets with many regions will be tested and used to evaluate the two groups, non-dyslexic and dyslexic subjects. Also, other possible wavelet-based features are going to be identified for differentiating them.
The authors wish to thank the Radiology department of Wake Forest University for providing dataset. The authors also thank to Rebecca Smith for valuable feedbacks and discussions for this work.
This article has been published as part of BMC Medical Informatics and Decision Making Volume 9, Supplement 1, 2009: 2008 International Workshop on Biomedical and Health Informatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1472-6947/9?issue=S1.
- Cohen L, Dehaene S: Specialization within the ventral stream: the case for the visual word form area. Neuroimage. 2004, 22 (1): 466-476. 10.1016/j.neuroimage.2003.12.049.View ArticlePubMedGoogle Scholar
- Cohen JD, Forman SD, Braver TS, Casey BJ, Schreiber DS, Noll DC: Activation of the Prefrontal Cortex in a Nonspatial Working Memory Task With Functional MRI. Human Brain Mapping. 1994, 1: 293-304. 10.1002/hbm.460010407.View ArticlePubMedGoogle Scholar
- DeYoe E, Bandettini P, Neitz J, Miller D, Winans P: Functional magnetic resonance imaging (fmri) of the human brain. Journal of Neurosci Methods. 1994, 54 (2): 171-187. 10.1016/0165-0270(94)90191-0.View ArticleGoogle Scholar
- Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A: Neurophysiological investigation of the basis of the fmri signal. Nature. 2001, 412: 150-157. 10.1038/35084005.View ArticlePubMedGoogle Scholar
- Friston K, William S, Howard R, Frackowiak R, Turner R: Movement related effects in fmri time series. Magn Reson Med. 1996, 35: 346-355.View ArticlePubMedGoogle Scholar
- Friston K, Ashburner J, Kiebel S, Nichols T, Penny W: Statistical Parametric Mapping: The Analysis of Functional Brain Images. 2006, Academic Press, Erewhon, NCGoogle Scholar
- Ford J, Makedon F, Owen C, Johnson S, Steinberg T, Savkin AJ: Stimuius trackingin functional magneticresonance imaging (fmrl). ACM international conference on Multimedia. 1988, 445-454.Google Scholar
- Ogawa S, Lee T, Kay AR, Tank DW: Brain magnetic resonance imaging with contrast dependent on blood oxygenation. National Academic Science. 1990, 87: 9868-9872. 10.1073/pnas.87.24.9868.View ArticleGoogle Scholar
- Bandettini PA, Jesmanowicz A, Wong EC, Hyde JS: Processing strategies for time course data set in functional mri of the human brain. Magnetic Resonance in Medicine. 1993, 30: 161-173. 10.1002/mrm.1910300204.View ArticlePubMedGoogle Scholar
- Friman O, Borga M, Lundberg P, Knutsson H: Detection of neural activity in fmri using maximum correlation modeling. NeuroImage. 2002, 15 (15): 386-395. 10.1006/nimg.2001.0972.View ArticlePubMedGoogle Scholar
- Ohnishi T, Matsuda H, Tabira T, Asada T, Uno M: Changes in brain morphology in alzheimer disease and normal aging: Is alzheimer disease an exaggerated aging process?. AJNR Am J Neuroradiol. 2001, 22 (9): 1680-1685.PubMedGoogle Scholar
- Bookheimer SY, Strojwas MH, Cohen MS, Saunders AM, Pericak-Vance MA, Mazziotta JC, Small GW: Patterns of brain activation in people at risk for alzheimer's disease. The New England Journal of Medicine. 2000, 343: 450-456. 10.1056/NEJM200008173430701.PubMed CentralView ArticlePubMedGoogle Scholar
- Gallaghera H, Happe F, Brunswicka N, Fletchera P, Frithc U, Fritha C: Reading the mind in cartoons and stories: an fmri study of'theory of mind' in verbal and nonverbal tasks. Neuropsychologia. 2000, 38: 11-21. 10.1016/S0028-3932(99)00053-6.View ArticleGoogle Scholar
- Logothetis NK: What we can do and what we cannot do with fMRI. Nature. 2008, 453: 869-878. 10.1038/nature06976.View ArticlePubMedGoogle Scholar
- Ramus F: Dyslexia. Talk of two theories. Nature. 2001, 412: 393-395. 10.1038/35086683.View ArticlePubMedGoogle Scholar
- Habib M: The neurological basis of developmental dyslexia brain. Brain. 2000, 123 (12): 2373-2399. 10.1093/brain/123.12.2373.View ArticlePubMedGoogle Scholar
- Critchley M: The dyslexic child/Macdonald Critchley. 1970, Heinemann Medical, 2Google Scholar
- Livingstone MS, Rosen GD, Drislane FW, Galaburda AM: Physiological and anatomical evidence for a magnocellular defect in developmental dyslexia. Proc Natl Acad Sci U S A. 1991, 88 (18): 7943-7947. 10.1073/pnas.88.18.7943.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruff S, Cardebat D, Marie N, Demonet JF: Enhanced response of the left frontal cortex to slowed down speech in dyslexia: an fmri study. Neuroreport. 2002, 13 (10): 1285-1289. 10.1097/00001756-200207190-00014.View ArticlePubMedGoogle Scholar
- Schulz E, Maurer U, Sanne van der Mark , Bucher K, Brem S, Martin E, Brandeisa D: Impaired semantic processing during sentence reading in children with dyslexia: Combined fMRI and ERP evidence. neuroImage. 2008, 41 (1): 153-168. 10.1016/j.neuroimage.2008.02.012.View ArticlePubMedGoogle Scholar
- Siok WT, Niu Z, Jin Z, Perfetti CA, Tan LH: A structural-functional basis for dyslexia in the cortex of Chinese readers. Neuroreport. 2000, 11 (11): 2443-2447.View ArticleGoogle Scholar
- Kronbichler M, Wimmer H, Staffen W, Hutzler F, Mair A, Ladurner Gr: Developmental Dyslexia: Gray Matter Abnormalities in the Occipitotemporal Cortex. Human Brain mapping. 2008, 29 (5): 613-625. 10.1002/hbm.20425.View ArticlePubMedGoogle Scholar
- Simos PG, Breier JI, Wheless JW, Foorman BR, Castillo EM, Papanicolaou AC: Brain Mechanisms for Reading Words and Pseudowords: an Integrated Approach. Cerebral Cortex. 2002, 12 (3): 297-305. 10.1093/cercor/12.3.297.View ArticlePubMedGoogle Scholar
- Ji SY, Kayvan N: A Modified Maximum Correlation Modeling Method for fMRI Brain. IEEE Bioinformatics and Biomedicine (BIBM). 2008, 64-69.Google Scholar
- Higuchi T: Approach to an irregular time series on the basis of the fractal theory. Physica. 1998, 31: 277-283.Google Scholar
- Rao RM, Bopardikar AS: Wavelet Transforms Introduction to theory and applications?. 1998, Addison WesleyGoogle Scholar
- Burrus CS, Gopinath RA, Guo H: Introduction to Wavelets and Wavelet Transforms: A Primer. 1997, Prentice Hall IncGoogle Scholar
- Veltman D, Hutton C: SPM2 user guide for fMRI. 2001, [http://www.fil.ion.ucl.ac.uk/spm/doc/]Google Scholar
- Wimmer K: fMRI time series analysis with the software spm99. 2003, Master's thesis, University of ViennaGoogle Scholar
- Mandelbrot BB: The Fractal Geometry of Nature. 1983, W. H. FreemanGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.