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

Emerging Imaging Modalities for Functional Assessment of Donor Lungs Ex Vivo

The severe shortage of functional donor lungs that can be offered to recipients has been a major challenge in lung transplantation. Innovative ex vivo lung perfusion (EVLP) and tissue engineering methodologies are now being developed to repair damaged donor lungs that are deemed unsuitable for transplantation. To assess the efficacy of donor lung reconditioning methods intended to rehabilitate rejected donor lungs, monitoring of lung function with improved spatiotemporal resolution is needed. Recent developments in live imaging are enabling non-destructive, direct, and longitudinal modalities for assessing local tissue and whole lung functions. In this review, we describe how emerging live imaging modalities can be coupled with lung tissue engineering approaches to promote functional recovery of ex vivo donor lungs.

Joint Deformable Image Registration and ADC Map Regularization: Application to DWI-based Lymphoma Classification

The Apparent Diffusion Coefficient (ADC) is considered an important imaging biomarker contributing to the assessment of tissue microstructure and pathophysiology. It is calculated from Diffusion-Weighted Magnetic Resonance Imaging (DWI) by means of a diffusion model, usually without considering any motion during image acquisition. We propose a method to improve the computation of the ADC by coping jointly with both motion artifacts in whole-body DWI (through group-wise registration) and possible instrumental noise in the diffusion model. The proposed deformable registration method yielded on average the lowest ADC reconstruction error on data with simulated motion and diffusion. Moreover, our approach was applied on whole-body diffusion weighted images obtained with five different b-values from a cohort of 38 patients with histologically confirmed lymphomas of three different types (Hodgkin, diffuse large B-cell lymphoma and follicular lymphoma). Evaluation on the real data showed that ADC-based features, extracted using our joint optimization approach classified lymphomas with an accuracy of approximately 78.6\% (yielding a 11\% increase in respect to the standard features extracted from unregistered diffusion-weighted images). Furthermore, the correlation between diffusion characteristics and histopathological findings was higher than any other previous approach of ADC computation.

Lung Respiratory Motion Estimation Based on Fast Kalman Filtering and 4D CT Image Registration

Respiratory motion estimation is an important part in image-guided radiation therapy and clinical diagnosis. However, most of the respiratory motion estimation methods rely on indirect measurements of external breathing indicators, which will not only introduce great estimation errors, but also bring invasive injury for patients. In this paper, we propose a method of lung respiratory motion estimation based on fast Kalman filtering and 4D CT image registration (LRME-4DCT). In order to perform dynamic motion estimation for continuous phases, a motion estimation model is constructed by combining two kinds of GPU-accelerated 4D CT image registration methods with fast Kalman filtering method. To address the high computational requirements of 4D CT image sequences, a multi-level processing strategy is adopted in the 4D CT image registration methods, and respiratory motion states are predicted from three independent directions. In the DIR-lab dataset and POPI dataset with 4D CT images, the average target registration error (TRE) of the LRME-4DCT method can reach 0.91 mm and 0.85 mm respectively. Compared with traditional estimation methods based on pair-wise image registration, the proposed LRME-4DCT method can estimate the physiological respiratory motion more accurately and quickly. Our proposed LRME-4DCT method fully meets the practical clinical requirements for rapid dynamic estimation of lung respiratory motion.

3D Magnetic Resonance Spirometry

Spirometry is today the gold standard technique for assessing pulmonary ventilatory function in humans. From the shape of a flow-volume loop measured while the patient is performing forced respiratory cycles, the Forced Vital Capacity (FVC) and the Forced Expiratory Volume in one second (FEV1) can be inferred, and the pulmonologist is able to detect and characterize common respiratory afflictions. This technique is non-invasive, simple, widely available, robust, repeatable and reproducible. Yet, its outcomes rely on the patient's cooperation and provide only global information over the lung. With 3D Magnetic Resonance (MR) Spirometry, local ventilation can be assessed by MRI anywhere in the lung while the patient is freely breathing. The larger dimensionality of 3D MR Spirometry advantageously allows the extraction of original metrics that characterize the anisotropic and hysteretic regional mechanical behavior of the lung. Here, we demonstrated the potential of this technique on a healthy human volunteer breathing along different respiratory patterns during the MR acquisition. These new results are discussed with lung physiology and recent pulmonary CT data. As respiratory mechanics inherently support lung ventilation, 3D MR Spirometry may open a new way to non-invasively explore lung function while providing improved diagnosis of localized pulmonary diseases.

Contour propagation for lung tumor delineation in 4D-CT using tensor-product surface of uniform and non-uniform closed cubic B-splines

A robust contour propagation method is proposed to help physicians delineate lung tumors on all phase images of four-dimensional computed tomography (4D-CT) by only manually delineating the contours on a reference phase. The proposed method models the trajectory surface swept by a contour in a respiratory cycle as a tensor-product surface of two closed cubic B-spline curves: a non-uniform B-spline curve which models the contour and a uniform B-spline curve which models the trajectory of a point on the contour. The surface is treated as a deformable entity, and is optimized from an initial surface by moving its control vertices such that the sum of the intensity similarities between the sampling points on the manually delineated contour and their corresponding ones on different phases is maximized. The initial surface is constructed by fitting the manually delineated contour on the reference phase with a closed B-spline curve. In this way, the proposed method can focus the registration on the contour instead of the entire image to prevent the deformation of the contour from being smoothed by its surrounding tissues, and greatly reduce the time consumption while keeping the accuracy of the contour propagation as well as the temporal consistency of the estimated respiratory motions across all phases in 4D-CT. Eighteen 4D-CT cases with 235 gross tumor volume (GTV) contours on the maximal inhale phase and 209 GTV contours on the maximal exhale phase are manually delineated slice by slice. The maximal inhale phase is used as the reference phase, which provides the initial contours. On the maximal exhale phase, the Jaccard similarity coefficient between the propagated GTV and the manually delineated GTV is 0.881 [Formula: see text] 0.026, and the Hausdorff distance is 3.07 [Formula: see text] 1.08 mm. The time for propagating the GTV to all phases is 5.55 [Formula: see text] 6.21 min. The results are better than those of the fast adaptive stochastic gradient descent B-spline method, the 3D  +  t B-spline method and the diffeomorphic demons method. The proposed method is useful for helping physicians delineate target volumes efficiently and accurately.

