Markerless tumor tracking using fast-kV switching dual-energy fluoroscopy on a benchtop system

Single-shot quantitative x-ray imaging using a primary modulator and dual-layer detector

BACKGROUND

Conventional x-ray imaging and fluoroscopy have limitations in quantitation due to several challenges, including scatter, beam hardening, and overlapping tissues. Dual-energy (DE) imaging, with its capability to quantify area density of specific materials, is well-suited to address such limitations, but only if the dual-energy projections are acquired with perfect spatial and temporal alignment and corrected for scatter.

PURPOSE

In this work, we propose single-shot quantitative imaging (SSQI) by combining the use of a primary modulator (PM) and dual-layer (DL) detector, which enables motion-free DE imaging with scatter correction in a single exposure.

METHODS

The key components of our SSQI setup include a PM and DL detector, where the former enables scatter correction for the latter while the latter enables beam hardening correction for the former. The SSQI algorithm allows simultaneous recovery of two material-specific images and two scatter images using four sub-measurements from the PM encoding. The concept was first demonstrated using simulation of chest x-ray imaging for a COVID patient. For validation, we set up SSQI geometry on our tabletop system and imaged acrylic and copper slabs with known thicknesses (acrylic: 0-22.5 cm; copper: 0-0.9 mm), estimated scatter with our SSQI algorithm, and compared the material decomposition (MD) for different combinations of the two materials with ground truth. Second, we imaged an anthropomorphic chest phantom containing contrast in the coronary arteries and compared the MD with and without SSQI. Lastly, to evaluate SSQI in dynamic applications, we constructed a flow phantom that enabled dynamic imaging of iodine contrast.

RESULTS

Our simulation study demonstrated that SSQI led to accurate scatter correction and MD, particularly for smaller focal blur and finer PM pitch. In the validation study, we found that the root mean squared error (RMSE) of SSQI estimation was 0.13 cm for acrylic and 0.04 mm for copper. For the anthropomorphic phantom, direct MD resulted in incorrect interpretation of contrast and soft tissue, while SSQI successfully distinguished them quantitatively, reducing RMSE in material-specific images by 38%-92%. For the flow phantom, SSQI was able to perform accurate dynamic quantitative imaging, separating contrast from the background.

CONCLUSIONS

We demonstrated the potential of SSQI for robust quantitative x-ray imaging. The integration of SSQI is straightforward with the addition of a PM and upgrade to a DL detector, which may enable its widespread adoption, including in techniques such as radiography and dynamic imaging (i.e., real-time image guidance and cone-beam CT).

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.

Effect of different noise reduction techniques and template matching parameters on markerless tumor tracking using dual-energy imaging

PURPOSE

To evaluate the impact of various noise reduction algorithms and template matching parameters on the accuracy of markerless tumor tracking (MTT) using dual-energy (DE) imaging.

METHODS

A Varian TrueBeam linear accelerator was used to acquire a series of alternating 60 and 120 kVp images (over a 180° arc) using fast kV switching, on five early-stage lung cancer patients. Subsequently, DE logarithmic weighted subtraction was performed offline on sequential images to remove bone. Various noise reduction techniques-simple smoothing, anticorrelated noise reduction (ACNR), noise clipping (NC), and NC-ACNR-were applied to the resultant DE images. Separately, tumor templates were generated from the individual planning CT scans, and band-pass parameter settings for template matching were varied. Template tracking was performed for each combination of noise reduction techniques and templates (based on band-pass filter settings). The tracking success rate (TSR), root mean square error (RMSE), and missing frames (percent unable to track) were evaluated against the estimated ground truth, which was obtained using Bayesian inference.

RESULTS

