Segmentation of skin lesions in 2-D and 3-D ultrasound images using a spatially coherent generalized Rayleigh mixture model

Computer-Aided Diagnosis Methods for High-Frequency Ultrasound Data Analysis: A Review

Over the last few decades, computer-aided diagnosis systems have become a part of clinical practice. They have the potential to assist clinicians in daily diagnostic tasks. The image processing techniques are fast, repeatable, and robust, which helps physicians to detect, classify, segment, and measure various structures. The recent rapid development of computer methods for high-frequency ultrasound image analysis opens up new diagnostic paths in dermatology, allergology, cosmetology, and aesthetic medicine. This paper, being the first in this area, presents a research overview of high-frequency ultrasound image processing techniques, which have the potential to be a part of computer-aided diagnosis systems. The reviewed methods are categorized concerning the application, utilized ultrasound device, and image data-processing type. We present the bridge between diagnostic needs and already developed solutions and discuss their limitations and future directions in high-frequency ultrasound image analysis. A search was conducted of the technical literature from 2005 to September 2022, and in total, 31 studies describing image processing methods were reviewed. The quantitative and qualitative analysis included 39 algorithms, which were selected as the most effective in this field. They were completed by 20 medical papers and define the needs and opportunities for high-frequency ultrasound application and CAD development.

Automated segmentation of epidermis in high-frequency ultrasound of pathological skin using a cascade of DeepLab v3+ networks and fuzzy connectedness

This study proposes a novel, fully automated framework for epidermal layer segmentation in different skin diseases based on 75 MHz high-frequency ultrasound (HFUS) image data. A robust epidermis segmentation is a vital first step to detect changes in thickness, shape, and intensity and therefore support diagnosis and treatment monitoring in inflammatory and neoplastic skin lesions. Our framework links deep learning and fuzzy connectedness for image analysis. It consists of a cascade of two DeepLab v3+ models with a ResNet-50 backbone and a fuzzy connectedness analysis module for fine segmentation. Both deep models are pre-trained on the ImageNet dataset and subjected to transfer learning using our HFUS database of 580 images with atopic dermatitis, psoriasis and non-melanocytic skin tumors. The first deep model is used to detect the appropriate region of interest, while the second stands for the main segmentation procedure. We use the softmax layer of the latter twofold to prepare the input data for fuzzy connectedness analysis: as a reservoir of seed points and a direct contribution to the input image. In the experiments, we analyze different configurations of the framework, including region of interest detection, deep model backbones and training loss functions, or fuzzy connectedness analysis with parameter settings. We also use the Dice index and epidermis thickness to compare our results to state-of-the-art approaches. The Dice index of 0.919 yielded by our model over the entire dataset (and exceeding 0.93 in inflammatory diseases) proves its superiority over the other methods.

Statistical modeling of retinal optical coherence tomography using the Weibull mixture model

In this paper, a novel statistical model is proposed for retinal optical coherence tomography (OCT) images. According to the layered structure of the retina, a mixture of six Weibull distributions is proposed to describe the main statistical features of OCT images. We apply Weibull distribution to establish a more comprehensive model but with fewer parameters that has better goodness of fit (GoF) than previous models. Our new model also takes care of features such as asymmetry and heavy-tailed nature of the intensity distribution of retinal OCT data. In order to test the effectiveness of this new model, we apply it to improve the low quality of the OCT images. For this purpose, the spatially constrained Gaussian mixture model (SCGMM) is implemented. Since SCGMM is designed for data with Gaussian distribution, we convert our Weibull mixture model to a Gaussian mixture model using histogram matching before applying SCGMM. The denoising results illustrate the remarkable performance in terms of the contrast to noise ratio (CNR) and texture preservation (TP) compared to other peer methods. In another test to evaluate the efficiency of our proposed model, the parameters and GoF criteria are considered as a feature vector for support vector machine (SVM) to classify the healthy retinal OCT images from pigment epithelial detachment (PED) disease. The confusion matrix demonstrates the impact of the proposed model in our preliminary study on the OCT classification problem.

Deep Learning-Based High-Frequency Ultrasound Skin Image Classification with Multicriteria Model Evaluation

