MRI simulation-based evaluation of image-processing and classification methods

Adaptive anatomical preservation optimal denoising for radiation therapy daily MRI

Low-field magnetic resonance imaging (MRI) has recently been integrated with radiation therapy systems to provide image guidance for daily cancer radiation treatments. The main benefit of the low-field strength is minimal electron return effects. The main disadvantage of low-field strength is increased image noise compared to diagnostic MRIs conducted at 1.5 T or higher. The increased image noise affects both the discernibility of soft tissues and the accuracy of further image processing tasks for both clinical and research applications, such as tumor tracking, feature analysis, image segmentation, and image registration. An innovative method, adaptive anatomical preservation optimal denoising (AAPOD), was developed for optimal image denoising, i.e., to maximally reduce noise while preserving the tissue boundaries. AAPOD employs a series of adaptive nonlocal mean (ANLM) denoising trials with increasing denoising filter strength (i.e., the block similarity filtering parameter in the ANLM algorithm), and then detects the tissue boundary losses on the differences of sequentially denoised images using a zero-crossing edge detection method. The optimal denoising filter strength per voxel is determined by identifying the denoising filter strength value at which boundary losses start to appear around the voxel. The final denoising result is generated by applying the ANLM denoising method with the optimal per-voxel denoising filter strengths. The experimental results demonstrated that AAPOD was capable of reducing noise adaptively and optimally while avoiding tissue boundary losses. AAPOD is useful for improving the quality of MRIs with low-contrast-to-noise ratios and could be applied to other medical imaging modalities, e.g., computed tomography.

BlochSolver: A GPU-optimized fast 3D MRI simulator for experimentally compatible pulse sequences

A magnetic resonance imaging (MRI) simulator, which reproduces MRI experiments using computers, has been developed using two graphic-processor-unit (GPU) boards (GTX 1080). The MRI simulator was developed to run according to pulse sequences used in experiments. Experiments and simulations were performed to demonstrate the usefulness of the MRI simulator for three types of pulse sequences, namely, three-dimensional (3D) gradient-echo, 3D radio-frequency spoiled gradient-echo, and gradient-echo multislice with practical matrix sizes. The results demonstrated that the calculation speed using two GPU boards was typically about 7 TFLOPS and about 14 times faster than the calculation speed using CPUs (two 18-core Xeons). We also found that MR images acquired by experiment could be reproduced using an appropriate number of subvoxels, and that 3D isotropic and two-dimensional multislice imaging experiments for practical matrix sizes could be simulated using the MRI simulator. Therefore, we concluded that such powerful MRI simulators are expected to become an indispensable tool for MRI research and development.

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.

Mixed model phase evolution for correction of magnetic field inhomogeneity effects in 3D quantitative gradient echo-based MRI

PURPOSE

In 3D gradient echo magnetic resonance imaging (MRI), strong field gradients B0_macro are visually observed at air/tissue interfaces. At low spatial resolution in particular, the respective field gradients lead to an apparent increase in intravoxel dephasing, and subsequently, to signal loss or inaccurate R2* estimates. If the strong field gradients are measured, their influence can be removed by postprocessing.

METHODS

Conventional corrections usually assume a linear phase evolution with time. For high macroscopic gradient inhomogeneities near the edge of the brain and at the paranasal sinuses, however, this assumption is often broken. Herein, we explored a novel model that considers both linear and stochastic dependences of the phase evolution with echo time in the presence of weak and strong macroscopic field inhomogeneities. We tested the performance of the model at large field gradients using simulation, phantom, and human in vivo studies.

RESULTS

The performance of the proposed approach was markedly better than the standard correction method, providing a correction equivalent to that of the conventional approach in regions with high signal to noise ratio (SNR > 10), but appearing more robust in regions with low SNR (SNR < 4).

CONCLUSION

The proposed technique shows promise to improve R2* measurements in regions of large susceptibilities. The clinical and research applications still require further investigation.

Technical Note: Development of a 3D printed subresolution sandwich phantom for validation of brain SPECT analysis

PURPOSE

To make an adaptable, head shaped radionuclide phantom to simulate molecular imaging of the brain using clinical acquisition and reconstruction protocols. This will allow the characterization and correction of scanner characteristics, and improve the accuracy of clinical image analysis, including the application of databases of normal subjects.

METHODS

A fused deposition modeling 3D printer was used to create a head shaped phantom made up of transaxial slabs, derived from a simulated MRI dataset. The attenuation of the printed polylactide (PLA), measured by means of the Hounsfield unit on CT scanning, was set to match that of the brain by adjusting the proportion of plastic filament and air (fill ratio). Transmission measurements were made to verify the attenuation of the printed slabs. The radionuclide distribution within the phantom was created by adding (99m)Tc pertechnetate to the ink cartridge of a paper printer and printing images of gray and white matter anatomy, segmented from the same MRI data. The complete subresolution sandwich phantom was assembled from alternate 3D printed slabs and radioactive paper sheets, and then imaged on a dual headed gamma camera to simulate an HMPAO SPECT scan.

RESULTS

