A cubic B-spline-based hybrid registration of lung CT images for a dynamic airway geometric model with large deformation

An automated methodology for whole-body, multimodality tracking of individual cancer lesions

Objective. Manual analysis of individual cancer lesions to assess disease response is clinically impractical and requires automated lesion tracking methodologies. However, no methodology has been developed for whole-body individual lesion tracking, across an arbitrary number of scans, and acquired with various imaging modalities.Approach. This study introduces a lesion tracking methodology and benchmarked it using 2368Ga-DOTATATE PET/CT and PET/MR images of eight neuroendocrine tumor patients. The methodology consists of six steps: (1) alignment of multiple scans via image registration, (2) body-part labeling, (3) automatic lesion-wise dilation, (4) clustering of lesions based on local lesion shape metrics, (5) assignment of lesion tracks, and (6) output of a lesion graph. Registration performance was evaluated via landmark distance, lesion matching accuracy was evaluated between each image pair, and lesion tracking accuracy was evaluated via identical track ratio. Sensitivity studies were performed to evaluate the impact of lesion dilation (fixed versus automatic dilation), anatomic location, image modalities (inter- versus intra-modality), registration mode (direct versus indirect registration), and track size (number of time-points and lesions) on lesion matching and tracking performance.Main results. Manual contouring yielded 956 lesions, 1570 lesion-matching decisions, and 493 lesion tracks. The median residual registration error was 2.5 mm. The automatic lesion dilation led to 0.90 overall lesion matching accuracy, and an 88% identical track ratio. The methodology is robust regarding anatomic locations, image modalities, and registration modes. The number of scans had a moderate negative impact on the identical track ratio (94% for 2 scans, 91% for 3 scans, and 81% for 4 scans). The number of lesions substantially impacted the identical track ratio (93% for 2 nodes versus 54% for ≥5 nodes).Significance. The developed methodology resulted in high lesion-matching accuracy and enables automated lesion tracking in PET/CT and PET/MR.

Deep learning-based lung image registration: A review

Lung image registration can effectively describe the relative motion of lung tissues, thereby helping to solve series problems in clinical applications. Since the lungs are soft and fairly passive organs, they are influenced by respiration and heartbeat, resulting in discontinuity of lung motion and large deformation of anatomic features. This poses great challenges for accurate registration of lung image and its applications. The recent application of deep learning (DL) methods in the field of medical image registration has brought promising results. However, a versatile registration framework has not yet emerged due to diverse challenges of registration for different regions of interest (ROI). DL-based image registration methods used for other ROI cannot achieve satisfactory results in lungs. In addition, there are few review articles available on DL-based lung image registration. In this review, the development of conventional methods for lung image registration is briefly described and a more comprehensive survey of DL-based methods for lung image registration is illustrated. The DL-based methods are classified according to different supervision types, including fully-supervised, weakly-supervised and unsupervised. The contributions of researchers in addressing various challenges are described, as well as the limitations of these approaches. This review also presents a comprehensive statistical analysis of the cited papers in terms of evaluation metrics and loss functions. In addition, publicly available datasets for lung image registration are also summarized. Finally, the remaining challenges and potential trends in DL-based lung image registration are discussed.

Automatic landmark detection and registration of brain cortical surfaces via quasi-conformal geometry and convolutional neural networks

In medical imaging, surface registration is extensively used for performing systematic comparisons between anatomical structures, with a prime example being the highly convoluted brain cortical surfaces. To obtain a meaningful registration, a common approach is to identify prominent features on the surfaces and establish a low-distortion mapping between them with the feature correspondence encoded as landmark constraints. Prior registration works have primarily focused on using manually labeled landmarks and solving highly nonlinear optimization problems, which are time-consuming and hence hinder practical applications. In this work, we propose a novel framework for the automatic landmark detection and registration of brain cortical surfaces using quasi-conformal geometry and convolutional neural networks. We first develop a landmark detection network (LD-Net) that allows for the automatic extraction of landmark curves given two prescribed starting and ending points based on the surface geometry. We then utilize the detected landmarks and quasi-conformal theory for achieving the surface registration. Specifically, we develop a coefficient prediction network (CP-Net) for predicting the Beltrami coefficients associated with the desired landmark-based registration and a mapping network called the disk Beltrami solver network (DBS-Net) for generating quasi-conformal mappings from the predicted Beltrami coefficients, with the bijectivity guaranteed by quasi-conformal theory. Experimental results are presented to demonstrate the effectiveness of our proposed framework. Altogether, our work paves a new way for surface-based morphometry and medical shape analysis.

DAART: a deep learning platform for deeply accelerated adaptive radiation therapy for lung cancer

Purpose

The study aimed to implement a novel, deeply accelerated adaptive radiation therapy (DAART) approach for lung cancer radiotherapy (RT). Lung cancer is the most common cause of cancer-related death, and RT is the preferred medically inoperable treatment for early stage non-small cell lung cancer (NSCLC). In the current lengthy workflow, it takes a median of four weeks from diagnosis to RT treatment, which can result in complete restaging and loss of local control with delay. We implemented the DAART approach, featuring a novel deepPERFECT system, to address unwanted delays between diagnosis and treatment initiation.

Materials and methods

We developed a deepPERFECT to adapt the initial diagnostic imaging to the treatment setup to allow initial RT planning and verification. We used data from 15 patients with NSCLC treated with RT to train the model and test its performance. We conducted a virtual clinical trial to evaluate the treatment quality of the proposed DAART for lung cancer radiotherapy.

Results

We found that deepPERFECT predicts planning CT with a mean high-intensity fidelity of 83 and 14 HU for the body and lungs, respectively. The shape of the body and lungs on the synthesized CT was highly conformal, with a dice similarity coefficient (DSC) of 0.91, 0.97, and Hausdorff distance (HD) of 7.9 mm, and 4.9 mm, respectively, compared with the planning CT scan. The tumor showed less conformality, which warrants acquisition of treatment Day1 CT and online adaptive RT. An initial plan was designed on synthesized CT and then adapted to treatment Day1 CT using the adapt to position (ATP) and adapt to shape (ATS) method. Non-inferior plan quality was achieved by the ATP scenario, while all ATS-adapted plans showed good plan quality.

