Use of a megavoltage electronic portal imaging device to identify prosthetic materials

Characterization and evaluation methods of fused deposition modeling and stereolithography additive manufacturing for clinical linear accelerator photon and electron radiotherapy applications

PURPOSE

To propose comprehensive characterization methods of additive manufacturing (AM) materials for MV photon and MeV electron radiotherapy.

METHODOLOGY

This study investigated 15 AM materials using CT machines. Geometrical accuracy, tissue-equivalence, uniformity, and fabrication parameters were considered. Selected soft tissue equivalent filaments were used to fabricate slab phantoms and compared with water equivalent RW3 phantom by delivering planar 6 & 10 MV photons and 6, 9, 12, 15, & 18 MeV electrons. Finally, a 3D printed CT-Electron Density characterization phantom was fabricated.

RESULTS

Materials used to print test objects can simulate tissues from adipose (relative electron density, ρe=0.72) up to near inner bone-equivalent (ρe=1.08). Lower densities such as breast and lung can be simulated using infills from 90 % down to 30 %, respectively. The gyroid infill pattern shows the lowest CT number variation and is recommended for low infill percentage printing. CT number uniformity can be observed from 40 % up to 100 % infill, while printing orientation does not significantly affect the CT number. The measured doses using the 3D printed phantoms show to have good agreement with TPS calculated dose for photon (< 1 % difference) and electron (< 5 % difference). Varying the printed slab thicknesses shows very similar response (< 3 % difference) compared with RW3 slabs except for 6 MeV electrons. Lastly, the fabricated CT-ED phantom generally matches the lung- up to the soft tissue- equivalence.

CONCLUSION

The proposed methods give the outline for characterization of AM materials as tissue-equivalent substitute. Printing parameters affect the radiological quality of 3D-printed object.

Characterisation of 3D-printable thermoplastics to be used as tissue-equivalent materials in photon and proton beam radiotherapy end-to-end quality assurance devices

Objective.To investigate the potential of 3D-printable thermoplastics as tissue-equivalent materials to be used in multimodal radiotherapy end-to-end quality assurance (QA) devices.Approach.Six thermoplastics were investigated: Polylactic Acid (PLA), Acrylonitrile Butadiene Styrene (ABS), Polyethylene Terephthalate Glycol (PETG), Polymethyl Methacrylate (PMMA), High Impact Polystyrene (HIPS) and StoneFil. Measurements of mass density (ρ), Relative Electron Density (RED), in a nominal 6 MV photon beam, and Relative Stopping Power (RSP), in a 210 MeV proton pencil-beam, were performed. Average Hounsfield Units (HU) were derived from CTs acquired with two independent scanners. The calibration curves of both scanners were used to predict averageρ,RED and RSP values and compared against the experimental data. Finally, measured data ofρ,RED and RSP was compared against theoretical values estimated for the thermoplastic materials and biological tissues.Main results.Overall, goodρand RSP CT predictions were made; only PMMA and PETG showed differences >5%. The differences between experimental and CT predicted RED values were also <5% for PLA, ABS, PETG and PMMA; for HIPS and StoneFil higher differences were found (6.94% and 9.42/15.34%, respectively). Small HU variations were obtained in the CTs for all materials indicating good uniform density distribution in the samples production. ABS, PLA, PETG and PMMA showed potential equivalency for a variety of soft tissues (adipose tissue, skeletal muscle, brain and lung tissues, differences within 0.19%-8.35% for all properties). StoneFil was the closest substitute to bone, but differences were >10%. Theoretical calculations of all properties agreed with experimental values within 5% difference for most thermoplastics.Significance.Several 3D-printed thermoplastics were promising tissue-equivalent materials to be used in devices for end-to-end multimodal radiotherapy QA and may not require corrections in treatment planning systems' dose calculations. Theoretical calculations showed promise in identifying thermoplastics matching target biological tissues before experiments are performed.

Evaluation of radiation detectors for the determination of field output factors in Leksell Gamma Knife dosimetry using 3D printed phantom inserts

