Single-pass metal artifact reduction using a dual-layer flat panel detector.

PURPOSE: Metal artifact remains a challenge in cone-beam CT images. Many image domain-based segmentation methods have been proposed for metal artifact reduction (MAR), which require two-pass reconstruction. Such methods first segment metal from a first-pass reconstruction and then forward-project the metal mask to identify them in projections. These methods work well in general but are limited when the metal is outside the scan field-of-view (FOV) or when the metal is moving during the scan. In the former, even reconstructing with a larger FOV does not guarantee a good estimate of metal location in the projections; and in the latter, the metal location in each projection is difficult to identify due to motion. Single-pass methods that detect metal in single-energy projections have also been developed, but often have imperfect metal detection that leads to residual artifacts. In this work, we develop a MAR method using a dual-layer (DL) flat panel detector, which improves performance for single-pass reconstruction. METHODS: In this work, we directly detect metal objects in projections using dual-energy (DE) imaging that generates material-specific images (e.g., soft tissue and bone), where the metal stands out in bone images when nonuniform soft tissue background is removed. Metal is detected via simple thresholding, and entropy filtration is further applied to remove false-positive detections. A DL detector provides DE images with superior temporal and spatial registration and was used to perform the task. Scatter correction was first performed on DE raw projections to improve the accuracy of material decomposition. One phantom mimicking a liver biopsy setup and a cadaver head were used to evaluate the metal reduction performance of the proposed method and compared with that of a standard two-pass reconstruction, a previously published sinogram-based method using a Markov random field (MRF) model, and a single-pass projection-domain method using single-energy imaging. The phantom has a liver steering setup placed in a hollow chest phantom, with embedded metal and a biopsy needle crossing the phantom boundary. The cadaver head has dental fillings and a metal tag attached to its surface. The identified metal regions in each projection were corrected by interpolation using surrounding pixels, and the images were reconstructed using filtered backprojection. RESULTS: Our current approach removes metal from the projections, which is robust to FOV truncation during imaging acquisition. In case of FOV truncation, the method outperformed the two-pass reconstruction method. The proposed method using DE renders better accuracy in metal segmentation than the MRF method and single-energy method, which were prone to false-positive errors that cause additional streaks. For the liver steering phantom, the average spatial nonuniformity was reduced from 0.127 in uncorrected images to 0.086 using a standard two-pass reconstruction and to 0.077 using the proposed method. For the cadaver head, the average standard deviation within selected soft tissue regions ( σ s ) was reduced from 209.1 HU in uncorrected images to 69.1 HU using a standard two-pass reconstruction and to 46.8 HU using our proposed method. The proposed method reduced the processing time by 31% as compared with the two-pass method. CONCLUSIONS: We proposed a MAR method that directly detects metal in the projection domain using DE imaging, which is robust to truncation and superior to that of single-energy imaging. The method requires only a single-pass reconstruction that substantially reduces processing time compared with the standard two-pass metal reduction method.

High-resolution model-based material decomposition in dual-layer flat-panel CBCT.

PURPOSE: Spectral CT uses energy-dependent measurements that enable material discrimination in addition to reconstruction of structural information. Flat-panel detectors (FPDs) have been widely used in dedicated and interventional systems to deliver high spatial resolution, volumetric cone-beam CT (CBCT) in compact and OR-friendly designs. In this work, we derive a model-based method that facilitates high-resolution material decomposition in a spectral CBCT system equipped with a prototype dual-layer FPD. Through high-fidelity modeling of multilayer detector, we seek to avoid resolution loss that is present in more traditional processing and decomposition approaches. METHOD: A physical model for spectral measurements in dual-layer flat-panel CBCT is developed including layer-dependent differences in system geometry, spectral sensitivities, and detector blur (e.g., due to varied scintillator thicknesses). This forward model is integrated into a model-based material decomposition (MBMD) method based on minimization of a penalized weighted least-squared (PWLS) objective function. The noise and resolution performance of this approach was compared with traditional projection-domain decomposition (PDD) and image-domain decomposition (IDD) approaches as well as one-step MBMD with lower-fidelity models that use approximated geometry, projection interpolation, or an idealized system geometry without system blur model. Physical studies using high-resolution three-dimensional (3D)-printed water-iodine phantoms were conducted to demonstrate the high-resolution imaging performance of the compared decomposition methods in iodine basis images and synthetic monoenergetic images. RESULTS: Physical experiments demonstrate that the MBMD methods incorporating an accurate geometry model can yield higher spatial resolution iodine basis images and synthetic monoenergetic images than PDD and IDD results at the same noise level. MBMD with blur modeling can further improve the spatial-resolution compared with the decomposition results obtained with IDD, PDD, and MBMD methods with lower-fidelity models. Using the MBMD without or with blur model can increase the absolute modulation at 1.75 lp/mm by 10% and 22% compared with IDD at the same noise level. CONCLUSION: The proposed model-based material decomposition method for a dual-layer flat-panel CBCT system has demonstrated an ability to extend high-resolution performance through sophisticated detector modeling including the layer-dependent blur. The proposed work has the potential to not only facilitate high-resolution spectral CT in interventional and dedicated CBCT systems, but may also provide the opportunity to evaluate different flat-panel design trade-offs including multilayer FPDs with mismatched geometries, scintillator thicknesses, and spectral sensitivities.

