Kinematic Analysis of the Distal Radioulnar Joint in Asymptomatic Wrists Using 4-Dimensional Computed Tomography-Motion Pattern and Interreader Reliability.

PURPOSE: The aim of this study was to determine the normal measurement values and interobserver performance of the distal radioulnar joint during wrist pronation-supination using 4-dimensional computed tomography (CT). METHODS: Four-dimensional CT examinations were performed on the asymptomatic contralateral wrists of 10 patients with unilateral chronic wrist pain. Measurements were conducted using the modified radioulnar (mRU) line and epicenter (Epi) methods. Volar subluxation of the ulnar head was demonstrated with negative values. Wilcoxon rank sum test was used to determine the measurement changes. Interobserver agreements were assessed using interclass correlation coefficients. RESULTS: In pronation, mRU line measurements (median, 0.09; interquartile range, 0-0.15) were significantly larger than in supination (median, -0.1; interquartile range, -0.18 to 0; P = 0.008).The Epi measurements were not significantly different in pronation (median, 0.03; interquartile range, 0.01-0.07) and supination (median, 0.06; interquartile range, 0.01-0.1; P = 0.799). There was an excellent inter-observer agreement between the two readers using mRU and Epi methods in pronation (0.982, 0.898), midpoint (0.994, 0.827) and supination (0.989, 0.972) positions, respectively. CONCLUSIONS: Using 4-dimensional CT examination, distal radioulnar joint kinematics in asymptomatic wrists demonstrate excellent interobserver agreements with increased volar ulnar subluxation with supination as detected using mRU, but not the Epi method.

Evidence-based recommendations for musculoskeletal kinematic 4D-CT studies using wide area-detector scanners: a phantom study with cadaveric correlation.

PROPOSE: To establish evidence-based recommendations for musculoskeletal kinematic 4D-CT on wide area-detector CT. MATERIALS AND METHODS: In order to assess factors influencing image quality in kinematic CT studies, a phantom consisting of a polymethylmethacrylate rotating disk with round wells of different sizes was imaged with various acquisition protocols. Cadaveric acquisitions were performed on the ankle joint during motion in two different axes and at different speeds to allow validation of phantom data. Images were acquired with a 320 detector-row CT scanner and were evaluated by two readers. RESULTS: Motion artefacts were significantly correlated with various parameters (movement axis, distance to centre, rotation speed and volume acquisition speed) (p < 0.0001). The relation between motion artefacts and distance to motion fulcrum was exponential (R2 0.99). Half reconstruction led to a 23 % increase in image noise and a 40 % decrease in motion artefacts. Cadaveric acquisitions confirmed phantom data. Based on these findings, high tube rotation speed and half reconstruction are recommended for kinematic CT. The axis of motion significantly influences image artefacts and should be considered in patient training and evaluation of acquisition protocol suitability. CONCLUSION: This study provides evidence-based recommendations for musculoskeletal kinematic 4D-CT. KEY POINTS: • Motion artefacts can hamper the quality and interpretation of dynamic joint studies • The recommendations presented here help increase image quality • Patient training and preparation can be improved • The artefact-free distance concept helps protocol adaptation and comparison.

Baseline correction of a correlation model for improving the prediction accuracy of infrared marker-based dynamic tumor tracking.

