Characterization and identification of spatial artifacts during 4D-CT imaging

4DCT image artifact detection using deep learning

BACKGROUND

Four-dimensional computed tomography (4DCT) is an es sential tool in radiation therapy. However, the 4D acquisition process may cause motion artifacts which can obscure anatomy and distort functional measurements from CT scans.

PURPOSE

We describe a deep learning algorithm to identify the location of artifacts within 4DCT images. Our method is flexible enough to handle different types of artifacts, including duplication, misalignment, truncation, and interpolation.

METHODS

We trained and validated a U-net convolutional neural network artifact detection model on more than 23 000 coronal slices extracted from 98 4DCT scans. The receiver operating characteristic (ROC) curve and precision-recall curve were used to evaluate the model's performance at identifying artifacts compared to a manually identified ground truth. The model was adjusted so that the sensitivity in identifying artifacts was equivalent to that of a human observer, as measured by computing the average ratio of artifact volume to lung volume in a given scan.

RESULTS

The model achieved a sensitivity, specificity, and precision of 0.78, 0.99, and 0.58, respectively. The ROC area-under-the-curve (AUC) was 0.99 and the precision-recall AUC was 0.73. Our model sensitivity is 8% higher than previously reported state-of-the-art artifact detection methods.

CONCLUSIONS

The model developed in this study is versatile, designed to handle duplication, misalignment, truncation, and interpolation artifacts within a single image, unlike earlier models that were designed for a single artifact type.

Robust quantification of CT-ventilation biomarker techniques and repeatability in a porcine model

BACKGROUND

Biomarkers estimating local lung ventilation have been derived from computed tomography (CT) imaging using various image acquisition and post-processing techniques. CT-ventilation biomarkers have potential clinical use in functional avoidance radiation therapy (RT), in which RT treatment plans are optimized to reduce dose delivered to highly ventilated lung. Widespread clinical implementation of CT-ventilation biomarkers necessitates understanding of biomarker repeatability. Performing imaging within a highly controlled experimental design enables quantification of error associated with remaining variables.

PURPOSE

To characterize CT-ventilation biomarker repeatability and dependence on image acquisition and post-processing methodology in anesthetized and mechanically ventilated pigs.

METHODS

Five mechanically ventilated Wisconsin Miniature Swine (WMS) received multiple consecutive four-dimensional CT (4DCT) and maximum inhale and exhale breath-hold CT (BH-CT) scans on five dates to generate CT-ventilation biomarkers. Breathing maneuvers were controlled with an average tidal volume difference <200 cc. As surrogates for ventilation, multiple local expansion ratios (LERs) were calculated from the acquired CT scans using Jacobian-based post-processing techniques.L E R 2 $LER_2$ measured local expansion between an image pair using either inhale and exhale BH-CT images or two 4DCT breathing phase images.L E R N $LER_N$ measured the maximum local expansion across the 4DCT breathing phase images. Breathing maneuver consistency, intra- and interday biomarker repeatability, image acquisition and post-processing technique dependence were quantitatively analyzed.

RESULTS

Biomarkers showed strong agreement with voxel-wise Spearman correlationρ > 0.9 $\rho > 0.9$ for intraday repeatability andρ > 0.8 $\rho > 0.8$ for all other comparisons, including between image acquisition techniques. Intra- and interday repeatability were significantly different (p < 0.01). LER2 and LERN post-processing did not significantly affect intraday repeatability.

CONCLUSIONS

4DCT and BH-CT ventilation biomarkers derived from consecutive scans show strong agreement in controlled experiments with nonhuman subjects.

Metrics of dose to highly ventilated lung are predictive of radiation-induced pneumonitis in lung cancer patients

PURPOSE

To identify metrics of radiation dose delivered to highly ventilated lung that are predictive of radiation-induced pneumonitis.

METHODS AND MATERIALS

A cohort of 90 patients with locally advanced non-small cell lung cancer treated with standard fractionated radiation therapy (RT) (60-66 Gy in 30-33 fractions) were evaluated. Regional lung ventilation was determined from pre-RT 4-dimensional computed tomography (4DCT) using the Jacobian determinant of a B-spline deformable image registration to estimate lung tissue expansion during respiration. Multiple voxel-wise population- and individual-based thresholds for defining high functioning lung were considered. Mean dose and volumes receiving dose ≥ 5-60 Gy were analyzed for both total lung-ITV (MLD,V5-V60) and highly ventilated functional lung-ITV (fMLD,fV5-fV60). The primary endpoint was symptomatic grade 2+ (G2+) pneumonitis. Receiver operator curve (ROC) analyses were used to identify predictors of pneumonitis.

