Respiratory-gated CT Research Articles

Selected advanced imaging techniques were unable to quantify in vivo parasitic burden in heartworm-infested dogs.

This prospective exploratory study aimed to determine whether certain noninvasive advanced imaging techniques could estimate parasitic burden in heartworm-infested dogs; a noninvasive method is needed for ethical considerations and permitting longitudinal drug studies. Three cardiac-gated and respiratory-gated 3T MRI techniques and CT pulmonary angiography (CTPA) were performed in three healthy beagles to optimize imaging techniques. Once the imaging techniques were established, a pilot study was performed to determine which one of the MRI techniques would be used in an observer comparison study. Ultimately, spoiled gradient recalled (SPGR)-cine-MRI and CTPA were performed in four and five heartworm-infested dogs, respectively. Heartworms were detected in the pulmonary arteries in all dogs during SPGR-cine-MRI and in no dog during CTPA. However, counting the number of worms was unsuccessful. In conclusion, CTPA and SPGR-cine-MRI were unable to replace necropsy for quantifying parasitic burden in heartworm-infested dogs.

Comparison of Different Attenuation Correction Methods for Dual Gating Myocardial Perfusion SPECT/CT

Dual respiratory-cardiac gating realizes motion-freeze myocardial perfusion single photon emission computed tomography (SPECT). The aim of this paper is to evaluate the performance of various attenuation maps for dual gating cardiac SPECT. We used extended cardiac-torso phantom with Tc-99m-sestamibi distribution to generate a dataset with both respiratory and cardiac motion and then grouped them based on six respiratory and eight cardiac phases. The corresponding attenuation correction (AC) maps were used to simulate different CT protocols: dual gating CT (DCT), respiratory gated CT (RCT), cine average CT (ACT), interpolated CT (ICT), helical CTs at end-inspiration (HCT-in), end-expiration (HCT-ex), and mid-respiration (HCT-mid), respectively. One hundred and twenty noise-free and realistic noisy projections were generated over 180° from RAO to LPO and reconstructed by OS-EM method using different AC maps. The reconstructed images in different respiratory phases were registered to one reference end-expiration phase and summed up as a cardiac image for each cardiac phase. The corresponding polar maps were generated and further analyzed using 17-segment plots. Relative difference (RD) of the average intensity was computed in each segment using the polar maps with DCT-AC as reference. Substantial artifacts could be found in HCT-in AC images in every cardiac phase. For all eight noise-free cardiac phases, the average RDmax for AC with ACT, HCT-ex, HCT-in, HCT-mid, RCT, and ICT comparing to DCT were 4.76%, 4.95%, 10.54%, 4.23%, 2.39%, and 2.76%, respectively. For all eight noisy cardiac phases, the average RDmax for AC with ACT, HCT-ex, HCT-in, HCT-mid, RCT, and ICT comparing to DCT were 4.09%, 5.04%, 9.24%, 4.02%, 3.22%, and 2.85%, respectively. The ICT and RCT have the best AC performance for dual gating SPECT, while HCT-ex, HCT-in, HCT-mid, and ACT show inferior results. Since RCT and DCT are not commonly used in clinics due to complex implementation and radiation concern, ICT is recommended, followed by HCT-mid or HCT-ex for AC in dual gating cardiac SPECT with good accuracy and relative low radiation dose.

Real-time respiratory phase matching between 2D fluoroscopic images and 3D CT images for precise percutaneous lung biopsy.

