An energy minimization method for MS lesion segmentation from T1-w and FLAIR images

In this paper, we extend the multiplicative intrinsic component optimization (MICO) algorithm to multichannel MR image segmentation, with focus on segmentation of multiple sclerosis (MS) lesions. The MICO algorithm was originally proposed by Li et al. in Ref. [1] for normal brain tissue segmentation and intensity inhomogeneity correction of a single channel MR image, which exhibits desirable advantages over other methods for MR image segmentation and intensity inhomogeneity correction in terms of segmentation accuracy and robustness. In this paper, we extend the MICO algorithm to multi-channel MR image segmentation and enable the segmentation of MS lesions. We assign different weights for different channels to control the impact of each channel. The weighted channels allow the enhancement of the impact of the FLAIR image on the segmentation of MS lesions by assigning a larger weight to the FLAIR image channel than the other channels. With the inherent mechanism of estimation of the bias field, our method is able to deal with the intensity inhomogeneity in the input multi-channel MR images. In the application of our method, we only use T1-w and FLAIR images as the input two channel MR images. Experimental results show promising result of our method.     

MR image-based measurement of rates of change in volumes of brain structures. Part I: method and validation

A detailed analysis procedure is described for evaluating rates of volumetric change in brain structures based on structural magnetic resonance (MR) images. In this procedure, a series of image processing tools have been employed to address the problems encountered in measuring rates of change based on structural MR images. These tools include an algorithm for intensity non-uniformity correction, a robust algorithm for three-dimensional image registration with sub-voxel precision and an algorithm for brain tissue segmentation. However, a unique feature in the procedure is the use of a fractional volume model that has been developed to provide a quantitative measure for the partial volume effect. With this model, the fractional constituent tissue volumes are evaluated for voxels at the tissue boundary that manifest partial volume effect, thus allowing tissue boundaries be defined at a sub-voxel level and in an automated fashion. Validation studies are presented on key algorithms including segmentation and registration. An overall assessment of the method is provided through the evaluation of the rates of brain atrophy in a group of normal elderly subjects for which the rate of brain atrophy due to normal aging is predictably small. An application of the method is given in Part II where the rates of brain atrophy in various brain regions are studied in relation to normal aging and Alzheimer's disease.     

Multicontext wavelet-based thresholding segmentation of brain tissues in magnetic resonance images

A novel segmentation method based on wavelet transform is presented for gray matter, white matter and cerebrospinal fluid in thin-sliced single-channel brain magnetic resonance (MR) scans. On the basis of the local image model, multicontext wavelet-based thresholding segmentation (MCWT) is proposed to classify 2D MR data into tissues automatically. In MCWT, the wavelet multiscale transform of local image gray histogram is done, and the gray threshold is gradually revealed from large-scale to small-scale coefficients. Image segmentation is independently performed in each local image to calculate the degree of membership of a pixel to each tissue class. Finally, a strategy is adopted to integrate the intersected outcomes from different local images. The result of the experiment indicates that MCWT outperforms other traditional segmentation methods in classifying brain MR images.     

Automatic segmentation for brain MR images via a convex optimized segmentation and bias field correction coupled model

Accurate segmentation of magnetic resonance (MR) images remains challenging mainly due to the intensity inhomogeneity, which is also commonly known as bias field. Recently active contour models with geometric information constraint have been applied, however, most of them deal with the bias field by using a necessary pre-processing step before segmentation of MR data. This paper presents a novel automatic variational method, which can segment brain MR images meanwhile correcting the bias field when segmenting images with high intensity inhomogeneities. We first define a function for clustering the image pixels in a smaller neighborhood. The cluster centers in this objective function have a multiplicative factor that estimates the bias within the neighborhood. In order to reduce the effect of the noise, the local intensity variations are described by the Gaussian distributions with different means and variances. Then, the objective functions are integrated over the entire domain. In order to obtain the global optimal and make the results independent of the initialization of the algorithm, we reconstructed the energy function to be convex and calculated it by using the Split Bregman theory. A salient advantage of our method is that its result is independent of initialization, which allows robust and fully automated application. Our method is able to estimate the bias of quite general profiles, even in 7T MR images. Moreover, our model can also distinguish regions with similar intensity distribution with different variances. The proposed method has been rigorously validated with images acquired on variety of imaging modalities with promising results.     

结合小波融合和深度学习的脑胶质瘤自动分割