Reconstructions of phantom scans successfully used automated ellipse fitting to apply attenuation correction. This removed the variability inherent in manual application of attenuation correction and registration inherent in existing cylindrical phantom designs. The resulting images were assessed visually and by count profiles and found to be similar to those from an existing elliptical PMMA phantom.

CONCLUSIONS

The authors have demonstrated the ability to create physically realistic HMPAO SPECT simulations using a novel head-shaped 3D printed subresolution sandwich method phantom. The phantom can be used to validate all neurological SPECT imaging applications. A simple modification of the phantom design to use thinner slabs would make it suitable for use in PET.

Longitudinal multiple sclerosis lesion segmentation: Resource and challenge

In conjunction with the ISBI 2015 conference, we organized a longitudinal lesion segmentation challenge providing training and test data to registered participants. The training data consisted of five subjects with a mean of 4.4 time-points, and test data of fourteen subjects with a mean of 4.4 time-points. All 82 data sets had the white matter lesions associated with multiple sclerosis delineated by two human expert raters. Eleven teams submitted results using state-of-the-art lesion segmentation algorithms to the challenge, with ten teams presenting their results at the conference. We present a quantitative evaluation comparing the consistency of the two raters as well as exploring the performance of the eleven submitted results in addition to three other lesion segmentation algorithms. The challenge presented three unique opportunities: (1) the sharing of a rich data set; (2) collaboration and comparison of the various avenues of research being pursued in the community; and (3) a review and refinement of the evaluation metrics currently in use. We report on the performance of the challenge participants, as well as the construction and evaluation of a consensus delineation. The image data and manual delineations will continue to be available for download, through an evaluation website2 as a resource for future researchers in the area. This data resource provides a platform to compare existing methods in a fair and consistent manner to each other and multiple manual raters.

Optimizing MR Scan Design for Model-Based ${T}_{1}$ , ${T}_{2}$ Estimation From Steady-State Sequences

Rapid, reliable quantification of MR relaxation parameters T1 and T2 is desirable for many clinical applications. Steady-state sequences such as Spoiled Gradient-Recalled Echo (SPGR) and Dual-Echo Steady-State (DESS) are fast and well-suited for relaxometry because the signals they produce are quite sensitive to T1 and T2 variation. However, T1, T2 estimation with these sequences typically requires multiple scans with varied sets of acquisition parameters. This paper describes a systematic framework for selecting scan types (e.g., combinations of SPGR and DESS scans) and optimizing their respective parameters (e.g., flip angles and repetition times). The method is based on a Cramér-Rao Bound (CRB)-inspired min-max optimization that finds scan parameter combinations that robustly enable precise object parameter estimation. We apply this technique to optimize combinations of SPGR and DESS scans for T1, T2 relaxometry in white matter (WM) and grey matter (GM) regions of the human brain at 3T field strength. Phantom accuracy experiments show that SPGR/DESS scan combinations are in excellent agreement with reference measurements. Phantom precision experiments show that trends in T1,T2 pooled sample standard deviations reflect CRB-based predictions. In vivo experiments show that in WM and GM, T1 and T2 estimates from a pair of optimized DESS scans exhibit precision (but not necessarily accuracy) comparable to that of optimized combinations of SPGR and DESS scans. To our knowledge, T1 maps from DESS acquisitions alone are new. This example application illustrates that scan optimization may help reveal new parameter mapping techniques from combinations of established pulse sequences.

Fast Realistic MRI Simulations Based on Generalized Multi-Pool Exchange Tissue Model

We present MRiLab, a new comprehensive simulator for large-scale realistic MRI simulations on a regular PC equipped with a modern graphical processing unit (GPU). MRiLab combines realistic tissue modeling with numerical virtualization of an MRI system and scanning experiment to enable assessment of a broad range of MRI approaches including advanced quantitative MRI methods inferring microstructure on a sub-voxel level. A flexible representation of tissue microstructure is achieved in MRiLab by employing the generalized tissue model with multiple exchanging water and macromolecular proton pools rather than a system of independent proton isochromats typically used in previous simulators. The computational power needed for simulation of the biologically relevant tissue models in large 3D objects is gained using parallelized execution on GPU. Three simulated and one actual MRI experiments were performed to demonstrate the ability of the new simulator to accommodate a wide variety of voxel composition scenarios and demonstrate detrimental effects of simplified treatment of tissue micro-organization adapted in previous simulators. GPU execution allowed  ∼ 200× improvement in computational speed over standard CPU. As a cross-platform, open-source, extensible environment for customizing virtual MRI experiments, MRiLab streamlines the development of new MRI methods, especially those aiming to infer quantitatively tissue composition and microstructure.

Dissipative particle dynamics of diffusion-NMR requires high Schmidt-numbers

