Model refinement increases confidence levels and clinical agreement when commissioning a three-dimensional secondary dose calculation system.

PURPOSE: Evaluate custom beam models for a second check dose calculation system using statistically verifiable passing criteria for film analysis, DVH, and 3D gamma metrics. METHODS: Custom beam models for nine linear accelerators for the Sun Nuclear Dose Calculator algorithm (SDC, Sun Nuclear) were evaluated using the AAPM-TG119 test suite (5 Intensity Modulated Radiation Therapy (IMRT) and 5 Volumetric Modulated Arc Therapy (VMAT) plans) and a set of clinical plans. Where deemed necessary, adjustments to Multileaf Collimator (MLC) parameters were made to improve results. Comparisons to the Analytic Anisotropic Algorithm (AAA), and gafchromic film measurements were performed. Confidence intervals were set to 95% per TG-119. Film gamma criteria were 3%/3 mm (conventional beams) or 3%/1 mm (Stereotactic Radiosurgery [SRS] beams). Dose distributions in solid water phantom were evaluated based on DVH metrics (e.g., D95, V20) and 3D gamma criteria (3%/3 mm or 3%/1 mm). Film passing rates, 3D gamma passing rates, and DVH metrics were reported for HD MLC machines and Millennium MLC Machines. RESULTS: For HD MLC machines, SDC gamma film agreement was 98.76% ± 2.30% (5.74% CL) for 6FFF/6srs (3%/1 mm), and 99.80% ± 0.32% (0.83% CL) for 6x (3%/3 mm). For Millennium MLC machines, film passing rates were 98.20% ± 3.14% (7.96% CL), 99.52% ± 1.14% (2.71% CL), and 99.69% ± 0.82% (1.91% CL) for 6FFF, 6x, and 10x, respectively. For SDC to AAA comparisons: HD MLC Linear Accelerators (LINACs); DVH point agreement was 0.97% ± 1.64% (4.18% CL) and 1.05% ± 2.12% (5.20% CL); 3D gamma agreement was 99.97% ± 0.14% (0.30% CL) and 100.00% ± 0.02% (0.05% CL), for 6FFF/6srs and 6x, respectively; Millennium MLC LINACs: DVH point agreement was 0.77% ± 2.40% (5.47% CL), 0.80% ± 3.40% (7.47% CL), and 0.07% ± 2.15% (4.30% CL); 3D gamma agreement was 99.97% ± 0.13% (0.29% CL), 99.97% ± 0.17% (0.36% CL), and 99.99% ± 0.06% (0.12% CL) for 6FFF, 6x, and 10x, respectively. CONCLUSION: SDC shows agreement well within TG119 CLs for film and redundant dose calculation comparisons with AAA. In some models (SRS), this was achieved using stricter criteria. TG119 plans can be used to help guide model adjustments and to establish clinical baselines for DVH and 3D gamma criteria.

A unified scatter rejection and correction method for cone beam computed tomography.

PURPOSE: Scattered radiation is a major cause of image quality degradation in flat panel detector-based cone beam CT (CBCT). While recently introduced 2D antiscatter grids reject the majority of scatter fluence, the small percentage of scatter fluence still transmitted to the detector remains a major challenge for implementation of quantitative imaging techniques such as dual energy imaging in CBCT. Additionally, this residual scatter is also a major source of grid-induced artifacts, which impedes implementation of 2D grids in CBCT. We therefore present a new method to achieve both robust scatter rejection and residual scatter correction using a 2D antiscatter grid; in doing so, we expand the role of 2D grids from mere scatter rejection devices to scatter measurement devices. METHOD: In our method, the radiopaque septa of the 2D grid emulate a micro array of beam-stops placed on the detector which introduce spatially periodic septal shadows. By selecting sufficiently thin grid septa, the primary intensity can be reduced while preserving the uniformity of scatter intensity. This enables us to correlate the modulated pixel signal intensity in septal shadows with local scatter intensity. Our method then exploits this correlation to measure and remove residual scatter intensity from projections. No assumptions are made about the object being imaged. We refer to this as Grid-based Scatter Sampling (GSS). In this work, we evaluate the principle of signal modulation with grid septa, the accuracy of scatter estimates, and the effect of the GSS method on image quality using simulations and measurements. We also implement the GSS method experimentally using a 2D grid prototype. RESULTS: Our results demonstrate that the GSS method increased CT number accuracy and reduced image artifacts associated with scatter. With 2D grid and residual scatter correction, HU nonuniformity was reduced from 65 HU to 30 HU in pelvis sized phantoms, and HU variations due to change in phantom size were reduced from 59 HU to 20 HU, when compared to use of only a 2D grid. With residual scatter correction via GSS method, grid-induced ring artifacts were suppressed, leading to a 41% reduction in noise. The shape of the modulation transfer function (MTF) was preserved before and after suppression of ring artifacts. CONCLUSIONS: Our grid-based scatter sampling method enables utilization of a 2D grid as a scatter measurement and correction device. This method significantly improves quantitative accuracy in CBCT, further reducing the image quality gap between CBCT and multi-detector CT. By correcting residual scatter with the proposed method, grid-induced line artifacts in projections and associated ring artifacts in CBCT images were also suppressed with no compromise of spatial resolution.