RESULTS

G2+ pneumonitis occurred in 22.2% of patients, with no differences between stage, smoking status, COPD, or chemo/immunotherapy use between G<2 and G2+ patients (P≥ 0.18). Highly ventilated lung was defined as voxels exceeding the population-wide median of 18% voxel-level expansion. All total and functional metrics were significantly different between patients with and without pneumonitis (P≤ 0.039). Optimal ROC points predicting pneumonitis from functional lung dose were fMLD ≤ 12.3 Gy, fV5 ≤ 54% and fV20 ≤ 19 %. Patients with fMLD ≤ 12.3 Gy had a 14% risk of developing G2+ pneumonitis whereas risk significantly increased to 35% for those with fMLD > 12.3 Gy (P = 0.035).

CONCLUSIONS

Dose to highly ventilated lung is associated with symptomatic pneumonitis and treatment planning strategies should focus on limiting dose to functional regions. These findings provide important metrics to be used in functional lung avoidance RT planning and designing clinical trials.

Quantifying robustness of CT-ventilation biomarkers to image noise

Purpose: To quantify the impact of image noise on CT-based lung ventilation biomarkers calculated using Jacobian determinant techniques. Methods: Five mechanically ventilated swine were imaged on a multi-row CT scanner with acquisition parameters of 120 kVp and 0.6 mm slice thickness in static and 4-dimensional CT (4DCT) modes with respective pitches of 1 and 0.09. A range of tube current time product (mAs) values were used to vary image dose. On two dates, subjects received two 4DCTs: one with 10 mAs/rotation (low-dose, high-noise) and one with CT simulation standard of care 100 mAs/rotation (high-dose, low-noise). Additionally, 10 intermediate noise level breath-hold (BHCT) scans were acquired with inspiratory and expiratory lung volumes. Images were reconstructed with and without iterative reconstruction (IR) using 1 mm slice thickness. The Jacobian determinant of an estimated transformation from a B-spline deformable image registration was used to create CT-ventilation biomarkers estimating lung tissue expansion. 24 CT-ventilation maps were generated per subject per scan date: four 4DCT ventilation maps (two noise levels each with and without IR) and 20 BHCT ventilation maps (10 noise levels each with and without IR). Biomarkers derived from reduced dose scans were registered to the reference full dose scan for comparison. Evaluation metrics were gamma pass rate (Γ) with 2 mm distance-to-agreement and 6% intensity criterion, voxel-wise Spearman correlation (ρ) and Jacobian ratio coefficient of variation (CoV JR ). Results: Comparing biomarkers derived from low (CTDI vol = 6.07 mGy) and high (CTDI vol = 60.7 mGy) dose 4DCT scans, mean Γ, ρ and CoV JR values were 93% ± 3%, 0.88 ± 0.03 and 0.04 ± 0.009, respectively. With IR applied, those values were 93% ± 4%, 0.90 ± 0.04 and 0.03 ± 0.003. Similarly, comparisons between BHCT-based biomarkers with variable dose (CTDI vol = 1.35-7.95 mGy) had mean Γ, ρ and CoV JR of 93% ± 4%, 0.97 ± 0.02 and 0.03 ± 0.006 without IR and 93% ± 4%, 0.97 ± 0.03 and 0.03 ± 0.007 with IR. Applying IR did not significantly change any metrics (p > 0.05). Discussion: This work demonstrated that CT-ventilation, calculated using the Jacobian determinant of an estimated transformation from a B-spline deformable image registration, is invariant to Hounsfield Unit (HU) variation caused by image noise. This advantageous finding may be leveraged clinically with potential applications including dose reduction and/or acquiring repeated low-dose acquisitions for improved ventilation characterization.

Using 4D dose accumulation to calculate organ-at-risk dose deviations from motion-synchronized liver and lung tomotherapy treatments