We present an efficient mesoscale model to simulate the diffusion measurement with nuclear magnetic resonance (NMR). On the level of mesoscopic thermal motion of fluid particles, we couple the Bloch equations with dissipative particle dynamics (DPD). Thereby we establish a physically consistent scaling relation between the diffusion constant measured for DPD-particles and the diffusion constant of a real fluid. The latter is based on a splitting into a centre-of-mass contribution represented by DPD, and an internal contribution which is not resolved in the DPD-level of description. As a consequence, simulating the centre-of-mass contribution with DPD requires high Schmidt numbers. After a verification for fundamental pulse sequences, we apply the NMR-DPD method to NMR diffusion measurements of anisotropic fluids, and of fluids restricted by walls of microfluidic channels. For the latter, the free diffusion and the localisation regime are considered.

A Rough Set Bounded Spatially Constrained Asymmetric Gaussian Mixture Model for Image Segmentation

Accurate image segmentation is an important issue in image processing, where Gaussian mixture models play an important part and have been proven effective. However, most Gaussian mixture model (GMM) based methods suffer from one or more limitations, such as limited noise robustness, over-smoothness for segmentations, and lack of flexibility to fit data. In order to address these issues, in this paper, we propose a rough set bounded asymmetric Gaussian mixture model with spatial constraint for image segmentation. First, based on our previous work where each cluster is characterized by three automatically determined rough-fuzzy regions, we partition the target image into three rough regions with two adaptively computed thresholds. Second, a new bounded indicator function is proposed to determine the bounded support regions of the observed data. The bounded indicator and posterior probability of a pixel that belongs to each sub-region is estimated with respect to the rough region where the pixel lies. Third, to further reduce over-smoothness for segmentations, two novel prior factors are proposed that incorporate the spatial information among neighborhood pixels, which are constructed based on the prior and posterior probabilities of the within- and between-clusters, and considers the spatial direction. We compare our algorithm to state-of-the-art segmentation approaches in both synthetic and real images to demonstrate the superior performance of the proposed algorithm.

Bayesian MRI denoising in complex domain

In recent years, several efforts have been done for producing Magnetic Resonance Image scanner with higher magnetic field strength mainly for increasing the Signal to Noise Ratio and the Contrast to Noise Ratio of the acquired images. However, denoising methodologies still play an important role for achieving images neatness. Several denoising algorithms have been presented in literature. Some of them exploit the statistical characteristics of the involved noise, some others project the image in a transformed domain, some others look for geometrical properties of the image. However, the common denominator consists in working in the amplitude domain, i.e. on the gray scale, real valued image. Within this manuscript we propose the idea of performing the noise filtering in the complex domain, i.e. on the real and on the imaginary parts of the acquired images. The advantage of the proposed methodology is that the statistical model of the involved signals is greatly simplified and no approximations are required, together with the full exploitation of the whole acquired signal. More in detail, a Maximum A Posteriori estimator developed for the handling complex data, which adopts Markov Random Fields for modeling the images, is proposed. First results and comparison with other widely adopted denoising filters confirm the validity of the method.

Extended Modality Propagation: Image Synthesis of Pathological Cases

This paper describes a novel generative model for the synthesis of multi-modal medical images of pathological cases based on a single label map. Our model builds upon i) a generative model commonly used for label fusion and multi-atlas patch-based segmentation of healthy anatomical structures, ii) the Modality Propagation iterative strategy used for a spatially-coherent synthesis of subject-specific scans of desired image modalities. The expression Extended Modality Propagation is coined to refer to the extension of Modality Propagation to the synthesis of images of pathological cases. Moreover, image synthesis uncertainty is estimated. An application to Magnetic Resonance Imaging synthesis of glioma-bearing brains is i) validated on the training dataset of a Multimodal Brain Tumor Image Segmentation challenge, ii) compared to the state-of-the-art in glioma image synthesis, and iii) illustrated using the output of two different tumor growth models. Such a generative model allows the generation of a large dataset of synthetic cases, which could prove useful for the training, validation, or benchmarking of image processing algorithms.

Simulation of MR angiography imaging for validation of cerebral arteries segmentation algorithms

BACKGROUND AND OBJECTIVE

Accurate vessel segmentation of magnetic resonance angiography (MRA) images is essential for computer-aided diagnosis of cerebrovascular diseases such as stenosis or aneurysm. The ability of a segmentation algorithm to correctly reproduce the geometry of the arterial system should be expressed quantitatively and observer-independently to ensure objectivism of the evaluation.

METHODS

This paper introduces a methodology for validating vessel segmentation algorithms using a custom-designed MRA simulation framework. For this purpose, a realistic reference model of an intracranial arterial tree was developed based on a real Time-of-Flight (TOF) MRA data set. With this specific geometry blood flow was simulated and a series of TOF images was synthesized using various acquisition protocol parameters and signal-to-noise ratios. The synthesized arterial tree was then reconstructed using a level-set segmentation algorithm available in the Vascular Modeling Toolkit (VMTK). Moreover, to present versatile application of the proposed methodology, validation was also performed for two alternative techniques: a multi-scale vessel enhancement filter and the Chan-Vese variant of the level-set-based approach, as implemented in the Insight Segmentation and Registration Toolkit (ITK). The segmentation results were compared against the reference model.

RESULTS

