Verification of radiosurgery target point alignment with an electronic portal imaging device (EPID)

An Optical Reusable 2D Radiochromic Gel-Based System for Ionising Radiation Measurements in Radiotherapy

This work describes the development of a reusable 2D detector based on radiochromic reaction for radiotherapy dosimetric measurements. It consists of a radiochromic gel dosimeter in a cuboidal plastic container, scanning with a flatbed scanner, and data processing using a dedicated software package. This tool is assessed using the example of the application of the coincidence test of radiation and mechanical isocenters for a medical accelerator. The following were examined: scanning repeatability and image homogeneity, the impact of image processing on data processing in coincidence tests, and irradiation conditions-monitor units per radiation beam and irradiation field are selected. Optimal conditions for carrying out the test are chosen: (i) the multi-leaf collimator gap should preferably be 5 mm for 2D star shot irradiation, (ii) it is recommended to apply ≥2500-≤5000 MU per beam to obtain a strong signal enabling easy data processing, (iii) Mean filter can be applied to the images to improve calculations. An approach to dosimeter reuse with the goal of reducing costs is presented; the number of reuses is related to the MUs per beam, which, in this study, is about 5-57 for 30,000-2500 MU per beam (four fields). The proposed reusable system was successfully applied to the coincidence tests, confirming its suitability as a new potential quality assurance tool in radiotherapy.

Fast Isocenter Determination Using 3D Polymer Gel Dosimetry with Kilovoltage Cone-Beam CT Reading and the PolyGeVero-CT Software Package for Linac Quality Assurance in Radiotherapy

This work presents an approach to the fast determination of a medical accelerator irradiation isocenter as a quality assurance (QA) procedure in radiotherapy. The isocenter determination tool is the tissue equivalent high-resolution 3D polymer gel dosimeter (PABIG^nx) in a dedicated container combined with kilovoltage imaging systems and the polyGeVero-CT software package (v. 1.2, GeVero Co., Poland). Two accelerators were employed: Halcyon and TrueBeam (Varian, USA), both equipped with cone beam computed tomography (CBCT) and iterative reconstruction CBCT (iCBCT) algorithms. The scope of this work includes: (i) the examination of factors influencing image quality (reconstruction algorithms and modes), radiation field parameters (dose and multi-leaf collimator (MLC) gaps), fiducial markers, signal averaging for reconstruction algorithms and the scanning time interval between consecutive scans, (ii) the examination of factors influencing the isocenter determination, image processing (signal averaging, background subtraction, image filtering) and (iii) an isocenter determination report using a 2D and 3D approach. An optimized protocol and isocenter determination conditions were found. The time and effort required to determine the isocenter are discussed.

Positions of radiation isocenter and the couch rotation center established by Winston-Lutz and optical measurements

Purpose

To compare x-ray and optical imaging methods for measuring the relative position of radiation isocenter and couch rotation center. To show the impact of radiation isocenter size and target motion on the margins for target contours.

Methods

Winston-Lutz measurements are made using EPID images. Image analysis was done with public domain software, ImageJ, and spreadsheets written in Microsoft Excel. A comparison between the center of a high density test object and center of the MLC collimated beam is used to judge the relative position of the radiation isocenter in space for gantry and couch rotation. Additionally, motion of the target with couch rotation is determined with an optical imaging system. Five different accelerators, two TrueBeams, a Trilogy, and two VersaHDs, were assessed by Winston-Lutz and optical methods.

Results

The shift in the radiation isocenter with gantry rotation is found to be a tri-axial ellipsoid. Shifts in the target position with respect to radiation isocenter with couch rotation were between 0.4 and 0.6 mm. The Winston-Lutz and optical method determination of couch rotation center agreed within measurement uncertainty.

Conclusions

Image analysis yields precise data on linear accelerator radiation isocenter and rotation centers of the couch. The Winston-Lutz and optical methods agreed within measurement uncertainty.

A two-layer cylinder phantom developed for film-based isocenter verification of radiotherapy machine

PURPOSE

A two-layer cylinder (TLC) phantom was developed for simplifying film-based isocenter verification of linear accelerators in radiotherapy.

METHODS AND MATERIALS

The phantom mainly consists of two parts: (1) two nested solid cylinders between which a radiochromic film can be inserted and irradiated; (2) a tungsten ball supported by a thin rod and located at the phantom center for alignment with the mechanical isocenter. In practice, the phantom was first positioned by the room laser to align the tungsten ball to the mechanical isocenter of the linear accelerator. Then, a radiochromic film was precisely inserted into the gap between the two cylinders of the phantom and irradiated by beams with preset gantry and couch angles. Later the irradiated film was scanned and processed by an in-house developed analysis software. Finally, the offset of the radiation isocenter from the mechanical isocenter was determined by the built-in three-dimensional (3D) reconstruction algorithms. The accuracy of this method was evaluated by positioning the phantom with a known couch shift, then checking the residual error after couch shift correction. The reliability of this method was evaluated by comparing the calculated offset with the corresponding result determined by the traditional star-shot method.