Conclusion

DAART reduces the common online ART (ART) treatment course by at least two weeks, resulting in a 50% shorter time to treatment to lower the chance of restaging and loss of local control.

deepPERFECT: Novel Deep Learning CT Synthesis Method for Expeditious Pancreatic Cancer Radiotherapy

Major sources of delay in the standard of care RT workflow are the need for multiple appointments and separate image acquisition. In this work, we addressed the question of how we can expedite the workflow by synthesizing planning CT from diagnostic CT. This idea is based on the theory that diagnostic CT can be used for RT planning, but in practice, due to the differences in patient setup and acquisition techniques, separate planning CT is required. We developed a generative deep learning model, deepPERFECT, that is trained to capture these differences and generate deformation vector fields to transform diagnostic CT into preliminary planning CT. We performed detailed analysis both from an image quality and a dosimetric point of view, and showed that deepPERFECT enabled the preliminary RT planning to be used for preliminary and early plan dosimetric assessment and evaluation.

Assessment of Heterogeneity in Lung Structure and Function During Mechanical Ventilation: A Review of Methodologies

The mammalian lung is characterized by heterogeneity in both its structure and function, by incorporating an asymmetric branching airway tree optimized for maintenance of efficient ventilation, perfusion, and gas exchange. Despite potential benefits of naturally occurring heterogeneity in the lungs, there may also be detrimental effects arising from pathologic processes, which may result in deficiencies in gas transport and exchange. Regardless of etiology, pathologic heterogeneity results in the maldistribution of regional ventilation and perfusion, impairments in gas exchange, and increased work of breathing. In extreme situations, heterogeneity may result in respiratory failure, necessitating support with a mechanical ventilator. This review will present a summary of measurement techniques for assessing and quantifying heterogeneity in respiratory system structure and function during mechanical ventilation. These methods have been grouped according to four broad categories: (1) inverse modeling of heterogeneous mechanical function; (2) capnography and washout techniques to measure heterogeneity of gas transport; (3) measurements of heterogeneous deformation on the surface of the lung; and finally (4) imaging techniques used to observe spatially-distributed ventilation or regional deformation. Each technique varies with regard to spatial and temporal resolution, degrees of invasiveness, risks posed to patients, as well as suitability for clinical implementation. Nonetheless, each technique provides a unique perspective on the manifestations and consequences of mechanical heterogeneity in the diseased lung.

Development and validation of a longitudinal soft-tissue metastatic lesion matching algorithm

Metastatic cancer presents with many, sometimes hundreds of metastatic lesions through the body, which often respond heterogeneously to treatment. Therefore, lesion-level assessment is necessary for a complete understanding of disease response. Lesion-level assessment typically requires manual matching of corresponding lesions, which is a tedious, subjective, and error-prone task. This study introduces a fully automated algorithm for matching of metastatic lesions in longitudinal medical images. The algorithm entails four steps: (1) image registration, (2) lesion dilation, (3) lesion clustering, and (4) linear assignment. In step (1), 3D deformable registration is used to register the scans. In step (2), lesion contours are conformally dilated. In step (3), lesion clustering is evaluated based on local metrics. In step (4), matching is assigned based on non-greedy cost minimization. The algorithm was optimized (e.g. choice of deformable registration algorithm, dilatation size) and validated on 140 scan-pairs of 32 metastatic cancer patients from two independent clinical trials, who received longitudinal PET/CT scans as part of their treatment response assessment. Registration error was evaluated using landmark distance. A sensitivity study was performed to evaluate the optimal lesion dilation magnitude. Lesion matching performance accuracy was evaluated for all patients and for a subset with high disease burden. Two investigated deformable registration approaches (whole body deformable and articulated deformable registrations) led to similar performance with the overall registration accuracy between 2.3 and 2.6 mm. The optimal dilation magnitude of 25 mm yielded almost a perfect matching accuracy of 0.98. No significant matching accuracy decrease was observed in the subset of patients with high lesion disease burden. In summary, lesion matching using our new algorithm was highly accurate and a significant improvement, when compared to previously established methods. The proposed method enables accurate automated metastatic lesion matching in whole-body longitudinal scans.

Improving deformable image registration with point metric and masking technique for postoperative breast cancer radiotherapy

Background

Deformable image registration (DIR) is increasingly used for target volume definition in radiotherapy. However, this method is challenging for postoperative breast cancer patients due to the large deformations and non-correspondence caused by tumor resection and clip insertion. In this study, an improved B-splines based DIR method was developed to address this issue for higher registration accuracy.

Methods

The conventional B-splines based DIR method was improved with the introduction of point metric and masking technique. The point metric minimizes the distance between 2 point sets with known correspondence for regularization of intensity-based B-splines registration. The masking technique reduces the influence of non-corresponding regions in breast computed tomography (CT) images. Two sets of CT images before and after breast surgery were used for image registration. One set was the diagnostic CT image acquired before surgery, and another set was the planning CT image acquired after surgery for breast cancer radiotherapy. A total of 26 sets of CT images from 13 patients were collected retrospectively for the test. The improved DIR method's registration accuracy was evaluated by target registration error (TRE), the Jacobian determinant, and visual assessment.

Results

For soft tissue, the difference in the median TRE between the improved DIR method and the conventional DIR method was statistically significant (2.27 vs. 5.88, P<0.05). The Jacobian determinant of the deformation field was positive for all patients. For visual assessment, the improved DIR method with point metric achieved better matching for soft tissue.

Conclusions

The improved DIR method's registration accuracy was higher than the conventional DIR method based on the preliminary results. With point metric and masking technique, the influence of large deformations and non-correspondence on registration between pre- and post-operative CT images can be effectively reduced. Therefore, this method provides a feasible way for target volume definition in postoperative breast cancer radiotherapy treatment planning.

An integrated 1D breathing lung simulation with relative hysteresis of airway structure and regional pressure for healthy and asthmatic human lungs

