A Review on MR Image Intensity Inhomogeneity Correction

Deep learning based correction of RF field induced inhomogeneities for T2w prostate imaging at 7 T

At ultrahigh field strengths images of the body are hampered by B1 -field inhomogeneities. These present themselves as inhomogeneous signal intensity and contrast, which is regarded as a "bias field" to the ideal image. Current bias field correction methods, such as the N4 algorithm, assume a low frequency bias field, which is not sufficiently valid for T2w images at 7 T. In this work we propose a deep learning based bias field correction method to address this issue for T2w prostate images at 7 T. By combining simulated B1 -field distributions of a multi-transmit setup at 7 T with T2w prostate images at 1.5 T, we generated artificial 7 T images for which the homogeneous counterpart was available. Using these paired data, we trained a neural network to correct the bias field. We predicted either a homogeneous image (t-Image neural network) or the bias field (t-Biasf neural network). In addition, we experimented with the single-channel images of the receive array and the corresponding sum of magnitudes of this array as the input image. Testing was carried out on four datasets: the test split of the synthetic training dataset, volunteer and patient images at 7 T, and patient images at 3 T. For the test split, the performance was evaluated using the structural similarity index measure, Wasserstein distance, and root mean squared error. For all other test data, the features Homogeneity and Energy derived from the gray level co-occurrence matrix (GLCM) were used to quantify the improvement. For each test dataset, the proposed method was compared with the current gold standard: the N4 algorithm. Additionally, a questionnaire was filled out by two clinical experts to assess the homogeneity and contrast preservation of the 7 T datasets. All four proposed neural networks were able to substantially reduce the B1 -field induced inhomogeneities in T2w 7 T prostate images. By visual inspection, the images clearly look more homogeneous, which is confirmed by the increase in Homogeneity and Energy in the GLCM, and the questionnaire scores from two clinical experts. Occasionally, changes in contrast within the prostate were observed, although much less for the t-Biasf network than for the t-Image network. Further, results on the 3 T dataset demonstrate that the proposed learning based approach is on par with the N4 algorithm. The results demonstrate that the trained networks were capable of reducing the B1 -field induced inhomogeneities for prostate imaging at 7 T. The quantitative evaluation showed that all proposed learning based correction techniques outperformed the N4 algorithm. Of the investigated methods, the single-channel t-Biasf neural network proves most reliable for bias field correction.

Integrating Radiomics With Clinicoradiological Scoring Can Predict High-Risk Patients Who Need Surgery in Crohn's Disease: A Pilot Study

BACKGROUND

Early identification of Crohn's disease (CD) patients at risk for complications could enable targeted surgical referral, but routine magnetic resonance enterography (MRE) has not been definitively correlated with need for surgery. Our objective was to identify computer-extracted image (radiomic) features from MRE associated with risk of surgery in CD and combine them with clinical and radiological assessments to predict time to intervention.

METHODS

This was a retrospective single-center pilot study of CD patients who had an MRE within 3 months prior to initiating medical therapy. Radiomic features were extracted from annotated terminal ileum regions on MRE and combined with clinical variables and radiological assessment (via Simplified Magnetic Resonance Index of Activity scoring for wall thickening, edema, fat stranding, ulcers) in a random forest classifier. The primary endpoint was high- and low-risk groups based on need for surgery within 1 year of MRE. The secondary endpoint was time to surgery after treatment.

RESULTS