The accuracy in determining the vessels centerline courses was very high for each tested segmentation algorithm (mean error rate = 5.6% if using VMTK). However, the estimated radii exhibited deviations from ground truth values with mean error rates ranging from 7% up to 79%, depending on the vessel size, image acquisition and segmentation method.

CONCLUSIONS

We demonstrated the practical application of the designed MRA simulator as a reliable tool for quantitative validation of MRA image processing algorithms that provides objective, reproducible results and is observer independent.

Comparison-Based Image Quality Assessment for Selecting Image Restoration Parameters

Image quality assessment (IQA) is traditionally classified into full-reference (FR) IQA, reduced-reference (RR) IQA, and no-reference (NR) IQA according to the amount of information required from the original image. Although NR-IQA and RR-IQA are widely used in practical applications, room for improvement still remains because of the lack of the reference image. Inspired by the fact that in many applications, such as parameter selection for image restoration algorithms, a series of distorted images are available, the authors propose a novel comparison-based IQA (C-IQA) framework. The new comparison-based framework parallels FR-IQA by requiring two input images and resembles NR-IQA by not using the original image. As a result, the new comparison-based approach has more application scenarios than FR-IQA does, and takes greater advantage of the accessible information than the traditional single-input NR-IQA does. Further, C-IQA is compared with other state-of-the-art NR-IQA methods and another RR-IQA method on two widely used IQA databases. Experimental results show that C-IQA outperforms the other methods for parameter selection, and the parameter trimming framework combined with C-IQA saves the computation of iterative image reconstruction up to 80%.

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.

Deep and Structured Robust Information Theoretic Learning for Image Analysis

This paper presents a robust information theoretic (RIT) model to reduce the uncertainties, i.e., missing and noisy labels, in general discriminative data representation tasks. The fundamental pursuit of our model is to simultaneously learn a transformation function and a discriminative classifier that maximize the mutual information of data and their labels in the latent space. In this general paradigm, we, respectively, discuss three types of the RIT implementations with linear subspace embedding, deep transformation, and structured sparse learning. In practice, the RIT and deep RIT are exploited to solve the image categorization task whose performances will be verified on various benchmark data sets. The structured sparse RIT is further applied to a medical image analysis task for brain magnetic resonance image segmentation that allows group-level feature selections on the brain tissues.

Robust Intensity Standardization in Brain Magnetic Resonance Images

The paper is focused on a tiSsue-Based Standardization Technique (SBST) of magnetic resonance (MR) brain images. Magnetic Resonance Imaging intensities have no fixed tissue-specific numeric meaning, even within the same MRI protocol, for the same body region, or even for images of the same patient obtained on the same scanner in different moments. This affects postprocessing tasks such as automatic segmentation or unsupervised/supervised classification methods, which strictly depend on the observed image intensities, compromising the accuracy and efficiency of many image analyses algorithms. A large number of MR images from public databases, belonging to healthy people and to patients with different degrees of neurodegenerative pathology, were employed together with synthetic MRIs. Combining both histogram and tissue-specific intensity information, a correspondence is obtained for each tissue across images. The novelty consists of computing three standardizing transformations for the three main brain tissues, for each tissue class separately. In order to create a continuous intensity mapping, spline smoothing of the overall slightly discontinuous piecewise-linear intensity transformation is performed. The robustness of the technique is assessed in a post hoc manner, by verifying that automatic segmentation of images before and after standardization gives a high overlapping (Dice index >0.9) for each tissue class, even across images coming from different sources. Furthermore, SBST efficacy is tested by evaluating if and how much it increases intertissue discrimination and by assessing gaussianity of tissue gray-level distributions before and after standardization. Some quantitative comparisons to already existing different approaches available in the literature are performed.

Postreconstruction filtering of 3D PET images by using weighted higher-order singular value decomposition

Background

Positron emission tomography (PET) always suffers from high levels of noise due to the constraints of the injected dose and acquisition time, especially in the studies of dynamic PET imaging. To improve the quality of PET image, several approaches have been introduced to suppress noise. However, traditional filters often blur the image edges, or erase small detail, or rely on multiple parameters. In order to solve such problems, nonlocal denoising methods have been adapted to denoise PET images.

Methods

In this paper, we propose to use the weighted higher-order singular value decomposition for PET image denoising. We first modeled the noise in the PET image as Poisson distribution. Then, we transformed the noise to an additive Gaussian noise by use of the anscombe root transformation. Finally, we denoised the transformed image using the proposed higher-order singular value decomposition (HOSVD)-based algorithms. The denoised results were compared with results from some general filters by performing physical phantom and mice studies.

Results

Compared to other commonly used filters, HOSVD-based denoising algorithms can preserve boundaries and quantitative accuracy better. The spatial resolution and the low activity features in PET image also can be preserved by use of HOSVD-based methods. Comparing with the standard HOSVD-based algorithm, the proposed weighted HOSVD algorithm can suppress the stair-step artifact, and the time-consumption is about half of that needed by the Wiener-augmented HOSVD algorithm.

Conclusions

The proposed weighted HOSVD denoising algorithm can suppress noise while better preserving of boundary and quantity in PET images.

Electronic supplementary material

