Table of contents

Polymer gel dosimeters are fabricated from radiation sensitive chemicals which, upon irradiation, polymerize as a function of the absorbed radiation dose. These gel dosimeters, with the capacity to uniquely record the radiation dose distribution in three-dimensions (3D), have specific advantages when compared to one-dimensional dosimeters, such as ion chambers, and two-dimensional dosimeters, such as film. These advantages are particularly significant in dosimetry situations where steep dose gradients exist such as in intensity-modulated radiation therapy (IMRT) and stereotactic radiosurgery. Polymer gel dosimeters also have specific advantages for brachytherapy dosimetry. Potential dosimetry applications include those for low-energy x-rays, high-linear energy transfer (LET) and proton therapy, radionuclide and boron capture neutron therapy dosimetries. These 3D dosimeters are radiologically soft-tissue equivalent with properties that may be modified depending on the application. The 3D radiation dose distribution in polymer gel dosimeters may be imaged using magnetic resonance imaging (MRI), optical-computerized tomography (optical-CT), x-ray CT or ultrasound. The fundamental science underpinning polymer gel dosimetry is reviewed along with the various evaluation techniques. Clinical dosimetry applications of polymer gel dosimetry are also presented.

State-of-the-art small-animal single photon emission computed tomography (SPECT) enables sub-half-mm resolution imaging of radio-labelled molecules. Due to severe photon penetration through pinhole edges, current multi-pinhole SPECT is not suitable for high-resolution imaging of photons with high energies, such as the annihilation photons emitted by positron emitting tracers (511 keV). To deal with this edge penetration, we introduce here clustered multi-pinhole SPECT (CMP): each pinhole in a cluster has a narrow opening angle to reduce photon penetration. Using simulations, CMP is compared with (i) a collimator with traditional pinholes that is currently used for sub-half-mm imaging of SPECT isotopes (U-SPECT-II), and (ii), like (i) but with collimator thickness adapted to image high-energy photons (traditional multi-pinhole SPECT, TMP). At 511 keV, U-SPECT-II is able to resolve the 0.9 mm rods of an iteratively reconstructed Jaszczak-like capillary hot rod phantom, and while TMP only leads to small improvements, CMP can resolve rods as small as 0.7 mm. Using a digital tumour phantom, we show that CMP resolves many details not assessable with standard USPECT-II and TMP collimators. Furthermore, CMP makes it possible to visualize uptake of positron emitting tracers in sub-compartments of a digital mouse striatal brain phantom. This may open up unique possibilities for analysing processes such as those underlying the function of neurotransmitter systems. Additional potential of CMP may include (i) the imaging of other high-energy single-photon emitters (e.g. I-131) and (ii) localized imaging of positron emitting tracers simultaneously with single photon emitters, with an even better resolution than coincidence PET.

Achieving accurate small field dosimetry is challenging. This study investigates the utility of a radiochromic plastic PRESAGE® read with optical-CT for the acquisition of radiosurgery field commissioning data from a Novalis Tx system with a high-definition multileaf collimator (HDMLC). Total scatter factors (S_{c, p}), beam profiles, and penumbrae were measured for five different radiosurgery fields (5, 10, 20, 30 and 40 mm) using a commercially available optical-CT scanner (OCTOPUS, MGS Research). The percent depth dose (PDD), beam profile and penumbra of the 10 mm field were also measured using a higher resolution in-house prototype CCD-based scanner. Gafchromic EBT® film was used for independent verification. Measurements of S_{c, p} made with PRESAGE® and film agreed with mini-ion chamber commissioning data to within 4% for every field (range 0.2–3.6% for PRESAGE®, and 1.6–3.6% for EBT). PDD, beam profile and penumbra measurements made with the two PRESAGE®/optical-CT systems and film showed good agreement with the high-resolution diode commissioning measurements with a competitive resolution (0.5 mm pixels). The in-house prototype optical-CT scanner allowed much finer resolution compared with previous applications of PRESAGE®. The advantages of the PRESAGE® system for small field dosimetry include 3D measurements, negligible volume averaging, directional insensitivity, an absence of beam perturbations, energy and dose rate independence.