High-Resolution Model-based Material Decomposition for Multi-Layer Flat-Panel Detectors.

In this work we compare a novel model-based material decomposition (MBMD) approach against a standard approach in high-resolution spectral CT using multi-layer flat-panel detectors. Physical experiments were conducted using a prototype dual-layer detector and a custom high-resolution iodine-enhanced line-pair phantom. Reconstructions were performed using three methods: traditional filtered back-projection (FBP) followed by image-domain decomposition, idealized MBMD with no blur modeling (iMBMD), and MBMD with system blur modeling (bMBMD). We find that both MBMD methods yielded higher resolution decompositions with lower noise than the FBP method, and that bMBMD further improves spatial resolution over iMBMD due to the additional blur modeling. These results demonstrate the advantages of MBMD in resolution performance and noise control over traditional methods for spectral CT. Model-based material decomposition hence has great potential in high-resolution spectral CT applications.

Model-based Material Decomposition with System Blur Modeling.

In this work, we present a novel model-based material decomposition (MBMD) approach for x-ray CT that includes system blur in the measurement model. Such processing has the potential to extend spatial resolution in material density estimates - particularly in systems where different spectral channels exhibit different spatial resolutions. We illustrate this new approach for a dual-layer detector x-ray CT and compare MBMD algorithms with and without blur in the reconstruction forward model. Both qualitative and quantitative comparisons of performance with and without blur modeling are reported. We find that blur modeling yields images with better recovery of high-resolution structures in an investigation of reconstructed line pairs as well as lower cross-talk bias between material bases that is ordinarily found due to mismatches in spatial resolution between spectral channels. The extended spatial resolution of the material decompositions has potential application in a range of high-resolution clinical tasks and spectral CT systems where spectral channels exhibit different spatial resolutions.

Characterization and potential applications of a dual-layer flat-panel detector.

PURPOSE: Dual-energy (DE) x-ray imaging has many clinical applications in radiography, fluoroscopy, and CT. This work characterizes a prototype dual-layer (DL) flat-panel detector (FPD) and investigates its DE imaging capabilities for applications in two-dimensional (2D) radiography/fluoroscopy and quantitative three-dimensional (3D) cone-beam CT. Unlike other DE methods like kV switching, a DL FPD obtains DE images from a single exposure, making it robust against patient and system motion. METHODS: The DL FPD consists of a top layer with a 200 µm-thick CsI scintillator coupled to an amorphous silicon (aSi) FPD of 150 µm pixel size and a bottom layer with a 550 µm thick CsI scintillator coupled to an identical aSi FPD. The two layers are separated by a 1-mm Cu filter to increase spectral separation. Images (43 × 43 cm2 active area) can be readout in 2 × 2 binning mode (300 µm pixels) at up to 15 frames per second. Detector performance was first characterized by measuring the MTF, NPS, and DQE for the top and bottom layers. For 2D applications, a qualitative study was conducted using an anthropomorphic thorax phantom containing a porcine heart with barium-filled coronary arteries (similar to iodine). Additionally, fluoroscopic lung tumor tracking was investigated by superimposing a moving tumor phantom on the thorax phantom. Tracking accuracies of single-energy (SE) and DE fluoroscopy were compared against the ground truth motion of the tumor. For 3D quantitative imaging, a phantom containing water, iodine, and calcium inserts was used to evaluate overall DE material decomposition capabilities. Virtual monoenergetic (VM) images ranging from 40 to 100 keV were generated, and the optimal VM image energy which achieved the highest image uniformity and maximum contrast-to-noise ratio (CNR) was determined. RESULTS: The spatial resolution of the top layer was substantially higher than that of the bottom layer (top layer 50% MTF = 2.2 mm-1 , bottom layer = 1.2 mm-1 ). A substantial increase in NNPS and reduction in DQE were observed for the bottom layer mainly due to photon loss within the top layer and Cu filter. For 2D radiographic and fluoroscopic applications, the DL FPD was capable of generating high-quality material-specific images separating soft tissue from bone and barium. For lung tumor tracking, DE fluoroscopy yielded more accurate results than SE fluoroscopy, with an average reduction in the root mean square error (RMSE) of over 10×. For the DE-CBCT studies, accurate basis material decompositions were obtained. The estimated material densities were 294.68  ±  17.41 and 92.14  ±  15.61 mg/ml for the 300 and 100 mg/ml calcium inserts, respectively, and 8.93  ±  1.45, 4.72  ±  1.44, and 2.11  ±  1.32 mg/ml for the 10, 5, and 2 mg/ml iodine inserts, respectively, with an average error of less than 5%. The optimal VM image energy was found to be 60 keV. CONCLUSIONS: We characterized a prototype DL FPD and demonstrated its ability to perform accurate single-exposure DE radiography/fluoroscopy and DE-CBCT. The merits of the DL detector approach include superior spatial and temporal registration between its constituent images, and less complicated acquisition sequences.

