High-DQE EPIDs based on thick, segmented BGO and CsI:Tl scintillators: performance evaluation at extremely low dose

Full modulation transfer functions of thick parallel- and focused-element scintillator arrays obtained by a Monte Carlo optical transport model

BACKGROUND

Arrays of thick segmented crystalline scintillators are useful x-ray converters for image-guided radiation therapy using electronic portal imaging (EPI) and megavoltage cone-beam computed tomography (MV-CBCT). Ionizing-radiation-only simulations previously showed relatively low modulation transfer function (MTF) in parallel-element arrays because of beam divergence. Hence, a focused-element geometry (matching the beam divergence) has been proposed. The "full" (ionizing and optical) MTF performance of such a focused geometry compared to its radiation-only MTF has, however, not been fully investigated.

PURPOSE

To study the full MTF performance of such arrays in a more realistic situation in which optical characteristics are also included using an in-house detector model that supports light transport, and quantify the errors in MTF estimation when the optical stage is ignored.

METHODS

First, radiation (x-ray and electron) transport was simulated. Then, transport of the generated optical photons was modeled using ScintSim2, an optical Monte Carlo (MC) code developed in MATLAB for simulation of two-dimensional (2D) parallel- and focused-element scintillator arrays. The full-MTF responses of focused- and parallel-element geometries, for a large array of 3 × 3 mm² CsI:Tl detector elements of 10, 40, and 60 mm thicknesses, were examined. For each configuration, a composite line spread function (LSF) was calculated to obtain the MTF.

RESULTS

At the Nyquist frequency, for 10 mm-thick central elements and 60 mm-thick peripheral parallel elements, full-MTF exhibited a drop of up to 15 and 79 times, respectively, compared with radiation-only MTF. This was found to be partly attributable to the angular distribution of the light emerging from the detector-element exit face and the dependence on its aspect ratio, since the light exiting thicker scintillators exhibited a more forward-directed distribution. Focused elements provided an increase of up to nine times in peripheral-area full MTF values.

CONCLUSIONS

Full MTF was up to 79 times lower than radiation-only MTF. Focused arrays preserved full MTF by up to nine times compared to parallel elements. The differences in the results obtained with and without inclusion of optical photons emphasize the need to include light transport when optimizing thick segmented scintillation detectors. Besides their application in detector optimization for radiotherapy megavoltage photon imaging, these findings can also be useful for other segmented-scintillator-based imaging systems, for example, in nuclear medicine, or in 2D detection systems for quality assurance of MR-linacs.

Development of a novel high quantum efficiency MV x-ray detector for image-guided radiotherapy: A feasibility study

Purpose

To develop a new scintillating fiber‐based electronic portal imaging device (EPID) with a high quantum efficiency (QE) while preserving an adequate spatial resolution.

Methods

Two prototypes were built: one with a single pixel readout and the other with an active matrix flat‐panel imager (AMFPI) for readout. The energy conversion layer of both prototypes was made of scintillating fiber layers interleaved with corrugated lead sheets to form a honeycomb pattern. The scintillating fibers have a diameter of 1 mm and the distance between the centers of neighboring fibers on the same layer is 1.35 mm. The layers have 1.22 mm spacing between them. The energy conversion layer has a thickness of 2 cm. The modulation transfer function (MTF), antiscatter properties and sensitivity of the detector with a single pixel readout were measured using a 6‐MV beam on a LINAC machine. In addition, a Monte Carlo simulation was conducted to calculate the zero‐frequency detective quantum efficiency (DQE(0)) of the proposed detector with an active matrix flat‐panel imager for readout.

Results

The DQE(0) of the proposed detector can be 11.5%, which is about an order of magnitude higher than that of current EPIDs. The frequency of 50% modulation (f50) of the measured MTF is 0.2 mm-1 at 6 MV, which is comparable to that of video‐based EPIDs. The scatter to primary ratio (SPR) measured with the detector at 10 cm air gap and 20 × 20 cm2 field size is approximately 30% lower than that of ionization chamber–based detectors with a comparable QE. The detector noise which includes the x‐ray quantum noise and absorption noise is much larger than the electronic noise per pixel of the flat‐panel imager at a dose of less than two LINAC pulses. Thus, the proposed detector is quantum noise limited down to very low doses (∼a couple of radiation pulses of the LINAC). A proof‐of‐concept image has been obtained using a 6‐MV beam.

