Monte Carlo treatment planning for photon and electron beams

PRIMO: a graphical environment for the Monte Carlo simulation of Varian and Elekta linacs

BACKGROUND

The accurate Monte Carlo simulation of a linac requires a detailed description of its geometry and the application of elaborate variance-reduction techniques for radiation transport. Both tasks entail a substantial coding effort and demand advanced knowledge of the intricacies of the Monte Carlo system being used.

METHODS

PRIMO, a new Monte Carlo system that allows the effortless simulation of most Varian and Elekta linacs, including their multileaf collimators and electron applicators, is introduced. PRIMO combines (1) accurate physics from the PENELOPE code, (2) dedicated variance-reduction techniques that significantly reduce the computation time, and (3) a user-friendly graphical interface with tools for the analysis of the generated data. PRIMO can tally dose distributions in phantoms and computerized tomographies, handle phase-space files in IAEA format, and import structures (planning target volumes, organs at risk) in the DICOM RT-STRUCT standard.

RESULTS

A prostate treatment, conformed with a high definition Millenium multileaf collimator (MLC 120HD) from a Varian Clinac 2100 C/D, is presented as an example. The computation of the dose distribution in 1.86×3.00×1.86 mm3 voxels with an average 2% standard statistical uncertainty, performed on an eight-core Intel Xeon at 2.67 GHz, took 1.8 h-excluding the patient-independent part of the linac, which required 3.8 h but it is simulated only once.

CONCLUSION

PRIMO is a self-contained user-friendly system that facilitates the Monte Carlo simulation of dose distributions produced by most currently available linacs. This opens the door for routine use of Monte Carlo in clinical research and quality assurance purposes. It is free software that can be downloaded from http://www.primoproject.net.

Comparison of different breast planning techniques and algorithms for radiation therapy treatment

This work aims at investigating the impact of treating breast cancer using different radiation therapy (RT) techniques--forwardly-planned intensity-modulated, f-IMRT, inversely-planned IMRT and dynamic conformal arc (DCART) RT--and their effects on the whole-breast irradiation and in the undesirable irradiation of the surrounding healthy tissues. Two algorithms of iPlan BrainLAB treatment planning system were compared: Pencil Beam Convolution (PBC) and commercial Monte Carlo (iMC). Seven left-sided breast patients submitted to breast-conserving surgery were enrolled in the study. For each patient, four RT techniques--f-IMRT, IMRT using 2-fields and 5-fields (IMRT2 and IMRT5, respectively) and DCART - were applied. The dose distributions in the planned target volume (PTV) and the dose to the organs at risk (OAR) were compared analyzing dose-volume histograms; further statistical analysis was performed using IBM SPSS v20 software. For PBC, all techniques provided adequate coverage of the PTV. However, statistically significant dose differences were observed between the techniques, in the PTV, OAR and also in the pattern of dose distribution spreading into normal tissues. IMRT5 and DCART spread low doses into greater volumes of normal tissue, right breast, right lung and heart than tangential techniques. However, IMRT5 plans improved distributions for the PTV, exhibiting better conformity and homogeneity in target and reduced high dose percentages in ipsilateral OAR. DCART did not present advantages over any of the techniques investigated. Differences were also found comparing the calculation algorithms: PBC estimated higher doses for the PTV, ipsilateral lung and heart than the iMC algorithm predicted.

Monte Carlo-based dose calculation engine for minibeam radiation therapy

Minibeam radiation therapy (MBRT) is an innovative radiotherapy approach based on the well-established tissue sparing effect of arrays of quasi-parallel micrometre-sized beams. In order to guide the preclinical trials in progress at the European Synchrotron Radiation Facility (ESRF), a Monte Carlo-based dose calculation engine has been developed and successfully benchmarked with experimental data in anthropomorphic phantoms. Additionally, a realistic example of treatment plan is presented. Despite the micron scale of the voxels used to tally dose distributions in MBRT, the combination of several efficiency optimisation methods allowed to achieve acceptable computation times for clinical settings (approximately 2 h). The calculation engine can be easily adapted with little or no programming effort to other synchrotron sources or for dose calculations in presence of contrast agents.

Impact of the number of discrete angles used during dose computation for TomoTherapy treatments

PURPOSE

To quantify systematically the effect on accuracy of discretizing gantry rotation during the dose calculation process of TomoTherapy treatments.

METHODS