RESULTS

For the accuracy test, the residual errors were -0.14 ± 0.03 mm, 0.05 ± 0.06 mm, and 0.05 ± 0.06 mm in the lateral, longitudinal, and vertical axes, respectively. For the reliability test, the differences between the calculated offset and the result determined by the star-shot method were -0.10 mm, 0.12 mm, and 0.12 mm in the lateral, longitudinal, and vertical axes, respectively.

CONCLUSION

The proposed method is able to reconstruct beams in 3D with one film, which is more time-saving and accurate. Additionally, with this design, the phantom positioning, film loading, beam delivery, and data analyzing are simpler. This phantom and analysis software provides an efficient and effective way to perform film-based isocenter verification of linear accelerators in radiotherapy.

A novel method to determine linac mechanical isocenter position and size and examples of specific QA applications

The most important geometric characteristic of stereotactic treatment is the accuracy of positioning the target at the treatment isocenter and the accuracy of directing the radiation beam at the treatment isocenter. Commonly, the radiation isocenter is used as the reference for the treatment isocenter, but its method of localization is not strictly defined, and it depends on the linac-specific beam steering parameters. A novel method is presented for determining the linac mechanical isocenter position and size based on the localization of the collimator axis of rotation at arbitrary gantry angle. The collimator axis of rotation position is determined from the radiation beam center position corrected for the focal spot offset. The focal spot offset is determined using the image center shift method with a custom-design rigid phantom with two sets of ball-bearings. Three specific quality assurance (QA) applications and assessment methods are also presented to demonstrate the functionality of linac mechanical isocenter position and size determination in clinical practice. The first is a mechanical and radiation isocenters coincidence test suitable for quick congruence assessment of these two isocenters for a selected energy, usually required after a nonroutine linac repair and/or energy adjustment. The second is a stereotactic beam isocentricity assessment suitable for pretreatment stereotactic QA. The third is a comprehensive linac geometrical performance test suitable for routine linac QA. The uncertainties of the method for determining mechanical isocenter position and size were measured to be 0.05 mm and 0.04 mm, respectively, using four available photon energies, and were significantly smaller than those of determining the radiation isocenter position and size, which were 0.36 mm and 0.12 mm respectively. It is therefore recommended that the mechanical isocenter position and size be used as the reference linac treatment isocenter and a linac mechanical characteristic parameter respectively.

Portal dosimetry scripting application programming interface (PDSAPI) for Winston-Lutz test employing ceramic balls

PURPOSE

Stereotactic radiosurgery (SRS) and stereotactic body radiation therapy (SBRT) treatments require a high degree of accuracy. Mechanical, imaging, and radiation isocenter coincidence is especially important. As a common method, the Winston-Lutz (WL) test plays an important role. However, weekly or daily WL test can be very time consuming. We developed novel methods using Portal Dosimetry Scripting Application Programming Interface (PDSAPI) to facilitate the test as well as documentation.

METHODS

Winston-Lutz PDSAPI was developed and tested on our routine weekly WL imaging. The results were compared against two commercially available software RIT (Radiological Imaging Technology, Colorado Springs, CO) and DoseLab (Varian Medical Systems, Inc. Palo Alto, CA). Two manual methods that served as ground truth were used to verify PDSAPI results. Twenty WL test image data sets (10 fields per tests, and 200 images in total) were analyzed by these five methods in this report.

RESULTS

More than 99.5% of WL PDSAPI 1D shifts agreed with each of four other methods within ±0.33 mm, which is roughly the pixel width of a-Si 1200 portal imager when source to imager distance (SID) is at 100 cm. 1D shifts agreement for ±0.22 mm and 0.11 mm were 96% and 63%, respectively. Same trend was observed for 2D displacement.

CONCLUSIONS

Winston-Lutz PDSAPI delivers similar accuracy as two commercial applications for WL test. This new application can save time spent transferring data and has the potential to implement daily WL test with reasonable test time. It also provides the data storage capability, and enables easy access to imaging and shift data.

3D isocentricity analysis for clinical linear accelerators

PURPOSE

To perform a three-dimensional (3D) concurrent isocentricity measurement of a clinical linear accelerator's (linac) using a single 3D dosimeter, PRESAGE.

METHODS

