Phantom Setups Research Articles

Simultaneous iodine and barium imaging with photon-counting CT

Photon-counting computed tomography (PCCT) could potentially offer significant improvements in image quality over conventional CT. This study leverages PCCT technology to achieve simultaneous multi-contrast imaging of barium (Z=56) and iodine (Z=53), addressing the limitations of conventional CT in differentiating materials with similar atomic numbers. Using a bench-top PCCT system with a cadmium zinc telluride (CZT) detector, experiments were conducted with various phantom setups and 1-5\% contrast agent concentrations and in a biological sample. Energy thresholds were adjusted to the K-edge absorption energies of barium (37.4 keV) and iodine (33.2 keV) to capture multi-energy CT images. K-edge decomposition was performed with K-edge subtraction and principal component analysis (PCA) techniques, demonstrating clear differentiation and accurate quantification of contrast agents in phantom combinations and the biological sample. The PCCT system successfully differentiated and quantified barium and iodine, with high correlations ($R^2 \approx 1$) between true and reconstructed concentrations. PCA outperformed K-edge subtraction by providing a superior capability to differentiate barium and iodine, particularly when calcium was present in the scanned object. These findings highlight the potential of PCCT for reliable, detailed imaging in clinical and research applications, particularly for contrast agents with close atomic numbers.

Characterizing the voxel-based approaches in radioembolization dosimetry with reDoseMC.

Yttrium-90 ( ) represents the primary radioisotope used in radioembolization procedures, while holmium-166 ( ) is hypothesized to serve as a viable substitute for due to its comparable therapeutic potential and improved quantitative imaging. Voxel-based dosimetry for these radioisotopes relies on activity images obtained through PET or SPECT and dosimetry methods, including the voxel S-value (VSV) and the local deposition method (LDM). However, the evaluation of the accuracy of absorbed dose calculations has been limited by the use of non-ideal reference standards and investigations restricted to the liver. The objective of this study was to expand upon these dosimetry characterizations by investigating the impact of image resolutions, voxel sizes, target volumes, and tissue materials on the accuracy of and dosimetry techniques. A specialized radiopharmaceutical dosimetry software called reDoseMC was developed using the Geant4 Monte Carlo toolkit and validated by benchmarking the generated kernels with published data. The decay spectra of both and were also compared. Multiple VSV kernels were generated for the liver, lungs, soft tissue, and bone for isotropic voxel sizes of 1 mm, 2 mm, and 4 mm. Three theoretical phantom setups were created with 20 or 40mm activity and mass density inserts for the same three voxel sizes. To replicate the limited spatial resolutions present in PET and SPECT images, image resolutions were modeled using a 3D Gaussian kernel with a Full Width at Half Maximum (FWHM) ranging from 0 to 16mm and with no added noise. The VSV and LDM dosimetry methods were evaluated by characterizing their respective kernels and analyzing their absorbed dose estimates calculated on theoretical phantoms. The ground truth for these estimations was calculated usingreDoseMC. The decay spectra obtained through reDoseMC showed less than a 1% difference when compared to previously published experimental data for energies below 1.9MeV in the case of and less than 1% for energies below 1.5MeV for . Additionally, the validation kernels for VSV exhibited results similar to those found in published Monte Carlo codes, with source dose depositions having less than a 3% error margin. Resolution thresholds ( ), defined as resolutions that resulted in similar dose estimates between the LDM and VSV methods, were observed for . They were 1.5mm for bone, 2.5 mm for soft tissue and liver, and 8.5 mm for lungs. For , the accuracy of absorbed dose deposition was found to be dependent on the contributions of absorbed dose from photons. Volume errors due to variations in voxel size impacted the final dose estimates. Larger target volumes yielded more accurate mean doses than smaller volumes. For both radioisotopes, the radial dose profiles for the VSV and LDM approximated but never matched the referencestandard. reDoseMC was developed and validated for radiopharmaceutical dosimetry. The accuracy of voxel-based dosimetry was found to vary widely with changes in image resolutions, voxel sizes, chosen target volumes, and tissue material; hence, the standardization of dosimetry protocols was found to be of great importance for comparable dosimetryanalysis.

Gamma camera imaging characteristics of 166Ho and 99mTc used in Selective Internal Radiation Therapy

BackgroundThe administration of a 166Ho scout dose is available as an alternative to 99mTc particles for pre-treatment imaging in Selective Internal Radiation Therapy (SIRT). It has been reported that the 166Ho scout dose may be more accurate for the prediction of microsphere distribution and the associated therapy planning. The aim of the current study is to compare the scintigraphic imaging characteristics of both isotopes, considering the objectives of the pre-treatment imaging using clinically geared phantoms.MethodsPlanar and SPECT/CT images were obtained using a NEMA image quality phantom in different phantom setups and another body-shaped phantom with several inserts. The influence of collimator type, count statistics, dead time effects, isotope properties and patient obesity on spatial resolution, contrast recovery and the detectability of small activity accumulations was investigated. Furthermore, the effects of the imaging characteristics on personalized dosimetry are discussed.ResultsThe images with 99mTc showed up to 3 mm better spatial resolution, up to two times higher contrast recovery and significantly lower image noise than those with 166Ho. The contrast-to-noise ratio was up to five times higher for 99mTc than for 166Ho. Only when using 99mTc all activity-filled spheres could be distinguished from the activity-filled background. The measurements mimicking an obese patient resulted in a degraded image quality for both isotopes.ConclusionsOur measurements demonstrate better scintigraphic imaging properties for 99mTc compared to 166Ho in terms of spatial resolution, contrast recovery, image noise, and lesion detectability. While the 166Ho scout dose promises better prediction of the microsphere distribution, it is important to consider the inferior imaging characteristics of 166Ho, which may affect individualized treatment planning in SIRT.

Prediction of proton stopping power ratios using dual-energy CT basis material decomposition.

Proton radiotherapy treatment plans are currently restricted by the range uncertainties originating from the stopping power ratio (SPR) prediction based on single-energy computed tomography (SECT). Various studies have shown that multi-energy CT (MECT) can reduce the range uncertainties due to medical implant materials and age-related variations in tissue composition. None of these has directly applied the basis material density (MD) images produced by projection-based MECT systems for SPR prediction. To present and evaluate a novel proton SPR prediction method based on MD images from dual-energy CT (DECT), which could reduce the range uncertainties currently associated with proton radiotherapy. A theoretical basis material decomposition into water and iodine material densities was performed for various pediatric and adult human reference tissues, as well as other non-tissue materials, by minimizing the root-mean-square relative attenuation error in the energy interval from 40 to 140keV. A model (here called MD-SPR) mapping predicted MDs to theoretically calculated reference SPRs was created with locally weighted scatterplot smoothing (LOWESS) data-fitting. The goodness of fit of the MD-SPR model was evaluated for the included reference tissues. MD images of two electron density phantoms, combined to form a head- and an abdomen-sized phantom setup, were acquired with a clinical projection-based fast-kV switching DECT scanner. The MD images were compared to the theoretically predicted MDs of the tissue surrogates and other non-tissue materials in the phantoms, as well as used for input to the MD-SPR model for generation of SPR images. The SPR images were subsequently compared to theoretical reference SPRs of the materials in the phantoms, as well as to SPR images from a commercial algorithm (DirectSPR, Siemens Healthineers, Forchheim, Germany) using image-based consecutive scan DECT for the same phantom setups. The predicted SPRs of the tissue surrogates were similar for MD-SPR and DirectSPR, where the adipose and bone tissue surrogates were within 1% difference to the reference SPRs, while other non-adipose soft tissue surrogates (breast, brain, liver, muscle) were all underestimated by between -0.7% and -1.8%. The SPRs of the non-tissue materials (polymethyl methacrylate (PMMA), polyether ether ketone (PEEK), graphite and Teflon) were within 2.8% for MD-SPR images, compared to 6.8% for DirectSPR. The MD-SPR model performed similar compared to other published methods for the human reference tissues. The SPR prediction for tissue surrogates was similar to DirectSPR and showed potential to improve SPR prediction for non-tissue materials.

Nonlinearly scaled prior image-controlled frequency split for high-frequency metal artifact reduction in computed tomography.

