Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

High resolution dual detector volume-of-interest cone beam breast CT - Demonstration with a bench top system

Purpose: In this study, we used a small field high resolution detector in conjunction with a full field flat panel detector to implement and investigate the dual detector volume-of-interest (VOI) cone beam breast computed tomography (CBCT) technique on a bench-top system. The potential of using this technique to image small calcifications without increasing the overall dose to the breast was demonstrated. Significant reduction of scatter components in the high resolution projection image data of the VOI was also shown. Methods: With the regular flat panel based CBCT technique, exposures were made at 80 kVp to generate an air kerma of 6 mGys at the isocenter. With the dual detector VOI CBCT technique, a high resolution small field CMOS detector was used to scan a cylindrical VOI (2.5 cm in diameter and height, 4.5 cm off-center) with collimated x-rays at four times of regular exposure level. A flat panel detector was used for full field scan with low x-ray exposures at half of the regular exposure level. The low exposure full field image data were used to fill in the truncated space in the VOI scan data and generate a complete projection image set. The Feldkamp-Davis-Kress (FDK) filtered backprojection algorithm was used to reconstruct high resolution images for the VOI. Two scanning techniques, one breast centered and the other VOI centered, were implemented and investigated. Paraffin cylinders with embedded thin aluminum (Al) wires were imaged and used in conjunction with optically stimulated luminescence (OSL) dose measurements to demonstrate the ability of this technique to image small calcifications without increasing the mean glandular dose (MGD). Results: Using exposures that produce an air kerma of 6 mGys at the isocenter, the regular CBCT technique was able to resolve the cross-sections of Al wires as thin as 254 {mu}m in diameter in the phantom. For the specific VOI studied, by increasing the exposure level by a factor of 4 for the VOI scan and reducing the exposure level by a factor of 2 for the full filed scan, the dual-detector CBCT technique was able to resolve the cross-sections of Al wires as thin as 152 {mu}m in diameter. The CNR evaluated for the entire Al wire cross-section was found to be improved from 5.5 in regular CBCT to 14.4 and 16.8 with the breast centered and VOI centered scanning techniques, respectively. Even inside VOI center, the VOI scan resulted in significant dose saving with the dose reduced by a factor of 1.6 at the VOI center. Dose saving outside the VOI was substantial with the dose reduced by a factor of 7.3 and 7.8 at the breast center for the breast centered and VOI centered scans, respectively, when compared to full field scan at the same exposure level. The differences between the two dual detector techniques in terms of dose saving and scatter reduction were small with VOI scan at 4x exposure level and full field scan at 0.5x exposure level. The MGDs were only 94% of that from the regular CBCT scan. Conclusions: For the specific VOI studied, the dual detector VOI CBCT technique has the potential to provide high quality images inside the VOI with MGD similar to or even lower than that of full field breast CBCT. It was also found that our results were compromised by the use of inadequate detectors for the VOI scan. An appropriately selected detector would better optimize the image quality improvement that can be achieved with the VOI CBCT technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Antiscatter grids in mobile C-arm cone-beam CT: Effect on image quality and dose

Purpose: X-ray scatter is a major detriment to image quality in cone-beam CT (CBCT). Existing geometries exhibit strong differences in scatter susceptibility with more compact geometries, e.g., dental or musculoskeletal, benefiting from antiscatter grids, whereas in more extended geometries, e.g., IGRT, grid use carries tradeoffs in image quality per unit dose. This work assesses the tradeoffs in dose and image quality for grids applied in the context of low-dose CBCT on a mobile C-arm for image-guided surgery. Methods: Studies were performed on a mobile C-arm equipped with a flat-panel detector for high-quality CBCT. Antiscatter grids of grid ratio (GR) 6:1-12:1, 40 lp/cm, were tested in ''body'' surgery, i.e., spine, using protocols for bone and soft-tissue visibility in the thoracic and abdominal spine. Studies focused on grid orientation, CT number accuracy, image noise, and contrast-to-noise ratio (CNR) in quantitative phantoms at constant dose. Results: There was no effect of grid orientation on possible gridline artifacts, given accurate angle-dependent gain calibration. Incorrect calibration was found to result in gridline shadows in the projection data that imparted high-frequency artifacts in 3D reconstructions. Increasing GR reduced errors in CT number from 31%, thorax, and 37%, abdomen, for gridless operation to 2% and 10%, respectively, with a 12:1 grid, while image noise increased by up to 70%. The CNR of high-contrast objects was largely unaffected by grids, but low-contrast soft-tissues suffered reduction in CNR, 2%-65%, across the investigated GR at constant dose. Conclusions: While grids improved CT number accuracy, soft-tissue CNR was reduced due to attenuation of primary radiation. CNR could be restored by increasing dose by factors of {approx}1.6-2.5 depending on GR, e.g., increase from 4.6 mGy for the thorax and 12.5 mGy for the abdomen without antiscatter grids to approximately 12 mGy and 30 mGy, respectively, with a high-GR grid. However, increasing the dose poses a significant impediment to repeat intraoperative CBCT and can cause the cumulative intraoperative dose to exceed that of a single diagnostic CT scan. This places the mobile C-arm in the category of extended CBCT geometries with sufficient air gap for which the tradeoffs between CNR and dose typically do not favor incorporation of an antiscatter grid.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Scatter correction in cone-beam CT via a half beam blocker technique allowing simultaneous acquisition of scatter and image information