X-ray scatter is a major cause of nonlinearity in densitometry measurements using digital mammography. Previous scatter correction techniques have primarily used a single scatter point spread function to estimate x-ray scatter. In this study, a new algorithm to correct x-ray scatter based on image convolution was implemented using a spatially variant scatter point spread function which is energy and thickness dependent. The scatter kernel was characterized in terms of its scattering fraction (SF) and scatter radial extent (k) on uniform Lucite phantoms with thickness of 0.8–8.0 cm. The algorithm operates on a pixel-by-pixel basis by grouping pixels of similar thicknesses into a series of mask images that are individually deconvolved using Fourier image analysis with a distinct kernel for each image. The algorithm was evaluated with three Lucite step phantoms and one anthropomorphic breast phantom using a full-field digital mammography system at energies of 24, 28, 31 and 49 kVp. The true primary signal was measured with a multi-hole collimator. The effect on image quality was also evaluated. For all 16 studies, the average mean percentage error in estimating the true primary signal was found to be −2.13% and the average rms percentage error was 2.60%. The image quality was seen to improve at every energy up to 25% at 49 kVp. The results indicate that a technique based on a spatially variant scatter point spread function can accurately estimate x-ray scatter.

Effective delivery of adaptive radiotherapy requires locating the target with high precision in real time. System latency caused by data acquisition, streaming, processing and delivery control necessitates prediction. Prediction is particularly challenging for highly mobile targets such as thoracic and abdominal tumors undergoing respiration-induced motion. The complexity of the respiratory motion makes it difficult to build and justify explicit models. In this study, we honor the intrinsic uncertainties in respiratory motion and propose a statistical treatment of the prediction problem. Instead of asking for a deterministic covariate–response map and a unique estimate value for future target position, we aim to obtain a distribution of the future target position (response variable) conditioned on the observed historical sample values (covariate variable). The key idea is to estimate the joint probability distribution (pdf) of the covariate and response variables using an efficient kernel density estimation method. Then, the problem of identifying the distribution of the future target position reduces to identifying the section in the joint pdf based on the observed covariate. Subsequently, estimators are derived based on this estimated conditional distribution. This probabilistic perspective has some distinctive advantages over existing deterministic schemes: (1) it is compatible with potentially inconsistent training samples, i.e., when close covariate variables correspond to dramatically different response values; (2) it is not restricted by any prior structural assumption on the map between the covariate and the response; (3) the two-stage setup allows much freedom in choosing statistical estimates and provides a full nonparametric description of the uncertainty for the resulting estimate. We evaluated the prediction performance on ten patient RPM traces, using the root mean squared difference between the prediction and the observed value normalized by the standard deviation of the observed data as the error metric. Furthermore, we compared the proposed method with two benchmark methods: most recent sample and an adaptive linear filter. The kernel density estimation-based prediction results demonstrate universally significant improvement over the alternatives and are especially valuable for long lookahead time, when the alternative methods fail to produce useful predictions.

Coherent bremsstrahlung denotes the process of bremsstrahlung emission by electrons traversing a single crystal, causing prominent peaks in the resulting photon spectra at low energies. While this phenomenon has been known for decades, little attempt has been made to exploit its potential for megavoltage imaging, where the quality of images is affected by low contrast due to the lack of sufficient photons at the energy range suited for diagnostic purposes. We provide a theoretical foundation of coherent bremsstrahlung in the first-order Born approximation without confinement to high energies. Based on this theory, first evidence is given that diamond crystals are capable of boosting the amount of diagnostic photons by about 10–20%. It is shown that this behaviour is largely conserved for polychromatic electrons hitting thick targets, where multiple scattering dominates the energy distribution of the emitted photons.

We discovered an empirical relationship between the logarithm of mean excitation energy (ln I_m) and the effective atomic number (EAN) of human tissues, which allows for computing patient-specific proton stopping power ratios (SPRs) using dual-energy CT (DECT) imaging. The accuracy of the DECT method was evaluated for 'standard' human tissues as well as their variance. The DECT method was compared to the existing standard clinical practice—a procedure introduced by Schneider et al at the Paul Scherrer Institute (the stoichiometric calibration method). In this simulation study, SPRs were derived from calculated CT numbers of known material compositions, rather than from measurement. For standard human tissues, both methods achieved good accuracy with the root-mean-square (RMS) error well below 1%. For human tissues with small perturbations from standard human tissue compositions, the DECT method was shown to be less sensitive than the stoichiometric calibration method. The RMS error remained below 1% for most cases using the DECT method, which implies that the DECT method might be more suitable for measuring patient-specific tissue compositions to improve the accuracy of treatment planning for charged particle therapy. In this study, the effects of CT imaging artifacts due to the beam hardening effect, scatter, noise, patient movement, etc were not analyzed. The true potential of the DECT method achieved in theoretical conditions may not be fully achievable in clinical settings. Further research and development may be needed to take advantage of the DECT method to characterize individual human tissues.