Eight radiomic features capturing localized texture heterogeneity within the terminal ileum were significantly associated with risk of surgery within 1 year of treatment (P < .05); yielding a discovery cohort area under the receiver-operating characteristic curve of 0.67 (n = 50) and validation cohort area under the receiver-operating characteristic curve of 0.74 (n = 23). Kaplan-Meier analysis of radiomic features together with clinical variables and Simplified Magnetic Resonance Index of Activity scores yielded the best hazard ratio of 4.13 (P = (7.6 × 10-6) and concordance index of 0.71 in predicting time to surgery after MRE.

CONCLUSIONS

Radiomic features on MRE may be associated with risk of surgery in CD, and in combination with clinicoradiological scoring can yield an accurate prognostic model for time to surgery.

SynthMorph: Learning Contrast-Invariant Registration Without Acquired Images

We introduce a strategy for learning image registration without acquired imaging data, producing powerful networks agnostic to contrast introduced by magnetic resonance imaging (MRI). While classical registration methods accurately estimate the spatial correspondence between images, they solve an optimization problem for every new image pair. Learning-based techniques are fast at test time but limited to registering images with contrasts and geometric content similar to those seen during training. We propose to remove this dependency on training data by leveraging a generative strategy for diverse synthetic label maps and images that exposes networks to a wide range of variability, forcing them to learn more invariant features. This approach results in powerful networks that accurately generalize to a broad array of MRI contrasts. We present extensive experiments with a focus on 3D neuroimaging, showing that this strategy enables robust and accurate registration of arbitrary MRI contrasts even if the target contrast is not seen by the networks during training. We demonstrate registration accuracy surpassing the state of the art both within and across contrasts, using a single model. Critically, training on arbitrary shapes synthesized from noise distributions results in competitive performance, removing the dependency on acquired data of any kind. Additionally, since anatomical label maps are often available for the anatomy of interest, we show that synthesizing images from these dramatically boosts performance, while still avoiding the need for real intensity images. Our code is available at doic https://w3id.org/synthmorph.

Preoperative Assessment of Language Dominance through Combined Resting-State and Task-Based Functional Magnetic Resonance Imaging

Brain lesions in language-related cortical areas remain a challenge in the clinical routine. In recent years, the resting-state fMRI (RS-fMRI) was shown to be a feasible method for preoperative language assessment. The aim of this study was to examine whether language-related resting-state components, which have been obtained using a data-driven independent-component-based identification algorithm, can be supportive in determining language dominance in the left or right hemisphere. Twenty patients suffering from brain lesions close to supposed language-relevant cortical areas were included. RS-fMRI and task-based (TB-fMRI) were performed for the purpose of preoperative language assessment. TB-fMRI included a verb generation task with an appropriate control condition (a syllable switching task) to decompose language-critical and language-supportive processes. Subsequently, the best fitting ICA component for the resting-state language network (RSLN) referential to general linear models (GLMs) of the TB-fMRI (including models with and without linguistic control conditions) was identified using an algorithm based on the Dice index. Thereby, the RSLNs associated with GLMs using a linguistic control condition led to significantly higher laterality indices than GLM baseline contrasts. LIs derived from GLM contrasts with and without control conditions alone did not differ significantly. In general, the results suggest that determining language dominance in the human brain is feasible both with TB-fMRI and RS-fMRI, and in particular, the combination of both approaches yields a higher specificity in preoperative language assessment. Moreover, we can conclude that the choice of the language mapping paradigm is crucial for the mentioned benefits.

Head CT: Toward Making Full Use of the Information the X-Rays Give

Although clinical head CT images are typically interpreted qualitatively, automated methods applied to routine clinical head CTs enable quantitative assessment of brain volume, brain parenchymal fraction, brain radiodensity, and brain radiomass. These metrics gain clinical meaning when viewed relative to a reference database and expressed as quantile regression values. Quantitative imaging data can aid in objective reporting and in the identification of outliers, with possible diagnostic implications. The comparison to a reference database necessitates standardization of head CT imaging parameters and protocols. Future research is needed to learn the effects of virtual monochromatic imaging on the quantitative characteristics of head CT images.

ADVIAN: Alzheimer's Disease VGG-Inspired Attention Network Based on Convolutional Block Attention Module and Multiple Way Data Augmentation

Aim: Alzheimer's disease is a neurodegenerative disease that causes 60-70% of all cases of dementia. This study is to provide a novel method that can identify AD more accurately. Methods: We first propose a VGG-inspired network (VIN) as the backbone network and investigate the use of attention mechanisms. We proposed an Alzheimer's Disease VGG-Inspired Attention Network (ADVIAN), where we integrate convolutional block attention modules on a VIN backbone. Also, 18-way data augmentation is proposed to avoid overfitting. Ten runs of 10-fold cross-validation are carried out to report the unbiased performance. Results: The sensitivity and specificity reach 97.65 ± 1.36 and 97.86 ± 1.55, respectively. Its precision and accuracy are 97.87 ± 1.53 and 97.76 ± 1.13, respectively. The F1 score, MCC, and FMI are obtained as 97.75 ± 1.13, 95.53 ± 2.27, and 97.76 ± 1.13, respectively. The AUC is 0.9852. Conclusion: The proposed ADVIAN gives better results than 11 state-of-the-art methods. Besides, experimental results demonstrate the effectiveness of 18-way data augmentation.

LEARNING MRI CONTRAST-AGNOSTIC REGISTRATION

We introduce a strategy for learning image registration without acquired imaging data, producing powerful networks agnostic to magnetic resonance imaging (MRI) contrast. While classical methods accurately estimate the spatial correspondence between images, they solve an optimization problem for every new image pair. Learning methods are fast at test time but limited to images with contrasts and geometric content similar to those seen during training. We propose to remove this dependency using a generative strategy that exposes networks to a wide range of images synthesized from segmentations during training, forcing them to generalize across contrasts. We show that networks trained within this framework generalize to a broad array of unseen MRI contrasts and surpass classical state-of-the-art brain registration accuracy by up to 12.4 Dice points for a variety of tested contrast combinations. Critically, training on arbitrary shapes synthesized from noise distributions results in competitive performance, removing the dependency on acquired data of any kind. Additionally, since anatomical label maps are often available for the anatomy of interest, we show that synthesizing images from these dramatically boosts performance, while still avoiding the need for real intensity images during training.

Intensity non-uniformity correction in MR imaging using residual cycle generative adversarial network

Correcting or reducing the effects of voxel intensity non-uniformity (INU) within a given tissue type is a crucial issue for quantitative magnetic resonance (MR) image analysis in daily clinical practice. Although having no severe impact on visual diagnosis, the INU can highly degrade the performance of automatic quantitative analysis such as segmentation, registration, feature extraction and radiomics. In this study, we present an advanced deep learning based INU correction algorithm called residual cycle generative adversarial network (res-cycle GAN), which integrates the residual block concept into a cycle-consistent GAN (cycle-GAN). In cycle-GAN, an inverse transformation was implemented between the INU uncorrected and corrected magnetic resonance imaging (MRI) images to constrain the model through forcing the calculation of both an INU corrected MRI and a synthetic corrected MRI. A fully convolution neural network integrating residual blocks was applied in the generator of cycle-GAN to enhance end-to-end raw MRI to INU corrected MRI transformation. A cohort of 55 abdominal patients with T1-weighted MR INU images and their corrections with a clinically established and commonly used method, namely, N4ITK were used as a pair to evaluate the proposed res-cycle GAN based INU correction algorithm. Quantitatively comparisons of normalized mean absolute error (NMAE), peak signal-to-noise ratio (PSNR), normalized cross-correlation (NCC) indices, and spatial non-uniformity (SNU) were made among the proposed method and other approaches. Our res-cycle GAN based method achieved an NMAE of 0.011 ± 0.002, a PSNR of 28.0 ± 1.9 dB, an NCC of 0.970 ± 0.017, and a SNU of 0.298 ± 0.085. Our proposed method has significant improvements (p < 0.05) in NMAE, PSNR, NCC and SNU over other algorithms including conventional GAN and U-net. Once the model is well trained, our approach can automatically generate the corrected MR images in a few minutes, eliminating the need for manual setting of parameters.

A new approach to automatic fetal brain extraction from MRI using a variational level set method

BACKGROUND AND PURPOSE

Appropriate images extracted from the MRI of mothers' wombs can be of great help in the medical diagnosis of fetal abnormalities. As maternal tissue may appear in such images, affecting visualization of myelination of the fetal brain, it is not possible to use methods routinely used for extraction of adult brains for fetal brains. The aim of the present study was to use a variational level set approach to extract fetal brain from T2-weighted MR images of the womb.

METHODS

Coronal T2-weighted images were acquired using fast MRI protocols (to avoid artifacts). The database includes 105 MR images from eight subjects. After correcting the inhomogeneity of the images, the fetal eyes were located, and from that information, the location of the fetus brain was automatically determined. Then, the variational level set was used for fetus brain extraction. The results were analyzed by a clinical specialist (radiologist) and the similarity (Dice and Jaccard coefficients), sensitivity and specificity were calculated.

RESULTS AND CONCLUSIONS

The means of the statistical analysis for the Dice and Jaccard coefficients, sensitivity and specificity, were 99.56%, 96.89%, 95.71%, and 97.96%, respectively. Thus, extraction of fetal brain from MR images was confirmed, both statistically and visually through cross-validation.

Classification of Alzheimer's Disease with and without Imagery using Gradient Boosted Machines and ResNet-50

Background. Alzheimer’s is a disease for which there is no cure. Diagnosing Alzheimer’s disease (AD) early facilitates family planning and cost control. The purpose of this study is to predict the presence of AD using socio-demographic, clinical, and magnetic resonance imaging (MRI) data. Early detection of AD enables family planning and may reduce costs by delaying long-term care. Accurate, non-imagery methods also reduce patient costs. The Open Access Series of Imaging Studies (OASIS-1) cross-sectional MRI data were analyzed. A gradient boosted machine (GBM) predicted the presence of AD as a function of gender, age, education, socioeconomic status (SES), and a mini-mental state exam (MMSE). A residual network with 50 layers (ResNet-50) predicted the clinical dementia rating (CDR) presence and severity from MRI’s (multi-class classification). The GBM achieved a mean 91.3% prediction accuracy (10-fold stratified cross validation) for dichotomous CDR using socio-demographic and MMSE variables. MMSE was the most important feature. ResNet-50 using image generation techniques based on an 80% training set resulted in 98.99% three class prediction accuracy on 4139 images (20% validation set) at Epoch 133 and nearly perfect multi-class predication accuracy on the training set (99.34%). Machine learning methods classify AD with high accuracy. GBM models may help provide initial detection based on non-imagery analysis, while ResNet-50 network models might help identify AD patients automatically prior to provider review.

Automated Estimation of Acute Infarct Volume from Noncontrast Head CT Using Image Intensity Inhomogeneity Correction

Identification of early ischemic changes (EIC) on noncontrast head CT scans performed within the first few hours of stroke onset may have important implications for subsequent treatment, though early stroke is poorly delimited on these studies. Lack of sharp lesion boundary delineation in early infarcts precludes manual volume measures, as well as measures using edge-detection or region-filling algorithms. We wished to test a hypothesis that image intensity inhomogeneity correction may provide a sensitive method for identifying the subtle regional hypodensity which is characteristic of early ischemic infarcts. A digital image analysis algorithm was developed using image intensity inhomogeneity correction (IIC) and intensity thresholding. Two different IIC algorithms (FSL and ITK) were compared. The method was evaluated using simulated infarcts and clinical cases. For synthetic infarcts, measured infarct volumes demonstrated strong correlation to the true lesion volume (for 20% decreased density “infarcts,” Pearson r = 0.998 for both algorithms); both algorithms demonstrated improved accuracy with increasing lesion size and decreasing lesion density. In clinical cases (41 acute infarcts in 30 patients), calculated infarct volumes using FSL IIC correlated with the ASPECTS scores (Pearson r = 0.680) and the admission NIHSS (Pearson r = 0.544). Calculated infarct volumes were highly correlated with the clinical decision to treat with IV-tPA. Image intensity inhomogeneity correction, when applied to noncontrast head CT, provides a tool for image analysis to aid in detection of EIC, as well as to evaluate and guide improvements in scan quality for optimal detection of EIC.

Tissue segmentation-based electron density mapping for MR-only radiotherapy treatment planning of brain using conventional T1-weighted MR images

Purpose

Magnetic resonance imaging (MRI) is the primary modality for targeting brain tumors in radiotherapy treatment planning (RTP). MRI is not directly used for dose calculation since image voxel intensities of MRI are not associated with EDs of tissues as those of computed tomography (CT). The purpose of the present study is to develop and evaluate a tissue segmentation‐based method to generate a synthetic‐CT (sCT) by mapping EDs to corresponding tissues using only T1‐weighted MR images for MR‐only RTP.

Methods

Air regions were contoured in several slices. Then, air, bone, brain, cerebrospinal fluid (CSF), and other soft tissues were automatically segmented with an in‐house algorithm based on edge detection and anatomical information and relative intensity distribution. The intensities of voxels in each segmented tissue were mapped into their CT number range to generate a sCT. Twenty‐five stereotactic radiosurgery and stereotactic ablative radiotherapy patients’ T1‐weighted MRI and coregistered CT images from two centers were retrospectively evaluated. The CT was used as ground truth. Distances between bone contours of the external skull of sCT and CT were measured. The mean error (ME) and mean absolute error (MAE) of electron density represented by standardized CT number was calculated in HU.

Results

The average distance between the contour of the external skull in sCT and the contour in coregistered CT is 1.0 ± 0.2 mm (mean ± 1SD). The ME and MAE differences for air, soft tissue and whole body voxels within external body contours are −4 HU/24 HU, 2 HU/26 HU, and −2 HU/125 HU, respectively.

Conclusions

A MR‐sCT generation technique was developed based on tissue segmentation and voxel‐based tissue ED mapping. The generated sCT is comparable to real CT in terms of anatomical position of tissues and similarity to the ED assignment. This method provides a feasible method to generate sCT for MR‐only radiotherapy treatment planning.

Brain tissue segmentation using improved kernelized rough-fuzzy C-means with spatio-contextual information from MRI

Segmentation of brain tissues from MRI often becomes crucial to properly investigate any region of the brain in order to detect abnormalities. However, the accurate segmentation of the brain tissues is a challenging task as the different tissue regions are usually imprecise, indiscernible, ambiguous, and overlapping. Additionally, different tissue regions are non-linearly separable. Noises and other artifacts may also present in the brain MRI. Therefore, conventional segmentation techniques may not often achieve desired accuracy. To deal those challenges, a robust kernelized rough fuzzy C-means clustering with spatial constraints (KRFCMSC) is proposed in this article for brain tissue segmentation. Here, the brain tissue segmentation from MRI is considered as a clustering of pixels problem. The basic idea behind the proposed technique is the judicious integration of the fuzzy set, rough set, and kernel trick along with spatial constraints (in the form of contextual information) to increase the clustering (segmentation) performance. The use of rough and fuzzy set theory in the clustering process handles the ambiguity, indiscernibility, vagueness and overlappingness of different brain tissue regions. While, the kernel trick increases the chance of linear separability of the complex regions which are otherwise not linearly separable in its original feature space. In order to deal the noisy pixels, here in the clustering process, the spatio-contextual information is introduced from the neighbouring pixels. Experiments are carried out on different real and synthetic benchmark brain MRI datasets (publicly available from Brainweb, and IBSR) without and with added noise. The performance of the proposed method is compared with five other counterpart clustering based segmentation techniques and evaluated using various supervised as well as unsupervised validity indices such as, overall accuracy, precision, recall, kappa, Jaccard, dice, and kernelized Xie-Beni index. Experimental results justify the superiority and robustness of the proposed method over other state-of-the-art methods on both benchmark real life and synthetic brain MRI datasets with and without added noise. Statistical significance of the better segmentation accuracy can be confirmed from the paired t-test results in favour of the proposed method compared to other counterpart methods.

A computer-aided diagnosis system for plus disease in retinopathy of prematurity with structure adaptive segmentation and vessel based features

Retinopathy of Prematurity (ROP) is a blinding disease affecting the retina of low birth-weight preterm infants. Accurate diagnosis of ROP is essential to identify treatment-requiring ROP, which would help to prevent childhood blindness. Plus disease, which characterizes abnormal twisting, widening and branching of the blood vessels, is a significant symptom of treatment requiring ROP. In this paper, we have developed and evaluated a computer-based analysis system for objective assessment of plus disease in ROP, which best mimics the clinical method of disease diagnosis by identifying unique vessel based features. The proposed system consists of an initial segmentation stage, which will efficiently extract blood vessels of varying width and length by utilizing structure adaptive filtering, connectivity analysis and image fusion. The paper proposes the usage of additional retinal features namely leaf node count and vessel density, to portray the abnormal growth and branching of the blood vessels and to complement the commonly used features namely tortuosity and width. The test results show a better classification of plus disease in terms of sensitivity (95%) and specificity (93%), emphasizing the superiority of the proposed segmentation algorithm and vessel-based features. An additional advantage of the proposed system is that the process of selection of relevant vessels for feature extraction is fully automated, which makes the system highly useful to the non-physician graders, owing to the unavailability of a sufficient number of ROP specialists.

A deep learning framework for efficient analysis of breast volume and fibroglandular tissue using MR data with strong artifacts

PURPOSE

The main purpose of this work is to develop, apply, and evaluate an efficient approach for breast density estimation in magnetic resonance imaging data, which contain strong artifacts including intensity inhomogeneities.

METHODS

We present a pipeline for breast density estimation, which consists of intensity inhomogeneity correction, breast volume segmentation, nipple extraction, and fibroglandular tissue segmentation. For the segmentation steps, a well-known deep learning architecture is employed.

RESULTS

The average Dice coefficient for the breast parenchyma is [Formula: see text], which outperforms the classical state-of-the-art approach by a margin of [Formula: see text].

CONCLUSION

The proposed solution is accurate and highly efficient and has potential to be applied for big epidemiological data with thousands of participants.

Normalized averaged range (nAR), a robust quantification method for MPIO-content

Micron-sized paramagnetic iron oxide particles (MPIO) are commonly used as contrast agents in magnetic resonance imaging (MRI) that produce negative contrast enhancement, i.e. darkening, on T2^*-weighted images. However, estimation and quantification of MPIO in vivo is still challenging. This limitation mainly arises from smearing and displacement of the negative contrast of the MPIO, so-called blooming, potentially leading to false-positive detection. Further, the bias field induced by the MR coils also hinders visualization and quantification of the MPIO. To mitigate these drawbacks, a positive contrast image can be generated, for example by using a frequency offset technique, which can significantly improve the accuracy of quantification methods. In this research, we introduce the normalized average range (nAR) as a new way to quantify the relative MPIO content within a study. The method compares the average value of test ROIs to that of a control ROI in range filtered images. The nAR can be used on both positive and negative contrast images. The nAR was tested on agar phantoms containing various MPIO concentrations, and on a rostral migration model for MPIO labeled stem cells in mice. The amount of MPIO was quantified for biased and unbiased data, and both for positive and negative contrast images. In addition, the presence of MPIOs in the olfactory bulb was verified by histology. The results show the nAR can indicate the presence and relative content of MPIO for both negative and positive images. However, the nAR showed slightly higher sensitivity in optimized positive contrast images compared to negative contrast images. In all cases, the bias field played a minor role in the quantification, making debiasing less of a concern. The dependency of the nAR values on the MPIO content in the ROI was further validated histologically. Thus, the nAR provides a robust and reliable tool for quantification of MPIO in mice.

LATE: A Level-Set Method Based on Local Approximation of Taylor Expansion for Segmenting Intensity Inhomogeneous Images

Intensity inhomogeneity is common in real-world images and inevitably leads to many difficulties for accurate image segmentation. Numerous level-set methods have been proposed to segment images with intensity inhomogeneity. However, most of these methods are based on linear approximation, such as locally weighted mean, which may cause problems when handling images with severe intensity inhomogeneities. In this paper, we view segmentation of such images as a nonconvex optimization problem, since the intensity variation in such an image follows a nonlinear distribution. Then, we propose a novel level-set method named local approximation of Taylor expansion (LATE), which is a nonlinear approximation method to solve the nonconvex optimization problem. In LATE, we use the statistical information of the local region as a fidelity term and the differentials of intensity inhomogeneity as an adjusting term to model the approximation function. In particular, since the first-order differential is represented by the variation degree of intensity inhomogeneity, LATE can improve the approximation quality and enhance the local intensity contrast of images with severe intensity inhomogeneity. Moreover, LATE solves the optimization of function fitting by relaxing the constraint condition. In addition, LATE can be viewed as a constraint relaxation of classical methods, such as the region-scalable fitting model and the local intensity clustering model. Finally, the level-set energy functional is constructed based on the Taylor expansion approximation. To validate the effectiveness of our method, we conduct thorough experiments on synthetic and real images. Experimental results show that the proposed method clearly outperforms other solutions in comparison.

Construction of population-specific Indian MRI brain template: Morphometric comparison with Chinese and Caucasian templates

OBJECTIVE

Spatial normalization of brain MR images is highly dependent on the choice of target brain template. Morphological differences caused by factors like genetic and environmental exposures, generates a necessity to construct population specific brain templates. Brain image analysis performed using brain templates from Caucasian population may not be appropriate for non-Caucasian population. In this study, our objective was to construct an Indian brain template from a large population (N = 157 subjects) and compare the morphometric parameters of this template with that of Chinese-56 and MNI-152 templates. In addition, using an independent MRI data of 15 Indian subjects, we also evaluated the potential registration accuracy differences using these three templates.

METHODS

Indian brain template was constructed using iterative routines as per established procedures. We compared our Indian template with standard MNI-152 template and Chinese template by measuring global brain features. We also examined accuracy of registration by aligning 15 new Indian brains to Indian, Chinese and MNI templates. Furthermore, we supported our measurement protocol with inter-rater and intra-rater reliability analysis.

RESULTS

Our results showed that there were significant differences in global brain features of Indian template in comparison with Chinese and MNI brain templates. The results of registration accuracy analysis revealed that fewer deformations are required when Indian brains are registered to Indian template as compared to Chinese and MNI templates.

CONCLUSION

This study concludes that population specific Indian template is likely to be more appropriate for structural and functional image analysis of Indian population.

Automatic Spine Tissue Segmentation from MRI Data Based on Cascade of Boosted Classifiers and Active Appearance Model

The study introduces a novel method for automatic segmentation of vertebral column tissue from MRI images. The paper describes a method that combines multiple stages of Machine Learning techniques to recognize and separate different tissues of the human spine. For the needs of this paper, 50 MRI examinations presenting lumbosacral spine of patients with low back pain were selected. After the initial filtration, automatic vertebrae recognition using Cascade Classifier takes place. Afterwards the main segmentation process using the patch based Active Appearance Model is performed. Obtained results are interpolated using centripetal Catmull–Rom splines. The method was tested on previously unseen vertebrae images segmented manually by 5 physicians. A test validating algorithm convergence per iteration was performed and the Intraclass Correlation Coefficient was calculated. Additionally, the 10-fold cross-validation analysis has been done. Presented method proved to be comparable to the physicians (FF = 90.19 ± 1.01%). Moreover results confirmed a proper algorithm convergence. Automatically segmented area correlated well with manual segmentation for single measurements (r¯=0.8336) and for average measurements (r¯=0.9068) with p = 0.05. The 10-fold cross-validation analysis (FF = 91.37 ± 1.13%) confirmed a good model generalization resulting in practical performance.

Baby brain atlases

The baby brain is constantly changing due to its active neurodevelopment, and research into the baby brain is one of the frontiers in neuroscience. To help guide neuroscientists and clinicians in their investigation of this frontier, maps of the baby brain, which contain a priori knowledge about neurodevelopment and anatomy, are essential. "Brain atlas" in this review refers to a 3D-brain image with a set of reference labels, such as a parcellation map, as the anatomical reference that guides the mapping of the brain. Recent advancements in scanners, sequences, and motion control methodologies enable the creation of various types of high-resolution baby brain atlases. What is becoming clear is that one atlas is not sufficient to characterize the existing knowledge about the anatomical variations, disease-related anatomical alterations, and the variations in time-dependent changes. In this review, the types and roles of the human baby brain MRI atlases that are currently available are described and discussed, and future directions in the field of developmental neuroscience and its clinical applications are proposed. The potential use of disease-based atlases to characterize clinically relevant information, such as clinical labels, in addition to conventional anatomical labels, is also discussed.

Intensity inhomogeneity correction of SD-OCT data using macular flatspace

Images of the retina acquired using optical coherence tomography (OCT) often suffer from intensity inhomogeneity problems that degrade both the quality of the images and the performance of automated algorithms utilized to measure structural changes. This intensity variation has many causes, including off-axis acquisition, signal attenuation, multi-frame averaging, and vignetting, making it difficult to correct the data in a fundamental way. This paper presents a method for inhomogeneity correction by acting to reduce the variability of intensities within each layer. In particular, the N3 algorithm, which is popular in neuroimage analysis, is adapted to work for OCT data. N3 works by sharpening the intensity histogram, which reduces the variation of intensities within different classes. To apply it here, the data are first converted to a standardized space called macular flat space (MFS). MFS allows the intensities within each layer to be more easily normalized by removing the natural curvature of the retina. N3 is then run on the MFS data using a modified smoothing model, which improves the efficiency of the original algorithm. We show that our method more accurately corrects gain fields on synthetic OCT data when compared to running N3 on non-flattened data. It also reduces the overall variability of the intensities within each layer, without sacrificing contrast between layers, and improves the performance of registration between OCT images.

A non-iterative multi-scale approach for intensity inhomogeneity correction in MRI

Intensity inhomogeneity is the prime obstacle for MR image processing like automatic segmentation, registration etc. This complication has strong dependence on the associated acquisition hardware and patient anatomy which recommends retrospective correction. In this paper, a new method is developed for correcting the intensity inhomogeneity using a non-iterative multi-scale approach that doesn't necessitate segmentation and any prior knowledge on the scanner or subject. The proposed algorithm extracts bias field at different scales using a Log-Gabor filter bank followed by smoothing operation. Later, they are combined to fit a third degree polynomial to estimate the bias field. Finally, the corrected image is estimated by performing pixel-wise division of original image and bias field. The performance of the same was tested on BrainWeb simulated data, HCP dataset and is found to provide better performance than the state-of-the-art method, N4. A good agreement between the extracted and ground truth bias field is observed through correlation coefficient on different MR modality images that include T1w, T2w and PD. Significant reduction in coefficient variation and coefficient of joint variation ratios in real data indicate an improved inter-class separation and reduced intra-class intensity variations across white and grey matter tissue regions.

A unified Bayesian-based compensated magnetic resonance imaging

Magnetic resonance (MR) images of higher quality is demanded for helping with more accurate and earlier diagnosis of different diseases. The overall quality of MR images is limited due to the existence of different degradation factors such as (1) MR aberrations due to intrinsic properties of the MR scanner, (2) magnetic field inhomogeneity, and (3) inherent MRI noise. Correcting each MRI degradation factor could be solely useful for the quality enhancement of MR imaging with a limited impact. Here, we propose a unified Bayesian based compensated MR imaging (CMRI) system which jointly corrects for the different aforementioned MR aberrations as well as MR noise and hence generates compensated MR (CMR) images with a higher quality. Testing the proposed CMRI system on both MR physical phantom as well as diffusion weighted and T2 weighted MR imaging data resulted in producing MR images with an overall higher quality that better represents different structures of tissue. The quantitative performance analysis shows a higher Signal to Noise (SNR) and Contrast to Noise (CNR) ratios as well as less Coefficient of Variation (CV) for reconstructed images using the proposed CMRI system compared to the Blind Deconvolution Compensation (BDC) method as state-of-the-art. As such, the proposed CMRI system has potential in improving MR image quality, which is important for accurate and consistent clinical interpretation.

Intensity standardisation of 7T MR images for intensity-based segmentation of the human hypothalamus

The high spatial resolution of 7T MRI enables us to identify subtle volume changes in brain structures, providing potential biomarkers of mental disorders. Most volumetric approaches require that similar intensity values represent similar tissue types across different persons. By applying colour-coding to T1-weighted MP2RAGE images, we found that the high measurement accuracy achieved by high-resolution imaging may be compromised by inter-individual variations in the image intensity. To address this issue, we analysed the performance of five intensity standardisation techniques in high-resolution T1-weighted MP2RAGE images. Twenty images with extreme intensities in the GM and WM were standardised to a representative reference image. We performed a multi-level evaluation with a focus on the hypothalamic region—analysing the intensity histograms as well as the actual MR images, and requiring that the correlation between the whole-brain tissue volumes and subject age be preserved during standardisation. The results were compared with T1 maps. Linear standardisation using subcortical ROIs of GM and WM provided good results for all evaluation criteria: it improved the histogram alignment within the ROIs and the average image intensity within the ROIs and the whole-brain GM and WM areas. This method reduced the inter-individual intensity variation of the hypothalamic boundary by more than half, outperforming all other methods, and kept the original correlation between the GM volume and subject age intact. Mixed results were obtained for the other four methods, which sometimes came at the expense of unwarranted changes in the age-related pattern of the GM volume. The mapping of the T1 relaxation time with the MP2RAGE sequence is advertised as being especially robust to bias field inhomogeneity. We found little evidence that substantiated the T1 map’s theoretical superiority over the T1-weighted images regarding the inter-individual image intensity homogeneity.

Comparative evaluation of performance measures for shading correction in time-lapse fluorescence microscopy

Time-lapse fluorescence microscopy is a valuable technology in cell biology, but it suffers from the inherent problem of intensity inhomogeneity due to uneven illumination or camera nonlinearity, known as shading artefacts. This will lead to inaccurate estimates of single-cell features such as average and total intensity. Numerous shading correction methods have been proposed to remove this effect. In order to compare the performance of different methods, many quantitative performance measures have been developed. However, there is little discussion about which performance measure should be generally applied for evaluation on real data, where the ground truth is absent. In this paper, the state-of-the-art shading correction methods and performance evaluation methods are reviewed. We implement 10 popular shading correction methods on two artificial datasets and four real ones. In order to make an objective comparison between those methods, we employ a number of quantitative performance measures. Extensive validation demonstrates that the coefficient of joint variation (CJV) is the most applicable measure in time-lapse fluorescence images. Based on this measure, we have proposed a novel shading correction method that performs better compared to well-established methods for a range of real data tested.

Quantitative Evaluation of Intensity Inhomogeneity Correction Methods for Structural MR Brain Images

The correction of intensity non-uniformity (INU) in magnetic resonance (MR) images is extremely important to ensure both within-subject and across-subject reliability. Here we tackled the problem of objectively comparing INU correction techniques for T1-weighted images, which are the most commonly used in structural brain imaging. We focused our investigations on the methods integrated in widely used software packages for MR data analysis: FreeSurfer, BrainVoyager, SPM and FSL. We used simulated data to assess the INU fields reconstructed by those methods for controlled inhomogeneity magnitudes and noise levels. For each method, we evaluated a wide range of input parameters and defined an enhanced configuration associated with best reconstruction performance. By comparing enhanced and default configurations, we found that the former often provide much more accurate results. Accordingly, we used enhanced configurations for a more objective comparison between methods. For different levels of INU magnitude and noise, SPM and FSL, which integrate INU correction with brain segmentation, generally outperformed FreeSurfer and BrainVoyager, whose methods are exclusively dedicated to INU correction. Nonetheless, accurate INU field reconstructions can be obtained with FreeSurfer on images with low noise and with BrainVoyager for slow and smooth inhomogeneity profiles. Our study may prove helpful for an accurate selection of the INU correction method to be used based on the characteristics of actual MR data.

Imaging biomarkers in multiple Sclerosis: From image analysis to population imaging

The production of imaging data in medicine increases more rapidly than the capacity of computing models to extract information from it. The grand challenges of better understanding the brain, offering better care for neurological disorders, and stimulating new drug design will not be achieved without significant advances in computational neuroscience. The road to success is to develop a new, generic, computational methodology and to confront and validate this methodology on relevant diseases with adapted computational infrastructures. This new concept sustains the need to build new research paradigms to better understand the natural history of the pathology at the early phase; to better aggregate data that will provide the most complete representation of the pathology in order to better correlate imaging with other relevant features such as clinical, biological or genetic data. In this context, one of the major challenges of neuroimaging in clinical neurosciences is to detect quantitative signs of pathological evolution as early as possible to prevent disease progression, evaluate therapeutic protocols or even better understand and model the natural history of a given neurological pathology. Many diseases encompass brain alterations often not visible on conventional MRI sequences, especially in normal appearing brain tissues (NABT). MRI has often a low specificity for differentiating between possible pathological changes which could help in discriminating between the different pathological stages or grades. The objective of medical image analysis procedures is to define new quantitative neuroimaging biomarkers to track the evolution of the pathology at different levels. This paper illustrates this issue in one acute neuro-inflammatory pathology: Multiple Sclerosis (MS). It exhibits the current medical image analysis approaches and explains how this field of research will evolve in the next decade to integrate larger scale of information at the temporal, cellular, structural and morphological levels.

A subspace-based coil combination method for phased-array magnetic resonance imaging

PURPOSE

Coil-by-coil reconstruction methods are followed by coil combination to obtain a single image representing a spin density map. Typical coil combination methods, such as square-root sum-of-squares and adaptive coil combining, yield images that exhibit spatially varying modulation of image intensity. Existing practice is to first combine coils according to a signal-to-noise criterion, then postprocess to correct intensity inhomogeneity. If inhomogeneity is severe, however, intensity correction methods can yield poor results. The purpose of this article is to present an alternative optimality criterion for coil combination; the resulting procedure yields reduced intensity inhomogeneity while preserving contrast.

THEORY AND METHODS

A minimum mean squared error criterion is adopted for combining coils via a subspace decomposition. Techniques are compared using both simulated and in vivo data.

RESULTS

Experimental results for simulated and in vivo data demonstrate lower bias, higher signal-to-noise ratio (about 7×) and contrast-to-noise ratio (about 2×), compared to existing coil combination techniques.

CONCLUSION

The proposed coil combination method is noniterative and does not require estimation of coil sensitivity maps or image mask; the method is particularly suited to cases where intensity inhomogeneity is too severe for existing approaches.

An efficient level set method for simultaneous intensity inhomogeneity correction and segmentation of MR images

Intensity inhomogeneity (bias field) is a common artefact in magnetic resonance (MR) images, which hinders successful automatic segmentation. In this work, a novel algorithm for simultaneous segmentation and bias field correction is presented. The proposed energy functional allows for explicit regularization of the bias field term, making the model more flexible, which is crucial in presence of strong inhomogeneities. An efficient minimization procedure, attempting to find the global minimum, is applied to the energy functional. The algorithm is evaluated qualitatively and quantitatively using a synthetic example and real MR images of different organs. Comparisons with several state-of-the-art methods demonstrate the superior performance of the proposed technique. Desirable results are obtained even for images with strong and complicated inhomogeneity fields and sparse tissue structures.

Bayesian model selection for pathological neuroimaging data applied to white matter lesion segmentation

In neuroimaging studies, pathologies can present themselves as abnormal intensity patterns. Thus, solutions for detecting abnormal intensities are currently under investigation. As each patient is unique, an unbiased and biologically plausible model of pathological data would have to be able to adapt to the subject's individual presentation. Such a model would provide the means for a better understanding of the underlying biological processes and improve one's ability to define pathologically meaningful imaging biomarkers. With this aim in mind, this work proposes a hierarchical fully unsupervised model selection framework for neuroimaging data which enables the distinction between different types of abnormal image patterns without pathological a priori knowledge. Its application on simulated and clinical data demonstrated the ability to detect abnormal intensity clusters, resulting in a competitive to improved behavior in white matter lesion segmentation when compared to three other freely-available automated methods.

Segmentation and quantification of adipose tissue by magnetic resonance imaging

In this brief review, introductory concepts in animal and human adipose tissue segmentation using proton magnetic resonance imaging (MRI) and computed tomography are summarized in the context of obesity research. Adipose tissue segmentation and quantification using spin relaxation-based (e.g., T1-weighted, T2-weighted), relaxometry-based (e.g., T1-, T2-, T2*-mapping), chemical-shift selective, and chemical-shift encoded water-fat MRI pulse sequences are briefly discussed. The continuing interest to classify subcutaneous and visceral adipose tissue depots into smaller sub-depot compartments is mentioned. The use of a single slice, a stack of slices across a limited anatomical region, or a whole body protocol is considered. Common image post-processing steps and emerging atlas-based automated segmentation techniques are noted. Finally, the article identifies some directions of future research, including a discussion on the growing topic of brown adipose tissue and related segmentation considerations.

Retrospective non-uniform illumination correction techniques in images of tuberculosis

Image pre-processing is highly significant in automated analysis of microscopy images. In this work, non-uniform illumination correction has been attempted using the surface fitting method (SFM), multiple regression method (MRM), and bidirectional empirical mode decomposition (BEMD) in digital microscopy images of tuberculosis (TB). The sputum smear positive and negative images recorded under a standard image acquisition protocol were subjected to illumination correction techniques and evaluated by error and statistical measures. Results show that SFM performs more efficiently than MRM or BEMD. The SFM produced sharp images of TB bacilli with better contrast. To further validate the results, multifractal analysis was performed that showed distinct variation before and after implementation of illumination correction by SFM. Results demonstrate that after illumination correction, there is a 26% increase in the number of bacilli, which aids in classification of the TB images into positive and negative, as TB positivity depends on the count of bacilli.

Quantitative 3D-OCT motion correction with tilt and illumination correction, robust similarity measure and regularization

Variability in illumination, signal quality, tilt and the amount of motion pose challenges for post-processing based 3D-OCT motion correction algorithms. We present an advanced 3D-OCT motion correction algorithm using image registration and orthogonal raster scan patterns aimed at addressing these challenges. An intensity similarity measure using the pseudo Huber norm and a regularization scheme based on a pseudo L0.5 norm are introduced. A two-stage registration approach was developed. In the first stage, only axial motion and axial tilt are coarsely corrected. This result is then used as the starting point for a second stage full optimization. In preprocessing, a bias field estimation based approach to correct illumination differences in the input volumes is employed. Quantitative evaluation was performed using a large set of data acquired from 73 healthy and glaucomatous eyes using SD-OCT systems. OCT volumes of both the optic nerve head and the macula region acquired with three independent orthogonal volume pairs for each location were used to assess reproducibility. The advanced motion correction algorithm using the techniques presented in this paper was compared to a basic algorithm corresponding to an earlier version and to performing no motion correction. Errors in segmentation-based measures such as layer positions, retinal and nerve fiber thickness, as well as the blood vessel pattern were evaluated. The quantitative results consistently show that reproducibility is improved considerably by using the advanced algorithm, which also significantly outperforms the basic algorithm. The mean of the mean absolute retinal thickness difference over all data was 9.9 um without motion correction, 7.1 um using the basic algorithm and 5.0 um using the advanced algorithm. Similarly, the blood vessel likelihood map error is reduced to 69% of the uncorrected error for the basic and to 47% of the uncorrected error for the advanced algorithm. These results demonstrate that our advanced motion correction algorithm has the potential to improve the reliability of quantitative measurements derived from 3D-OCT data substantially.

A hybrid method based on fuzzy clustering and local region-based level set for segmentation of inhomogeneous medical images

medical images are more affected by intensity inhomogeneity rather than noise and outliers. This has a great impact on the efficiency of region-based image segmentation methods, because they rely on homogeneity of intensities in the regions of interest. Meanwhile, initialization and configuration of controlling parameters affect the performance of level set segmentation. To address these problems, this paper proposes a new hybrid method that integrates a local region-based level set method with a variation of fuzzy clustering. Specifically it takes an information fusion approach based on a coarse-to-fine framework that seamlessly fuses local spatial information and gray level information with the information of the local region-based level set method. Also, the controlling parameters of level set are directly computed from fuzzy clustering result. This approach has valuable benefits such as automation, no need to prior knowledge about the region of interest (ROI), robustness on intensity inhomogeneity, automatic adjustment of controlling parameters, insensitivity to initialization, and satisfactory accuracy. So, the contribution of this paper is to provide these advantages together which have not been proposed yet for inhomogeneous medical images. Proposed method was tested on several medical images from different modalities for performance evaluation. Experimental results approve its effectiveness in segmenting medical images in comparison with similar methods.

Axonal transport rate decreased at the onset of optic neuritis in EAE mice

Optic neuritis is frequently the first symptom of multiple sclerosis (MS), an inflammatory demyelinating neurodegenerative disease. Impaired axonal transport has been considered as an early event of neurodegenerative diseases. However, few studies have assessed the integrity of axonal transport in MS or its animal models. We hypothesize that axonal transport impairment occurs at the onset of optic neuritis in experimental autoimmune encephalomyelitis (EAE) mice. In this study, we employed manganese-enhanced MRI (MEMRI) to assess axonal transport in optic nerves in EAE mice at the onset of optic neuritis. Axonal transport was assessed as (a) optic nerve Mn(2+) accumulation rate (in % signal change/h) by measuring the rate of increased total optic nerve signal enhancement, and (b) Mn(2+) transport rate (in mm/h) by measuring the rate of change in optic nerve length enhanced by Mn(2+). Compared to sham-treated healthy mice, Mn(2+) accumulation rate was significantly decreased by 19% and 38% for EAE mice with moderate and severe optic neuritis, respectively. The axonal transport rate of Mn(2+) was significantly decreased by 43% and 65% for EAE mice with moderate and severe optic neuritis, respectively. The degree of axonal transport deficit correlated with the extent of impaired visual function and diminished microtubule-associated tubulins, as well as the severity of inflammation, demyelination, and axonal injury at the onset of optic neuritis.

A new kernel-based fuzzy level set method for automated segmentation of medical images in the presence of intensity inhomogeneity

Researchers recently apply an integrative approach to automate medical image segmentation for benefiting available methods and eliminating their disadvantages. Intensity inhomogeneity is a challenging and open problem in this area, which has received less attention by this approach. It has considerable effects on segmentation accuracy. This paper proposes a new kernel-based fuzzy level set algorithm by an integrative approach to deal with this problem. It can directly evolve from the initial level set obtained by Gaussian Kernel-Based Fuzzy C-Means (GKFCM). The controlling parameters of level set evolution are also estimated from the results of GKFCM. Moreover the proposed algorithm is enhanced with locally regularized evolution based on an image model that describes the composition of real-world images, in which intensity inhomogeneity is assumed as a component of an image. Such improvements make level set manipulation easier and lead to more robust segmentation in intensity inhomogeneity. The proposed algorithm has valuable benefits including automation, invariant of intensity inhomogeneity, and high accuracy. Performance evaluation of the proposed algorithm was carried on medical images from different modalities. The results confirm its effectiveness for medical image segmentation.