3D EPID-based in vivo dosimetry for IMRT and VMAT

IPEM topical report: guidance for the clinical implementation of online treatment monitoring solutions for IMRT/VMAT

This report provides guidance for the implementation of online treatment monitoring (OTM) solutions in radiotherapy (RT), with a focus on modulated treatments. Support is provided covering the implementation process, from identification of an OTM solution to local implementation strategy. Guidance has been developed by a RT special interest group (RTSIG) working party (WP) on behalf of the Institute of Physics and Engineering in Medicine (IPEM). Recommendations within the report are derived from the experience of the WP members (in consultation with manufacturers, vendors and user groups), existing guidance or legislation and a UK survey conducted in 2020 (Stevenset al2021). OTM is an inclusive term representing any system capable of providing a direct or inferred measurement of the delivered dose to a RT patient. Information on each type of OTM is provided but, commensurate with UK demand, guidance is largely influenced byin vivodosimetry methods utilising the electronic portal imager device (EPID). Sections are included on the choice of OTM solutions, acceptance and commissioning methods with recommendations on routine quality control, analytical methods and tolerance setting, clinical introduction and staffing/resource requirements. The guidance aims to give a practical solution to sensitivity and specificity testing. Functionality is provided for the user to introduce known errors into treatment plans for local testing. Receiver operating characteristic analysis is discussed as a tool to performance assess OTM systems. OTM solutions can help verify the correct delivery of radiotherapy treatment. Furthermore, modern systems are increasingly capable of providing clinical decision-making information which can impact the course of a patient's treatment. However, technical limitations persist. It is not within the scope of this guidance to critique each available solution, but the user is encouraged to carefully consider workflow and engage with manufacturers in resolving compatibility issues.

Assessment of a Therapeutic X-ray Radiation Dose Measurement System Based on a Flexible Copper Indium Gallium Selenide Solar Cell

Several detectors have been developed to measure radiation doses during radiotherapy. However, most detectors are not flexible. Consequently, the airgaps between the patient surface and detector could reduce the measurement accuracy. Thus, this study proposes a dose measurement system based on a flexible copper indium gallium selenide (CIGS) solar cell. Our system comprises a customized CIGS solar cell (with a size 10 × 10 cm2 and thickness 0.33 mm), voltage amplifier, data acquisition module, and laptop with in-house software. In the study, the dosimetric characteristics, such as dose linearity, dose rate independence, energy independence, and field size output, of the dose measurement system in therapeutic X-ray radiation were quantified. For dose linearity, the slope of the linear fitted curve and the R-square value were 1.00 and 0.9999, respectively. The differences in the measured signals according to changes in the dose rates and photon energies were <2% and <3%, respectively. The field size output measured using our system exhibited a substantial increase as the field size increased, contrary to that measured using the ion chamber/film. Our findings demonstrate that our system has good dosimetric characteristics as a flexible in vivo dosimeter. Furthermore, the size and shape of the solar cell can be easily customized, which is an advantage over other flexible dosimeters based on an a-Si solar cell.

A regression-based approach to compute the pixels sensitivity map of linear accelerator portal imaging devices

PURPOSE

To determine the pixel sensitivity map (PSM) for amorphous silicon electronic portal imaging devices (EPIDs) using a single flood field signal.

METHOD AND MATERIALS

A raw EPID signal results from the incident particle energy fluence, the inherent pixels response, and the background signal. In large open fields, particle energy fluence is a slow-varying signal that is locally considered spatially constant. Pixels response is a fast and abrupt varying behavior. The background signal is due to the EPID panel electronics, which is determined during radiation absence. To determine the PSM, after correcting for the background signal, we apply a model that captures the underlying smooth particle energy fluence-induced signal. This fluence signal-fitted model is then used to determine the PSM. Here, we use a polynomial-based regression surface model in both x and y dimensions. To validate the generated PSM, we measure beams and compute PSMs for multiple beam energies with and without flattening filters and for multiple source-to-imager distances. Since the PSM is a detector characteristic, it should be independent of those variables. We also intercompare measurements of fixed slit fields with the EPID being shifted between measurements.

RESULTS