The online version of this article (doi:10.1186/s12938-016-0221-y) contains supplementary material, which is available to authorized users.

Dual-domain denoising in three dimensional magnetic resonance imaging

Denoising is a crucial preprocessing procedure for three dimensional magnetic resonance imaging (3D MRI). Existing denoising methods are predominantly implemented in a single domain, ignoring information in other domains. However, denoising methods are becoming increasingly complex, making analysis and implementation challenging. The present study aimed to develop a dual-domain image denoising (DDID) algorithm for 3D MRI that encapsulates information from the spatial and transform domains. In the present study, the DDID method was used to distinguish signal from noise in the spatial and frequency domains, after which robust accurate noise estimation was introduced for iterative filtering, which is simple and beneficial for computation. In addition, the proposed method was compared quantitatively and qualitatively with existing methods for synthetic and in vivo MRI datasets. The results of the present study suggested that the novel DDID algorithm performed well and provided competitive results, as compared with existing MRI denoising filters.

Adaptive Dimensionality Reduction with Semi-Supervision (AdDReSS): Classifying Multi-Attribute Biomedical Data

Medical diagnostics is often a multi-attribute problem, necessitating sophisticated tools for analyzing high-dimensional biomedical data. Mining this data often results in two crucial bottlenecks: 1) high dimensionality of features used to represent rich biological data and 2) small amounts of labelled training data due to the expense of consulting highly specific medical expertise necessary to assess each study. Currently, no approach that we are aware of has attempted to use active learning in the context of dimensionality reduction approaches for improving the construction of low dimensional representations. We present our novel methodology, AdDReSS (Adaptive Dimensionality Reduction with Semi-Supervision), to demonstrate that fewer labeled instances identified via AL in embedding space are needed for creating a more discriminative embedding representation compared to randomly selected instances. We tested our methodology on a wide variety of domains ranging from prostate gene expression, ovarian proteomic spectra, brain magnetic resonance imaging, and breast histopathology. Across these various high dimensional biomedical datasets with 100+ observations each and all parameters considered, the median classification accuracy across all experiments showed AdDReSS (88.7%) to outperform SSAGE, a SSDR method using random sampling (85.5%), and Graph Embedding (81.5%). Furthermore, we found that embeddings generated via AdDReSS achieved a mean 35.95% improvement in Raghavan efficiency, a measure of learning rate, over SSAGE. Our results demonstrate the value of AdDReSS to provide low dimensional representations of high dimensional biomedical data while achieving higher classification rates with fewer labelled examples as compared to without active learning.

Resolution Improvement in Positron Emission Tomography Using Anatomical Magnetic Resonance Imaging

An ideal imaging system should provide information with high-sensitivity, high spatial, and temporal resolution. Unfortunately, it is not possible to satisfy all of these desired features in a single modality. In this paper, we discuss methods to improve the spatial resolution in positron emission imaging (PET) using a priori information from Magnetic Resonance Imaging (MRI). Our approach uses an image restoration algorithm based on the maximization of mutual information (MMI), which has found significant success for optimizing multimodal image registration. The MMI criterion is used to estimate the parameters in the Sharpness-Constrained Wiener filter. The generated filter is then applied to restore PET images of a realistic digital brain phantom. The resulting restored images show improved resolution and better signal-to-noise ratio compared to the interpolated PET images. We conclude that a Sharpness-Constrained Wiener filter having parameters optimized from a MMI criterion may be useful for restoring spatial resolution in PET based on a priori information from correlated MRI.

Global denoising for 3D MRI

BACKGROUND

Denoising is the primary preprocessing step for subsequent application of MRI. However, most commonly-used patch-based denoising methods are heavily dependent on the degree of patch matching. Due to the large number of voxels in the 3D MRI dataset, the procedure of searching sufficient similarity patches was limited by the empirical compromising between computational efficiency and estimation accuracy, and cannot fulfill the application in multimodal MRI dataset with different SNR and resolutions.

METHODS

In this study, we propose a modified global filtering framework for 3D MRI. For each denoising voxel, the similarity weighting matrix is computed using the reference patch and other patches from the whole dataset. This large weighting matrix is then approximated using the k-means clustering Nyström method to achieve computational viability.

RESULTS

Experiments on both synthetic and in vivo MRI datasets demonstrated that the proposed adaptive Nyström low-rank approximation could achieve competitive estimation compared with exact global filter while reducing the sampling rate by four orders of magnitude. In addition, the corresponding global filter improved patches-based method in both spatial and transform domain.

CONCLUSION

We propose a global denoising framework for 3D MRI which extracts information from the entire dataset to restore each voxel. This large weighting matrix of the global filter is approximated using Nyström low-rank approximation with an adaptive k-means clustering sampling scheme, which significantly reduce the sampling rate as well as the running time. The proposed method is capable of denoising in multimodal MRI dataset and can be used to improve currently used patch-based methods.

Automatic Region-Based Brain Classification of MRI-T1 Data