Purpose: X-ray scatter incurred to detectors degrades the quality of cone-beam computed tomography (CBCT) and represents a problem in volumetric image guided and adaptive radiation therapy. Several methods using a beam blocker for the estimation and subtraction of scatter have been proposed. However, due to missing information resulting from the obstruction of the blocker, such methods require dual scanning or dynamically moving blocker to obtain a complete volumetric image. Here, we propose a half beam blocker-based approach, in conjunction with a total variation (TV) regularized Feldkamp-Davis-Kress (FDK) algorithm, to correct scatter-induced artifacts by simultaneously acquiring image and scatter information from a single-rotation CBCT scan. Methods: A half beam blocker, comprising lead strips, is used to simultaneously acquire image data on one side of the projection data and scatter data on the other half side. One-dimensional cubic B-Spline interpolation/extrapolation is applied to derive patient specific scatter information by using the scatter distributions on strips. The estimated scatter is subtracted from the projection image acquired at the opposite view. With scatter-corrected projections where this subtraction is completed, the FDK algorithm based on a cosine weighting function is performed to reconstruct CBCT volume. To suppress the noise in the reconstructed CBCT images produced by geometric errors between two opposed projections and interpolated scatter information, total variation regularization is applied by a minimization using a steepest gradient descent optimization method. The experimental studies using Catphan504 and anthropomorphic phantoms were carried out to evaluate the performance of the proposed scheme. Results: The scatter-induced shading artifacts were markedly suppressed in CBCT using the proposed scheme. Compared with CBCT without a blocker, the nonuniformity value was reduced from 39.3% to 3.1%. The root mean square error relative to values inside the regions of interest selected from a benchmark scatter free image was reduced from 50 to 11.3. The TV regularization also led to a better contrast-to-noise ratio. Conclusions: An asymmetric half beam blocker-based FDK acquisition and reconstruction technique has been established. The proposed scheme enables simultaneous detection of patient specific scatter and complete volumetric CBCT reconstruction without additional requirements such as prior images, dual scans, or moving strips.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantitative cone-beam CT imaging in radiation therapy using planning CT as a prior: First patient studies

Purpose: Quantitative cone-beam CT (CBCT) imaging is on increasing demand for high-performance image guided radiation therapy (IGRT). However, the current CBCT has poor image qualities mainly due to scatter contamination. Its current clinical application is therefore limited to patient setup based on only bony structures. To improve CBCT imaging for quantitative use, we recently proposed a correction method using planning CT (pCT) as the prior knowledge. Promising phantom results have been obtained on a tabletop CBCT system, using a correction scheme with rigid registration and without iterations. More challenges arise in clinical implementations of our method, especially because patients have large organ deformation in different scans. In this paper, we propose an improved framework to extend our method from bench to bedside by including several new components. Methods: The basic principle of our correction algorithm is to estimate the primary signals of CBCT projections via forward projection on the pCT image, and then to obtain the low-frequency errors in CBCT raw projections by subtracting the estimated primary signals and low-pass filtering. We improve the algorithm by using deformable registration to minimize the geometry difference between the pCT and the CBCT images. Since the registration performance relies on the accuracy of the CBCT image, we design an optional iterative scheme to update the CBCT image used in the registration. Large correction errors result from the mismatched objects in the pCT and the CBCT scans. Another optional step of gas pocket and couch matching is added into the framework to reduce these effects. Results: The proposed method is evaluated on four prostate patients, of which two cases are presented in detail to investigate the method performance for a large variety of patient geometry in clinical practice. The first patient has small anatomical changes from the planning to the treatment room. Our algorithm works well even without the optional iterations and the gas pocket and couch matching. The image correction on the second patient is more challenging due to the effects of gas pockets and attenuating couch. The improved framework with all new components is used to fully evaluate the correction performance. The enhanced image quality has been evaluated using mean CT number and spatial nonuniformity (SNU) error as well as contrast improvement factor. If the pCT image is considered as the ground truth, on the four patients, the overall mean CT number error is reduced from over 300 HU to below 16 HU in the selected regions of interest (ROIs), and the SNU error is suppressed from over 18% to below 2%. The average soft-tissue contrast is improved by an average factor of 2.6. Conclusions: We further improve our pCT-based CBCT correction algorithm for clinical use. Superior correction performance has been demonstrated on four patient studies. By providing quantitative CBCT images, our approach significantly increases the accuracy of advanced CBCT-based clinical applications for IGRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Feasibility study of a synchronized-moving-grid (SMOG) system to improve image quality in cone-beam computed tomography (CBCT)