The Leksell Gamma Knife® Perfexion™ and Icon™ have a unique geometry, containing 192 60Co sources with collimation for field sizes of 4 mm, 8 mm, and 16 mm. 4 mm and 8 mm collimated fields lack lateral charged particle equilibrium, so accurate field output factors are essential. This study performs field output factor measurements for the microDiamond, microSilicon, and RAZOR™ Nano detectors. 3D printed inserts for the spherical Solid Water® Phantom were fabricated for microDiamond detector, the microSilicon unshielded diode and the RAZOR™ Nano micro-ionisation chamber. Detectors were moved iteratively to identify the peak detector signal for each collimator, representing the effective point of measurement of the chamber. In addition, field output correction factors were calculated for each detector relative to vendor supplied Monte Carlo simulated field output factors and field output factors measured with a W2 scintillator. All field output factors where within 1.1 % for the 4 mm collimator and within 2.3 % for the 8 mm collimator. The 3D printed phantom inserts were suitable for routine measurements if the user identifies the effective point of measurement, and ensures a reproducible setup by marking the rotational alignment of the cylindrical print. Measurements with the microDiamond and microSilicon can be performed faster compared to the RAZOR™ Nano due to differences in the signal to noise ratio. All detectors are suitable for field output factor measurements for the Leksell Gamma Knife® Perfexion™ and Icon™.

Regression fitting megavoltage depth dose curves to determine material relative electron density in radiotherapy

The relative electron density (RED) parameter is ubiquitous throughout radiotherapy for clinical dosimetry and treatment planning purposes as it provides a more accurate description of the relevant radiological properties over mass density alone. RED is in practice determined indirectly from calibrated CT Hounsfield Units (HU). While CT images provide useful 3D information, the spectral differences between CT and clinical LINAC beams may impact the validity of the CT-ED calibration, especially in the context of novel tissue-mimicking materials where deviations from biologically typical atomic number to atomic weight ratios 〈Z/A〉 occur and/or high-Z materials are present. A theoretical basis for determining material properties directly in a clinical beam spectrum via an electron-density equivalent pathlength (eEPL) method has been previously established. An experimental implementation of this approach is introduced whereby material-specific measured percentage depth dose curves (PDDs) are regressed to a PDD measured in a reference material (water), providing an inference of 〈Z/A〉, which when combined with the physical density provides a determination of RED. This method is validated over a range of tissue-mimicking materials and compared against the standard CT output, as well as compositional information obtained from the manufacturer's specifications. The measured PDD regression method shows consistent results against both manufacturer-provided and CT-derived values between 0.9 and 1.15 RED. Outside of this soft-tissue range a trend was observed whereby the 〈Z/A〉 determined becomes unrealistic indicating the method is no longer reporting RED alone and the assumptions around the eEPL model are constrained. Within the soft-tissue RED range of validity, the regression method provides a practical and robust characterisation for unknown materials in the clinical setting and may be used to improve on the CT derived RED where high Z material components are suspected.

A markerless beam's eye view tumor tracking algorithm based on unsupervised deformable registration learning framework of VoxelMorph in medical image with partial data

PURPOSE

To propose an unsupervised deformable registration learning framework-based markerless beam's eye view (BEV) tumor tracking algorithm for the inferior quality megavolt (MV) images with occlusion and deformation.

METHODS

Quality assurance (QA) plans for thorax phantom were delivered to the linear accelerator with artificially treatment offsets. Electronic portal imaging device (EPID) images (682 in total) and corresponding digitally reconstructed radiograph (DRR) were gathered as the moving and fixed image pairs, which were randomly divided into training and testing set in a ratio of 0.7:0.3 to train a non-rigid registration model with Voxelmorph. Simultaneously, 533 pairs of patient images from 21 lung tumor plans were acquired for tumor tracking investigation to offer quantifiable tumor motion data. Tracking error and image similarity measures were employed to evaluate the algorithm's accuracy.

RESULTS

The tracking algorithm can handle image registration with non-rigid deformation and losses ranging from 10 % to 80 %. The tracking errors in the phantom study were below 3 mm in about 86.8 % of cases, and below 2 mm in about 80 % of cases. The normalized mutual information (NMI) changes from 1.182 ± 0.024 to 1.198 ± 0.024 (p < 0.005). The patient target had an average translation of 3.784 mm and a maximum comprehensive displacement value of 7.455 mm. NMI of patient images changes from 1.209 ± 0.027 to 1.217 ± 0.026 (p < 0.005), indicating the presence of non-negligible non-rigid deformation.

CONCLUSIONS