This study aims to develop a one-dimensional (1D) computational fluid dynamics (CFD) model with dynamic airway geometry that considers airway wall compliance and acinar dynamics. The proposed 1D model evaluates the pressure distribution and the hysteresis between the pressure and tidal volume (V_tidal) in the central and terminal airways for healthy and asthmatic subjects. Four-dimensional CT images were captured at 11-14 time points during the breathing cycle. The airway diameter and length were reconstructed using a volume-filling method and a stochastic model at respective time points. The obtained values for the airway diameter and length were interpolated via the Akima spline to avoid unboundedness. A 1D energy balance equation considering the effects of wall compliance and parenchymal inertance was solved using the efficient aggregation-based algebraic multigrid solver, a sparse matrix solver, reducing the computational costs by around 90% when compared with the generalized minimal residual solver. In the V_tidal versus displacement in the basal direction (z-coordinate), the inspiration curve was lower than the expiration curve, leading to relative hysteresis. The dynamic deformation model was the major factor influencing the difference in the workload in the central and terminal airways. In contrast, wall compliance and parenchymal inertance appeared only marginally to affect the pressure and workload. The integrated 1D model mimicked dynamic deformation by predicting airway diameter and length at each time point, describing the effects of wall compliance and parenchymal inertance. This computationally efficient model could be utilized to assess breathing mechanism as an alternative to pulmonary function tests.NEW & NOTEWORTHY This study introduces a one-dimensional (1D) computational fluid dynamics (CFD) model mimicking the realistic changes in diameter and length in whole airways and reveals differences in lung deformation between healthy and asthmatic subjects. Utilizing computational models, the effects of parenchymal inertance and airway wall compliance are investigated by changing ventilation frequency and airway wall elastance, respectively.

Imaging the Injured Lung: Mechanisms of Action and Clinical Use

Acute respiratory distress syndrome (ARDS) consists of acute hypoxemic respiratory failure characterized by massive and heterogeneously distributed loss of lung aeration caused by diffuse inflammation and edema present in interstitial and alveolar spaces. It is defined by consensus criteria, which include diffuse infiltrates on chest imaging-either plain radiography or computed tomography. This review will summarize how imaging sciences can inform modern respiratory management of ARDS and continue to increase the understanding of the acutely injured lung. This review also describes newer imaging methodologies that are likely to inform future clinical decision-making and potentially improve outcome. For each imaging modality, this review systematically describes the underlying principles, technology involved, measurements obtained, insights gained by the technique, emerging approaches, limitations, and future developments. Finally, integrated approaches are considered whereby multimodal imaging may impact management of ARDS.

A robust deformable image registration enhancement method based on radial basis function

Background

To develop and evaluate a robust deformable image registration (DIR) enhancement method based on radial basis function (RBF) expansion.

Methods

To improve DIR accuracy using sparsely available measured displacements, it is crucial to estimate the motion correlation between the voxels. In the proposed method, we chose to derive this correlation from the initial displacement vector fields (DVFs), and represent it in the form of RBF expansion coefficients of the voxels. The method consists of three steps: (I) convert an initial DVF to a coefficient matrix comprising expansion coefficients of the Wendland's RBF; (II) modify the coefficient matrix under the guidance of sparely distributed landmarks to generate the post-enhancement coefficient matrix; and (III) convert the post-enhancement coefficient matrix to the post-enhancement DVF. The method was tested on five DIR algorithms using a digital phantom. 3D registration errors were calculated for comparisons between the pre-/post-enhancement DVFs and the ground-truth DVFs. Effects of the number and locations of landmarks on DIR enhancement were evaluated.

Results

After applying the DIR enhancement method, the 3D registration errors per voxel (unit: mm) were reduced from pre-enhancement to post-enhancement by 1.3 (2.4 to 1.1, 54.2%), 0.0 (0.9 to 0.9, 0.0%), 6.1 (8.2 to 2.1, 74.4%), 3.2 (4.7 to 1.5, 68.1%), and 1.7 (2.9 to 1.2, 58.6%) for the five tested DIR algorithms respectively. The average DIR error reduction was 2.5±2.3 mm (percentage error reduction: 51.1%±29.1%). 3D registration errors decreased inverse-exponentially as the number of landmarks increased, and were insensitive to the landmarks' locations in relation to the down-sampling DVF grids.

Conclusions

We demonstrated the feasibility of a robust RBF-based method for enhancing DIR accuracy using sparsely distributed landmarks. This method has been shown robust and effective in reducing DVF errors using different numbers and distributions of landmarks for various DIR algorithms.

In vivo dynamics of the tracheal airway and its influences on respiratory airflows

Respiration is a dynamic process accompanied by morphological changes in the airways. Although deformation of large airways is expected to exacerbate pulmonary disease symptoms by obstructing airflow during increased minute ventilation, its quantitative effects on airflow characteristics remain unclear. Here, we used an exemplar case derived from in vivo dynamic imaging and examined the effects of tracheal deformation on airflow characteristics under different conditions. First, we measured tracheal deformation profiles of a healthy lung using magnetic resonance imaging during forced exhalation, which we simulated to characterize subject-specific airflow patterns. Subsequently, for both inhalation and exhalation, we compared the airflows when the maximal deformation in tracheal cross-sectional area was 0% (rigid), 33% (mild), 50% (moderate), or 75% (severe). We quantified differences in airflow patterns between deformable and rigid airways by computing the correlation coefficients (R) and the root-mean-square of differences (Drms) between their velocity contours. For both inhalation and exhalation, airflow patterns were similar in all branches between the rigid and mild conditions (R > 0.9; Drms < 32%). However, airflow characteristics in the moderate and severe conditions differed markedly from those in the rigid and mild conditions in all lung branches, particularly for inhalation (moderate: R > 0.1, Drms < 76%; severe: R > 0.2, Drms < 96%). Our exemplar case supports the use of a rigid airway assumption to compute flows for mild deformation. For moderate or severe deformation, however, dynamic contraction should be considered, especially during inhalation, to accurately predict airflow and elucidate the underlying pulmonary pathology.

A novel method to generate dynamic boundary conditions for airway CFD by mapping upper airway movement with non-rigid registration of dynamic and static MRI