We previously found that the baseline drift of external and internal respiratory motion reduced the prediction accuracy of infrared (IR) marker-based dynamic tumor tracking irradiation (IR Tracking) using the Vero4DRT system. Here, we proposed a baseline correction method, applied immediately before beam delivery, to improve the prediction accuracy of IR Tracking. To perform IR Tracking, a four-dimensional (4D) model was constructed at the beginning of treatment to correlate the internal and external respiratory signals, and the model was expressed using a quadratic function involving the IR marker position (x) and its velocity (v), namely function F(x,v). First, the first 4D model, F1st(x,v), was adjusted by the baseline drift of IR markers (BDIR) along the x-axis, as function F'(x,v). Next, BDdetect, that defined as the difference between the target positions indicated by the implanted fiducial markers (Pdetect) and the predicted target positions with F'(x,v) (Ppredict) was determined using orthogonal kV X-ray images at the peaks of the Pdetect of the end-inhale and end-exhale phases for 10 s just before irradiation. F'(x,v) was corrected with BDdetect to compensate for the residual error. The final corrected 4D model was expressed as Fcor(x,v) = F1st{(x-BDIR),v}-BDdetect. We retrospectively applied this function to 53 paired log files of the 4D model for 12 lung cancer patients who underwent IR Tracking. The 95th percentile of the absolute differences between Pdetect and Ppredict (|Ep|) was compared between F1st(x,v) and Fcor(x,v). The median 95th percentile of |Ep| (units: mm) was 1.0, 1.7, and 3.5 for F1st(x,v), and 0.6, 1.1, and 2.1 for Fcor(x,v) in the left-right, anterior-posterior, and superior-inferior directions, respectively. Over all treatment sessions, the 95th percentile of |Ep| peaked at 3.2 mm using Fcor(x,v) compared with 8.4 mm using F1st(x,v). Our proposed method improved the prediction accuracy of IR Tracking by correcting the baseline drift immediately before irradiation.

Intrafractional prostate motion during external beam radiotherapy monitored by a real-time target localization system.

This paper investigates the clinical significance of real-time monitoring of intrafractional prostate motion during external beam radiotherapy using a commercial 4D localization system. Intrafractional prostate motion was tracked during 8,660 treatment fractions for 236 patients. The following statistics were analyzed: 1) the percentage of fractions in which the prostate shifted 2-7 mm for a certain duration; 2) the proportion of the entire tracking time during which the prostate shifted 2-7mm; and 3) the proportion of each minute in which the shift exceeded 2-7 mm. The ten patients exhibiting maximum intrafractional-motion patterns were analyzed separately. Our results showed that the percentage of fractions in which the prostate shifted by > 2, 3, 5, and 7 mm off the baseline in any direction for > 30 s was 56.8%, 27.2%, 4.6%, and 0.7% for intact prostate and 68.7%, 35.6%, 10.1%, and 1.8% for postprostatectomy patients, respectively. For the ten patients, these percentages were 91.3%, 72.4%, 36.3%, and 6%, respectively. The percentage of tracking time during which the prostate shifted > 2, 3, 5, and 7 mm was 27.8%, 10.7%, 1.6%, and 0.3%, respectively, and it was 56.2%, 33.7%, 11.2%, and 2.1%, respectively, for the ten patients. The percentage of tracking time for a > 3 mm posterior motion was four to five times higher than that in other directions. For treatments completed in 5 min (VMAT) and 10 min (IMRT), the proportion for the prostate to shift by > 3mm was 4% and 12%, respectively. Although intrafractional prostate motion was generally small, caution should be taken for patients who exhibit frequent large intrafractional motion. For those patients, adjustment of patient positioning may be necessary or a larger treatment margin may be used. After the initial alignment, the likelihood of prostate motion increases with time. Therefore, it is favorable to use advanced techniques (e.g., VMAT) that require less delivery time in order to reduce the treatment uncertainty resulting from intrafractional prostate motion.

Combination effects of tissue heterogeneity and geometric targeting error in stereotactic body radiotherapy for lung cancer using CyberKnife.

We have investigated the combined effect of tissue heterogeneity and its variation associated with geometric error in stereotactic body radiotherapy (SBRT) for lung cancer. The treatment plans for eight lung cancer patients were calculated using effective path length (EPL) correction and Monte Carlo (MC) algorithms, with both having the same beam configuration for each patient. These two kinds of plans for individual patients were then subsequently recalculated with adding systematic and random geometric errors. In the ordinary treatment plans calculated with no geometric offset, the EPL calculations, compared with the MC calculations, largely overestimated the doses to PTV by ~ 21%, whereas the overestimation were markedly lower in GTV by ~ 12% due to relatively higher density of GTV than of PTV. When recalculating the plans for individual patients with assigning the systematic and random geometric errors, no significant changes in the relative dose distribution, except for overall shift, were observed in the EPL calculations, whereas largely altered in the MC calculations with a consistent increase in dose to GTV. Considering the better accuracy of MC than EPL algorithms, the present results demonstrated the strong coupling of tissue heterogeneity and geometric error, thereby emphasizing the essential need for simultaneous correction for tissue heterogeneity and geometric targeting error in SBRT of lung cancer.