Conclusions

This work indicates that by using scintillating fibers and lead layers it is possible to increase the thickness of the detecting materials, and therefore the QE or the DQE(0) of the detector, while maintaining an adequate spatial resolution for MV x‐ray imaging. Due to the use of lead as the spacing material, the new detector also has antiscatter property, which will help improve the signal‐to‐noise ratio of the images. Further investigation to optimize the design of the detector and achieve a better combination of DQE and spatial resolution is warranted.

Characterizing a novel scintillating glass for application to megavoltage cone-beam computed tomography

PURPOSE

The purpose of this study was to evaluate the performance of a prototype electric portal imaging device (EPID) with a high detective quantum efficiency (DQE) scintillator, LKH-5. Specifically, image quality in context of both planar and megavoltage (MV) cone-beam computed tomography (CBCT) is analyzed.

METHODS

Planar image quality in terms of modulation transfer function (MTF), noise power spectrum (NPS), and DQE are measured and compared to an existing EPID (AS-1200) using the 6 MV beamline for a Varian TrueBeam linac. Imager performance is contextualized for three-dimensional (3D), MV-CBCT performance by measuring imager lag and analyzing the expected degradation of the DQE as a function of dose. Finally, comparisons between reconstructed images of the Catphan phantom in terms of qualitative quality and signal-difference-to-noise ratio (SDNR) are made for 6 MV images using both conventional and LKH-5 EPIDs as well as for the kilovoltage (kV) on-board imager (OBI).

RESULTS

Analysis of the NPS reveals linearity at all measured doses using the prototype LKH-5 detector. While the first zero of the MTF is much lower for the LKH-5 detector than the conventional EPID (0.6 cycles/mm vs 1.6 cycles/mm), the normalized NPS (NNPS) multiplied by total quanta (qNNPS) of the LKH-5 detector is roughly a factor of seven to eight times lower, yielding a DQE(0) of approximately 8%. First, second, and third frame lag were measured at approximately 23%, 5%, and 1%, respectively, although no noticeable image artifacts were apparent in reconstructed volumes. Analysis of low-dose performance reveals that DQE(0) remains at 80% of its maximum value at a dose as low as 7.5 × 10^-6  MU. For a 400 projection technique, this represents a total scan dose of 0.0030 MU, suggesting that if imaging doses are increased to a value typical of kV-CBCT scans (~2.7 cGy), the LKH-5 detector will retain quantum noise limited performance. Finally, comparing Catphan scans, the prototype detector exhibits much lower image noise than the conventional EPID, resulting in improved small object representation. Furthermore, SDNR of H₂ O and polystyrene cylinders improved from -1.95 and 2.94 to -15 and 18.7, respectively.

CONCLUSIONS

Imaging performance of the prototype LKH-5 detector was measured and analyzed for both planar and 3D contexts. Improving noise transfer of the detector results in concurrent improvement of DQE(0). For 3D imaging, temporal characteristics were adequate for artifact-free performance and at relevant doses, the detector retained quantum noise limited performance. Although quantitative MTF measurements suggest poorer resolution, small object representation of the prototype imager is qualitatively improved over the conventional detector due to the measured reduction in noise.

A novel multilayer MV imager computational model for component optimization

PURPOSE

A novel Megavoltage (MV) multilayer imager (MLI) design featuring higher detective quantum efficiency and lower noise than current conventional MV imagers in clinical use has been recently reported. Optimization of the MLI design for multiple applications including tumor tracking, MV-CBCT and portal dosimetry requires a computational model that will provide insight into the physics processes that affect the overall and individual components' performance. The purpose of the current work was to develop and validate a comprehensive computational model that can be used for MLI optimization.

METHODS

The MLI model was built using the Geant4 Application for Tomographic Emission (GATE) application. The model includes x-ray and charged-particle interactions as well as the optical transfer within the phosphor. A first prototype MLI device featuring a stack of four detection layers was used for model validation. Each layer of the prototype contains a copper buildup plate, a phosphor screen and photodiode array. The model was validated against measured data of Modulation Transfer Function (MTF), Noise-Power Spectrum (NPS), and Detective Quantum Efficiency (DQE). MTF was computed using a slanted slit with 2.3^° angle and 0.1 mm width. NPS was obtained using the autocorrelation function technique. DQE was calculated from MTF and NPS data. The comparison metrics between simulated and measured data were the Pearson's correlation coefficient (r) and the normalized root-mean-square error (NRMSE).