Computational fluid dynamics (CFD) simulations of airflow in the human airways have the potential to provide a great deal of information that can aid clinicians in case management and surgical decision making, such as airway resistance, energy expenditure, airflow distribution, heat and moisture transfer, and particle deposition, as well as the change in each of these due to surgical interventions. However, the clinical relevance of CFD simulations has been limited to date, as previous models either did not incorporate neuromuscular motion or any motion at all. Many common airway pathologies, such as obstructive sleep apnea (OSA) and tracheomalacia, involve large movements of the structures surrounding the airway, such as the tongue and soft palate. Airway wall motion may be due to many factors including neuromuscular motion, internal aerodynamic forces, and external forces such as gravity. Therefore, to realistically model these airway diseases, a method is required to derive the airway wall motion, whatever the cause, and apply it as a boundary condition to CFD simulations. This paper presents and validates a novel method of capturing in vivo motion of airway walls from magnetic resonance images with high spatiotemporal resolution, through a novel combination of non-rigid image, surface, and surface-normal-vector registration. Coupled with image-synchronous pneumotachography, this technique provides the necessary boundary conditions for dynamic CFD simulations of breathing, allowing the effect of the airway's complex motion to be calculated for the first time, in both normal subjects and those with conditions such as OSA.

Quantitative computed tomography determined regional lung mechanics in normal nonsmokers, normal smokers and metastatic sarcoma subjects

Objectives

Extra-thoracic tumors send out pilot cells that attach to the pulmonary endothelium. We hypothesized that this could alter regional lung mechanics (tissue stiffening or accumulation of fluid and inflammatory cells) through interactions with host cells. We explored this with serial inspiratory computed tomography (CT) and image matching to assess regional changes in lung expansion.

Materials and methods

We retrospectively assessed 44 pairs of two serial CT scans on 21 sarcoma patients: 12 without lung metastases and 9 with lung metastases. For each subject, two or more serial inspiratory clinically-derived CT scans were retrospectively collected. Two research-derived control groups were included: 7 normal nonsmokers and 12 asymptomatic smokers with two inspiratory scans taken the same day or one year apart respectively. We performed image registration for local-to-local matching scans to baseline, and derived local expansion and density changes at an acinar scale. Welch two sample t test was used for comparison between groups. Statistical significance was determined with a p value < 0.05.

Results

Lung regions of metastatic sarcoma patients (but not the normal control group) demonstrated an increased proportion of normalized lung expansion between the first and second CT. These hyper-expanded regions were associated with, but not limited to, visible metastatic lung lesions. Compared with the normal control group, the percent of increased normalized hyper-expanded lung in sarcoma subjects was significantly increased (p < 0.05). There was also evidence of increased lung “tissue” volume (non-air components) in the hyper-expanded regions of the cancer subjects relative to non-hyper-expanded regions. “Tissue” volume increase was present in the hyper-expanded regions of metastatic and non-metastatic sarcoma subjects. This putatively could represent regional inflammation related to the presence of tumor pilot cell-host related interactions.

Conclusions

This new quantitative CT (QCT) method for linking serial acquired inspiratory CT images may provide a diagnostic and prognostic means to objectively characterize regional responses in the lung following oncological treatment and monitoring for lung metastases.

Patient-specific respiratory motion estimation using Sparse Motion Field Presentation

Respiratory motion estimation plays a significant role in radiation therapy. Previous motion estimation approaches usually depended on 4DCT, which introduced extra radio dose for patients, and the local motion details were ignored in the statistical model. In this paper, we propose a novel estimation framework, which employs the Sparse Motion Field Presentation (SMFP) method to obtain a coarse motion estimation which preserves patient-specific respiratory motion details and an Adaptive Variable Coefficient (AVC) motion prior registration approach is applied for the accurate estimation. The experimental results show that the proposed framework effectively preserved the local motion details and achieved more accurate motion estimations compared to the Mean Motion Model (MMM) and the Principal Component Analysis (PCA) model. We achieved motion estimations for diaphragmatic breathing type, thoracic breathing type and mixed type, respectively. The accuracy measured in the average symmetric surface distance (standard deviation) were 1.9(0.9) mm, 2.4(1.1) mm and 2.2(1.0) mm, when the sum of squared intensity difference (SSD) were 5.0, 6.1 and 5.6, respectively.

Automated identification and reduction of artifacts in cine four-dimensional computed tomography (4DCT) images using respiratory motion model

PURPOSE

Four-dimensional computed tomography (4DCT) images are often marred by artifacts that substantially degrade image quality and confound image interpretation. Human observation remains the standard method of 4DCT artifact evaluation, which is time-consuming and subjective. We develop a method to automatically identify and reduce artifacts in cine 4DCT images.

METHODS

We proposed an algorithm that consisted of two main stages: deformable image registration and respiratory motion simulation. Specifically, each 4DCT phase image was registered to the breath-holding CT image using the block-matching method, with erroneous spatial matches removed by the least median of squares filter and the full displacement vector field generated by the moving least squares interpolation. The lung's respiratory motion trajectory was then recovered from the displacement vector field using the parameterized polynomial function, with fitting parameters estimated by combinatorial optimization. In this way, artifacts were located according to deviations between image points and their motion trajectories and further corrected based on position prediction.

RESULTS

The mean spatial error (standard deviation) was 1.00 (0.85) mm after registration as opposed to 6.96 (4.61) mm before registration. In addition, we took human observation conducted by medical experts as the gold standard. The average sensitivity, specificity, and accuracy of the proposed method in artifact identification were 0.97, 0.84, and 0.89, respectively.

CONCLUSIONS

The proposed method identified and reduced artifacts accurately and automatically, providing an alternative way to analyze 4DCT image quality and to correct problematic images for radiation therapy.

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry.

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry.

Selecting registration schemes in case of interstitial lung disease follow-up in CT

PURPOSE

Primary goal of this study is to select optimal registration schemes in the framework of interstitial lung disease (ILD) follow-up analysis in CT.

METHODS