This study presents the first application of convolutional neural networks to high-frequency ultrasound skin image classification. This type of imaging opens up new opportunities in dermatology, showing inflammatory diseases such as atopic dermatitis, psoriasis, or skin lesions. We collected a database of 631 images with healthy skin and different skin pathologies to train and assess all stages of the methodology. The proposed framework starts with the segmentation of the epidermal layer using a DeepLab v3+ model with a pre-trained Xception backbone. We employ transfer learning to train the segmentation model for two purposes: to extract the region of interest for classification and to prepare the skin layer map for classification confidence estimation. For classification, we train five models in different input data modes and data augmentation setups. We also introduce a classification confidence level to evaluate the deep model’s reliability. The measure combines our skin layer map with the heatmap produced by the Grad-CAM technique designed to indicate image regions used by the deep model to make a classification decision. Moreover, we propose a multicriteria model evaluation measure to select the optimal model in terms of classification accuracy, confidence, and test dataset size. The experiments described in the paper show that the DenseNet-201 model fed with the extracted region of interest produces the most reliable and accurate results.

Automated seeding for ultrasound skin lesion segmentation

The segmentation of cancer-suspicious skin lesions using ultrasound may help their differential diagnosis and treatment planning. Active contour models (ACM) require an initial seed, which when manually chosen may cause variations in segmentation accuracy. Fully-automated skin segmentation typically employs layer-by-layer segmentation using a combination of methods; however, such segmentation has not yet been applied on cancerous lesions. In the current work, fully automated segmentation is achieved in two steps: an automated seeding (AS) step using a layer-by-layer method followed by a growing step using an ACM. The method was tested on images of nevi, melanomas, and basal cell carcinomas from two ultrasound imaging systems (N=60), with all lesions being successfully located. For the seeding step, manual seeding (MS) was used as a reference. AS approached the accuracy of MS when the latter used an optimal bounding rectangle based on the ground truth (Sørensen-Dice coefficient (SDC) of 72.3 vs 74.6, respectively). The effect of varying the manual seed was also investigated; a 0.7 decrease in seed height and width caused a mean SDC of 54.6. The results show the robustness of automated seeding for skin lesion segmentation.

Fast and accurate segmentation method of active shape model with Rayleigh mixture model clustering for prostate ultrasound images

BACKGROUND AND OBJECTIVE

The prostate cancer interventions, which need an accurate prostate segmentation, are performed under ultrasound imaging guidance. However, prostate ultrasound segmentation is facing two challenges. The first is the low signal-to-noise ratio and inhomogeneity of the ultrasound image. The second is the non-standardized shape and size of the prostate.

METHODS

For prostate ultrasound image segmentation, this paper proposed an accurate and efficient method of Active shape model (ASM) with Rayleigh mixture model Clustering (ASM-RMMC). Firstly, Rayleigh mixture model (RMM) is adopted for clustering the image regions which present similar speckle distributions. These content-based clustered images are then used to initialize and guide the deformation of an ASM model.

RESULTS

The performance of the proposed method is assessed on 30 prostate ultrasound images using four metrics as Mean Average Distance (MAD), Dice Similarity Coefficient (DSC), False Positive Error (FPE) and False Negative Error (FNE). The proposed ASM-RMMC reaches high segmentation accuracy with 95% ± 0.81% for DSC, 1.86 ± 0.02 pixels for MAD, 2.10% ± 0.36% for FPE and 2.78% ± 0.71% for FNE, respectively. Moreover, the average segmentation time is less than 8 s when treating a single prostate ultrasound image through ASM-RMMC.

CONCLUSIONS

This paper presents a method for prostate ultrasound image segmentation, which achieves high accuracy with less computational complexity and meets the clinical requirements.

Multi-Grid Phase Field Skin Tumor Segmentation in 3D Ultrasound Images

The aim of this paper is to present a new method for skin tumor segmentation in the 3D ultrasound images. We consider a variational formulation, the energy of which combines a diffuse interface phase field model (regularization term) and a log-likelihood computed using nonparametric estimates (data attachment term). We propose a multi-grid implementation with the exact solutions which has the advantage to avoid space discretization and numerical instabilities. The resulting algorithm is simple and easy to implement in multi-dimensions. Concerning applications, we focus on skin tumor segmentation. The clinical dataset used for the experiments is composed of 12 images with the ground truth given by a dermatologist. Comparisons with the reference methods show that the proposed method is more robust to the choice of the volume initialization. Moreover, thanks to the flexibility introduced by the diffuse interface, the sensitivity increases by 12% if the initialization is inside the lesion, and the Dice index increases by 59%, if the initialization covers the entire lesion. These results show that this new method is well designed to tackle the problem of underestimation of tumor volumes.

