3D Multi-branch Tubular Surface and Centerline Extraction with 4D Iterative Key Points

Novel Radiomic Measurements of Tumor-Associated Vasculature Morphology on Clinical Imaging as a Biomarker of Treatment Response in Multiple Cancers

PURPOSE

The tumor-associated vasculature (TAV) differs from healthy blood vessels by its convolutedness, leakiness, and chaotic architecture, and these attributes facilitate the creation of a treatment-resistant tumor microenvironment. Measurable differences in these attributes might also help stratify patients by likely benefit of systemic therapy (e.g., chemotherapy). In this work, we present a new category of computational image-based biomarkers called quantitative tumor-associated vasculature (QuanTAV) features, and demonstrate their ability to predict response and survival across multiple cancer types, imaging modalities, and treatment regimens involving chemotherapy.

EXPERIMENTAL DESIGN

We isolated tumor vasculature and extracted mathematical measurements of twistedness and organization from routine pretreatment radiology (CT or contrast-enhanced MRI) of a total of 558 patients, who received one of four first-line chemotherapy-based therapeutic intervention strategies for breast (n = 371) or non-small cell lung cancer (NSCLC, n = 187).

RESULTS

Across four chemotherapy-based treatment strategies, classifiers of QuanTAV measurements significantly (P < 0.05) predicted response in held out testing cohorts alone (AUC = 0.63-0.71) and increased AUC by 0.06-0.12 when added to models of significant clinical variables alone. Similarly, we derived QuanTAV risk scores that were prognostic of recurrence-free survival in treatment cohorts who received surgery following chemotherapy for breast cancer [P = 0.0022; HR = 1.25; 95% confidence interval (CI), 1.08-1.44; concordance index (C-index) = 0.66] and chemoradiation for NSCLC (P = 0.039; HR = 1.28; 95% CI, 1.01-1.62; C-index = 0.66). From vessel-based risk scores, we further derived categorical QuanTAV high/low risk groups that were independently prognostic among all treatment groups, including patients with NSCLC who received chemotherapy only (P = 0.034; HR = 2.29; 95% CI, 1.07-4.94; C-index = 0.62). QuanTAV response and risk scores were independent of clinicopathologic risk factors and matched or exceeded models of clinical variables including posttreatment response.

CONCLUSIONS

Across these domains, we observed an association of vascular morphology on CT and MRI-as captured by metrics of vessel curvature, torsion, and organizational heterogeneity-and treatment outcome. Our findings suggest the potential of shape and structure of the TAV in developing prognostic and predictive biomarkers for multiple cancers and different treatment strategies.

Trajectory Grouping With Curvature Regularization for Tubular Structure Tracking

Tubular structure tracking is a crucial task in the fields of computer vision and medical image analysis. The minimal paths-based approaches have exhibited their strong ability in tracing tubular structures, by which a tubular structure can be naturally modeled as a minimal geodesic path computed with a suitable geodesic metric. However, existing minimal paths-based tracing approaches still suffer from difficulties such as the shortcuts and short branches combination problems, especially when dealing with the images involving complicated tubular tree structures or background. In this paper, we introduce a new minimal paths-based model for minimally interactive tubular structure centerline extraction in conjunction with a perceptual grouping scheme. Basically, we take into account the prescribed tubular trajectories and curvature-penalized geodesic paths to seek suitable shortest paths. The proposed approach can benefit from the local smoothness prior on tubular structures and the global optimality of the used graph-based path searching scheme. Experimental results on both synthetic and real images prove that the proposed model indeed obtains outperformance comparing with the state-of-the-art minimal paths-based tubular structure tracing algorithms.

Minimal Paths for Tubular Structure Segmentation With Coherence Penalty and Adaptive Anisotropy

The minimal path method has proven to be particularly useful and efficient in tubular structure segmentation applications. In this paper, we propose a new minimal path model associated with a dynamic Riemannian metric embedded with an appearance feature coherence penalty and an adaptive anisotropy enhancement term. The features that characterize the appearance and anisotropy properties of a tubular structure are extracted through the associated orientation score. The proposed the dynamic Riemannian metric is updated in the course of the geodesic distance computation carried out by the efficient single-pass fast marching method. Compared to the state-of-the-art minimal path models, the proposed minimal path model is able to extract the desired tubular structures from a complicated vessel tree structure. In addition, we propose an efficient prior path-based method to search for vessel radius value at each centerline position of the target. Finally, we perform the numerical experiments on both synthetic and real images. The quantitive validation is carried out on retinal vessel images. The results indicate that the proposed model indeed achieves a promising performance.