A 3D CT image is used along with real-time 2D fluoroscopic images in the state-of-the-art cone-beam CT system to guide percutaneous lung biopsy (PLB). To improve the guiding accuracy by compensating for respiratory motion, we propose an algorithm for real-time matching of 2D fluoroscopic images to multiple 3D CT images of different respiratory phases that is robust to the small movement and deformation due to cardiac motion. Based on the transformations obtained from nonrigid registration between two 3D CT images acquired at expiratory and inspiratory phases, we first generate sequential 3D CT images (or a 4D CT image) and the corresponding 2D digitally reconstructed radiographs (DRRs) of vessels. We then determine 3D CT images corresponding to each real-time 2D fluoroscopic image, by matching the 2D fluoroscopic image to a 2D DRR. Quantitative evaluations performed with 20 clinical datasets show that registration errors of anatomical features between a 2D fluoroscopic image and its matched 2D DRR are less than 3mm on average. Registration errors of a target lesion are determined to be roughly 3mm on average for 10 datasets. We propose a real-time matching algorithm to compensate for respiratory motion between a 2D fluoroscopic image and 3D CT images of the lung, regardless of cardiac motion, based on a newly improved matching measure. The proposed algorithm can improve the accuracy of a guiding system for the PLB by providing 3D images precisely registered to 2D fluoroscopic images in real-time, without time-consuming respiratory-gated or cardiac-gated CT images.

Multi-modal anatomical Optical Coherence Tomography and CT for in vivo Dynamic Upper Airway Imaging.

We describe a novel, multi-modal imaging protocol for validating quantitative dynamic airway imaging performed using anatomical Optical Coherence Tomography (aOCT). The aOCT system consists of a catheter-based aOCT probe that is deployed via a bronchoscope, while a programmable ventilator is used to control airway pressure. This setup is employed on the bed of a Siemens Biograph CT system capable of performing respiratory-gated acquisitions. In this arrangement the position of the aOCT catheter may be visualized with CT to aid in co-registration. Utilizing this setup we investigate multiple respiratory pressure parameters with aOCT, and respiratory-gated CT, on both ex vivo porcine trachea and live, anesthetized pigs. This acquisition protocol has enabled real-time measurement of airway deformation with simultaneous measurement of pressure under physiologically relevant static and dynamic conditions- inspiratory peak or peak positive airway pressures of 10-40 cm H2O, and 20-30 breaths per minute for dynamic studies. We subsequently compare the airway cross sectional areas (CSA) obtained from aOCT and CT, including the change in CSA at different stages of the breathing cycle for dynamic studies, and the CSA at different peak positive airway pressures for static studies. This approach has allowed us to improve our acquisition methodology and to validate aOCT measurements of the dynamic airway for the first time. We believe that this protocol will prove invaluable for aOCT system development and greatly facilitate translation of OCT systems for airway imaging into the clinical setting.

Dynamic CT imaging of volumetric changes in pulmonary nodules correlates with physical measurements of stiffness

Background and purposeA major challenge in CT screening for lung cancer is limited specificity when distinguishing between malignant and non-malignant pulmonary nodules (PN). Malignant nodules have different mechanical properties and tissue characteristics (‘stiffness’) from non-malignant nodules. This study seeks to improve CT specificity by demonstrating in rats that measurements of volumetric ratios in PNs with varying composition can be determined by respiratory-gated dynamic CT imaging and that these ratios correlate with direct physical measurements of PN stiffness. Methods and materialsRespiratory-gated MicroCT images acquired at extreme tidal volumes of 9 rats with PNs from talc, matrigel and A549 human lung carcinoma were analyzed and their volumetric ratios (δ) derived. PN stiffness was determined by measuring the Young’s modulus using atomic force microscopy (AFM) for each nodule excised immediately after MicroCT imaging. ResultsThere was significant correlation (p=0.0002) between PN volumetric ratios determined by respiratory-gated CT imaging and the physical stiffness of the PNs determined from AFM measurements. ConclusionWe demonstrated proof of concept that PN volume changes measured non-invasively correlate with direct physical measurements of stiffness. These results may translate clinically into a means of improving the specificity of CT screening for lung cancer and/or improving individual prognostic assessments based on lung tumor stiffness.

Impact of treatment planning with deformable image registration on dose distribution for carbon-ion beam lung treatment using a fixed irradiation port and rotating couch.