DE-ACNR, combined with template band-pass filter settings of σ_low  = 0.4 mm and σ_high  = 1.6 mm resulted in the highest TSR (87.5%), RMSE (1.40 mm), and a reasonable amount of missing frames (3.1%). In comparison to unprocessed DE images, with optimized band-pass filter settings of σ_low  = 0.6 mm and σ_high  = 1.2 mm, the TSR, RMSE, and missing frames were 85.3%, 1.62 mm, and 2.7%, respectively. Optimized band-pass filter settings resulted in improved TSR values and a lower missing frame rate for both unprocessed DE and DE-ACNR as compared to the use previously published band-pass parameters based on single energy kV images.

CONCLUSION

Noise reduction strategies combined with the optimal selection of band-pass filter parameters can improve the accuracy and TSR of MTT for lung tumors when using DE imaging.

Improved Tumor Image Estimation in X-Ray Fluoroscopic Images by Augmenting 4DCT Data for Radiotherapy

Measurement of tumor position is important for the radiotherapy of lung tumors with respiratory motion. Although tumors can be observed using X-ray fluoroscopy during radiotherapy, it is often difficult to measure tumor position from X-ray image sequences accurately because of overlapping organs. To measure tumor position accurately, a method for extracting tumor intensities from X-ray image sequences using a hidden Markov model (HMM) has been proposed. However, the performance of tumor intensity extraction depends on limited knowledge regarding the tumor motion observed in the four-dimensional computed tomography (4DCT) data used to construct the HMM. In this study, we attempted to improve the performance of tumor intensity extraction by augmenting 4DCT data. The proposed method was tested using simulated datasets of X-ray image sequences. The experimental results indicated that the HMM using the augmentation method could improve tumor-tracking performance when the range of tumor movement during treatment differed from that in the 4DCT data.

Adaptive dual-energy algorithm based on pre-calibrated weighting factors for chest radiography

Objective. To develop a dual-energy (DE) algorithm with spatially varying weighting factors for material selection and noise suppression. Approach. Calibration step-phantoms, with overlapping slabs of solid water and bone with different thicknesses, were used to obtain the pre-calibrated material selection and noise reduction weighting factors. The Material selection weighting factors were calculated by finding a zero of contrast-to-noise-ratio (CNR) between regions with two overlapping materials and regions of only target material, while noise suppression weighting factors were determined by maximizing signal-to-noise ratio for overlapping regions. The pre-calibrated weighting factors were fitted with low and high energy radiograph of two Rando phantoms to create maps of material selection and noise suppression weighting factors, which used with DE algorithm and anti-correlated noise reduction (ACNR) algorithm to generate DE images. Three different implementations, including two different sizes of Rando phantoms and two different orientations (oblique and anterior-posterior), were investigated. Soft-tissue and bone only images of Rando phantoms were obtained with five combinations of DE algorithms and CNR, contrast, and noise values of selected regions of interest were compared to evaluate the performance of the novel method: simple log subtraction (SLS), SLS with uniform ACNR, adaptive DE (aDE), aDE with uniform ACNR, and aDE and adaptive ACNR (aACNR). Main results. Compared to SLS, the aDE algorithm demonstrated improved image quality in all three orientations. CNR increased with better contrast for both soft-tissue and bone images. Implementation of aACNR algorithm resulted in further reduction of image noise and improvements in CNR at the cost of contrast. However, aACNR algorithm showed better contrast compared to ACNR method. Significance. A novel DE algorithm was proposed, which showed improved material selection and noise suppression as compared to the conventional DE techniques and can be easily implemented in a clinical environment for real-time DE image generation.

The markerless lung target tracking AAPM Grand Challenge (MATCH) results

PURPOSE

Lung stereotactic ablative body radiotherapy (SABR) is a radiation therapy success story with level 1 evidence demonstrating its efficacy. To provide real-time respiratory motion management for lung SABR, several commercial and preclinical markerless lung target tracking (MLTT) approaches have been developed. However, these approaches have yet to be benchmarked using a common measurement methodology. This knowledge gap motivated the MArkerless lung target Tracking CHallenge (MATCH). The aim was to localize lung targets accurately and precisely in a retrospective in silico study and a prospective experimental study.