Tracking systems such as Radixact Synchrony change the planned delivery of radiation during treatment to follow the target. This is typically achieved without considering the location changes of organs at risk (OARs). The goal of this work was to develop a novel 4D dose accumulation framework to quantify OAR dose deviations due to the motion and tracked treatment. The framework obtains deformation information and the target motion pattern from a four-dimensional computed tomography dataset. The helical tomotherapy treatment plan is split into 10 plans and motion correction is applied separately to the jaw pattern and multi-leaf collimator (MLC) sinogram for each phase based on the location of the target in each phase. Deformable image registration (DIR) is calculated from each phase to the references phase using a commercial algorithm, and doses are accumulated according to the DIR. The effect of motion synchronization on OAR dose was analyzed for five lung and five liver subjects by comparing planned versus synchrony-accumulated dose. The motion was compensated by an average of 1.6 cm of jaw sway and by an average of 5.7% of leaf openings modified, indicating that most of the motion compensation was from jaw sway and not MLC changes. OAR dose deviations as large as 19 Gy were observed, and for all 10 cases, dose deviations greater than 7 Gy were observed. Target dose remained relatively constant (D95% within 3 Gy), confirming that motion-synchronization achieved the goal of maintaining target dose. Dose deviations provided by the framework can be leveraged during the treatment planning process by identifying cases where OAR doses may change significantly from their planned values with respect to the critical constraints. The framework is specific to synchronized helical tomotherapy treatments, but the OAR dose deviations apply to any real-time tracking technique that does not consider location changes of OARs.

On the evaluation of mobile target trajectory between four-dimensional computer tomography and four-dimensional cone-beam computer tomography

Purpose

For mobile lung tumors, four‐dimensional computer tomography (4D CT) is often used for simulation and treatment planning. Localization accuracy remains a challenge in lung stereotactic body radiation therapy (SBRT) treatments. An attractive image guidance method to increase localization accuracy is 4D cone‐beam CT (CBCT) as it allows for visualization of tumor motion with reduced motion artifacts. However, acquisition and reconstruction of 4D CBCT differ from that of 4D CT. This study evaluates the discrepancies between the reconstructed motion of 4D CBCT and 4D CT imaging over a wide range of sine target motion parameters and patient waveforms.

Methods

A thorax motion phantom was used to examine 24 sine motions with varying amplitudes and cycle times and seven patient waveforms. Each programmed motion was imaged using 4D CT and 4D CBCT. The images were processed to auto segment the target. For sine motion, the target centroid at each phase was fitted to a sinusoidal curve to evaluate equivalence in amplitude between the two imaging modalities. The patient waveform motion was evaluated based on the average 4D data sets.

Results

The mean difference and root‐mean‐square‐error between the two modalities for sine motion were −0.35 ± 0.22 and 0.60 mm, respectively, with 4D CBCT slightly overestimating amplitude compared with 4D CT. The two imaging methods were determined to be significantly equivalent within ±1 mm based on two one‐sided t tests (p < 0.001). For patient‐specific motion, the mean difference was 1.5 ± 2.1 (0.8 ± 0.6 without outlier), 0.4 ± 0.3, and 0.8 ± 0.6 mm for superior/inferior (SI), anterior/posterior (AP), and left/right (LR), respectively.

Conclusion

In cases where 4D CT is used to image mobile tumors, 4D CBCT is an attractive localization method due to its assessment of motion with respect to 4D CT, particularly for lung SBRT treatments where accuracy is paramount.

Variation in target volume and centroid position due to breath holding during four-dimensional computed tomography scanning: A phantom study

This study investigated the effects of respiratory motion, including unwanted breath holding, on the target volume and centroid position on four‐dimensional computed tomography (4DCT) imaging. Cine 4DCT images were reconstructed based on a time‐based sorting algorithm, and helical 4DCT images were reconstructed based on both the time‐based sorting algorithm and an amplitude‐based sorting algorithm. A spherical object 20 mm in diameter was moved according to several simulated respiratory motions, with a motion period of 4.0 s and maximum amplitude of 5 mm. The object was extracted automatically, and the target volume and centroid position in the craniocaudal direction were measured using a treatment planning system. When the respiratory motion included unwanted breath‐holding times shorter than the breathing cycle, the root mean square errors (RSME) between the reference and imaged target volumes were 18.8%, 14.0%, and 5.5% in time‐based images in cine mode, time‐based images in helical mode, and amplitude‐based images in helical mode, respectively. In helical mode, the RSME between the reference and imaged centroid position was reduced from 1.42 to 0.50 mm by changing the reconstruction method from time‐ to amplitude‐based sorting. When the respiratory motion included unwanted breath‐holding times equal to the breathing cycle, the RSME between the reference and imaged target volumes were 19.1%, 24.3%, and 15.6% in time‐based images in cine mode, time‐based images in helical mode, and amplitude‐based images in helical mode, respectively. In helical mode, the RSME between the reference and imaged centroid position was reduced from 1.61 to 0.83 mm by changing the reconstruction method from time‐ to amplitude‐based sorting. With respiratory motion including breath holding of shorter duration than the breathing cycle, the accuracies of the target volume and centroid position were improved by amplitude‐based sorting, particularly in helical 4DCT.

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Background