Up to version 4.0.x included, TomoTherapy treatment planning system (TPS) approximates gantry rotation by computing dose from 51 discrete angles corresponding to the center of the projections used to control the binary multileaf collimator. Potential effects on dose computation accuracy for off-axis targets and low modulation factors have been shown previously for a few treatment configurations. In versions 4.1.x and later, TomoTherapy oversamples the projections to better account for gantry rotation, but only during full scatter optimization and final calculation (i.e., not during optimization in "beamlet" mode). The effect on accuracy of changing the number of angles was quantified with the following framework: (1) predict the impact of the discretization of gantry rotation for various modulation factors, target sizes, and off-axis positions using a simplified analytical algorithm; (2) perform regular quality assurance using measurements with EDR2 radiographic films; (3) isolating the effect of changing the number of discretized angles only (51, 153, and 459) using a previously validated Monte Carlo model (TomoPen). The diameters of the targets were 2, 3, and 5 cm; off-axis central positions of target volumes were 5, 10 and 15, and 17 cm (when accepted by the treatment unit); planned modulation factors were 1.3 and 2.0.

RESULTS

For extreme configurations (3 cm tumor, 1.3 modulation factor, 15 cm off-axis position), effects on dose distributions were significant with 89.3% and 95.4% of the points passing gamma tests with 2%∕2 mm and 3%∕3 mm criteria, respectively, for TPS software version 4.0.x (51 gantry angles). The passing rate was 100% for both gamma criteria for the 4.1.x version (153 gantry angles). Those differences could be attributed almost completely to gantry motion discretization using TomoPen. Using 51 gantry angles for dose computation, TomoPen reproduced within statistical uncertainties (<1% standard deviation) dose distributions computed with version 4.0.x. Using 153 and 459 gantry angles, TomoPen reproduced within statistical uncertainties measurements and dose distributions computed with version 4.1.x.

CONCLUSIONS

When low modulation factors and significant off-axis positions are used, accounting for gantry rotation during dose computation using at least 153 gantry angles is required to ensure optimal accuracy.

Helical tomotherapy for SIB and hypo-fractionated treatments in lung carcinomas: a 4D Monte Carlo treatment planning study

PURPOSE

To evaluate the impact of intra-fraction motion induced by regular breathing on treatment quality for helical tomotherapy treatments.

MATERIAL AND METHODS

Four patients treated by simultaneous-integrated boost (SIB) and three by hypo-fractionated stereotactic treatments (hypo-fractionated, 18 Gy/fraction) were included. All patients were coached to ensure regular breathing. For the SIB group, the tumor volume was delineated using CT information only (CTV(CT)) and the boost region was based on PET information (GTV(PET), no CTV extension). In the hypo-fractionated group, a GTV based on CT information was contoured. In both groups, ITVs were defined according to 4D data. The PTV included the ITV plus a setup error margin. The treatment was planned using the tomotherapy TPS on 3D CT images. In order to verify the impact of intra-fraction motion and interplay effects, dose calculations were performed using a previously validated Monte Carlo model of tomotherapy (TomoPen): first on the planning 3D CT ("planned dose") and second, on the 10 phases of the 4D scan. For the latter, two dose distributions, termed "interplay simulated" or "no interplay" were computed with and without beamlet-phase correlation over the 10 phases and combined using deformable dose registration.

RESULTS

In all cases, DVHs of "interplay simulated" dose distributions complied within 1% of the original clinical objectives used for planning, defined according to ICRU (report 83) and RTOG (trials 0236 and 0618) recommendations, for SIB and hypo-fractionated groups, respectively. For one patient in the hypo-fractionated group, D(mean) to the CTV(CT) was 2.6% and 2.5% higher than "planned" for "interplay simulated" and "no interplay", respectively.

CONCLUSION

For the patients included in this study, assuming regular breathing, the results showed that interplay of breathing and tomotherapy delivery motions did not affect significantly plan delivery accuracy. Hence, accounting for intra-fraction motion through the definition of an ITV volume was sufficient to ensure tumor coverage.

Investigation of pulsed low dose rate radiotherapy using dynamic arc delivery techniques