Deformation vector fields (DVF)-driven image reconstruction for 4D-CBCT.

BACKGROUND: High quality 4D-CBCT can be obtained by deforming a planning CT (pCT), where the deformation vector fields (DVF) are estimated by matching the forward projections of pCT and 4D-CBCT projections. The matching metric used in the previous study is the sum of squared intensity differences (SSID). The scatter signal level in CBCT projections is much higher than pCT, the SSID metric may not lead to optimal DVF. OBJECTIVE: To improve the DVF estimation accuracy, we develop a new matching metric that is less sensitive to the intensity level difference caused by the scatter signal. METHODS: The negative logarithm of correlation coefficient (NLCC) is used as the matching metric. A non-linear conjugate gradient optimization algorithm is used to estimate the DVF. A 4D NCAT phantom and an anthropomorphic thoracic phantom were used to evaluate the NLCC-based algorithm. RESULTS: In the NCAT phantom study, the relative reconstruction error is reduced from 18.0% in SSID to 14.13% in NLCC. In the thoracic phantom study, the root mean square error of the tumor motion is reduced from 1.16 mm in SSID to 0.43 mm in NLCC. CONCLUSION: NLCC metric can improve the image reconstruction and motion estimation accuracy of DVF-driven image reconstruction for 4D-CBCT.

Length of occlusion predicts recanalization and outcome after intravenous thrombolysis in middle cerebral artery stroke.

BACKGROUND AND PURPOSE: The length of large vessel occlusion is considered a major factor for therapy in patients with ischemic stroke. We used 4D-CT angiography evaluation of middle cerebral artery occlusion in prediction of recanalization and favorable clinical outcome and after intravenous thrombolysis (IV-tPA). METHODS: In 80 patients treated with IV-tPA for acute complete middle cerebral artery/M1 occlusion determined using CT angiography and temporal maximum intensity projection, calculated from 4D-CT angiography, the length of middle cerebral artery proximal stump, occlusion in M1 or M1 and M2 segment were measured. Univariate and multivariate analyses were performed to define independent predictors of successful recanalization after 24 hours and favorable outcome after 3 months. RESULTS: The length of occlusion was measureable in all patients using temporal maximum intensity projection. Recanalization thrombolysis in myocardial infarction 2 to 3 was achieved in 37 individuals (46%). The extension to M2 segment as a category (odds ratio, 4.58; 95% confidence interval, 1.39-15.05; P=0.012) and the length of M1 segment occlusion (odds ratio, 0.82; 95% confidence interval, 0.73-0.92; P=0.0007) with an optimal cutoff value of 12 mm (sensitivity 0.67; specificity 0.71) were significant independent predictors of recanalization. Favorable outcome (modified Rankin scale 0-2) was achieved in 25 patients (31%), baseline National Institutes of Health Stroke Scale (odds ratio, 0.82; 95% confidence interval, 0.72-0.93; P=0.003) and the length of occlusion M1 in segment (odds ratio, 0.79; 95% confidence interval, 0.69-0.91; P=0.0008) with an optimal cutoff value of 11 mm (sensitivity 0.74; specificity 0.76) were significant independent predictors of favorable outcome. CONCLUSIONS: The length of middle cerebral artery occlusion is an independent predictor of successful IV-tPA treatment.

Influence of trigger type, tube voltage and heart rate on calcified plaque imaging in dual source cardiac computed tomography: phantom study.