RESULTS

Good agreement between measured and simulated MTF and NPS values was observed. Pearson's correlation coefficient for the combined signal from all layers of the MLI was equal to 0.9991 for MTF and 0.9992 for NPS; NRMSE was 0.0121 for MTF and 0.0194 for NPS. Similarly, the DQE correlation coefficient for the combined signal was 0.9888 and the NRMSE was 0.0686.

CONCLUSIONS

A comprehensive model of the novel MLI design was developed using the GATE toolkit and validated against measured MTF, NPS, and DQE data acquired with a prototype device featuring four layers. This model will be used for further optimization of the imager components and configuration for clinical radiotherapy applications.

IGRT and motion management during lung SBRT delivery

Patient motion can cause misalignment of the tumour and toxicities to the healthy lung tissue during lung stereotactic body radiation therapy (SBRT). Any deviations from the reference setup can miss the target and have acute toxic effects on the patient with consequences onto its quality of life and survival outcomes. Correction for motion, either immediately prior to treatment or intra-treatment, can be realized with image-guided radiation therapy (IGRT) and motion management devices. The use of these techniques has demonstrated the feasibility of integrating complex technology with clinical linear accelerator to provide a higher standard of care for the patients and increase their quality of life.

Super-resolution imaging in a multiple layer EPID

A new portal imager consisting of four vertically stacked conventional electronic portal imaging device (EPID) layers has been constructed in pursuit of improved detective quantum efficiency (DQE). We hypothesize that super-resolution (SR) imaging can also be achieved in such a system by shifting each layer laterally by half a pixel relative to the layer above. Super-resolution imaging will improve resolution and contrast-to-noise ratio (CNR) in megavoltage (MV) planar and cone beam computed tomography (MV-CBCT) applications. Simulations are carried out to test this hypothesis with digital phantoms. To assess planar resolution, 2 mm long iron rods with 0.3 × 0.3 mm² square cross-section are arranged in a grid pattern at the center of a 1 cm thick solid water. For measuring CNR in MV-CBCT, a 20 cm diameter digital phantom with 8 inserts of different electron densities is used. For measuring resolution in MV-CBCT, a digital phantom featuring a bar pattern similar to the Gammex™ phantom is used. A 6 MV beam is attenuated through each phantom and detected by each of the four detector layers. Fill factor of the detector is explicitly considered. Projections are blurred with an estimated point spread function (PSF) before super-resolution reconstruction. When projections from multiple shifted layers are used in SR reconstruction, even a simple shift-add fusion can significantly improve the resolution in reconstructed images. In the reconstructed planar image, the grid pattern becomes visually clearer. In MV-CBCT, combining projections from multiple layers results in increased CNR and resolution. The inclusion of two, three and four layers increases CNR by 40%, 70% and 99%, respectively. Shifting adjacent layers by half a pixel almost doubles resolution. In comparison, using four perfectly aligned layers does not improve resolution relative to a single layer.

Optimization of the design of thick, segmented scintillators for megavoltage cone-beam CT using a novel, hybrid modeling technique

PURPOSE

Active matrix flat-panel imagers (AMFPIs) incorporating thick, segmented scintillators have demonstrated order-of-magnitude improvements in detective quantum efficiency (DQE) at radiotherapy energies compared to systems based on conventional phosphor screens. Such improved DQE values facilitate megavoltage cone-beam CT (MV CBCT) imaging at clinically practical doses. However, the MV CBCT performance of such AMFPIs is highly dependent on the design parameters of the scintillators. In this paper, optimization of the design of segmented scintillators was explored using a hybrid modeling technique which encompasses both radiation and optical effects.

METHODS

Imaging performance in terms of the contrast-to-noise ratio (CNR) and spatial resolution of various hypothetical scintillator designs was examined through a hybrid technique involving Monte Carlo simulation of radiation transport in combination with simulation of optical gain distributions and optical point spread functions. The optical simulations employed optical parameters extracted from a best fit to measurement results reported in a previous investigation of a 1.13 cm thick, 1016 μm pitch prototype BGO segmented scintillator. All hypothetical designs employed BGO material with a thickness and element-to-element pitch ranging from 0.5 to 6 cm and from 0.508 to 1.524 mm, respectively. In the CNR study, for each design, full tomographic scans of a contrast phantom incorporating various soft-tissue inserts were simulated at a total dose of 4 cGy.

