Low‐dose megavoltage cone‐beam computed tomography using a novel multi‐layer imager (MLI)

MV-based relative electron density estimation (iMREDe) for MR-LINAC dose calculation

BACKGROUND

MR-LINAC systems have been increasingly utilized for real-time imaging in adaptive treatments worldwide. Challenges in MR representation of air cavities and subsequent estimation of electron density maps impede planning efficiency and may lead to dose calculation uncertainties.

PURPOSE

To demonstrate the generation of accurate electron density maps using the primary MV beam with a flat-panel imager.

METHODS

The ViewRay MRIdian MR-LINAC system was modeled digitally for Monte Carlo simulations. Iron shimming, the magnetic field, and the proposed flat panel detector were included in the model. The effect of the magnetic field on the detector response was investigated. Acquisition of projections over 360 degrees was simulated for digital phantoms of the Catphan 505 phantom and a patient treated for Head and Neck cancer. Shim patterns on the projections were removed and detector noise linearity was assessed. Electron density maps were generated for the digital patient phantom using the flat-panel detector and compared with actual treatment planning CT generated electron density maps of the same patient.

RESULTS

The effect of the magnetic field on the detector point-spread function (PSF) was found to be substantial for field strengths above 50 mT. Shims correction in the projection images using air normalization and in-painting effectively removed reconstruction artifacts without affecting noise linearity. The relative difference between reconstructed electron density maps from the proposed method and electron density maps generated from the treatment planning CT was 11% on average along all slices included in the iMREDe reconstruction.

CONCLUSIONS

The proposed iMREDe technique demonstrated the feasibility of generating accurate electron densities for the ViewRay MRIdian MR-LINAC system with a flat-panel imager and the primary MV beam. This work is a step towards reducing the time and effort required for adaptive radiotherapy in the current ViewRay MR-LINAC systems.

Monte Carlo model of a prototype flat-panel detector for multi-energy applications in radiotherapy

BACKGROUND

The incorporation of multi-energy capabilities into radiotherapy flat-panel detectors offers advantages including enhanced soft tissue visualization by reduction of signal from overlapping anatomy such as bone in 2D image projections; creation of virtual monoenergetic images for 3D contrast enhancement, metal artefact reduction and direct acquisition of relative electron density. A novel dual-layer on-board imager offering dual energy processing capabilities is being designed. As opposed to other dual-energy implementation techniques which require separate acquisition with two different x-ray spectra, the dual-layer detector design enables simultaneous acquisition of high and low energy images with a single exposure. A computational framework is required to optimize the design parameters and evaluate detector performance for specific clinical applications.

PURPOSE

In this study, we report on the development of a Monte Carlo (MC) model of the imager including model validation.

METHODS

The stack-up of the dual-layer imager (DLI) was implemented in GEANT4 Application for Tomographic Emission (GATE). The DLI model has an active area of 43×43 cm2 , with top and bottom Cesium Iodide (CsI) scintillators of 600 and 800 μm thickness, respectively. Measurement of spatial resolution and imaging of dedicated multi-material dual-energy (DE) phantoms were used to validate the model. The modulation transfer function (MTF) of the detector was calculated for a 120 kVp x-ray spectrum using a 0.5 mm thick tantalum edge rotated by 2.5o . For imaging validation, the DE phantom was imaged using a 140 kVp x-ray spectrum. For both validation simulations, corresponding measurements were done using an initial prototype of the imager. Agreement between simulations and measurement was assessed using normalized root mean square error (NRMSE) and 1D profile difference for the MTF and phantom images respectively. Further comparison between measurement and simulation was made using virtual monoenergetic images (VMIs) generated from basis material images derived using precomputed look-up tables.

RESULTS

The MTF of the bottom layer of the dual-layer model shows values decreasing more quickly with spatial frequency, compared to the top layer, due to the thicker bottom scintillator thickness and scatter from the top layer. A comparison with measurement shows NRMSE of 0.013 and 0.015 as well as identical MTF50 of 0.8 mm1 and 1.0 mm1 for the top and bottom layer respectively. For the DE imaging of the DE-phantom, although a maximum deviation of 3.3% is observed for the 10 mm aluminum and Teflon inserts at the top layer, the agreement for all other inserts is less than 2.2% of the measured value at both layers. Material decomposition of simulated scatter-free DE images gives an average accuracy in PMMA and aluminum composition of 4.9% and 10.3% for 11-30 mm PMMA and 1-10 mm aluminum objects respectively. A comparison of decomposed values using scatter containing measured and simulated DE images shows good agreement within statistical uncertainty.