针对水肿区域边界模糊和瘤内结构复杂多变导致的脑胶质瘤分割不精确问题,本文提出了一种基于小波融合和3D-UNet网络的脑胶质瘤磁共振图像自动分割算法.首先,对脑胶质瘤磁共振图像的T1、T1ce、T2、Flair四种模态进行小波融合以及偏置场校正;然后,提取待分类的图像块;再利用提取的图像块训练3D-UNet网络以对图像块中的像素进行分类;最后加载损失率较小的网络模型进行分割,并采用基于连通区域的轮廓提取方法,以降低假阳性率.对57组Brats2018（Brain Tumor Segmentation 2018）磁共振图像测试集进行分割的结果显示,肿瘤的整体、核心和水肿部分的平均分割准确率（DSC）分别达到90.64%、80.74%和86.37%,这表明该算法分割脑胶质瘤准确率较高,与金标准相近.相比多模态图像融合前,该算法在减少输入网络数据量和图像冗余信息的同时,还一定程度上解决了胶质瘤边界模糊、分割不精确的问题,提高了分割的准确度和鲁棒性.     

一种确定极端事件阈值的新方法:随机重排去趋势波动分析方法

将去趋势波动分析法(detrended fluctuation analysis,DFA)和替代数据法相结合,同时引入启发式分割算法和卡方检验,提出了一种确定极端气候事件阈值的新方法,称为随机重排去趋势波动分析(stochastic re-sort detrended fluctuation analysis, S-DFA)方法. 同百分位阈值方法相比,S-DFA方法明确指出了极端事件和非极端事件之间的临界值. 基于中国气象局公布的中国165个国际交换站1961-2006年无缺测的逐日日平均气温资料,利用随机重排去趋势波动分析(S-DFA)方法计算并分析了中国极端低温事件阈值的空间分布特征,并对S-DFA方法在实际资料中的应用进行了检验. 从可预报性的角度给出了极端低温事件综合指标的定义. 这一综合指标将极端低温事件的发生频次和强度综合起来,且兼顾了不同地区各自特有的区域气候背景,进一步说明了综合指标定义的合理性. 基于极端低温事件综合指标的空间分布规律,将中国1961-2006年间极端低温事件分为四个不同等级的地区. 极端低温综合指标整体表现出下降趋势,在20世纪80年代初期之前综合指标的变化具有两个明显的准10年周期,而在这之后则一直处于下降趋势且大大低于平均值,直到90年代中期以后才再次上升至平均值附近. 关键词： 去趋势波动分析法 替代数据法 极端气候事件 阈值     

Application of FT-based MMSE deconvolution method for cerebral blood flow measurement in patients with leukoaraiosis

Introduction

The bolus-tracking (BT) technique is the most popular perfusion-weighted (PW) dynamic susceptibility contrast MRI method used for estimating cerebral blood flow (CBF), cerebral blood volume and mean transit time. The BT technique uses a convolution model that establishes the input–output relationship between blood flow and the vascular tracer concentration. Singular value decomposition (SVD)- and Fourier transform (FT)-based deconvolution methods are popular and widely used for estimating PW MRI parameters. However, from the published literature, it appears that SVD is more widely accepted than other methods. In a previous article, an FT-based minimum mean-squared error (MMSE) technique was proposed and simulation experiments were performed to compare it with the well-established circular SVD (oSVD) method. In this study, the FT-based MMSE method has been used to estimate relative CBF (rCBF) in 13 patients with white matter lesions (WMLs) (leukoaraiosis), and results are compared with the widely used oSVD method.

Materials and Methods

Thirteen patients with leukoaraiosis were imaged on a 1.5-T Siemens whole-body scanner. After acquiring the localizer and structural scans consisting of FLAIR (fluid attenuated with inversion recovery), T₁-weighted and T₂-weighted images, perfusion study was implemented as part of the MRI protocol. For each patient and method, two values were calculated: (a) rCBF for normal white matter (NWM) ROI, obtained by dividing the average CBF value in NWM ROI with average CBF in gray matter (GM) ROI, and (b) rCBF for WML ROI, obtained by dividing the average CBF value in WML ROI with average CBF in GM ROI. Results for the two deconvolution methods were computed.

Results and Discussion

A significant (P<.05) decrease in estimated rCBF was observed in the WML in all the patients using the MMSE method, while for the oSVD method, the decrease was observed in all but one patient. Initial results suggest that the MMSE method is comparable to the oSVD method for estimating rCBF in NMW while it may be better than oSVD for estimating rCBF in lesions of low flow. Studies involving a larger patient population may be required to further validate the findings of this work.     