The study provides a robust markerless tumor tracking algorithm for multi-modal, partial data and inferior quality image processing.

3D-printed patient-specific pelvis phantom for dosimetry measurements for prostate stereotactic radiotherapy with dominant intraprostatic lesion boost

AIM

Developing and assessing the feasibility of using a three-dimensional (3D) printed patient-specific anthropomorphic pelvis phantom for dose calculation and verification for stereotactic ablative radiation therapy (SABR) with dose escalation to the dominant intraprostatic lesions.

MATERIAL AND METHODS

A 3D-printed pelvis phantom, including bone-mimicking material, was fabricated based on the computed tomography (CT) images of a prostate cancer patient. To compare the extent to which patient and phantom body and bones overlapped, the similarity Dice coefficient was calculated. Modular cylindrical inserts were created to encapsulate radiochromic films and ionization chamber for absolute dosimetry measurements at the location of prostate and at the boost region. Gamma analysis evaluation with 2%/2mm criteria was performed to compare treatment planning system calculations and measured dose when delivering a 10 flattening filter free (FFF) SABR plan and a 10FFF boost SABR plan.

RESULTS

Dice coefficients of 0.98 and 0.91 were measured for body and bones, respectively, demonstrating agreement between patient and phantom outlines. For the boost plans the gamma analysis yielded 97.0% of pixels passing 2%/2mm criteria and these results were supported by the chamber average dose difference of 0.47 ± 0.03%. These results were further improved when overriding the bone relative electron density: 97.3% for the 2%/2mm gamma analysis, and 0.05 ± 0.03% for the ionization chamber average dose difference.

CONCLUSIONS

The modular patient-specific 3D-printed pelvis phantom has proven to be a highly attractive and versatile tool to validate prostate SABR boost plans using multiple detectors.

Investigation of the effects of spinal surgical implants on radiotherapy dosimetry: A study of 3D printed phantoms

PURPOSE

The use of three-dimensional (3D) printing to develop custom phantoms for dosimetric studies in radiotherapy is increasing. The process allows production of phantoms designed to evaluated specific geometries, patients, or patient groups with a defining feature. The ability to print bone-equivalent phantoms has, however, proved challenging. The purpose of this work was to 3D print a series of three similar spine phantoms containing no surgical implants, implants made of titanium, and implants made of carbon fiber, for future dosimetric and imaging studies. Phantoms were evaluated for (a) tissue and bone equivalence, (b) geometric accuracy compared to design, and (c) similarity to one another.

METHODS

Sample blocks of PLA, HIPS, and StoneFil PLA-concrete with different infill densities were printed to evaluate tissue and bone equivalence. The samples were used to develop CT to physical (PD) and effective relative electron density (RED_eff ) conversion curves and define the settings for printing the phantoms. CT scans of the printed phantoms were obtained to assess the geometry and densities achieved. Mean distance to agreement (MDA) and DICE coefficient (DSC) values were calculated between contours defining the different materials, obtained from design and like phantom modules. HU values were used to determine PD and RED_eff and subsequently evaluate tissue and bone equivalence.

RESULTS

Sample objects showed linear relationships between HU and both PD and RED_eff for both PLA and StoneFil. The PD and RED_eff of the objects calculated using clinical CT conversion curves were not accurate and custom conversion curves were required. PLA printed with 90% infill density was found to have a PD of 1.11 ± 0.03 g.cm^-3 and RED_eff of 1.04 ± 0.02 and selected for tissue- equivalent phantom elements. StoneFil printed with 100% infill density showed a PD of 1.35 ± 0.03 g.cm^-3 and RED_eff of 1.24 ± 0.04 and was selected for bone-equivalent elements. Upon evaluation of the final phantoms, the PLA elements displayed PD in the range of 1.10 ± 0.03 g.cm^-3 -1.13 ± 0.03 g.cm^-3 and RED_eff in the range of 1.02 ± 0.03-1.06 ± 0.03. The StoneFil elements showed PD in the range of 1.43 ± 0.04 g.cm^-3 -1.46 ± 0.04 g.cm^-3 and RED_eff in the range of 1.31 ± 0.04-1.33 ± 0.04. The PLA phantom elements were shown to have MDA of ≤1.00 mm and DSC of ≥0.95 compared to design, and ≤0.48 mm and ≥0.91 compared like modules. The StoneFil elements displayed MDA values of ≤0.44 mm and DSC of ≥0.98 compared to design and ≤0.43 mm and ≥0.92 compared like modules.

