Dataset Information

Automatic quantitative analysis of pulmonary vascular morphology in CT images.

ABSTRACT:

Purpose

Vascular remodeling is a significant pathological feature of various pulmonary diseases, which may be assessed by quantitative computed tomography (CT) imaging. The purpose of this study was therefore to develop and validate an automatic method for quantifying pulmonary vascular morphology in CT images.

Methods

The proposed method consists of pulmonary vessel extraction and quantification. For extracting pulmonary vessels, a graph-cuts-based method is proposed which considers appearance (CT intensity) and shape (vesselness from a Hessian-based filter) features, and incorporates distance to the airways into the cost function to prevent false detection of airway walls. For quantifying the extracted pulmonary vessels, a radius histogram is generated by counting the occurrence of vessel radii, calculated from a distance transform-based method. Subsequently, two biomarkers, slope ? and intercept ?, are calculated by linear regression on the radius histogram. A public data set from the VESSEL12 challenge was used to independently evaluate the vessel extraction. The quantitative analysis method was validated using images of a three-dimensional (3D) printed vessel phantom, scanned by a clinical CT scanner and a micro-CT scanner (to obtain a gold standard). To confirm the association between imaging biomarkers and pulmonary function, 77 scleroderma patients were investigated with the proposed method.

Results

In the independent evaluation with the public data set, our vessel segmentation method obtained an area under the receiver operating characteristic (ROC) curve of 0.976. The median radius difference between clinical and micro-CT scans of a 3D printed vessel phantom was 0.062 ± 0.020 mm, with interquartile range of 0.199 ± 0.050 mm. In the studied patient group, a significant correlation between diffusion capacity for carbon monoxide and the biomarkers, ? (R = -0.27, P = 0.018) and ? (R = 0.321, P = 0.004), was obtained.

Conclusion

In conclusion, the proposed method was validated independently using a public data set resulting in an area under the ROC curve of 0.976 and using a 3D printed vessel phantom data set, showing a vessel sizing error of 0.062 mm (0.16 in-plane pixel units). The correlation between imaging biomarkers and diffusion capacity in a clinical data set confirmed an association between lung structure and function. This quantification of pulmonary vascular morphology may be helpful in understanding the pathophysiology of pulmonary vascular diseases.

SUBMITTER: Zhai Z

PROVIDER: S-EPMC6852650 | biostudies-literature | 2019 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Automatic quantitative analysis of pulmonary vascular morphology in CT images.

Zhai Zhiwei Z Staring Marius M Hernández Girón Irene I Veldkamp Wouter J H WJH Kroft Lucia J LJ Ninaber Maarten K MK Stoel Berend C BC

Medical physics 20190709 9

<h4>Purpose</h4>Vascular remodeling is a significant pathological feature of various pulmonary diseases, which may be assessed by quantitative computed tomography (CT) imaging. The purpose of this study was therefore to develop and validate an automatic method for quantifying pulmonary vascular morphology in CT images.<h4>Methods</h4>The proposed method consists of pulmonary vessel extraction and quantification. For extracting pulmonary vessels, a graph-cuts-based method is proposed which consid ...[more]

PMID: 31206181

Similar Datasets

Project description:X-ray computed tomography (CT) imaging of the brain is commonly used in diagnostic settings. Although CT scans are primarily used in clinical practice, they are increasingly used in research. A fundamental processing step in brain imaging research is brain extraction - the process of separating the brain tissue from all other tissues. Methods for brain extraction have either been 1) validated but not fully automated, or 2) fully automated and informally proposed, but never formally validated.To systematically analyze and validate the performance of FSL's brain extraction tool (BET) on head CT images of patients with intracranial hemorrhage. This was done by comparing the manual gold standard with the results of several versions of automatic brain extraction and by estimating the reliability of automated segmentation of longitudinal scans. The effects of the choice of BET parameters and data smoothing is studied and reported.All images were thresholded using a 0-100 Hounsfield unit (HU) range. In one variant of the pipeline, data were smoothed using a 3-dimensional Gaussian kernel (?=1mm(3)) and re-thresholded to 0-100HU; in the other, data were not smoothed. BET was applied using 1 of 3 fractional intensity (FI) thresholds: 0.01, 0.1, or 0.35 and any holes in the brain mask were filled. For validation against a manual segmentation, 36 images from patients with intracranial hemorrhage were selected from 19 different centers from the MISTIE (Minimally Invasive Surgery plus recombinant-tissue plasminogen activator for Intracerebral Evacuation) stroke trial. Intracranial masks of the brain were manually created by one expert CT reader. The resulting brain tissue masks were quantitatively compared to the manual segmentations using sensitivity, specificity, accuracy, and the Dice Similarity Index (DSI). Brain extraction performance across smoothing and FI thresholds was compared using the Wilcoxon signed-rank test. The intracranial volume (ICV) of each scan was estimated by multiplying the number of voxels in the brain mask by the dimensions of each voxel for that scan. From this, we calculated the ICV ratio comparing manual and automated segmentation: ICVautomated/ICVmanual. To estimate the performance in a large number of scans, brain masks were generated from the 6 BET pipelines for 1095 longitudinal scans from 129 patients. Failure rates were estimated from visual inspection. ICV of each scan was estimated and an intraclass correlation (ICC) was estimated using a one-way ANOVA.Smoothing images improves brain extraction results using BET for all measures except specificity (all p<0.01, uncorrected), irrespective of the FI threshold. Using an FI of 0.01 or 0.1 performed better than 0.35. Thus, all reported results refer only to smoothed data using an FI of 0.01 or 0.1. Using an FI of 0.01 had a higher median sensitivity (0.9901) than an FI of 0.1 (0.9884, median difference: 0.0014, p<0.001), accuracy (0.9971 vs. 0.9971; median difference: 0.0001, p<0.001), and DSI (0.9895 vs. 0.9894; median difference: 0.0004, p<0.001) and lower specificity (0.9981 vs. 0.9982; median difference: -0.0001, p<0.001). These measures are all very high indicating that a range of FI values may produce visually indistinguishable brain extractions. Using smoothed data and an FI of 0.01, the mean (SD) ICV ratio was 1.002 (0.008); the mean being close to 1 indicates the ICV estimates are similar for automated and manual segmentation. In the 1095 longitudinal scans, this pipeline had a low failure rate (5.2%) and the ICC estimate was high (0.929, 95% CI: 0.91, 0.945) for successfully extracted brains.BET performs well at brain extraction on thresholded, 1mm(3) smoothed CT images with an FI of 0.01 or 0.1. Smoothing before applying BET is an important step not previously discussed in the literature. Analysis code is provided.