The fluence signal of the flattening filter-free (FFF) beams was optimally modeled as a 12th degree polynomial surfaces, which had ≤ 0.1% residuals near the central axis. The 6 and 10 MV FFF PSMs were within ˜0.1%, and independent of the EPID SID, suggesting that the PSM is energy independent. The 6, 10, and 15 MV flattened-beam PSMs were well modeled as 12th degree polynomial surfaces, which were equivalent within ˜0.24% but differed from the FFF PSM by up to 0.5% near the beam central axis. Applying the FFF PSMs to the flattened-beam measurements reduced the central-axis deviation between the raw and corrected signal to < 0.1%, confirming the PSM energy independence hypothesis. When the FFF PSM is utilized, output verification with shifted slit deliveries agreed within ˜0.5% for all beam energies, which is within the radiation delivery uncertainty of ˜0.57%.

CONCLUSION

PSM for MV EPIDs can be determined by separating out the slowly varying, well-behaved fluence signal from the pixel-to-pixel sensitivity variations. The quality of the PSM is found to be dependent on the quality of the surface fit, which is best for the 6 MV FFF beam measured at SID equal to 180 cm. Within fitting errors, the PSM is independent of beam energy for 6, 10, and 15 MV beams with and without flattening filters. The PSM generation does not require shifting the EPID panel nor multiple EPID panel irradiations and should be usable for linacs with fixed geometry EPIDs.

Portal dosimetry in radiotherapy repeatability evaluation

The accuracy of radiotherapy is the subject of continuous discussion, and dosimetry methods, particularly in dynamic techniques, are being developed. At the same time, many oncology centers develop quality procedures, including pretreatment and online dose verification and proper patient tracking methods. This work aims to present the possibility of using portal dosimetry in the assessment of radiotherapy repeatability. The analysis was conducted on 74 cases treated with dynamic techniques. Transit dosimetry was made for each collision‐free radiation beam. It allowed the comparison of summary fluence maps, obtained for fractions with the corresponding summary maps from all other treatment fractions. For evaluation of the compatibility in the fluence map pairs (6798), the gamma coefficient was calculated. The results were considered in four groups, depending on the used radiotherapy technique: stereotactic fractionated radiotherapy, breath‐hold, free‐breathing, and conventionally fractionated other cases. The chi² or Fisher's exact test was made depending on the size of the analyzed set and also Mann–Whitney U‐test was used to compare treatment repeatability of different techniques. The aim was to test whether the null hypothesis of error‐free therapy was met. The patient is treated repeatedly if the P‐value in all the fluence maps sets is higher than the level of 0.01. The best compatibility between treatment fractions was obtained for the stereotactic technique. The technique with breath‐holding gave the lowest percentage of compliance of the analyzed fluence pairs. The results indicate that the repeatability of the treatment is associated with the radiotherapy technique. Treated volume location is also an essential factor found in the evaluation of treatment accuracy. The EPID device is a useful tool in assessing the repeatability of radiotherapy. The proposed method of fluence maps comparison also allows us to assess in which therapeutic session the patient was treated differently from the other fractions.

A tri-hybrid method to estimate the patient-generated scattered photon fluence components to the EPID image plane

In vivo dosimetry methods can verify the prescription dose is delivered to the patient during treatment. Unfortunately, in exit dosimetry, the megavoltage image is contaminated with patient-generated scattered photons. However, estimation and removal of the effect of this fluence improves accuracy of in vivo dosimetry methods. This work develops a 'tri-hybrid' algorithm combining analytical, Monte Carlo (MC) and pencil-beam scatter kernel methods to provide accurate estimates of the total patient-generated scattered photon fluence entering the MV imager. For the multiply-scattered photon fluence, a modified MC simulation method was applied, using only a few histories. From each second- and higher-order interaction site in the simulation, energy fluence entering all pixels of the imager was calculated using analytical methods. For photon fluence generated by electron interactions in the patient (i.e. bremsstrahlung and positron annihilation), a convolution/superposition approach was employed using pencil-beam scatter fluence kernels as a function of patient thickness and air gap distance, superposed on the incident fluence distribution. The total patient-scattered photon fluence entering the imager was compared with a corresponding full MC simulation (EGSnrc) for several test cases. These included three geometric phantoms (water, half-water/half-lung, computed tomography thorax) using monoenergetic (1.5, 5.5 and 12.5 MeV) and polyenergetic (6 and 18 MV) photon beams, 10 × 10 cm² field, source-to-surface distance 100 cm, source-to-imager distance 150 cm and 40 × 40 cm² imager. The proposed tri-hybrid method is demonstrated to agree well with full MC simulation, with the average fluence differences and standard deviations found to be within 0.5% and 1%, respectively, for test cases examined here. The method, as implemented here with a single CPU (non-parallelized), takes ∼80 s, which is considerably shorter compared to full MC simulation (∼30 h). This is a promising method for fast yet accurate calculation of patient-scattered fluence at the imaging plane for in vivo dosimetry applications.

Simulation of x-ray-induced acoustic imaging for absolute dosimetry: Accuracy of image reconstruction methods

PURPOSE

Three-dimensional in-vivo dose verification is one of the standing challenges in radiation therapy. X-ray-induced acoustic tomography has recently been proposed as an imaging method for use in in-vivo dosimetry. The aim of this study was to investigate the accuracy of reconstructing three-dimensional (3D) absolute dose using x-ray-induced acoustic tomography. We performed this investigation using two different tomographic dose reconstruction techniques.

METHODS

Two examples of 3D dose reconstruction techniques for x-ray acoustic imaging are investigated. Dose distributions are calculated for varying field sizes using a clinical treatment planning system. The induced acoustic pressure waves which are generated by the increase in temperature due to the absorption of pulsed MV x-rays are simulated using an advanced numerical modeling package for acoustic wave propagation in the time domain. Two imaging techniques, back projection and iterative time reversal, are used to reconstruct the 3D dose distribution in a water phantom with open fields. Image analysis is performed and reconstructed depth dose curves from x-ray acoustic imaging are compared to the depth dose curves calculated from the treatment planning system. Calculated field sizes from the reconstructed dose profiles by back projection and time reversal are compared to the planned field size to determine their accuracy. The iterative time reversal imaging technique is also used to reconstruct dose in an example clinical dose distribution. Image analysis of this clinical test case is performed using the gamma passing rate. In addition, gamma passing rates are used to validate the stopping criteria in the iterative time reversal method.

RESULTS

Water phantom simulations showed that back projection does not adequately reconstruct the shape and intensity of the depth dose. When compared to the depth of maximum dose calculated by a treatment planning system, the maximum dose depth by back projection is shifted deeper by 55 and 75 mm for 4 × 4 cm and 10 × 10 cm field sizes, respectively. The reconstructed depth dose by iterative time reversal accurately agrees with the planned depth dose for a 4 × 4 cm field size and is shifted deeper by 12 mm for the 10 × 10 cm field size. When reconstructing field sizes, the back projection method leads to 18% and 35% larger sizes for the 4 × 4 cm and 10 × 10 cm fields, respectively, whereas the iterative time reversal method reconstructs both field sizes with < 2% error. For the clinical dose distribution, we were able to reconstruct the dose delivered by a 1 degree sub-arc with a good accuracy. The reconstructed and planned doses were compared using gamma analysis, with> 96% gamma passing rate at 3%/2 mm.

CONCLUSIONS

Our results show that the 3D x-ray acoustic reconstructed dose by iterative time reversal is considerably more accurate than the dose reconstructed by back projection. Iterative time reversal imaging has a potential for use in 3D absolute dosimetry.

EPID – a useful interfraction QC tool

Biomedical accelerators used in radiotherapy are equipped with detector arrays which are commonly used to obtain the image of patient position during the treatment session. These devices use both kilovolt and megavolt x-ray beams. The advantage of EPID (Electronic Portal Imaging Device) megavolt panels is the correlation of the measured signal with the calibrated dose. The EPID gives a possibility to verify delivered dose. The aim of the study is to answer the question whether EPID can be useful as a tool for interfraction QC (quality control) of dose and geometry repeatability.

The EPID system has been calibrated according to the manufacturer’s recommendations to obtain a signal and dose values correlation. Initially, the uncertainty of the EPID matrix measurement was estimated. According to that, the detecting sensitivity of two parameters was checked: discrepancies between the planned and measured dose and field geometry variance. Moreover, the linearity of measured signal-dose function was evaluated.

In the second part of the work, an analysis of several dose distributions was performed. In this study, the analysis of clinical cases was limited to stereotactic dynamic radiotherapy. Fluence maps were obtained as a result of the dose distribution measurements with the EPID during treatment sessions. The compatibility of fluence maps was analyzed using the gamma index. The fluence map acquired during the first fraction was the reference one. The obtained results show that EPID system can be used for interfraction control of dose and geometry repeatability.