When using a fixed irradiation port, treatment couch rotation is necessary to increase beam angle selection. We evaluated dose variations associated with positional morphological changes to organs. We retrospectively chose the data sets of ten patients with lung cancer who underwent respiratory-gated CT at three different couch rotation angles (0°, 20° and -20°). The respective CT data sets are referred to as CT0, CT20 and CT-20. Three treatment plans were generated as follows: in Plan 1, all compensating bolus designs and dose distributions were calculated using CT0. To evaluate the rotation effect without considering morphology changes, in Plan 2, the compensating boli designed using CT0 were applied to the CT±20 images. Plan 3 involved compensating boli designed using the CT±20 images. The accumulated dose distributions were calculated using deformable image registration (DIR). A sufficient prescribed dose was calculated for the planning target volume (PTV) in Plan 1 [minimum dose received by a volume ≥95% (D95) > 95.8%]. By contrast, Plan 2 showed degraded dose conformation to the PTV (D95 > 90%) owing to mismatch of the bolus design to the morphological positional changes in the respective CT. The dose assessment results of Plan 3 were very close to those of Plan 1. Dose distribution is significantly affected by whether or not positional organ morphology changes are factored into dose planning. In treatment planning using multiple CT scans with different couch positions, it is mandatory to calculate the accumulated dose using DIR.

Monitoring internal organ motion with continuous wave radar in CT

To avoid motion artifacts in medical imaging or to minimize the exposure of healthy tissues in radiation therapy, medical devices are often synchronized with the patient's respiratory motion. Today's respiratory motion monitors require additional effort to prepare the patients, e.g., mounting a motion belt or placing an optical reflector on the patient's breast. Furthermore, they are not able to measure internal organ motion without implanting markers. An interesting alternative to assess the patient's organ motion is continuous wave radar. The aim of this work is to design, implement, and evaluate such a radar system focusing on application in CT. The authors designed a radar system operating in the 860 MHz band to monitor the patient motion. In the intended application of the radar system, the antennas are located close to the patient's body inside the table of a CT system. One receive and four transmitting antennas are used to avoid the requirement of exact patient positioning. The radar waves propagate into the patient's body and are reflected at tissue boundaries, for example at the borderline between muscle and adipose tissue, or at the boundaries of organs. At present, the authors focus on the detection of respiratory motion. The radar system consists of the hardware mentioned above as well as of dedicated signal processing software to extract the desired information from the radar signal. The system was evaluated using simulations and measurements. To simulate the radar system, a simulation model based on radar and wave field equations was designed and 4D respiratory-gated CT data sets were used as input. The simulated radar signals and the measured data were processed in the same way. The radar system hardware and the signal processing algorithms were tested with data from ten volunteers. As a reference, the respiratory motion signal was recorded using a breast belt simultaneously with the radar measurements. Concerning the measurements of the test persons, there is a very good correlation (ρ = 0.917) between the respiratory motion phases received by the radar system and the external motion monitor. Our concept of using an array of transmitting antennas turned out to be widely insensitive to the positioning of the test persons. A time shift between the respiratory motion curves recorded with the radar system and the motion curves from the external respiratory monitor was observed which indicates a slight difference between internal organ motion and motion detected by the external respiratory monitor. The simulations were in good accordance with the measurements. A continuous wave radar operating in the near field of the antennas can be used to determine the respiratory motion of humans accurately. In contrast to trigger systems used today, the radar system is able to measure motion inside the body. If such a monitor was routinely available in clinical CT, it would be possible optimizing the scan start with respect to the respiratory state of the patient. Breathing commands would potentially widely be avoided, and as far as uncooperative patients or children are concerned, less sedation might be necessary. Further applications of the radar system could be in radiation therapy or interventional imaging for instance.

Reproducibility of intensity‐based estimates of lung ventilation