A 3D dosimeter, PRESAGE, set up on the treatment couch of a Varian TrueBeam LINAC using the setup lasers, was irradiated under gantry angles of 0 ∘ , 50 ∘ , 160 ∘ , and 270 ∘ with the couch fixed at 0 ∘ and subsequently, under couch angles of 10 ∘ , 330 ∘ , 300 ∘ , and 265 ∘ with the gantry fixed at 270 ∘ . The 1 cm 2 (at 100 cm SAD) square fields were delivered at 6 MV with 800 MU/field. After irradiation, the dosimeter was scanned using a single-beam optical scanner and images were reconstructed with submillimeter resolution using filtered back-projection. Postprocessing was used to extract views parallel to the star-shot planes from which beam trajectories and the smallest circles enclosing these were drawn and extracted. These circles and information from the view orthogonal to both star-shots were used to represent the rotational centers as spheroids. The linac isocenter was defined by the distribution of midpoints between any, randomly selected, points lying inside the center spheroids defined by the table and gantry rotations; isocenter location and size were defined by the average midpoint and the distribution's semi-axes. Collimator rotations were not included in this study.

RESULTS

Relative to the setup center defined by lasers, the table and gantry rotation center coordinates (lat., long., vert.) were measured in units of millimeters, to be (-0.24, 0.18, -0.52) and (0.10, 0.53, -0.52), respectively. Displacements from the setup center were 0.60 and 0.75 mm for the table and gantry centers, while the distance between them measured 0.49 mm. The linac's radiation isocenter was calculated to be at (-0.07, -0.17, 0.51) relative to the setup lasers and its size was found to be most easily described by a spheroid prolate in vertical direction with semi-axis lengths of 0.13 and 0.23 mm for the lateral-longitudinal and vertical directions, respectively.

CONCLUSIONS

This study demonstrates how to measure the location and sizes of rotational centers in 3D with one setup. The proposed method provides a more comprehensive view on the isocentricity of LINAC than the conventional two-dimensional film measurements. Additionally, a new definition of isocenter and its size was proposed.

Beam focal spot intrafraction motion and gantry angle dependence: A study of Varian linac focal spot alignment

The characteristics of the focal spot of the linear accelerator (linac) play a role in determining the resulting dose distribution within the patient, and hence probability of treatment success. A direct measurement of focal spot position is not recommended by AAPM Task Group 142, but factors influenced by focal spot position, such as beam symmetry and isocentre position, are. Traditional methods of measuring focal spot position are time consuming and can only be performed at gantry 0°. The presented method has been proposed using a phantom of novel design to accurately measure the position of the focal spot relative to the collimator's axis of rotation (CAX) at any gantry angle, and to measure the intra-fraction movement of the focal spot relative to the mean position during treatment. The method was reproducible to within 0.012 mm/0.029 mm (mean/max) for the three Varian linacs tested. The focal spot position was shown to deviate from the CAX by up to 0.386 mm during gantry rotation. The focal spot position was more unstable at the start of treatment, with the worst performing linac having an initial displacement of up to 0.15 mm from its mean position before stabilizing to within 0.01 mm after 3 s. The method proposed is a beneficial addition to the quality assurance (QA) schedule of any clinic, allowing quick determination of source position and movement at any gantry angle. Measurement of focal spot allows the possibility of fine-tuning the electron beam steering system to improve the standard of the photon beam and of stereotactic treatments.

A study of Winston-Lutz test on two different electronic portal imaging devices and with low energy imaging

Stereotactic radiosurgery requires sub-millimetre accuracy in patient positioning and target localization. Therefore, verification of the linear accelerator (linac) isocentre and the laser alignment to the isocentre is performed in some clinics prior to the treatment using the Winston-Lutz (W-L) test with films and more recently with images obtained using the electronic portal imaging devices (EPID). The W-L test is performed by acquiring EPID images of a radio-opaque ball of 6 mm diameter (the W-L phantom) placed at the isocentre of the linac at various gantry and table angles, with a predefined small square or circular radiation beam. In this study, the W-L test was performed on two linacs having EPIDs of different size and resolution, viz, a TrueBeam™ linac with aS1000 EPID of size 40 × 30 cm(2) with 1024 × 768 pixel resolution and an EDGE™ linac having an EPID of size 43 × 43 cm(2) with pixel resolution of 1280 × 1280. In order to determine the displacement of the radio-opaque ball centre from the radiation beam centre of the W-L test, an in-house MATLAB™ image processing code was developed using morphological operations. The displacement in radiation beam centre at each gantry and couch position was obtained by determining the distance between the radiation field centre and the radio-opaque ball centre for every image. Since the MATLAB code was based on image processing that was dependent on the image contrast and resolution, the W-L test was also compared for images obtained with different beam energies. The W-L tests were performed for 6 and 8 MV beams on the TrueBeam™ linac and for 2.5 and 6 MV beams on the EDGE™ linac with a higher resolution EPID. It was observed that the images obtained with the EPID of higher resolution resulted in same accuracy in the determination of the displacement between the centres of the radio-opaque ball and the radiation beam, and significant difference was not observed with images acquired with different energies. It is concluded that the software based on morphological operations provided an accurate estimation of the displacement of the ball centre from the radiation beam center.