RESULTS

Theoretical values for contrast, noise, and CNR were found to be in close agreement with empirical results from the BGO prototype, strongly supporting the validity of the modeling technique. CNR and spatial resolution for the various scintillator designs demonstrate complex behavior as scintillator thickness and element pitch are varied--with a clear trade-off between these two imaging metrics up to a thickness of ~3 cm. Based on these results, an optimization map indicating the regions of design that provide a balance between these metrics was obtained. The map shows that, for a given set of optical parameters, scintillator thickness and pixel pitch can be judiciously chosen to maximize performance without resorting to thicker, more costly scintillators.

CONCLUSIONS

Modeling radiation and optical effects in thick, segmented scintillators through use of a hybrid technique can provide a practical way to gain insight as to how to optimize the performance of such devices in radiotherapy imaging. Assisted by such modeling, the development of practical designs should greatly facilitate low-dose, soft tissue visualization employing MV CBCT imaging in external beam radiotherapy.

Metal artifact correction for x-ray computed tomography using kV and selective MV imaging

PURPOSE

The overall goal of this work is to improve the computed tomography (CT) image quality for patients with metal implants or fillings by completing the missing kilovoltage (kV) projection data with selectively acquired megavoltage (MV) data that do not suffer from photon starvation. When both of these imaging systems, which are available on current radiotherapy devices, are used, metal streak artifacts are avoided, and the soft-tissue contrast is restored, even for regions in which the kV data cannot contribute any information.

METHODS

Three image-reconstruction methods, including two filtered back-projection (FBP)-based analytic methods and one iterative method, for combining kV and MV projection data from the two on-board imaging systems of a radiotherapy device are presented in this work. The analytic reconstruction methods modify the MV data based on the information in the projection or image domains and then patch the data onto the kV projections for a FBP reconstruction. In the iterative reconstruction, the authors used dual-energy (DE) penalized weighted least-squares (PWLS) methods to simultaneously combine the kV/MV data and perform the reconstruction.

RESULTS

The authors compared kV/MV reconstructions to kV-only reconstructions using a dental phantom with fillings and a hip-implant numerical phantom. Simulation results indicated that dual-energy sinogram patch FBP and the modified dual-energy PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in the kV projections. The root-mean-square errors of soft-tissue patterns obtained using combined kV/MV data are 10-15 Hounsfield units smaller than those of the kV-only images, and the structural similarity index measure also indicates a 5%-10% improvement in the image quality. The added dose from the MV scan is much less than the dose from the kV scan if a high efficiency MV detector is assumed.

CONCLUSIONS

The authors have shown that it is possible to improve the image quality of kV CTs for patients with metal implants or fillings by completing the missing kV projection data with selectively acquired MV data that do not suffer from photon starvation. Numerical simulations demonstrated that dual-energy sinogram patch FBP and a modified kV/MV PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in kV projections. Combined kV/MV images may permit the improved delineation of structures of interest in CT images for patients with metal implants or fillings.

Monte Carlo simulation of a quantum noise limited Čerenkov detector based on air-spaced light guiding taper for megavoltage x-ray imaging

PURPOSE

Electronic Portal Imaging Devices (EPIDs) have been widely used in radiation therapy and are still needed on linear accelerators (Linacs) equipped with kilovoltage cone beam CT (kV-CBCT) or MRI systems. Our aim is to develop a new high quantum efficiency (QE) Čerenkov Portal Imaging Device (CPID) that is quantum noise limited at dose levels corresponding to a single Linac pulse.

METHODS

Recently a new concept of CPID for MV x-ray imaging in radiation therapy was introduced. It relies on Čerenkov effect for x-ray detection. The proposed design consisted of a matrix of optical fibers aligned with the incident x-rays and coupled to an active matrix flat panel imager (AMFPI) for image readout. A weakness of such design is that too few Čerenkov light photons reach the AMFPI for each incident x-ray and an AMFPI with an avalanche gain is required in order to overcome the readout noise for portal imaging application. In this work the authors propose to replace the optical fibers in the CPID with light guides without a cladding layer that are suspended in air. The air between the light guides takes on the role of the cladding layer found in a regular optical fiber. Since air has a significantly lower refractive index (∼ 1 versus 1.38 in a typical cladding layer), a much superior light collection efficiency is achieved.