Projection-domain metal artifact correction using a dual layer detector.

Metal artifact remains a challenge in cone-beam CT images. Many two-pass metal artifact reduction methods have been proposed, which work fairly well, but are limited when the metal is outside the scan field-of-view (FOV) or when the metal is moving during the scan. In the former, even reconstructing with a larger FOV does not guarantee a good estimate of metal location in the projections; and in the latter, the metal location in each projection is difficult to identify due to motion. Furthermore, two-pass methods increase the total reconstruction time. In this study, a projection-based metal detection and correction method with a dual layer detector is investigated. The dual layer detector provides dual energy images with perfect temporal and spatial registration in each projection, which aid in the identification of metal. A simple phantom with metal wires (copper) and a needle (steel) is used to evaluate the projection-based metal artifact reduction method from a dual layer scan and compared with that of a single layer scan. Preliminary results showed enhanced ability to identify metal regions, leading to substantially reduced metal artifact in reconstructed images. In summary, an effective single-pass, projection-domain method using a dual layer detector has been demonstrated, and it is expected to be robust against truncation and motion.

Comparative Study of Dual Energy Cone-Beam CT using a Dual-Layer Detector and kVp Switching for Material Decomposition.

Cone-beam CT (CBCT) is widely used in diagnostic imaging and image-guided procedures, leading to an increasing need for advanced CBCT techniques, such as dual energy (DE) imaging. Previous studies have shown that DE-CBCT can perform quantitative material decomposition, including quantification of contrast agents, electron density, and virtual monoenergetic images. Currently, most CBCT systems perform DE imaging using a kVp switching technique. However, the disadvantages of this method are spatial and temporal misregistration as well as total scan time increase, leading to errors in the material decomposition. DE-CBCT with a dual layer flat panel detector potentially overcomes these limitations by acquiring the dual energy images simultaneously. In this work, we investigate the DE imaging performance of a prototype dual layer detector by evaluating its material decomposition capability and comparing its performance to that of the kVp switching method. Two sets of x-ray spectra were used for kVp switching: 80/120 kVp and 80/120 kVp + 1 mm Cu filtration. Our results show the dual layer detector outperforms kVp switching at 80/120 kVp with matched dose. The performance of kVp switching was better by adding 1 mm copper filtration to the high energy images (80/120 kVp + 1 mm Cu), though the dual layer detector still provided comparable performance for material decomposition tasks. Overall, both the dual layer detector and kVp switching methods provided quantitative material decomposition images in DE-CBCT, with the dual layer detector having additional potential advantages.

A fast, linear Boltzmann transport equation solver for computed tomography dose calculation (Acuros CTD).