METHODS

MATCH was an American Association of Physicists in Medicine sponsored Grand Challenge. Common materials for the in silico and experimental studies were the experiment setup including an anthropomorphic thorax phantom with two targets within the lungs, and a lung SABR planning protocol. The phantom was moved rigidly with patient-measured lung target motion traces, which also acted as ground truth motion. In the retrospective in silico study a volumetric modulated arc therapy treatment was simulated and a dataset consisting of treatment planning data and intra-treatment kilovoltage (kV) and megavoltage (MV) images for four blinded lung motion traces was provided to the participants. The participants used their MLTT approach to localize the moving target based on the dataset. In the experimental study, the participants received the phantom experiment setup and five patient-measured lung motion traces. The participants used their MLTT approach to localize the moving target during an experimental SABR phantom treatment. The challenge was open to any participant, and participants could complete either one or both parts of the challenge. For both the in silico and experimental studies the MLTT results were analyzed and ranked using the prospectively defined metric of the percentage of the tracked target position being within 2 mm of the ground truth.

RESULTS

A total of 30 institutions registered and 15 result submissions were received, four for the in silico study and 11 for the experimental study. The participating MLTT approaches were: Accuray CyberKnife (2), Accuray Radixact (2), BrainLab Vero, C-RAD, and preclinical MLTT (5) on a conventional linear accelerator (Varian TrueBeam). For the in silico study the percentage of the 3D tracking error within 2 mm ranged from 50% to 92%. For the experimental study, the percentage of the 3D tracking error within 2 mm ranged from 39% to 96%.

CONCLUSIONS

A common methodology for measuring the accuracy of MLTT approaches has been developed and used to benchmark preclinical and commercial approaches retrospectively and prospectively. Several MLTT approaches were able to track the target with sub-millimeter accuracy and precision. The study outcome paves the way for broader clinical implementation of MLTT. MATCH is live, with datasets and analysis software being available online at https://www.aapm.org/GrandChallenge/MATCH/ to support future research.

Markerless 3D tumor tracking during single-fraction free-breathing 10MV flattening-filter-free stereotactic lung radiotherapy

BACKGROUND AND PURPOSE

Positional verification during single fraction lung SBRT could increase confidence and reduce the chance of geographic miss. As planar 2DkV imaging during VMAT irradiation is already available on current linear accelerators, markerless tracking based on these images could offer widely available and low-cost verification. We evaluated treatment delivery data and template matching and triangulation for 3D-positional verification during free-breathing, single fraction (34 Gy), 10 MV flattening-filter-free VMAT lung SBRT.

METHODS AND MATERIALS

Tumor tracking based on kV imaging at 7 frames/second was performed during irradiation in 6 consecutive patients (7 lesions). Tumor characteristics, tracking ability, comparison of tracking displacements with CBCT-based shifts, tumor position relative to the PTV margin, and treatment times are reported.

RESULTS

For all 7 lesions combined, 3D tumor position could be determined for, on average, 71% (51-84%) of the total irradiation time. Visually estimated tracked and automated match +/- manually-corrected CBCT-derived displacements generally agreed within 1 mm. During the tracked period, the longitudinal, lateral and vertical position of the tumor was within a 5 mm/3 mm PTV margin 95.5/85.3% of the time. The PTV was derived from the ITV including all tumor motion. The total time from first set-up imaging to end of the last arc was 18.3-31.4 min (mean = 23.4, SD = 4.1).

CONCLUSION

3D positional verification during irradiation of small lung targets with limited motion, was feasible. However, tumor position could not be determined for on average 29% of the time. Improvements are needed. Margin reduction may be feasible. Imaging and delivery of a single 34 Gy fraction was fast.

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.

Markerless Real-Time 3-Dimensional kV Tracking of Lung Tumors During Free Breathing Stereotactic Radiation Therapy

Purpose