Increased pixel noise and streak artifact reduce CT image quality and limit the potential for radiation dose reduction during CT of the thoracic inlet. We propose to quantify the pixel noise of mediastinal structures in the thoracic inlet, during low-dose (LDCT) and ultralow-dose (uLDCT) thoracic CT, and assess the utility of new software (quantum denoising system and BOOST3D) in addressing these limitations. Twelve patients had LDCT (120 kV, 25 mAs) and uLDCT (120 kV, 10 mAs) images reconstructed initially using standard mediastinal and lung filters followed by the quantum denoising system (QDS) to reduce pixel noise and BOOST3D (B3D) software to correct photon starvation noise as follows: group 1 no QDS, no B3D; group 2 B3D alone; group 3 QDS alone and group 4 both QDS and B3D. Nine regions of interest (ROIs) were replicated on mediastinal anatomy in the thoracic inlet, for each patient resulting in 3456 data points to calculate pixel noise and attenuation. QDS reduced pixel noise by 18.4% (lung images) and 15.8% (mediastinal images) at 25 mAs. B3D reduced pixel noise by ∼8% in the posterior thorax and in combination there was a 35.5% reduction in effective radiation dose (E) for LDCT (1.63–1.05 mSv) in lung images and 32.2% (1.55–1.05 mSv) in mediastinal images. The same combination produced 20.7% reduction (0.53–0.42 mSv) in E for uLDCT, for lung images and 17.3% (0.51–0.42) for mediastinal images. This quantitative analysis of image quality confirms the utility of dedicated processing software in targeting image noise and streak artifact in thoracic LDCT and uLDCT images taken in the thoracic inlet. This processing software potentiates substantial reductions in radiation dose during thoracic LDCT and uLDCT.

The purpose of this study is to investigate the response of transverse relaxation time (T2) and apparent diffusion coefficient (ADC) in human glioma tumor xenografts during and after fractionated radiotherapy. Tumor-bearing mice were divided into four treatment groups (n = 6 per group) that received a total dose of 800 cGy of 200 kVp x-rays, given over two or three fractions, with a fraction spacing of either 24 or 72 h. A fifth treatment group received 800 cGy in a single fraction, and a sixth group of mice served as an untreated control. All mice were scanned pretreatment, before each fraction and at multiple points after treatment using a 9.4 T magnetic resonance imaging (MRI) system. Quantitative T2 and ADC maps were produced. All treated groups showed an increase in mean tumor ADC, though the time for this response to reach a maximum and return toward baseline was delayed in the fractionated groups. The highest ADC was measured 7 days after the final fraction of treatment for all groups. There were no significant differences in the maximum measured change in ADC between any of the treated groups, with the average measured maximum value being 20.5% above baseline. After treatment, all groups showed an increase in mean tumor T2, with the average measured maximum T2 being 4.7% above baseline. This increase was followed by a transition to mean T2 values below baseline values, with the average measured tumor T2 being 92.4% of the pretreatment value. The transition between elevated and depressed T2 values was delayed in the cases of fractionated therapies and occurred between 3.6 and 7.3 days after the last fraction of treatment. These results further the understanding of the temporal evolution of T2 and ADC during fractionated radiotherapy and support their potential use as time-sensitive biomarkers for tumor response.

Flow patterns may affect the potential of thrombus formation following plaque rupture. Computational fluid dynamics (CFD) were employed to assess hemodynamic conditions, and particularly flow recirculation and vortex formation in reconstructed arterial models associated with ST-elevation myocardial infraction (STEMI) or stable coronary stenosis (SCS) in the left anterior descending coronary artery (LAD). Results indicate that in the arterial models associated with STEMI, a 50% diameter stenosis immediately before or after a bifurcation creates a recirculation zone and vortex formation at the orifice of the bifurcation branch, for most of the cardiac cycle, thus allowing the creation of stagnating flow. These flow patterns are not seen in the SCS model with an identical stenosis. Post-stenotic recirculation in the presence of a 90% stenosis was evident at both the STEMI and SCS models. The presence of 90% diameter stenosis resulted in flow reduction in the LAD of 51.5% and 35.9% in the STEMI models and 37.6% in the SCS model, for a 10 mmHg pressure drop. CFD simulations in a reconstructed model of stenotic LAD segments indicate that specific anatomic characteristics create zones of vortices and flow recirculation that promote thrombus formation and potentially myocardial infarction.