Purpose: To evaluate the feasibility of a synchronized moving grid (SMOG) system to remove scatter artifacts, improve the contrast-to-noise ratio (CNR), and reduce image lag artifacts in cone-beam CT (CBCT). Methods: The SMOG system proposed here uses a rapidly oscillating, synchronized moving grid attached to the kV source. Multiple partial projections are taken at different grid positions to form a complete projection in each gantry position, before the gantry moves to the next position during a scan. The grid has a low transmission factor, and it is used for both scatter reduction and scatter measurement for postscan scatter correction. Experimental studies using a static grid and an enlarged CATphan phantom were performed to evaluate the potential CNR enhancement for different SMOG exposure numbers (1, 2, and 4). Simulation studies were performed to evaluate the image lag correction for different exposure numbers (2, 3, and 4) and grid interspace widths in SMOG using the data from an anthropomorphic pelvis phantom scan. Imaging dose of SMOG was also estimated by measuring the imaging dose in a CIRS CT dose phantom using a static grid. Results: SMOG can enhance the CNR by 16% and 13% when increasing exposure number from 1 to 2 and from 2 to 4, respectively. This enhancement was more dramatic for larger phantoms and smaller initial exposure numbers. Simulation results indicated that SMOG could reduce the lag to less than 4.3% for 2-exposure mode and to less than 0.8% for 3-exposure mode when the grid interspace width was 1.4 cm. Increasing the number of exposures in SMOG dramatically reduced the residual lag in the image. Reducing the grid interspace width somewhat reduced the residual lag. Skin line artifacts were removed entirely in SMOG. Point dose measurement showed that imaging dose of SMOG at isocenter was similar as that of a conventional CBCT. Conclusions: Compared to our previously developed static-grid dual-rotation method, the proposed SMOG technique has the advantages of enhancing the CNR, correcting the image lag, and reducing the delivery time. Once implemented, SMOG has the potential to remove scatter and image lag artifacts, and significantly enhance CNR for CBCT using the same scanning time as conventional CBCT.

60 APPLIED LIFE SCIENCES↗

Noise suppression in reconstruction of low-Z target megavoltage cone-beam CT images

Purpose: To improve the image contrast-to-noise (CNR) ratio for low-Z target megavoltage cone-beam CT (MV CBCT) using a statistical projection noise suppression algorithm based on the penalized weighted least-squares (PWLS) criterion. Methods: Projection images of a contrast phantom, a CatPhan{sup Registered-Sign} 600 phantom and a head phantom were acquired by a Varian 2100EX LINAC with a low-Z (Al) target and low energy x-ray beam (2.5 MeV) at a low-dose level and at a high-dose level. The projections were then processed by minimizing the PWLS objective function. The weighted least square (WLS) term models the noise of measured projection and the penalty term enforces the smoothing constraints of the projection image. The variance of projection data was chosen as the weight for the PWLS objective function and it determined the contribution of each measurement. An anisotropic quadratic form penalty that incorporates the gradient information of projection image was used to preserve edges during noise reduction. Low-Z target MV CBCT images were reconstructed by the FDK algorithm after each projection was processed by the PWLS smoothing. Results: Noise in low-Z target MV CBCT images were greatly suppressed after the PWLS projection smoothing, without noticeable sacrifice of the spatial resolution. Depending on the choice of smoothing parameter, the CNR of selected regions of interest in the PWLS processed low-dose low-Z target MV CBCT image can be higher than the corresponding high-dose image.Conclusion: The CNR of low-Z target MV CBCT images was substantially improved by using PWLS projection smoothing. The PWLS projection smoothing algorithm allows the reconstruction of high contrast low-Z target MV CBCT image with a total dose of as low as 2.3 cGy.

60 APPLIED LIFE SCIENCES↗

Four-dimensional cone beam CT reconstruction and enhancement using a temporal nonlocal means method

Purpose: Four-dimensional cone beam computed tomography (4D-CBCT) has been developed to provide respiratory phase-resolved volumetric imaging in image guided radiation therapy. Conventionally, it is reconstructed by first sorting the x-ray projections into multiple respiratory phase bins according to a breathing signal extracted either from the projection images or some external surrogates, and then reconstructing a 3D CBCT image in each phase bin independently using FDK algorithm. This method requires adequate number of projections for each phase, which can be achieved using a low gantry rotation or multiple gantry rotations. Inadequate number of projections in each phase bin results in low quality 4D-CBCT images with obvious streaking artifacts. 4D-CBCT images at different breathing phases share a lot of redundant information, because they represent the same anatomy captured at slightly different temporal points. Taking this redundancy along the temporal dimension into account can in principle facilitate the reconstruction in the situation of inadequate number of projection images. In this work, the authors propose two novel 4D-CBCT algorithms: an iterative reconstruction algorithm and an enhancement algorithm, utilizing a temporal nonlocal means (TNLM) method. Methods: The authors define a TNLM energy term for a given set of 4D-CBCT images. Minimization of this term favors those 4D-CBCT images such that any anatomical features at one spatial point at one phase can be found in a nearby spatial point at neighboring phases. 4D-CBCT reconstruction is achieved by minimizing a total energy containing a data fidelity term and the TNLM energy term. As for the image enhancement, 4D-CBCT images generated by the FDK algorithm are enhanced by minimizing the TNLM function while keeping the enhanced images close to the FDK results. A forward-backward splitting algorithm and a Gauss-Jacobi iteration method are employed to solve the problems. The algorithms implementation on GPU is designed to avoid redundant and uncoalesced memory access, in order to ensure a high computational efficiency. Our algorithms have been tested on a digital NURBS-based cardiac-torso phantom and a clinical patient case. Results: The reconstruction algorithm and the enhancement algorithm generate visually similar 4D-CBCT images, both better than the FDK results. Quantitative evaluations indicate that, compared with the FDK results, our reconstruction method improves contrast-to-noise-ratio (CNR) by a factor of 2.56-3.13 and our enhancement method increases the CNR by 2.75-3.33 times. The enhancement method also removes over 80% of the streak artifacts from the FDK results. The total computation time is 509-683 s for the reconstruction algorithm and 524-540 s for the enhancement algorithm on an NVIDIA Tesla C1060 GPU card. Conclusions: By innovatively taking the temporal redundancy among 4D-CBCT images into consideration, the proposed algorithms can produce high quality 4D-CBCT images with much less streak artifacts than the FDK results, in the situation of inadequate number of projections.