Motion due to respiration may challenge dose prediction of dose calculation algorithms during treatment planning.

Materials and methods

The two versions of depth dose calculation algorithms, namely, Anisotropic Analytical Algorithm (AAA) version 10.0 (AAAv10.0), AAA version 13.6 (AAAv13.6) and Acuros XB dose calculation (AXB) algorithm version 10.0 (AXBv10.0), AXB version 13.6 (AXBv13.6), were compared against a full MC simulated 6X photon beam using QUASAR respiratory motion phantom with a moving chest wall. To simulate the moving chest wall, a 4 cm thick wax mould was attached to the lung insert of the phantom. Depth doses along the central axis were compared in the anterior and lateral beam direction for field sizes 2 × 2 cm², 4 × 4 cm² and 10 × 10 cm².

Results

For the lateral beam direction, the moving chest wall highlighted differences of up to 105% for AAAv10.0 and 40% for AXBv10.0 from MC calculations in the surface and buildup doses. AAAv13.6 and AXBv13.6 agrees with MC predictions to within 10% at similar depth. For anterior beam doses, dose differences predicted for both versions of AAA and AXB algorithm were within 7% and results were consistent with static heterogeneous studies.

Conclusions

The presence of the moving chest wall was capable of identifying depth dose differences between the two versions of the algorithms. These differences could not be identified in the static chest wall as shown in the anterior beam depth dose calculations.

Deformable image registration applied to lung SBRT: Usefulness and limitations

Radiation therapy (RT) of the lung requires deformation analysis. Deformable image registration (DIR) is the fundamental method to quantify deformations for various applications: motion compensation, contour propagation, dose accumulation, etc. DIR is therefore unavoidable in lung RT. DIR algorithms have been studied for decades and are now available both within commercial and academic packages. However, they are complex and have limitations that every user must be aware of before clinical implementation. In this paper, the main applications of DIR for lung RT with their associated uncertainties and their limitations are reviewed.

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.

Geometric validation of self-gating k-space-sorted 4D-MRI vs 4D-CT using a respiratory motion phantom

PURPOSE

MRI is increasingly being used for radiotherapy planning, simulation, and in-treatment-room motion monitoring. To provide more detailed temporal and spatial MR data for these tasks, we have recently developed a novel self-gated (SG) MRI technique with advantage of k-space phase sorting, high isotropic spatial resolution, and high temporal resolution. The current work describes the validation of this 4D-MRI technique using a MRI- and CT-compatible respiratory motion phantom and comparison to 4D-CT.

METHODS

The 4D-MRI sequence is based on a spoiled gradient echo-based 3D projection reconstruction sequence with self-gating for 4D-MRI at 3 T. Respiratory phase is resolved by using SG k-space lines as the motion surrogate. 4D-MRI images are reconstructed into ten temporal bins with spatial resolution 1.56 × 1.56 × 1.56 mm(3). A MRI-CT compatible phantom was designed to validate the performance of the 4D-MRI sequence and 4D-CT imaging. A spherical target (diameter 23 mm, volume 6.37 ml) filled with high-concentration gadolinium (Gd) gel is embedded into a plastic box (35 × 40 × 63 mm(3)) and stabilized with low-concentration Gd gel. The phantom, driven by an air pump, is able to produce human-type breathing patterns between 4 and 30 respiratory cycles/min. 4D-CT of the phantom has been acquired in cine mode, and reconstructed into ten phases with slice thickness 1.25 mm. The 4D images sets were imported into a treatment planning software for target contouring. The geometrical accuracy of the 4D MRI and CT images has been quantified using target volume, flattening, and eccentricity. The target motion was measured by tracking the centroids of the spheres in each individual phase. Motion ground-truth was obtained from input signals and real-time video recordings.