Image segmentation of medical images is a challenging problem with several still not totally solved issues, such as noise interference and image artifacts. Region-based and histogram-based segmentation methods have been widely used in image segmentation. Problems arise when we use these methods, such as the selection of a suitable threshold value for the histogram-based method and the over-segmentation followed by the time-consuming merge processing in the region-based algorithm. To provide an efficient approach that not only produce better results, but also maintain low computational complexity, a new region dividing based technique is developed for image segmentation, which combines the advantages of both regions-based and histogram-based methods. The proposed method is applied to the challenging applications: Gray matter (GM), White matter (WM) and cerebro-spinal fluid (CSF) segmentation in brain MR Images. The method is evaluated on both simulated and real data, and compared with other segmentation techniques. The obtained results have demonstrated its improved performance and robustness.

Intensity Inhomogeneity Correction of Structural MR Images: A Data-Driven Approach to Define Input Algorithm Parameters

Intensity non-uniformity (INU) in magnetic resonance (MR) imaging is a major issue when conducting analyses of brain structural properties. An inaccurate INU correction may result in qualitative and quantitative misinterpretations. Several INU correction methods exist, whose performance largely depend on the specific parameter settings that need to be chosen by the user. Here we addressed the question of how to select the best input parameters for a specific INU correction algorithm. Our investigation was based on the INU correction algorithm implemented in SPM, but this can be in principle extended to any other algorithm requiring the selection of input parameters. We conducted a comprehensive comparison of indirect metrics for the assessment of INU correction performance, namely the coefficient of variation of white matter (CVWM), the coefficient of variation of gray matter (CVGM), and the coefficient of joint variation between white matter and gray matter (CJV). Using simulated MR data, we observed the CJV to be more accurate than CVWM and CVGM, provided that the noise level in the INU-corrected image was controlled by means of spatial smoothing. Based on the CJV, we developed a data-driven approach for selecting INU correction parameters, which could effectively work on actual MR images. To this end, we implemented an enhanced procedure for the definition of white and gray matter masks, based on which the CJV was calculated. Our approach was validated using actual T1-weighted images collected with 1.5 T, 3 T, and 7 T MR scanners. We found that our procedure can reliably assist the selection of valid INU correction algorithm parameters, thereby contributing to an enhanced inhomogeneity correction in MR images.

3D multimodal spatial fuzzy segmentation of intramuscular connective and adipose tissue from ultrashort TE MR images of calf muscle

PURPOSE

To develop and evaluate an automated algorithm to segment intramuscular adipose (IMAT) and connective (IMCT) tissue from musculoskeletal MRI images acquired with a dual echo Ultrashort TE (UTE) sequence.

THEORY AND METHODS

The dual echo images and calculated structure tensor images are the inputs to the multichannel fuzzy cluster mean (MCFCM) algorithm. Modifications to the basic multichannel fuzzy cluster mean include an adaptive spatial term and bias shading correction. The algorithm was tested on digital phantoms simulating IMAT/IMCT tissue under varying conditions of image noise and bias and on ten subjects with varying amounts of IMAT/IMCT.

RESULTS

The MCFCM including the adaptive spatial term and bias shading correction performed better than the original MCFCM and adaptive spatial MCFCM algorithms. IMAT/IMCT was segmented from the unsmoothed simulated phantom data with a mean Dice coefficient of 0.933 ±0.001 when contrast-to-noise (CNR) was 140 and bias was varied between 30% and 65%. The algorithm yielded accurate in vivo segmentations of IMAT/IMCT with a mean Dice coefficient of 0.977 ±0.066.

CONCLUSION

The proposed algorithm is completely automated and yielded accurate segmentation of intramuscular adipose and connective tissue in the digital phantom and in human calf data. Magn Reson Med 77:870-883, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

A Novel Statistical Approach for Brain MR Images Segmentation Based on Relaxation Times

Brain tissue segmentation in Magnetic Resonance Imaging is useful for a wide range of applications. Classical approaches exploit the gray levels image and implement criteria for differentiating regions. Within this paper a novel approach for brain tissue joint segmentation and classification is presented. Starting from the estimation of proton density and relaxation times, we propose a novel method for identifying the optimal decision regions. The approach exploits the statistical distribution of the involved signals in the complex domain. The technique, compared to classical threshold based ones, is able to globally improve the classification rate. The effectiveness of the approach is evaluated on both simulated and real datasets.

Challenges and Benchmarks in Bioimage Analysis

Similar to the medical imaging community, the bioimaging community has recently realized the need to benchmark various image analysis methods to compare their performance and assess their suitability for specific applications. Challenges sponsored by prestigious conferences have proven to be an effective means of encouraging benchmarking and new algorithm development for a particular type of image data. Bioimage analysis challenges have recently complemented medical image analysis challenges, especially in the case of the International Symposium on Biomedical Imaging (ISBI). This review summarizes recent progress in this respect and describes the general process of designing a bioimage analysis benchmark or challenge, including the proper selection of datasets and evaluation metrics. It also presents examples of specific target applications and biological research tasks that have benefited from these challenges with respect to the performance of automatic image analysis methods that are crucial for the given task. Finally, available benchmarks and challenges in terms of common features, possible classification and implications drawn from the results are analysed.