An Efficient Melanoma Diagnosis Approach Using Integrated HMF Multi-Atlas Map Based Segmentation

Melanoma is a life threading disease when it grows outside the corium layer of the skin. Mortality rates of the Melanoma cases are maximum among the skin cancer patients. The cost required for the treatment of advanced melanoma cases is very high and the survival rate is low. Numerous computerized dermoscopy systems are developed based on the combination of shape, texture and color features to facilitate early diagnosis of melanoma. The availability and cost of the dermoscopic imaging system is still an issue. To mitigate this issue, this paper presented an integrated segmentation and Third Dimensional (3D) feature extraction approach for the accurate diagnosis of melanoma. A multi-atlas method is applied for the image segmentation. The patch-based label fusion model is expressed in a Bayesian framework to improve the segmentation accuracy. A depth map is obtained from the Two-dimensional (2D) dermoscopic image for reconstructing the 3D skin lesion represented as structure tensors. The 3D shape features including the relative depth features are obtained. Streaks are the significant morphological terms of the melanoma in the radial growth phase. The proposed method yields maximum segmentation accuracy, sensibility, specificity and minimum cost function than the existing segmentation technique and classifier.

Generalized Bayesian Model Selection for Speckle on Remote Sensing Images

Synthetic aperture radar (SAR) and ultrasound (US) are two important active imaging techniques for remote sensing, both of which are subject to speckle noise caused by coherent summation of back-scattered waves and subsequent nonlinear envelope transformations. Estimating the characteristics of this multiplicative noise is crucial to develop denoising methods and to improve statistical inference from remote sensing images. In this paper, reversible jump Markov chain Monte Carlo (RJMCMC) algorithm has been used with a wider interpretation and a recently proposed RJMCMC-based Bayesian approach, trans-space RJMCMC, has been utilized. The proposed method provides an automatic model class selection mechanism for remote sensing images of SAR and US where the model class space consists of popular envelope distribution families. The proposed method estimates the correct distribution family, as well as the shape and the scale parameters, avoiding performing an exhaustive search. For the experimental analysis, different SAR images of urban, forest and agricultural scenes, and two different US images of a human heart have been used. Simulation results show the efficiency of the proposed method in finding statistical models for speckle.

Semi-Automatic Algorithms for Estimation and Tracking of AP-Diameter of the IVC in Ultrasound Images

Acutely ill patients presenting with conditions such as sepsis, trauma, and congestive heart failure require judicious resuscitation in order to achieve and maintain optimal circulating blood volume. Increasingly, emergency and critical care physicians are using portable ultrasound to approximate the temporal changes of the anterior–posterior (AP)-diameter of the inferior vena cava (IVC) in order to guide fluid administration or removal. This paper proposes semi-automatic active ellipse and rectangle algorithms capable of improved and quantified measurement of the AP-diameter. The proposed algorithms are compared to manual measurement and a previously published active circle model. Results demonstrate that the rectangle model outperforms both active circle and ellipse irrespective of IVC shape and closely approximates tedious expert assessment.

Application of automatic statistical post-processing method for analysis of ultrasonic and digital dermatoscopy images

Ultrasonic and digital dermatoscopy diagnostic methods are used in order to estimate the changes of structure, as well as to non-invasively measure the changes of parameters of lesions of human tissue. These days, it is very actual to perform the quantitative analysis of medical data, which allows to achieve the reliable early-stage diagnosis of lesions and help to save more lives. The proposed automatic statistical post-processing method based on integration of ultrasonic and digital dermatoscopy measurements is intended to estimate the parameters of malignant tumours, measure spatial dimensions (e.g. thickness) and shape, and perform faster diagnostics by increasing the accuracy of tumours differentiation. It leads to optimization of time-consuming analysis procedures of medical images and could be used as a reliable decision support tool in the field of dermatology.

Estimation and tracking of AP-diameter of the inferior vena cava in ultrasound images using a novel active circle algorithm

