Lung tumor motion change during stereotactic body radiotherapy (SBRT): an evaluation using MRI.

The purpose of this study is to investigate changes in lung tumor internal target volume during stereotactic body radiotherapy treatment (SBRT) using magnetic resonance imaging (MRI). Ten lung cancer patients (13 tumors) undergoing SBRT (48 Gy over four consecutive days) were evaluated. Each patient underwent three lung MRI evaluations: before SBRT (MRI-1), after fraction 3 of SBRT (MRI-3), and three months after completion of SBRT (MRI-3m). Each MRI consisted of T1-weighted images in axial plane through the entire lung. A cone-beam CT (CBCT) was taken before each fraction. On MRI and CBCT taken before fractions 1 and 3, gross tumor volume (GTV) was contoured and differences between the two volumes were compared. Median tumor size on CBCT before fractions1 (CBCT-1) and 3 (CBCT-3) was 8.68 and 11.10 cm3, respectively. In 12 tumors, the GTV was larger on CBCT-3 compared to CBCT-1 (median enlargement, 1.56 cm3). Median tumor size on MRI-1, MRI-3, and MRI-3m was 7.91, 11.60, and 3.33 cm3, respectively. In all patients, the GTV was larger on MRI-3 compared to MRI-1 (median enlargement, 1.54 cm3). In all patients, GTV was smaller on MRI-3m compared to MRI-1 (median shrinkage, 5.44 cm3). On CBCT and MRI, all patients showed enlargement of the GTV during the treatment week of SBRT, except for one patient who showed minimal shrinkage (0.86 cm3). Changes in tumor volume are unpredictable; therefore, motion and breathing must be taken into account during treatment planning, and image-guided methods should be used, when treating with large fraction sizes.

Outcome measures for use in trials of paediatric otorrhoea: A systematic review.

INTRODUCTION: Paediatric otorrhoea (PO) describes a middle ear infection that results in a perforation of the tympanic membrane and ear discharge, in children and young people (CYP). Prolonged infection may be associated with hearing loss and developmental delay. The current management of paediatric otorrhoea is variable, including non-invasive treatments (conservative, oral antibiotics, topical antibiotics) and surgery, reflecting the lack of a sufficiently strong evidence base. Outcome reporting is fundamental to producing reliable and meaningful evidence to inform best practice. OBJECTIVES: Primary objective: to determine which outcome measures are currently used to evaluate treatment success in studies of non-surgical treatments for paediatric otorrhoea. SECONDARY OBJECTIVES: to identify outcome measurement instruments used in the literature and assess their applicability for use in clinical trials of PO. METHODS: This systematic review was registered with PROSPERO (CRD42023407976). Database searches of EMBASE, MEDLINE and Cochrane was performed on June 6, 2023, covering from Jan 1995 to May 2023. Randomised controlled trials or study protocols involving CYP with PO were included following PRISMA guidelines. Risk of bias was assessed with Cochrane's tool. RESULTS: Of the 377 papers identified, six were included in the systematic review. The primary outcome of five of the studies related to otorrhoea cessation; both time to cessation and proportion recovered at various time points were used as measures. Two measurement instruments were identified: Otitis Media-6 Questionnaire and the Institute for Medical Technology Assessment Productivity Cost Questionnaire. Both were shown to be applicable measurement instruments when used in clinical trials of PO. CONCLUSIONS: To promote homogeneity and facilitate meaningful comparison and combination of studies, we propose that time to cessation of otorrhoea from onset of otorrhoea should be used as the primary outcome in future studies. Further research is needed to establish if this is the most important outcome to children and their caregivers.

Protocol for the Paediatric Otorrhoea Study (POSt): a multi-methods study to understand the burden of paediatric otorrhoea in the UK.