A set of 128 multiresolution schemes composed of multiresolution nonrigid and combinations of rigid and nonrigid registration schemes are evaluated, utilizing ten artificially warped ILD follow-up volumes, originating from ten clinical volumetric CT scans of ILD affected patients, to select candidate optimal schemes. Specifically, all combinations of four transformation models (three rigid: rigid, similarity, affine and one nonrigid: third order B-spline), four cost functions (sum-of-square distances, normalized correlation coefficient, mutual information, and normalized mutual information), four gradient descent optimizers (standard, regular step, adaptive stochastic, and finite difference), and two types of pyramids (recursive and Gaussian-smoothing) were considered. The selection process involves two stages. The first stage involves identification of schemes with deformation field singularities, according to the determinant of the Jacobian matrix. In the second stage, evaluation methodology is based on distance between corresponding landmark points in both normal lung parenchyma (NLP) and ILD affected regions. Statistical analysis was performed in order to select near optimal registration schemes per evaluation metric. Performance of the candidate registration schemes was verified on a case sample of ten clinical follow-up CT scans to obtain the selected registration schemes.

RESULTS

By considering near optimal schemes common to all ranking lists, 16 out of 128 registration schemes were initially selected. These schemes obtained submillimeter registration accuracies in terms of average distance errors 0.18 ± 0.01 mm for NLP and 0.20 ± 0.01 mm for ILD, in case of artificially generated follow-up data. Registration accuracy in terms of average distance error in clinical follow-up data was in the range of 1.985-2.156 mm and 1.966-2.234 mm, for NLP and ILD affected regions, respectively, excluding schemes with statistically significant lower performance (Wilcoxon signed-ranks test, p < 0.05), resulting in 13 finally selected registration schemes.

CONCLUSIONS

Selected registration schemes in case of ILD CT follow-up analysis indicate the significance of adaptive stochastic gradient descent optimizer, as well as the importance of combined rigid and nonrigid schemes providing high accuracy and time efficiency. The selected optimal deformable registration schemes are equivalent in terms of their accuracy and thus compatible in terms of their clinical outcome.

Efficient methods for implementation of multi-level nonrigid mass-preserving image registration on GPUs and multi-threaded CPUs

BACKGROUND AND OBJECTIVE

Faster and more accurate methods for registration of images are important for research involved in conducting population-based studies that utilize medical imaging, as well as improvements for use in clinical applications. We present a novel computation- and memory-efficient multi-level method on graphics processing units (GPU) for performing registration of two computed tomography (CT) volumetric lung images.

METHODS

We developed a computation- and memory-efficient Diffeomorphic Multi-level B-Spline Transform Composite (DMTC) method to implement nonrigid mass-preserving registration of two CT lung images on GPU. The framework consists of a hierarchy of B-Spline control grids of increasing resolution. A similarity criterion known as the sum of squared tissue volume difference (SSTVD) was adopted to preserve lung tissue mass. The use of SSTVD consists of the calculation of the tissue volume, the Jacobian, and their derivatives, which makes its implementation on GPU challenging due to memory constraints. The use of the DMTC method enabled reduced computation and memory storage of variables with minimal communication between GPU and Central Processing Unit (CPU) due to ability to pre-compute values. The method was assessed on six healthy human subjects.

RESULTS

Resultant GPU-generated displacement fields were compared against the previously validated CPU counterpart fields, showing good agreement with an average normalized root mean square error (nRMS) of 0.044±0.015. Runtime and performance speedup are compared between single-threaded CPU, multi-threaded CPU, and GPU algorithms. Best performance speedup occurs at the highest resolution in the GPU implementation for the SSTVD cost and cost gradient computations, with a speedup of 112 times that of the single-threaded CPU version and 11 times over the twelve-threaded version when considering average time per iteration using a Nvidia Tesla K20X GPU.

CONCLUSIONS

The proposed GPU-based DMTC method outperforms its multi-threaded CPU version in terms of runtime. Total registration time reduced runtime to 2.9min on the GPU version, compared to 12.8min on twelve-threaded CPU version and 112.5min on a single-threaded CPU. Furthermore, the GPU implementation discussed in this work can be adapted for use of other cost functions that require calculation of the first derivatives.

Pulmonary CT and MRI phenotypes that help explain chronic pulmonary obstruction disease pathophysiology and outcomes

Pulmonary x-ray computed tomographic (CT) and magnetic resonance imaging (MRI) research and development has been motivated, in part, by the quest to subphenotype common chronic lung diseases such as chronic obstructive pulmonary disease (COPD). For thoracic CT and MRI, the main COPD research tools, disease biomarkers are being validated that go beyond anatomy and structure to include pulmonary functional measurements such as regional ventilation, perfusion, and inflammation. In addition, there has also been a drive to improve spatial and contrast resolution while at the same time reducing or eliminating radiation exposure. Therefore, this review focuses on our evolving understanding of patient-relevant and clinically important COPD endpoints and how current and emerging MRI and CT tools and measurements may be exploited for their identification, quantification, and utilization. Since reviews of the imaging physics of pulmonary CT and MRI and reviews of other COPD imaging methods were previously published and well-summarized, we focus on the current clinical challenges in COPD and the potential of newly emerging MR and CT imaging measurements to address them. Here we summarize MRI and CT imaging methods and their clinical translation for generating reproducible and sensitive measurements of COPD related to pulmonary ventilation and perfusion as well as parenchyma morphology. The key clinical problems in COPD provide an important framework in which pulmonary imaging needs to rapidly move in order to address the staggering burden, costs, as well as the mortality and morbidity associated with COPD.

Use of registration-based contour propagation in texture analysis for esophageal cancer pathologic response prediction