Project description:PURPOSE:We evaluated the feasibility of using an automatic segmentation tool to delineate cardiac substructures from noncontrast computed tomography (CT) images for cardiac dosimetry and toxicity analyses for patients with nonsmall cell lung cancer (NSCLC) after radiotherapy. MATERIAL AND METHODS:We used an in-house developed multi-atlas segmentation tool to delineate 11cardiac substructures, including the whole heart, four heart chambers, and six greater vessels, automatically from the averaged 4D-CT planning images of 49 patients with NSCLC. Two experienced radiation oncologists edited the auto-segmented contours. Times for automatic segmentation and modification were recorded. The modified contours were compared with the auto-segmented contours in terms of Dice similarity coefficient (DSC) and mean surface distance (MSD) to evaluate the extent of modification. Differences in dose-volume histogram (DVH) characteristics were also evaluated for the modified versus auto-segmented contours. RESULTS:The mean automatic segmentation time for all 11 structures was 7-9?min. For the 49 patients, the mean DSC values (±SD) ranged from .73?±?.08 to .95?±?.04, and the mean MSD values ranged from 1.3?±?.6?mm to 2.9?±?5.1?mm. Overall, the modifications were small; the largest modifications were in the pulmonary vein and the inferior vena cava. The heart V30 (volume receiving dose ?30?Gy) and the mean dose to the whole heart and the four heart chambers were not different for the modified versus the auto-segmented contours based on the statistically significant condition of p?<?.05. Also, the maximum dose to the great vessels was no different except for the pulmonary vein. CONCLUSIONS:Automatic segmentation of cardiac substructures did not require substantial modifications. Dosimetric evaluation showed no significant difference between the auto-segmented and modified contours for most structures, which suggests that the auto-segmented contours can be used to study cardiac dose-responses in clinical practice.

Project description:BackgroundIn patients with pulmonary embolism (PE), there is a growing interest in quantifying the pulmonary vascular obtruction index (PVOI), which may be an independent risk factor for PE recurrence. Perfusion SPECT/CT is a very attractive tool to provide an accurate quantification of the PVOI. However, there is currently no reliable method to automatically delineate and quantify it. The aim of this phantom study was to assess and compare 3 segmentation methods for PVOI quantification with perfusion SPECT/CT imaging.MethodsThree hundred ninety-six SPECT/CT scans, with various PE scenarios (n = 44), anterior to posterior perfusion gradients (n = 3), and lung volumes (n = 3) were simulated using Simind software. Three segmentation methods were assesssed: (1) using an intensity threshold expressed as a percentage of the maximal voxel value (MaxTh), (2) using a Z-score threshold (ZTh) after building a Z-score parametric lung map, and (3) using a relative difference threshold (RelDiffTh) after building a relative difference parametric map. Ninety randomly selected simulations were used to define the optimal threshold, and 306 simulations were used for the complete analysis. Spacial correlation between PE volumes from the phantom data and the delineated PE volumes was assessed by computing DICEPE indices. Bland-Altman statistics were used to calculate agreement for PVOI between the phantom data and the segmentation methods.ResultsMean DICEPE index was higher with the RelDiffTh method (0.85 ± 0.08), as compared with the MaxTh method (0.78 ± 0.16) and the ZTh method (0.67 ± 0.15). Using the RelDiffTh method, mean DICEPE index remained high (> 0.81) regardless of the perfusion gradient and the lung volumes. Using the RelDiffTh method, mean relative difference in PVOI was - 12%, and the limits of agreement were - 40% to 16%. Values were 3% (- 75% to 81%) for MaxTh method and 0% (- 120% to 120%) for ZTh method. Graphycal analysis of the Bland-Altman graph for the RelDiffTh method showed very close estimation of the PVOI for small and medium PE, and a trend toward an underestimation of large PE.ConclusionIn this phantom study, a delineation method based on a relative difference parametric map provided a good estimation of the PVOI, regardless of the extent of PE, the intensity of the anterior to posterior gradient, and the whole lung volumes.

Dataset Information

Automatic quantitative analysis of pulmonary vascular morphology in CT images.

Purpose

Methods

Results

Conclusion

Publications

Automatic quantitative analysis of pulmonary vascular morphology in CT images.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets