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.     

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.     

Image-guided regularization level set evolution for MR image segmentation and bias field correction

Magnetic resonance (MR) image segmentation is a crucial step in surgical and treatment planning. In this paper, we propose a level-set-based segmentation method for MR images with intensity inhomogeneous problem. To tackle the initialization sensitivity problem, we propose a new image-guided regularization to restrict the level set function. The maximum a posteriori inference is adopted to unify segmentation and bias field correction within a single framework. Under this framework, both the contour prior and the bias field prior are fully used. As a result, the image intensity inhomogeneity can be well solved. Extensive experiments are provided to evaluate the proposed method, showing significant improvements in both segmentation and bias field correction accuracies as compared with other state-of-the-art approaches.     

An anisotropic images segmentation and bias correction method

Intensity inhomogeneities cause considerable difficulty in the quantitative analysis of magnetic resonance (MR) images. Thus, bias field correction is a necessary step before quantitative analysis of MR data can be undertaken. This paper presents an anisotropic approach to bias correction and segmentation for images with intensity inhomogeneities and noise. Intensity-based methods are usually applied to estimate the bias field; however, most of them only concern the intensity information. When the images have noise or slender topological objects, these methods cannot obtain accurate results or bias fields. We use structure information to construct an anisotropic Gibbs field and combine the anisotropic Gibbs field with the Bayesian framework to segment images while estimating the bias fields. Our method is able to capture bias of quite general profiles. Moreover, it is robust to noise and slender topological objects. The proposed method has been used for images of various modalities with promising results.     

Bias field reduction by localized Lloyd-Max quantization

Bias field reduction is a common problem in medical imaging. A bias field usually manifests itself as a smooth intensity variation across the image. The resulting image inhomogeneity is a severe problem for posterior image processing and analysis techniques such as registration or segmentation. In this article, we present a novel debiasing technique based on localized Lloyd-Max quantization (LMQ). The local bias is modeled as a multiplicative field and is assumed to be slowly varying. The method is based on the assumption that the global, undegraded histogram is characterized by a limited number of gray values. The goal is then to find the discrete intensity values such that spreading those values according to the local bias field reproduces the global histogram as good as possible. We show that our method is capable of efficiently reducing (even strong) bias fields in 3D volumes.     

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.     

A novel active contour model driven by local and global intensity fitting energies

In this paper, we propose a novel hybrid active contour model for image segmentation. In our model, we define a new region-scalable fitting (RSF) energy functional which combines the local and the global image information. The RSF energy functional can not only attract the contour toward object boundaries, but also improve the robustness to initialization of the contours. In order to segment the image fast and accurately, the length term and regularization term is incorporated into the variational level set formulation. Finally, by adopting gradient descent method, the minimization of the energy equation can be given. Due to the new kernel function we defined, our model can cope with intensity inhomogeneity images and less sensitive to the initialization of the contour when compared with the other models. Experimental results demonstrated that the proposed model can also segment both the real and medical images accurately.     

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

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

基于非凸正则化项的合成孔径雷达图像分割新算法

合成孔径雷达图像中乘性噪音的存在使合成孔径雷达图像分割变得非常困难.针对这一难题,本文以提高分割准确度,保护图像的几何结构边缘和提高算法的鲁棒性为目的,提出了一种适用于处理合成孔径雷达图像分割的新模型.新模型结合合成孔径雷达图像的区域和边缘信息,首先通过引入非凸的正则化项,定义了能量泛函;然后极小化能量泛函,建立了水平集函数演化的偏微分方程;最后对水平集演化方程的数值求解,实现了对合成孔径雷达图像感兴趣区域的分割.分别采用仿真图像和实测合成孔径雷达图像对新模型进行验证,结果表明,新模型对合成孔径雷达图像具有很强的边缘定位能力,能使目标区域分割更完整.     

基于非凸正则化项的合成孔径雷达图像分割新算法

合成孔径雷达图像中乘性噪音的存在使合成孔径雷达图像分割变得非常困难.针对这一难题,本文以提高分割准确度,保护图像的几何结构边缘和提高算法的鲁棒性为目的,提出了一种适用于处理合成孔径雷达图像分割的新模型.新模型结合合成孔径雷达图像的区域和边缘信息,首先通过引入非凸的正则化项,定义了能量泛函;然后极小化能量泛函,建立了水平集函数演化的偏微分方程;最后对水平集演化方程的数值求解,实现了对合成孔径雷达图像感兴趣区域的分割.分别采用仿真图像和实测合成孔径雷达图像对新模型进行验证,结果表明,新模型对合成孔径雷达图像具有很强的边缘定位能力,能使目标区域分割更完整.     

Adaptive segmentation of magnetic resonance images with intensity inhomogeneity using level set method

It is a big challenge to segment magnetic resonance (MR) images with intensity inhomogeneity. The widely used segmentation algorithms are region based, which mostly rely on the intensity homogeneity, and could bring inaccurate results. In this paper, we propose a novel region-based active contour model in a variational level set formulation. Based on the fact that intensities in a relatively small local region are separable, a local intensity clustering criterion function is defined. Then, the local function is integrated around the neighborhood center to formulate a global intensity criterion function, which defines the energy term to drive the evolution of the active contour locally. Simultaneously, an intensity fitting term that drives the motion of the active contour globally is added to the energy. In order to segment the image fast and accurately, we utilize a coefficient to make the segmentation adaptive. Finally, the energy is incorporated into a level set formulation with a level set regularization term, and the energy minimization is conducted by a level set evolution process. Experiments on synthetic and real MR images show the effectiveness of our method.     