There has been no consensus standard of care to treat recurrent cancer patients who have previously been irradiated. Pulsed low dose rate (PLDR) external beam radiotherapy has the potential to reduce normal tissue toxicities while still providing significant tumor control for recurrent cancers. This work investigates the dosimetry feasibility of PLDR treatment using dynamic arc delivery techniques. Five treatment sites were investigated in this study including breast, pancreas, prostate, head and neck, and lung. Dynamic arc plans were generated using the Varian Eclipse system and the RapidArc delivery technique with 6 and 10 MV photon beams. Each RapidArc plan consisted of two full arcs and the plan was delivered five times to achieve a daily dose of 200 cGy. The dosimetry requirement was to deliver approximately 20 cGy/arc with a 3 min interval to achieve an effective dose rate of 6.7 cGy min⁻¹. Monte Carlo simulations were performed to calculate the actual dose delivered to the planning target volume (PTV) per arc taking into account beam attenuation/scattering and intensity modulation. The maximum, minimum and mean doses to the PTV were analyzed together with the dose volume histograms and isodose distributions. The dose delivery for the five plans was validated using solid water phantoms inserted with an ionization chamber and film, and a cylindrical detector array. Two intensity-modulated arcs were used to efficiently deliver the PLDR plans that provided conformal dose distributions for treating complex recurrent cancers. For the five treatment sites, the mean PTV dose ranged from 18.9 to 22.6 cGy/arc. For breast, the minimum and maximum PTV dose was 8.3 and 35.2 cGy/arc, respectively. The PTV dose varied between 12.9 and 27.5 cGy/arc for pancreas, 12.6 and 28.3 cGy/arc for prostate, 12.1 and 30.4 cGy/arc for H&N, and 16.2 and 27.6 cGy/arc for lung. Advanced radiation therapy can provide superior target coverage and normal tissue sparing for PLDR reirradiation of recurrent cancers, which can be delivered using dynamic arc delivery techniques with ten full arcs and an effective dose rate of 6.7 ± 4.0 cGy min⁻¹.

X-Ray Irradiation as a Microbial Intervention Strategy for Food

First recognized in 1895, X-ray irradiation soon became a breakthrough diagnostic tool for the dental and medical professions. However, the food industry remained slow to adopt X-ray irradiation as a means for controlling insects and microbial contaminants in food, instead using gamma and electron beam (E-beam) irradiation. However, the reinvention of X-ray machines with increased efficiency, combined with recent developments in legislation and engineering, is now allowing X-ray to actively compete with gamma irradiation and E-beam as a microbial reduction strategy for foods. This review summarizes the historical developments of X-rays and discusses the key technological advances over the past two decades that now have led to the development of several different X-ray irradiators capable of enhancing the safety and shelf life of many heat-sensitive products, including lettuce, spinach, tomatoes, and raw almonds, all of which have been linked to high profile outbreaks of foodborne illness.

A PENELOPE-based system for the automated Monte Carlo simulation of clinacs and voxelized geometries-application to far-from-axis fields

PURPOSE

Two new codes, PENEASY and PENEASYLINAC, which automate the Monte Carlo simulation of Varian Clinacs of the 600, 1800, 2100, and 2300 series, together with their electron applicators and multileaf collimators, are introduced. The challenging case of a relatively small and far-from-axis field has been studied with these tools.

METHODS

PENEASY is a modular, general-purpose main program for the PENELOPE Monte Carlo system that includes various source models, tallies and variance-reduction techniques (VRT). The code includes a new geometry model that allows the superposition of voxels and objects limited by quadric surfaces. A variant of the VRT known as particle splitting, called fan splitting, is also introduced. PENEASYLINAC, in turn, automatically generates detailed geometry and configuration files to simulate linacs with PENEASY. These tools are applied to the generation of phase-space files, and of the corresponding absorbed dose distributions in water, for two 6 MV photon beams from a Varian Clinac 2100 C∕D: a 40 × 40 cm(2) centered field; and a 3 × 5 cm(2) field centered at (4.5, -11.5) cm from the beam central axis. This latter configuration implies the largest possible over-traveling values of two of the jaws. Simulation results for the depth dose and lateral profiles at various depths are compared, by using the gamma index, with experimental values obtained with a PTW 31002 ionization chamber. The contribution of several VRTs to the computing speed of the more demanding off-axis case is analyzed.

RESULTS