On the selection of gantry and collimator angles for isocenter localization using Winston-Lutz tests

In Winston-Lutz (WL) tests, the isocenter of a linear accelerator (linac) is determined as the intersection of radiation central axes (CAX) from multiple gantry, collimator, and couch angles. It is well known that the CAX can wobble due to mechanical imperfections of the linac. Previous studies suggested that the wobble varies with gantry and collimator angles. Therefore, the isocenter determined in the WL tests has a profound dependence on the gantry and collimator angles at which CAX are sampled. In this study, we evaluated the systematic and random errors in the iso-centers determined with different CAX sampling schemes. Digital WL tests were performed on six linacs. For each WL test, 63 CAX were sampled at nine gantry angles and seven collimator angles. Subsets of these data were used to simulate the effects of various CAX sampling schemes. An isocenter was calculated from each subset of CAX and compared against the reference isocenter, which was calculated from 48 opposing CAX. The differences between the calculated isocenters and the reference isocenters ranged from 0 to 0.8 mm. The differences diminished to less than 0.2 mm when 24 or more CAX were sampled. Isocenters determined with collimator 0° were vertically lower than those determined with collimator 90° and 270°. Isocenter localization errors in the longitudinal direction (along the axis of gantry rotation) showed a strong dependence on the collimator angle selected. The errors in all directions were significantly reduced when opposing collimator angles and opposing gantry angles were employed. The isocenter localization errors were less than 0.2 mm with the common CAX sampling scheme, which used four cardinal gantry angles and two opposing collimator angles. Reproducibility stud-ies on one linac showed that the mean and maximum variations of CAX during the WL tests were 0.053 mm and 0.30 mm, respectively. The maximal variation in the resulting isocenters was 0.068 mm if 48 CAX were used, or 0.13 mm if four CAX were used. Quantitative results from this study are useful for understanding and minimizing the isocenter uncertainty in WL tests.

Isocenter verification for linac-based stereotactic radiation therapy: review of principles and techniques

There have been several manual, semi-automatic and fully-automatic methods proposed for verification of the position of mechanical isocenter as part of comprehensive quality assurance programs required for linear accelerator-based stereotactic radiosurgery/radiotherapy (SRS/SRT) treatments. In this paper, a systematic review has been carried out to discuss the present methods for isocenter verification and compare their characteristics, to help physicists in making a decision on selection of their quality assurance routine.

Initial application of a geometric QA tool for integrated MV and kV imaging systems on three image guided radiotherapy systems

PURPOSE

Several linacs with integrated kilovoltage (kV) imaging have been developed for delivery of image guided radiation therapy (IGRT). High geometric accuracy and coincidence of kV imaging systems and megavoltage (MV) beam delivery are essential for successful image guidance. A geometric QA tool has been adapted for routine QA for evaluating and characterizing the geometric accuracy of kV and MV cone-beam imaging systems. The purpose of this work is to demonstrate the application of methodology to routine QA across three IGRT-dedicated linac platforms.

METHODS

It has been applied to a Varian Trilogy (Varian Medical Systems, Palo Alto, CA), an Elekta SynergyS (Elekta, Stockholm, Sweden), and a Brainlab Vero (Brainlab AG, Feldkirchen, Germany). Both the Trilogy and SynergyS linacs are equipped with a retractable kV x-ray tube and a flat panel detector. The Vero utilizes a rotating, rigid ring structure integrating a MV x-ray head mounted on orthogonal gimbals, an electronic portal imaging device (EPID), two kV x-ray tubes, and two fixed flat panel detectors. This dual kV imaging system provides orthogonal radiographs, CBCT images, and real-time fluoroscopic monitoring. Two QA phantoms were built to suit different field sizes. Projection images of a QA phantom were acquired using MV and kV imaging systems at a series of gantry angles. Software developed for this study was used to analyze the projection images and calculate nine geometric parameters for each projection. The Trilogy was characterized five times over one year, while the SynergyS was characterized four times and the Vero once. Over 6500 individual projections were acquired and analyzed. Quantitative geometric parameters of both MV and kV imaging systems, as well as the isocenter consistency of the imaging systems, were successfully evaluated.

RESULTS

A geometric tool has been successfully implemented for calibration and QA of integrated kV and MV across a variety of radiotherapy platforms. X-ray source angle deviations up to 0.8 degrees, and detector center offsets up to 3 mm, were observed for three linacs, with the exception of the Vero, for which a significant center offset of one kV detector (prior to machine commissioning) was observed. In contrast, the gimbal-based MV source positioning of the Vero demonstrated differences between observed and expected source positions of less than 0.2 mm, both with and without gimbal rotation.

CONCLUSIONS