The effect of the choice of patient model on the performance of in vivo 3D EPID dosimetry to detect variations in patient position and anatomy

PURPOSE

In vivo EPID dosimetry is meant to trigger on relevant differences between delivered and planned dose distributions and should therefore be sensitive to changes in patient position and patient anatomy. Three-dimensional (3D) EPID back-projection algorithms can use either the planning computed tomography (CT) or the daily patient anatomy as patient model for dose reconstruction. The purpose of this study is to quantify the effect of the choice of patient model on the performance of in vivo 3D EPID dosimetry to detect patient-related variations.

METHODS

Variations in patient position and patient anatomy were simulated by transforming the reference planning CT images (pCT) into synthetic daily CT images (dCT) representing a variation of a given magnitude in patient position or in patient anatomy. For each variation, synthetic in vivo EPID data were also generated to simulate the reconstruction of in vivo EPID dose distributions. Both the planning CT images and the synthetic daily CT images could be used as patient model in the reconstructions yielding e D pCT and e D dCT EPID reconstructed dose distributions respectively. The accuracy of e D pCT and e D dCT reconstructions was evaluated against absolute dose measurements made in different phantom setups, and against dose distributions calculated by the treatment planning system (TPS). The comparison was performed by γ-analysis (3% local dose/2 mm). The difference in sensitivity between e D pCT and e D dCT reconstructions to detect variations in patient position and in patient anatomy was investigated using receiver operating characteristic analysis and the number of triggered alerts for 100 volumetric modulated arc therapy plans and 12 variations.

RESULTS

e D dCT showed good agreement with both absolute point dose measurements (<0.5%) and TPS data (γ-mean = 0.52 ± 0.11). The agreement degraded with e D pCT , with the magnitude of the deviation varying with each specific case. e D dCT readily detected combined 3 mm translation setup errors in all directions (AUC = 1.0) and combined 3° rotation setup errors around all axes (AUC = 0.86) whereas e D pCT showed good detectability only for 12 mm translations (AUC = 0.85) and 9° rotations (AUC = 0.80). Conversely, e D pCT manifested a higher sensitivity to patient anatomical changes resulting in AUC values of 0.92/0.95 for a 6 mm patient contour expansion/contraction compared to 0.70/0.64 with e D dCT . Using |ΔPTV_D50 | > 3% as clinical tolerance level, the percentage of alerts for 6 mm changes in patient contour were 85%/27% with e D pCT / e D dCT .

CONCLUSIONS

With planning CT images as patient model, EPID dose reconstructions underestimate the dosimetric effects caused by errors in patient positioning and overestimate the dosimetric effects caused by changes in patient anatomy. The use of the daily patient position and anatomy as patient model for in vivo 3D EPID transit dosimetry improves the ability of the system to detect uncorrected errors in patient position and it reduces the likelihood of false positives due to patient anatomical changes.

A Monte Carlo study of the impact of phosphor optical properties on EPID imaging performance

We have developed a Monte Carlo computational model of a clinically employed electronic portal imaging device (EPID), and demonstrated the impact of phosphor optical properties on the imaging performance. The EPID model was built with Geant4 application for tomographic emission. Both radiative and optical transport were included in the model. Modulation transfer function (MTF), normalized noise-power spectrum times the incident x-ray fluence (qNNPS), and detective quantum efficiency (DQE) were calculated for simulated and measured data, and their agreement was quantified by the normalized root-mean-square error (NRMSE). MTF was computed using a 100 µm wide slit tilted by 1.5° and qNNPS was estimated using the Fujita-Lubberts-Swank method. DQE was calculated from MTF and qNNPS data. The NRMSE value was 0.0467 for MTF, 0.0217 for qNNPS, and 0.0885 for DQE, showing good agreement between measurement and simulation. Five major optical properties, phosphor grain size, phosphor thickness, phosphor refractive index, binder refractive index, and packing ratio were tested for their influence on the qNNPS, MTF, and DQE(0) of the model. Generally, the effect on the qNNPS is greater than MTF, and no impact on DQE(0), except from phosphor thickness, was observed. Multiple applications, such as imager design optimization and investigations of the dosimetric performance, are expected to benefit from the validated model.

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.