Change in PET-based textural features has shown promise in predicting cancer response to treatment. However, contouring tumour volumes on longitudinal scans is time-consuming. This study investigated the usefulness of contour propagation in texture analysis for the purpose of pathologic response prediction in esophageal cancer. Forty-five esophageal cancer patients underwent PET/CT scans before and after chemo-radiotherapy. Patients were classified into responders and non-responders after the surgery. Physician-defined tumour ROIs on pre-treatment PET were propagated onto the post-treatment PET using rigid and ten deformable registration algorithms. PET images were converted into 256 discrete values. Co-occurrence, run-length, and size zone matrix textures were computed within all ROIs. The relative difference of each texture at different treatment time-points was used to predict the pathologic responders. Their predictive value was assessed using the area under the receiver-operating-characteristic curve (AUC). Propagated ROIs from different algorithms were compared using Dice similarity index (DSI). Contours propagated by the fast-demons, fast-free-form and rigid algorithms did not fully capture the high FDG uptake regions of tumours. Fast-demons propagated ROIs had the least agreement with other contours (DSI = 58%). Moderate to substantial overlap were found in the ROIs propagated by all other algorithms (DSI = 69%-79%). Rigidly propagated ROIs with co-occurrence texture failed to significantly differentiate between responders and non-responders (AUC = 0.58, q-value = 0.33), while the differentiation was significant with other textures (AUC = 0.71-0.73, p < 0.009). Among the deformable algorithms, fast-demons (AUC = 0.68-0.70, q-value < 0.03) and fast-free-form (AUC = 0.69-0.74, q-value < 0.04) were the least predictive. ROIs propagated by all other deformable algorithms with any texture significantly predicted pathologic responders (AUC = 0.72-0.78, q-value < 0.01). Propagated ROIs using deformable registration for all textures can lead to accurate prediction of pathologic response, potentially expediting the temporal texture analysis process. However, fast-demons, fast-free-form, and rigid algorithms should be applied with care due to their inferior performance compared to other algorithms.

Extent of Spine Deformity Predicts Lung Growth and Function in Rabbit Model of Early Onset Scoliosis

Early onset deformity of the spine and chest wall (initiated <8 years of age) is associated with increased morbidity at adulthood relative to adolescent onset deformity of comparable severity. Presumably, inhibition of thoracic growth during late stage alveolarization leads to an irreversible loss of pulmonary growth and thoracic function; however the natural history of this disease from onset to adulthood has not been well characterized. In this study we establish a rabbit model of early onset scoliosis to establish the extent that thoracic deformity affects structural and functional respiratory development. Using a surgical right unilateral rib-tethering procedure, rib fusion with early onset scoliosis was induced in 10 young New Zealand white rabbits (3 weeks old). Progression of spine deformity, functional residual capacity, total lung capacity, and lung mass was tracked through longitudinal breath-hold computed tomography imaging up to skeletal maturity (28 weeks old). Additionally at maturity forced vital capacity and regional specific volume were calculated as functional measurements and histo-morphometry performed with the radial alveolar count as a measure of acinar complexity. Data from tethered rib rabbits were compared to age matched healthy control rabbits (N = 8). Results show unilateral rib-tethering created a progressive spinal deformity ranging from 30° to 120° curvature, the severity of which was strongly associated with pulmonary growth and functional outcomes. At maturity rabbits with deformity greater than the median (55°) had decreased body weight (89%), right (59%) and left (86%) lung mass, right (74%) and left (69%) radial alveolar count, right lung volume at total lung capacity (60%), and forced vital capacity (75%). Early treatment of spinal deformity in children may prevent pulmonary complications in adulthood and these results provide a basis for the prediction of pulmonary development from thoracic structure. This model may also have future use as a platform to evaluate treatment effectiveness.

Towards local progression estimation of pulmonary emphysema using CT

PURPOSE

Whole lung densitometry on chest CT images is an accepted method for measuring tissue destruction in patients with pulmonary emphysema in clinical trials. Progression measurement is required for evaluation of change in health condition and the effect of drug treatment. Information about the location of emphysema progression within the lung may be important for the correct interpretation of drug efficacy, or for determining a treatment plan. The purpose of this study is therefore to develop and validate methods that enable the local measurement of lung density changes, which requires proper modeling of the effect of respiration on density.

METHODS

Four methods, all based on registration of baseline and follow-up chest CT scans, are compared. The first naïve method subtracts registered images. The second employs the so-called dry sponge model, where volume correction is performed using the determinant of the Jacobian of the transformation. The third and the fourth introduce a novel adaptation of the dry sponge model that circumvents its constant-mass assumption, which is shown to be invalid. The latter two methods require a third CT scan at a different inspiration level to estimate the patient-specific density-volume slope, where one method employs a global and the other a local slope. The methods were validated on CT scans of a phantom mimicking the lung, where mass and volume could be controlled. In addition, validation was performed on data of 21 patients with pulmonary emphysema.

RESULTS

The image registration method was optimized leaving a registration error below half the slice increment (median 1.0 mm). The phantom study showed that the locally adapted slope model most accurately measured local progression. The systematic error in estimating progression, as measured on the phantom data, was below 2 gr/l for a 70 ml (6%) volume difference, and 5 gr/l for a 210 ml (19%) difference, if volume correction was applied. On the patient data an underlying linearity assumption relating lung volume change with density change was shown to hold (fitR(2) = 0.94), and globalized versions of the local models are consistent with global results (R(2) of 0.865 and 0.882 for the two adapted slope models, respectively).

CONCLUSIONS

In conclusion, image matching and subsequent analysis of differences according to the proposed lung models (i) has good local registration accuracy on patient data, (ii) effectively eliminates a dependency on inspiration level at acquisition time, (iii) accurately predicts progression in phantom data, and (iv) is reasonably consistent with global results in patient data. It is therefore a potential future tool for assessing local emphysema progression in drug evaluation trials and in clinical practice.

Analysis of three-dimensional aerosol deposition in pharmacologically relevant terms: beyond black or white ROIs

BACKGROUND

This article presents a novel methodological approach to evaluate images of aerosol deposition taken with PET-CT cameras. Traditionally, Black-or-White (BW) Regions of Interest (ROIs) are created to cover Anatomical Regions (ARs) segmented from the high-resolution CT. Such ROIs do not usually consider blurring effects due to limited spatial resolution or breathing motion, and do not consider uncertainty in the AR position within the PET image. The new methodology presented here (Grayscale) addresses these issues, allows estimates of aerosol deposition within ARs, and expresses the deposition in terms of Tissue Dosing (in the lung periphery) and Inner Surface Concentration (in the larger airways).

METHODS

Imaging data included a PET deposition image acquired during breathing and two CT scans acquired during breath holds at different lung volumes. The lungs were segmented into anatomically consistent ARs to allow unbiased comparisons across subjects and across lobes. The Grayscale method involves defining Voxel Influence Matrices (VIMs) to consider how average activity within each AR influences the measured activity within each voxel. The BW and Grayscale methods were used to analyze aerosol deposition in 14 bronchoconstricted asthmatics.