This initial application of this geometric QA tool shows promise as a universal, independent tool for quantitative evaluation of geometric accuracies of both MV and integrated kV imaging systems across a range of platforms. It provides nine geometric parameters of any imaging system at every gantry angle as well as the isocenter coincidence of the MV and kV image systems.

A fast double template convolution isocenter evaluation algorithm with subpixel accuracy

PURPOSE

To design a fast Winston Lutz (fWL) algorithm for accurate analysis of radiation isocenter from images without edge detection or center of mass calculations.

METHODS

An algorithm has been developed to implement the Winston Lutz test for mechanical/ radiation isocenter agreement using an electronic portal imaging device (EPID). The algorithm detects the position of the radiation shadow of a tungsten ball within a stereotactic cone. The fWL algorithm employs a double convolution to independently find the position of the sphere and cone centers. Subpixel estimation is used to achieve high accuracy. Results of the algorithm were compared to (1) a human observer with template guidance and (2) an edge detection/center of mass (edCOM) algorithm. Testing was performed with high resolution (0.05 mm/px, film) and low resolution (0.78 mm/px, EPID) image sets.

RESULTS

Sphere and cone center relative positions were calculated with the fWL algorithm for high resolution test images with an accuracy of 0.002 +/- 0.061 mm compared to 0.042 +/- 0.294 mm for the human observer, and 0.003 +/- 0.038 mm for the edCOM algorithm. The fWL algorithm required 0.01 s per image compared to 5 s for the edCOM algorithm and 20 s for the human observer. For lower resolution images the fWL algorithm localized the centers with an accuracy of 0.083 +/- 0.12 mm compared to 0.03 +/- 0.5514 mm for the edCOM algorithm.

CONCLUSIONS

A fast (subsecond) subpixel algorithm has been developed that can accurately determine the center locations of the ball and cone in Winston Lutz test images without edge detection or COM calculations.

An automatic method of the isocentre position verification for micromultileaf collimator based radiosurgery system

An efficient procedure has been developed using an electronic portal imaging device (EPID) and in-house written software for verification of a target simulator alignment with the radiation isocentre. A 5 mm tungsten ball is aligned to a linac isocentre based on a lasers intersection point. The BrainLab(®) m3™ add-on multileaf collimator (MLC) forms a rectangular open field of 1.8 × 1.8 cm(2). At five different gantry and couch positions, EPID images are acquired. A computer search algorithm determines the centres of both a radiation field and a tungsten ball for each image. Based on the geometric differences between those centres, the optimum three-dimensional shift of a tungsten ball is calculated in order to minimise the misalignment error between a target simulator and a radiation isocentre. A decision can then be made whether or not the tungsten ball and lasers intersection point should be corrected. The accuracy and precision of the procedure has been tested and found to be 0.04 and 0.24 mm respectively at 95% confidence interval. The procedure is also quicker, easier and more reliable to perform compared to the previous method based on irradiating a radiographic film.

A quality assurance method with submillimeter accuracy for stereotactic linear accelerators

The Stereotactic Alignment for Linear Accelerator (S. A. Linac) system is developed to conveniently improve the alignment accuracy of a conventional linac equipped with stereotactic cones. From the Winston‐Lutz test, the SAlinac system performs three‐dimensional (3D) reconstruction of the quality assurance (QA) ball coordinates with respect to the radiation isocenter, and combines this information with digital images of the laser target to determine the absolute position of the room lasers. A handheld device provides near‐real‐time repositioning advice to enable the user to align the QA ball and room lasers to within 0.25 mm of the centroid of the radiation isocenter. The results of 37 Winston‐Lutz tests over 68 days showed that the median 3D QA ball alignment error was 0.09 mm, and 97% of the time the 3D error was ≤0.25 mm. All 3D isocentric errors in the study were 0.3 mm or less. The median x and y laser alignment coordinate error was 0.09 mm, and 94% of the time the x and y laser error was ≤0.25 mm. A phantom test showed that the system can make submillimeter end‐to‐end accuracy achievable, making a conventional linac a “Submillimeter Knife”.

PACS numbers: 87.53.Ly, 87.55.Qr

A simple method to quantify the coincidence between portal image graticules and radiation field centers or radiation isocenter

PURPOSE

The aim of this study was to develop a computerized method to quantify the coincidence between portal image graticules and radiation field centers or radiation isocenter. Three types of graticules were included in this study: Megavoltage (MV) mechanical graticule, MV electronic portal imaging device digital graticule, and kilovoltage (kV) on-board imaging digital graticule.

METHODS

A metal ball bearing (BB) was imaged with MV and kV x-ray beams in a procedure similar to a Winston-Lutz test. The radiation fields, graticules, and BB were localized in eight portal images using Hough transform-based computer algorithms. The center of the BB served as a static reference point in the 3D space so that the distances between the graticule centers and the radiation field centers were calculated. The radiation isocenter was determined from the radiation field centers at different gantry angles.