For the 40 × 40 cm(2) field, the percentages γ(1) and γ(1.2) of voxels with gamma indices (using 0.2 cm and 2% criteria) larger than unity and larger than 1.2 are 0.2% and 0%, respectively. For the 3 × 5 cm(2) field, γ(1) = 0%. These figures indicate an excellent agreement between simulation and experiment. The dose distribution for the off-axis case with voxels of 2.5 × 2.5 × 2.5 mm(3) and an average standard statistical uncertainty of 2% (1σ) is computed in 3.1 h on a single core of a 2.8 GHz Intel Core 2 Duo processor. This result is obtained with the optimal combination of the tested VRTs. In particular, fan splitting for the off-axis case accelerates execution by a factor of 240 with respect to standard particle splitting.

CONCLUSIONS

PENEASY and PENEASYLINAC can simulate the considered Varian Clinacs both in an accurate and efficient manner. Fan splitting is crucial to achieve simulation results for the off-axis field in an affordable amount of CPU time. Work to include Elekta linacs and to develop a graphical interface that will facilitate user input is underway.

First steps towards a fast-neutron therapy planning program

BACKGROUND

The Monte Carlo code GEANT4 was used to implement first steps towards a treatment planning program for fast-neutron therapy at the FRM II research reactor in Garching, Germany. Depth dose curves were calculated inside a water phantom using measured primary neutron and simulated primary photon spectra and compared with depth dose curves measured earlier. The calculations were performed with GEANT4 in two different ways, simulating a simple box geometry and splitting this box into millions of small voxels (this was done to validate the voxelisation procedure that was also used to voxelise the human body).

RESULTS

In both cases, the dose distributions were very similar to those measured in the water phantom, up to a depth of 30 cm. In order to model the situation of patients treated at the FRM II MEDAPP therapy beamline for salivary gland tumors, a human voxel phantom was implemented in GEANT4 and irradiated with the implemented MEDAPP neutron and photon spectra. The 3D dose distribution calculated inside the head of the phantom was similar to the depth dose curves in the water phantom, with some differences that are explained by differences in elementary composition. The lateral dose distribution was studied at various depths. The calculated cumulative dose volume histograms for the voxel phantom show the exposure of organs at risk surrounding the tumor.

CONCLUSIONS

In order to minimize the dose to healthy tissue, a conformal treatment is necessary. This can only be accomplished with the help of an advanced treatment planning system like the one developed here. Although all calculations were done for absorbed dose only, any biological dose weighting can be implemented easily, to take into account the increased radiobiological effectiveness of neutrons compared to photons.

3D visualisation of the stochastic patterns of the radial dose in nano-volumes by a Monte Carlo simulation of HZE ion track structure

The description of energy deposition by high charge and energy (HZE) nuclei is of importance for space radiation risk assessment and due to their use in hadrontherapy. Such ions deposit a large fraction of their energy within the so-called core of the track and a smaller proportion in the penumbra (or track periphery). We study the stochastic patterns of the radial dependence of energy deposition using Monte Carlo track structure codes RITRACKS and RETRACKS, that were used to simulate HZE tracks and calculate energy deposition in voxels of 40 nm. The simulation of a (56)Fe(26+) ion of 1 GeV u(-1) revealed zones of high-energy deposition which maybe found as far as a few millimetres away from the track core in some simulations. The calculation also showed that ∼43 % of the energy was deposited in the penumbra. These 3D stochastic simulations combined with a visualisation interface are a powerful tool for biophysicists which may be used to study radiation-induced biological effects such as double strand breaks and oxidative damage and the subsequent cellular and tissue damage processing and signalling.

Fast, accurate photon beam accelerator modeling using BEAMnrc: a systematic investigation of efficiency enhancing methods and cross-section data

