Markov-Gibbs Random Field Research Articles

Advanced OCTA imaging segmentation: Unsupervised, non-linear retinal vessel detection using modified self-organizing maps and joint MGRF modeling

Background and Objective:This paper proposes a fully automated and unsupervised stochastic segmentation approach using two-level joint Markov-Gibbs Random Field (MGRF) to detect the vascular system from retinal Optical Coherence Tomography Angiography (OCTA) images, which is a critical step in developing Computer-Aided Diagnosis (CAD) systems for detecting retinal diseases. Methods:Using a new probabilistic model based on a Linear Combination of Discrete Gaussian (LCDG), the first level models the appearance of OCTA images and their spatially smoothed images. The parameters of the LCDG model are estimated using a modified Expectation Maximization (EM) algorithm. The second level models the maps of OCTA images, including the vascular system and other retina tissues, using MGRF with analytically estimated parameters from the input images. The proposed segmentation approach employs modified self-organizing maps as a MAP-based optimizer maximizing the joint likelihood and handles the Joint MGRF model in a new, unsupervised way. This approach deviates from traditional stochastic optimization approaches and leverages non-linear optimization to achieve more accurate segmentation results. Results:The proposed segmentation framework is evaluated quantitatively on a dataset of 204 subjects. Achieving 0.92 ± 0.03 Dice similarity coefficient, 0.69 ± 0.25 95-percentile bidirectional Hausdorff distance, and 0.93 ± 0.03 accuracy, confirms the superior performance of the proposed approach. Conclusions:The conclusions drawn from the study highlight the superior performance of the proposed unsupervised and fully automated segmentation approach in detecting the vascular system from OCTA images. This approach not only deviates from traditional methods but also achieves more accurate segmentation results, demonstrating its potential in aiding the development of CAD systems for detecting retinal diseases.

An AI-based novel system for predicting respiratory support in COVID-19 patients through CT imaging analysis

The proposed AI-based diagnostic system aims to predict the respiratory support required for COVID-19 patients by analyzing the correlation between COVID-19 lesions and the level of respiratory support provided to the patients. Computed tomography (CT) imaging will be used to analyze the three levels of respiratory support received by the patient: Level 0 (minimum support), Level 1 (non-invasive support such as soft oxygen), and Level 2 (invasive support such as mechanical ventilation). The system will begin by segmenting the COVID-19 lesions from the CT images and creating an appearance model for each lesion using a 2D, rotation-invariant, Markov–Gibbs random field (MGRF) model. Three MGRF-based models will be created, one for each level of respiratory support. This suggests that the system will be able to differentiate between different levels of severity in COVID-19 patients. The system will decide for each patient using a neural network-based fusion system, which combines the estimates of the Gibbs energy from the three MGRF-based models. The proposed system were assessed using 307 COVID-19-infected patients, achieving an accuracy of 97.72%pm 1.57, a sensitivity of 97.76%pm 4.08, and a specificity of 98.87%pm 2.09, indicating a high level of prediction accuracy.

Spatial data modeling by means of Gibbs–Markov random fields based on a generalized planar rotator model

We introduce a Gibbs–Markov random field for spatial data on Cartesian grids, which is based on the generalized planar rotator (GPR) model. The GPR model generalizes the recently proposed modified planar rotator (MPR) model by including in the Hamiltonian additional terms that better capture realistic features of spatial data, such as smoothness, non-Gaussianity, and geometric anisotropy. In particular, the GPR model can include an infinite number of higher-order harmonics with exponentially vanishing interaction strength, directional dependence of the bilinear interaction term between nearest neighbors, longer-distance neighbor interactions, and two types of an external bias field. Hence, in contrast with the single-parameter MPR, the GPR model features five additional parameters: the number n of higher-order terms and the decay-rate parameter α, the exchange anisotropy parameter Jnn, the further-neighbor interaction coupling Jfn, and the external field (bias) parameters K (or K′). We present numerical tests on various synthetic and real data which demonstrate the effects of the respective terms on the model’s prediction performance, and we discuss these results in connection with the data properties.

Segmentation of Infant Brain Using Nonnegative Matrix Factorization