RESULTS

Misalignments of MV and kV portal imaging graticules varied with the gantry or x-ray source angle as a result of mechanical imperfections of the linear accelerator and its imaging system. While the three graticules in this study were aligned to the radiation field centers and the radiation isocenter within 2.0 mm, misalignments of 1.5-2.0 mm were found at certain gantry angles. These misalignments were highly reproducible with the gantry rotation.

CONCLUSIONS

A simple method was developed to quantify the alignments of portal image graticules directly against the radiation field centers or the radiation isocenter. The advantage of this method is that it does not require the BB to be placed exactly at the radiation isocenter through a precalibrated surrogating device such as room lasers or light field crosshairs. The present method is useful for radiation therapy modalities that require high-precision portal imaging such as image-guided stereotactic radiotherapy.

A fiducial detection algorithm for real-time image guided IMRT based on simultaneous MV and kV imaging

The advantage of highly conformal dose techniques such as 3DCRT and IMRT is limited by intrafraction organ motion. A new approach to gain near real-time 3D positions of internally implanted fiducial markers is to analyze simultaneous onboard kV beam and treatment MV beam images (from fluoroscopic or electronic portal image devices). Before we can use this real-time image guidance for clinical 3DCRT and IMRT treatments, four outstanding issues need to be addressed. (1) How will fiducial motion blur the image and hinder tracking fiducials? kV and MV images are acquired while the tumor is moving at various speeds. We find that a fiducial can be successfully detected at a maximum linear speed of 1.6 cm/s. (2) How does MV beam scattering affect kV imaging? We investigate this by varying MV field size and kV source to imager distance, and find that common treatment MV beams do not hinder fiducial detection in simultaneous kV images. (3) How can one detect fiducials on images from 3DCRT and IMRT treatment beams when the MV fields are modified by a multileaf collimator (MLC)? The presented analysis is capable of segmenting a MV field from the blocking MLC and detecting visible fiducials. This enables the calculation of nearly real-time 3D positions of markers during a real treatment. (4) Is the analysis fast enough to track fiducials in nearly real time? Multiple methods are adopted to predict marker positions and reduce search regions. The average detection time per frame for three markers in a 1024 x 768 image was reduced to 0.1 s or less. Solving these four issues paves the way to tracking moving fiducial markers throughout a 3DCRT or IMRT treatment. Altogether, these four studies demonstrate that our algorithm can track fiducials in real time, on degraded kV images (MV scatter), in rapidly moving tumors (fiducial blurring), and even provide useful information in the case when some fiducials are blocked from view by the MLC. This technique can provide a gating signal or be used for intra-fractional tumor tracking on a Linac equipped with a kV imaging system. Any motion exceeding a preset threshold can warn the therapist to suspend a treatment session and reposition the patient.

Development of a QA phantom and automated analysis tool for geometric quality assurance of on-board MV and kV x-ray imaging systems

The medical linear accelerator (linac) integrated with a kilovoltage (kV) flat-panel imager has been emerging as an important piece of equipment for image-guided radiation therapy. Due to the sagging of the linac head and the flexing of the robotic arms that mount the x-ray tube and flat-panel detector, geometric nonidealities generally exist in the imaging geometry no matter whether it is for the two-dimensional projection image or three-dimensional cone-beam computed tomography. Normally, the geometric parameters are established during the commissioning and incorporated in correction software in respective image formation or reconstruction. A prudent use of an on-board imaging system necessitates a routine surveillance of the geometric accuracy of the system like the position of the x-ray source, imager position and orientation, isocenter, rotation trajectory, and source-to-imager distance. Here we describe a purposely built phantom and a data analysis software for monitoring these important parameters of the system in an efficient and automated way. The developed tool works equally well for the megavoltage (MV) electronic portal imaging device and hence allows us to measure the coincidence of the isocenters of the MV and kV beams of the linac. This QA tool can detect an angular uncertainty of 0.1 degrees of the x-ray source. For spatial uncertainties, such as the source position, the imager position, or the kV/MV isocenter misalignment, the demonstrated accuracy of this tool was better than 1.6 mm. The developed tool provides us with a simple, robust, and objective way to probe and monitor the geometric status of an imaging system in a fully automatic process and facilitate routine QA workflow in a clinic.

Fast internal marker tracking algorithm for onboard MV and kV imaging systems