CONCLUSION

Validation using both MTF and phantom imaging shows good agreement between simulation and measurements. With the present configuration of the digital prototype, the model can generate material decomposed images and virtual monoenergetic images.

Impact of a novel multilayer imager on metal artifacts in MV-CBCT

Objective. Megavoltage cone-beam computed tomography (MV-CBCT) imaging offers several advantages including reduced metal artifacts and accurate electron density mapping for adaptive or emergent situations. However, MV-CBCT imaging is limited by the poor efficiency of current detectors. Here we examine a new MV imager and compare CBCT reconstructions under clinically relevant scenarios.Approach. A multilayer imager (MLI), consisting of four vertically stacked standard flat-panel imagers, was mounted to a clinical linear accelerator. A custom anthropomorphic pelvis phantom with replaceable femoral heads was imaged using MV-CBCT and kilovoltage CBCT (kV-CBCT). Bone, aluminum, and titanium were used as femoral head inserts. 8 MU 2.5 MV scans were acquired for all four layers and (as reference) the top layer. Prostate and bladder were contoured on a reference CT and transferred to the other scans after rigid registration, from which the structural similarity index measure (SSIM) was calculated. Prostate and bladder were also contoured on CBCT scans without guidance, and Dice coefficients were compared to CT contours.Main results. kV-CBCT demonstrated the highest SSIMs with bone inserts (prostate: 0.86, bladder: 0.94) and lowest with titanium inserts (0.32, 0.37). Four-layer MV-CBCT SSIMs were preserved with bone (0.75, 0.80) as compared to titanium (0.67, 0.74), outperforming kV-CBCT when metal is present. One-layer MV-CBCT consistently underperformed four-layer results across all phantom configurations. Unilateral titanium inserts and bilateral aluminum insert results fell between the bone and bilateral titanium results. Dice coefficients trended similarly, with four-layer MV-CBCT reducing metal artifact impact relative to KV-CBCT to provide better soft-tissue identification.Significance. MV-CBCT with a four-layer MLI showed improvement over single-layer MV scans, approaching kV-CBCT quality for soft-tissue contrast. In the presence of artifact-producing metal implants, four-layer MV-CBCT scans outperformed kV-CBCT by eliminating artifacts and single-layer MV-CBCT by reducing noise. MV-CBCT with a novel multi-layer imager may be a valuable alternative to kV-CBCT, particularly in the presence of metal.

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.

Improving image quality and reducing dose with 2.5 MV diamond target volume-of-interest cone beam CT imaging

PURPOSE

Over the past decades, continuous efforts have been made to improve megavoltage (MV) image quality versus dose characteristics, including the implementation of low atomic number (Z) targets in MV beamlines and the development of more efficient detectors. Recently, a diamond target beam within a commercial radiotherapy treatment platform demonstrated improved planar contrast-to-noise-ratio (CNR) per unit dose using a novel 2.5 MV sintered diamond target beam, which enabled image acquisition on the order of mGy. The present work assesses cone beam CT (CBCT) image quality characteristics for the novel 2.5 MV diamond target beam and the effects of volume-of-interest (VOI) collimation on the image quality and imaging dose distribution.

METHODS

A sintered diamond target was incorporated into the target arm of the linear accelerator, replacing the 2.5 MV commercial copper imaging target. CBCT image quality was evaluated against the commercial imaging beam with regard to spatial resolution and CNR versus dose. In addition to full-field acquisitions, we investigated VOI techniques that collimate the imaging beam to preselected anatomy, to determine potential image quality improvements and dose sparing capacity. Using an anthropomorphic phantom, VOI regions were defined to encompass the maxillary and ethmoid sinuses and ranged in dimension from 3 cm to 4.85 cm equivalent radius. The MLC was fit to each VOI structure throughout a full CBCT arc and the corresponding MLC sequences were produced as XML scripts for acquisition. Calibrated radiochromic film was used in phantom to measure cumulative axial dose distributions during each CBCT acquisition.