This study develops an atlas-based automated framework for segmenting infants’ brains from magnetic resonance imaging (MRI). For the accurate segmentation of different structures of an infant’s brain at the isointense age (6–12 months), our framework integrates features of diffusion tensor imaging (DTI) (e.g., the fractional anisotropy (FA)). A brain diffusion tensor (DT) image and its region map are considered samples of a Markov–Gibbs random field (MGRF) that jointly models visual appearance, shape, and spatial homogeneity of a goal structure. The visual appearance is modeled with an empirical distribution of the probability of the DTI features, fused by their nonnegative matrix factorization (NMF) and allocation to data clusters. Projecting an initial high-dimensional feature space onto a low-dimensional space of the significant fused features with the NMF allows for better separation of the goal structure and its background. The cluster centers in the latter space are determined at the training stage by the K-means clustering. In order to adapt to large infant brain inhomogeneities and segment the brain images more accurately, appearance descriptors of both the first-order and second-order are taken into account in the fused NMF feature space. Additionally, a second-order MGRF model is used to describe the appearance based on the voxel intensities and their pairwise spatial dependencies. An adaptive shape prior that is spatially variant is constructed from a training set of co-aligned images, forming an atlas database. Moreover, the spatial homogeneity of the shape is described with a spatially uniform 3D MGRF of the second-order for region labels. In vivo experiments on nine infant datasets showed promising results in terms of the accuracy, which was computed using three metrics: the 95-percentile modified Hausdorff distance (MHD), the Dice similarity coefficient (DSC), and the absolute volume difference (AVD). Both the quantitative and visual assessments confirm that integrating the proposed NMF-fused DTI feature and intensity MGRF models of visual appearance, the adaptive shape prior, and the shape homogeneity MGRF model is promising in segmenting the infant brain DTI.

Automated Diagnosis of Optical Coherence Tomography Angiography (OCTA) Based on Machine Learning Techniques.

Diabetic retinopathy (DR) refers to the ophthalmological complications of diabetes mellitus. It is primarily a disease of the retinal vasculature that can lead to vision loss. Optical coherence tomography angiography (OCTA) demonstrates the ability to detect the changes in the retinal vascular system, which can help in the early detection of DR. In this paper, we describe a novel framework that can detect DR from OCTA based on capturing the appearance and morphological markers of the retinal vascular system. This new framework consists of the following main steps: (1) extracting retinal vascular system from OCTA images based on using joint Markov-Gibbs Random Field (MGRF) model to model the appearance of OCTA images and (2) estimating the distance map inside the extracted vascular system to be used as imaging markers that describe the morphology of the retinal vascular (RV) system. The OCTA images, extracted vascular system, and the RV-estimated distance map is then composed into a three-dimensional matrix to be used as an input to a convolutional neural network (CNN). The main motivation for using this data representation is that it combines the low-level data as well as high-level processed data to allow the CNN to capture significant features to increase its ability to distinguish DR from the normal retina. This has been applied on multi-scale levels to include the original full dimension images as well as sub-images extracted from the original OCTA images. The proposed approach was tested on in-vivo data using about 91 patients, which were qualitatively graded by retinal experts. In addition, it was quantitatively validated using datasets based on three metrics: sensitivity, specificity, and overall accuracy. Results showed the capability of the proposed approach, outperforming the current deep learning as well as features-based detecting DR approaches.

The Role of 3D CT Imaging in the Accurate Diagnosis of Lung Function in Coronavirus Patients.

Early grading of coronavirus disease 2019 (COVID-19), as well as ventilator support machines, are prime ways to help the world fight this virus and reduce the mortality rate. To reduce the burden on physicians, we developed an automatic Computer-Aided Diagnostic (CAD) system to grade COVID-19 from Computed Tomography (CT) images. This system segments the lung region from chest CT scans using an unsupervised approach based on an appearance model, followed by 3D rotation invariant Markov–Gibbs Random Field (MGRF)-based morphological constraints. This system analyzes the segmented lung and generates precise, analytical imaging markers by estimating the MGRF-based analytical potentials. Three Gibbs energy markers were extracted from each CT scan by tuning the MGRF parameters on each lesion separately. The latter were healthy/mild, moderate, and severe lesions. To represent these markers more reliably, a Cumulative Distribution Function (CDF) was generated, then statistical markers were extracted from it, namely, 10th through 90th CDF percentiles with 10% increments. Subsequently, the three extracted markers were combined together and fed into a backpropagation neural network to make the diagnosis. The developed system was assessed on 76 COVID-19-infected patients using two metrics, namely, accuracy and Kappa. In this paper, the proposed system was trained and tested by three approaches. In the first approach, the MGRF model was trained and tested on the lungs. This approach achieved % accuracy and % kappa. In the second approach, we trained the MGRF model on the lesions and tested it on the lungs. This approach achieved % accuracy and % kappa. Finally, we trained and tested the MGRF model on lesions. It achieved 100% accuracy and 100% kappa. The results reported in this paper show the ability of the developed system to accurately grade COVID-19 lesions compared to other machine learning classifiers, such as k-Nearest Neighbor (KNN), decision tree, naïve Bayes, and random forest.