The purpose of this paper is to assess the feasibility of delivering intensity- and energy-modulated electron radiation treatment (MERT) by a photon multileaf collimator (xMLC) and to evaluate the improvements obtained in shallow head and neck (HN) tumors. Four HN patient cases covering different clinical situations were planned by MERT, which used an in-house treatment planning system that utilized Monte Carlo dose calculation. The cases included one oronasal, two parotid and one middle ear tumors. The resulting dose-volume histograms were compared with those obtained from conventional photon and electron treatment techniques in our clinic, which included IMRT, electron beam and mixed beams, most of them using fixed-thickness bolus. Experimental verification was performed with plane-parallel ionization chambers for absolute dose verification, and a PTW ionization chamber array and radiochromic film for relative dosimetry. A MC-based treatment planning system for target with compromised volumes in depth and laterally has been validated. A quality assurance protocol for individual MERT plans was launched. Relative MC dose distributions showed a high agreement with film measurements and absolute ion chamber dose measurements performed at a reference point agreed with MC calculations within 2% in all cases. Clinically acceptable PTV coverage and organ-at-risk sparing were achieved by using the proposed MERT approach. MERT treatment plans, based on delivery of intensity-modulated electron beam using the xMLC, for superficial head and neck tumors, demonstrated comparable or improved PTV dose homogeneity with significantly lower dose to normal tissues. The clinical implementation of this technique will be able to offer a viable alternative for the treatment of shallow head and neck tumors.

In this study we compared two approaches that have recently been used to minimize precision errors in 3D quantitative computed tomography (QCT) images of the hip and the spine in order to optimize the detection of longitudinal changes in bone mineral density (BMD). In 30 subjects we obtained baseline and 1 year follow-up 3D CT scans of the proximal femur and the spine. QCT analysis was applied to a variety of volumes of interest (VOIs) automatically positioned relative to anatomic coordinate systems (ACS). In the first approach (A1) baseline and follow-up scans were analyzed independently. In the second approach (A2) a 3D versor-based rigid intensity registration method was applied to match baseline and follow-up images, and the baseline ACS was mapped on the follow-up image using the registration transformation. Afterwards, the analysis VOIs were again independently calculated for baseline and follow-up images. There were no significant differences of percent BMD changes between baseline and follow-up images between A1 and A2 for any of the VOIs investigated. With advanced image processing methods a time-consuming 3D registration between baseline and follow-up images before the analysis does not improve analysis precision compared to the use of anatomical coordinate systems.

The rapidly increasing number of diagnostic computed tomography (CT) procedures in the recent decades has spurred heightened concern over the potential risk to patients. Although an accurate organ dose assessment tool has now become highly desirable, existing software packages depend on stylized computational phantoms that were originally developed more than 40 years ago, exhibiting very large discrepancies when compared with phantoms that are anatomically realistic. However, past comparative studies did not focus on CT protocols for adult patients. This study was designed to quantitatively compare two types of phantoms, the stylized phantoms and a pair of recently developed RPI-adult male and adult female (RPI-AM and RPI-AF) phantoms, for various CT scanning protocols involving the chest, abdomen–pelvis and chest–abdomen–pelvis. Organ doses were based on Monte Carlo simulations using the MCNPX code and a detailed CT scanner model for the GE LightSpeed 16. Results are presented as ratios of organ doses from the stylized phantoms to those from the RPI phantoms. It is found that, for most organs contained in the scan volume, the ratios were within the range of 0.75–1.16. However, the stomach doses are significantly different and the ratio is found to be up to 1.86 in male phantoms and 2.29 in the female phantoms due to the anatomical differences between the two types of phantoms. Organs that lie near a scan boundary also exhibit a significant relative difference in organ doses between the two types of phantoms. This study concludes that, due to relatively low x-ray energies, CT doses are very sensitive to organ shape, size and position, and thus anatomically realistic phantoms should be used to avoid the dose uncertainties caused by the lack of anatomical realism. The new phantoms, such as the RPI-AM and AF phantoms that are designed using advanced surface meshes, are deformable and will make it possible to match the anatomy of a specific patient leading to further improvement in dose and risk assessments for patients undergoing CT examinations.

The addition of accurate system modeling in PET image reconstruction results in images with distinct noise texture and characteristics. In particular, the incorporation of point spread functions (PSF) into the system model has been shown to visually reduce image noise, but the noise properties have not been thoroughly studied. This work offers a systematic evaluation of noise and signal properties in different combinations of reconstruction methods and parameters. We evaluate two fully 3D PET reconstruction algorithms: (1) OSEM with exact scanner line of response modeled (OSEM+LOR), (2) OSEM with line of response and a measured point spread function incorporated (OSEM+LOR+PSF), in combination with the effects of four post-reconstruction filtering parameters and 1–10 iterations, representing a range of clinically acceptable settings. We used a modified NEMA image quality (IQ) phantom, which was filled with ⁶⁸Ge and consisted of six hot spheres of different sizes with a target/background ratio of 4:1. The phantom was scanned 50 times in 3D mode on a clinical system to provide independent noise realizations. Data were reconstructed with OSEM+LOR and OSEM+LOR+PSF using different reconstruction parameters, and our implementations of the algorithms match the vendor's product algorithms. With access to multiple realizations, background noise characteristics were quantified with four metrics. Image roughness and the standard deviation image measured the pixel-to-pixel variation; background variability and ensemble noise quantified the region-to-region variation. Image roughness is the image noise perceived when viewing an individual image. At matched iterations, the addition of PSF leads to images with less noise defined as image roughness (reduced by 35% for unfiltered data) and as the standard deviation image, while it has no effect on background variability or ensemble noise. In terms of signal to noise performance, PSF-based reconstruction has a 7% improvement in contrast recovery at matched ensemble noise levels and 20% improvement of quantitation SNR in unfiltered data. In addition, the relations between different metrics are studied. A linear correlation is observed between background variability and ensemble noise for all different combinations of reconstruction methods and parameters, suggesting that background variability is a reasonable surrogate for ensemble noise when multiple realizations of scans are not available.

The purpose of this work was to create a computational platform for studying motion in intensity modulated radiotherapy (IMRT). Specifically, the non-uniform rational B-spline (NURB) cardiac and torso (NCAT) phantom was modified for use in a four-dimensional Monte Carlo (4D-MC) simulation system to investigate the effect of respiratory-induced intra-fraction organ motion on IMRT dose distributions as a function of diaphragm motion, lesion size and lung density. Treatment plans for four clinical scenarios were designed: diaphragm peak-to-peak amplitude of 1 cm and 3 cm, and two lesion sizes—2 cm and 4 cm diameter placed in the lower lobe of the right lung. Lung density was changed for each phase using a conservation of mass calculation. Further, a new heterogeneous lung model was implemented and tested. Each lesion had an internal target volume (ITV) subsequently expanded by 15 mm isotropically to give the planning target volume (PTV). The PTV was prescribed to receive 72 Gy in 40 fractions. The MLC leaf sequence defined by the planning system for each patient was exported and used as input into the MC system. MC simulations using the dose planning method (DPM) code together with deformable image registration based on the NCAT deformation field were used to find a composite dose distribution for each phantom. These composite distributions were subsequently analyzed using information from the dose volume histograms (DVH). Lesion motion amplitude has the largest effect on the dose distribution. Tumor size was found to have a smaller effect and can be mitigated by ensuring the planning constraints are optimized for the tumor size. The use of a dynamic or heterogeneous lung density model over a respiratory cycle does not appear to be an important factor with a ≤ 0.6% change in the mean dose received by the ITV, PTV and right lung. The heterogeneous model increases the realism of the NCAT phantom and may provide more accurate simulations in radiation therapy investigations that use the phantom. This work further evaluates the NCAT phantom for use as a tool in radiation therapy research in addition to its extensive use in diagnostic imaging and nuclear medicine research. Our results indicate that the NCAT phantom, combined with 4D-MC simulations, is a useful tool in radiation therapy investigations and may allow the study of relative effects in many clinically relevant situations.