60 APPLIED LIFE SCIENCES↗

Impact of flat panel-imager veiling glare on scatter-estimation accuracy and image quality of a commercial on-board cone-beam CT imaging system

Purpose: The purposes of this study is to measure the low frequency drop (LFD) of the modulation transfer function (MTF), associated with the long tails of the detector point spread function (PSF) of an on-board flat panel imager and study its impact on cone-beam CT (CBCT) image quality and scatter measurement accuracy. Methods: Two different experimental methods were used to characterize LFD and its associated PSF of a Varian OBI flat-panel detector system: the edge response function (ERF) method and the disk transfer function (DTF) method. PSF was estimated by fitting parametric models to these measurements for four values of the applied voltage (kVp). The resultant PSF was used to demonstrate the effect of LFD on image contrast and CT number accuracy in CBCT images reconstructed from synthetic datasets, as well as, accuracy of scatter measurements with the beam-stop method. Results: The MTFs derived from the measured ERF data revealed LFDs varying from 8% (at 60 kVp) to 10.5% (at 120 kVp), while the intensity of the long PSF tails was found to increase with increasing kVp. The veiling glare line spread functions derived from the ERF and DTF methods were in excellent agreement. Uncorrected veiling glare reduced contrast and the image intensity in CBCT reconstruction, near the phantom periphery (by 67 Hounsfield units in a 20 cm-in-diameter water phantom) and (to a smaller degree) near inhomogeneities. Use of the bow-tie filter mitigated these effects. Veiling glare also resulted in about 10%-15% overestimation of the scatter-to-primary ratio when measured with the beam-stop or beam-stop array method. Conclusions: The long tails of the detector PSF were found to have a modest dependence of beam spectrum, which is reflected on the MTF curve LFD. Our findings show that uncorrected veiling glare can affect quantitative accuracy and contrast in CBCT imaging, based on flat panel imager. In addition, it results in overestimation of the scatter-to-primary ratio, measured with the beam-stop methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Target-specific optimization of four-dimensional cone beam computed tomography

Purpose: Under-sampling artifacts are a major problem in four-dimensional cone-beam computed tomography (4D-CBCT) and may compromise evaluation of target motion. Several studies have addressed scan parameter selection for minimizing under-sampling artifacts; however, the role of the target characteristics in scan parameter selection has not been investigated. In this work, the authors evaluated 4D-CBCT performance by assessing the accuracy of target motion measurements for various target sizes and motions. The results may serve as patient-specific guidelines for scan parameters selection in 4D-CBCT. Methods: The authors acquired 4D-CBCT scans of a moving phantom consisting of six water-filled sphere targets of sizes 10-37 mm, using various scan times ranging from 30 s to 3 min., setting the motion to 3-s and 6-s periods. The authors used automatic image registration to extract the target motion trajectories and evaluated these measurements for various target sizes and motions over various combinations of scan parameters including scan time, detector configuration, number of respiration phases, and reconstruction filters. Results: The most important object parameter to 4D-CBCT performance was the period of motion. Measurements for the 6-s motion were always systematically less accurate than measurements for the 3-s motion for 34 of 36 objects of various sizes and periods of motion. The 6-s motion required a greater scan time than the 3-s motion for equivalent measurement accuracy. The second most influential parameter was the target size. For the 3-s period of motion, objects larger than 13 mm were tracked with sub-millimeter accuracy with a 1-min scan time. For the 6-s period of motion, objects larger than 22-mm were tracked with sub-millimeter accuracy with a 1.5-min scan time. For all sizes and motions, temporal blurring was observed when the number of phases was fewer than 8. Offset detector configuration provided the same performance as centered detector except for small targets (20 mm.) at short scan times ({<=}1 min). Finally, reconstruction filtering and number of respiratory phases did not affect performance. Conclusions: Scan time should be set according to target size and motion. The authors have provided figures that provide the minimum scan time needed to achieve the particular motion measurement accuracy for a particular size and motion period.

60 APPLIED LIFE SCIENCES↗

Robust methods for automatic image-to-world registration in cone-beam CT interventional guidance