Computer-Assisted Image Processing System for Early Assessment of Lung Nodule Malignancy.

Simple SummaryLung cancer is the second most common cancer in men after prostate cancer and in women after breast cancer, but it is the leading cause of cancer death among both genders. This manuscript proposes a new computer-aided diagnosis system that uses only a single computed tomography scan to diagnose the pulmonary nodule as benign or malignant. This system helps in the early detection of the pulmonary nodules and shows its ability to identify the pulmonary nodules precisely.Lung cancer is one of the most dreadful cancers, and its detection in the early stage is very important and challenging. This manuscript proposes a new computer-aided diagnosis system for lung cancer diagnosis from chest computed tomography scans. The proposed system extracts two different kinds of features, namely, appearance features and shape features. For the appearance features, a Histogram of oriented gradients, a Multi-view analytical Local Binary Pattern, and a Markov Gibbs Random Field are developed to give a good description of the lung nodule texture, which is one of the main distinguishing characteristics between benign and malignant nodules. For the shape features, Multi-view Peripheral Sum Curvature Scale Space, Spherical Harmonics Expansion, and a group of some fundamental morphological features are implemented to describe the outer contour complexity of the nodules, which is main factor in lung nodule diagnosis. Each feature is fed into a stacked auto-encoder followed by a soft-max classifier to generate the initial malignancy probability. Finally, all these probabilities are combined together and fed to the last network to give the final diagnosis. The system is validated using 727 nodules which are subset from the Lung Image Database Consortium (LIDC) dataset. The system shows very high performance measures and achieves , , and for the accuracy, sensitivity, and specificity, respectively. This high performance shows the ability of the system to distinguish between the malignant and benign nodules precisely.

A Novel Computer-Aided Diagnostic System for Early Detection of Diabetic Retinopathy Using 3D-OCT Higher-Order Spatial Appearance Model.

Early diagnosis of diabetic retinopathy (DR) is of critical importance to suppress severe damage to the retina and/or vision loss. In this study, an optical coherence tomography (OCT)-based computer-aided diagnosis (CAD) method is proposed to detect DR early using structural 3D retinal scans. This system uses prior shape knowledge to automatically segment all retinal layers of the 3D-OCT scans using an adaptive, appearance-based method. After the segmentation step, novel texture features are extracted from the segmented layers of the OCT B-scans volume for DR diagnosis. For every layer, Markov–Gibbs random field (MGRF) model is used to extract the 2nd-order reflectivity. In order to represent the extracted image-derived features, we employ cumulative distribution function (CDF) descriptors. For layer-wise classification in 3D volume, using the extracted Gibbs energy feature, an artificial neural network (ANN) is fed the extracted feature for every layer. Finally, the classification outputs for all twelve layers are fused using a majority voting schema for global subject diagnosis. A cohort of 188 3D-OCT subjects are used for system evaluation using different k-fold validation techniques and different validation metrics. Accuracy of , , and are achieved using 4-, 5-, and 10-fold cross-validation, respectively. Additional comparison with deep learning networks, which represent the state-of-the-art, documented the promise of our system’s ability to diagnose the DR early.

Precise Segmentation of COVID-19 Infected Lung from CT Images Based on Adaptive First-Order Appearance Model with Morphological/Anatomical Constraints.