Accurate verification of tumor position during irradiation could reduce the probability of target miss. We investigated whether a commercial gantry-mounted 2-dimensional (2D) kilo-voltage (kV) imaging system could be used for real-time 3D tumor tracking during volumetric modulated arc therapy (VMAT) lung stereotactic body radiation therapy (SBRT). Markerless tumor tracking on kV fluoroscopic images was validated using a life-like moving thorax phantom and subsequently performed on kV images continuously acquired before and during free-breathing VMAT lung SBRT.

Methods and Materials

The 3D-printed/molded phantom containing 3 lung tumors was moved in 3D in TrueBeam developer mode, using simulated regular/irregular breathing patterns. Planar kV images were acquired at 7 frames/s during 11 Gy/fraction 10 MV flattening filter free VMAT. 2D reference templates were created for each gantry angle using the planning 4D computed tomography inspiration phase. kV images and templates were matched using normalized cross correlation to determine 2D tumor position, and triangulation of 2D matched projections determined the third dimension. 3D target tracking performed on cone beam computed tomography projection data from 18 patients (20 tumors) and real-time online tracking data from 2 of the 18 patients who underwent free-breathing VMAT lung SBRT are presented.

Results

For target 1 and 2 of the phantom (upper lung and middle/medial lung, mean density –130 Hounsfield units), 3D results within 2 mm of the known position were present in 92% and 96% of the kV projections, respectively. For target 3 (inferior lung, mean density –478 Hounsfield units) this dropped to 80%. Benchmarking against the respiratory signal, 13/20 (65%) tumors (10.5 ± 11.1 cm³) were considered successfully tracked on the cone beam computed tomography data. Tracking was less successful (≤50% of the time) in 7/20 (1.2 ± 1.5 cm³). Successful online tracking during lung SBRT was demonstrated.

Conclusions

3D markerless tumor tracking on a standard linear accelerator using template matching and triangulation of free-breathing kV fluoroscopic images was possible in 65% of small lung tumors. The smallest tumors were most challenging.

Densely sampled spectral modulation for x-ray CT using a stationary modulator with flying focal spot: a conceptual and feasibility study of scatter and spectral correction

PURPOSE

Modulation of the x-ray source in computed tomography (CT) by a designated filter to achieve a desired distribution of photon flux has been greatly advanced in recent years. In this work, we present a densely sampled spectral modulation (DSSM) as a promising low-cost solution to quantitative CT imaging in the presence of scatter. By leveraging a special stationary filter (namely a spectral modulator) and a flying focal spot, DSSM features a strong correlation in the scatter distributions across focal spot positions and sees no substantial projection sparsity or misalignment in data sampling, making it possible to simultaneously correct for scatter and spectral effects in a unified framework.

METHODS

The concept of DSSM is first introduced, followed by an analysis of the design and benefits of using the stationary spectral modulator with a flying focal spot (SMFFS) that dramatically changes the data sampling and its associated data processing. With an assumption that the scatter distributions across focal spot positions have strong correlation, a scatter estimation and spectral correction algorithm from DSSM is then developed, where a dual-energy modulator along with two flying focal spot positions is of interest. Finally, a phantom study on a tabletop cone-beam CT system is conducted to understand the feasibility of DSSM by SMFFS, using a copper modulator and by moving the x-ray tube position in the X direction to mimic the flying focal spot.

RESULTS

Based on our analytical analysis of the DSSM by SMFFS, the misalignment of low- and high-energy projection rays can be reduced by a factor of more than 10 when compared with a stationary modulator only. With respect to modulator design, metal materials such as copper, molybdenum, silver, and tin could be good candidates in terms of energy separation at a given attenuation of photon flux. Physical experiments using a Catphan phantom as well as an anthropomorphic chest phantom demonstrate the effectiveness of DSSM by SMFFS with much better CT number accuracy and less image artifacts. The root mean squared error was reduced from 297.9 to 6.5 Hounsfield units (HU) for the Catphan phantom and from 409.3 to 39.2 HU for the chest phantom.