Purpose: Real-time surgical navigation relies on accurate image-to-world registration to align the coordinate systems of the image and patient. Conventional manual registration can present a workflow bottleneck and is prone to manual error and intraoperator variability. This work reports alternative means of automatic image-to-world registration, each method involving an automatic registration marker (ARM) used in conjunction with C-arm cone-beam CT (CBCT). The first involves a Known-Model registration method in which the ARM is a predefined tool, and the second is a Free-Form method in which the ARM is freely configurable. Methods: Studies were performed using a prototype C-arm for CBCT and a surgical tracking system. A simple ARM was designed with markers comprising a tungsten sphere within infrared reflectors to permit detection of markers in both x-ray projections and by an infrared tracker. The Known-Model method exercised a predefined specification of the ARM in combination with 3D-2D registration to estimate the transformation that yields the optimal match between forward projection of the ARM and the measured projection images. The Free-Form method localizes markers individually in projection data by a robust Hough transform approach extended from previous work, backprojected to 3D image coordinates based on C-arm geometric calibration. Image-domain point sets were transformed to world coordinates by rigid-body point-based registration. The robustness and registration accuracy of each method was tested in comparison to manual registration across a range of body sites (head, thorax, and abdomen) of interest in CBCT-guided surgery, including cases with interventional tools in the radiographic scene. Results: The automatic methods exhibited similar target registration error (TRE) and were comparable or superior to manual registration for placement of the ARM within {approx}200 mm of C-arm isocenter. Marker localization in projection data was robust across all anatomical sites, including challenging scenarios involving the presence of interventional tools. The reprojection error of marker localization was independent of the distance of the ARM from isocenter, and the overall TRE was dominated by the configuration of individual fiducials and distance from the target as predicted by theory. The median TRE increased with greater ARM-to-isocenter distance (e.g., for the Free-Form method, TRE increasing from 0.78 mm to 2.04 mm at distances of {approx}75 mm and 370 mm, respectively). The median TRE within {approx}200 mm distance was consistently lower than that of the manual method (TRE = 0.82 mm). Registration performance was independent of anatomical site (head, thorax, and abdomen). The Free-Form method demonstrated a statistically significant improvement (p= 0.0044) in reproducibility compared to manual registration (0.22 mm versus 0.30 mm, respectively). Conclusions: Automatic image-to-world registration methods demonstrate the potential for improved accuracy, reproducibility, and workflow in CBCT-guided procedures. A Free-Form method was shown to exhibit robustness against anatomical site, with comparable or improved TRE compared to manual registration. It was also comparable or superior in performance to a Known-Model method in which the ARM configuration is specified as a predefined tool, thereby allowing configuration of fiducials on the fly or attachment to the patient.

60 APPLIED LIFE SCIENCES↗

4D cone beam CT via spatiotemporal tensor framelet

Purpose: On-board 4D cone beam CT (4DCBCT) offers respiratory phase-resolved volumetric imaging, and improves the accuracy of target localization in image guided radiation therapy. However, the clinical utility of this technique has been greatly impeded by its degraded image quality, prolonged imaging time, and increased imaging dose. The purpose of this letter is to develop a novel iterative 4DCBCT reconstruction method for improved image quality, increased imaging speed, and reduced imaging dose. Methods: The essence of this work is to introduce the spatiotemporal tensor framelet (STF), a high-dimensional tensor generalization of the 1D framelet for 4DCBCT, to effectively take into account of highly correlated and redundant features of the patient anatomy during respiration, in a multilevel fashion with multibasis sparsifying transform. The STF-based algorithm is implemented on a GPU platform for improved computational efficiency. To evaluate the method, 4DCBCT full-fan scans were acquired within 30 s, with a gantry rotation of 200°; STF is also compared with a state-of-art reconstruction method via spatiotemporal total variation regularization. Results: Both the simulation and experimental results demonstrate that STF-based reconstruction achieved superior image quality. The reconstruction of 20 respiratory phases took less than 10 min on an NVIDIA Tesla C2070 GPU card. The STF codes are available at https://sites.google.com/site/spatiotemporaltensorframelet . Conclusions: By effectively utilizing the spatiotemporal coherence of the patient anatomy among different respiratory phases in a multilevel fashion with multibasis sparsifying transform, the proposed STF method potentially enables fast and low-dose 4DCBCT with improved image quality.

60 APPLIED LIFE SCIENCES↗

The effect of irregular breathing patterns on internal target volumes in four-dimensional CT and cone-beam CT images in the context of stereotactic lung radiotherapy