INTRODUCTION: Paediatric otorrhoea (PO) refers to the leakage of fluid through a perforation in the ear drum, resulting from an infection of the middle ear of a child or young person (CYP). PO frequently results in hearing loss which may lead to developmental delay, restricted communication and reduced educational attainment.Epidemiological information for PO is largely derived from low-income countries. The aim of this study will be to establish the incidence of PO within the UK and to understand the impact of PO on CYP and their families' everyday lives. It will build the foundations for a randomised controlled trial investigating the best antibiotic treatment for PO. METHODS AND ANALYSIS: The study will consist of two work packages. (1) Data from the Clinical Practice Research Datalink (CPRD), January 2005 to July 2021, will be used to determine the incidence of patient presentations with PO to primary care in the UK. It will also explore the current antimicrobial prescribing practice for PO in primary care. (2) Thirty semi-structured interviews will be conducted from 13 July to 31 October 2023 with CYP and their parents/carers to help identify the impact of PO on everyday life, the patient journey and how service users define treatment success. Three medical professional focus groups will be used to understand the current management practice, how treatment success is measured and acceptability to randomise patients. Thematic analysis will be used. ETHICS AND DISSEMINATION: The Health Research Authority, The Health and Social Care Research Ethics Committee (23/NI/0082) and the CPRD's research data governance panel (22_002508) reviewed this study. Results will be disseminated at medical conferences, in peer-reviewed journals and via social media. The study will cocreate a webpage on healthtalk.org, with the Dipex Charity, about PO to ensure members of the public can learn more about the condition. TRIAL REGISTRATION NUMBER: ISRCTN46071200.

Dosimetry tools and techniques for IMRT.

Intensity modulated radiation therapy (IMRT) poses a number of challenges for properly measuring commissioning data and quality assurance (QA) radiation dose distributions. This report provides a comprehensive overview of how dosimeters, phantoms, and dose distribution analysis techniques should be used to support the commissioning and quality assurance requirements of an IMRT program. The proper applications of each dosimeter are described along with the limitations of each system. Point detectors, arrays, film, and electronic portal imagers are discussed with respect to their proper use, along with potential applications of 3D dosimetry. Regardless of the IMRT technique utilized, some situations require the use of multiple detectors for the acquisition of accurate commissioning data. The overall goal of this task group report is to provide a document that aids the physicist in the proper selection and use of the dosimetry tools available for IMRT QA and to provide a resource for physicists that describes dosimetry measurement techniques for purposes of IMRT commissioning and measurement-based characterization or verification of IMRT treatment plans. This report is not intended to provide a comprehensive review of commissioning and QA procedures for IMRT. Instead, this report focuses on the aspects of metrology, particularly the practical aspects of measurements that are unique to IMRT. The metrology of IMRT concerns the application of measurement instruments and their suitability, calibration, and quality control of measurements. Each of the dosimetry measurement tools has limitations that need to be considered when incorporating them into a commissioning process or a comprehensive QA program. For example, routine quality assurance procedures require the use of robust field dosimetry systems. These often exhibit limitations with respect to spatial resolution or energy response and need to themselves be commissioned against more established dosimeters. A chain of dosimeters, from secondary standards to field instruments, is established to assure the quantitative nature of the tests. This report is intended to describe the characteristics of the components of these systems; dosimeters, phantoms, and dose evaluation algorithms. This work is the report of AAPM Task Group 120.

Assessment of the setup dependence of detector response functions for mega-voltage linear accelerators.

PURPOSE: Accurate modeling of beam profiles is important for precise treatment planning dosimetry. Calculated beam profiles need to precisely replicate profiles measured during machine commissioning. Finite detector size introduces perturbations into the measured profiles, which, in turn, impact the resulting modeled profiles. The authors investigate a method for extracting the unperturbed beam profiles from those measured during linear accelerator commissioning. METHODS: In-plane and cross-plane data were collected for an Elekta Synergy linac at 6 MV using ionization chambers of volume 0.01, 0.04, 0.13, and 0.65 cm3 and a diode of surface area 0.64 mm2. The detectors were orientated with the stem perpendicular to the beam and pointing away from the gantry. Profiles were measured for a 10 x 10 cm2 field at depths ranging from 0.8 to 25.0 cm and SSDs from 90 to 110 cm. Shaping parameters of a Gaussian response function were obtained relative to the Edge detector. The Gaussian function was deconvolved from the measured ionization chamber data. The Edge detector profile was taken as an approximation to the true profile, to which deconvolved data were compared. Data were also collected with CC13 and Edge detectors for additional fields and energies on an Elekta Synergy, Varian Trilogy, and Siemens Oncor linear accelerator and response functions obtained. Response functions were compared as a function of depth, SSD, and detector scan direction. Variations in the shaping parameter were introduced and the effect on the resulting deconvolution profiles assessed. RESULTS: Up to 10% setup dependence in the Gaussian shaping parameter occurred, for each detector for a particular plane. This translated to less than a +/- 0.7 mm variation in the 80%-20% penumbral width. For large volume ionization chambers such as the FC65 Farmer type, where the cavity length to diameter ratio is far from 1, the scan direction produced up to a 40% difference in the shaping parameter between in-plane and cross-plane measurements. This is primarily due to the directional difference in penumbral width measured by the FC65 chamber, which can more than double in profiles obtained with the detector stem parallel compared to perpendicular to the scan direction. For the more symmetric CC13 chamber the variation was only 3% between in-plane and cross-plane measurements. CONCLUSIONS: The authors have shown that the detector response varies with detector type, depth, SSD, and detector scan direction. In-plane vs. cross-plane scanning can require calculation of a direction dependent response function. The effect of a 10% overall variation in the response function, for an ionization chamber, translates to a small deviation in the penumbra from that of the Edge detector measured profile when deconvolved. Due to the uncertainties introduced by deconvolution the Edge detector would be preferable in obtaining an approximation of the true profile, particularly for field sizes where the energy dependence of the diode can be neglected. However, an averaged response function could be utilized to provide a good approximation of the true profile for large ionization chambers and for larger fields for which diode detectors are not recommended.