A new segmentation technique is introduced for delineating the lung region in 3D computed tomography (CT) images. To accurately model the distribution of Hounsfield scale values within both chest and lung regions, a new probabilistic model is developed that depends on a linear combination of Gaussian (LCG). Moreover, we modified the conventional expectation-maximization (EM) algorithm to be run in a sequential way to estimate both the dominant Gaussian components (one for the lung region and one for the chest region) and the subdominant Gaussian components, which are used to refine the final estimated joint density. To estimate the marginal density from the mixed density, a modified k-means clustering approach is employed to classify the Gaussian subdominant components to determine which components belong properly to a lung and which components belong to a chest. The initial segmentation, based on the LCG-model, is then refined by the imposition of 3D morphological constraints based on a 3D Markov–Gibbs random field (MGRF) with analytically estimated potentials. The proposed approach was tested on CT data from 32 coronavirus disease 2019 (COVID-19) patients. Segmentation quality was quantitatively evaluated using four metrics: Dice similarity coefficient (DSC), overlap coefficient, 95th-percentile bidirectional Hausdorff distance (BHD), and absolute lung volume difference (ALVD), and it achieved , , , and , respectively. The reported results showed the capability of the proposed approach to accurately segment healthy lung tissues in addition to pathological lung tissues caused by COVID-19, outperforming four current, state-of-the-art deep learning-based lung segmentation approaches.

Early assessment of lung function in coronavirus patients using invariant markers from chest X-rays images

The primary goal of this manuscript is to develop a computer assisted diagnostic (CAD) system to assess pulmonary function and risk of mortality in patients with coronavirus disease 2019 (COVID-19). The CAD system processes chest X-ray data and provides accurate, objective imaging markers to assist in the determination of patients with a higher risk of death and thus are more likely to require mechanical ventilation and/or more intensive clinical care.To obtain an accurate stochastic model that has the ability to detect the severity of lung infection, we develop a second-order Markov-Gibbs random field (MGRF) invariant under rigid transformation (translation or rotation of the image) as well as scale (i.e., pixel size). The parameters of the MGRF model are learned automatically, given a training set of X-ray images with affected lung regions labeled. An X-ray input to the system undergoes pre-processing to correct for non-uniformity of illumination and to delimit the boundary of the lung, using either a fully-automated segmentation routine or manual delineation provided by the radiologist, prior to the diagnosis. The steps of the proposed methodology are: (i) estimate the Gibbs energy at several different radii to describe the inhomogeneity in lung infection; (ii) compute the cumulative distribution function (CDF) as a new representation to describe the local inhomogeneity in the infected region of lung; and (iii) input the CDFs to a new neural network-based fusion system to determine whether the severity of lung infection is low or high. This approach is tested on 200 clinical X-rays from 200 COVID-19 positive patients, 100 of whom died and 100 who recovered using multiple training/testing processes including leave-one-subject-out (LOSO), tenfold, fourfold, and twofold cross-validation tests. The Gibbs energy for lung pathology was estimated at three concentric rings of increasing radii. The accuracy and Dice similarity coefficient (DSC) of the system steadily improved as the radius increased. The overall CAD system combined the estimated Gibbs energy information from all radii and achieved a sensitivity, specificity, accuracy, and DSC of 100%, 97% ± 3%, 98% ± 2%, and 98% ± 2%, respectively, by twofold cross validation. Alternative classification algorithms, including support vector machine, random forest, naive Bayes classifier, K-nearest neighbors, and decision trees all produced inferior results compared to the proposed neural network used in this CAD system. The experiments demonstrate the feasibility of the proposed system as a novel tool to objectively assess disease severity and predict mortality in COVID-19 patients. The proposed tool can assist physicians to determine which patients might require more intensive clinical care, such a mechanical respiratory support.

Precise higher-order reflectivity and morphology models for early diagnosis of diabetic retinopathy using OCT images