PURPOSE: To improve dose reporting of CT scans, patient-specific organ doses are highly desired. However, estimating the dose distribution in a fast and accurate manner remains challenging, despite advances in Monte Carlo methods. In this work, we present an alternative method that deterministically solves the linear Boltzmann transport equation (LBTE), which governs the behavior of x-ray photon transport through an object. METHODS: Our deterministic solver for CT dose (Acuros CTD) is based on the same approach used to estimate scatter in projection images of a CT scan (Acuros CTS). A deterministic method is used to compute photon fluence within the object, which is then converted to deposited energy by multiplying by known, material-specific conversion factors. To benchmark Acuros CTD, we used the AAPM Task Group 195 test for CT dose, which models an axial, fan beam scan (10 mm thick beam) and calculates energy deposited in each organ of an anthropomorphic phantom. We also validated our own Monte Carlo implementation of Geant4 to use as a reference to compare Acuros against for other common geometries like an axial, cone beam scan (160 mm thick beam) and a helical scan (40 mm thick beam with table motion for a pitch of 1). RESULTS: For the fan beam scan, Acuros CTD accurately estimated organ dose, with a maximum error of 2.7% and RMSE of 1.4% when excluding organs with <0.1% of the total energy deposited. The cone beam and helical scans yielded similar levels of accuracy compared to Geant4. Increasing the number of source positions beyond 18 or decreasing the voxel size below 5 × 5 × 5 mm3 provided marginal improvement to the accuracy for the cone beam scan but came at the expense of increased run time. Across the different scan geometries, run time of Acuros CTD ranged from 8 to 23 s. CONCLUSIONS: In this digital phantom study, a deterministic LBTE solver was capable of fast and accurate organ dose estimates.

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

Feasibility of closed-MLC tracking using high sensitivity and multi-layer electronic portal imagers.

In radiation therapy, improvements in treatment conformality are often limited by movement of target tissue. To better treat the target, tumor tracking strategies involving beam's-eye-view (BEV) have been explored. However, localization surrogates like implanted fiducial markers may sometimes leave the field-of-view (FOV), as defined by the linear accelerator (LINAC) multi-leaf collimator (MLC). Radiation leakage through the MLC has been measured previously at approximately 1%-2%. High sensitivity prototype detectors imagers may improve the ability to visualize objects outside of the MLC FOV during treatment. The present study presents a proof-of-concept for tracking fiducial markers outside the MLC FOV by employing high sensitivity detectors using a high-efficiency, prototype scintillating glass called LKH-5 and also investigates the impact of multi-layer imager (MLI) architecture. It was found that by improving the detector efficiency, using either of these methods results in a reduction of dose required for fiducial marker visibility. Further, image correction by a rectangular median filter will improve fiducial marker representation in the MLC blocked images. Quantified by measuring the peak-to-sidelobe ratio (PSR) of the normalized cross correlation (NCC) between a template of the fiducial marker with the blocked MLC acquisition, visibility has been found at a threshold of roughly 5 for all configurations with a 3 × 3 cm2 ROI. For typical gadolinium oxysulfide (GOS) detectors in single and simulated 4-layer configurations, the minimum dose required for visualization was 20 and 10 MU, respectively. For LKH-5 detectors in single and simulated 4-layer configurations, this minimum dose was reduced to 4 and 2 MU, respectively. With a 6 MV flattening filter free (FFF) beam dose rate of 1400 MU min-1, the maximum detector frame rate while maintaining fiducial visibility is approximately 12 fps for a 4-layer LKH-5 configuration.

Dual edge apparatus and algorithm for measurement of x-ray beam spot parameters.

PURPOSE: The shape, size, and location of the x-ray beam spot (where the electron beam strikes the target) in a linac-based radiation therapy machine are of potential clinical significance. Established techniques to measure the beam spot parameters involve specialized hardware and typically assess size and shape of the beam spot or its position, but not both. A simple apparatus and algorithm for measuring all beam spot parameters simultaneously is proposed here. METHODS: The apparatus is composed of two partially transmitting edge plates mounted at different vertical positions. The mount for the apparatus slides into the accessory tray of the linac treatment head so that it rotates with the collimator, and it is imaged by the existing electronic portal imaging device (EPID) over multiple collimator angles. A software algorithm takes the acquired images and uses a parallel-beam CT reconstruction technique to compute beam spot size, shape, and position in one computation. In addition, the wobble of the collimator assembly can be estimated. The overall method was validated with both Monte Carlo simulation and with comparison to in-house spot camera measurements on a radiation therapy system. RESULTS: The algorithm correctly predicted the beam spot parameters used for the Monte Carlo simulation to better than 50 µm accuracy in all cases. Furthermore, results from the dual edge method matched spot camera results with 30 µm accuracy for beam spot size and shape, with 80 µm average accuracy for beam spot position, and better than 200 µm accuracy for collimator assembly wobble. CONCLUSIONS: We have developed a combination dual edge apparatus and image processing algorithm that, when used on a radiotherapy linac with an EPID, can accurately determine the size and shape of the electron beam spot, its position relative to collimator rotation axis, and the wobble of the collimator assembly.