Accurate heterogeneous dose calculation for lung cancer patients without high-resolution CT densities.

The aim of this study was to investigate the relative accuracy of megavoltage photon-beam dose calculations employing either 5 bulk densities or independent voxel densities determined by calibration of the CT Houndsfield number. Full-resolution CT and bulk density treatment plans were generated for 70 lung or esophageal cancer tumors (66 cases) using a commercial treatment planning system with an adaptive convolution dose calculation algorithm (Pinnacle3, Philips Medicals Systems). Bulk densities were applied to segmented regions. Individual and population average densities were compared to the full-resolution plan for each case. Monitor units were kept constant and no normalizations were employed. Dose volume histograms (DVH) and dose difference distributions were examined for all cases. The average densities of the segmented air, lung, fat, soft tissue, and bone for the entire set were found to be 0.14, 0.26, 0.89, 1.02, and 1.12 g/cc, respectively. In all cases, the normal tissue DVH agreed to better than 2% in dose. In 62 of 70 target DVHs, agreement to better than 3% in dose was observed. Six cases demonstrated emphysema, one with bullous formations and one with a hiatus hernia having a large volume of gas. They required the additional assignment of density to the emphysemic lung and inflammatory changes to the lung, the regions of collapsed lung, the bullous formations, and the hernia gas. Bulk tissue density dose calculation provides an accurate method of heterogeneous dose calculation. However, patients with advanced emphysema may require high-resolution CT studies for accurate treatment planning.

Evaluating target cold spots by the use of tail EUDs.

PURPOSE: To propose a new measure of target underdose that can be used in the evaluation and optimization of radiotherapy dose distributions. METHODS AND MATERIALS: We compare various formulations of the equivalent uniform dose (EUD) and introduce a modification of existing EUD definitions, which we call tail EUD. Tail EUD is a measure of "cold spots" below the prescription dose in the target dose distribution, using units of gray (Gy). We investigate the mathematical properties of various target EUD concepts, including tail EUD. We apply the tail EUD measure retrospectively to intensity modulated radiation therapy (IMRT) treatment plans from our plan database. We also use tail EUD as an optimization objective in the optimization of prostate, pancreas, and head-and-neck plans. RESULTS: Tail EUD has desirable mathematical properties. In particular, it is convex and it leads to convex level sets (i.e., no local minima) if the EUD from which it is derived is concave. The tail EUD value is correlated with the subjective degree of target coverage. Constraining tail EUDs to a certain level in plan optimization leads to comparable target coverage in different plans and treatment sites. CONCLUSIONS: The newly introduced concept of tail EUD appears to be useful for both plan evaluation and optimization. In addition it can potentially be applied in the design of new clinical protocols.

A prototype quantitative film scanner for radiochromic film dosimetry.