BACKGROUND: To investigate the impact of high pitch cardiac CT vs. retrospective ECG gated CT on the quantification of calcified vessel stenoses, with assessment of the influence of tube voltage, reconstruction kernel and heart rate. METHODS: A 4D cardiac movement phantom equipped with three different plaque phantoms (12.5%, 25% and 50% stenosis at different calcification levels), was scanned with a 128-row dual source CT scanner, applying different trigger types (gated vs. prospectively triggered high pitch), tube voltages (100-120 kV) and heart rates (50-90 beats per minute, bpm). Images were reconstructed using different standard (B26f, B46f, B70f) and iterative (I26f, I70f) convolution kernels. Absolute and relative plaque sizes were measured and statistically compared. Radiation dose associated with the different methods (gated vs. high pitch, 100 kV vs. 120 kV) were compared. RESULTS: Compared to the known diameters of the phantom plaques and vessels both CT-examination techniques overestimated the degrees of stenoses. Using the high pitch CT-protocol plaques appeared larger (0.09 ± 0.31 mm, 2 ± 8 percent points, PP) in comparison to the ECG-gated CT-scans. Reducing tube voltage had a similar effect, resulting in higher grading of the same stenoses by 3 ± 8 PP. In turn, sharper convolution kernels lead to a lower grading of stenoses (differences of up to 5%). Pairwise comparison of B26f and I26f, B46f and B70f, and B70f and I70f showed differences of 0-1 ± 6-8 PP of the plaque depiction. Motion artifacts were present only at 90 bpm high pitch experiments. High-pitch protocols were associated with significantly lower radiation doses compared with the ECG-gated protocols (258.0 mGy vs. 2829.8 mGy CTDIvol, p ≤ 0.0001). CONCLUSION: Prospectively triggered high-pitch cardiac CT led to an overestimation of plaque diameter and degree of stenoses in a coronary phantom. This overestimation is only slight and probably negligible in a clinical situation. Even at higher heart rates high pitch CT-scanning allowed reliable measurements of plaque and vessel diameters with only slight differences compared ECG-gated protocols, although motion artifacts were present at 90 bpm using the high pitch protocols.

Stereotactic body radiotherapy for pulmonary metastases. Prognostic factors and adverse respiratory events.

PURPOSE: The aim of this retrospective study was to evaluate the feasibility, safety, and effectiveness of stereotactic body radiotherapy (SBRT) for pulmonary metastases. PATIENTS AND METHODS: Between April 2007 and March 2011, 87 patients underwent SBRT for pulmonary metastases using the in-house Air-Bag System(TM) to obtain the four-dimensional image for treatment planning and to reduce intrafractional intrathoracic organ motion with abdominal compression to reduce the risk of radiation pneumonitis. Survival and respiratory adverse events were analyzed. RESULTS: The 2- and 3-year overall survival (OS) rates were 47 and 32 %, and the corresponding cause-specific survivals were 52 and 36 %. The 2- and 3-year OS rates were 57 and 49 % for patients in group 1, respectively, while the corresponding OS rates were 48 and 21 %, and 40 and 32 % for patients in groups 2 and 3, respectively. The 2- and 3-year local control (LC) rates were 80 and 80 %, respectively. The corresponding intrathoracic progression-free survival rates were 40 and 32 %, respectively. Concerning adverse respiratory events after SBRT for pulmonary metastases, 14 % were grade 0 (G0), 66 % G1, 13 % G2, 6 % G3, and 1 % G4. Concerning the adverse respiratory events (NCI-CTC) by grade scale, 1- and 2-year cumulative probabilities of radiation pneumonitis were 12 and 20 % for G2 and 4 and 10 % for G3/4, respectively. The mean values for cumulative V20 were 11.6 ± 8.5 %, 29.8 ± 18.6 %, and 25.7 ± 12.8 % in G0/1, G2, and G3/4, respectively. The number of pulmonary metastases that could be safely treated with SBRT was 6 PTVs (or seven gross tumor volumes) within a cumulative V20 of 30 % under the restricted intrafractional respiratory tumor motion using the Air-Bag System(TM). CONCLUSION: We propose that the number of pulmonary metastases that can be safely treated with SBRT is 6 PTVs with a cumulative V20 of 30 % under the restricted respiratory tumor motion using the Air-Bag System(TM). SBRT for pulmonary metastases offers locally effective treatment for recurrent or residual lesions after first line chemotherapy.