RESULTS

Grayscale resulted in a closer description of the PET image than BW (p<0.0001) and exposed a seven-fold underestimation in measures of specific deposition. The Average Tissue Dosing was 2.11×10(-6) Total Lung Dose/mg. The average Inner Surface Concentration was 45×10(-6) Total Lung Dose/mm(2), with the left lower lobe having a lower ISC than lobes of the right lung (p<0.05). There was a strong lobar heterogeneity in these measures (COV=0.3).

CONCLUSION

The Grayscale approach is an improvement over the BW approach and provides a closer description of the PET image. It can be used to characterize heterogeneous concentrations throughout the lung and may be important in translational research and in the evaluation of aerosol delivery systems.

A fast deformable registration method for 4D lung CT in hybrid framework

PURPOSE

A pulmonary respiration model for deformable registration of lung CT for the surgery path planning and surgical navigation is an important, difficult, and time-consuming task. This paper presents a new fast deformable registration method for 4D lung CT in a hybrid framework incorporating point set registration with mutual information registration.

METHOD

The point sets of the lung surface and vessels are automatically extracted. Their displacement vectors are obtained by point set registration. The sum of squared Euclidean distance between the displacement vectors of these point sets and the displacement vectors based on the B-spline transformation model is minimized as a novel similarity measure to derive the rough transformation function. Finally, the rough transformation function is refined by using the mutual information-based registration method. To evaluate the effectiveness of the proposed method, the authors performed registrations on 20 4D lung volume cases from two different CT scanners. The proposed method was compared with the point set-based method, the mutual information-based method, and the ANTS method, which is a state-of-the-art deformable registration technique.

RESULTS

The results show that the landmark distance errors and computation time of the proposed method decreased an average of 5 and 70 %, respectively, when compared to the mutual information-alone-based method. The proposed method results in an average of 28 % lower landmark distance error than registration method based on point sets in spite of increase in computation time. Moreover, compared with ANTS, the computation time of the proposed method is reduced by an average of 93 % in the case of comparable landmark distance errors.

CONCLUSION

The accuracy and speed of the proposed deformable registration method indicate that the method is suitable for use in a clinical image-guided intervention system.

Deformable registration of the inflated and deflated lung in cone-beam CT-guided thoracic surgery: initial investigation of a combined model- and image-driven approach

PURPOSE

Surgical resection is the preferred modality for curative treatment of early stage lung cancer, but localization of small tumors (<10 mm diameter) during surgery presents a major challenge that is likely to increase as more early-stage disease is detected incidentally and in low-dose CT screening. To overcome the difficulty of manual localization (fingers inserted through intercostal ports) and the cost, logistics, and morbidity of preoperative tagging (coil or dye placement under CT-fluoroscopy), the authors propose the use of intraoperative cone-beam CT (CBCT) and deformable image registration to guide targeting of small tumors in video-assisted thoracic surgery (VATS). A novel algorithm is reported for registration of the lung from its inflated state (prior to pleural breach) to the deflated state (during resection) to localize surgical targets and adjacent critical anatomy.

METHODS

The registration approach geometrically resolves images of the inflated and deflated lung using a coarse model-driven stage followed by a finer image-driven stage. The model-driven stage uses image features derived from the lung surfaces and airways: triangular surface meshes are morphed to capture bulk motion; concurrently, the airways generate graph structures from which corresponding nodes are identified. Interpolation of the sparse motion fields computed from the bounding surface and interior airways provides a 3D motion field that coarsely registers the lung and initializes the subsequent image-driven stage. The image-driven stage employs an intensity-corrected, symmetric form of the Demons method. The algorithm was validated over 12 datasets, obtained from porcine specimen experiments emulating CBCT-guided VATS. Geometric accuracy was quantified in terms of target registration error (TRE) in anatomical targets throughout the lung, and normalized cross-correlation. Variations of the algorithm were investigated to study the behavior of the model- and image-driven stages by modifying individual algorithmic steps and examining the effect in comparison to the nominal process.

RESULTS

The combined model- and image-driven registration process demonstrated accuracy consistent with the requirements of minimally invasive VATS in both target localization (∼3-5 mm within the target wedge) and critical structure avoidance (∼1-2 mm). The model-driven stage initialized the registration to within a median TRE of 1.9 mm (95% confidence interval (CI) maximum = 5.0 mm), while the subsequent image-driven stage yielded higher accuracy localization with 0.6 mm median TRE (95% CI maximum = 4.1 mm). The variations assessing the individual algorithmic steps elucidated the role of each step and in some cases identified opportunities for further simplification and improvement in computational speed.

CONCLUSIONS

The initial studies show the proposed registration method to successfully register CBCT images of the inflated and deflated lung. Accuracy appears sufficient to localize the target and adjacent critical anatomy within ∼1-2 mm and guide localization under conditions in which the target cannot be discerned directly in CBCT (e.g., subtle, nonsolid tumors). The ability to directly localize tumors in the operating room could provide a valuable addition to the VATS arsenal, obviate the cost, logistics, and morbidity of preoperative tagging, and improve patient safety. Future work includes in vivo testing, optimization of workflow, and integration with a CBCT image guidance system.

Multiscale image-based modeling and simulation of gas flow and particle transport in the human lungs

Improved understanding of structure and function relationships in the human lungs in individuals and subpopulations is fundamentally important to the future of pulmonary medicine. Image-based measures of the lungs can provide sensitive indicators of localized features, however to provide a better prediction of lung response to disease, treatment, and environment, it is desirable to integrate quantifiable regional features from imaging with associated value-added high-level modeling. With this objective in mind, recent advances in computational fluid dynamics (CFD) of the bronchial airways-from a single bifurcation symmetric model to a multiscale image-based subject-specific lung model-will be reviewed. The interaction of CFD models with local parenchymal tissue expansion-assessed by image registration-allows new understanding of the interplay between environment, hot spots where inhaled aerosols could accumulate, and inflammation. To bridge ventilation function with image-derived central airway structure in CFD, an airway geometrical modeling method that spans from the model 'entrance' to the terminal bronchioles will be introduced. Finally, the effects of turbulent flows and CFD turbulence models on aerosol transport and deposition will be discussed.