RESULTS

The dynamic phantom has been operated in four respiratory rate (RR) settings, 6, 10, 15, and 20/min, and was scanned with 4D-MRI and 4D-CT. 4D-CT images have target-stretching, partial-missing, and other motion artifacts in various phases, whereas the 4D-MRI images are visually free of those artifacts. Volume percentage difference for the 6.37 ml target ranged from 5.3% ± 4.3% to 10.3% ± 5.9% for 4D-CT, and 1.47 ± 0.52 to 2.12 ± 1.60 for 4D-MRI. With an increase of respiratory rate, the target volumetric and geometric deviations increase for 4D-CT images while remaining stable for the 4D-MRI images. Target motion amplitude errors at different RRs were measured with a range of 0.66-1.25 mm for 4D-CT and 0.2-0.42 mm for 4D-MRI. The results of Mann-Whitney tests indicated that 4D-MRI significantly outperforms 4D-CT in phase-based target volumetric (p = 0.027) and geometric (p < 0.001) measures. Both modalities achieve equivalent accuracy in measuring motion amplitude (p = 0.828).

CONCLUSIONS

The k-space self-gated 4D-MRI technique provides a robust method for accurately imaging phase-based target motion and geometry. Compared to 4D-CT, the current 4D-MRI technique demonstrates superior spatiotemporal resolution, and robust resistance to motion artifacts caused by fast target motion and irregular breathing patterns. The technique can be used extensively in abdominal targeting, motion gating, and toward implementing MRI-based adaptive radiotherapy.

Maximum-Likelihood Joint Image Reconstruction/Motion Estimation in Attenuation-Corrected Respiratory Gated PET/CT Using a Single Attenuation Map

This work provides an insight into positron emission tomography (PET) joint image reconstruction/motion estimation (JRM) by maximization of the likelihood, where the probabilistic model accounts for warped attenuation. Our analysis shows that maximum-likelihood (ML) JRM returns the same reconstructed gates for any attenuation map (μ-map) that is a deformation of a given μ-map, regardless of its alignment with the PET gates. We derived a joint optimization algorithm accordingly, and applied it to simulated and patient gated PET data. We first evaluated the proposed algorithm on simulations of respiratory gated PET/CT data based on the XCAT phantom. Our results show that independently of which μ-map is used as input to JRM: (i) the warped μ-maps correspond to the gated μ-maps, (ii) JRM outperforms the traditional post-registration reconstruction and consolidation (PRRC) for hot lesion quantification and (iii) reconstructed gated PET images are similar to those obtained with gated μ-maps. This suggests that a breath-held μ-map can be used. We then applied JRM on patient data with a μ-map derived from a breath-held high resolution CT (HRCT), and compared the results with PRRC, where each reconstructed PET image was obtained with a corresponding cine-CT gated μ-map. Results show that JRM with breath-held HRCT achieves similar reconstruction to that using PRRC with cine-CT. This suggests a practical low-dose solution for implementation of motion-corrected respiratory gated PET/CT.

Respiratory effort correction strategies to improve the reproducibility of lung expansion measurements

PURPOSE

Four-dimensional computed tomography (4DCT) can be used to make measurements of pulmonary function longitudinally. The sensitivity of such measurements to identify change depends on measurement uncertainty. Previously, intrasubject reproducibility of Jacobian-based measures of lung tissue expansion was studied in two repeat prior-RT 4DCT human acquisitions. Difference in respiratory effort such as breathing amplitude and frequency may affect longitudinal function assessment. In this study, the authors present normalization schemes that correct ventilation images for variations in respiratory effort and assess the reproducibility improvement after effort correction.

METHODS