Purpose: Stereotactic lung radiotherapy is complicated by tumor motion from patient respiration. Four-dimensional CT (4DCT) imaging is a motion compensation method used in treatment planning to generate a maximum intensity projection (MIP) internal target volume (ITV). Image guided radiotherapy during treatment may involve acquiring a volumetric cone-beam CT (CBCT) image and visually aligning the tumor to the planning 4DCT MIP ITV contour. Moving targets imaged with CBCT can appear blurred and currently there are no studies reporting on the effect that irregular breathing patterns have on CBCT volumes and their alignment to 4DCT MIP ITV contours. The objective of this work was therefore to image a phantom moving with irregular breathing patterns to determine whether any configurations resulted in errors in volume contouring or alignment. Methods: A Perspex thorax phantom was used to simulate a patient. Three wooden 'lung' inserts with embedded Perspex 'lesions' were moved up to 4 cm with computer-generated motion patterns, and up to 1 cm with patient-specific breathing patterns. The phantom was imaged on 4DCT and CBCT with the same acquisition settings used for stereotactic lung patients in the clinic and the volumes on all phantom images were contoured. This project assessed the volumes for qualitative and quantitative changes including volume, length of the volume, and errors in alignment between CBCT volumes and 4DCT MIP ITV contours. Results: When motion was introduced 4DCT and CBCT volumes were reduced by up to 20% and 30% and shortened by up to 7 and 11 mm, respectively, indicating that volume was being under-represented at the extremes of motion. Banding artifacts were present in 4DCT MIP images, while CBCT volumes were largely reduced in contrast. When variable amplitudes from patient traces were used and CBCT ITVs were compared to 4DCT MIP ITVs there was a distinct trend in reduced ITV with increasing amplitude that was not seen when compared to true ITVs. Breathing patterns with a rest period following expiration resulted in well-defined superior edges and were better aligned using an edge-to-edge alignment technique. In most cases, sinusoidal motion patterns resulted in the closest agreements to true values and the smallest misalignments. Conclusions: Strategies are needed to compensate for volume losses at the extremes of motion for both 4DCT MIP and CBCT images for larger and varied amplitudes, and for patterns with rest periods following expiration. Lesions moving greater than 2 cm would warrant larger treatment margins added to the 4DCT MIP ITV to account for the volume being under-represented at the extremes of motion. Lesions moving with a rest period following expiration would be better aligned using an edge-to-edge alignment technique. Sinusoidal patterns represented the ideal clinical scenario, reinforcing the importance of investigating clinically relevant motions and their effects on 4DCT MIP and CBCT volumes. Since most patients do not breathe sinusoidally this may lead to misinterpretation of previous studies using only sinusoidal motion.

60 APPLIED LIFE SCIENCES↗

A dual cone-beam CT system for image guided radiotherapy: Initial performance characterization

Purpose: The purpose of this study is to evaluate the performance of a recently developed benchtop dual cone-beam computed tomography (CBCT) system with two orthogonally placed tube/detector sets. Methods: The benchtop dual CBCT system consists of two orthogonally placed 40 Multiplication-Sign 30 cm flat-panel detectors and two conventional x-ray tubes with two individual high-voltage generators sharing the same rotational axis. The x-ray source to detector distance is 150 cm and x-ray source to rotational axis distance is 100 cm for both subsystems. The objects are scanned through 200 Degree-Sign of rotation. The dual CBCT system utilized 110 Degree-Sign of projection data from one detector and 90 Degree-Sign from the other while the two individual single CBCTs utilized 200 Degree-Sign data from each detector. The system performance was characterized in terms of uniformity, contrast, spatial resolution, noise power spectrum, and CT number linearity. The uniformities, within the axial slice and along the longitudinal direction, and noise power spectrum were assessed by scanning a water bucket; the contrast and CT number linearity were measured using the Catphan phantom; and the spatial resolution was evaluated using a tungsten wire phantom. A skull phantom and a ham were also scanned to provide qualitative evaluation of high- and low-contrast resolution. Each measurement was compared between dual and single CBCT systems. Results: Compared to single CBCT, the dual CBCT presented: (1) a decrease in uniformity by 1.9% in axial view and 1.1% in the longitudinal view, as averaged for four energies (80, 100, 125, and 150 kVp); (2) comparable or slightly better contrast (0{approx}25 HU) for low-contrast objects and comparable contrast for high-contrast objects; (3) comparable spatial resolution; (4) comparable CT number linearity with R{sup 2}{>=} 0.99 for all four tested energies; (5) lower noise power spectrum in magnitude. Dual CBCT images of the skull phantom and the ham demonstrated both high-contrast resolution and good soft-tissue contrast. Conclusions: The performance of a benchtop dual CBCT imaging system has been characterized and is comparable to that of a single CBCT.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Respiratory triggered 4D cone-beam computed tomography: A novel method to reduce imaging dose

Purpose: A novel method called respiratory triggered 4D cone-beam computed tomography (RT 4D CBCT) is described whereby imaging dose can be reduced without degrading image quality. RT 4D CBCT utilizes a respiratory signal to trigger projections such that only a single projection is assigned to a given respiratory bin for each breathing cycle. In contrast, commercial 4D CBCT does not actively use the respiratory signal to minimize image dose. Methods: To compare RT 4D CBCT with conventional 4D CBCT, 3600 CBCT projections of a thorax phantom were gathered and reconstructed to generate a ground truth CBCT dataset. Simulation pairs of conventional 4D CBCT acquisitions and RT 4D CBCT acquisitions were developed assuming a sinusoidal respiratory signal which governs the selection of projections from the pool of 3600 original projections. The RT 4D CBCT acquisition triggers a single projection when the respiratory signal enters a desired acquisition bin; the conventional acquisition does not use a respiratory trigger and projections are acquired at a constant frequency. Acquisition parameters studied were breathing period, acquisition time, and imager frequency. The performance of RT 4D CBCT using phase based and displacement based sorting was also studied. Image quality was quantified by calculating difference images of the test dataset from the ground truth dataset. Imaging dose was calculated by counting projections. Results: Using phase based sorting RT 4D CBCT results in 47% less imaging dose on average compared to conventional 4D CBCT. Image quality differences were less than 4% at worst. Using displacement based sorting RT 4D CBCT results in 57% less imaging dose on average, than conventional 4D CBCT methods; however, image quality was 26% worse with RT 4D CBCT. Conclusions: Simulation studies have shown that RT 4D CBCT reduces imaging dose while maintaining comparable image quality for phase based 4D CBCT; image quality is degraded for displacement based RT 4D CBCT in its current implementation.