Medical research suggests that the anterior-posterior (AP)-diameter of the inferior vena cava (IVC) and its associated temporal variation as imaged by bedside ultrasound is useful in guiding fluid resuscitation of the critically-ill patient. Unfortunately, indistinct edges and gaps in vessel walls are frequently present which impede accurate estimation of the IVC AP-diameter for both human operators and segmentation algorithms. The majority of research involving use of the IVC to guide fluid resuscitation involves manual measurement of the maximum and minimum AP-diameter as it varies over time. This effort proposes using a time-varying circle fitted inside the typically ellipsoid IVC as an efficient, consistent and novel approach to tracking and approximating the AP-diameter even in the context of poor image quality. In this active-circle algorithm, a novel evolution functional is proposed and shown to be a useful tool for ultrasound image processing. The proposed algorithm is compared with an expert manual measurement, and state-of-the-art relevant algorithms. It is shown that the algorithm outperforms other techniques and performs very close to manual measurement.

Fast Unsupervised Bayesian Image Segmentation With Adaptive Spatial Regularisation

This paper presents a new Bayesian estimation technique for hidden Potts-Markov random fields with unknown regularisation parameters, with application to fast unsupervised K -class image segmentation. The technique is derived by first removing the regularisation parameter from the Bayesian model by marginalisation, followed by a small-variance-asymptotic (SVA) analysis in which the spatial regularisation and the integer-constrained terms of the Potts model are decoupled. The evaluation of this SVA Bayesian estimator is then relaxed into a problem that can be computed efficiently by iteratively solving a convex total-variation denoising problem and a least-squares clustering ( K -means) problem, both of which can be solved straightforwardly, even in high-dimensions, and with parallel computing techniques. This leads to a fast fully unsupervised Bayesian image segmentation methodology in which the strength of the spatial regularisation is adapted automatically to the observed image during the inference procedure, and that can be easily applied in large 2D and 3D scenarios or in applications requiring low computing times. Experimental results on synthetic and real images, as well as extensive comparisons with state-of-the-art algorithms, confirm that the proposed methodology offer extremely fast convergence and produces accurate segmentation results, with the important additional advantage of self-adjusting regularisation parameters.

Semi-automatic bowel wall thickness measurements on MR enterography in patients with Crohn's disease

OBJECTIVE

To evaluate a semi-automatic method for delineation of the bowel wall and measurement of the wall thickness in patients with Crohn's disease.

METHODS

53 patients with suspected or proven Crohn's disease were selected. Two radiologists independently supervised the delineation of regions with active Crohn's disease on MRI, yielding manual annotations (Ano1, Ano2). Three observers manually measured the maximal bowel wall thickness of each annotated segment. An active contour segmentation approach semi-automatically delineated the bowel wall. For each active region, two segmentations (Seg1, Seg2) were obtained by independent observers, in which the maximum wall thickness was automatically determined. The overlap between (Seg1, Seg2) was compared with the overlap of (Ano1, Ano2) using Wilcoxon's signed rank test. The corresponding variances were compared using the Brown-Forsythe test. The variance of the semi-automatic thickness measurements was compared with the overall variance of manual measurements through an F-test. Furthermore, the intraclass correlation coefficient (ICC) of semi-automatic thickness measurements was compared with the ICC of manual measurements through a likelihood-ratio test.

RESULTS

Patient demographics: median age, 30 years; interquartile range, 25-38 years; 33 females. The median overlap of the semi-automatic segmentations (Seg1 vs Seg2: 0.89) was significantly larger than the median overlap of the manual annotations (Ano1 vs Ano2: 0.72); p = 1.4 × 10^-5. The variance in overlap of the semi-automatic segmentations was significantly smaller than the variance in overlap of the manual annotations (p = 1.1 × 10^-9). The variance of the semi-automated measurements (0.46 mm²) was significantly smaller than the variance of the manual measurements (2.90 mm², p = 1.1 × 10^-7). The ICC of semi-automatic measurement (0.88) was significantly higher than the ICC of manual measurement (0.45); p = 0.005.

CONCLUSION

The semi-automatic technique facilitates reproducible delineation of regions with active Crohn's disease. The semi-automatic thickness measurement sustains significantly improved interobserver agreement. Advances in knowledge: Automation of bowel wall thickness measurements strongly increases reproducibility of these measurements, which are commonly used in MRI scoring systems of Crohn's disease activity.

Melanoma Is Skin Deep: A 3D Reconstruction Technique for Computerized Dermoscopic Skin Lesion Classification