RESULTS

A Monte Carlo simulation of the new design has been conducted to investigate its feasibility. Detector quantities such as quantum efficiency (QE), spatial resolution (MTF), and frequency dependent detective quantum efficiency (DQE) have been evaluated. The detector signal and the quantum noise have been compared to the readout noise.

CONCLUSIONS

Our studies show that the modified new CPID has a QE and DQE more than an order of magnitude greater than that of current clinical systems and yet a spatial resolution similar to that of current low-QE flat-panel based EPIDs. Furthermore it was demonstrated that the new CPID does not require an avalanche gain in the AMFPI and is quantum noise limited at dose levels corresponding to a single Linac pulse.

Rapid Monte Carlo simulation of detector DQE(f)

PURPOSE

Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10(7) - 10(9) detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS.

METHODS

The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels.

RESULTS

On the order of only 10-100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10(7) reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10(6) CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU.

CONCLUSIONS

The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.

ScintSim1: A new Monte Carlo simulation code for transport of optical photons in 2D arrays of scintillation detectors

Two-dimensional (2D) arrays of thick segmented scintillators are of interest as X-ray detectors for both 2D and 3D image-guided radiotherapy (IGRT). Their detection process involves ionizing radiation energy deposition followed by production and transport of optical photons. Only a very limited number of optical Monte Carlo simulation models exist, which has limited the number of modeling studies that have considered both stages of the detection process. We present ScintSim1, an in-house optical Monte Carlo simulation code for 2D arrays of scintillation crystals, developed in the MATLAB programming environment. The code was rewritten and revised based on an existing program for single-element detectors, with the additional capability to model 2D arrays of elements with configurable dimensions, material, etc., The code generates and follows each optical photon history through the detector element (and, in case of cross-talk, the surrounding ones) until it reaches a configurable receptor, or is attenuated. The new model was verified by testing against relevant theoretically known behaviors or quantities and the results of a validated single-element model. For both sets of comparisons, the discrepancies in the calculated quantities were all <1%. The results validate the accuracy of the new code, which is a useful tool in scintillation detector optimization.

Optimization of the performance of segmented scintillators for radiotherapy imaging through novel binning techniques

Thick, segmented crystalline scintillators have shown increasing promise as replacement x-ray converters for the phosphor screens currently used in active matrix flat-panel imagers (AMFPIs) in radiotherapy, by virtue of providing over an order of magnitude improvement in the detective quantum efficiency (DQE). However, element-to-element misalignment in current segmented scintillator prototypes creates a challenge for optimal registration with underlying AMFPI arrays, resulting in degradation of spatial resolution. To overcome this challenge, a methodology involving the use of a relatively high resolution AMFPI array in combination with novel binning techniques is presented. The array, which has a pixel pitch of 0.127 mm, was coupled to prototype segmented scintillators based on BGO, LYSO and CsI:Tl materials, each having a nominal element-to-element pitch of 1.016 mm and thickness of ∼ 1 cm. The AMFPI systems incorporating these prototypes were characterized at a radiotherapy energy of 6 MV in terms of modulation transfer function, noise power spectrum, DQE, and reconstructed images of a resolution phantom acquired using a cone-beam CT geometry. For each prototype, the application of 8 × 8 pixel binning to achieve a sampling pitch of 1.016 mm was optimized through use of an alignment metric which minimized misregistration and thereby improved spatial resolution. In addition, the application of alternative binning techniques that exclude the collection of signal near septal walls resulted in further significant improvement in spatial resolution for the BGO and LYSO prototypes, though not for the CsI:Tl prototype due to the large amount of optical cross-talk resulting from significant light spread between scintillator elements in that device. The efficacy of these techniques for improving spatial resolution appears to be enhanced for scintillator materials that exhibit mechanical hardness, high density and high refractive index, such as BGO. Moreover, materials that exhibit these properties as well as offer significantly higher light output than BGO, such as CdWO4, should provide the additional benefit of preserving DQE performance.

Countering beam divergence effects with focused segmented scintillators for high DQE megavoltage active matrix imagers