We have developed a high resolution, quantitative, two-dimensional optical film scanner for use with a commercial high sensitivity radiochromic film (RCF) for measuring single fraction external-beam radiotherapy dose distributions. The film scanner was designed to eliminate artifacts commonly observed in RCF dosimetry. The scanner employed a stationary light source and detector with a moving antireflective glass film platen attached to a high precision computerized X-Y translation stage. An ultrabright red light emitting diode (LED) with a peak output at 633 nm and full width at half maximum (FWHM) of 16 nm was selected as the scanner light source to match the RCF absorption peak. A dual detector system was created using two silicon photodiode detectors to simultaneously measure incident and transmitted light. The LED light output was focused to a submillimeter (FWHM 0.67 mm) spot size, which was determined from a scanning knife-edge technique for measuring Gaussian optical beams. Data acquisition was performed with a 16-bit A/D card in conjunction with commercial software. The linearity of the measured densities on the scanner was tested using a calibrated neutral-density step filter. Sensitometric curves and three IMRT field scans were acquired with a spatial resolution of 1 mm for both radiographic film and RCF. The results were compared with measurements taken with a commercial diode array under identical delivery conditions. The RCF was rotated by 90 deg and rescanned to study orientation effects. Comparison between the RCF and the diode array measurements using percent dose difference and distance-to-agreement criteria produced average passing rates of 99.0% using 3%/3 mm criteria and 96.7% using 2%/2 mm criteria. The same comparison between the radiographic film and diode array measurements resulted in average passing rates 96.6% and 91.6% for the above two criteria, respectively. No measurable light-scatter or interference scanner artifacts were observed. The RCF rotated by 90 deg showed no measurable orientation effect. A scan of a 15 x 15 cm2 area with 1 mm resolution required 22 min to acquire. The LED densitometer provides an accurate film dosimetry system with 1 mm or better resolution. The scanner eliminates the orientation dependence of RCF dosimetry that was previously reported with commercial flatbed scanners.

A computational implementation and comparison of several intensity modulated proton therapy treatment planning algorithms.

The authors present a comparative study of intensity modulated proton therapy (IMPT) treatment planning employing algorithms of three-dimensional (3D) modulation, and 2.5-dimensional (2.5D) modulation, and intensity modulated distal edge tracking (DET) [A. Lomax, Phys. Med. Biol. 44, 185-205 (1999)] applied to the treatment of head-and-neck cancer radiotherapy. These three approaches were also compared with 6 MV photon intensity modulated radiation therapy (IMRT). All algorithms were implemented in the University of Florida Optimized Radiation Therapy system using a finite sized pencil beam dose model and a convex fluence map optimization model. The 3D IMPT and the DET algorithms showed considerable advantages over the photon IMRT in terms of dose conformity and sparing of organs at risk when the beam number was not constrained. The 2.5D algorithm did not show an advantage over the photon IMRT except in the dose reduction to the distant healthy tissues, which is inherent in proton beam delivery. The influences of proton beam number and pencil beam size on the IMPT plan quality were also studied. Out of 24 cases studied, three cases could be adequately planned with one beam and 12 cases could be adequately planned with two beams, but the dose uniformity was often marginally acceptable. Adding one or two more beams in each case dramatically improved the dose uniformity. The finite pencil beam size had more influence on the plan quality of the 2.5D and DET algorithms than that of the 3D IMPT. To obtain a satisfactory plan quality, a 0.5 cm pencil beam size was required for the 3D IMPT and a 0.3 cm size was required for the 2.5D and the DET algorithms. Delivery of the IMPT plans produced in this study would require a proton beam spot scanning technique that has yet to be developed clinically.

Comparative analysis of 60Co intensity-modulated radiation therapy.

