Daniel Frauchiger

Optoacoustic imaging using a three-dimensional reconstruction algorithm

A novel three-dimensional (3-D) computational algorithm to reconstruct 3-D optoacoustic images from two-dimensional (2-D) pressure distributions generated at the sample surface is presented. The 2-D pressure distributions were measured as images at different delay times after the excitation laser pulse. The pressure images were captured with a gated CCD camera as the local pressure induces intensity changes of a reflected probe beam at the surface of the irradiated sample. The illumination time was 10 ns and the resolution of the surface pressure image was 20 /spl mu/m. The reconstruction algorithm is based on the decomposition into plane waves. The algorithm was tested in the back projection of simulated pressure transients of three sources, and applied to different biological systems. Furthermore the algorithm was compared with the time of flight back projection algorithm. Optoacoustic images with a depth resolution of 15 /spl mu/m and a lateral resolution of 100 /spl mu/m are presented.

Real-time optoacoustic imaging using a Schlieren transducer

Optoacoustics, which combines the merits and most compelling features of light and sound, is a technique for diagnosing optical tissue properties and for performing noninvasive medical imaging. We present a detection method based on Schlieren optical imaging used to record optoacoustic waves generated with 6 ns long laser pulses. These waves result from irradiation of optically absorbing targets hidden deep inside highly scattering media. Since they are generated from one single shot, Schlieren images contain all the information necessary for the reconstruction of the absorbing structures. The detection method allows on-line imaging with a spatial resolution of about 30 μm with high contrast. The proposed real-time imaging system not only enhances the level of visible detail but it is ideally suited for monitoring dynamic tissue changes.

Monte Carlo implementation, validation, and characterization of a 120 leaf MLC

PURPOSE Recently, the new high definition multileaf collimator (HD120 MLC) was commercialized by Varian Medical Systems providing high resolution in the center section of the treatment field. The aim of this work is to investigate the characteristics of the HD120 MLC using Monte Carlo (MC) methods. METHODS Based on the information of the manufacturer, the HD120 MLC was implemented into the already existing Swiss MC Plan (SMCP). The implementation has been configured by adjusting the physical density and the air gap between adjacent leaves in order to match transmission profile measurements for 6 and 15 MV beams of a Novalis TX. These measurements have been performed in water using gafchromic films and an ionization chamber at an SSD of 95 cm and a depth of 5 cm. The implementation was validated by comparing diamond measured and calculated penumbra values (80%-20%) for different field sizes and water depths. Additionally, measured and calculated dose distributions for a head and neck IMRT case using the DELTA(4) phantom have been compared. The validated HD120 MLC implementation has been used for its physical characterization. For this purpose, phase space (PS) files have been generated below the fully closed multileaf collimator (MLC) of a 40 × 22 cm(2) field size for 6 and 15 MV. The PS files have been analyzed in terms of energy spectra, mean energy, fluence, and energy fluence in the direction perpendicular to the MLC leaves and have been compared with the corresponding data using the well established Varian 80 leaf (MLC80) and Millennium M120 (M120 MLC) MLCs. Additionally, the impact of the tongue and groove design of the MLCs on dose has been characterized. RESULTS Calculated transmission values for the HD120 MLC are 1.25% and 1.34% in the central part of the field for the 6 and 15 MV beam, respectively. The corresponding ionization chamber measurements result in a transmission of 1.20% and 1.35%. Good agreement has been found for the comparison between transmission profiles resulting from MC simulations and film measurements. The simulated and measured values for the penumbra agreed within <0.5 mm for all field sizes, depths, and beam energies, and a good agreement has been found between the measured and the calculated dose distributions for the IMRT case. The total energy spectra are almost identical for the three MLCs. However, the mean energy, fluence and energy fluence are significantly different. Due to the different leaf widths of the MLCs, the shape of these distributions is different, each representing its leave structure. Due to the increase in width from the inner to the outer HD120 MLC leaves, the fluence and energy fluence clearly decrease below the outer leaves. The MLC80 and the M120 MLC resulted in an increase of the fluence and energy fluence compared with those resulted for the HD120 MLC. The dose reduction can exceed 20% compared with the dose of the open field due to the tongue and groove design of the HD120 MLC. CONCLUSIONS The HD120 MLC has been successfully implemented into the SMCP. Comparisons between MC calculations and measurements show very good agreement. The SMCP is now able to calculate accurate dose distributions for treatment plans using the HD120 MLC.

Spectral optoacoustic imaging using a scanning transducer

A scanning optoacoustic transducer was used to generate images of biological tissue in vivo. The transducer consisted of an optical fiber for delivery of laser pulses and an annular piezoelectric film. Optoacoustic waves were generated by irradiating the tissue with 5 ns long pulses from an optical parametric oscillator (OPO) which was tuned between 500 and 600 nm. The motorized scanner allowed to generate a 200x200 pixel image within 20 seconds. In images of human skin in vivo, certain structures such as the epidermis, blood capillaries and deeper blood vessels could be localized and identified by their spectral characteristics. From the spectral dependence of blood containing structures the oxygenation level of hemoglobin could be estimated.