Melanoma mortality rates are the highest amongst skin cancer patients. Melanoma is life threating when it grows beyond the dermis of the skin. Hence, depth is an important factor to diagnose melanoma. This paper introduces a non-invasive computerized dermoscopy system that considers the estimated depth of skin lesions for diagnosis. A 3-D skin lesion reconstruction technique using the estimated depth obtained from regular dermoscopic images is presented. On basis of the 3-D reconstruction, depth and 3-D shape features are extracted. In addition to 3-D features, regular color, texture, and 2-D shape features are also extracted. Feature extraction is critical to achieve accurate results. Apart from melanoma, in-situ melanoma the proposed system is designed to diagnose basal cell carcinoma, blue nevus, dermatofibroma, haemangioma, seborrhoeic keratosis, and normal mole lesions. For experimental evaluations, the PH2, ISIC: Melanoma Project, and ATLAS dermoscopy data sets is considered. Different feature set combinations is considered and performance is evaluated. Significant performance improvement is reported the post inclusion of estimated depth and 3-D features. The good classification scores of sensitivity = 96%, specificity = 97% on PH2 data set and sensitivity = 98%, specificity = 99% on the ATLAS data set is achieved. Experiments conducted to estimate tumor depth from 3-D lesion reconstruction is presented. Experimental results achieved prove that the proposed computerized dermoscopy system is efficient and can be used to diagnose varied skin lesion dermoscopy images.

Segmentation of Skin Tumors in High-Frequency 3-D Ultrasound Images

High-frequency 3-D ultrasound imaging is an informative tool for diagnosis, surgery planning and skin lesion examination. The purpose of this article was to describe a semi-automated segmentation tool providing easy access to the extent, shape and volume of a lesion. We propose an adaptive log-likelihood level-set segmentation procedure using non-parametric estimates of the intensity distribution. The algorithm has a single parameter to control the smoothness of the contour, and we describe how a fixed value yields satisfactory segmentation results with an average Dice coefficient of D = 0.76. The algorithm is implemented on a grid, which increases the speed by a factor of 100 compared with a standard pixelwise segmentation. We compare the method with parametric methods making the hypothesis of Rayleigh or Nakagami distributed signals, and illustrate that our method has greater robustness with similar computational speed. Benchmarks are made on realistic synthetic ultrasound images and a data set of nine clinical 3-D images acquired with a 50-MHz imaging system. The proposed algorithm is suitable for use in a clinical context as a post-processing tool.

Local Transverse-Slice-Based Level-Set Method for Segmentation of 3-D High-Frequency Ultrasonic Backscatter From Dissected Human Lymph Nodes

OBJECTIVE

To detect metastases in freshly excised human lymph nodes (LNs) using three-dimensional (3-D), high-frequency, quantitative ultrasound (QUS) methods, the LN parenchyma (LNP) must be segmented to preclude QUS analysis of data in regions outside the LNP and to compensate ultrasound attenuation effects due to overlying layers of LNP and residual perinodal fat (PNF).

METHODS

After restoring the saturated radio-frequency signals from PNF using an approach based on smoothing cubic splines, the three regions, i.e., LNP, PNF, and normal saline (NS), in the LN envelope data are segmented using a new, automatic, 3-D, three-phase, statistical transverseslice-based level-set (STS-LS) method that amends Lankton's method. Due to ultrasound attenuation and focusing effects, the speckle statistics of the envelope data vary with imaged depth. Thus, to mitigate depth-related inhomogeneity effects, the STS-LS method employs gamma probabilitydensity functions to locally model the speckle statistics within consecutive transverse slices.

RESULTS

Accurate results were obtained on simulated data. On a representative dataset of 54 LNs acquired from colorectal-cancer patients, the Dice similarity coefficient for LNP, PNF, and NS were 0.938 ± 0.025, 0.832 ± 0.086, and 0.968 ± 0.008, respectively, when compared to expert manual segmentation.

CONCLUSION

The STS-LS outperforms the established methods based on global and local statistics in our datasets and is capable of accurately handling the depth-dependent effects due to attenuation and focusing.

SIGNIFICANCE

This advance permits the automatic QUS-based cancer detection in the LNs. Furthermore, the STS-LS method could potentially be used in a wide range of ultrasound-imaging applications suffering from depth-dependent effects.

Joint Segmentation and Deconvolution of Ultrasound Images Using a Hierarchical Bayesian Model Based on Generalized Gaussian Priors