Lung function depends on lung expansion and contraction during the respiratory cycle. Respiratory-gated CT imaging and image registration can be used to estimate the regional lung volume change by observing CT voxel density changes during inspiration or expiration. In this study, the authors examine the reproducibility of intensity-based estimates of lung tissue expansion and contraction in three mechanically ventilated sheep and ten spontaneously breathing humans. The intensity-based estimates are compared to the estimates of lung function derived from image registration deformation field. 4DCT data set was acquired for a cohort of spontaneously breathing humans and anesthetized and mechanically ventilated sheep. For each subject, two 4DCT scans were performed with a short time interval between acquisitions. 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. The CT density change in the registered image pair was used to compute intensity-based specific air volume change (SAC) and the intensity-based Jacobian (IJAC), while the transformation-based Jacobian (TJAC) was computed directly from the image registration deformation field. IJAC is introduced to make the intensity-based and transformation-based methods comparable since SAC and Jacobian may not be associated with the same physiological phenomenon and have different units. Scan-to-scan variations in respiratory effort were corrected using a global scaling factor for normalization. A gamma index metric was introduced to quantify voxel-by-voxel reproducibility considering both differences in ventilation and distance between matching voxels. The authors also tested how different CT prefiltering levels affected intensity-based ventilation reproducibility. Higher reproducibility was found for anesthetized mechanically ventilated animals than for the humans for both the intensity-based (IJAC) and transformation-based (TJAC) ventilation estimates. The human IJAC maps had scan-to-scan correlation coefficients of 0.45 ± 0.14, a gamma pass rate 70 ± 8 without normalization and 75 ± 5 with normalization. The human TJAC maps had correlation coefficients 0.81 ± 0.10, a gamma pass rate 86 ± 11 without normalization and 93 ± 4 with normalization. The gamma pass rate and correlation coefficient of the IJAC maps gradually increased with increased smoothing, but were still much lower than those of the TJAC maps. The transformation-based ventilation maps show better reproducibility than the intensity-based maps, especially in human subjects. Reproducibility was also found to depend on variations in respiratory effort; all techniques were better when applied to images from mechanically ventilated sheep compared to spontaneously breathing human subjects. Nevertheless, intensity-based techniques applied to mechanically ventilated sheep were less reproducible than the transformation-based applied to spontaneously breathing humans, suggesting the method used to determine ventilation maps is important. Prefiltering of the CT images may help to improve the reproducibility of the intensity-based ventilation estimates, but even with filtering the reproducibility of the intensity-based ventilation estimates is not as good as that of transformation-based ventilation estimates.

Reproducibility of registration-based measures of lung tissue expansion

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. 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. 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. 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.

Current concepts in F18 FDG PET/CT-based radiation therapy planning for lung cancer

Radiation therapy is an important component of cancer therapy for early stage as well as locally advanced lung cancer. The use of F18 FDG PET/CT has come to the forefront of lung cancer staging and overall treatment decision-making. FDG PET/CT parameters such as standard uptake value and metabolic tumor volume provide important prognostic and predictive information in lung cancer. Importantly, FDG PET/CT for radiation planning has added biological information in defining the gross tumor volume as well as involved nodal disease. For example, accurate target delineation between tumor and atelectasis is facilitated by utilizing PET and CT imaging. Furthermore, there has been meaningful progress in incorporating metabolic information from FDG PET/CT imaging in radiation treatment planning strategies such as radiation dose escalation based on standard uptake value thresholds as well as using respiratory-gated PET and CT planning for improved target delineation of moving targets. In addition, PET/CT-based follow-up after radiation therapy has provided the possibility of early detection of local as well as distant recurrences after treatment. More research is needed to incorporate other biomarkers such as proliferative and hypoxia biomarkers in PET as well as integrating metabolic information in adaptive, patient-centered, tailored radiation therapy.

4D XCAT phantom for multimodality imaging research