Tissue characterization using a phantom to validate four-dimensional tissue deformation.

PURPOSE: This project proposes using a real tissue phantom for 4D tissue deformation reconstruction (4DTDR) and 4D deformable image registration (DIR) validation, which allows for the complete verification of the motion path rather than limited end-point to end-point of motion. METHODS: Three electro-magnetic-tracking (EMT) fiducials were implanted into fresh porcine liver that was subsequently animated in a clinically realistic phantom. The animation was previously shown to be similar to organ motion, including hysteresis, when driven using a real patient's breathing pattern. For this experiment, 4DCTs and EMT traces were acquired when the phantom was animated using both sinusoidal and recorded patient-breathing traces. Fiducial were masked prior to 4DTDR for reconstruction. The original 4DCT data (with fiducials) were sampled into 20 CT phase sets and fiducials' coordinates were recorded, resulting in time-resolved fiducial motion paths. Measured values of fiducial location were compared to EMT measured traces and the result calculated by 4DTDR. RESULTS: For the sinusoidal breathing trace, 95% of EMT measured locations were within 1.2 mm of the measured 4DCT motion path, allowing for repeatable accurate motion characterization. The 4DTDR traces matched 95% of the EMT trace within 1.6 mm. Using the more irregular (in amplitude and frequency) patient trace, 95% of the EMT trace points fitted both 4DCT and 4DTDR motion path within 4.5 mm. The average match of the 4DTDR estimation of the tissue hysteresis over all CT phases was 0.9 mm using a sinusoidal signal for animation and 1.0 mm using the patient trace. CONCLUSIONS: The real tissue phantom is a tool which can be used to accurately characterize tissue deformation, helping to validate or evaluate a DIR or 4DTDR algorithm over a complete motion path. The phantom is capable of validating, evaluating, and quantifying tissue hysteresis, thereby allowing for full motion path validation.

The effect of breathing irregularities on quantitative accuracy of respiratory gated PET∕CT.

PURPOSE: 4D positron emission tomography and computed tomography (PET∕CT) can be used to reduce motion artifacts by correlating the raw PET data with the respiratory cycle. The accuracy of each PET phase is dependent on the reproducibility and consistency of the breathing cycle during acquisition. The objective of this study is to evaluate the impact of breathing amplitude and phase irregularities on the quantitative accuracy of 4D PET standardized uptake value (SUV) measurements. In addition, the magnitude of quantitative errors due to respiratory motion and partial volume error are compared. METHODS: Phantom studies were performed using spheres filled with (18)F ranging from 9 to 47 mm in diameter with background activity. Motion was simulated using patient breathing data. The authors compared the accuracy of SUVs derived from gated PET (4 bins and 8 bins, phase-based) for ideal, average, and highly irregular breathing patterns. RESULTS: Under ideal conditions, gated PET produced SUVs that were within (-5.4 ± 5.3)% of the static phantom measurements averaged across all sphere sizes. With breathing irregularities, the quantitative accuracy of gated PET decreased. Gated PET SUVs (best of 4 bins) were (-9.6 ± 13.0)% of the actual value for an average breather and decreased to (-17.1 ± 10.8)% for a highly irregular breather. Without gating, the differences in the SUV from actual value were (-28.5 ± 18.2)%, (-25.9 ± 14.4)%, and (-27.9 ± 18.2)% for the ideal, average, and highly irregular breather, respectively. CONCLUSIONS: Breathing irregularities reduce the quantitative accuracy of gated PET∕CT. Current gated PET techniques may underestimate the actual lesion SUV due to phase assignment errors. Evaluation of respiratory trace is necessary to assess accuracy of data binning and its effect on 4D PET SUVs.

4D patient dose reconstruction using online measured EPID cine images for lung SBRT treatment validation.