60 APPLIED LIFE SCIENCES↗

Actively triggered 4d cone-beam CT acquisition

Purpose: 4d cone-beam computed tomography (CBCT) scans are usually reconstructed by extracting the motion information from the 2d projections or an external surrogate signal, and binning the individual projections into multiple respiratory phases. In this “after-the-fact” binning approach, however, projections are unevenly distributed over respiratory phases resulting in inefficient utilization of imaging dose. To avoid excess dose in certain respiratory phases, and poor image quality due to a lack of projections in others, the authors have developed a novel 4d CBCT acquisition framework which actively triggers 2d projections based on the forward-predicted position of the tumor.Methods: The forward-prediction of the tumor position was independently established using either (i) an electromagnetic (EM) tracking system based on implanted EM-transponders which act as a surrogate for the tumor position, or (ii) an external motion sensor measuring the chest-wall displacement and correlating this external motion to the phase-shifted diaphragm motion derived from the acquired images. In order to avoid EM-induced artifacts in the imaging detector, the authors devised a simple but effective “Faraday” shielding cage. The authors demonstrated the feasibility of their acquisition strategy by scanning an anthropomorphic lung phantom moving on 1d or 2d sinusoidal trajectories.Results: With both tumor position devices, the authors were able to acquire 4d CBCTs free of motion blurring. For scans based on the EM tracking system, reconstruction artifacts stemming from the presence of the EM-array and the EM-transponders were greatly reduced using newly developed correction algorithms. By tuning the imaging frequency independently for each respiratory phase prior to acquisition, it was possible to harmonize the number of projections over respiratory phases. Depending on the breathing period (3.5 or 5 s) and the gantry rotation time (4 or 5 min), between ∼90 and 145 projections were acquired per respiratory phase resulting in a dose of ∼1.7–2.6 mGy per respiratory phase. Further dose savings and decreases in the scanning time are possible by acquiring only a subset of all respiratory phases, for example, peak-exhale and peak-inhale only scans.Conclusions: This study is the first experimental demonstration of a new 4d CBCT acquisition paradigm in which imaging dose is efficiently utilized by actively triggering only those projections that are desired for the reconstruction process.

61 RADIATION PROTECTION AND DOSIMETRY↗

Simultaneous motion estimation and image reconstruction (SMEIR) for 4D cone-beam CT

Purpose: Image reconstruction and motion model estimation in four-dimensional cone-beam CT (4D-CBCT) are conventionally handled as two sequential steps. Due to the limited number of projections at each phase, the image quality of 4D-CBCT is degraded by view aliasing artifacts, and the accuracy of subsequent motion modeling is decreased by the inferior 4D-CBCT. The objective of this work is to enhance both the image quality of 4D-CBCT and the accuracy of motion model estimation with a novel strategy enabling simultaneous motion estimation and image reconstruction (SMEIR).Methods: The proposed SMEIR algorithm consists of two alternating steps: (1) model-based iterative image reconstruction to obtain a motion-compensated primary CBCT (m-pCBCT) and (2) motion model estimation to obtain an optimal set of deformation vector fields (DVFs) between the m-pCBCT and other 4D-CBCT phases. The motion-compensated image reconstruction is based on the simultaneous algebraic reconstruction technique (SART) coupled with total variation minimization. During the forward- and backprojection of SART, measured projections from an entire set of 4D-CBCT are used for reconstruction of the m-pCBCT by utilizing the updated DVF. The DVF is estimated by matching the forward projection of the deformed m-pCBCT and measured projections of other phases of 4D-CBCT. The performance of the SMEIR algorithm is quantitatively evaluated on a 4D NCAT phantom. The quality of reconstructed 4D images and the accuracy of tumor motion trajectory are assessed by comparing with those resulting from conventional sequential 4D-CBCT reconstructions (FDK and total variation minimization) and motion estimation (demons algorithm). The performance of the SMEIR algorithm is further evaluated by reconstructing a lung cancer patient 4D-CBCT.Results: Image quality of 4D-CBCT is greatly improved by the SMEIR algorithm in both phantom and patient studies. When all projections are used to reconstruct a 3D-CBCT by FDK, motion-blurring artifacts are present, leading to a 24.4% relative reconstruction error in the NACT phantom. View aliasing artifacts are present in 4D-CBCT reconstructed by FDK from 20 projections, with a relative error of 32.1%. When total variation minimization is used to reconstruct 4D-CBCT, the relative error is 18.9%. Image quality of 4D-CBCT is substantially improved by using the SMEIR algorithm and relative error is reduced to 7.6%. The maximum error (MaxE) of tumor motion determined from the DVF obtained by demons registration on a FDK-reconstructed 4D-CBCT is 3.0, 2.3, and 7.1 mm along left–right (L-R), anterior–posterior (A-P), and superior–inferior (S-I) directions, respectively. From the DVF obtained by demons registration on 4D-CBCT reconstructed by total variation minimization, the MaxE of tumor motion is reduced to 1.5, 0.5, and 5.5 mm along L-R, A-P, and S-I directions. From the DVF estimated by SMEIR algorithm, the MaxE of tumor motion is further reduced to 0.8, 0.4, and 1.5 mm along L-R, A-P, and S-I directions, respectively.Conclusions: The proposed SMEIR algorithm is able to estimate a motion model and reconstruct motion-compensated 4D-CBCT. The SMEIR algorithm improves image reconstruction accuracy of 4D-CBCT and tumor motion trajectory estimation accuracy as compared to conventional sequential 4D-CBCT reconstruction and motion estimation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Development and validation of a measurement-based source model for kilovoltage cone-beam CT Monte Carlo dosimetry simulations