Magnetic-resonance image segmentation based on improved variable weight multi-resolution Markov random field in undecimated complex wavelet domain

To solve the problem that the magnetic resonance (MR) image has weak boundaries, large amount of information, and low signal-to-noise ratio, we propose an image segmentation method based on the multi-resolution Markov random field (MRMRF) model. The algorithm uses undecimated dual-tree complex wavelet transformation to transform the image into multiple scales. The transformed low-frequency scale histogram is used to improve the initial clustering center of the K-means algorithm, and then other cluster centers are selected according to the maximum distance rule to obtain the coarse-scale segmentation. The results are then segmented by the improved MRMRF model. In order to solve the problem of fuzzy edge segmentation caused by the gray level inhomogeneity of MR image segmentation under the MRMRF model, it is proposed to introduce variable weight parameters in the segmentation process of each scale. Furthermore, the final segmentation results are optimized. We name this algorithm the variable-weight multi-resolution Markov random field (VWMRMRF). The simulation and clinical MR image segmentation verification show that the VWMRMRF algorithm has high segmentation accuracy and robustness, and can accurately and stably achieve low signal-to-noise ratio, weak boundary MR image segmentation.     

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.     

A robust level set method based on local statistical information for noisy image segmentation

The paper presents an improved local region-based active contour model for image segmentation, which is robust to noise. A data fitting energy functional is defined in terms of curves and the energy terms which are based on the differences between the local average intensity and the global intensity means. Then the energy is incorporated into a level set variational formulation, from which a curve evolution equation is derived for energy minimization. And then the level set function is regularized by Gaussian filter to keep smooth and eliminate the re-initialization. By using the local statistical information, the proposed model can handle with noisy images. The proposed model is first presented as a two-phase level set formulation and then extended to a multi-phase one. Experimental results show desirable performances of the proposed model for both noisy synthetic and real images, especially with high level noise.     

Brain tissue classification based on DTI using an improved Fuzzy C-means algorithm with spatial constraints

We present an effective method for brain tissue classification based on diffusion tensor imaging (DTI) data. The method accounts for two main DTI segmentation obstacles: random noise and magnetic field inhomogeneities. In the proposed method, DTI parametric maps were used to resolve intensity inhomogeneities of brain tissue segmentation because they could provide complementary information for tissues and define accurate tissue maps. An improved fuzzy c-means with spatial constraints proposal was used to enhance the noise and artifact robustness of DTI segmentation. Fuzzy c-means clustering with spatial constraints (FCM_S) could effectively segment images corrupted by noise, outliers, and other imaging artifacts. Its effectiveness contributes not only to the introduction of fuzziness for belongingness of each pixel but also to the exploitation of spatial contextual information. We proposed an improved FCM_S applied on DTI parametric maps, which explores the mean and covariance of the feature spatial information for automated segmentation of DTI. The experiments on synthetic images and real-world datasets showed that our proposed algorithms, especially with new spatial constraints, were more effective.     

A fully automated algorithm under modified FCM framework for improved brain MR image segmentation

Automated brain magnetic resonance image (MRI) segmentation is a complex problem especially if accompanied by quality depreciating factors such as intensity inhomogeneity and noise. This article presents a new algorithm for automated segmentation of both normal and diseased brain MRI. An entropy driven homomorphic filtering technique has been employed in this work to remove the bias field. The initial cluster centers are estimated using a proposed algorithm called histogram-based local peak merger using adaptive window. Subsequently, a modified fuzzy c-mean (MFCM) technique using the neighborhood pixel considerations is applied. Finally, a new technique called neighborhood-based membership ambiguity correction (NMAC) has been used for smoothing the boundaries between different tissue classes as well as to remove small pixel level noise, which appear as misclassified pixels even after the MFCM approach. NMAC leads to much sharper boundaries between tissues and, hence, has been found to be highly effective in prominently estimating the tissue and tumor areas in a brain MR scan. The algorithm has been validated against MFCM and FMRIB software library using MRI scans from BrainWeb. Superior results to those achieved with MFCM technique have been observed along with the collateral advantages of fully automatic segmentation, faster computation and faster convergence of the objective function.     

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.     

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.     

基于伪魏格纳分布分解的合成孔径雷达图像目标与阴影分割

在二维魏格纳分布的框架内,针对魏格纳变换的交叉项问题和计算量大的问题,提出了合成孔径雷达图像局部伪魏格纳变换的目标和目标阴影的分割方法.首先,将合成孔径雷达图像进行二维伪魏格纳变换,得到各像素点的二维能量谱图|然后提取各像素点的二维能量谱图对应位置值形成多个不同频段的与原图像同大小的能量谱图|最后,对不同频段的能量谱图采用不同的处理方法后,将各能量谱图相加处理后形成区域标识图像,最终得到原图像的目标和目标阴影分割图像.本文利用该方法对MSTAR切片图像进行了分割试验,并对分割图像与频谱最大值距离或方位分割算法和基于双参量CFAR与隐马尔科夫联合分割算法进行了分割图像对比度对比.实验结果表明,采用本文算法的合成孔径雷达分割图像,对比度明显提高,且保留了目标图像细节.