Towards a theory of a solution space for the biplane imaging geometry problem.

Biplane angiographic imaging is a primary method for visual and quantitative assessment of the vasculature. In order to reliably reconstruct the three-dimensional (3D) position, orientation, and shape of the vessel structure, a key problem is to determine the rotation matrix R and the translation vector t which relate the two coordinate systems. This so-called Imaging Geometry Determination problem is well studied in the medical imaging and computer vision communities and a number of interesting approaches have been reported. Each such technique determines a solution which yields 3D vasculature reconstructions with errors comparable to other techniques. From the literature, we see that different techniques with different optimization strategies yield reconstructions with equivalent errors. We have investigated this behavior, and it appears that the error in the input data leads to this equivalence effectively yielding what we call the solution space of feasible geometries, i.e., geometries which could be solutions given the error or uncertainty in the input image data. In this paper, we lay the theoretical framework for this concept of a solution space of feasible geometries using simple schematic constructions, deriving the underlying mathematical relationships, presenting implementation details, and discussing implications and applications of the proposed idea. Because the solution space of feasible geometries encompasses equivalent solutions given the input error, the solution space approach can be used to evaluate the precision of calculated geometries or 3D data based on known or estimated uncertainties in the input image data. We also use the solution space approach to calculate an imaging geometry, i.e., a solution.

Vessel size measurements in angiograms: manual measurements.

Vessel size measurement is perhaps the most often performed quantitative analysis in diagnostic and interventional angiography. Although automated vessel sizing techniques are generally considered to have good accuracy and precision, we have observed that clinicians rarely use these techniques in standard clinical practice, choosing to indicate the edges of vessels and catheters to determine sizes and calibrate magnifications, i.e., manual measurements. Thus, we undertook an investigation of the accuracy and precision of vessel sizes calculated from manually indicated edges of vessels. Manual measurements were performed by three neuroradiologists and three physicists. Vessel sizes ranged from 0.1-3.0 mm in simulation studies and 0.3-6.4 mm in phantom studies. Simulation resolution functions had full-widths-at-half-maximum (FWHM) ranging from 0.0 to 0.5 mm. Phantom studies were performed with 4.5 in., 6 in., 9 in., and 12 in. image intensifier modes, magnification factor = 1, with and without zooming. The accuracy and reproducibility of the measurements ranged from 0.1 to 0.2 mm, depending on vessel size, resolution, and pixel size, and zoom. These results indicate that manual measurements may have accuracies comparable to automated techniques for vessels with sizes greater than 1 mm, but that automated techniques which take into account the resolution function should be used for vessels with sizes smaller than 1 mm.

Vessel size measurements in angiograms: a comparison of techniques.

As interventional procedures become more complicated, the need for accurate quantitative vascular information increases. In response to this need, many commercial vendors provide techniques for measurement of vessel sizes, usually based on derivative techniques. In this study, we investigate the accuracy of several techniques used in the measurement of vessel size. Simulated images of vessels having circular cross sections were generated and convolved with various focal spot distributions taking into account the magnification. These vessel images were then convolved with Gaussian image detector line spread functions (LSFs). Additionally, images of a phantom containing vessels with a range of diameters were acquired for the 4.5", 6", 9", and 12" modes of an image intensifier-TV (II-TV) system. Vessel sizes in the images were determined using a first-derivative technique, a second-derivative technique, a linear combination of these two measured sizes, a thresholding technique, a densitometric technique, and a model-based technique. For the same focal spot size, the shape of the focal spot distribution does not affect measured vessel sizes except at large magnifications. For vessels with diameters larger than the full-width-at-half-maximum (FWHM) of the LSF, accurate vessel sizes (errors approximately 0.1 mm) could be obtained by using an average of sizes determined by the first and second derivatives. For vessels with diameters smaller than the FWHM of the LSF, the densitometric and model-based techniques can provide accurate vessel sizes when these techniques are properly calibrated.