PURPOSE: This study aims to develop an EPID-guided 4D patient dose reconstruction framework and to investigate its feasibility for lung SBRT treatment validation. METHODS: Both the beam apertures and tumor movements were detected based on the continuously acquired EPID images during the treatment. Instead of directly using the transit photon fluence measured by the EPID, this method reconstructed the entrance fluence with the measured beam apertures and the delivered MUs. The entrance fluence distributions were sorted into their corresponding phases based on the detected tumor motion pattern and then accumulated for each phase. Together with the in-room 4DCT taken before every treatment to consider the interfractional-motion, the entrance fluence was then used for the patient dose calculation. Deformable registration was performed to sum up the phase doses for final treatment assessment. The feasibility of using the transit EPID images for entrance fluence reconstruction was evaluated against EPID in-air measurements. The accuracy of 3D- and 4D-dose reconstruction was validated by experiments with a motor-driven cylindrical diode array for six clinical-SBRT plans. RESULTS: The average difference between the measured and reconstructed fluence maps was within 0.16%. The reconstructed 3D-dose showed a less than 1.4% difference for the CAX-dose and at least a 98.3% gamma-passing-rate (2%∕2 mm) for the peripheral dose. Distorted dose distributions were observed in the measurement with the moving phantom. The comparison between the measured and the reconstructed 4D-dose without considering temporal information failed the gamma-evaluation for most cases. In contrast, when temporal information was considered, the dose distortion phenomena were successfully represented in the reconstructed dose (97.6%-99.7% gamma-passing rate). CONCLUSIONS: The proposed method considered uncertainties of the beam delivery system, the interfractional- and intrafractional-motion, and the interplay effect. The experimental validation demonstrates that this method is practical and accurate for online or offline SBRT patient dose verification.

A phantom for testing of 4D-CT for radiotherapy of small lesions.

PURPOSE: The use of time-resolved four-dimensional computed tomography (4D-CT) in radiotherapy requires strict quality assurance to ensure the accuracy of motion management protocols. The aim of this work was to design and test a phantom capable of large amplitude motion for use in 4D-CT, with particular interest in small lesions typical for stereotactic body radiotherapy. METHODS: The phantom of "see-saw" design is light weight, capable of including various sample materials and compatible with several surrogate marker signal acquisition systems. It is constructed of polymethylmethacrylate (Perspex) and its movement is controlled via a dc motor and drive wheel. It was tested using two CT scanners with different 4D acquisition methods: the Philips Brilliance Big Bore CT (helical scan, pressure belt) and a General Electric Discovery STE PET∕CT (axial scan, infrared marker). Amplitudes ranging from 1.5 to 6.0 cm and frequencies of up to 40 cycles per minute were used to study the effect of motion on image quality. Maximum intensity projections (MIPs), as well as average intensity projections (AIPs) of moving objects were investigated and their quality dependence on the number of phase reconstruction bins assessed. RESULTS: CT number discrepancies between moving and stationary objects were found to have no systematic dependence on amplitude, frequency, or specific interphase variability. MIP-delineated amplitudes of motion were found to match physical phantom amplitudes to within 2 mm for all motion scenarios tested. Objects undergoing large amplitude motions (>3.0 cm) were shown to cause artefacts in MIP and AIP projections when ten phase bins were assigned. This problem can be mitigated by increasing the number of phase bins in a 4D-CT scan. CONCLUSIONS: The phantom was found to be a suitable tool for evaluating the image quality of 4D-CT motion management technology, as well as providing a quality assurance tool for intercenter∕intervendor testing of commercial 4D-CT systems. When imaging objects with large amplitudes, the completeness criterion described here indicates the number of phase bins required to prevent missing data in MIPs and AIPs. This is most relevant for small lesions undergoing large motions.

Four dimensional CT imaging: a review of current technologies and modalities.

Organ motion is a substantial concern in the treatment of thoracic tumours using radiotherapy. A number of technologies have evolved in order to address this both during computed tomography (CT) imaging and radiation delivery. This review paper investigates the various technologies which have been developed in the field of CT scanning as well as their accuracy, cost and the implications of their clinical implementation. The scanning modalities covered include: slow CT, breath hold CT, gated CT and retrospectively correlated CT (4DCT). It was found that there are advantages and drawbacks to each of the mentioned techniques relating to patient dose, scan time, extra equipment and workload. Also some scanning techniques are only compatible with certain treatment modalities which would further influence the decision as to which technologies to implement.