Biomagnetic instruments usually employ sensors with approximately radial normal vectors arranged on a near-spherical surface. The multipole expansion employed in the recently introduced signal space separation method (SSS) reveals limitations in this traditional sensor array design. Specifically, we show that the excessive symmetry of the sensor array impedes separation of multipole components arising from inside and outside of the array. This motivates consideration of novel instrument designs that would sample the field in a more comprehensive way. We evaluated several simulated sensor arrays that employ vector sensors in one or two layers, giving information on multiple field components and the radial dependence of the field. Our results indicate that this kind of sensor array design could significantly improve SSS performance, leading to enhanced shielding against external interference and reduced noise after signal reconstruction. The best two-layer array evaluated here attains a shielding factor of nearly 1000 or 60 dB with about 400 sensors. Due to limited spatial coverage, a traditional reference array geometry does not give the same level of improvement. In addition to improved software shielding, enhanced detection of different multipole components increases the information obtained about the magnetic field, which has fundamental importance.

Parametric imaging using dynamic positron emission tomography (PET) provides important information for biological research and clinical diagnosis. Indirect and direct methods have been developed for reconstructing linear parametric images from dynamic PET data. Indirect methods are relatively simple and easy to implement because the image reconstruction and kinetic modeling are performed in two separate steps. Direct methods estimate parametric images directly from raw PET data and are statistically more efficient. However, the convergence rate of direct algorithms can be slow due to the coupling between the reconstruction and kinetic modeling. Here we present two fast gradient-type algorithms for direct reconstruction of linear parametric images. The new algorithms decouple the reconstruction and linear parametric modeling at each iteration by employing the principle of optimization transfer. Convergence speed is accelerated by running more sub-iterations of linear parametric estimation because the computation cost of the linear parametric modeling is much less than that of the image reconstruction. Computer simulation studies demonstrated that the new algorithms converge much faster than the traditional expectation maximization (EM) and the preconditioned conjugate gradient algorithms for dynamic PET.

This paper presents finite-difference time-domain (FDTD) calculations of SAR in the University of Florida newborn female model. The newborn model is based upon a surface representation of the organs of the body, using non-uniform rational B-spline surfaces (NURBS). The surface model can then be converted into voxels at any resolution required. This flexibility allows the preparation of voxel models at 2, 1 and 0.5 mm to investigate the effect of resolution on dispersion and the choice of algorithms to calculate SAR in the Yee cell as the frequency increases up to 6 GHz. The added advantage of the newborn model is that it is relatively small and so FDTD calculations can be made tractable at a very fine resolution of 0.5 mm. A comparison is made between the calculated external electric fields required to produce the basic restriction on whole-body-averaged SAR and the ICNIRP reference levels for public exposure. At 250 MHz, the whole body resonance, the ICNIRP reference level does not provide a conservative estimate of the whole-body-averaged SAR restriction. The reference level is also breached in the range 700–2450 MHz by all of the irradiation geometries considered.

X-ray microbeam radiation therapy (MRT) is a novel form of treatment, currently in its preclinical stage, which uses microplanar x-ray beams from a synchrotron radiation source. It is important to perform accurate dosimetry on these microbeams, but, to date, there has been no accurate enough method available for making 3D dose measurements with isotropic, high spatial resolution to verify the results of Monte Carlo dose simulations. Here, we investigate the potential of optical computed tomography for satisfying these requirements. The construction of a simple optical CT microscopy (optical projection tomography) system from standard commercially available hardware is described. The measurement of optical densities in projection data is shown to be highly linear (r² = 0.999). The depth-of-field (DOF) of the imaging system is calculated based on the previous literature and measured experimentally using a commercial DOF target. It is shown that high quality images can be acquired despite the evident lack of telecentricity and despite DOF of the system being much lower than the sample diameter. Possible reasons for this are discussed. Results are presented for a complex irradiation of a 22 mm diameter cylinder of the radiochromic polymer PRESAGE™, demonstrating the exquisite 'dose-painting' abilities available in the MRT hutch of beamline ID-17 at the European Synchrotron Radiation Facility. Dose distributions in this initial experiment are equally well resolved on both an optical CT scan and a corresponding transmission image of radiochromic film, down to a line width of 83 µm (6 lp mm⁻¹) with an MTF value of 0.40. A group of 33 µm wide lines was poorly resolved on both the optical CT and film images, and this is attributed to an incorrect exposure time calculation, leading to under-delivery of dose. Image artefacts in the optical CT scan are discussed. PRESAGE™ irradiated using the microbeam facility is proposed as a suitable material for producing phantom samples for quantitative characterization of optical CT microscopy systems.