In this work, an investigation of efficiency enhancing methods and cross-section data in the BEAMnrc Monte Carlo (MC) code system is presented. Additionally, BEAMnrc was compared with VMC++, another special-purpose MC code system that has recently been enhanced for the simulation of the entire treatment head. BEAMnrc and VMC++ were used to simulate a 6 MV photon beam from a Siemens Primus linear accelerator (linac) and phase space (PHSP) files were generated at 100 cm source-to-surface distance for the 10 x 10 and 40 x 40 cm2 field sizes. The BEAMnrc parameters/techniques under investigation were grouped by (i) photon and bremsstrahlung cross sections, (ii) approximate efficiency improving techniques (AEITs), (iii) variance reduction techniques (VRTs), and (iv) a VRT (bremsstrahlung photon splitting) in combination with an AEIT (charged particle range rejection). The BEAMnrc PHSP file obtained without the efficiency enhancing techniques under study or, when not possible, with their default values (e.g., EXACT algorithm for the boundary crossing algorithm) and with the default cross-section data (PEGS4 and Bethe-Heitler) was used as the "base line" for accuracy verification of the PHSP files generated from the different groups described previously. Subsequently, a selection of the PHSP files was used as input for DOSXYZnrc-based water phantom dose calculations, which were verified against measurements. The performance of the different VRTs and AEITs available in BEAMnrc and of VMC++ was specified by the relative efficiency, i.e., by the efficiency of the MC simulation relative to that of the BEAMnrc base-line calculation. The highest relative efficiencies were approximately 935 (approximately 111 min on a single 2.6 GHz processor) and approximately 200 (approximately 45 min on a single processor) for the 10 x 10 field size with 50 million histories and 40 x 40 cm2 field size with 100 million histories, respectively, using the VRT directional bremsstrahlung splitting (DBS) with no electron splitting. When DBS was used with electron splitting and combined with augmented charged particle range rejection, a technique recently introduced in BEAMnrc, relative efficiencies were approximately 420 (approximately 253 min on a single processor) and approximately 175 (approximately 58 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. Calculations of the Siemens Primus treatment head with VMC++ produced relative efficiencies of approximately 1400 (approximately 6 min on a single processor) and approximately 60 (approximately 4 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. BEAMnrc PHSP calculations with DBS alone or DBS in combination with charged particle range rejection were more efficient than the other efficiency enhancing techniques used. Using VMC++, accurate simulations of the entire linac treatment head were performed within minutes on a single processor. Noteworthy differences (+/- 1%-3%) in the mean energy, planar fluence, and angular and spectral distributions were observed with the NIST bremsstrahlung cross sections compared with those of Bethe-Heitler (BEAMnrc default bremsstrahlung cross section). However, MC calculated dose distributions in water phantoms (using combinations of VRTs/AEITs and cross-section data) agreed within 2% of measurements. Furthermore, MC calculated dose distributions in a simulated water/air/water phantom, using NIST cross sections, were within 2% agreement with the BEAMnrc Bethe-Heitler default case.

Fast Monte Carlo simulation on a voxelized human phantom deformed to a patient

PURPOSE

A method for performing fast simulations of absorbed dose using a patient's computerized tomography (CT) scan without explicitly relying on a calibration curve is presented.

METHODS

The method is based on geometrical deformations performed on a standard voxelized human phantom. This involves spatially transforming the human phantom to align it with the patient CT image. Since the chemical composition and density of each voxel are given in the phantom data, a calibration curve is not used in the proposed method. For this study, the Monte Carlo (MC) code PENELOPE has been used as the simulation of reference. The results obtained with PENELOPE simulations are compared to those obtained with PENFAST and with the collapsed cone convolution algorithm implemented in a commercial treatment planning system.

RESULTS

The comparisons of the absorbed doses calculated with the different algorithms on two patient CTs and the corresponding deformed phantoms show a maximum distance to agreement of 2 mm, and in general, the obtained absorbed dose distributions are compatible within the reached statistical uncertainty. The validity of the deformation method for a broad range of patients is shown using MC simulations in random density phantoms. A PENFAST simulation of a 6 MV photon beam impinging on a patient CT reaches 2% statistical uncertainty in the absorbed dose, in a 0.1 cm3 voxel along the central axis, in 10 min running on a single core of a 2.8 GHz CPU.

CONCLUSIONS

The proposed method of the absorbed dose calculation in a deformed voxelized phantom allows for dosimetric studies in the geometry of a patient CT scan. This is due to the fact that the chemical composition and material density of the phantom are known. Furthermore, simulation using the phantom geometry can provide dosimetric information for each organ. The method can be used for quality assurance procedures. In relation to PENFAST, it is shown that a purely condensed-history algorithm (class I) can be used for absorbed dose estimation in patient CTs.

Low-coupled parallel strategy for Monte Carlo radiation dose calculation

In this work a parallelization technique for Monte Carlo calculation of the radiation absorbed dose is presented. The proposed method reduces the overall time of calculation exploiting the probabilistic nature of the problem. Reducing dose calculation times to meet practical application of use of results is a topic of great interest in the Monte Carlo simulation, dosimetry evaluation and treatment planning. The general problem of dose calculation parallelization was analyzed. The optimal number of calculation units for a given context was also determined.

A preliminary study of in-house Monte Carlo simulations: an integrated Monte Carlo verification system

PURPOSE

To develop an infrastructure for the integrated Monte Carlo verification system (MCVS) to verify the accuracy of conventional dose calculations, which often fail to accurately predict dose distributions, mainly due to inhomogeneities in the patient's anatomy, for example, in lung and bone.

METHODS AND MATERIALS

The MCVS consists of the graphical user interface (GUI) based on a computational environment for radiotherapy research (CERR) with MATLAB language. The MCVS GUI acts as an interface between the MCVS and a commercial treatment planning system to import the treatment plan, create MC input files, and analyze MC output dose files. The MCVS consists of the EGSnrc MC codes, which include EGSnrc/BEAMnrc to simulate the treatment head and EGSnrc/DOSXYZnrc to calculate the dose distributions in the patient/phantom. In order to improve computation time without approximations, an in-house cluster system was constructed.

RESULTS

The phase-space data of a 6-MV photon beam from a Varian Clinac unit was developed and used to establish several benchmarks under homogeneous conditions. The MC results agreed with the ionization chamber measurements to within 1%. The MCVS GUI could import and display the radiotherapy treatment plan created by the MC method and various treatment planning systems, such as RTOG and DICOM-RT formats. Dose distributions could be analyzed by using dose profiles and dose volume histograms and compared on the same platform. With the cluster system, calculation time was improved in line with the increase in the number of central processing units (CPUs) at a computation efficiency of more than 98%.

CONCLUSIONS

Development of the MCVS was successful for performing MC simulations and analyzing dose distributions.

Evaluation of PENFAST--a fast Monte Carlo code for dose calculations in photon and electron radiotherapy treatment planning

The aim of the present study is to demonstrate the potential of accelerated dose calculations, using the fast Monte Carlo (MC) code referred to as PENFAST, rather than the conventional MC code PENELOPE, without losing accuracy in the computed dose. For this purpose, experimental measurements of dose distributions in homogeneous and inhomogeneous phantoms were compared with simulated results using both PENELOPE and PENFAST. The simulations and experiments were performed using a Saturne 43 linac operated at 12 MV (photons), and at 18 MeV (electrons). Pre-calculated phase space files (PSFs) were used as input data to both the PENELOPE and PENFAST dose simulations. Since depth-dose and dose profile comparisons between simulations and measurements in water were found to be in good agreement (within +/-1% to 1 mm), the PSF calculation is considered to have been validated. In addition, measured dose distributions were compared to simulated results in a set of clinically relevant, inhomogeneous phantoms, consisting of lung and bone heterogeneities in a water tank. In general, the PENFAST results agree to within a 1% to 1 mm difference with those produced by PENELOPE, and to within a 2% to 2 mm difference with measured values. Our study thus provides a pre-clinical validation of the PENFAST code. It also demonstrates that PENFAST provides accurate results for both photon and electron beams, equivalent to those obtained with PENELOPE. CPU time comparisons between both MC codes show that PENFAST is generally about 9-21 times faster than PENELOPE.

An overview of Monte Carlo treatment planning for radiotherapy

The implementation of Monte Carlo dose calculation algorithms in clinical radiotherapy treatment planning systems has been anticipated for many years. Despite a continuous increase of interest in Monte Carlo Treatment Planning (MCTP), its introduction into clinical practice has been delayed by the extent of calculation time required. The development of newer and faster MC codes is behind the commercialisation of the first MC-based treatment planning systems. The intended scope of this article is to provide the reader with a compact 'primer' on different approaches to MCTP with particular attention to the latest developments in the field.

Anatomical imaging for radiotherapy

The goal of radiation therapy is to achieve maximal therapeutic benefit expressed in terms of a high probability of local control of disease with minimal side effects. Physically this often equates to the delivery of a high dose of radiation to the tumour or target region whilst maintaining an acceptably low dose to other tissues, particularly those adjacent to the target. Techniques such as intensity modulated radiotherapy (IMRT), stereotactic radiosurgery and computer planned brachytherapy provide the means to calculate the radiation dose delivery to achieve the desired dose distribution. Imaging is an essential tool in all state of the art planning and delivery techniques: (i) to enable planning of the desired treatment, (ii) to verify the treatment is delivered as planned and (iii) to follow-up treatment outcome to monitor that the treatment has had the desired effect. Clinical imaging techniques can be loosely classified into anatomic methods which measure the basic physical characteristics of tissue such as their density and biological imaging techniques which measure functional characteristics such as metabolism. In this review we consider anatomical imaging techniques. Biological imaging is considered in another article. Anatomical imaging is generally used for goals (i) and (ii) above. Computed tomography (CT) has been the mainstay of anatomical treatment planning for many years, enabling some delineation of soft tissue as well as radiation attenuation estimation for dose prediction. Magnetic resonance imaging is fast becoming widespread alongside CT, enabling superior soft-tissue visualization. Traditionally scanning for treatment planning has relied on the use of a single snapshot scan. Recent years have seen the development of techniques such as 4D CT and adaptive radiotherapy (ART). In 4D CT raw data are encoded with phase information and reconstructed to yield a set of scans detailing motion through the breathing, or cardiac, cycle. In ART a set of scans is taken on different days. Both allow planning to account for variability intrinsic to the patient. Treatment verification has been carried out using a variety of technologies including: MV portal imaging, kV portal/fluoroscopy, MVCT, conebeam kVCT, ultrasound and optical surface imaging. The various methods have their pros and cons. The four x-ray methods involve an extra radiation dose to normal tissue. The portal methods may not generally be used to visualize soft tissue, consequently they are often used in conjunction with implanted fiducial markers. The two CT-based methods allow measurement of inter-fraction variation only. Ultrasound allows soft-tissue measurement with zero dose but requires skilled interpretation, and there is evidence of systematic differences between ultrasound and other data sources, perhaps due to the effects of the probe pressure. Optical imaging also involves zero dose but requires good correlation between the target and the external measurement and thus is often used in conjunction with an x-ray method. The use of anatomical imaging in radiotherapy allows treatment uncertainties to be determined. These include errors between the mean position at treatment and that at planning (the systematic error) and the day-to-day variation in treatment set-up (the random error). Positional variations may also be categorized in terms of inter- and intra-fraction errors. Various empirical treatment margin formulae and intervention approaches exist to determine the optimum strategies for treatment in the presence of these known errors. Other methods exist to try to minimize error margins drastically including the currently available breath-hold techniques and the tracking methods which are largely in development. This paper will review anatomical imaging techniques in radiotherapy and how they are used to boost the therapeutic benefit of the treatment.

Monte Carlo simulation of a realistic anatomical phantom described by triangle meshes: application to prostate brachytherapy imaging

PURPOSE

Monte Carlo codes can simulate the transport of radiation within matter with high accuracy and can be used to study medical applications of ionising radiations. The aim of our work was to develop a Monte Carlo code capable of generating projection images of the human body. In order to obtain clinically realistic images a detailed anthropomorphic phantom was prepared. These two simulation tools are intended to study the multiple applications of imaging in radiotherapy, from image guided treatments to portal imaging.

METHODS

We adapted the general-purpose code PENELOPE 2006 to simulate a radiation source, an ideal digital detector, and a realistic model of the patient anatomy. The anthropomorphic phantom was developed using computer-aided design tools, and is based on the NCAT phantom. The surface of each organ is modelled using a closed triangle mesh, and the full phantom contains 330 organs and more than 5 million triangles. A novel object-oriented geometry package, which includes an octree structure to sort the triangles, has been developed to use this complex geometry with PENELOPE.

RESULTS

As an example of the capabilities of the new code, projection images of the human pelvis region were simulated. Radioactive seeds were included inside the phantom's prostate. Therefore, the resulting simulated images resemble what would be obtained in a clinical procedure to assess the positioning of the seeds in a prostate brachytherapy treatment.

CONCLUSIONS

The new code can produce projection images of the human body that are comparable to those obtained by a real imaging system (within the limitations of the anatomical phantom and the detector model). The simulated images can be used to study and optimise an imaging task (i.e., maximise the object detectability, minimise the delivered dose, find the optimum beam energy, etc.). Since PENELOPE can simulate radiation from 50 eV to 1 GeV, the code can also be used to simulate radiotherapy treatments and portal imaging. Using the octree data structure, the new geometry model does not significantly increase the computing time when compared to the simulation of a much simpler quadric geometry. In conclusion, we have shown that it is feasible to use PENELOPE and a complex triangle mesh geometry to simulate real medical physics applications.