The imaging performance of active matrix flat-panel imagers designed for megavoltage imaging (MV AMFPIs) is severely constrained by relatively low x-ray detection efficiency, which leads to a detective quantum efficiency (DQE) of only ∼1%. Previous theoretical and empirical studies by our group have demonstrated the potential for addressing this constraint through the utilization of thick, two-dimensional, segmented scintillators with optically isolated crystals. However, this strategy is constrained by the degradation of high-frequency DQE resulting from spatial resolution loss at locations away from the central beam axis due to oblique incidence of radiation. To address this challenge, segmented scintillators constructed so that the crystals are individually focused toward the radiation source are proposed and theoretically investigated. The study was performed using Monte Carlo simulations of radiation transport to examine the modulation transfer function and DQE of focused segmented scintillators with thicknesses ranging from 5 to 60 mm. The results demonstrate that, independent of scintillator thickness, the introduction of focusing largely restores spatial resolution and DQE performance otherwise lost in thick, unfocused segmented scintillators. For the case of a 60 mm thick BGO scintillator and at a location 20 cm off the central beam axis, use of focusing improves DQE by up to a factor of ∼130 at non-zero spatial frequencies. The results also indicate relatively robust tolerance of such scintillators to positional displacements, of up to 10 cm in the source-to-detector direction and 2 cm in the lateral direction, from their optimal focusing position, which could potentially enhance practical clinical use of focused segmented scintillators in MV AMFPIs.

Does kV-MV dual-energy computed tomography have an advantage in determining proton stopping power ratios in patients?

Conventional kilovoltage (kV) x-ray-based dual-energy CT (DECT) imaging using two different x-ray energy spectra is sensitive to image noise and beam hardening effects. The purpose of this study was to evaluate the theoretical advantage of the DECT method for determining proton stopping power ratios (SPRs) using a combination of kV and megavoltage (MV) x-ray energies. We investigated three representative x-ray energy pairs: 100 and 140 kVp comprised the kV-kV pair, 100 kVp and 1 MV comprised the kV-MV pair, and two 1 MV x-ray beams-one with and one without external filtration-comprised the MV-MV pair. The SPRs of 34 human tissues were determined using the DECT method with these three x-ray energy pairs. Small perturbations were introduced into the CT numbers and x-ray spectra used for the DECT calculation to simulate the effects of random noise and beam hardening. An error propagation analysis was performed on the DECT calculation algorithm to investigate the propagation of CT number uncertainty to final SPR estimation and to suggest the best x-ray energy combination. We found that the DECT method using each of the three beam pairs achieved similar accuracy in determining the SPRs of human tissues in ideal conditions. However, when CT number uncertainties and artifacts such as imaging noise and beam hardening effects were considered, the kV-MV DECT improved the accuracy of SPR estimation substantially over the kV-kV or MV-MV DECT methods. Furthermore, our error propagation analysis showed that the combination of 100 kVp and 1 MV beams was close to the optimal selection when using the DECT method to determine SPRs. Overall, the kV-MV combination makes the DECT method more robust in resolving the effective atomic numbers for biological tissues than the traditional kV-kV DECT method.

Low-dose megavoltage cone-beam CT imaging using thick, segmented scintillators

Megavoltage, cone-beam computed tomography (MV CBCT) employing an electronic portal imaging device (EPID) is a highly promising technique for providing soft-tissue visualization in image-guided radiotherapy. However, current EPIDs based on active matrix flat-panel imagers (AMFPIs), which are regarded as the gold standard for portal imaging and referred to as conventional MV AMFPIs, require high radiation doses to achieve this goal due to poor x-ray detection efficiency (∼2% at 6 MV). To overcome this limitation, the incorporation of thick, segmented, crystalline scintillators, as a replacement for the phosphor screens used in these AMFPIs, has been shown to significantly improve the detective quantum efficiency (DQE) performance, leading to improved image quality for projection imaging at low dose. Toward the realization of practical AMFPIs capable of low dose, soft-tissue visualization using MV CBCT imaging, two prototype AMFPIs incorporating segmented scintillators with ∼11 mm thick CsI:Tl and Bi(4)Ge(3)O(12) (BGO) crystals were evaluated. Each scintillator consists of 120 × 60 crystalline elements separated by reflective septal walls, with an element-to-element pitch of 1.016 mm. The prototypes were evaluated using a bench-top CBCT system, allowing the acquisition of 180 projection, 360° tomographic scans with a 6 MV radiotherapy photon beam. Reconstructed images of a spatial resolution phantom, as well as of a water-equivalent phantom, embedded with tissue equivalent objects having electron densities (relative to water) varying from ∼0.28 to ∼1.70, were obtained down to one beam pulse per projection image, corresponding to a scan dose of ∼4 cGy--a dose similar to that required for a single portal image obtained from a conventional MV AMFPI. By virtue of their significantly improved DQE, the prototypes provided low contrast visualization, allowing clear delineation of an object with an electron density difference of ∼2.76%. Results of contrast, noise and contrast-to-noise ratio are presented as a function of dose and compared to those from a conventional MV AMFPI.