Proton radiography and fluoroscopy of lung tumors: a Monte Carlo study using patient-specific 4DCT phantoms.

PURPOSE: Monte Carlo methods are used to simulate and optimize a time-resolved proton range telescope (TRRT) in localization of intrafractional and interfractional motions of lung tumor and in quantification of proton range variations. METHODS: The Monte Carlo N-Particle eXtended (MCNPX) code with a particle tracking feature was employed to evaluate the TRRT performance, especially in visualizing and quantifying proton range variations during respiration. Protons of 230 MeV were tracked one by one as they pass through position detectors, patient 4DCT phantom, and finally scintillator detectors that measured residual ranges. The energy response of the scintillator telescope was investigated. Mass density and elemental composition of tissues were defined for 4DCT data. RESULTS: Proton water equivalent length (WEL) was deduced by a reconstruction algorithm that incorporates linear proton track and lateral spatial discrimination to improve the image quality. 4DCT data for three patients were used to visualize and measure tumor motion and WEL variations. The tumor trajectories extracted from the WEL map were found to be within 1 mm agreement with direct 4DCT measurement. Quantitative WEL variation studies showed that the proton radiograph is a good representation of WEL changes from entrance to distal of the target. CONCLUSIONS: MCNPX simulation results showed that TRRT can accurately track the motion of the tumor and detect the WEL variations. Image quality was optimized by choosing proton energy, testing parameters of image reconstruction algorithm, and comparing to ground truth 4DCT. The future study will demonstrate the feasibility of using the time resolved proton radiography as an imaging tool for proton treatments of lung tumors.

An investigation of 4D cone-beam CT algorithms for slowly rotating scanners.

PURPOSE: To evaluate several algorithms for 4D cone-beam computed tomography (4D CBCT) with slow rotating devices. 4D CBCT is used to perform phase-correlated (PC) reconstructions of moving objects, such as breathing patients, for example. Such motion phase-dependent reconstructions are especially useful for updating treatment plans in radiation therapy. The treatment plan can be registered more precisely to the motion of the tumor and, in consequence, the irradiation margins for the treatment, the so-called planning target volume, can be reduced significantly METHODS: In the study, several algorithms were evaluated for kilovoltage cone-beam CT units attached to linear particle accelerators. The reconstruction algorithms were the conventional PC reconstruction, the McKinnon-Bates (MKB) algorithm, the prior image constrained compressed sensing (PICCS) approach, a total variation minimization (ASD-POCS) algorithm, and the auto-adaptive phase correlation (AAPC) algorithm. For each algorithm, the same motion-affected raw data were used, i.e., one simulated and one measured data set. The reconstruction results from the authors' implementation of these algorithms were evaluated regarding their noise and artifact levels, their residual motion blur, and their computational complexity and convergence. RESULTS: In general, it turned out that the residual motion blur was lowest in those algorithms which exclusively use data from a single motion phase. Algorithms using the image from the full data set as initialization or as a reference for the reconstruction were not capable of fully removing the motion blurring. The iterative algorithms, especially approaches based on total variation minimization, showed lower noise and artifact levels but were computationally complex. The conventional methods based on a single filtered backprojection were computationally inexpensive but suffered from higher noise and streak artifacts which limit the usability. In contrast, these methods showed to be less demanding and more predictable in their outcome than the total variation minimization based approaches. CONCLUSIONS: The reconstruction algorithms including at least one iterative step can reduce the 4 CBCT specific artifacts. Nevertheless, the algorithms that use the full data set, at least for initialization, such as MKB and PICCS in the authors' implementation, are only a trade-off and may not fully achieve the optimal temporal resolution. A predictable image quality as seen in conventional reconstruction methods, i.e., without total variation minimization, is a desirable property for reconstruction algorithms. Fast, iterative approaches such as the MKB can therefore be seen as a suitable tradeoff.

Evaluation and commissioning of a surface based system for respiratory sensing in 4D CT.