Repeat 4DCT image data acquired within a short time interval from 24 patients prior to radiation therapy (RT) were used for this analysis. Using a tissue volume preserving deformable image registration algorithm, Jacobian ventilation maps in two scanning sessions were computed and compared on the same coordinate for reproducibility analysis. In addition to computing the ventilation maps from end expiration to end inspiration, the authors investigated the effort normalization strategies using other intermediated inspiration phases upon the principles of equivalent tidal volume (ETV) and equivalent lung volume (ELV). Scatter plots and mean square error of the repeat ventilation maps and the Jacobian ratio map were generated for four conditions: no effort correction, global normalization, ETV, and ELV. In addition, gamma pass rate was calculated from a modified gamma index evaluation between two ventilation maps, using acceptance criterions of 2 mm distance-to-agreement and 5% ventilation difference.

RESULTS

The pattern of regional pulmonary ventilation changes as lung volume changes. All effort correction strategies improved reproducibility when changes in respiratory effort were greater than 150 cc (p < 0.005 with regard to the gamma pass rate). Improvement of reproducibility was correlated with respiratory effort difference (R = 0.744 for ELV in the cohort with tidal volume difference greater than 100 cc). In general for all subjects, global normalization, ETV and ELV significantly improved reproducibility compared to no effort correction (p = 0.009, 0.002, 0.005 respectively). When tidal volume difference was small (less than 100 cc), none of the three effort correction strategies improved reproducibility significantly (p = 0.52, 0.46, 0.46 respectively). For the cohort (N = 13) with tidal volume difference greater than 100 cc, the average gamma pass rate improves from 57.3% before correction to 66.3% after global normalization, and 76.3% after ELV. ELV was found to be significantly better than global normalization (p = 0.04 for all subjects, and p = 0.003 for the cohort with tidal volume difference greater than 100 cc).

CONCLUSIONS

All effort correction strategies improve the reproducibility of the authors' pulmonary ventilation measures, and the improvement of reproducibility is highly correlated with the changes in respiratory effort. ELV gives better results as effort difference increase, followed by ETV, then global. However, based on the spatial and temporal heterogeneity in the lung expansion rate, a single scaling factor (e.g., global normalization) appears to be less accurate to correct the ventilation map when changes in respiratory effort are large.

Assessment of a quantitative metric for 4D CT artifact evaluation by observer consensus

The benefits of four-dimensional computed tomography (4D CT) are limited by the presence of artifacts that remain difficult to quantify. A correlation-based metric previously proposed for ciné 4D CT artifact identification was further validated as an independent artifact evaluator by using a novel qualitative assessment featuring a group of observers reaching a consensus decision on artifact location and magnitude. The consensus group evaluated ten ciné 4D CT scans for artifacts over each breathing phase of coronal lung views assuming one artifact per couch location. Each artifact was assigned a magnitude score of 1-5, 1 indicating lowest severity and 5 indicating highest severity. Consensus group results served as the ground truth for assessment of the correlation metric. The ten patients were split into two cohorts; cohort 1 generated an artifact identification threshold derived from receiver operating characteristic analysis using the Youden Index, while cohort 2 generated sensitivity and specificity values from application of the artifact threshold. The Pearson correlation coefficient was calculated between the correlation metric values and the consensus group scores for both cohorts. The average sensitivity and specificity values found with application of the artifact threshold were 0.703 and 0.476, respectively. The correlation coefficients of artifact magnitudes for cohort 1 and 2 were 0.80 and 0.61, respectively, (p < 0.001 for both); these correlation coefficients included a few scans with only two of the five possible magnitude scores. Artifact incidence was associated with breathing phase (p < 0.002), with presentation less likely near maximum exhale. Overall, the correlation metric allowed accurate and automated artifact identification. The consensus group evaluation resulted in efficient qualitative scoring, reduced interobserver variation, and provided consistent identification of artifact location and magnitudes.

Estimating the 4D respiratory lung motion by spatiotemporal registration and super-resolution image reconstruction

PURPOSE

One of the main challenges in lung cancer radiation therapy is how to reduce the treatment margin but accommodate the geometric uncertainty of moving tumor. 4D-CT is able to provide the full range of motion information for the lung and tumor. However, accurate estimation of lung motion with respect to the respiratory phase is difficult due to various challenges in image registration, e.g., motion artifacts and large interslice thickness in 4D-CT. Meanwhile, the temporal coherence across respiration phases is usually not guaranteed in the conventional registration methods which consider each phase image in 4D-CT independently. To address these challenges, the authors present a unified approach to estimate the respiratory lung motion with two iterative steps.

METHODS