The authors develop the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. Highly detailed whole-body anatomies for the adult male and female were defined in the XCAT using nonuniform rational B-spline (NURBS) and subdivision surfaces based on segmentation of the Visible Male and Female anatomical datasets from the National Library of Medicine as well as patient datasets. Using the flexibility of these surfaces, the Visible Human anatomies were transformed to match body measurements and organ volumes for a 50th percentile (height and weight) male and female. The desired body measurements for the models were obtained using the PEOPLESIZE program that contains anthropometric dimensions categorized from 1st to the 99th percentile for US adults. The desired organ volumes were determined from ICRP Publication 89 [ICRP, "Basic anatomical and physiological data for use in radiological protection: reference values," ICRP Publication 89 (International Commission on Radiological Protection, New York, NY, 2002)]. The male and female anatomies serve as standard templates upon which anatomical variations may be modeled in the XCAT through user-defined parameters. Parametrized models for the cardiac and respiratory motions were also incorporated into the XCAT based on high-resolution cardiac- and respiratory-gated multislice CT data. To demonstrate the usefulness of the phantom, the authors show example simulation studies in PET, SPECT, and CT using publicly available simulation packages. As demonstrated in the pilot studies, the 4D XCAT (which includes thousands of anatomical structures) can produce realistic imaging data when combined with accurate models of the imaging process. With the flexibility of the NURBS surface primitives, any number of different anatomies, cardiac or respiratory motions or patterns, and spatial resolutions can be simulated to perform imaging research. With the ability to produce realistic, predictive 3D and 4D imaging data from populations of normal and abnormal patients under various imaging parameters, the authors conclude that the XCAT provides an important tool in imaging research to evaluate and improve imaging devices and techniques. In the field of x-ray CT, the phantom may also provide the necessary foundation with which to optimize clinical CT applications in terms of image quality versus radiation dose, an area of research that is becoming more significant with the growing use of CT.

Tensor grid based image registration with application to ventilation estimation on 4D CT lung data

For many image registration tasks, the information contained in the original resolution of the image data is crucial for a subsequent medical analysis, e.g. accurate assessment of local pulmonary ventilation. However, the complexity of a non-parametric registration scheme is directly connected to the resolution of the images. Therefore, the registration is often performed on a downsampled version in order to meet runtime demands and thereby producing suboptimal results. To enable the application of the highest resolution at least in regions of high clinical importance, an approach is presented replacing the usually taken equidistant grids by tensor grids for image representation. We employ a non-parametric approach for the registration of a respiratory-gated 4D CT thorax scan. Tensor grids are introduced for the registration setting and compared to equidistant grids. For ventilation assessment, the Jacobian metric is explored. The application of the tensor grid approach makes the local usage of the original resolution feasible; thereby a smaller registration error is achieved in regions of higher resolution using the tensor grids, while the two types of grids perform similar in regions of equal resolution. Concerning the ventilation assessment, the Jacobian metric yields reasonable results, showing more detail using the tensor grids due to the higher resolution. The proposed approach using tensor grids preserves registration accuracy, while reducing computational demands. The application of the Jacobian metric for ventilation assessment in conjunction with tensor grids is promising; however, due to a missing ground-truth the medical relevance could not be established for the ventilation estimation so far.

Estimation of Cardiac Respiratory-Motion by Semi-Automatic Segmentation and Registration of Non-Contrast-Enhanced 4D-CT Cardiac Datasets