In this study, we perform a scientific comparative analysis of using (60)Co beams in intensity-modulated radiation therapy (IMRT). In particular, we evaluate the treatment plan quality obtained with (i) 6 MV, 18 MV and (60)Co IMRT; (ii) different numbers of static multileaf collimator (MLC) delivered (60)Co beams and (iii) a helical tomotherapy (60)Co beam geometry. We employ a convex fluence map optimization (FMO) model, which allows for the comparison of plan quality between different beam energies and configurations for a given case. A total of 25 clinical patient cases that each contain volumetric CT studies, primary and secondary delineated targets, and contoured structures were studied: 5 head-and-neck (H&N), 5 prostate, 5 central nervous system (CNS), 5 breast and 5 lung cases. The DICOM plan data were anonymized and exported to the University of Florida optimized radiation therapy (UFORT) treatment planning system. The FMO problem was solved for each case for 5-71 equidistant beams as well as a helical geometry for H&N, prostate, CNS and lung cases, and for 3-7 equidistant beams in the upper hemisphere for breast cases, all with 6 MV, 18 MV and (60)Co dose models. In all cases, 95% of the target volumes received at least the prescribed dose with clinical sparing criteria for critical organs being met for all structures that were not wholly or partially contained within the target volume. Improvements in critical organ sparing were found with an increasing number of equidistant (60)Co beams, yet were marginal above 9 beams for H&N, prostate, CNS and lung. Breast cases produced similar plans for 3-7 beams. A helical (60)Co beam geometry achieved similar plan quality as static plans with 11 equidistant (60)Co beams. Furthermore, 18 MV plans were initially found not to provide the same target coverage as 6 MV and (60)Co plans; however, adjusting the trade-offs in the optimization model allowed equivalent target coverage for 18 MV. For plans with comparable target coverage, critical structure sparing was best achieved with 6 MV beams followed closely by (60)Co beams, with 18 MV beams requiring significantly increased dose to critical structures. In this paper, we report in detail on a representative set of results from these experiments. The results of the investigation demonstrate the potential for IMRT radiotherapy employing commercially available (60)Co sources and a double-focused MLC. Increasing the number of equidistant beams beyond 9 was not observed to significantly improve target coverage or critical organ sparing and static plans were found to produce comparable plans to those obtained using a helical tomotherapy treatment delivery when optimized using the same well-tuned convex FMO model. While previous studies have shown that 18 MV plans are equivalent to 6 MV for prostate IMRT, we found that the 18 MV beams actually required more fluence to provide similar quality target coverage.

Performance of bilingual speakers on the English and Spanish versions of the Hearing in Noise Test (HINT).

This study compared the performance of bilingual participants on the English and Spanish versions of the Hearing in Noise Test (HINT). The participants were divided into an early bilingual (EB) group and a late bilingual (LB) group based on age of second-language acquisition. All participants acquired Spanish as their first language (L1) and English as a second language (L2). Care was taken to ensure that all participants demonstrated at least a "good competence level" for self-rated speaking, understanding, reading, and writing skills in both English and Spanish. Results revealed superior performance on the Spanish HINT versus the English HINT in both quiet and in noise for both groups of participants. Significant differences in performance were noted for the EB versus the LB participants. A number of possible explanations for superior performance in L1 are provided, and implications for educating students in their L2 are discussed.

TG-69: radiographic film for megavoltage beam dosimetry.

TG-69 is a task group report of the AAPM on the use of radiographic film for dosimetry. Radiographic films have been used for radiation dosimetry since the discovery of x-rays and have become an integral part of dose verification for both routine quality assurance and for complex treatments such as soft wedges (dynamic and virtual), intensity modulated radiation therapy (IMRT), image guided radiation therapy (IGRT), and small field dosimetry like stereotactic radiosurgery. Film is convenient to use, spatially accurate, and provides a permanent record of the integrated two dimensional dose distributions. However, there are several challenges to obtaining high quality dosimetric results with film, namely, the dependence of optical density on photon energy, field size, depth, film batch sensitivity differences, film orientation, processing conditions, and scanner performance. Prior to the clinical implementation of a film dosimetry program, the film, processor, and scanner need to be tested to characterize them with respect to these variables. Also, the physicist must understand the basic characteristics of all components of film dosimetry systems. The primary mission of this task group report is to provide guidelines for film selection, irradiation, processing, scanning, and interpretation to allow the physicist to accurately and precisely measure dose with film. Additionally, we present the basic principles and characteristics of film, processors, and scanners. Procedural recommendations are made for each of the steps required for film dosimetry and guidance is given regarding expected levels of accuracy. Finally, some clinical applications of film dosimetry are discussed.

Optimization of equivalent uniform dose using the L-curve criterion.

Optimization of equivalent uniform dose (EUD) in inverse planning for intensity-modulated radiation therapy (IMRT) prevents variation in radiobiological effect between different radiotherapy treatment plans, which is due to variation in the pattern of dose nonuniformity. For instance, the survival fraction of clonogens would be consistent with the prescription when the optimized EUD is equal to the prescribed EUD. One of the problems in the practical implementation of this approach is that the spatial dose distribution in EUD-based inverse planning would be underdetermined because an unlimited number of nonuniform dose distributions can be computed for a prescribed value of EUD. Together with ill-posedness of the underlying integral equation, this may significantly increase the dose nonuniformity. To optimize EUD and keep dose nonuniformity within reasonable limits, we implemented into an EUD-based objective function an additional criterion which ensures the smoothness of beam intensity functions. This approach is similar to the variational regularization technique which was previously studied for the dose-based least-squares optimization. We show that the variational regularization together with the L-curve criterion for the regularization parameter can significantly reduce dose nonuniformity in EUD-based inverse planning.