Intrafraction organ motion can limit the advantage of highly conformal dose techniques such as intensity modulated radiation therapy (IMRT) due to target position uncertainty. To ensure high accuracy in beam targeting, real-time knowledge of the target location is highly desired throughout the beam delivery process. This knowledge can be gained through imaging of internally implanted radio-opaque markers with fluoroscopic or electronic portal imaging devices (EPID). In the case of MV based images, marker detection can be problematic due to the significantly lower contrast between different materials in comparison to their kV-based counterparts. This work presents a fully automated algorithm capable of detecting implanted metallic markers in both kV and MV images with high consistency. Using prior CT information, the algorithm predefines the volumetric search space without manual region-of-interest (ROI) selection by the user. Depending on the template selected, both spherical and cylindrical markers can be detected. Multiple markers can be simultaneously tracked without indexing confusion. Phantom studies show detection success rates of 100% for both kV and MV image data. In addition, application of the algorithm to real patient image data results in successful detection of all implanted markers for MV images. Near real-time operational speeds of approximately 10 frames/sec for the detection of five markers in a 1024 x 768 image are accomplished using an ordinary PC workstation.

An MLC-based linac QA procedure for the characterization of radiation isocenter and room lasers’ position

We have designed and implemented a new stereotactic linac QA test with stereotactic precision. The test is used to characterize gantry sag, couch wobble, cone placement, MLC offsets, and room lasers' positions relative to the radiation isocenter. Two MLC star patterns, a cone pattern, and the laser line patterns are recorded on the same imaging medium. Phosphor plates are used as imaging medium due to their sensitivity to red light. The red light of room lasers erases some of the irradiation information stored on the phosphor plates enabling accurate and direct measurements for the position of room lasers and radiation isocenter. Using film instead of the phosphor plate as imaging medium is possible, however, it is less practical. The QA method consists of irradiating four phosphor plates that record the gantry sag between the 0 degrees and 180 degrees gantry angles, the position and stability of couch rotational axis, the sag between the 90 degrees and 270 degrees gantry angles, the accuracy of cone placement on the collimator, the MLC offsets from the collimator rotational axis, and the position of laser lines relative to the radiation isocenter. The estimated accuracy of the method is +/- 0.2 mm. The observed reproducibility of the method is about +/- 0.1 mm. The total irradiation/ illumination time is about 10 min per image. Data analysis, including the phosphor plate scanning, takes less than 5 min for each image. The method characterizes the radiation isocenter geometry with the high accuracy required for the stereotactic radiosurgery. In this respect, it is similar to the standard ball test for stereotactic machines. However, due to the usage of the MLC instead of the cross-hair/ball, it does not depend on the cross-hair/ball placement errors with respect to the lasers and it provides more information on the mechanical integrity of the linac/couch/laser system. Alternatively, it can be used as a highly accurate QA procedure for the nonstereotactic machines. Noteworthy is its ability to characterize the MLC position accuracy, which is an important factor in IMRT delivery.

Theoretical foundation for real-time prostate localization using an inductively coupled transmitter and a superconducting quantum interference device (SQUID) magnetometer system

Real‐time, 3D localization of the prostate for intensity‐modulated radiotherapy can be accomplished with passively charged radio frequency transmitters and superconducting quantum interference device (SQUID) magnetometers. The overall system design consists of an external dipole antenna as a power source for charging a microchip implant transmitter and SQUID magnetometers for signal detection. An external dipole antenna charges an on‐chip capacitor through inductive coupling in the near field region through a small implant inductor. The charge and discharge sequence between the external antenna and the implant circuit can be defined by half duplex, full duplex, or sequential operations. The resulting implant discharge current creates an alternating magnetic field through the inductor. The field is detected by the surrounding magnetometers, and the location of the implant transmitter can be calculated. Problems associated with this system design are interrelated with the signal strength at the detectors, detector sensitivity, and charge time of the implant capacitor. The physical parameters required for optimizing the system for real‐time applications are the operating frequency, implant inductance and capacitance, the external dipole current and loop radius, the detector distance, and mutual inductance. Consequently, the sequential operating mode is the best choice for real‐time localization for constraints requiring positioning within 1 s due to the mutual inductance and detector sensitivity. We present the theoretical foundation for designing a real‐time, 3D prostate localization system including the associated physical parameters and demonstrate the feasibility and physical limitations for such a system.

PACS number: 87.53.‐j

Introducing a system for automated control of rotation axes, collimator and laser adjustment for a medical linear accelerator

Mechanical stability and precise adjustment of rotation axes, collimator and room lasers are essential for the success of radiotherapy and particularly stereotactic radiosurgery with a linear accelerator. Quality assurance procedures, at present mainly based on visual tests and radiographic film evaluations, should desirably be little time consuming and highly accurate. We present a method based on segmentation and analysis of digital images acquired with an electronic portal imaging device (EPID) that meets these objectives. The method can be employed for routine quality assurance with a square field formed by the built-in collimator jaws as well as with a circular field using an external drill hole collimator. A number of tests, performed to evaluate accuracy and reproducibility of the algorithm, yielded very satisfying results. Studies performed over a period of 18 months prove the applicability of the inspected accelerator for stereotactic radiosurgery.