A modified McKinnon-Bates (MKB) algorithm for improved 4D cone-beam computed tomography (CBCT) of the lung.

PURPOSE: Four-dimensional (4D) cone-beam computed tomography (CBCT) of the lung is an effective tool for motion management in radiotherapy but presents a challenge because of slow gantry rotation times. Sorting the individual projections by breathing phase and using an established technique such as Feldkamp-Davis-Kress (FDK) to generate corresponding phase-correlated (PC) three-dimensional (3D) images results in reconstructions (FDK-PC) that often contain severe streaking artifacts due to the sparse angular sampling distributions. These can be reduced by further slowing down the gantry at the expense of incurring unwanted increases in scan times and dose. A computationally efficient alternative is the McKinnon-Bates (MKB) reconstruction algorithm that has shown promise in reducing view aliasing-induced streaking but can produce ghosting artifacts that reduce contrast and impede the determination of motion trajectories. The purpose of this work was to identify and correct shortcomings in the MKB algorithm. METHODS: In the general MKB approach, a time-averaged 3D prior image is first reconstructed. The prior is then forward-projected at the same angles as the original projection data creating time-averaged reprojections. These reprojections are subsequently subtracted from the original (unblurred) projections to create motion-encoded difference projections. The difference projections are reconstructed into PC difference images that are added to the well-sampled 3D prior to create the higher quality 4D image. The cause of the ghosting in the traditional 4D MKB images was studied and traced to motion-induced streaking in the prior that, when reprojected, has the undesirable effect of re-encoding for motion in what should be a purely time-averaged reprojection. A new method, designated as the modified McKinnon-Bates (mMKB) algorithm, was developed based on destreaking the prior. This was coupled with a postprocessing 4D bilateral filter for noise suppression and edge preservation (mMKBbf ). The algorithms were tested with the 4D XCAT phantom using four simulated scan times (57, 60, 120, 180 s) and with two in vivo thorax studies (acquisition time of 60 and 90 s). Contrast-to-noise ratios (CNRs) of the target lesions and overall visual quality of the images were assessed. RESULTS: Prior destreaking (mMKB algorithm) reduced ghosting artifacts and increased CNRs for all cases, with the biggest impacts seen in the end inhale (EI) and end exhale (EE) phases of the respiratory cycle. For the XCAT phantom, mMKB lesion CNR was 44% higher than the MKB lesion CNR and was 81% higher than the FDK-PC lesion CNR (EI and EE phases). The bilateral filter provided a further average CNR improvement of 87% with the highest increases associated with longer scan times. Across all phases and scan times, the maximum mMKBbf -to-FDK-PC CNR improvement was over 300%. In vivo results agreed with XCAT results. Significantly less ghosting was observed throughout the mMKB images including near the lesions-of-interest and the diaphragm allowing for, in one case, visualization of a small tumor with nearly 30 mm of motion. The maximum FDK-PC-to-MKBbf CNR improvement for Patient 1's lesion was 261% and for Patient 2's lesion was 318%. CONCLUSIONS: The 4D mMKB algorithm yields good quality coronal and sagittal images in the thorax that may provide sufficient information for patient verification.

Physics considerations in MV-CBCT multi-layer imager design.