The goal of this work is to investigate, for a large set of patients, the motion of the heart with respiration during free-breathing supine medical imaging. For this purpose we analyzed the motion of the heart in 32 non-contrast enhanced respiratory-gated 4D-CT datasets acquired during quiet unconstrained breathing. The respiratory-gated CT images covered the cardiac region and were acquired at each of 10 stages of the respiratory cycle, with the first stage being end-inspiration. We devised a 3-D semi-automated segmentation algorithm that segments the heart in the 4D-CT datasets acquired without contrast enhancement for use in estimating respiratory motion of the heart. Our semi-automated segmentation results were compared against interactive hand segmentations of the coronal slices by a cardiologist and a radiologist. The pairwise difference in segmentation among the algorithm and the physicians was on the average 11% and 10% of the total average segmented volume across the patient, with a couple of patients as outliers above the 95% agreement limit. The mean difference among the two physicians was 8% with an outlier above the 95% agreement limit. The 3-D segmentation was an order of magnitude faster than the Physicians' manual segmentation and represents significant reduction of Physicians' time. The segmented first stages of respiration were used in 12 degree-of-freedom (DOF) affine registration to estimate the motion at each subsequent stage of respiration. The registration results from the 32 patients indicate that the translation in the superior-inferior direction was the largest component motion, with a maximum of 10.7 mm, mean of 6.4 mm, and standard deviation of 2.2 mm. Translation in the anterior-posterior direction was the next largest component of motion, with a maximum of 4.0 mm, mean of 1.7 mm, and standard deviation of 1.0 mm. Rotation about the right-left axis was on average the largest component of rotation observed, with a maximum of 4.6 degrees, mean of 1.6 degrees, and standard deviation of 2.1 degrees. The other rotation and shear parameters were all close to zero on average indicting the motion could be reasonably well approximated by rigid-body motion. However, the product of the three scale factors averaged about 0.97 indicating the possibility of a small decrease in heart volume with expiration. The motion results were similar whether we used the Physician's segmentation or the 3-D algorithm.

Registration-based estimates of local lung tissue expansion compared to xenon CT measures of specific ventilation

The main function of the respiratory system is gas exchange. Since many disease or injury conditions can cause biomechanical or material property changes that can alter lung function, there is a great interest in measuring regional lung ventilation and regional specific volume change. We describe a registration-based technique for estimating local lung expansion from multiple respiratory-gated CT images of the thorax. The degree of regional lung expansion is measured using the Jacobian (a function of local partial derivatives) of the registration displacement field, which we show is directly related to specific volume change. We compare the ventral-dorsal patterns of lung expansion estimated across five pressure changes to a xenon CT based measure of specific ventilation in five anesthetized sheep studied in the supine orientation. Using 3D image registration to match images acquired at 10 cm H(2)O and 15 cm H(2)O airway pressures gave the best match between the average Jacobian and the xenon CT specific ventilation (linear regression, average r(2)=0.73).

Respiratory-gated CT as a tool for the simulation of breathing artifacts in PET and PET/CT

Respiratory motion in PET and PET/CT blurs the images and can cause attenuation-related errors in quantitative parameters such as standard uptake values. In rare instances, this problem even causes localization errors and the disappearance of tumors that should be detectable. Attenuation errors are severe near the diaphragm and can be enhanced when the attenuation correction is based on a CT series acquired during a breath-hold. To quantify the errors and identify the parameters associated with them, the authors performed a simulated PET scan based on respiratory-gated CT studies of five lung cancer patients. Diaphragmatic motion ranged from 8 to 25 mm in the five patients. The CT series were converted to 511-keV attenuation maps which were forward-projected and exponentiated to form sinograms of PET attenuation factors at each phase of respiration. The CT images were also segmented to form a PET object, moving with the same motion as the CT series. In the moving PET object, spherical 20 mm mobile tumors were created in the vicinity of the dome of the liver and immobile 20 mm tumors in the midchest region. The moving PET objects were forward-projected and attenuated, then reconstructed in several ways: phase-matched PET and CT, gated PET with ungated CT, ungated PET with gated CT, and conventional PET. Spatial resolution and statistical noise were not modeled. In each case, tumor uptake recovery factor was defined by comparing the maximum reconstructed pixel value with the known correct value. Mobile 10 and 30 mm tumors were also simulated in the case of a patient with 11 mm of breathing motion. Phase-matched gated PET and CT gave essentially perfect PET reconstructions in the simulation. Gated PET with ungated CT gave tumors of the correct shape, but recovery was too large by an amount that depended on the extent of the motion, as much as 90% for mobile tumors and 60% for immobile tumors. Gated CT with ungated PET resulted in blurred tumors and caused recovery errors between -50% and +75%. Recovery in clinical scans would be 0%-20% lower than stated because spatial resolution was not included in the simulation. Mobile tumors near the dome of the liver were subject to the largest errors in either case. Conventional PET for 20 mm tumors was quantitative in cases of motion less than 15 mm because of canceling errors in blurring and attenuation, but the recovery factors were too low by as much as 30% in cases of motion greater than 15 mm. The 10 mm tumors were blurred by motion to a greater extent, causing a greater SUV underestimation than in the case of 20 mm tumors, and the 30 mm tumors were blurred less. Quantitative PET imaging near the diaphragm requires proper matching of attenuation information to the emission information. The problem of missed tumors near the diaphragm can be reduced by acquiring attenuation-correction information near end expiration. A simple PET/CT protocol requiring no gating equipment also addresses this problem.