New methods for the geometrical analysis of tubular organs

This paper presents new methods to study the shape of tubular organs. Determining precise cross-sections is of major importance to perform geometrical measurements, such as diameter, wall-thickness estimation or area measurement. Our first contribution is a robust method to estimate orthogonal planes based on the Voronoi Covariance Measure. Our method is not relying on a curve-skeleton computation beforehand. This means our orthogonal plane estimator can be used either on the skeleton or on the volume. Another important step towards tubular organ characterization is achieved through curve-skeletonization, as skeletons allow to compare two tubular organs, and to perform virtual endoscopy. Our second contribution is dedicated to correcting common defects of the skeleton by new pruning and recentering methods. Finally, we propose a new method for curve-skeleton extraction. Various results are shown on different types of segmented tubular organs, such as neurons, airway-tree and blood vessels.

Use of image registration and fusion algorithms and techniques in radiotherapy: Report of the AAPM Radiation Therapy Committee Task Group No. 132

Image registration and fusion algorithms exist in almost every software system that creates or uses images in radiotherapy. Most treatment planning systems support some form of image registration and fusion to allow the use of multimodality and time-series image data and even anatomical atlases to assist in target volume and normal tissue delineation. Treatment delivery systems perform registration and fusion between the planning images and the in-room images acquired during the treatment to assist patient positioning. Advanced applications are beginning to support daily dose assessment and enable adaptive radiotherapy using image registration and fusion to propagate contours and accumulate dose between image data taken over the course of therapy to provide up-to-date estimates of anatomical changes and delivered dose. This information aids in the detection of anatomical and functional changes that might elicit changes in the treatment plan or prescription. As the output of the image registration process is always used as the input of another process for planning or delivery, it is important to understand and communicate the uncertainty associated with the software in general and the result of a specific registration. Unfortunately, there is no standard mathematical formalism to perform this for real-world situations where noise, distortion, and complex anatomical variations can occur. Validation of the software systems performance is also complicated by the lack of documentation available from commercial systems leading to use of these systems in undesirable 'black-box' fashion. In view of this situation and the central role that image registration and fusion play in treatment planning and delivery, the Therapy Physics Committee of the American Association of Physicists in Medicine commissioned Task Group 132 to review current approaches and solutions for image registration (both rigid and deformable) in radiotherapy and to provide recommendations for quality assurance and quality control of these clinical processes.

Lung motion estimation by robust point matching and spatiotemporal tracking for 4D CT

We propose a deformable registration approach to estimate patient-specific lung motion during free breathing for four-dimensional (4D) computed tomography (CT) based on point matching and tracking between images in different phases. First, a robust point matching (RPM) algorithm coarsely aligns the source phase image onto all other target phase images of 4D CT. Scale-invariant feature transform (SIFT) is introduced into the cost function in order to accelerate and stabilize the convergence of the point matching. Next, the temporal consistency of the estimated lung motion model is preserved by fitting the trajectories of the points in the respiratory phase using L₁ norm regularization. Then, the fitted positions of a point along the trajectory are used as the initial positions for the point tracking. Spatial mean-shift iteration is employed to track points in all phase images. The tracked positions in all phases are used to perform RPM again. These steps are repeated until the number of updated points is smaller than a given threshold σ. With this method, the correspondence between the source phase image and other target phase image is established more accurately. Trajectory fitting ensures the estimated trajectory does not fluctuate violently. We evaluated our method by using the public DIR-lab, POPI-model, CREATIS and COPDgene lung datasets. In the experimental results, the proposed method achieved satisfied accuracy for image registration. Our method also preserved the topology of the deformation fields well for image registration with large deformation.

Sequential Monte Carlo tracking of the marginal artery by multiple cue fusion and random forest regression

Given the potential importance of marginal artery localization in automated registration in computed tomography colonography (CTC), we have devised a semi-automated method of marginal vessel detection employing sequential Monte Carlo tracking (also known as particle filtering tracking) by multiple cue fusion based on intensity, vesselness, organ detection, and minimum spanning tree information for poorly enhanced vessel segments. We then employed a random forest algorithm for intelligent cue fusion and decision making which achieved high sensitivity and robustness. After applying a vessel pruning procedure to the tracking results, we achieved statistically significantly improved precision compared to a baseline Hessian detection method (2.7% versus 75.2%, p<0.001). This method also showed statistically significantly improved recall rate compared to a 2-cue baseline method using fewer vessel cues (30.7% versus 67.7%, p<0.001). These results demonstrate that marginal artery localization on CTC is feasible by combining a discriminative classifier (i.e., random forest) with a sequential Monte Carlo tracking mechanism. In so doing, we present the effective application of an anatomical probability map to vessel pruning as well as a supplementary spatial coordinate system for colonic segmentation and registration when this task has been confounded by colon lumen collapse.