Improved segmentation of low-contrast lesions using sigmoid edge model

PURPOSE

The intensity profile of an image in the vicinity of a tissue's boundary is modeled by a step/ramp function. However, this assumption does not hold in cases of low-contrast images, heterogeneous tissue textures, and where partial volume effect exists. We propose a hybrid algorithm for segmentation of CT/MR tumors in low-contrast, noisy images having heterogeneous/homogeneous or hyper-/hypo-intense abnormalities. We also model a smoothed noisy intensity profile by a sigmoid function and employ it to find the true location of boundary more accurately.

METHODS

A novel combination of the SVM, watershed, and scattered data approximation algorithms is employed to initially segment a tumor. Small and large abnormalities are treated distinctly. Next, the proposed sigmoid edge model is fitted to the normal profile of the border. The estimated parameters of the model are then utilized to find true boundary of a tissue.

RESULTS

We extensively evaluated our method using synthetic images (contaminated with varying levels of noise) and clinical CT/MR data. Clinical images included 57 CT/MR volumes consisting of small/large tumors, very low-/high-contrast images, liver/brain tumors, and hyper-/hypo-intense abnormalities. We achieved a Dice measure of [Formula: see text] and average symmetric surface distance of [Formula: see text] mm. Regarding IBSR dataset, we fulfilled Jaccard index of [Formula: see text]. The average run-time of our code was [Formula: see text] s.

CONCLUSION

Individual treatment of small and large tumors and boundary correction using the proposed sigmoid edge model can be used to develop a robust tumor segmentation algorithm which deals with any types of tumors.

Realistic simulation of artefacts in diffusion MRI for validating post-processing correction techniques

In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. In the literature, improved alignment is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We apply this framework in three ways. Firstly we assess the most commonly used method for artefact correction, FSL's eddy_correct, and compare it to a recently proposed alternative, eddy. We demonstrate quantitatively that using eddy_correct leads to significant errors in the corrected data, whilst eddy is able to provide much improved correction. Secondly we investigate the datasets required to achieve good correction with eddy, by looking at the minimum number of directions required and comparing the recommended full-sphere acquisitions to equivalent half-sphere protocols. Finally, we investigate the impact of correction quality by examining the fits from microstructure models to real and simulated data.

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.

A Model of Population and Subject (MOPS) Intensities With Application to Multiple Sclerosis Lesion Segmentation

White matter (WM) lesions are thought to play an important role in multiple sclerosis (MS) disease burden. Recent work in the automated segmentation of white matter lesions from magnetic resonance imaging has utilized a model in which lesions are outliers in the distribution of tissue signal intensities across the entire brain of each patient. However, the sensitivity and specificity of lesion detection and segmentation with these approaches have been inadequate. In our analysis, we determined this is due to the substantial overlap between the whole brain signal intensity distribution of lesions and normal tissue. Inspired by the ability of experts to detect lesions based on their local signal intensity characteristics, we propose a new algorithm that achieves lesion and brain tissue segmentation through simultaneous estimation of a spatially global within-the-subject intensity distribution and a spatially local intensity distribution derived from a healthy reference population. We demonstrate that MS lesions can be segmented as outliers from this intensity model of population and subject. We carried out extensive experiments with both synthetic and clinical data, and compared the performance of our new algorithm to those of state-of-the art techniques. We found this new approach leads to a substantial improvement in the sensitivity and specificity of lesion detection and segmentation.

MR image synthesis by contrast learning on neighborhood ensembles

Automatic processing of magnetic resonance images is a vital part of neuroscience research. Yet even the best and most widely used medical image processing methods will not produce consistent results when their input images are acquired with different pulse sequences. Although intensity standardization and image synthesis methods have been introduced to address this problem, their performance remains dependent on knowledge and consistency of the pulse sequences used to acquire the images. In this paper, an image synthesis approach that first estimates the pulse sequence parameters of the subject image is presented. The estimated parameters are then used with a collection of atlas or training images to generate a new atlas image having the same contrast as the subject image. This additional image provides an ideal source from which to synthesize any other target pulse sequence image contained in the atlas. In particular, a nonlinear regression intensity mapping is trained from the new atlas image to the target atlas image and then applied to the subject image to yield the particular target pulse sequence within the atlas. Both intensity standardization and synthesis of missing tissue contrasts can be achieved using this framework. The approach was evaluated on both simulated and real data, and shown to be superior in both intensity standardization and synthesis to other established methods.

An Example-Based Brain MRI Simulation Framework

The simulation of magnetic resonance (MR) images plays an important role in the validation of image analysis algorithms such as image segmentation, due to lack of sufficient ground truth in real MR images. Previous work on MRI simulation has focused on explicitly modeling the MR image formation process. However, because of the overwhelming complexity of MR acquisition these simulations must involve simplifications and approximations that can result in visually unrealistic simulated images. In this work, we describe an example-based simulation framework, which uses an "atlas" consisting of an MR image and its anatomical models derived from the hard segmentation. The relationships between the MR image intensities and its anatomical models are learned using a patch-based regression that implicitly models the physics of the MR image formation. Given the anatomical models of a new brain, a new MR image can be simulated using the learned regression. This approach has been extended to also simulate intensity inhomogeneity artifacts based on the statistical model of training data. Results show that the example based MRI simulation method is capable of simulating different image contrasts and is robust to different choices of atlas. The simulated images resemble real MR images more than simulations produced by a physics-based model.