This study proposes a novel computer assisted diagnostic (CAD) system for early diagnosis of diabetic retinopathy (DR) using optical coherence tomography (OCT) B-scans. The CAD system is based on fusing novel OCT markers that describe both the morphology/anatomy and the reflectivity of retinal layers to improve DR diagnosis. This system separates retinal layers automatically using a segmentation approach based on an adaptive appearance and their prior shape information. High-order morphological and novel reflectivity markers are extracted from individual segmented layers. Namely, the morphological markers are layer thickness and tortuosity while the reflectivity markers are the 1st-order reflectivity of the layer in addition to local and global high-order reflectivity based on Markov-Gibbs random field (MGRF) and gray-level co-occurrence matrix (GLCM), respectively. The extracted image-derived markers are represented using cumulative distribution function (CDF) descriptors. The constructed CDFs are then described using their statistical measures, i.e., the 10th through 90th percentiles with a 10% increment. For individual layer classification, each extracted descriptor of a given layer is fed to a support vector machine (SVM) classifier with a linear kernel. The results of the four classifiers are then fused using a backpropagation neural network (BNN) to diagnose each retinal layer. For global subject diagnosis, classification outputs (probabilities) of the twelve layers are fused using another BNN to make the final diagnosis of the B-scan. This system is validated and tested on 130 patients, with two scans for both eyes (i.e. 260 OCT images), with a balanced number of normal and DR subjects using different validation metrics: 2-folds, 4-folds, 10-folds, and leave-one-subject-out (LOSO) cross-validation approaches. The performance of the proposed system was evaluated using sensitivity, specificity, F1-score, and accuracy metrics. The system’s performance after the fusion of these different markers showed better performance compared with individual markers and other machine learning fusion methods. Namely, it achieved 96.15%, 99.23%, 97.66%, and 97.69%, respectively, using the LOSO cross-validation technique. The reported results, based on the integration of morphology and reflectivity markers and by using state-of-the-art machine learning classifications, demonstrate the ability of the proposed system to diagnose the DR early.

A novel technology to integrate imaging and clinical markers for non-invasive diagnosis of lung cancer

This study presents a non-invasive, automated, clinical diagnostic system for early diagnosis of lung cancer that integrates imaging data from a single computed tomography scan and breath bio-markers obtained from a single exhaled breath to quickly and accurately classify lung nodules. CT imaging and breath volatile organic compounds data were collected from 47 patients. Spherical Harmonics-based shape features to quantify the shape complexity of the pulmonary nodules, 7th-Order Markov Gibbs Random Field based appearance model to describe the spatial non-homogeneities in the pulmonary nodule, and volumetric features (size) of pulmonary nodules were calculated from CT images. 27 VOCs in exhaled breath were captured by a micro-reactor approach and quantied using mass spectrometry. CT and breath markers were input into a deep-learning autoencoder classifier with a leave-one-subject-out cross validation for nodule classification. To mitigate the limitation of a small sample size and validate the methodology for individual markers, retrospective CT scans from 467 patients with 727 pulmonary nodules, and breath samples from 504 patients were analyzed. The CAD system achieved 97.8% accuracy, 97.3% sensitivity, 100% specificity, and 99.1% area under curve in classifying pulmonary nodules.

A novel 3D segmentation approach for extracting retinal layers from optical coherence tomography images.

Accurate segmentation of retinal layers of the eye in 3D Optical Coherence Tomography (OCT) data provides relevant information for clinical diagnosis. This manuscript describes a 3D segmentation approach that uses an adaptive patient-specific retinal atlas, as well as an appearance model for 3D OCT data. To reconstruct the atlas of 3D retinal scan, the central area of the macula (macula mid-area) where the fovea could be clearly identified, was segmented initially. Markov Gibbs Random Field (MGRF) including intensity, spatial information, and shape of 12 retinal layers were used to segment the selected area of retinal fovea. A set of coregistered OCT scans that were gathered from 200 different individuals were used to build a 2D shape prior. This shape prior was adapted subsequently to the first order appearance and second order spatial interaction MGRF model. After segmenting the center of the macula "foveal area", the labels and appearances of the layers that were segmented were utilized to segment the adjacent slices. The final step was repeated recursively until a 3D OCT scan of the patient was segmented. This approach was tested in 50 patients with normal and with ocular pathological conditions. The segmentation was compared to a manually segmented ground truth. The results were verified by clinical retinal experts. Dice Similarity Coefficient (DSC), 95% bidirectional modified Hausdorff Distance (HD), Unsigned Mean Surface Position Error (MSPE), and Average Volume Difference (AVD) metrics were used to quantify the performance of the proposed approach. The proposed approach was proved to be more accurate than the current state-of-the-art 3D OCT approaches. The proposed approach has the advantage of segmenting all the 12 retinal layers rapidly and more accurately than current state-of-the-art 3D OCT approaches.

Efficient and Scalable Approach to Equilibrium Conditional Simulation of Gibbs Markov Random Fields