Online monitoring of the stopping distribution of particle beams used for radiotherapy provides the possibility of detecting possible errors in dose deposition early during a given treatment session, and may therefore help to improve the quality of the therapy. Antiproton annihilation events produce several long-range secondary particles which can be detected in real time by standard high energy particle physics detector systems. In this note, Monte Carlo calculations are performed in order to study the feasibility of real-time imaging by detecting charged pions produced during antiproton irradiation of typical biological targets. A simple treatment plan in a water phantom is simulated and the results show that by detecting π^± the position and the size of the planned target volume can be located with precision in the order of 1 mm.

In this study, we demonstrate the proof of principle of the University Medical Center Utrecht (UMCU) robot dedicated to magnetic resonance imaging (MRI)-guided interventions in patients. The UMCU robot consists of polymers and non-ferromagnetic materials. For transperineal prostate interventions, it can be placed between the patient's legs inside a closed bore 1.5T MR scanner. The robot can manually be translated and rotated resulting in five degrees of freedom. It contains a pneumatically driven tapping device to automatically insert a needle stepwise into the prostate using a controller unit outside the scanning room. To define the target positions and to verify the needle insertion point and the needle trajectory, a high-resolution 3D balanced steady state free precession (bSSFP) scan that provides a T2/T1-weighted contrast is acquired. During the needle insertion fast 2D bSSFP images are generated to track the needle on-line. When the target position is reached, the radiation oncologist manually places a fiducial gold marker (small seed) at this location. In total two needle trajectories are used to place all markers. Afterwards, a high-resolution 3D bSSFP scan is acquired to visualize the fiducial gold markers. Four fiducial gold markers were placed transperineally into the prostate of a patient with a clinical stage T3 prostate cancer. In the generated scans, it was possible to discriminate the patient's anatomy, the needle and the markers. All markers were delivered inside the prostate. The procedure time was 1.5 h. This study proves that MRI-guided needle placement and seed delivery in the prostate with the UMCU robot are feasible.

The authors of a recent paper (Wang and Rogers 2009 Phys. Med. Biol.54 1609) have used the Monte Carlo method to simulate the 'classical' experiment made more than 30 years ago by Johansson et al (1978 National and International Standardization of Radiation Dosimetry (Atlanta 1977) vol 2 (Vienna: IAEA) pp 243–70) on the displacement (or replacement) perturbation correction factor p_dis for cylindrical chambers in ⁶⁰Co and high-energy photon beams. They conclude that an 'unreasonable normalization at d_max' of the ionization chambers response led to incorrect results, and for the IAEA TRS-398 Code of Practice, which uses ratios of those results, 'the difference in the correction factors can lead to a beam calibration deviation of more than 0.5% for Farmer-like chambers'. The present work critically examines and questions some of the claims and generalized conclusions of the paper. It is demonstrated that for real, commercial Farmer-like chambers, the possible deviations in absorbed dose would be much smaller (typically 0.13%) than those stated by Wang and Rogers, making the impact of their proposed values negligible on practical high-energy photon dosimetry. Differences of the order of 0.4% would only appear at the upper extreme of the energies potentially available for clinical use (around 25 MV) and, because lower energies are more frequently used, the number of radiotherapy photon beams for which the deviations would be larger than say 0.2% is extremely small. This work also raises concerns on the proposed value of p_dis for Farmer chambers at the reference quality of ⁶⁰Co in relation to their impact on electron beam dosimetry, both for direct dose determination using these chambers and for the cross-calibration of plane-parallel chambers. The proposed increase of about 1% in p_dis (compared with TRS-398) would lower the k_Q factors and therefore D_w in electron beams by the same amount. This would yield a severe discrepancy with the current good agreement between electron dosimetry based on an electron cross-calibrated plane-parallel chamber (against a Farmer) or on a directly ⁶⁰Co calibrated plane-parallel chamber, which is not likely to be in error by 1%. It is suggested that the influence of the ⁶⁰Co source spectrum used in the simulations may not be negligible for calculations aimed at an uncertainty level of 0.1%.

In this letter we reply to the letter by Professor Andreo (Andreo 2010 Phys. Med. Biol.55 L9–16) regarding the values of the replacement correction factor for cylindrical chambers, which in turn commented on our recent paper on this (Wang and Rogers 2009 Phys. Med. Biol.54 1609–20).