A Simulation Framework for Quantitative Validation of Artefact Correction in Diffusion MRI

In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. The ultimate aim of these techniques is improved image alignment, yet in the literature this is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We demonstrate the application of this framework by testing one of the most commonly used methods for EC correction, registration of DWIs to b = 0, and reveal the systematic bias this introduces into corrected datasets.

The Metabolic Epicenter of Supratentorial Gliomas: A 1H-MRSI Study

Background:

Assessing the impact of glioma location on prognosis remains elusive. We approached the problem using multivoxel proton magnetic resonance spectroscopic imaging (1H-MRSI) to define a tumor “metabolic epicenter”, and examined the relationship of metabolic epicenter location to survival and histopathological grade.

Methods:

We studied 54 consecutive patients with a supratentorial glioma (astrocytoma or oligodendroglioma, WHO grades II-IV). The metabolic epicenter in each tumor was defined as the 1H-MRSI voxel containing maximum intra-tumoral choline on preoperative imaging. Tumor location was considered the X-Y-Z coordinate position, in a standardized stereotactic space, of the metabolic epicenter. Correlation between epicenter location and survival or grade was assessed.

Results:

Metabolic epicenter location correlated significantly with patient survival for all tumors (r2 = 0.30, p = 0.0002) and astrocytomas alone (r2 = 0.32, p = 0.005). A predictive model based on both metabolic epicenter location and histopathological grade accounted for 70% of the variability in survival, substantially improving on histology alone to predict survival. Location also correlated significantly with grade (r2 = 0.25, p = 0.001): higher grade tumors had a metabolic epicenter closer to the midpoint of the brain.

Conclusions:

The concept of the metabolic epicenter eliminates several problems related to existing methods of classifying glioma location. The location of the metabolic epicenter is strongly correlated with overall survival and histopathological grade, suggesting that it reflects biological factors underlying glioma growth and malignant dedifferentiation. These findings may be clinically relevant to predicting patterns of local glioma recurrence, and in planning resective surgery or radiotherapy.

Partial volume estimation in brain MRI revisited

We propose a fast algorithm to estimate brain tissue concentrations from conventional T1-weighted images based on a Bayesian maximum a posteriori formulation that extends the "mixel" model developed in the 90's. A key observation is the necessity to incorporate additional prior constraints to the "mixel" model for the estimation of plausible concentration maps. Experiments on the ADNI standardized dataset show that global and local brain atrophy measures from the proposed algorithm yield enhanced diagnosis testing value than with several widely used soft tissue labeling methods.

In silico modeling of magnetic resonance flow imaging in complex vascular networks

The paper presents a computational model of magnetic resonance (MR) flow imaging. The model consists of three components. The first component is used to generate complex vascular structures, while the second one provides blood flow characteristics in the generated vascular structures by the lattice Boltzmann method. The third component makes use of the generated vascular structures and flow characteristics to simulate MR flow imaging. To meet computational demands, parallel algorithms are applied in all the components. The proposed approach is verified in three stages. In the first stage, experimental validation is performed by an in vitro phantom. Then, the simulation possibilities of the model are shown. Flow and MR flow imaging in complex vascular structures are presented and evaluated. Finally, the computational performance is tested. Results show that the model is able to reproduce flow behavior in large vascular networks in a relatively short time. Moreover, simulated MR flow images are in accordance with the theoretical considerations and experimental images. The proposed approach is the first such an integrative solution in literature. Moreover, compared to previous works on flow and MR flow imaging, this approach distinguishes itself by its computational efficiency. Such a connection of anatomy, physiology and image formation in a single computer tool could provide an in silico solution to improving our understanding of the processes involved, either considered together or separately.

Regularized estimation of magnitude and phase of multi-coil b1 field via Bloch-Siegert B1 mapping and coil combination optimizations

Parallel excitation requires fast and accurate B1 map estimation. Bloch-Siegert (BS) B1 mapping is very fast and accurate over a large dynamic range. When applied to multi-coil systems, however, this phase-based method may produce low signal-to-noise ratio estimates in low magnitude regions due to localized excitation patterns of parallel excitation systems. Also, the imaging time increases with the number of coils. In this work, we first propose to modify the standard BS B1 mapping sequence so that it avoids the scans required by previous B1 phase estimation methods. A regularized method is then proposed to jointly estimate the magnitude and phase of multi-coil B1 maps from BS B1 mapping data, improving estimation quality by using the prior knowledge of the smoothness of B1 magnitude and phase. Lastly, we use Cramer-Rao lower bound analysis to optimize the coil combinations, to improve the quality of the raw data for B1 estimation. The proposed methods are demonstrated by simulations and phantom experiments.