An isocenter position verification device for electronic portal imaging: physical and dosimetrical characteristics

The physical and dosimetrical characteristics of a device, designed to visualize the isocenter position on electronic portal images, were examined. The device, to be mounted on the gantry head of the accelerator, containing five spheric lead markers, was designed in order to visualize the isocenter position on portal images. A quality control device was designed to check the reliability of this technique. The disturbance of the dose distribution by the markers was studied with gel dosimetry. The use of markers resulted in a precise and accurate method to visualize the isocenter on portal images. A maximum underdosage of 11%, due to attenuation by the markers, was observed. The use of markers to visualize the isocenter position on portal images, is a fast and reliable method when analyzing patient setup errors with online electronic portal imaging.

Clinical use of electronic portal imaging: report of AAPM Radiation Therapy Committee Task Group 58

AAPM Task Group 58 was created to provide materials to help the medical physicist and colleagues succeed in the clinical implementation of electronic portal imaging devices (EPIDs) in radiation oncology. This complex technology has matured over the past decade and is capable of being integrated into routine practice. However, the difficulties encountered during the specification, installation, and implementation process can be overwhelming. TG58 was charged with providing sufficient information to allow the users to overcome these difficulties and put EPIDs into routine clinical practice. In answering the charge, this report provides; comprehensive information about the physics and technology of currently available EPID systems; a detailed discussion of the steps required for successful clinical implementation, based on accumulated experience; a review of software tools available and clinical use protocols to enhance EPID utilization; and specific quality assurance requirements for initial and continuing clinical use of the systems. Specific recommendations are summarized to assist the reader with successful implementation and continuing use of an EPID.

Guide to clinical use of electronic portal imaging

The Electronic Portal Imaging Device (EPID) provides localization quality images and computer-aided analysis, which should in principal, replace portal film imaging. Modern EPIDs deliver superior image quality and an array of analysis tools that improve clinical decision making. It has been demonstrated that the EPID can be a powerful tool in the reduction of treatment setup errors and the quality assurance and verification of complex treatments. However, in many radiation therapy clinics EPID technology is not in routine clinical use. This low utilization suggests that the capability and potential of the technology alone do not guarantee its full adoption. This paper addresses basic considerations required to facilitate clinical implementation of the EPID technology and gives specific examples of successful implementations.

Setup verification in linac-based radiosurgery

A semi-automatic technique for the direct setup alignment of radiosurgical circular fields from an isocentric linac to treatment room laser cross-hairs is described. Alignment is achieved by acquiring images of the treatment room positioning laser cross-hairs superimposed on the radiosurgical circular field image. An alignment algorithm calculates the center of the radiosurgical field image as well as the intersection of the laser cross-hairs. This determines any alignment deviations and the information is then used to translate the radiosurgical collimator to its correct aligned position. Two detectors, each being sensitive to the lasers and ionizing radiation, were used to acquire the radiation/laser images. The first detector consists of a 0.3-mm-thick layer of photoconducting a-Se deposited on a 1.5-mm-thick copper plate and the second is film. The algorithm and detector system can detect deviations with a precision of approximately 0.04 mm. A device with gyroscopic degrees of freedom was built in order to firmly hold the detector at any orientation perpendicular to the radiosurgical beam axis. This device was used in conjunction with our alignment algorithm to quantify the isocentric sphere relative to the treatment room lasers over all gantry and couch angles used in dynamic stereotactic radiosurgery.

Linear accelerator output variations and their consequences for megavoltage imaging

An experimental study of radiation output intensity fluctuations of a Philips SL25 linear accelerator is presented. Measurements are obtained using an electronic portal imaging device, and the consequences of the measured fluctuations for various different applications of megavoltage imaging including portal imaging, transit dosimetry and megavoltage computed tomography (MVCT) are discussed with examples. Fluctuations in output of +/- 0.7% (1 SD) are seen on every radiation pulse after photon noise and uncertainties caused by the detection system have been accounted for. Large fluctuations are also seen during the initial beam stabilization period (15%), during normal accelerator operation after the beam has been on for more than 1 min (4.5%) and during are therapy as a repeatable function of gantry angle (9%). Such output intensity fluctuations are shown to produce image artifacts in portal imaging devices with scanned detector readout and can also produce systematic errors in detector calibration that would lead to uncertainty in transit dose calculations. The propagation of these intensity fluctuations through MVCT image reconstruction is shown to produce ring artifacts in the reconstructed image. Sample portal and MVCT images are presented. All observed fluctuations in accelerator output are well within the manufacturer's specifications and do not affect the total dose delivered during normal treatment. Finally, megavoltage imaging is shown to be a powerful tool for accelerator quality assurance and treatment verification.