Optoacoustic tomography using a two-dimensional optical pressure transducer and two different reconstruction algorithms

Optoacoustic pressure waves were generated by irradiating blood vessels embedded into calf cartilage with 5 ns long pulses from an optical parametric oscillator (OPO). A two dimensional pressure transducer consisting of a glass prism and a liquid layer was used to measure the laser-induced pressure waves. The pressure dependent changes of reflection at the prism-liquid interface were recorded with a time-gated CCD-camera. On the basis of the 2D pressure measurements the spatial distribution of the blood vessels was reconstructed using two different algorithms, one based on radial back projection and one using Fourier transformation. High lateral (100 micrometers ) and depth resolution (20 micrometers ) was achieved.

Dose calculation of dynamic trajectory radiotherapy using Monte Carlo

PURPOSE Using volumetric modulated arc therapy (VMAT) delivery technique gantry position, multi-leaf collimator (MLC) as well as dose rate change dynamically during the application. However, additional components can be dynamically altered throughout the dose delivery such as the collimator or the couch. Thus, the degrees of freedom increase allowing almost arbitrary dynamic trajectories for the beam. While the dose delivery of such dynamic trajectories for linear accelerators is technically possible, there is currently no dose calculation and validation tool available. Thus, the aim of this work is to develop a dose calculation and verification tool for dynamic trajectories using Monte Carlo (MC) methods. METHODS The dose calculation for dynamic trajectories is implemented in the previously developed Swiss Monte Carlo Plan (SMCP). SMCP interfaces the treatment planning system Eclipse with a MC dose calculation algorithm and is already able to handle dynamic MLC and gantry rotations. Hence, the additional dynamic components, namely the collimator and the couch, are described similarly to the dynamic MLC by defining data pairs of positions of the dynamic component and the corresponding MU-fractions. For validation purposes, measurements are performed with the Delta4 phantom and film measurements using the developer mode on a TrueBeam linear accelerator. These measured dose distributions are then compared with the corresponding calculations using SMCP. First, simple academic cases applying one-dimensional movements are investigated and second, more complex dynamic trajectories with several simultaneously moving components are compared considering academic cases as well as a clinically motivated prostate case. RESULTS The dose calculation for dynamic trajectories is successfully implemented into SMCP. The comparisons between the measured and calculated dose distributions for the simple as well as for the more complex situations show an agreement which is generally within 3% of the maximum dose or 3mm. The required computation time for the dose calculation remains the same when the additional dynamic moving components are included. CONCLUSION The results obtained for the dose comparisons for simple and complex situations suggest that the extended SMCP is an accurate dose calculation and efficient verification tool for dynamic trajectory radiotherapy. This work was supported by Varian Medical Systems.

Electron beam collimation with a photon MLC for standard electron treatments

Standard electron treatments are currently still performed using standard or molded patient-specific cut-outs placed in the electron applicator. Replacing cut-outs and electron applicators with a photon multileaf collimator (pMLC) for electron beam collimation would make standard electron treatments more efficient and would facilitate advanced treatment techniques like modulated electron radiotherapy (MERT) and mixed beam radiotherapy (MBRT). In this work, a multiple source Monte Carlo beam model for pMLC shaped electron beams commissioned at a source-to-surface distance (SSD) of 70 cm is extended for SSDs of up to 100 cm and validated for several Varian treatment units with field sizes typically used for standard electron treatments. Measurements and dose calculations agree generally within 3% of the maximal dose or 2 mm distance to agreement. To evaluate the dosimetric consequences of using pMLC collimated electron beams for standard electron treatments, pMLC-based and cut-out-based treatment plans are created for a left and a right breast boost, a sternum, a testis and a parotid gland case. The treatment plans consist of a single electron field, either alone (1E) or in combination with two 3D conformal tangential photon fields (1E2X). For each case, a pMLC plan with similar treatment plan quality in terms of dose homogeneity to the target and absolute mean dose values to the organs at risk (OARs) compared to a cut-out plan is found. The absolute mean dose to an OAR is slightly increased for pMLC-based compared to cut-out-based 1E plans if the OAR is located laterally close to the target with respect to beam direction, or if a 6 MeV electron beam is used at an extended SSD. In conclusion, treatment plans using cut-out collimation can be replaced by plans of similar treatment plan quality using pMLC collimation with accurately calculated dose distributions.

Assessing dose rate distributions in VMAT plans.

Dose rate is an essential factor in radiobiology. As modern radiotherapy delivery techniques such as volumetric modulated arc therapy (VMAT) introduce dynamic modulation of the dose rate, it is important to assess the changes in dose rate. Both the rate of monitor units per minute (MU rate) and collimation are varied over the course of a fraction, leading to different dose rates in every voxel of the calculation volume at any point in time during dose delivery. Given the radiotherapy plan and machine specific limitations, a VMAT treatment plan can be split into arc sectors between Digital Imaging and Communications in Medicine control points (CPs) of constant and known MU rate. By calculating dose distributions in each of these arc sectors independently and multiplying them with the MU rate, the dose rate in every single voxel at every time point during the fraction can be calculated. Independently calculated and then summed dose distributions per arc sector were compared to the whole arc dose calculation for validation. Dose measurements and video analysis were performed to validate the calculated datasets. A clinical head and neck, cranial and liver case were analyzed using the tool developed. Measurement validation of synthetic test cases showed linac agreement to precalculated arc sector times within ±0.4 s and doses ±0.1 MU (one standard deviation). Two methods for the visualization of dose rate datasets were developed: the first method plots a two-dimensional (2D) histogram of the number of voxels receiving a given dose rate over the course of the arc treatment delivery. In similarity to treatment planning system display of dose, the second method displays the dose rate as color wash on top of the corresponding computed tomography image, allowing the user to scroll through the variation over time. Examining clinical cases showed dose rates spread over a continuous spectrum, with mean dose rates hardly exceeding 100 cGy min(-1) for conventional fractionation. A tool to analyze dose rate distributions in VMAT plans with sub-second accuracy was successfully developed and validated. Dose rates encountered in clinical VMAT test cases show a continuous spectrum with a mean less than or near 100 cGy min(-1) for conventional fractionation.

SU-F-T-412: Spatiotemporal Distribution of Dose Rates in VMAT Plans

PURPOSE To assess dose rates in volumetric modulated arc therapy (VMAT) plans. METHODS For dose rate analysis, using Monte Carlo methods, dose distributions are calculated in sectors of consecutive DICOM control points and each multiplied with the MU-rate obtaining dose rate distributions for every voxel at every time-point during a treatment fraction. Dose rate distributions were calculated for a clinical head and neck, cranial (2 Gy fraction dose) and a SBRT liver treatment case (5 Gy per fraction). To assess optimizer influence on dose rate distributions, five conventionally fractionated head and neck cases were optimized using the progressive resolution optimizer (PRO) versions 2 and 3 (Varian Medical Systems, Palo Alto, CA). RESULTS Firstly, histograms of dose rates over the course of treatment delivery in volumes of interest were generated for the PTV and spine or brainstem of the clinical treatment plans. Secondly, dose rates were assessed as color wash overlay on top of CT slices, allowing the user to explore the dose rate variation over time. For the cases considered, dose rates were spread out over a continuous spectrum with a mean value in the lower half. PRO3 generated plans showed higher maximum dose rates compared to PRO2 (699.80 cGy/min vs. 543.46 cGy/min) and modulation of MU rates in blocks as opposed to gradual modulation. CONCLUSION Calculation and visualization of dose rate distributions in VMAT have been performed. Dose rates in clinical relevant VMAT cases are spread out over a continuous spectrum with a mean in the lower part for both conventional and SBRT test cases. Different optimization algorithms may lead to substantially different dose rate distributions in the same cases. The development of the Swiss Monte Carlo Plan framework has been partially supported by Varian Medical Systems.

SU-F-T-586: Pre-Treatment QA of InCise2 MLC Plans On a Cyberknife-M6 Using the Delta4 System in SBRT

PURPOSE Performing pre-treatment quality assurance (QA) with the Delta4 system (ScandiDos Inc., Madison, WI) is well established for linac-based radiotherapy. This is not true when using a Cyberknife (Accuray Inc., Sunnyvale, CA) where, typically film-based QA is applied. The goal of this work was to test the feasibility to use the Delta4 system for pre-treatment QA for stereotactic body radiation therapy (SBRT) using a Cyberknife-M6 equipped with the InCise2 multileaf collimator (MLC). METHODS In order to perform measurements without accelerator pulse signal, the Tomotherapy option within the Delta4 software was used. Absolute calibration of the Delta4 phantom was performed using a 10×10 cm2 field shaped by the InCise2 MLC of the Cyberknife-M6. Five fiducials were attached to the Delta4 phantom in order to be able to track the phantom before and during measurements. For eight SBRT treatment plans (two liver, two prostate, one lung, three bone metastases) additional verification plans were recalculated on the Delta4 phantom using MultiPlan. Dicom data was exported from MultiPlan and was adapted in order to be compatible with the Delta4 software. The measured and calculated dose distributions were compared using the gamma analysis of the Delta4 system. RESULTS All eight SBRT plans were successfully measured with the aid of the Delta4 system. In the mean, 98.0±1.9%, 95.8±4.1% and 88.40±11.4% of measured dose points passed the gamma analysis using a global dose deviation criterion of 3% (100% corresponds to the dose maximum) and a distance-to-agreement criterion of 3 mm, 2 mm and 1 mm, respectively, and a threshold of 20%. CONCLUSION Pre-treatment QA of SBRT plans using the Delta4 system on a Cyberknife-M6 is feasible. Measured dose distributions of SBRT plans showed clinically acceptable agreement with the corresponding calculated dose distributions.