CONCLUSIONS

The concept of DSSM using a SMFFS is proposed. Phantom results on its scatter estimation and spectral correction performance validate our main ideas and key assumptions, demonstrating its potential and feasibility for quantitative CT imaging.

Clinical practice vs. state-of-the-art research and future visions: Report on the 4D treatment planning workshop for particle therapy - Edition 2018 and 2019

The 4D Treatment Planning Workshop for Particle Therapy, a workshop dedicated to the treatment of moving targets with scanned particle beams, started in 2009 and since then has been organized annually. The mission of the workshop is to create an informal ground for clinical medical physicists, medical physics researchers and medical doctors interested in the development of the 4D technology, protocols and their translation into clinical practice. The 10th and 11th editions of the workshop took place in Sapporo, Japan in 2018 and Krakow, Poland in 2019, respectively. This review report from the Sapporo and Krakow workshops is structured in two parts, according to the workshop programs. The first part comprises clinicians and physicists review of the status of 4D clinical implementations. Corresponding talks were given by speakers from five centers around the world: Maastro Clinic (The Netherlands), University Medical Center Groningen (The Netherlands), MD Anderson Cancer Center (United States), University of Pennsylvania (United States) and The Proton Beam Therapy Center of Hokkaido University Hospital (Japan). The second part is dedicated to novelties in 4D research, i.e. motion modelling, artificial intelligence and new technologies which are currently being investigated in the radiotherapy field.

Adaptive noise reduction for dual-energy x-ray imaging based on spatial variations in beam attenuation

PURPOSE

The main goal of this work is to improve the previously proposed patient-specific pixel-based dual-energy (PP-DE) algorithm by developing an adaptive anti-correlated noise reduction (ACNR) method, resulting in reduced image noise.

METHODS

Theoretical models of contrast-to-noise (CNR) and signal-to-noise (SNR) ratio were developed as functions of weighting factors for DE bone ω _Bn or soft tissue ω _ST cancellation. These analytical expressions describe CNR and SNR properties of dual-energy (DE) images, obtained with both simple log subtraction (SLS) and ACNR algorithms, and allow for a direct comparison between experimental and theoretical results. The theoretical models demonstrate the importance of ACNR weighting factor (ω _A ) optimization leading to the maximization of the SNR of the final image. A step phantom was constructed, which consisted of overlapping slabs of solid water (0-30 cm) and bone-mimicking material (0-6 cm), resulting in a total of 7 × 7 regions. High-energy (HE) and low-energy (LE) images were acquired at 140 kVp and 60 kVp with a clinical ExacTrac imaging system. The CNR and SNR were obtained for the DE images as functions of ω _Bn,ST and noise reduction weighting factor ω _A for different combinations of thicknesses. Weighting factors for bone cancellation were optimized for each region of interest (ROI) by finding zeros of CNR function for DE images between soft tissue only and soft tissue plus bone regions (and vice versa for soft tissue cancellation). The weighting factor for the ACNR algorithm ω _A was then optimized by maximizing the SNR function for each ROI. HE and LE images for an anthropomorphic Rando phantom were obtained with the same acquisition parameters as for the step phantom. DE images for bone only and soft tissue only were obtained with three algorithms: SLS and PP-DE with conventional ACNR (uniform ω _A ), and PP-DE with adaptive ACNR (region-varying ω _A ). Weighting factor maps for PP-DE and adaptive ACNR methods were obtained for Rando phantom geometry (which was determined from its CT scans) by interpolation (or extrapolation) of weighting factors for the step phantom. CNR values were calculated for different regions.

RESULTS