CONCLUSIONS

Phantoms which were radiologically equivalent to tissue and bone were produced with a high level of similarity to design and even higher level of similarity of one another. When used in conjunction with the derived CT to PD or RED_eff conversion curves they are suitable for evaluating the effects of spinal surgical implants of varying material of construction.

Addressing the dosimetric impact of bone cement and vertebroplasty in stereotactic body radiation therapy

PURPOSE

Bone cement used for vertebroplasty can affect the accuracy on the dose calculation of the radiation therapy treatment. In addition the CT values of high density objects themselves can be misrepresented in kVCT images. The aim of our study is then to propose a streamlined approach for estimating the real density of cement implants used in stereotactic body radiation therapy.

METHODS

Several samples of cement were manufactured and irradiated in order to investigate the impact of their composition on the radiation dose. The validity of the CT conversion method for a range of photon energies was investigated, for the studied samples and on six patients. Calculations and measurements were carried out with various overridden densities and dose prediction algorithms (AXB with dose-to-medium reporting or AAA) in order to find the effective density override.

RESULTS

Relative dose differences of several percent were found between the dose measured and calculated downstream of the implant using an ion chamber and TPS or EPID dosimetry. If the correct density is assigned to the implant, calculations can provide clinically acceptable accuracy (gamma criteria of 3%/2 mm). The use of MV imaging significantly favors the attribution of a correct equivalent density to the implants compared to the use of kVCT images.

CONCLUSION

The porosity and relative density of the various studied implants vary significantly. Bone cement density estimations can be characterized using MV imaging or planar in vivo dosimetry, which could help determining whether errors in dose calculations are due to incorrect densities.

Impact of radiopacified bone cement on radiotherapy dose calculation

Background and purpose

Radiopacifiers are introduced to bone cements to provide the appearance of bone in kilovoltage (kV) radiographic images. For higher energy megavoltage (MV) radiotherapy treatment beams, however, these radiopacifiers do not cause a bone-like perturbation of dose. This study therefore aimed to determine the impact of the barium-contrasted plastic-based cement materials on radiotherapy dose calculations.

Materials and methods

The radiological properties of a physical sample of bone cement were characterised by computed tomography (CT) imaging and transmission measurements. Monte Carlo simulations of percentage depth-dose profiles were performed to determine the possible dose error for MV treatment beams. Dose differences were then investigated for clinical volumetric modulated radiotherapy treatment plans, with and without density overrides applied.

Results

Differences of up to 7% were observed at the downstream interface of a 0.6 cm thick bone cement layer, compared to bone. Differences in planning target volume dose-volume metrics varied between −0.5% and 2.0%.

Conclusion

Before planning radiotherapy treatments for patients who have undergone cranioplasty, every effort should be made to identify whether a radiopacified bone cement has been implanted. Density overrides should be applied to minimise dose calculation errors, whenever bone cement is used.

Verification of the Elekta Monaco TPS Monte Carlo in modelling radiation transmission through metals in a water equivalent phantom