Performance evaluation of polycrystalline HgI2 photoconductors for radiation therapy imaging

PURPOSE

Electronic portal imaging devices based on megavoltage (MV), active matrix, flat-panel imagers (AMFPIs) are presently regarded as the gold standard in portal imaging for external beam radiation therapy. These devices, employing indirect detection of incident radiation by means of a metal plate plus phosphor screen combination, offer a quantum efficiency of only approximately 2% at 6 MV, leading to a detective quantum efficiency (DQE) of only approximately 1%. In order to significantly improve the DQE performance of MV AMFPIs, a strategy based on the development of direct detection imagers incorporating thick films of polycrystalline mercuric iodide (HgI2) photoconductor was undertaken and is reported.

METHODS

Two MV AMFPI prototypes, one incorporating an approximately 300 microm thick HgI2 layer created through physical vapor deposition (PVD) and a second incorporating an approximately 460 microm thick HgI2 layer created through screen-printing of particle-in-binder (PIB) material, were quantitatively evaluated using a 6 MV photon beam. The reported measurements include empirical determination of x-ray sensitivity, lag, modulation transfer function (MTF), noise power spectrum, and DQE.

RESULTS

For both prototypes, MTF and DQE results were found to be consistent with theoretical expectations and the MTFs were also found to be higher than that measured from a conventional MV AMFPI. In addition, the DQE results exhibit input-quantum-limited behavior, even at extremely low doses. Compared to PVD, the PIB prototype exhibits much lower dark current, slightly higher lag, and similar DQE. Finally, the challenges associated with this approach, as well as strategies for achieving considerably higher DQE through thicker HgI2 layers, are discussed.

CONCLUSIONS

The DQE of each of the prototypes is found to be comparable to that of conventional MV AMFPIs, commensurate with the modest photoconductor thicknesses of these early samples. It is anticipated that thicker layers of HgI2 based on PIB deposition can provide higher DQE while maintaining good material properties.

Monte Carlo investigations of the effect of beam divergence on thick, segmented crystalline scintillators for radiotherapy imaging

The use of thick, segmented scintillators in electronic portal imagers offers the potential for significant improvement in x-ray detection efficiency compared to conventional phosphor screens. Such improvement substantially increases the detective quantum efficiency (DQE), leading to the possibility of achieving soft-tissue visualization at clinically practical (i.e. low) doses using megavoltage (MV) cone-beam computed tomography. While these DQE increases are greatest at zero spatial frequency, they are diminished at higher frequencies as a result of degradation of spatial resolution due to lateral spreading of secondary radiation within the scintillator--an effect that is more pronounced for thicker scintillators. The extent of this spreading is even more accentuated for radiation impinging the scintillator at oblique angles of incidence due to beam divergence. In this paper, Monte Carlo simulations of radiation transport, performed to investigate and quantify the effects of beam divergence on the imaging performance of MV imagers based on two promising scintillators (BGO and CsI:Tl), are reported. In these studies, 10-40 mm thick scintillators, incorporating low-density polymer, or high-density tungsten septal walls, were examined for incident angles corresponding to that encountered at locations up to approximately 15 cm from the central beam axis (for an imager located 130 cm from a radiotherapy x-ray source). The simulations demonstrate progressively more severe spatial resolution degradation (quantified in terms of the effect on the modulation transfer function) as a function of increasing angle of incidence (as well as of the scintillator thickness). Since the noise power behavior was found to be largely independent of the incident angle, the dependence of the DQE on the incident angle is therefore primarily determined by the spatial resolution. The observed DQE degradation suggests that 10 mm thick scintillators are not strongly affected by beam divergence for detector areas up to approximately 30x30 cm2. For thicker scintillators, the area that is relatively unaffected is significantly reduced, requiring a focused scintillator geometry in order to preserve spatial resolution, and thus DQE.