This paper introduces a new approach for the dedicated reduction of high-frequency metal artifacts, which applies a nonlinear scaling (NLS) transfer function on the high-frequency projection domain to reduce artifacts, while preserving edge information and anatomic detail by incorporating prior image information. An NLS function is applied to suppress high-frequency streak artifacts, but to restrict the correction to metal projections only, scaling is performed in the sinogram domain. Anatomic information should be preserved and is excluded from scaling by incorporating a prior image from tissue classification. The corrected high-frequency sinogram is reconstructed and combined with the low-frequency component of a normalized metal artifact reduction (NMAR) image. Scans of different anthropomorphic phantoms were acquired (unilateral hip, bilateral hip, dental implants, and embolization coil). Multiple regions of interest (ROIs) were drawn around the metal implants and hounsfield unit (HU) deviations were analyzed. Clinical data sets including single image slices of dental fillings, a bilateral hip implant, spinal fixation screws, and an aneurysm coil were reconstructed and assessed. The prior image-controlled NLS can remove streak artifacts while preserving anatomic detail within the bone and soft tissue. The qualitative analysis of clinical cases showed a tremendous enhancement within dental fillings and neuro coils, and a significant enhancement within spinal screws or hip implants. The phantom scan measurements support this observation. In all phantom setups, the NLS-corrected result showed lowest HU derivation and the best visualization of thedata. The prior image-controlled NLS provides a method to reduce high-frequency streaks in metal-corrupted computed tomography (CT)data.

Validation and practical use of Plan2Heat hyperthermia treatment planning for capacitive heating

Background Capacitive devices are used for hyperthermia delivery, initially mainly in Asia, but nowadays also increasingly in Europe. Treatment planning can be very useful to determine the most effective patient-specific treatment set-up. This paper provides a validation of GPU-based simulations using Plan2Heat for capacitive hyperthermia devices. Methods Validation was first performed by comparing simulations with an analytical solution for a spherical object placed inside a uniform electric field. Resolution was 5, 2.5 or 1 mm. Next, simulations for homogeneous and inhomogeneous phantom setups were performed for Thermotron RF8 and Celsius TCS capacitive heating devices at 2.5 mm resolution. Also different combinations of electrode sizes were evaluated. Normalized SAR profiles were compared to phantom measurements from the literature. Possible clinical use of treatment planning was demonstrated for an anal cancer patient, evaluating different treatment set-ups in prone and supine position. Results Numerical and analytical solutions showed excellent agreement. At the center of the sphere, the error was 5.1%, 2.9% and 0.2% for a resolution of 5, 2.5 and 1 mm, respectively. Comparison of measurements and simulations for both Thermotron RF8 and Celsius TCS showed very good agreement within 5% for all phantom set-ups. Simulations were capable of accurately predicting the penetration depth; a very relevant parameter for clinical application. The patient case illustrated that planning can provide insight by comparing effectiveness of different treatment strategies. Conclusion Plan2Heat can rapidly and accurately predict heating patterns generated by capacitive devices. Thus, Plan2Heat is suitable for patient-specific treatment planning for capacitive hyperthermia.

Application of a novel diamond detector for commissioning of FLASH radiotherapy electron beams.