The CNR and SNR characteristics as functions of material cancellation and noise reduction weighting factors were obtained from theoretical models and experimental data from the step phantom. This showed a good qualitative validation of the models. For the ANCR algorithm, both the theory and experiment demonstrated that the material cancellation weighting factors (ω _Bn,ST ) can be optimized independently of the noise cancellation weighting factors (ω _A ), which can be optimized by maximizing SNR. For each ROI (with different overlapping bone and soft tissue thicknesses) the weighting factors ω _Bn,ST were determined as well as corresponding optimal weighting factors ω _A for noise reduction. For the Rando phantom, CNR values for regions representing different anatomical structures (ribs, spine, and tumor) were evaluated. It was shown that the proposed adaptive ACNR further improves image quality, compared to the conventional ACNR algorithm. The improvement is maximized for regions with bones (ribs or spine), where the largest attenuation is observed.

CONCLUSION

The ACNR weighting factors are dependent on the material thicknesses due to varying beam attenuation leading to different levels of quantum noise. This was shown with the derived theoretical expressions of the CNR and SNR functions and was validated by experimental data. The adaptive ANCR DE algorithm was developed, which allows for an increase in image quality by spatially varying weighting factors for noise reduction. This algorithm complements the previously developed PP-DE algorithm to obtain better quality DE images.

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.

The first prospective implementation of markerless lung target tracking in an experimental quality assurance procedure on a standard linear accelerator

The ability to track tumour motion without implanted markers on a standard linear accelerator (linac) could enable wide access to real-time adaptive radiotherapy for cancer patients. We previously have retrospectively validated a method for 3D markerless target tracking using intra-fractional kilovoltage (kV) projections acquired on a standard linac. This paper presents the first prospective implementation of markerless lung target tracking on a standard linac and its quality assurance (QA) procedure. The workflow and the algorithm developed to track the 3D target position during volumetric modulated arc therapy treatment delivery were optimised. The linac was operated in clinical QA mode, while kV projections were streamed to a dedicated computer using a frame-grabber software. The markerless target tracking accuracy and precision were measured in a lung phantom experiment under the following conditions: static localisation of seven distinct positions, dynamic localisation of five patient-measured motion traces, and dynamic localisation with treatment interruption. The QA guidelines were developed following the AAPM Task Group 147 report with the requirement that the tracking margin components, the margins required to account for tracking errors, did not exceed 5 mm in any direction. The mean tracking error ranged from 0.0 to 0.9 mm (left-right), -0.6 to  -0.1 mm (superior-inferior) and  -0.7 to 0.1 mm (anterior-posterior) over the three tests. Larger errors were found in cases with large left-right or anterior-posterior and small superior-inferior motion. The tracking margin components did not exceed 5 mm in any direction and ranged from 0.4 to 3.2 mm (left-right), 0.7 to 1.6 mm (superior-inferior) and 0.8 to 1.5 mm (anterior-posterior). This study presents the first prospective implementation of markerless lung target tracking on a standard linac and provides a QA procedure for its safe clinical implementation, potentially enabling real-time adaptive radiotherapy for a large population of lung cancer patients.

Fast-switching dual energy cone beam computed tomography using the on-board imager of a commercial linear accelerator

To evaluate fast-kV switching (FS) dual energy (DE) cone beam computed tomography (CBCT) using the on-board imager (OBI) of a commercial linear accelerator to produce virtual monoenergetic (VM) and relative electron density (RED) images. Using an polynomial attenuation mapping model, CBCT phantom projections obtained at 80 and 140 kVp with FS imaging, were decomposed into equivalent thicknesses of aluminum (Al) and polymethyl methacrylate (PMMA). All projections were obtained with the titanium foil and bowtie filter in place. Basis material projections were then recombined to create VM images by using the linear attenuation coefficients at the specified energy for each material. Similarly, RED images were produced by replacing the linear attenuation values of Al and PMMA by their respective RED values in the projection space. VM and RED images were reconstructed using Feldkamp-Davis-Kress (FDK) and an iterative algorithm (iCBCT, Varian Medical Systems). Hounsfield units (HU), contrast-to-noise ratio (CNR) and RED values were compared against known values. The results after VM-CBCT production showed good material decomposition and consistent HU_VM values, with measured root mean square errors (RMSE) from theoretical values, after FDK reconstruction, of 20.5, 5.7, 12.8 and 21.7 HU for 50, 80, 100 and 150 keV, respectively. The largest CNR improvements, when compared to polychromatic images, were observed for the 50 keV VM images. Image noise was reduced up to 28% in the VM-CBCT images after iterative image reconstruction. RED values measured for our method resulted in a mean percentage error of 0.0%  ±  1.8%. This study describes a method to generate VM-CBCT and RED images using FS-DE scans obtained using the OBI of a linac, including the effects of the bowtie filter. The creation of VM and RED images increases the dynamic range of CBCT images, and provides additional data that may be used for adaptive radiotherapy, and on table verification for radiotherapy treatments.