An exact approach to direct aperture optimization in IMRT treatment planning.

We consider the problem of intensity-modulated radiation therapy (IMRT) treatment planning using direct aperture optimization. While this problem has been relatively well studied in recent years, most approaches employ a heuristic approach to the generation of apertures. In contrast, we use an exact approach that explicitly formulates the fluence map optimization (FMO) problem as a convex optimization problem in terms of all multileaf collimator (MLC) deliverable apertures and their associated intensities. However, the number of deliverable apertures, and therefore the number of decision variables and constraints in the new problem formulation, is typically enormous. To overcome this, we use an iterative approach that employs a subproblem whose optimal solution either provides a suitable aperture to add to a given pool of allowable apertures or concludes that the current solution is optimal. We are able to handle standard consecutiveness, interdigitation and connectedness constraints that may be imposed by the particular MLC system used, as well as jaws-only delivery. Our approach has the additional advantage that it can explicitly account for transmission of dose through the part of an aperture that is blocked by the MLC system, yielding a more precise assessment of the treatment plan than what is possible using a traditional beamlet-based FMO problem. Finally, we develop and test two stopping rules that can be used to identify treatment plans of high clinical quality that are deliverable very efficiently. Tests on clinical head-and-neck cancer cases showed the efficacy of our approach, yielding treatment plans comparable in quality to plans obtained by the traditional method with a reduction of more than 75% in the number of apertures and a reduction of more than 50% in beam-on time, with only a modest increase in computational effort. The results also show that delivery efficiency is very insensitive to the addition of traditional MLC constraints; however, jaws-only treatment requires about a doubling in beam-on time and number of apertures used. Finally, we showed the importance of accounting for transmission effects when assessing or, preferably, optimizing treatment plan quality.

The impact of heterogeneity correction on dosimetric parameters that predict for radiation pneumonitis.

PURPOSE: To determine if heterogeneity correction significantly affects commonly measured dosimetric parameters predicting pulmonary toxicity in patients receiving radiation for lung cancer. METHODS AND MATERIALS: Sixty-eight patients treated for lung cancer were evaluated. The conformal treatment technique mostly employed anteroposterior/posterior-anterior fields and off-cord obliques. The percent total lung volume receiving 20 Gy or higher (V20) and mean lung dose (MLD) were correlated with the incidence of radiation pneumonitis. Parameters from both heterogeneity-corrected and heterogeneity-uncorrected plans were used to assess this risk. RESULTS: Univariate analysis revealed a significant correlation between the development of radiation pneumonitis and both V20 and MLD. A best-fit line to a plot of V20 from the homogeneous plan against the corresponding V20 heterogeneous value produced a slope of 1.00 and zero offset, indicating no difference between the two parameters. For MLD, a similarly significant correlation is seen between the heterogeneous and homogeneous parameters, indicating a 4% difference when correcting for heterogeneity. A significant correlation was also observed between the MLD and V20 parameters (p < 0.0001). CONCLUSIONS: A high degree of correlation exists between heterogeneity-corrected and heterogeneity-uncorrected dosimetric parameters for lung and the risk of developing pneumonitis. Either V20 or MLD predicts the pneumonitis risk with similar effect.

A fourier analysis on the maximum acceptable grid size for discrete proton beam dose calculation.