The purpose of this study is to assess the temporal and reconstruction accuracy of a surface imaging system, the GateCT under ideal conditions, and compare the device with a commonly used respiratory surrogate: the Varian RPM. A clinical CT scanner, run in cine mode, was used with two optical devices, GateCT and RPM, to detect respiratory motion. A radiation detector, GM-10, triggers the X-ray on/off to GateCT system, while the RPM is directly synchronized with the CT scanner through an electronic connection. Two phantoms were imaged: the first phantom translated on a rigid plate along the anterior-posterior (AP) direction, and was used to assess the temporal synchronization of each optical system with the CT scanner. The second phantom, consisting of five spheres translating 3 cm peak-to-peak in the superior-inferior direction, was used to assess the quality of rebinned images created by GateCT and RPM. Calibration assessment showed a nearly perfect synchronization with the scanner for both the RPM and GateCT systems, thus demonstrating the good performance of the radiation detector. Results for the volume rebinning test showed discrepancies in volumes for the 3D reconstruction (compared to ground truth) of up to 36% for GateCT and up to 40% for RPM. No statistical difference was proven between the two systems in volume sorting. Errors are mainly due to phase detection inaccuracies and to the large motion of the phantom. This feasibility study assessed the consistency of two optical systems in synchronizing the respiratory signal with the image acquisition. A new patient protocol based on both RPM and GateCT will be soon started.

Development of a deformable lung phantom for the evaluation of deformable registration.

The deformable lung phantom was developed to account for the patient breathing motion, and to evaluate for a deformable image registration algorithm. The phantom consisted of an acryl cylinder filled with water and a latex balloon located in the inner space of the cylinder. A silicon membrane was attached to the inferior end of the phantom. This silicon membrane was designed to simulate a real lung diaphragm and to reduce motor workload. This specific design was able to reduce the metal use which may prevent infrared sensing of the real position management (RPM) gating system on 4D CT image acquisition. Verification of intensity based 3D demon deformable registration was based on peak exhale and peak inhale breathing phases. The registration differences ranged from 0.85 mm to 1.47 mm, and accuracy was determined according to inner target deformation. This phantom was able to simulate the features and deformation of real human lung and has the potential for wide application in 4D radiation treatment planning.

Consistent and invertible deformation vector fields for a breathing anthropomorphic phantom: a post-processing framework for the XCAT phantom.

Breathing motion is challenging for radiotherapy planning and delivery. This requires advanced four-dimensional (4D) imaging and motion mitigation strategies and associated validation tools with known deformations. Numerical phantoms such as the XCAT provide reproducible and realistic data for simulation-based validation. However, the XCAT generates partially inconsistent and non-invertible deformations where tumours remain rigid and structures can move through each other. We address these limitations by post-processing the XCAT deformation vector fields (DVF) to generate a breathing phantom with realistic motion and quantifiable deformation. An open-source post-processing framework was developed that corrects and inverts the XCAT-DVFs while preserving sliding motion between organs. Those post-processed DVFs are used to warp the first XCAT-generated image to consecutive time points providing a 4D phantom with a tumour that moves consistently with the anatomy, the ability to scale lung density as well as consistent and invertible DVFs. For a regularly breathing case, the inverse consistency of the DVFs was verified and the tumour motion was compared to the original XCAT. The generated phantom and DVFs were used to validate a motion-including dose reconstruction (MIDR) method using isocenter shifts to emulate rigid motion. Differences between the reconstructed doses with and without lung density scaling were evaluated. The post-processing framework produced DVFs with a maximum [Formula: see text]-percentile inverse-consistency error of 0.02 mm. The generated phantom preserved the dominant sliding motion between the chest wall and inner organs. The tumour of the original XCAT phantom preserved its trajectory while deforming consistently with the underlying tissue. The MIDR was compared to the ground truth dose reconstruction illustrating its limitations. MIDR with and without lung density scaling resulted in small dose differences up to 1 Gy (prescription 54 Gy). The proposed open-source post-processing framework overcomes important limitations of the original XCAT phantom and makes it applicable to a wider range of validation applications within radiotherapy.