Many studies have performed dosimetric studies using various metal implants however these are difficult to translate to other implants of a different geometry or material (Rijken and Colyer, J Appl Clin Med Phys 18:5:301-306, 2017; Ade and du Plessis, J Appl Clin Med Phys 18:5:162-173, 2017; Prabhakar et al. Rep Pract Oncol Radiother 18:209-213, 2013; Ng et al. Rep Pract Oncol Radiother 20:273-277, 2015; Reft et al. Med Phys 30:1162-1182, 2003; Sasaki et al., Nihon Hoshasen Gijutsu Gakkai Zasshi 72(9):735-745, 2016). In this study, the ability of the Monaco Monte Carlo algorithm (Elekta AB, Stockholm, Sweden) to model radiation transport through different types of metals was evaluated. Investigation of the capabilities and limitations of the algorithm is required for the potential use of Monaco for planning radiotherapy treatments when avoidance of metal implants is clinically undesirable. A MapCHECK 2 diode array (Sun Nuclear Corp, Melbourne, USA) and a PTW 30013 Farmer chamber was used to measure the dose at depth, downstream of 1 cm × 5 cm × 5 cm metal blocks of three known compositions; stainless steel, aluminium and MCP96. The setup was imaged using a CT scanner and imported into the Monaco TPS where the beam arrangement was replicated. The density of the metals was overridden using the known electron density of each (IMPAC Medical Systems Inc, Monaco dose calculation technical reference. IMPAC Medical Systems, Sunnydale, CA, 2013). The differences between the dose measured using the ion chamber and calculated using Monaco downstream of the 1 cm metal blocks were respectively: - 1.2%, - 2.2% and 9.5% when irradiated using a 6 MV beam, and - 0.9%, - 1.3% and 14%, when irradiated using a 15 MV beam. This was then repeated using 2 cm and 3 cm of each metal type giving similar results for aluminium and stainless steel and increased discrepancy for MCP96. Discrepancies between treatment planning software and measurements at depth have been shown to give uncertainties between 5 and 23% in previous studies (Rijken and Colyer, J Appl Clin Med Phys 18:5:301-306, 2017; Ade and du Plessis, J Appl Clin Med Phys 18:5:162-173, 2017; Prabhakar et al. Rep Pract Oncol Radiother 18:209-213, 2013; Ng et al. Rep Pract Oncol Radiother 20:273-277, 2015; Reft et al. Med Phys 30:1162-1182, 2003; Sasaki et al., Nihon Hoshasen Gijutsu Gakkai Zasshi 72(9):735-745, 2016). This study uses basic shapes providing results that remove the uncertainties in geometry and can therefore be applied to any shape. This will help determine whether errors in dose calculations are due to the TPS particle transport algorithms or due to other effects, such as inaccurate contouring or incorrect densities. Thus giving the planner an additional degree of freedom in their planning and decision making process.

Web of Science, Scopus, and Google Scholar citation rates: a case study of medical physics and biomedical engineering: what gets cited and what doesn't?

There are often differences in a publication's citation count, depending on the database accessed. Here, aspects of citation counts for medical physics and biomedical engineering papers are studied using papers published in the journal Australasian physical and engineering sciences in medicine. Comparison is made between the Web of Science, Scopus, and Google Scholar. Papers are categorised into subject matter, and citation trends are examined. It is shown that review papers as a group tend to receive more citations on average; however the highest cited individual papers are more likely to be research papers.

Use of electronic portal imaging devices for electron treatment verification

This study aims to help broaden the use of electronic portal imaging devices (EPIDs) for pre-treatment patient positioning verification, from photon-beam radiotherapy to photon- and electron-beam radiotherapy, by proposing and testing a method for acquiring clinically-useful EPID images of patient anatomy using electron beams, with a view to enabling and encouraging further research in this area. EPID images used in this study were acquired using all available beams from a linac configured to deliver electron beams with nominal energies of 6, 9, 12, 16 and 20 MeV, as well as photon beams with nominal energies of 6 and 10 MV. A widely-available heterogeneous, approximately-humanoid, thorax phantom was used, to provide an indication of the contrast and noise produced when imaging different types of tissue with comparatively realistic thicknesses. The acquired images were automatically calibrated, corrected for the effects of variations in the sensitivity of individual photodiodes, using a flood field image. For electron beam imaging, flood field EPID calibration images were acquired with and without the placement of blocks of water-equivalent plastic (with thicknesses approximately equal to the practical range of electrons in the plastic) placed upstream of the EPID, to filter out the primary electron beam, leaving only the bremsstrahlung photon signal. While the electron beam images acquired using a standard (unfiltered) flood field calibration were observed to be noisy and difficult to interpret, the electron beam images acquired using the filtered flood field calibration showed tissues and bony anatomy with levels of contrast and noise that were similar to the contrast and noise levels seen in the clinically acceptable photon beam EPID images. The best electron beam imaging results (highest contrast, signal-to-noise and contrast-to-noise ratios) were achieved when the images were acquired using the higher energy electron beams (16 and 20 MeV) when the EPID was calibrated using an intermediate (12 MeV) electron beam energy. These results demonstrate the feasibility of acquiring clinically-useful EPID images of patient anatomy using electron beams and suggest important avenues for future investigation, thus enabling and encouraging further research in this area. There is manifest potential for the EPID imaging method proposed in this work to lead to the clinical use of electron beam imaging for geometric verification of electron treatments in the future.