We developed an analytical method for determining the maximum acceptable grid size for discrete dose calculation in proton therapy treatment plan optimization, so that the accuracy of the optimized dose distribution is guaranteed in the phase of dose sampling and the superfluous computational work is avoided. The accuracy of dose sampling was judged by the criterion that the continuous dose distribution could be reconstructed from the discrete dose within a 2% error limit. To keep the error caused by the discrete dose sampling under a 2% limit, the dose grid size cannot exceed a maximum acceptable value. The method was based on Fourier analysis and the Shannon-Nyquist sampling theorem as an extension of our previous analysis for photon beam intensity modulated radiation therapy [J. F. Dempsey, H. E. Romeijn, J. G. Li, D. A. Low, and J. R. Palta, Med. Phys. 32, 380-388 (2005)]. The proton beam model used for the analysis was a near monoenergetic (of width about 1% the incident energy) and monodirectional infinitesimal (nonintegrated) pencil beam in water medium. By monodirection, we mean that the proton particles are in the same direction before entering the water medium and the various scattering prior to entrance to water is not taken into account. In intensity modulated proton therapy, the elementary intensity modulation entity for proton therapy is either an infinitesimal or finite sized beamlet. Since a finite sized beamlet is the superposition of infinitesimal pencil beams, the result of the maximum acceptable grid size obtained with infinitesimal pencil beam also applies to finite sized beamlet. The analytic Bragg curve function proposed by Bortfeld [T. Bortfeld, Med. Phys. 24, 2024-2033 (1997)] was employed. The lateral profile was approximated by a depth dependent Gaussian distribution. The model included the spreads of the Bragg peak and the lateral profiles due to multiple Coulomb scattering. The dependence of the maximum acceptable dose grid size on the orientation of the beam with respect to the dose grid was also investigated. The maximum acceptable dose grid size depends on the gradient of dose profile and in turn the range of proton beam. In the case that only the phantom scattering was considered and that the beam was aligned with the dose grid, grid sizes from 0.4 to 6.8 mm were required for proton beams with ranges from 2 to 30 cm for 2% error limit at the Bragg peak point. A near linear relation between the maximum acceptable grid size and beam range was observed. For this analysis model, the resolution requirement was not significantly related to the orientation of the beam with respect to the grid.

Fast voxel and polygon ray-tracing algorithms in intensity modulated radiation therapy treatment planning.

We present work on combining three algorithms to improve ray-tracing efficiency in radiation therapy dose computation. The three algorithms include: An improved point-in-polygon algorithm, incremental voxel ray tracing algorithm, and stereographic projection of beamlets for voxel truncation. The point-in-polygon and incremental voxel ray-tracing algorithms have been used in computer graphics and nuclear medicine applications while the stereographic projection algorithm was developed by our group. These algorithms demonstrate significant improvements over the current standard algorithms in peer reviewed literature, i.e., the polygon and voxel ray-tracing algorithms of Siddon for voxel classification (point-in-polygon testing) and dose computation, respectively, and radius testing for voxel truncation. The presented polygon ray-tracing technique was tested on 10 intensity modulated radiation therapy (IMRT) treatment planning cases that required the classification of between 0.58 and 2.0 million voxels on a 2.5 mm isotropic dose grid into 1-4 targets and 5-14 structures represented as extruded polygons (a.k.a. Siddon prisms). Incremental voxel ray tracing and voxel truncation employing virtual stereographic projection was tested on the same IMRT treatment planning cases where voxel dose was required for 230-2400 beamlets using a finite-size pencil-beam algorithm. Between a 100 and 360 fold cpu time improvement over Siddon's method was observed for the polygon ray-tracing algorithm to perform classification of voxels for target and structure membership. Between a 2.6 and 3.1 fold reduction in cpu time over current algorithms was found for the implementation of incremental ray tracing. Additionally, voxel truncation via stereographic projection was observed to be 11-25 times faster than the radial-testing beamlet extent approach and was further improved 1.7-2.0 fold through point-classification using the method of translation over the cross product technique.

A high-speed scintillation-based electronic portal imaging device to quantitatively characterize IMRT delivery.

We have developed an electronic portal imaging device (EPID) employing a fast scintillator and a high-speed camera. The device is designed to accurately and independently characterize the fluence delivered by a linear accelerator during intensity modulated radiation therapy (IMRT) with either step-and-shoot or dynamic multileaf collimator (MLC) delivery. Our aim is to accurately obtain the beam shape and fluence of all segments delivered during IMRT, in order to study the nature of discrepancies between the plan and the delivered doses. A commercial high-speed camera was combined with a terbium-doped gadolinium-oxy-sulfide (Gd2O2S:Tb) scintillator to form an EPID for the unaliased capture of two-dimensional fluence distributions of each beam in an IMRT delivery. The high speed EPID was synchronized to the accelerator pulse-forming network and gated to capture every possible pulse emitted from the accelerator, with an approximate frame rate of 360 frames-per-second (fps). A 62-segment beam from a head-and-neck IMRT treatment plan requiring 68 s to deliver was recorded with our high speed EPID producing approximately 6 Gbytes of imaging data. The EPID data were compared with the MLC instruction files and the MLC controller log files. The frames were binned to provide a frame rate of 72 fps with a signal-to-noise ratio that was sufficient to resolve leaf positions and segment fluence. The fractional fluence from the log files and EPID data agreed well. An ambiguity in the motion of the MLC during beam on was resolved. The log files reported leaf motions at the end of 33 of the 42 segments, while the EPID observed leaf motions in only 7 of the 42 segments. The static IMRT segment shapes observed by the high speed EPID were in good agreement with the shapes reported in the log files. The leaf motions observed during beam-on for step-and-shoot delivery were not temporally resolved by the log files.