First, the authors propose a novel spatiotemporal registration algorithm to align all phase images of 4D-CT (in low-resolution) to a high-resolution group-mean image in the common space. The temporal coherence of registration is maintained by a set of temporal fibers that delineate temporal correspondences across different respiratory phases. Second, a super-resolution technique is utilized to build the high-resolution group-mean image with more anatomical details than any individual phase image, thus largely alleviating the registration uncertainty especially in correspondence detection. In particular, the authors use the concept of sparse representation to keep the group-mean image as sharp as possible.

RESULTS

The performance of our 4D motion estimation method has been extensively evaluated on both the simulated datasets and real lung 4D-CT datasets. In all experiments, our method achieves more accurate and consistent results in lung motion estimation than all other state-of-the-art approaches under comparison.

CONCLUSIONS

The authors have proposed a novel spatiotemporal registration method to estimate the lung motion in 4D-CT. Promising results have been obtained, which indicates the high applicability of our method in clinical lung cancer radiation therapy.

Reproducibility of registration-based measures of lung tissue expansion

PURPOSE

Lung function depends on lung expansion and contraction during the respiratory cycle. Respiratory-gated CT imaging and 3D image registration can be used to locally estimate lung tissue expansion and contraction (regional lung volume change) by computing the determinant of the Jacobian matrix of the image registration deformation field. In this study, the authors examine the reproducibility of Jacobian-based measures of lung tissue expansion in two repeat 4DCT acquisitions of mechanically ventilated sheep and free-breathing humans.

METHODS

4DCT image data from three white sheep and nine human subjects were used for this analysis. In each case, two 4DCT studies were acquired for each subject within a short time interval. The animal subjects were anesthetized and mechanically ventilated, while the humans were awake and spontaneously breathing based on respiratory pacing audio cues. From each 4DCT data set, an image pair consisting of a volume reconstructed near end inspiration and a volume reconstructed near end exhalation was selected. The end inspiration and end exhalation images were registered using a tissue volume preserving deformable registration algorithm and the Jacobian of the registration deformation field was used to measure regional lung expansion. The Jacobian map from the baseline data set was compared to the Jacobian map from the followup data by measuring the voxel-by-voxel Jacobian ratio.

RESULTS

In the animal subjects, the mean Jacobian ratio (baseline scan Jacobian divided by followup scan Jacobian, voxel-by-voxel) was 0.9984±0.021 (mean ± standard deviation, averaged over the entire lung region). The mean Jacobian ratio was 1.0224±0.058 in the human subjects. The reproducibility of the Jacobian values was found to be strongly dependent on the reproducibility of the subject's respiratory effort and breathing pattern.

CONCLUSIONS

Lung expansion, a surrogate for lung function, can be assessed using two or more respiratory-gated CT image acquisitions. The results show that good reproducibility can be obtained in anesthetized, mechanically ventilated animals, but variations in respiratory effort and breathing patterns reduce reproducibility in spontaneously-breathing humans. The global linear normalization can globally compensate for breathing effort differences, but a homogeneous scaling does not account for differences in regional lung expansion rates. Additional work is needed to develop compensation procedures or normalization schemes that can account for local variations in lung expansion during respiration.

Feature Guided Motion Artifact Reduction with Structure-Awareness in 4D CT Images

In this paper, we propose a novel method to reduce the magnitude of 4D CT artifacts by stitching two images with a data-driven regularization constrain, which helps preserve the local anatomy structures. Our method first computes an interface seam for the stitching in the overlapping region of the first image, which passes through the "smoothest" region, to reduce the structure complexity along the stitching interface. Then, we compute the displacements of the seam by matching the corresponding interface seam in the second image. We use sparse 3D features as the structure cues to guide the seam matching, in which a regularization term is incorporated to keep the structure consistency. The energy function is minimized by solving a multiple-label problem in Markov Random Fields with an anatomical structure preserving regularization term. The displacements are propagated to the rest of second image and the two image are stitched along the interface seams based on the computed displacement field. The method was tested on both simulated data and clinical 4D CT images. The experiments on simulated data demonstrated that the proposed method was able to reduce the landmark distance error on average from 2.9 mm to 1.3 mm, outperforming the registration-based method by about 55%. For clinical 4D CT image data, the image quality was evaluated by three medical experts, and all identified much fewer artifacts from the resulting images by our method than from those by the compared method.