PurposeA diamond detector prototype was recently proposed by Marinelli et al. (Medical Physics 2022, https://doi.org/10.1002/mp.15473) for applications in ultrahigh‐dose‐per‐pulse (UH‐DPP) and ultrahigh‐dose‐rate (UH‐DR) beams, as used in FLASH radiotherapy (FLASH‐RT). In the present study, such so‐called flashDiamond (fD) was investigated from the dosimetric point of view, under pulsed electron beam irradiation. It was then used for the commissioning of an ElectronFlash linac (SIT S.p.A., Italy) both in conventional and UH‐DPP modalities.MethodsDetector calibration was performed in reference conditions, under 60Co and electron beam irradiation. Its response linearity was investigated in UH‐DPP conditions. For this purpose, the DPP was varied in the 1.2–11.9 Gy range, by changing either the beam applicator or the pulse duration from 1 to 4 μs. Dosimetric validation of the fD detector prototype was then performed in conventional modality, by measuring percentage depth dose (PDD) curves, beam profiles, and output factors (OFs). All such measurements were carried out in a motorized water phantom. The obtained results were compared with the ones from commercially available dosimeters, namely, a microDiamond, an Advanced Markus ionization chamber, a silicon diode detector, and EBT‐XD GAFchromic films. Finally, the fD detector was used to fully characterize the 7 and 9 MeV UH‐DPP electron beams delivered by the ElectronFlash linac. In particular, PDDs, beam profiles, and OFs were measured, for both energies and all the applicators, and compared with the ones from EBT‐XD films irradiated in the same experimental conditions.ResultsThe fD calibration coefficient resulted to be independent from the investigated beam qualities. The detector response was found to be linear in the whole investigated DPP range. A very good agreement was observed among PDDs, beam profiles, and OFs measured by the fD prototype and reference detectors, both in conventional and UH‐DPP irradiation modalities.ConclusionsThe fD detector prototype was validated from the dosimetric point of view against several commercial dosimeters in conventional beams. It was proved to be suitable in UH‐DPP and UH‐DR conditions, for which no other commercial real‐time active detector is available to date. It was shown to be a very useful tool to perform fast and reproducible beam characterizations in standard clinical motorized water phantom setups. All of the previously mentioned demonstrate the suitability of the proposed detector for the commissioning of UH‐DR linac beams for preclinical FLASH‐RT applications.

Experimental and computational evaluation of capacitive hyperthermia

Objective Hyperthermia as an enhancer of radio- and/or chemotherapy has been confirmed by various trials. Quite a few positive randomized trials have been carried out with capacitive hyperthermia systems (CHS), even though specific absorption rates (SAR) in deep regions are known to be inferior to the established annular-phased array techniques. Due to a lack of systematic SAR measurements for current capacitive technology, we performed phantom measurements in combination with simulation studies. Materials and Methods According to the current guidelines, homogeneous and inhomogeneous agarose phantoms were manufactured for the commercial CHS Celsius42. Temperature/time curves were registered, and specific absorption rate (SAR) profiles and distributions were derived using the temperature gradient method. We implemented models for electrodes and phantom setups for simulation studies using Sim4Life. Results For a standard total power of 200 W, we measured effective SAR until depths of 6–8 cm in a homogeneous phantom, which indicates fair heating conditions for tumor diseases in superficial and intermediate depths. A fat layer of 1 cm strongly weakens the SAR, but 10–20 W/kg are still achieved in intermediate to deep regions (2–10 cm). In the phantom setup with integrated bone, we measured low SAR of 5–10 W/kg in the cancellous bone. Our simulations could fairly describe the measured SAR distributions, but predict tendentially higher SAR than measured. Additional simulations suggest that we would achieve higher SAR with vital fatty tissue and bone metastases in clinical situations. Conclusion Capacitive systems are suitable to heat superficial and medium-deep tumors as well as some bone metastases, and CHS application is feasible for a specific class of patients with pelvic and abdominal tumors. These findings are consistent with positive clinical studies.

Characterization of an advanced cone beam CT (CBCT) reconstruction algorithm used for dose calculation on Varian Halcyon linear accelerators

In this study, the performance of a new iterative reconstruction algorithm, the pre-clinical AcurosXB iCBCT algorithm, has been characterized on Varian Halcyon linear accelerators with respect to the potential of radiotherapy dose calculations on CBCT images. The study utilized various phantom setups to verify the accuracy of the pre-clinical algorithm under different scatter conditions and compared dose calculations performed on CBCT images reconstructed with the pre-clinical algorithm to those performed on typical planning CT images. The results indicated that despite showing improvements compared to the existing iCBCT protocol, certain restrictions should be introduced when the pre-clinical AcurosXB iCBCT algorithm was used for dose calculations. Changes in the scatter condition exhibited a larger effect on CBCTs than on planning CTs. Therefore, users should be careful in offsetting the patient and positioning the patient’s arms if the resultant images will be used for dose calculations. In addition, protocols with different kV settings should be approached with caution, where 100 kV protocols should only be used to scan the head and neck area, while the rest of the body should be scanned with the 125 kV and 140 kV protocols. When the patient is set up properly and the appropriate energy is selected for the anatomical area, the uncertainty of using the novel AcurosXB iCBCT algorithm for treatment planning dose calculation is within ±2.0%.

Experimental comparison of clinically used ion beams for imaging applications using a range telescope

In particle therapy, the x-ray based treatment planning converting photon attenuation values to relative stopping power ratio (RSP) introduces clinically relevant range uncertainties. Recently, novel imaging technologies using transmission ion beams have been investigated to directly assess the water equivalent thickness (WET) of tissue, showing improved accuracy in RSP reconstruction, while potentially reducing the imaging dose. Due to their greater availability, protons have been mostly used for ion imaging. To this end, in this work, the influence of three ion species (protons, helium and carbon ions) on the image quality of radiographic WET retrieval has been explored with a dedicated experimental setup and compared to Monte Carlo (MC) simulations. Three phantom setups with different tissue interfaces and features have been irradiated with clinically validated proton, helium and carbon ion pencil beams under comparable imaging dose and beam settings at the Heidelberg Ion-Beam Therapy Center. Ion radiographies (iRADs) were acquired with an integration mode detector, that functions as a range telescope with 61 parallel plate ionization chambers. For comparison, experiments were reproduced in-silico with FLUKA MC simulations. Carbon ions provide iRADs with highest image quality in terms of normalized root mean square error, followed by helium ions and protons. All ions show similar capabilities of resolving WET for the considered phantoms, as shown by the similar average relative error < 3%. Besides for the slab phantom, MC simulations yielded better results than the experiment, indicating potential improvement of the experimental setup. Our results showed that the ability to resolve the WET is similar for all particles, intrinsically limited by the granularity of the detector system. While carbon ions are best suited for acquiring iRADs with the investigated integration mode detector, helium ions are put forward as a less technical challenging alternative.

An integrated quality assurance phantom for frameless single-isocenter multitarget stereotactic radiosurgery

Brain stereotactic radiosurgery (SRS) treatments require multiple quality assurance (QA) procedures to ensure accurate and precise treatment delivery. As single-isocenter multitarget SRS treatments become more popular, the quantification of off-axis accuracy of the linear accelerator is crucial. In this study, a novel brain SRS integrated phantom was developed and validated to enable SRS QA with a single phantom to facilitate implementation of a frameless single-isocenter, multitarget SRS program. This phantom combines the independent verification of each positioning system, the Winston-Lutz, off-axis accuracy evaluation (i.e. off-axis Winston-Lutz), and the dosimetric accuracy utilizing both point dose measurements as well as film measurement, without moving the phantom. A novel 3D printed phantom, coined OneIso, was designed with a movable insert which can switch between the Winston-Lutz test target and dose measurement without moving the phantom itself. For dose verification, ten brain SRS clinical treatment plans with 10 MV flattening-filter-free beams were delivered on a Varian TrueBeam with a high-definition multileaf collimator (HD-MLC). Radiochromic film and pinpoint ion chamber comparison measurements were made between the OneIso and solid water (SW) phantom setups. For the off-axis Winston-Lutz measurements, a row of off-axis ball bearings (BBs) was integrated into the OneIso. To quantify the spatial accuracy versus distance from the isocenter, two-dimensional displacements were calculated between the planned and delivered BB locations relative to their respective MLC defined field border. OneIso and the SW phantoms agree within 1%, for both film and point dose measurements. OneIso identified a reduction in spatial accuracy further away from the isocenter. Differences increased as distance from the isocenter increased, exceeding recommended SRS accuracy tolerances at 7 cm away from the isocenter. OneIso provides a streamlined, single-setup workflow for single-isocenter multitarget frameless linac-based SRS QA. Additionally, with the ability to quantify off-axis spatial discrepancies, we can determine limitations on the maximum distance between targets to ensure a single-isocenter multitarget SRS program meets recommended guidelines.

Technical Note: An anthropomorphic phantom with implanted neurostimulator for investigation of MRI safety.

The objective of this work was to design and construct an improved anthropomorphic phantom for use in studying magnetic resonance imaging (MRI) radiofrequency (RF) safety at 3T related to deep brain stimulation (DBS), and especially the role of DBS lead trajectories. Based on a computer-aided design including reasonable representation of human features, the phantom was fabricated by three-dimensional (3D) printing and then fully assembled with a human skull, a commercial DBS device implanted using the surgical standard at our institution, and fiber-optic temperature sensors embedded in two tissue mimicking solutions (e.g., the heterogeneous setup). Preliminary MRI safety experiments were conducted using turbo spin-echo (TSE) imaging with the device powered on and powered off. These results were then compared to analogous results for a homogeneous phantom setup that filled the structure with a standard body average solution. Both phantom setups produced temperature increases of ~1.0°C, with a maximum increase of 1.1±0.2°C recorded during imaging of the heterogeneous phantom setup. The preliminary experimental results suggest that improved phantom structures capable of replicating actual DBS lead trajectories may be advisable when conducting DBS-related MRI safety studies. An anthropomorphic phantom was constructed with promising initial results indicating different DBS lead trajectories and phantom setups may impact temperature elevations along an implanted DBS lead. Although additional work will be necessary to validate its efficacy over conventional phantoms, the anthropomorphic phantom can likely be used in the future to assess different procedures for DBS lead placement, the RF power deposition of MRI protocols applicable to DBS patients, and to validate novel methods to reduce localized heating effects associated with DBS devices, such as parallel RF transmission.

Geometric inaccuracy and co-registration errors for CT, DynaCT and MRI images used in robotic stereotactic radiosurgery treatment planning

PurposeTo measure the combined errors due to geometric inaccuracy and image co-registration on secondary images (dynamic CT angiography (dCTA), 3D DynaCT angiography (DynaCTA), and magnetic resonance images (MRI)) that are routinely used to aid in target delineation and planning for stereotactic radiosurgery (SRS). MethodsThree phantoms (one commercial and two in-house built) and two different analysis approaches (commercial and MATLAB based) were used to quantify the magnitude of geometric image distortion and co-registration errors for different imaging modalities within CyberKnife’s MultiPlan treatment planning software. For each phantom, the combined errors were reported as a mean target registration error (TRE). The mean TRE’s for different intramodality imaging parameters (e.g., mAs, kVp, and phantom set-ups) and for dCTA, DynaCTA, and MRI systems were measured. ResultsOnly X-ray based imaging can be performed with the commercial phantom, and the mean TRE ± standard deviation values were large compared to the in-house analysis using MATLAB. With the 3D printed phantom, even drastic changes in treatment planning CT imaging protocols did not greatly influence the mean TRE (<0.5 mm for a 1 mm slice thickness CT). For all imaging modalities, the largest mean TRE was found on DynaCT, followed by T2-weighted MR images (albeit all <1 mm). ConclusionsThe user may overestimate the mean TRE if the commercial phantom and MultiPlan were used solely. The 3D printed phantom design is a sensitive and suitable quality assurance tool for measuring 3D geometric inaccuracy and co-registration errors across all imaging modalities.

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.

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

Effects of metal implants and a metal artifact reduction tool on calculation accuracy of AAA and Acuros XB algorithms in small fields.

In this study, the dosimetric accuracy of analytical anisotropic algorithm (AAA) and Acuros XB (AXB) dose calculation algorithms (Varian Medical Systems, Palo Alto, CA) was investigated for small radiation fields incident on phantoms of various metals that include stainless steel grade 316L (SS316L) and titanium alloy grade 5 (Ti5) implants. In addition, the effects of using metal artifact reduction for orthopedic implants (O-MAR, Philips Healthcare, Cleveland, OH) were evaluated. The evaluations of AAA and AXB were performed by comparing the crossline profiles calculated by AAA and AXB with GafChromicTM EBT3 film measurements at the phantom-implant interfaces and in close vicinity of implant materials for small field sizes (1×1cm2 , 2×2cm2 , 3×3cm2 , and 4×4cm2 ) of a 6MV flattening filter free photon beam. O-MAR corrected and uncorrected (UC) computed tomography (CT) images were used for dose calculations. The values of average and standard deviations (SD) of Hounsfield unit (HU) for selected regions of each case were evaluated. The differences in average dose percentages in defined regions were calculated to quantify the relative dosimetric changes between doses calculated on UC and O-MAR corrected CT images. Compared to UC images, the values of SD were reduced, and the average HU became closer to its reference value in the O-MAR images. There was some discrepancy in average dose percentage differences between calculations using UC and O-MAR images at 1cm above the SS316L implant (average dose percentage differences were AXB/UC=5.9% and AXB/O-MAR=-1.2%; AAA/UC=2.2%, and AAA/O-MAR=-0.8%). Neither AAA nor AXB algorithms predict increase in dose at upper phantom-implant interface (4.9%, 9.9%. 13.5%, and 13.8% for the fields from 1×1cm2 to 4×4cm2 , respectively). At the side of the SS316L implant (where dark streak artifacts exist), dose difference averages were estimated as-1.1% and 22.3% when AXB/O-MAR and AXB/UC calculations are compared with EBT3 measurements, respectively. Dose predictions at 1cm below the SS316L implant were underestimated by AXB/O-MAR (average -0.5%) and AXB/UC (average 2.0%). The O-MAR tool was shown to have a favorable dosimetric effect or no effect on the calculations in the upper proximity of the implant materials. The dose differences between EBT3 film measurements and calculations at upper phantom-implant interfaces were smaller when they were calculated using O-MAR images. However, the dose differences increased when O-MAR corrected images were used for AAA calculations at lower phantom-implant interfaces. Use of O-MAR enabled closer agreement for the AXB algorithm, especially in the dark streak artifact regions. The O-MAR algorithm should be used when the dose is calculated with the AXB algorithm in cases of patients with the metal implants. The estimations using AAA and AXB algorithms, in phantom setups, with Ti5 implant material were found to be closer to the EBT3 film measurements, when compared with the same estimations using SS316L implant material.

Two-dimensional EPID dosimetry for an MR-linac: Proof of concept.

At our institute, in vivo patient dose distributions are reconstructed for all treatments delivered using conventional linacs from electronic portal imaging device (EPID) transit images acquired during treatment using a simple back-projection model. Currently, the clinical implementation of MRI-guided radiotherapy systems, which aims for online and real-time adaptation of the treatment plan, is progressing. In our department, the MR-linac (Unity, Elekta AB, Stockholm, Sweden) is now in clinical use. The aim of this work is to demonstrate the feasibility of two-dimensional (2D) EPID dosimetric verification for the magnetic resonance (MR)-linac by comparing back-projected EPID doses to ionization chamber (IC) array dose distributions. Our conventional back-projection algorithm was adapted for the MR-linac. The most important changes involve modeling of the attenuation by and scatter from the cryostat. The commissioning process involved the acquisition of square field EPID measurements using various phantom setups (varying SSD, phantom thickness, and field size). Commissioning models were created for gantry 0°, 90°, and 180° and verified by comparing EPID-reconstructed 2D dose distributions to measurements made with the OCTAVIUS 1500 IC array (PTW, Freiburg, Germany) for two prostate and one rectum IMRT plans (25 beams total). The average of the γ parameters (y-mean and y-pass rate) and the dose difference at a reference point were reported. Due to their construction, the attenuation of couch, bridge, and cryostat shows a much stronger dependence on gantry angle in the MR-linac compared to conventional linacs. We present a method to correct for these effects. This method is validated by dose reconstruction of the 25 intensity-modulated radiation therapy beams recorded at a certain gantry angle using the model of another gantry angle, combined with the correction method. For dose verification performed at a gantry angle identical to the commissioned model, the average y-mean and y-pass rate values (3% global dose, 2mm, 10% isodose) were 0.37±0.07 and 98.1, 95% CI [98.1±2.4], respectively. The average dose difference at the reference point was -0.5%±1.8%. Verification at gantry angles different from the commissioned model (i.e., using the gantry angle dependent correction) reported 0.39±0.08 and 97.6, 95% CI [96.9, 98.3] average y-mean and y-pass rate values. The average dose difference at the reference point was -0.1%±1.8%. The EPID dosimetry back-projection model was successfully adapted for the MR-linac at gantry 0°, 90°, and 180°, accounting for the presence of the MRI housing between phantom (or patient) and the EPID. A method to account for the gantry angle dependence was also tested reporting similar results.

Impact of backscatter material thickness on the depth dose of orthovoltage irradiators for radiobiology research

The orthovoltage x-ray energy frequently used in radiation research is prone to dosimetry errors due to insufficient backscatter conditions. In many radiobiology studies, especially for cell irradiations, precise dose calculation algorithms such as Convolution-Superposition or Monte Carlo are impractical and as such, less accurate hand calculation methods are used for dose estimation. These dose estimation methods typically assume full backscatter conditions. The purpose of this study is to demonstrate the magnitude of the dose error that results from insufficient backscatter, and to provide lookup tables to account this issue. The beam spectra of several widely used commercial systems (XRAD-225, XRAD-320, SARRP) were used in Monte Carlo (MC) simulations on a series of phantom setups to investigate the impact of varying backscatter conditions on dosimetry. The depth dose curves for different field sizes, water phantom thicknesses and beam qualities were generated. In addition, depth dependent backscatter factors for different field sizes and different beam qualities were calculated. It is demonstrated that as much as a 50% dose difference exists for different backscatter conditions at the beam qualities studied. The choice of cell dish size as well as other changes in the experiment setup can have more than 10% impact on the dose. The impact of backscatter is reduced with a decrease in field size. Further, the thickness needed to provide full backscatter can be approximated as being equal to the field size. It is imperative to ensure full backscatter conditions during system and dosimeter calibration, or to use the look-up table provided in this study.

Part I: Out-of-field dose mapping for 6X and 6X-flattening-filter-free beams on the TrueBeam for extended distances.

With increasing cancer treatment success rates, many patients go on to live long, productive lives following recovery. Therefore, minimizing potential side effects due to dose outside the treated field is becoming a significant consideration in radiation therapy. With many potential treatment configurations available, it is important to quantify how out-of-field dose varies with common variables such as distance from isocenter, couch angle, jaw size, and flattening-filter setting. The accurate quantification of out-of-field dose at extended distances could also benefit researchers and detector developers. While data exist for out-of-field dose from older linear accelerator (Linac) models, the phenomenon has not been described for the latest generation of machines, such as the Varian TrueBeam. The purpose of this study was to comprehensively quantify out-of-field dose for the Varian TrueBeam Linac low energy photons in a wide range of positions and treatment geometries. Out-of-field doses were measured using two phantom setups: (a) A large volume ion chamber with a buildup sleeve to quantify head leakage and collimator scatter background dose; and (b) A farmer ion chamber in solid water to incorporate phantom scatter in addition to collimator scatter, and head leakage background dose. In both cases, the ion chamber was positioned with its length along the slowly varying transverse direction (perpendicular to the radial from isocenter). Doses were measured for four symmetric jaw settings (2×2cm2 , 4×4cm2 , 10×10cm2 , and 20×20cm2 ) for a range of distances from the isocenter (0-100cm). The angular dependence of the out-of-field dose was measured using four different angles: 0°, 45°, 90°, and 135° with respect to the in-plane direction. All measurements were performed for both 6X and 6X-flattening-filter-free (FFF) beams. The lowest out-of-field doses were observed at 60cm away from isocenter in both in-plane and cross-plane directions for fields smaller than 10×10cm2 . Out-of-field dose decreased with decreasing jaw size (a factor of 4.7 for 6X-FFF and a factor of 3.1 for 6X going from 20×20cm2 to 2×2cm2 at 60cm from isocenter in the in-plane direction). The 6X-FFF beam produced out-of-field doses as low as 64% of the 6X beam. This study presents a comprehensive description of 6X and 6X-FFF out-of-field doses on a Varian TrueBeam Linac including measurements at a range of positions, angles, and jaw settings and with and without phantom scatter.

Structure delineation in the presence of metal - A comparative phantom study using single and dual-energy computed tomography with and without metal artefact reduction.

Background and purposeMetal artefacts in computed tomography (CT) images impairs structure delineation. These artefacts can potentially be reduced with dual-energy CT (DECT) with or without using metal artefact reduction (MAR). The purpose was to investigate how structure delineation in DECT with or without MAR and single-energy CT (SECT) images were affected by metals.Materials and methodsA phantom with known irregular structures was developed. Reference structures were determined from a low-noise scan without metal. Bilateral hip prostheses were simulated with steel or titanium inserts. The phantom was scanned with SECT and fast-kV switching DECT with optional MAR. Four radiation oncologists delineated the structures in two phantom set-ups. Delineated structures were evaluated with Dice similarity coefficient (DSC) and Hausdorff distance relative to the reference structures.ResultsWith titanium inserts, more structures were detected for non-MAR DECT compared to SECT while the same or less were detected with steel inserts. MAR improved delineation in DECT images. For steel inserts, three structures in the region of artefacts, were delineated by at least two oncologists with MAR-DECT compared to none with non-MAR DECT or SECT. The highest values of DSC for MAR-DECT were 0.69, 0.81 and 0.77 for those structures.ConclusionsDelineation was improved with non-MAR DECT compared to SECT, especially for titanium inserts. A larger improvement was seen with the use of MAR for both steel and titanium inserts. The improvement was dependent on the location of the structure relative to the inserts, and the structure contrast relative to the background.

Navigating surgical fluorescence cameras using near-infrared optical tracking.

Fluorescence guidance facilitates real-time intraoperative visualization of the tissue of interest. However, due to attenuation, the application of fluorescence guidance is restricted to superficial lesions. To overcome this shortcoming, we have previously applied three-dimensional surgical navigation to position the fluorescence camera in reach of the superficial fluorescent signal. Unfortunately, in open surgery, the near-infrared (NIR) optical tracking system (OTS) used for navigation also induced an interference during NIR fluorescence imaging. In an attempt to support future implementation of navigated fluorescence cameras, different aspects of this interference were characterized and solutions were sought after. Two commercial fluorescence cameras for open surgery were studied in (surgical) phantom and human tissue setups using two different NIR OTSs and one OTS simulating light-emitting diode setup. Following the outcome of these measurements, OTS settings were optimized. Measurements indicated the OTS interference was caused by: (1)spectral overlap between the OTS light and camera, (2)OTS light intensity, (3)OTS duty cycle, (4)OTS frequency, (5)fluorescence camera frequency, and (6)fluorescence camera sensitivity. By optimizing points 2 to 4, navigation of fluorescence cameras during open surgery could be facilitated. Optimization of the OTS and camera compatibility can be used to support navigated fluorescence guidance concepts.