Technical Principles of Dual-Energy Cone Beam Computed Tomography and Clinical Applications for Radiation Therapy

PURPOSE

Medical imaging is an indispensable tool in radiotherapy for dose planning, image guidance and treatment monitoring. Cone beam CT (CBCT) is a low dose imaging technique with high spatial resolution capability as a direct by-product of using flat-panel detectors. However, certain issues such as x-ray scatter, beam hardening and other artifacts limit its utility to the verification of patient positioning using image-guided radiotherapy.

METHODS AND MATERIALS

Dual-energy (DE)-CBCT has recently demonstrated promise as an improved tool for tumor visualization in benchtop applications. It has the potential to improve soft-tissue contrast and reduce artifacts caused by beam hardening and metal. In this review, the practical aspects of developing a DE-CBCT based clinical and technical workflow are presented based on existing DE-CBCT literature and concepts adapted from the well-established library of work in DE-CT. Furthermore, the potential applications of DE-CBCT on its future role in radiotherapy are discussed.

RESULTS AND CONCLUSIONS

Based on current literature and an investigation of future applications, there is a clear potential for DE-CBCT technologies to be incorporated into radiotherapy. The applications of DE-CBCT include (but are not limited to): adaptive radiotherapy, brachytherapy, proton therapy, radiomics and theranostics.

Adaptive weighted log subtraction based on neural networks for markerless tumor tracking using dual-energy fluoroscopy

PURPOSE

To present a novel method, based on convolutional neural networks (CNN), to automate weighted log subtraction (WLS) for dual-energy (DE) fluoroscopy to be used in conjunction with markerless tumor tracking (MTT).

METHODS

A CNN was developed to automate WLS (aWLS) of DE fluoroscopy to enhance soft tissue visibility. Briefly, this algorithm consists of two phases: training a CNN architecture to predict pixel-wise weighting factors followed by application of WLS subtraction to reduce anatomical noise. To train the CNN, a custom phantom was built consisting of aluminum (Al) and acrylic (PMMA) step wedges. Per-pixel ground truth (GT) weighting factors were calculated by minimizing the contrast of Al in the step wedge phantom to train the CNN. The pretrained model was then utilized to predict pixel-wise weighting factors for use in WLS. For comparison, the weighting factor was manually determined in each projection (mWLS). A thorax phantom with five simulated spherical targets (5-25 mm) embedded in a lung cavity, was utilized to assess aWLS performance. The phantom was imaged with fast-kV dual-energy (120 and 60 kVp) fluoroscopy using the on-board imager of a commercial linear accelerator. DE images were processed offline to produce soft tissue images using both WLS methods. MTT was compared using soft tissue images produced with both mWLS and aWLS techniques.

RESULTS

Qualitative evaluation demonstrated that both methods achieved soft tissue images with similar quality. The use of aWLS increased the number of tracked frames by 1-5% compared to mWLS, with the largest increase observed for the smallest simulated tumors. The tracking errors for both methods produced agreement to within 0.1 mm.

CONCLUSIONS

A novel method to perform automated WLS for DE fluoroscopy was developed. Having similar soft tissue quality as well as bone suppression capability as mWLS, this method allows for real-time processing of DE images for MTT.