Design of a compensating bolus by use of exhalation CT data for covering residual motion in respiratory-gated charged-particle lung therapy: four-dimensional carbon beam dose calculation

We developed an algorithm which we used to design a compensating bolus by using respiratory-gated CT data for respiratory-gated carbon beam lung therapy and evaluated it by calculating dose distributions as a function of time. Four-dimensional CT (4DCT) images were obtained for seven lung cancer patients under free breathing conditions. The internal target volume (ITV) was calculated by maximum intensity projection processing which use of three types of gross tumor volumes (GTVs): at peak exhalation and with a 5 mm shift of the GTV to both superior and inferior sides. Then a compensating bolus was designed which use of the ITV and applied to 4DCT data at the gating window (around exhalation phase). The carbon beam distribution was calculated by a pencil-beam algorithm as a function of time. The compensating bolus provides a sufficient prescribed dose to the target in the gating window and minimizes any excessive dose to the normal tissues. The metric of dosimetric assessment metrics of D95 in all patients is greater than 96% of the prescribed dose in the gating window. Our results will be beneficial for improving the accuracy of charged-particle radiotherapy for hospitals where 4DCT cannot be used.

Determination of Respiratory Motion for Distal Esophagus Cancer Using Four-Dimensional Computed Tomography

To investigate the motion characteristics of distal esophagus cancer primary tumors using four-dimensional computed tomography (4D CT). Thirty-one consecutive patients treated for esophagus cancer who received respiratory-gated 4D CT imaging for treatment planning were selected. Deformable image registration was used to map the full expiratory motion gross tumor volume (GTV) to the full-inspiratory CT image, allowing quantitative assessment of each voxel's displacement. These displacements were correlated with patient tumor and respiratory characteristics. The mean (SE) tidal volume was 608 (73) mL. The mean GTV volume was 64.3 (10.7) mL on expiration and 64.1 (10.7) mL on inspiration (no significant difference). The mean tumor motion in the x-direction was 0.13 (0.006) cm (average of absolute values), in the y-direction 0.23 (0.01) cm (anteriorly), and in the z-direction 0.71 (0.02) cm (inferiorly). Tumor motion correlated with tidal volume. Comparison of tumor motion above vs. below the diaphragm was significant for the average net displacement (p = 0.014), motion below the diaphragm was greater than above. From the cumulative distribution 95% of the tumors moved less than 0.80 cm radially and 1.75 cm inferiorly. Primary esophagus tumor motion was evaluated with 4D CT. According to the results of this study, when 4D CT is not available, a radial margin of 0.8 cm and axial margin of +/-1.8 cm would provide tumor motion coverage for 95% of the cases in our study population.

Improvement in banding artefacts in four-dimensional computed tomography for radiotherapy planning