Purpose: The purpose of this study is to adapt an equivalent source model originally developed for conventional CT Monte Carlo dose quantification to the radiation oncology context and validate its application for evaluating concomitant dose incurred by a kilovoltage (kV) cone-beam CT (CBCT) system integrated into a linear accelerator.Methods: In order to properly characterize beams from the integrated kV CBCT system, the authors have adapted a previously developed equivalent source model consisting of an equivalent spectrum module that takes into account intrinsic filtration and an equivalent filter module characterizing the added bowtie filtration. An equivalent spectrum was generated for an 80, 100, and 125 kVp beam with beam energy characterized by half-value layer measurements. An equivalent filter description was generated from bowtie profile measurements for both the full- and half-bowtie. Equivalent source models for each combination of equivalent spectrum and filter were incorporated into the Monte Carlo software package MCNPX. Monte Carlo simulations were then validated against in-phantom measurements for both the radiographic and CBCT mode of operation of the kV CBCT system. Radiographic and CBCT imaging dose was measured for a variety of protocols at various locations within a body (32 cm in diameter) and head (16 cm in diameter) CTDI phantom. The in-phantom radiographic and CBCT dose was simulated at all measurement locations and converted to absolute dose using normalization factors calculated from air scan measurements and corresponding simulations. The simulated results were compared with the physical measurements and their discrepancies were assessed quantitatively.Results: Strong agreement was observed between in-phantom simulations and measurements. For the radiographic protocols, simulations uniformly underestimated measurements by 0.54%–5.14% (mean difference =−3.07%, SD = 1.60%). For the CBCT protocols, simulations uniformly underestimated measurements by 1.35%–5.31% (mean difference =−3.42%, SD = 1.09%).Conclusions: This work demonstrates the feasibility of using a measurement-based kV CBCT source model to facilitate dose calculations with Monte Carlo methods for both the radiographic and CBCT mode of operation. While this initial work validates simulations against measurements for simple geometries, future work will involve utilizing the source model to investigate kV CBCT dosimetry with more complex anthropomorphic phantoms and patient specific models.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Motion-map constrained image reconstruction (MCIR): Application to four-dimensional cone-beam computed tomography

Purpose: Utilization of respiratory correlated four-dimensional cone-beam computed tomography (4DCBCT) has enabled verification of internal target motion and volume immediately prior to treatment. However, with current standard CBCT scan, 4DCBCT poses challenge for reconstruction due to the fact that multiple phase binning leads to insufficient number of projection data to reconstruct and thus cause streaking artifacts. The purpose of this study is to develop a novel 4DCBCT reconstruction algorithm framework called motion-map constrained image reconstruction (MCIR), that allows reconstruction of high quality and high phase resolution 4DCBCT images with no more than the imaging dose as well as projections used in a standard free breathing 3DCBCT (FB-3DCBCT) scan.Methods: The unknown 4DCBCT volume at each phase was mathematically modeled as a combination of FB-3DCBCT and phase-specific update vector which has an associated motion-map matrix. The motion-map matrix, which is the key innovation of the MCIR algorithm, was defined as the matrix that distinguishes voxels that are moving from stationary ones. This 4DCBCT model was then reconstructed with compressed sensing (CS) reconstruction framework such that the voxels with high motion would be aggressively updated by the phase-wise sorted projections and the voxels with less motion would be minimally updated to preserve the FB-3DCBCT. To evaluate the performance of our proposed MCIR algorithm, we evaluated both numerical phantoms and a lung cancer patient. The results were then compared with the (1) clinical FB-3DCBCT reconstructed using the FDK, (2) 4DCBCT reconstructed using the FDK, and (3) 4DCBCT reconstructed using the well-known prior image constrained compressed sensing (PICCS).Results: Examination of the MCIR algorithm showed that high phase-resolved 4DCBCT with sets of up to 20 phases using a typical FB-3DCBCT scan could be reconstructed without compromising the image quality. Moreover, in comparison with other published algorithms, the image quality of the MCIR algorithm is shown to be excellent.Conclusions: This work demonstrates the potential for providing high-quality 4DCBCT during on-line image-guided radiation therapy (IGRT), without increasing the imaging dose. The results showed that (at least) 20 phase images could be reconstructed using the same projections data, used to reconstruct a single FB-3DCBCT, without streak artifacts that are caused by insufficient projections.

62 RADIOLOGY AND NUCLEAR MEDICINE↗