Megavoltage (MV) cone-beam computed tomography (CBCT) using an electronic portal imaging (EPID) offers advantageous features, including 3D mapping, treatment beam registration, high-z artifact suppression, and direct radiation dose calculation. Adoption has been slowed by image quality limitations and concerns about imaging dose. Developments in imager design, including pixelated scintillators, structured phosphors, inexpensive scintillation materials, and multi-layer imager (MLI) architecture have been explored to improve EPID image quality and reduce imaging dose. The present study employs a hybrid Monte Carlo and linear systems model to determine the effect of detector design elements, such as multi-layer architecture and scintillation materials. We follow metrics of image quality including modulation transfer function (MTF) and noise power spectrum (NPS) from projection images to 3D reconstructions to in-plane slices and apply a task based figure-of-merit, the ideal observer signal-to-noise ratio (d') to determine the effect of detector design on object detectability. Generally, detectability was limited by detector noise performance. Deploying an MLI imager with a single scintillation material for all layers yields improvement in noise performance and d' linear with the number of layers. In general, improving x-ray absorption using thicker scintillators results in improved DQE(0). However, if light yield is low, performance will be affected by electronic noise at relatively high doses, resulting in rapid image quality degradation. Maximizing image quality in a heterogenous MLI detector (i.e. multiple different scintillation materials) is most affected by limiting total noise. However, while a second-order effect, maximizing total spatial resolution of the MLI detector is a balance between the intensity contribution of each layer against its individual MTF. So, while a thinner scintillator may yield a maximal individual-layer MTF, its quantum efficiency will be relatively low in comparison to a thicker scintillator and thus, intensity contribution may be insufficient to noticeably improve the total detector MTF.

Multi-layer imager design for mega-voltage spectral imaging.

The architecture of multi-layer imagers (MLIs) can be exploited to provide megavoltage spectral imaging (MVSPI) for specific imaging tasks. In the current work, we investigated bone suppression and gold fiducial contrast enhancement as two clinical tasks which could be improved with spectral imaging. A method based on analytical calculations that enables rapid investigation of MLI component materials and thicknesses was developed and validated against Monte Carlo computations. The figure of merit for task-specific imaging performance was the contrast-to-noise ratio (CNR) of the gold fiducial when the CNR of bone was equal to zero after a weighted subtraction of the signals obtained from each MLI layer. Results demonstrated a sharp increase in the CNR of gold when the build-up component or scintillation materials and thicknesses were modified. The potential for low-cost, prompt implementation of specific modifications (e.g. composition of the build-up component) could accelerate clinical translation of MVSPI.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part I: Core algorithms and validation.

PURPOSE: To describe Acuros® CTS, a new software tool for rapidly and accurately estimating scatter in x-ray projection images by deterministically solving the linear Boltzmann transport equation (LBTE). METHODS: The LBTE describes the behavior of particles as they interact with an object across spatial, energy, and directional (propagation) domains. Acuros CTS deterministically solves the LBTE by modeling photon transport associated with an x-ray projection in three main steps: (a) Ray tracing photons from the x-ray source into the object where they experience their first scattering event and form scattering sources. (b) Propagating photons from their first scattering sources across the object in all directions to form second scattering sources, then repeating this process until all high-order scattering sources are computed using the source iteration method. (c) Ray-tracing photons from scattering sources within the object to the detector, accounting for the detector's energy and anti-scatter grid responses. To make this process computationally tractable, a combination of analytical and discrete methods is applied. The three domains are discretized using the Linear Discontinuous Finite Elements, Multigroup, and Discrete Ordinates methods, respectively, which confer the ability to maintain the accuracy of a continuous solution. Furthermore, through the implementation in CUDA, we sought to exploit the parallel computing capabilities of graphics processing units (GPUs) to achieve the speeds required for clinical utilization. Acuros CTS was validated against Geant4 Monte Carlo simulations using two digital phantoms: (a) a water phantom containing lung, air, and bone inserts (WLAB phantom) and (b) a pelvis phantom derived from a clinical CT dataset. For these studies, we modeled the TrueBeam® (Varian Medical Systems, Palo Alto, CA) kV imaging system with a source energy of 125 kVp. The imager comprised a 600 µm-thick Cesium Iodide (CsI) scintillator and a 10:1 one-dimensional anti-scatter grid. For the WLAB studies, the full-fan geometry without a bowtie filter was used (with and without the anti-scatter grid). For the pelvis phantom studies, a half-fan geometry with bowtie was used (with the anti-scatter grid). Scattered and primary photon fluences and energies deposited in the detector were recorded. RESULTS: The Acuros CTS and Monte Carlo results demonstrated excellent agreement. For the WLAB studies, the average percent difference between the Monte Carlo- and Acuros-generated scattered photon fluences at the face of the detector was -0.7%. After including the detector response, the average percent differences between the Monte Carlo- and Acuros-generated scatter fractions (SF) were -0.1% without the grid and 0.6% with the grid. For the digital pelvis simulation, the Monte Carlo- and Acuros-generated SFs agreed to within 0.1% on average, despite the scatter-to-primary ratios (SPRs) being as high as 5.5. The Acuros CTS computation time for each scatter image was ~1 s using a single GPU. CONCLUSIONS: Acuros CTS enables a fast and accurate calculation of scatter images by deterministically solving the LBTE thus offering a computationally attractive alternative to Monte Carlo methods. Part II describes the application of Acuros CTS to scatter correction of CBCT scans on the TrueBeam system.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part II: System modeling, scatter correction, and optimization.

PURPOSE: To correct for scatter in kV cone-beam CT (CBCT) projection data on a clinical system using a new tool, Acuros® CTS, that estimates scatter images rapidly and accurately by deterministically solving the linear Boltzmann transport equation. METHODS: Phantom and patient CBCT scans were acquired on TrueBeam® radiotherapy machines. A first-pass reconstruction was used to create water and bone density maps of the imaged object, which was updated to include a more accurate representation of the patient couch. The imaging system model accounted for the TrueBeam x-ray source (polychromatic spectrum, beam filtration, bowtie filter, and collimation hardware) and x-ray detection system (antiscatter grid, flat-panel imager). Acuros CTS then used the system and object models to estimate the scatter component of each projection image, which was subtracted from the measured projections. The corrected projections were then reconstructed to produce the final result. We examined the tradeoff between run time and accuracy using a Pareto optimization of key parameters, including the voxel size of the down-sampled object model, the number of pixels in the down-sampled detector, and the number of scatter images (angular down-sampling). All computations and reconstructions were performed on a research workstation containing two graphics processing units (GPUs). In addition, we established a method for selecting a subset of projections for which scatter images were calculated. The projections were selected to minimize interpolation errors in the remaining projections. Image quality improvement was assessed by measuring the accuracy of the reconstructed phantom and patient images. RESULTS: The Pareto optimization yielded a set of parameters with an average run time of 26 seconds for scatter correction while maintaining high accuracy of scatter estimation. This was achieved in part by means of optimizing the projection angles that were processed, thus favoring the use of more angles in the lateral (i.e., horizontal) direction and fewer angles in the AP direction. In a 40 cm solid water phantom reconstruction, nonuniformities were decreased from 217 HU without scatter correction to 51 HU with conventional (kernel-based) scatter correction to 17 HU with Acuros CTS-based scatter correction. In clinical pelvis scans, nonuniformities in the bladder were reduced from 85 HU with conventional scatter correction to 14 HU with Acuros CTS. CONCLUSIONS: Acuros CTS is a promising new tool for fast and accurate scatter correction for CBCT imaging. By carefully modeling the imaging chain and optimizing several parameters, we achieved high correction accuracies with computation times compatible with the clinical workflow. The improvement in image quality enables better soft-tissue visualization and potentially enables applications such as adaptive radiotherapy.

Leveraging multi-layer imager detector design to improve low-dose performance for megavoltage cone-beam computed tomography.

While megavoltage cone-beam computed tomography (CBCT) using an electronic portal imaging device (EPID) provides many advantages over kilovoltage (kV) CBCT, clinical adoption is limited by its high doses. Multi-layer imager (MLI) EPIDs increase DQE(0) while maintaining high resolution. However, even well-designed, high-performance MLIs suffer from increased electronic noise from each readout, degrading low-dose image quality. To improve low-dose performance, shift-and-bin addition (ShiBA) imaging is proposed, leveraging the unique architecture of the MLI. ShiBA combines hardware readout-binning and super-resolution concepts, reducing electronic noise while maintaining native image sampling. The imaging performance of full-resolution (FR); standard, aligned binned (BIN); and ShiBA images in terms of noise power spectrum (NPS), electronic NPS, modulation transfer function (MTF), and the ideal observer signal-to-noise ratio (SNR)-the detectability index (d')-are compared. The FR 4-layer readout of the prototype MLI exhibits an electronic NPS magnitude 6-times higher than a state-of-the-art single layer (SLI) EPID. Although the MLI is built on the same readout platform as the SLI, with each layer exhibiting equivalent electronic noise, the multi-stage readout of the MLI results in electronic noise 50% higher than simple summation. Electronic noise is mitigated in both BIN and ShiBA imaging, reducing its total by ~12 times. ShiBA further reduces the NPS, effectively upsampling the image, resulting in a multiplication by a sinc2 function. Normalized NPS show that neither ShiBA nor BIN otherwise affects image noise. The LSF shows that ShiBA removes the pixilation artifact of BIN images and mitigates the effect of detector shift, but does not quantifiably improve the MTF. ShiBA provides a pre-sampled representation of the images, mitigating phase dependence. Hardware binning strategies lower the quantum noise floor, with 2 × 2 implementation reducing the dose at which DQE(0) degrades by 10% from 0.01 MU to 0.004 MU, representing 20% improvement in d'.

Spectral imaging using clinical megavoltage beams and a novel multi-layer imager.

We assess the feasibility of clinical megavoltage (MV) spectral imaging for material and bone separation with a novel multi-layer imager (MLI) prototype. The MLI provides higher detective quantum efficiency and lower noise than conventional electronic portal imagers. Simulated experiments were performed using a validated Monte Carlo model of the MLI to estimate energy absorption and energy separation between the MLI components. Material separation was evaluated experimentally using solid water and aluminum (Al), copper (Cu) and gold (Au) for 2.5 MV, 6 MV and 6 MV flattening filter free (FFF) clinical photon beams. An anthropomorphic phantom with implanted gold fiducials was utilized to further demonstrate bone/gold separation. Weighted subtraction imaging was employed for material and bone separation. The weighting factor (w) was iteratively estimated, with the optimal w value determined by minimization of the relative signal difference ([Formula: see text]) and signal-difference-to-noise ratio (SDNR) between material (or bone) and the background. Energy separation between layers of the MLI was mainly the result of beam hardening between components with an average energy separation between 34 and 47 keV depending on the x-ray beam energy. The minimum average energy of the detected spectrum in the phosphor layer was 123 keV in the top layer of the MLI with the 2.5 MV beam. The w values that minimized [Formula: see text] and SDNR for Al, Cu and Au were 0.89, 0.76 and 0.64 for 2.5 MV; for 6 MV FFF, w was 0.98, 0.93 and 0.77 respectively. Bone suppression in the anthropomorphic phantom resulted in improved visibility of the gold fiducials with the 2.5 MV beam. Optimization of the MLI design is required to achieve optimal separation at clinical MV beam energies.

A novel method for quantification of beam's-eye-view tumor tracking performance.

PURPOSE: In-treatment imaging using an electronic portal imaging device (EPID) can be used to confirm patient and tumor positioning. Real-time tumor tracking performance using current digital megavolt (MV) imagers is hindered by poor image quality. Novel EPID designs may help to improve quantum noise response, while also preserving the high spatial resolution of the current clinical detector. Recently investigated EPID design improvements include but are not limited to multi-layer imager (MLI) architecture, thick crystalline and amorphous scintillators, and phosphor pixilation and focusing. The goal of the present study was to provide a method of quantitating improvement in tracking performance as well as to reveal the physical underpinnings of detector design that impact tracking quality. The study employs a generalizable ideal observer methodology for the quantification of tumor tracking performance. The analysis is applied to study both the effect of increasing scintillator thickness on a standard, single-layer imager (SLI) design as well as the effect of MLI architecture on tracking performance. METHODS: The present study uses the ideal observer signal-to-noise ratio (d') as a surrogate for tracking performance. We employ functions which model clinically relevant tasks and generalized frequency-domain imaging metrics to connect image quality with tumor tracking. A detection task for relevant Cartesian shapes (i.e., spheres and cylinders) was used to quantitate trackability of cases employing fiducial markers. Automated lung tumor tracking algorithms often leverage the differences in benign and malignant lung tissue textures. These types of algorithms (e.g., soft-tissue localization - STiL) were simulated by designing a discrimination task, which quantifies the differentiation of tissue textures, measured experimentally and fit as a power-law in trend (with exponent ß) using a cohort of MV images of patient lungs. The modeled MTF and NPS were used to investigate the effect of scintillator thickness and MLI architecture on tumor tracking performance. RESULTS: Quantification of MV images of lung tissue as an inverse power-law with respect to frequency yields exponent values of ß = 3.11 and 3.29 for benign and malignant tissues, respectively. Tracking performance with and without fiducials was found to be generally limited by quantum noise, a factor dominated by quantum detective efficiency (QDE). For generic SLI construction, increasing the scintillator thickness (gadolinium oxysulfide - GOS) from a standard 290 µm to 1720 µm reduces noise to about 10%. However, 81% of this reduction is appreciated between 290 and 1000 µm. In comparing MLI and SLI detectors of equivalent individual GOS layer thickness, the improvement in noise is equal to the number of layers in the detector (i.e., 4) with almost no difference in MTF. Further, improvement in tracking performance was slightly less than the square-root of the reduction in noise, approximately 84-90%. In comparing an MLI detector with an SLI with a GOS scintillator of equivalent total thickness, improvement in object detectability is approximately 34-39%. CONCLUSIONS: We have presented a novel method for quantification of tumor tracking quality and have applied this model to evaluate the performance of SLI and MLI EPID designs. We showed that improved tracking quality is primarily limited by improvements in NPS. When compared to very thick scintillator SLI, employing MLI architecture exhibits the same gains in QDE, but by mitigating the effect of optical Swank noise, results in more dramatic improvements in tracking performance.

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