RESULTS

In full-field CBCT, the 2.5 MV diamond target beam demonstrated improved CNR versus dose compared to the commercial imaging beam, by factors of up to 1.7. The calculated modulation transfer function (MTF) displayed an increase of nearly 30% in f₅₀ for the 2.5 MV diamond target beam compared to the commercial beam. Using VOI techniques, CNR increased monotonically as a function of equivalent radius at the bone-tissue interface. At the bone-sinus interface, the CNR for the full-field case was slightly decreased compared to the largest VOI case. Imaging dose in the anteroposterior direction increased with increasing VOI equivalent radius.

CONCLUSION

The novel 2.5 MV sintered diamond target beam presents a simple modification to the commercial imaging beam which provides improved image quality in full-field CBCT and the potential for simultaneous dose sparing and CNR improvement at high-contrast interfaces using VOI acquisition techniques.

Abbreviated on-treatment CBCT using roughness penalized mono-energization of kV-MV data and a multi-layer MV imager

Simultaneous acquisition of cone beam CT (CBCT) projections using both the kV and MV imagers of an image guided radiotherapy system reduces set-up scan times-a benefit to lung cancer radiation oncology patients-but increases noise in the 3D reconstruction. In this article, we present a kV-MV scan time reduction technique that uses two noise-reducing measures to achieve superior performance. The first is a high-DQE multi-layer MV imager prototype. The second is a beam hardening correction algorithm which combines poly-energetic modeling with edge-preserving, regularized smoothing of the projections. Performance was tested in real acquisitions of the Catphan 604 and a thorax phantom. Percent noise was quantified from voxel values in a soft tissue volume of interest (VOI) while edge blur was quantified from a VOI straddling a boundary between air and soft material. Comparisons in noise/resolution performance trade-off were made between our proposed approach, a dose-equivalent kV-only scan, and a kV-MV reconstruction technique previously published by Yinet al(2005Med. Phys.329). The proposed technique demonstrated lower noise as a function of spatial resolution than the baseline kV-MV method, notably a 50% noise reduction at typical edge blur levels. Our proposed method also exhibited fainter non-uniformity artifacts and in some cases superior contrast. Overall, we find that the combination of a multi-layer MV imager, acquiring at a LINAC source energy of 2.5 MV, and a denoised beam hardening correction algorithm enables noise, resolution, and dose performance comparable to standard kV-imager only set-up CBCT, but with nearly half the gantry rotation time.

Clinical translation of a new flat-panel detector for beam's-eye-view imaging

Electronic portal imaging devices (EPIDs) lend themselves to beams-eye view clinical applications, such as tumor tracking, but are limited by low contrast and detective quantum efficiency (DQE). We characterize a novel EPID prototype consisting of multiple layers and investigate its suitability for use under clinical conditions. A prototype multi-layer imager (MLI) was constructed utilizing four conventional EPID layers, each consisting of a copper plate, a Gd2O2S:Tb phosphor scintillator, and an amorphous silicon flat panel array detector. We measured the detector's response to a 6 MV photon beam with regards to modulation transfer function, noise power spectrum, DQE, contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), and the linearity of the detector's response to dose. Additionally, we compared MLI performance to the single top layer of the MLI and the standard Varian AS-1200 detector. Pre-clinical imaging was done on an anthropomorphic phantom, and the detector's CNR, SNR and spatial resolution were assessed in a clinical environment. Images obtained from spine and liver patient treatment deliveries were analyzed to verify CNR and SNR improvements. The MLI has a DQE(0) of 9.7%, about 5.7 times the reference AS-1200 detector. Improved noise performance largely drives the increase. CNR and SNR of clinical images improved three-fold compared to reference. A novel MLI was characterized and prepared for clinical translation. The MLI substantially improved DQE and CNR performance while maintaining the same resolution. Pre-clinical tests on an anthropomorphic phantom demonstrated improved performance as predicted theoretically. Preliminary patient data were analyzed, confirming improved CNR and SNR. Clinical applications are anticipated to include more accurate soft tissue tracking.