We study the performance of an automated hybrid Monte Carlo (HMC) approach for conditional simulation of a recently proposed, single-parameter Gibbs Markov random field. This is based on a modified version of the planar rotator (MPR) model and is used for efficient gap filling in gridded data. HMC combines the deterministic over-relaxation method and the stochastic Metropolis update with dynamically adjusted restriction and performs automatic detection of the crossover to the targeted equilibrium state. We focus on the ability of the algorithm to efficiently drive the system to equilibrium at very low temperatures even with sparse conditioning data. These conditions are the most challenging computationally, requiring extremely long relaxation times if simulated by means of the standard Metropolis algorithm. We demonstrate that HMC has considerable benefits in terms of both computational efficiency and prediction performance of the MPR method.

Accurate Segmentation of Cerebrovasculature From TOF-MRA Images Using Appearance Descriptors

Analyzing cerebrovascular changes can significantly lead to not only detecting the presence of serious diseases e.g., hypertension and dementia, but also tracking their progress. Such analysis could be better performed using Time-of-Flight Magnetic Resonance Angiography (ToF–MRA) images, but this requires accurate segmentation of the cerebral vasculature from the surroundings. To achieve this goal, we propose a fully automated cerebral vasculature segmentation approach based on extracting both prior and current appearance features that have the ability to capture the appearance of macro and micro-vessels in ToF–MRA. The appearance prior is modeled with a novel translation and rotation invariant Markov-Gibbs Random Field (MGRF) of voxel intensities with pairwise interaction analytically identified from a set of training data sets. The appearance of the cerebral vasculature is also represented with a marginal probability distribution of voxel intensities by using a Linear Combination of Discrete Gaussians (LCDG) that its parameters are estimated by using a modified Expectation-Maximization (EM) algorithm. The extracted appearance features are separable and can be classified by any classifier, as demonstrated by our segmentation results. To validate the accuracy of our algorithm, we tested the proposed approach on in-vivo data using 270 data sets, which were qualitatively validated by a neuroradiology expert. The results were quantitatively validated using the three commonly used metrics for segmentation evaluation: the Dice coefficient, the modified Hausdorff distance, and the absolute volume difference. The proposed approach showed a higher accuracy compared to two of the existing segmentation approaches.

GPU-Accelerated Simulation of Massive Spatial Data Based on the Modified Planar Rotator Model

A novel Gibbs Markov random field for spatial data on Cartesian grids based on the modified planar rotator (MPR) model of statistical physics has been recently introduced for efficient and automatic interpolation of big data sets, such as satellite and radar images. The MPR model does not rely on Gaussian assumptions. Spatial correlations are captured via nearest-neighbor interactions between transformed variables. This allows vectorization of the model which, along with an efficient hybrid Monte Carlo algorithm, leads to fast execution times that scale approximately linearly with system size. The present study takes advantage of the short-range nature of the interactions between the MPR variables to parallelize the algorithm on graphics processing units (GPUs) in the Compute Unified Device Architecture programming environment. It is shown that, for the processors employed, the GPU implementation can lead to impressive computational speedups, up to almost 500 times on large grids, compared to single-processor calculations. Consequently, massive data sets comprising millions of data points can be automatically processed in less than one second on an ordinary GPU.

Novel stochastic framework for automatic segmentation of human thigh MRI volumes and its applications in spinal cord injured individuals.

Severe spinal cord injury (SCI) leads to skeletal muscle atrophy and adipose tissue infiltration in the skeletal muscle, which can result in compromised muscle mechanical output and lead to health-related complications. In this study, we developed a novel automatic 3-D approach for volumetric segmentation and quantitative assessment of thigh Magnetic Resonance Imaging (MRI) volumes in individuals with chronic SCI as well as non-disabled individuals. In this framework, subcutaneous adipose tissue, inter-muscular adipose tissue and total muscle tissue are segmented using linear combination of discrete Gaussians algorithm. Also, three thigh muscle groups were segmented utilizing the proposed 3-D Joint Markov Gibbs Random Field model that integrates first order appearance model, spatial information, and shape model to localize the muscle groups. The accuracy of the automatic segmentation method was tested both on SCI (N = 16) and on non-disabled (N = 14) individuals, showing an overall 0.93±0.06 accuracy for adipose tissue and muscle compartments segmentation based on Dice Similarity Coefficient. The proposed framework for muscle compartment segmentation showed an overall higher accuracy compared to ANTs and STAPLE, two previously validated atlas-based segmentation methods. Also, the framework proposed in this study showed similar Dice accuracy and better Hausdorff distance measure to that obtained using DeepMedic Convolutional Neural Network structure, a well-known deep learning network for 3-D medical image segmentation. The automatic segmentation method proposed in this study can provide fast and accurate quantification of adipose and muscle tissues, which have important health and functional implications in the SCI population.