大鼠创伤性脑水肿模型中近红外光有效检测深度研究

采用近红外光谱技术实现颅脑损伤的无损监测过程中,存在着检测深度不明确的问题.利用Monte Carlo模拟光子在生物组织中的传输过程,建立有效检测深度模型,对光纤探头在大鼠创伤性脑水肿模型中的有效检测深度规律进行了研究.采用不依赖于模板的阈值分割和窄带水平集分割方法,将大鼠头部MRI图像分为头皮、头骨、脑脊液、灰质和白质五部分,建立真实的大鼠脑组织三维模型,使Monte Carlo仿真结果更加准确.改进复杂组织光场分布仿真的tMCimg软件,使其能够实时记录光子在组织中的位置和光子被检测器接收时的能量,从而计算出探头在组织中的有效检测深度.分析了不同光源和检测器的中心距、光源芯径对有效检测深度的影响,结果表明光在大鼠脑组织中的有效检测深度小于或者等于光源和检测器中心距的一半,并随光源芯径的增大逐渐增大.建立大鼠脑水肿模型,验证了仿真结果的正确性.研究结果对于无创脑水肿模型的光纤探头的设计和脑水肿区域的判定有着重要的意义.     

Increased differentiation of intracranial white matter lesions by multispectral 3D-tissue segmentation: preliminary results

MRI is a very sensitive imaging modality, however with relatively low specificity. The aim of this work was to determine the potential of image post-processing using 3D-tissue segmentation technique for identification and quantitative characterization of intracranial lesions primarily in the white matter. Forty subjects participated in this study: 28 patients with brain multiple sclerosis (MS), 6 patients with subcortical ischemic vascular dementia (SIVD), and 6 patients with lacunar white matter infarcts (LI). In routine MR imaging these pathologies may be almost indistinguishable. The 3D-tissue segmentation technique used in this study was based on three input MR images (T(1), T(2)-weighted, and proton density). A modified k-Nearest-Neighbor (k-NN) algorithm optimized for maximum computation speed and high quality segmentation was utilized. In MS lesions, two very distinct subsets were classified using this procedure. Based on the results of segmentation one subset probably represent gliosis, and the other edema and demyelination. In SIVD, the segmented images demonstrated homogeneity, which differentiates SIVD from the heterogeneity observed in MS. This homogeneity was in agreement with the general histological findings. The LI changes pathophysiologically from subacute to chronic. The segmented images closely correlated with these changes, showing a central area of necrosis with cyst formation surrounded by an area that appears like reactive gliosis. In the chronic state, the cyst intensity was similar to that of CSF, while in the subacute stage, the peripheral rim was more prominent. Regional brain lesion load were also obtained on one MS patient to demonstrate the potential use of this technique for lesion load measurements. The majority of lesions were identified in the parietal and occipital lobes. The follow-up study showed qualitatively and quantitatively that the calculated MS load increase was associated with brain atrophy represented by an increase in CSF volume as well as decrease in "normal" brain tissue volumes. Importantly, these results were consistent with the patient's clinical evolution of the disease after a six-month period. In conclusion, these results show there is a potential application for a 3D tissue segmentation technique to characterize white matter lesions with similar intensities on T(2)-weighted MR images. The proposed methodology warrants further clinical investigation and evaluation in a large patient population.     

Multispectral MR images segmentation based on fuzzy knowledge and modified seeded region growing

Magnetic resonance imaging (MRI) is a valuable diagnostic tool in medical science due to its capability for soft-tissue characterization and three-dimensional visualization. One potential application of MRI in clinical practice is brain parenchyma classification and segmentation. Based on fuzzy knowledge and modified seeded region growing, this work proposes a novel image segmentation method, called Fuzzy Knowledge-Based Seeded Region Growing (FKSRG), for multispectral MR images. In this work, fuzzy knowledge includes the fuzzy edge, fuzzy similarity and fuzzy distance, which are obtained from relationships between pixels in multispectral MR images and are applied to the modified seeded regions growing process. In conventional regions merging, the final number of regions is unknown. Therefore, a Target Generation Process is proposed and applied to support conventional regions merging, such that the FKSRG method does not over- or undersegment images. Finally, two image sets, namely, computer-generated phantom images and real MR images, are used in experiments to assess the effectiveness of the proposed FKSRG method. Experimental results demonstrate that the FKSRG method segments multispectral MR images much more effectively than the Functional MRI of the Brain Automated Segmentation Tool, K-means and Support Vector Machine methods.     