This paper proposes a joint segmentation and deconvolution Bayesian method for medical ultrasound (US) images. Contrary to piecewise homogeneous images, US images exhibit heavy characteristic speckle patterns correlated with the tissue structures. The generalized Gaussian distribution (GGD) has been shown to be one of the most relevant distributions for characterizing the speckle in US images. Thus, we propose a GGD-Potts model defined by a label map coupling US image segmentation and deconvolution. The Bayesian estimators of the unknown model parameters, including the US image, the label map, and all the hyperparameters are difficult to be expressed in a closed form. Thus, we investigate a Gibbs sampler to generate samples distributed according to the posterior of interest. These generated samples are finally used to compute the Bayesian estimators of the unknown parameters. The performance of the proposed Bayesian model is compared with the existing approaches via several experiments conducted on realistic synthetic data and in vivo US images.

ULTRASOUND SPECTROSCOPY

We introduce the concept of "Ultrasound Spectroscopy". The premise of ultrasound spectroscopy is that by acquiring ultrasound RF data at multiple power and frequency settings, a rich set of features can be extracted from that RF data and used to characterize the underlying tissues. This is beneficial for a variety of problems, such as accurate tissue classification, application-specific image generation, and numerous other quantitative tasks. These capabilities are particularly relevant to point-of-care ultrasound (POCUS) applications, where operator experience with ultrasound may be limited. Instead of displaying B-mode images, a POCUS application using ultrasound spectroscopy can, for example, automatically detect internal abdominal bleeding. In this paper, we present ex vivo tissue phantom studies to demonstrate the accuracy of ultrasound spectroscopy over previous approaches. Our studies suggest that ultrasound spectroscopy provides exceptional accuracy and informative features for classifying blood versus other tissues across image locations and body habitus.

Active learning based segmentation of Crohns disease from abdominal MRI

This paper proposes a novel active learning (AL) framework, and combines it with semi supervised learning (SSL) for segmenting Crohns disease (CD) tissues from abdominal magnetic resonance (MR) images. Robust fully supervised learning (FSL) based classifiers require lots of labeled data of different disease severities. Obtaining such data is time consuming and requires considerable expertise. SSL methods use a few labeled samples, and leverage the information from many unlabeled samples to train an accurate classifier. AL queries labels of most informative samples and maximizes gain from the labeling effort. Our primary contribution is in designing a query strategy that combines novel context information with classification uncertainty and feature similarity. Combining SSL and AL gives a robust segmentation method that: (1) optimally uses few labeled samples and many unlabeled samples; and (2) requires lower training time. Experimental results show our method achieves higher segmentation accuracy than FSL methods with fewer samples and reduced training effort.

Motion estimation-based image enhancement in ultrasound imaging

High resolution medical ultrasound (US) imaging is an ongoing challenge in many diagnosis applications and can be achieved either by instrumentation or by post-processing. Though many works have considered the issue of resolution enhancement in optical imaging, very few works have investigated this issue in US imaging. In optics, several algorithms have been proposed to achieve super-resolution (SR) image reconstruction, which consists of merging several low resolution images to create a higher resolution image. However, the straightforward implementation of such techniques for US imaging is unsuccessful, due to the interaction of ultrasound with tissue and speckle. We show how to overcome the limit of SR in this framework by refining the registration part of common multiframe techniques. For this purpose, we investigate motion estimation methods adapted to US imaging. Performance of the proposed technique is evaluated on both realistic simulated US images (providing an estimated best-case performance) and real US sequences of phantom and in-vivo thyroid images. Compared to classical SR methods, our technique brings both quantitative and qualitative improvements. Resolution gain was found to be 1.41 for the phantom sequence and 1.12 for the thyroid sequence and a quantitative study using the phantom further confirmed the spatial resolution enhancement. Furthermore, the contrast-to-noise ratio was increased by 27% and 13% for simulated and experimental US images, respectively.

Curvelet based automatic segmentation of supraspinatus tendon from ultrasound image: a focused assistive diagnostic method

Background

Disorders of rotator cuff tendons results in acute pain limiting the normal range of motion for shoulder. Of all the tendons in rotator cuff, supraspinatus (SSP) tendon is affected first of any pathological changes. Diagnosis of SSP tendon using ultrasound is considered to be operator dependent with its accuracy being related to operator’s level of experience.

Methods