Gibbs Markov random fields with continuous values based on the modified planar rotator model

We introduce a novel Gibbs Markov random field for spatial data on Cartesian grids based on the modified planar rotator (MPR) model of statistical physics. The MPR captures spatial correlations using nearest-neighbor interactions of continuously-valued spins and does not rely on Gaussian assumptions. The only model parameter is the reduced temperature, which we estimate by means of an ergodic specific energy matching principle. We propose an efficient hybrid Monte Carlo simulation algorithm that leads to fast relaxation of the MPR model and allows vectorization. Consequently, the MPR computational time for inference and simulation scales approximately linearly with system size. This makes it more suitable for big data sets, such as satellite and radar images, than conventional geostatistical approaches. The performance (accuracy and computational speed) of the MPR model is validated with conditional simulation of Gaussian synthetic and non-Gaussian real data (atmospheric heat release measurements and Walker-lake DEM-based concentrations) and comparisons with standard gap-filling methods.

Early diabetic retinopathy diagnosis based on local retinal blood vessel analysis in optical coherence tomography angiography (OCTA) images.

This paper introduces a new computer-aided diagnosis (CAD) system for detecting early-stage diabetic retinopathy (DR) using optical coherence tomography angiography (OCTA) images. The proposed DR-CAD system is based on the analysis of new local features that describe both the appearance and retinal structure in OCTA images. It starts with a new segmentation approach that has the ability to extract the blood vessels from superficial and deep retinal OCTA maps. The high capability of our segmentation approach stems from using a joint Markov-Gibbs random field stochastic model integrating a 3D spatial statistical model with a first-order appearance model of the blood vessels. Following the segmentation step, three new local features are estimated from the segmented vessels and the foveal avascular zone (FAZ): (a) vessels density, (b) blood vessel calibre, and (c) width of the FAZ. To distinguish mild DR patients from normal cases, the estimated three features are used to train and test a support vector machine (SVM) classifier with the radial basis function (RBF) kernel. On a cohort of 105 subjects, the presented DR-CAD system demonstrated an overall accuracy (ACC) of 94.3%, a sensitivity of 97.9%, a specificity of 87.0%, the area under the curve (AUC) of 92.4%, and a Dice similarity coefficient (DSC) of 95.8%. This in turn demonstrates the promise of the proposed CAD system as a supplemental tool for early detection of DR. We developed a new DR-CAD system that is capable of diagnosing DR in its early stage. The proposed system is based on extracting three different features from the segmented OCTA images, which reflect the changes in the retinal vasculature network.

3D kidney segmentation from abdominal diffusion MRI using an appearance-guided deformable boundary.

A new technique for more accurate automatic segmentation of the kidney from its surrounding abdominal structures in diffusion-weighted magnetic resonance imaging (DW-MRI) is presented. This approach combines a new 3D probabilistic shape model of the kidney with a first-order appearance model and fourth-order spatial model of the diffusion-weighted signal intensity to guide the evolution of a 3D geometric deformable model. The probabilistic shape model was built from labeled training datasets to produce a spatially variant, independent random field of region labels. A Markov-Gibbs random field spatial model with up to fourth-order interactions was adequate to capture the inhomogeneity of renal tissues in the DW-MRI signal. A new analytical approach estimated the Gibbs potentials directly from the DW-MRI data to be segmented, in order that the segmentation procedure would be fully automatic. Finally, to better distinguish the kidney object from the surrounding tissues, marginal gray level distributions inside and outside of the deformable boundary were modeled with adaptive linear combinations of discrete Gaussians (first-order appearance model). The approach was tested on a cohort of 64 DW-MRI datasets with b-values ranging from 50 to 1000 s/mm2. The performance of the presented approach was evaluated using leave-one-subject-out cross validation and compared against three other well-known segmentation methods applied to the same DW-MRI data using the following evaluation metrics: 1) the Dice similarity coefficient (DSC); 2) the 95-percentile modified Hausdorff distance (MHD); and 3) the percentage kidney volume difference (PKVD). High performance of the new approach was confirmed by the high DSC (0.95±0.01), low MHD (3.9±0.76) mm, and low PKVD (9.5±2.2)% relative to manual segmentation by an MR expert (a board certified radiologist).