Optimization of tissue segmentation of brain MR images based on multispectral 3D feature maps

The purpose of this work was to optimize and increase the accuracy of tissue segmentation of the brain magnetic resonance (MR) images based on multispectral 3D feature maps. We used three sets of MR images as input to the in-house developed semi-automated 3D tissue segmentation algorithm: proton density (PD) and T2-weighted fast spin echo and, T1-weighted spin echo. First, to eliminate the random noise, non-linear anisotropic diffusion type filtering was applied to all the images. Second, to reduce the nonuniformity of the images, we devised and applied a correction algorithm based on uniform phantoms. Following these steps, the qualified observer "seeded" (identified training points) the tissue of interest. To reduce the operator dependent errors, cluster optimization was also used; this clustering algorithm identifies the densest clusters pertaining to the tissues. Finally, the images were segmented using k-NN (k-Nearest Neighborhood) algorithm and a stack of color-coded segmented images were created along with the connectivity algorithm to generate the entire surface of the brain. The application of pre-processing optimization steps substantially improved the 3D tissue segmentation methodology.     

Handcrafted fuzzy rules for tissue classification

This article proposes a handcrafted fuzzy rule-based system for segmentation and identification of different tissue types in magnetic resonance (MR) brain images. The proposed fuzzy system uses a combination of histogram and spatial neighborhood-based features. The intensity variation from one type of tissue to another is gradual at the boundaries due to the inherent nature of the MR signal (MR physics). A fuzzy rule-based approach is expected to better handle these variations and variability in features corresponding to different types of tissues. The proposed segmentation is tested to classify the pixels of the T2-weighted axial MR images of the brain into three primary tissue types: white matter, gray matter and cerebral-spinal fluid. The results are compared with those from manual segmentation by an expert, demonstrating good agreement between them.     

Automatic segmentation of cerebral white matter hyperintensities using only 3D FLAIR images

Magnetic Resonance (MR) white matter hyperintensities have been shown to predict an increased risk of developing cognitive decline. However, their actual role in the conversion to dementia is still not fully understood. Automatic segmentation methods can help in the screening and monitoring of Mild Cognitive Impairment patients who take part in large population-based studies. Most existing segmentation approaches use multimodal MR images. However, multiple acquisitions represent a limitation in terms of both patient comfort and computational complexity of the algorithms. In this work, we propose an automatic lesion segmentation method that uses only three-dimensional fluid-attenuation inversion recovery (FLAIR) images. We use a modified context-sensitive Gaussian mixture model to determine voxel class probabilities, followed by correction of FLAIR artifacts. We evaluate the method against the manual segmentation performed by an experienced neuroradiologist and compare the results with other unimodal segmentation approaches. Finally, we apply our method to the segmentation of multiple sclerosis lesions by using a publicly available benchmark dataset. Results show a similar performance to other state-of-the-art multimodal methods, as well as to the human rater.     

A soft-segmentation visualization scheme for magnetic resonance images

Prevalent visualization tools exploit gray value distribution in images through modified histogram equalization and matching technique, referred to as the window width/window level-based method, to improve visibility and enhance diagnostic value. The window width/window level tool is extensively used in magnetic resonance (MR) images to highlight tissue boundaries during image interpretation. However, the identification of different regions and distinct boundaries between them based on gray-level distribution and displayed intensity levels is extremely difficult because of the large dynamic range of tissue intensities inherent in MR images. We propose a soft-segmentation visualization scheme to generate pixel partitions from the histogram of MR image data using a connectionist approach and then generate selective visual depictions of pixel partitions using pseudo color based on an appropriate fuzzy membership function. By applying the display scheme in clinical examples in this study, we could demonstrate additional overlapping regions between distinct tissue types in healthy and diseased areas (in the brain) that could help improve the tissue characterization ability of MR images.     

Automatic white matter lesion segmentation using an adaptive outlier detection method