Quantitative vessel tortuosity: A potential CT imaging biomarker for distinguishing lung granulomas from adenocarcinomas

Adenocarcinomas and active granulomas can both have a spiculated appearance on computed tomography (CT) and both are often fluorodeoxyglucose (FDG) avid on positron emission tomography (PET) scan, making them difficult to distinguish. Consequently, patients with benign granulomas are often subjected to invasive surgical biopsies or resections. In this study, quantitative vessel tortuosity (QVT), a novel CT imaging biomarker to distinguish between benign granulomas and adenocarcinomas on routine non-contrast lung CT scans is introduced. Our study comprised of CT scans of 290 patients from two different institutions, one cohort for training (N = 145) and the other (N = 145) for independent validation. In conjunction with a machine learning classifier, the top informative and stable QVT features yielded an area under receiver operating characteristic curve (ROC AUC) of 0.85 in the independent validation set. On the same cohort, the corresponding AUCs for two human experts including a radiologist and a pulmonologist were found to be 0.61 and 0.60, respectively. QVT features also outperformed well known shape and textural radiomic features which had a maximum AUC of 0.73 (p-value = 0.002), as well as features learned using a convolutional neural network AUC = 0.76 (p-value = 0.028). Our results suggest that QVT features could potentially serve as a non-invasive imaging biomarker to distinguish granulomas from adenocarcinomas on non-contrast CT scans.

Identification of dental root canals and their medial line from micro-CT and cone-beam CT records

Background

Shape of the dental root canal is highly patient specific. Automated identification methods of the medial line of dental root canals and the reproduction of their 3D shape can be beneficial for planning endodontic interventions as severely curved root canals or multi-rooted teeth may pose treatment challenges. Accurate shape information of the root canals may also be used by manufacturers of endodontic instruments in order to make more efficient clinical tools.

Method

Novel image processing procedures dedicated to the automated detection of the medial axis of the root canal from dental micro-CT and cone-beam CT records are developed. For micro-CT, the 3D model of the root canal is built up from several hundred parallel cross sections, using image enhancement, histogram based fuzzy c-means clustering, center point detection in the segmented slice, three dimensional inner surface reconstruction, and potential field driven curve skeleton extraction in three dimensions. Cone-beam CT records are processed with image enhancement filters and fuzzy chain based regional segmentation, followed by the reconstruction of the root canal surface and detecting its skeleton via a mesh contraction algorithm.

Results

The proposed medial line identification and root canal detection algorithms are validated on clinical data sets. 25 micro-CT and 36 cone-beam-CT records are used in the validation procedure. The overall success rate of the automatic dental root canal identification was about 92% in both procedures. The algorithms proved to be accurate enough for endodontic therapy planning.

Conclusions

Accurate medial line identification and shape detection algorithms of dental root canal have been developed. Different procedures are defined for micro-CT and cone-beam CT records. The automated execution of the subsequent processing steps allows easy application of the algorithms in the dental care. The output data of the image processing procedures is suitable for mathematical modeling of the central line. The proposed methods can help automate the preparation and design of several kinds of endodontic interventions.

Semi-automated segmentation of the sigmoid and descending colon for radiotherapy planning using the fast marching method

A fast and accurate segmentation of organs at risk, such as the healthy colon, would be of benefit for planning of radiotherapy, in particular in an adaptive scenario. For the treatment of pelvic tumours, a great challenge is the segmentation of the most adjacent and sensitive parts of the gastrointestinal tract, the sigmoid and descending colon. We propose a semi-automated method to segment these bowel parts using the fast marching (FM) method. Standard 3D computed tomography (CT) image data obtained from routine radiotherapy planning were used. Our pre-processing steps distinguish the intestine, muscles and air from connective tissue. The core part of our method separates the sigmoid and descending colon from the muscles and other segments of the intestine. This is done by utilizing the ability of the FM method to compute a specified minimal energy functional integrated along a path, and thereby extracting the colon centre line between user-defined control points in the sigmoid and descending colon. Further, we reconstruct the tube-shaped geometry of the sigmoid and descending colon by fitting ellipsoids to points on the path and by adding adjacent voxels that are likely voxels belonging to these bowel parts. Our results were compared to manually outlined sigmoid and descending colon, and evaluated using the Dice coefficient (DC). Tests on 11 patients gave an average DC of 0.83 (+/-0.07) with little user interaction. We conclude that the proposed method makes it possible to fast and accurately segment the sigmoid and descending colon from routine CT image data.