The automatic segmentation of SSP tendon ultrasound image was performed to provide focused and more accurate diagnosis. The image processing techniques were employed for automatic segmentation of SSP tendon. The image processing techniques combines curvelet transform and mathematical concepts of logical and morphological operators along with area filtering. The segmentation assessment was performed using true positives rate, false positives rate and also accuracy of segmentation. The specificity and sensitivity of the algorithm was tested for diagnosis of partial thickness tears (PTTs) and full thickness tears (FTTs). The ultrasound images of SSP tendon were taken from medical center with the help of experienced radiologists. The algorithm was tested on 116 images taken from 51 different patients.

Results

The accuracy of segmentation of SSP tendon was calculated to be 95.61% in accordance with the segmentation performed by radiologists, with true positives rate of 91.37% and false positives rate of 8.62%. The specificity and sensitivity was found to be 93.6%, 94% and 95%, 95.6% for partial thickness tears and full thickness tears respectively. The proposed methodology was successfully tested over a database of more than 116 US images, for which radiologist assessment and validation was performed.

Conclusions

The segmentation of SSP tendon from ultrasound images helps in focused, accurate and more reliable diagnosis which has been verified with the help of two experienced radiologists. The specificity and sensitivity for accurate detection of partial and full thickness tears has been considerably increased after segmentation when compared with existing literature.

Automatic Detection and Segmentation of Crohn's Disease Tissues From Abdominal MRI

We propose an information processing pipeline for segmenting parts of the bowel in abdominal magnetic resonance images that are affected with Crohn's disease. Given a magnetic resonance imaging test volume, it is first oversegmented into supervoxels and each supervoxel is analyzed to detect presence of Crohn's disease using random forest (RF) classifiers. The supervoxels identified as containing diseased tissues define the volume of interest (VOI). All voxels within the VOI are further investigated to segment the diseased region. Probability maps are generated for each voxel using a second set of RF classifiers which give the probabilities of each voxel being diseased, normal or background. The negative log-likelihood of these maps are used as penalty costs in a graph cut segmentation framework. Low level features like intensity statistics, texture anisotropy and curvature asymmetry, and high level context features are used at different stages. Smoothness constraints are imposed based on semantic information (importance of each feature to the classification task) derived from the second set of learned RF classifiers. Experimental results show that our method achieves high segmentation accuracy with Dice metric values of 0.90 ± 0.04 and Hausdorff distance of 7.3 ± 0.8 mm. Semantic information and context features are an integral part of our method and are robust to different levels of added noise.

Estimating the granularity coefficient of a Potts-Markov random field within a Markov chain Monte Carlo algorithm

This paper addresses the problem of estimating the Potts parameter β jointly with the unknown parameters of a Bayesian model within a Markov chain Monte Carlo (MCMC) algorithm. Standard MCMC methods cannot be applied to this problem because performing inference on β requires computing the intractable normalizing constant of the Potts model. In the proposed MCMC method, the estimation of β is conducted using a likelihood-free Metropolis-Hastings algorithm. Experimental results obtained for synthetic data show that estimating β jointly with the other unknown parameters leads to estimation results that are as good as those obtained with the actual value of β. On the other hand, choosing an incorrect value of β can degrade estimation performance significantly. To illustrate the interest of this method, the proposed algorithm is successfully applied to real bidimensional SAR and tridimensional ultrasound images.

Segmentation of skin lesions in 2-D and 3-D ultrasound images using a spatially coherent generalized Rayleigh mixture model

This paper addresses the problem of jointly estimating the statistical distribution and segmenting lesions in multiple-tissue high-frequency skin ultrasound images. The distribution of multiple-tissue images is modeled as a spatially coherent finite mixture of heavy-tailed Rayleigh distributions. Spatial coherence inherent to biological tissues is modeled by enforcing local dependence between the mixture components. An original Bayesian algorithm combined with a Markov chain Monte Carlo method is then proposed to jointly estimate the mixture parameters and a label-vector associating each voxel to a tissue. More precisely, a hybrid Metropolis-within-Gibbs sampler is used to draw samples that are asymptotically distributed according to the posterior distribution of the Bayesian model. The Bayesian estimators of the model parameters are then computed from the generated samples. Simulation results are conducted on synthetic data to illustrate the performance of the proposed estimation strategy. The method is then successfully applied to the segmentation of in vivo skin tumors in high-frequency 2-D and 3-D ultrasound images.