Prospects for quantitative two-dimensional radiochromic film dosimetry for low dose-rate brachytherapy sources.

Radiochromic film (RCF) has been shown to be a precise and accurate two-dimensional dosimeter for acute exposure radiation fields. However, "temporal history" mismatch between calibration and brachytherapy films due to RCF dose-rate effects could introduce potentially large uncertainties in low dose-rate (LDR) brachytherapy absolute dose measurement. This article presents a quantitative evaluation of the precision and accuracy of a laser scanner-based RCF-dosimetry system and the effect of the temporal history mismatch in LDR absolute dose measurement. MD-55-2 RCF was used to measure absolute dose for a low dose-rate 137Cs brachytherapy source using both single- and double-exposure techniques. Dose-measurement accuracy was evaluated by comparing RCF to Monte Carlo photon-transport simulation. The temporal history mismatch effect was investigated by examining dependence of RCF accuracy on irradiation-to-densitometry time interval. The predictions of the empirical cumulative dose superposition model (CDSM) were compared with measurements. For the double-exposure technique, the agreement between measurement and Monte Carlo simulation was better than 4% in the 3-60 Gy dose range with measurement precisions (coverage factor k = 1) of <2% and <6% for the doses greater or less than 3 Gy, respectively. The overall uncertainty (k = 1) of dose rate/air-kerma strength measurements achievable by this dosimetry system for a spatial resolution of 0.1 mm is less than 4% for doses greater than 5 Gy. The measured temporal history mismatch systematic error is about 1.8% for a 48 h postexposure time when using the double exposure technique and agrees with CDSM's prediction qualitatively. This work demonstrates that the model MD-55-2 RCF detector has the potential to support quantitative dose measurements about LDR brachytherapy sources with precision and accuracy better than that of previously described dosimeters. The impacts of this work on the future use of new type of RCF were also discussed.

Important considerations for radiochromic film dosimetry with flatbed CCD scanners and EBT GAFCHROMIC film.

In this study, we present three significant artifacts that have the potential to negatively impact the accuracy and precision of film dosimetry measurements made using GAFCHROMIC EBT radiochromic film when read out with CCD flatbed scanners. Films were scanned using three commonly employed instruments: a Macbeth TD932 spot densitometer, an Epson Expression 1680 CCD array scanner, and a Microtek ScanMaker i900 CCD array scanner. For the two scanners we assessed the variation in optical density (OD) of GAFCHROMIC EBT film with scanning bed position, angular rotation of the film with respect to the scan line direction, and temperature inside the scanner due to repeated scanning. Scanning uniform radiochromic films demonstrated a distinct bowing effect in profiles in the direction of the CCD array with a nonuniformity of up to 17%. Profiles along a direction orthogonal to the CCD array demonstrated a 7% variation. A strong angular dependence was found in measurements made with the flatbed scanners; the effect could not be reproduced with the spot densitometer. An IMRT quality assurance film was scanned twice rotating the film 90' between the scans. For films scanned on the Epson scanner, up to 12% variation was observed in unirradiated EBT films rotated between 0 degrees and 90 degrees, which decreased to approximately 8% for EBT films irradiated to 300 cGy. Variations of up to 80% were observed for films scanned with the Microtek scanner. The scanners were found to significantly increase the film temperature with repeated scanning. Film temperature between 18 and 33 degrees C caused OD changes of approximately 7%. Considering these effects, we recommend adherence to a strict scanning protocol that includes: maintaining the orientation of films scanned on flatbed scanners, limiting scanning to the central portion of the scanner bed, and limiting the number of consecutive scans to minimize changes in OD caused by film heating.