White matter (WM) lesions are diffuse WM abnormalities that appear as hyperintense (bright) regions in cranial magnetic resonance imaging (MRI). WM lesions are often observed in older populations and are important indicators of stroke, multiple sclerosis, dementia and other brain-related disorders. In this paper, a new automated method for WM lesions segmentation is presented. In the proposed method, the presence of WM lesions is detected as outliers in the intensity distribution of the fluid-attenuated inversion recovery (FLAIR) MR images using an adaptive outlier detection approach. Outliers are detected using a novel adaptive trimmed mean algorithm and box-whisker plot. In addition, pre- and postprocessing steps are implemented to reduce false positives attributed to MRI artifacts commonly observed in FLAIR sequences. The approach is validated using the cranial MRI sequences of 38 subjects. A significant correlation (R=0.9641, P value=3.12×10(-3)) is observed between the automated approach and manual segmentation by radiologist. The accuracy of the proposed approach was further validated by comparing the lesion volumes computed using the automated approach and lesions manually segmented by an expert radiologist. Finally, the proposed approach is compared against leading lesion segmentation algorithms using a benchmark dataset.     

An automatic cerebellum extraction method in T1-weighted brain MR images using an active contour model with a shape prior

Purpose

The objective of this paper was to automatically segment the cerebellum from T1-weighted human brain magnetic resonance (MR) images.

Materials and Methods

The proposed method constructs a cerebellum template using five sets of 3-T MR imaging (MRI) data, which are used to determine the initial position and the shape prior of the cerebellum for the active contour model. Our formulation includes the active contour model with shape prior, which thereby maintains the shape of the template. The proposed active contour model is sequentially applied to sagittal-, coronal- and transverse-view images. To evaluate the proposed method, it is applied to BrainWeb data and a 3-T MRI data set and compared with FreeSurfer with respect to performance assessment metrics.

Results

The segmented cerebellum was compared with the results from FreeSurfer. Using the manually segmented cerebellum as reference, we measured the average Jaccard coefficients of the proposed method, which were 0.882 and 0.885 for the BrainWeb data and 3-T MRI data set, respectively.

Conclusion

We presented the active contour model with shape prior for extracting the cerebellum from T1-weighted brain MR images. The proposed method yielded a robust and accurate segmentation result.     

A fully automatic three-step liver segmentation method on LDA-based probability maps for multiple contrast MR images

Automatic 3D liver segmentation in magnetic resonance (MR) data sets has proven to be a very challenging task in the domain of medical image analysis. There exist numerous approaches for automatic 3D liver segmentation on computer tomography data sets that have influenced the segmentation of MR images. In contrast to previous approaches to liver segmentation in MR data sets, we use all available MR channel information of different weightings and formulate liver tissue and position probabilities in a probabilistic framework. We apply multiclass linear discriminant analysis as a fast and efficient dimensionality reduction technique and generate probability maps then used for segmentation. We develop a fully automatic three-step 3D segmentation approach based upon a modified region growing approach and a further threshold technique. Finally, we incorporate characteristic prior knowledge to improve the segmentation results. This novel 3D segmentation approach is modularized and can be applied for normal and fat accumulated liver tissue properties.     

Non-locally regularized segmentation of multiple sclerosis lesion from multi-channel MRI data

Segmentation of multiple sclerosis (MS) lesion is important for many neuroimaging studies. In this paper, we propose a novel algorithm for automatic segmentation of MS lesions from multi-channel MR images (T1W, T2W and FLAIR images). The proposed method is an extension of Li et al.'s algorithm in [1], which only segments the normal tissues from T1W images. The proposed method is aimed to segment MS lesions, while normal tissues are also segmented and bias field is estimated to handle intensity inhomogeneities in the images. Another contribution of this paper is the introduction of a nonlocal means technique to achieve spatially regularized segmentation, which overcomes the influence of noise. Experimental results have demonstrated the effectiveness and advantages of the proposed algorithm.     

Automatic measurement of changes in brain volume on consecutive 3D MR images by segmentation propagation

This article presents a technique to automatically measure changes in the volume of a structure of interest in successive 3D magnetic resonance (MR) images and its application in the study of the brain and lateral cerebral ventricles. The only manual step is a segmentation of the structure of interest in the first image. The analysis comprises, first, precise rigid co-registration of the time series of images; second, computation of residual deformations between pairs of images; third, automatic quantification of the volume change, obtained by propagation of the segmentation of the structure of interest through the series of MR images. This approach has been applied to monitor changes in the volume of the brain and lateral cerebral ventricles in a healthy subject and a patient with primary progressive aphasia (PPA). Results are consistent with those obtained by application of the boundary shift integral (BSI) and by stereology in the same subjects.