Respiratory-gated CT (RGCT) and four-dimensional CT (4DCT) scan techniques cover consecutive segments of the respiratory cycle. However, motion artefacts may occur in fast respiratory phases such as mid-inhalation and -exhalation. CT imaging involves the use of a number of x-ray tube positions for each couch position. We investigated the fundamental nature of motion artefacts using a constant-velocity moving phantom in motion in the CT plane or perpendicular to the CT plane, and in pigs to simulate a human model. Artefacts and movement distance were evaluated in a moving phantom and artificially ventilated pigs with a 256-multi-detector row CT (256MDCT). The phantom moved in the CT plane or perpendicular to the CT plane with a constant velocity. Backprojection used variable initial backprojection angles (IBAs). The phantom length for motion perpendicular to the CT plane was independent of IBA but was represented by phantom diameter plus the distance of movement per gantry rotation. In contrast, that for the motion in the CT plane was dependent on IBA, as represented by phantom diameter plus the distance of movement per rotation for IBA perpendicular to the phantom movement direction, and phantom diameter plus half the distance of movement per gantry rotation for other IBAs. Results for volumetric CT images with different IBAs showed the presence of banding artefacts. Similar findings were seen in artificially ventilated pigs. Motion artefacts are unavoidable in both conventional CT and 256MDCT. Banding artefacts will be improved if the same IBAs at each couch position are accounted for during image reconstruction. This improvement will be beneficial in respiratory gated and 4D radiation therapies.

Improved motion compensation in 3D-CT using respiratory-correlated segment reconstruction: diagnostic and radiotherapy applications

Conventional respiratory-gated CT and four-dimensional CT (4DCT) are disadvantaged by their low temporal resolution, which results in the inclusion of anatomic motion-induced artefacts. These represent a significant source of error both in radiotherapy treatment planning for the thorax and upper abdomen and in diagnostic procedures. In particular, temporal resolution and image quality are vitally important to accurate diagnosis and the minimization of planning target volume margin due to respiratory motion. To improve both temporal resolution and signal-to-noise ratio (SNR), we developed a respiratory-correlated segment reconstruction method (RS) and adapted it to the Feldkamp-Davis-Kress algorithm (FDK) with a 256 multidetector row CT (256MDCT). The 256MDCT scans approximately 100 mm in the craniocaudal direction with a 0.5 mm slice thickness in one rotation. Data acquisition for the RS-FDK relies on the assistance of a respiratory sensing system operating in cine scan mode (continuous axial scan with the table stationary). We evaluated the RS-FDK for volume accuracy and image noise in a phantom study with the 256MDCT and compared results with those for a full scan (FS-FDK), which is usually employed in conventional 4DCT and in half scan (HS-FDK). Results showed that the RS-FDK gave a more accurate volume than the others and had the same SNR as the FS-FDK. In a subsequent animal study, we demonstrated a practical sorting process for projection data which was unaffected by variations in respiratory period, and found that the RS-FDK gave the clearest visualization among the three algorithms of the margins of the liver and pulmonary vessels. In summary, the RS-FDK algorithm provides multi-phase images with higher temporal resolution and better SNR. This method should prove useful when combined with new radiotherapeutic and diagnostic techniques.

Volumetric cine imaging for four-dimensional radiation therapy planning using the second model of the 256-detector row CT-scanner: Initial experience in lung cancer

This is an initial demonstration of volumetric cine imaging under free breathing in lung cancer patients using the second model of the 256-detector row CT (256MDCT). Scan range is 128 mm in the cranio-caudal direction with a 0.5 mm slice thickness and effective temporal resolution of 250 ms. Volumetric cine images are satisfactorily obtained for the thin sections which can be used to create cine loops in multiple planes. Thus, the 256MDCT overcomes some of the limitations of current respiratory-gated CT or 4D-CT scan techniques and has the potential for use in 4D-IMRT planning and 4D radiation therapy.