A multiscale MDCT image-based breathing lung model with time-varying regional ventilation

A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C₁ continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.

Deformable medical image registration: a survey

Deformable image registration is a fundamental task in medical image processing. Among its most important applications, one may cite: 1) multi-modality fusion, where information acquired by different imaging devices or protocols is fused to facilitate diagnosis and treatment planning; 2) longitudinal studies, where temporal structural or anatomical changes are investigated; and 3) population modeling and statistical atlases used to study normal anatomical variability. In this paper, we attempt to give an overview of deformable registration methods, putting emphasis on the most recent advances in the domain. Additional emphasis has been given to techniques applied to medical images. In order to study image registration methods in depth, their main components are identified and studied independently. The most recent techniques are presented in a systematic fashion. The contribution of this paper is to provide an extensive account of registration techniques in a systematic manner.

A PERFECT MATCH CONDITION FOR POINT-SET MATCHING PROBLEMS USING THE OPTIMAL MASS TRANSPORT APPROACH

We study the performance of optimal mass transport-based methods applied to point-set matching problems. The present study, which is based on the L2 mass transport cost, states that perfect matches always occur when the product of the point-set cardinality and the norm of the curl of the non-rigid deformation field does not exceed some constant. This analytic result is justified by a numerical study of matching two sets of pulmonary vascular tree branch points whose displacement is caused by the lung volume changes in the same human subject. The nearly perfect match performance verifies the effectiveness of this mass transport-based approach.

Multiscale imaging and registration-driven model for pulmonary acinar mechanics in the mouse

A registration-based multiscale method to obtain a deforming geometric model of mouse acinus is presented. An intact mouse lung was fixed by means of vascular perfusion at a hydrostatic inflation pressure of 20 cmH(2)O. Microcomputed tomography (μCT) scans were obtained at multiple resolutions. Substructural morphometric analysis of a complete acinus was performed by computing a surface-to-volume (S/V) ratio directly from the 3D reconstruction of the acinar geometry. A geometric similarity is observed to exist in the acinus where S/V is approximately preserved anywhere in the model. Using multiscale registration, the shape of the acinus at an elevated inflation pressure of 25 cmH(2)O is estimated. Changes in the alveolar geometry suggest that the deformation within the acinus is not isotropic. In particular, the expansion of the acinus (from 20 to 25 cmH(2)O) is accompanied by an increase in both surface area and volume in such a way that the S/V ratio is not significantly altered. The developed method forms a useful tool in registration-driven fluid and solid mechanics studies as displacement of the alveolar wall becomes available in a discrete sense.

Improving Intensity-Based Lung CT Registration Accuracy Utilizing Vascular Information

Accurate pulmonary image registration is a challenging problem when the lungs have a deformation with large distance. In this work, we present a nonrigid volumetric registration algorithm to track lung motion between a pair of intrasubject CT images acquired at different inflation levels and introduce a new vesselness similarity cost that improves intensity-only registration. Volumetric CT datasets from six human subjects were used in this study. The performance of four intensity-only registration algorithms was compared with and without adding the vesselness similarity cost function. Matching accuracy was evaluated using landmarks, vessel tree, and fissure planes. The Jacobian determinant of the transformation was used to reveal the deformation pattern of local parenchymal tissue. The average matching error for intensity-only registration methods was on the order of 1 mm at landmarks and 1.5 mm on fissure planes. After adding the vesselness preserving cost function, the landmark and fissure positioning errors decreased approximately by 25% and 30%, respectively. The vesselness cost function effectively helped improve the registration accuracy in regions near thoracic cage and near the diaphragm for all the intensity-only registration algorithms tested and also helped produce more consistent and more reliable patterns of regional tissue deformation.

Tracking the motion trajectories of junction structures in 4D CT images of the lung

Respiratory motion poses a major challenge in lung radiotherapy. Based on 4D CT images, a variety of intensity-based deformable registration techniques have been proposed to study the pulmonary motion. However, the accuracy achievable with these approaches can be sub-optimal because the deformation is defined globally in space. Therefore, the accuracy of the alignment of local structures may be compromised. In this work, we propose a novel method to detect a large collection of natural junction structures in the lung and use them as the reliable markers to track the lung motion. Specifically, detection of the junction centers and sizes is achieved by analysis of local shape profiles on one segmented image. To track the temporal trajectory of a junction, the image intensities within a small region of interest surrounding the center are selected as its signature. Under the assumption of the cyclic motion, we describe the trajectory by a closed B-spline curve and search for the control points by maximizing a metric of combined correlation coefficients. Local extrema are suppressed by improving the initial conditions using random walks from pair-wise optimizations. Several descriptors are introduced to analyze the motion trajectories. Our method was applied to 13 real 4D CT images. More than 700 junctions in each case are detected with an average positive predictive value of greater than 90%. The average tracking error between automated and manual tracking is sub-voxel and smaller than the published results using the same set of data.

Introducing the Jacobian-volume-histogram of deforming organs: application to parotid shrinkage evaluation

The Jacobian of the deformation field of elastic registration between images taken during radiotherapy is a measure of inter-fraction local deformation. The histogram of the Jacobian values (Jac) within an organ was introduced (JVH-Jacobian-volume-histogram) and first applied in quantifying parotid shrinkage. MVCTs of 32 patients previously treated with helical tomotherapy for head-neck cancers were collected. Parotid deformation was evaluated through elastic registration between MVCTs taken at the first and last fractions. Jac was calculated for each voxel of all parotids, and integral JVHs were calculated for each parotid; the correlation between the JVH and the planning dose-volume histogram (DVH) was investigated. On average, 82% (±17%) of the voxels shrinks (Jac < 1) and 14% (±17%) shows a local compression >50% (Jac < 0.5). The best correlation between the DVH and the JVH was found between V10 and V15, and Jac < 0.4-0.6 (p < 0.01). The best constraint predicting a higher number of largely compressing voxels (Jac0.5<7.5%, median value) was V15 ≥ 75% (OR: 7.6, p = 0.002). Jac and the JVH are promising tools for scoring/modelling toxicity and for evaluating organ/contour variations with potential applications in adaptive radiotherapy.