A robust Hough transform algorithm for determining the radiation centers of circular and rectangular fields with subpixel accuracy

A customizable, open-source Winston-Lutz system for multi-target, single isocentre radiotherapy

OBJECTIVE

To present and share an open-source system (phantom and software) for verifying the targeting accuracy of linac-based, single-isocenter, multi-target radiotherapy. This quality assurance test extends the traditional Winston-Lutz test, which considers a single target located at isocentre.

APPROACH

Plans for a 3D-printed phantom are provided, which can be customized to accommodate various target (BB) positions. Given BB positions and gantry/collimator/couch combinations, the software generates multi-leaf collimator positions to facilitate multi-target Winston-Lutz (MTWL) plan creation. The software determines deviations between detected and expected BB positions on MV images resulting from MTWL plan delivery. BBs are located using a Hough circle detection algorithm, which is modified to favour the detection of circles: (1) having reasonable size, (2) that are contained within the radiation field, and (3) having reasonable pixel intensities. Validation was performed in two ways: (1) using synthetic data with zero targeting errors and (2) by measuring real linac targeting errors and comparing against results obtained using a commercial system.

MAIN RESULTS

Validation using the synthetic data yielded a mean (maximum) absolute discrepancy of 0.11 mm (0.21 mm), which is comparable to the synthetic phantom resolution (0.2 mm). The mean (maximum) absolute discrepancy compared to the commercial system is 0.13 mm (0.43 mm). These values are similar to results obtained with repeated deliveries of the same MTWL plan with the same phantom setup. Both validation tests yield reasonable results and are therefore considered successful. The MTWL test was performed independently by three physicists on two linacs to investigate repeatability, resulting in a mean (maximum) absolute discrepancy of 0.14 mm (0.51 mm) among the various attempts.

SIGNIFICANCE

Successful completion of this quality assurance test, using our customizable and open-source system, provides confidence that multi-target, single isocentre radiotherapy treatments can be delivered with sufficient geometric accuracy according to the chosen tolerance level.

Design and automation of specific geometric quality controls for cone-radiosurgery treatments

PURPOSE

The Varian TrueBeam STx linac can be equipped with BrainLAB stereotactic cones and ExacTrac imaging system for SRS treatments. However, these two third-party systems lack integration in a self-performance diagnosis tool dedicated to the SRS platform. The aim of this work was to design and automate essential geometric tests considering the complete set of cones with diameter range from 4 to 15 mm.

METHODS

EPID-based tests were focused on the cone alignment, the radiation isocentricity and the isocenters congruence. Images acquired with or without the BrainLAB pointer were analysed using the Hough transform and morphological filtering operations, after assessment of the algorithm accuracy using simulated images. The new Machine QA program was experimented over one year.

RESULTS

A subpixel resolution of 0.02 mm was found for the circular-field center detection algorithm. The tests results did not depend on the pointer location. The maximum deviations reported were in accordance with the AAPM SRS guidelines. The accurate localization of the linac radiation isocenter allowed for guidance of the ExacTrac calibration. A misalignment reaching 0.2 mm was measured for all cones but one, highlighting the benefit of systematizing this control before each patient SRS treatment.

CONCLUSION

An effective in-house QA program dedicated to SRS cones was developed to supplement the standard machine performance check on our mixed SRS platform. Specific geometric tests even include the smallest 4-mm cone, which could be of great interest for future clinical indications such as the radiosurgery of functional disorders.

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.

Investigation of Elekta AQUA software for kilovoltage to megavoltage radiation isocenter coincidence

Elekta AQUA v2.02 software (Gantry Runout isocenter test) was investigated as a tool for verification of kilovoltage to megavoltage gantry radiation isocenter coincidence. AQUA reported megavoltage (6 MV) isocenter was independent of field size over the range 5 cm × 5 cm to 20 cm × 20 cm. For the 10 cm × 10 cm field size, standard deviation in AQUA reported 3D megavoltage (6MV) isocenter over ten consecutive deliveries was less than 0.04 mm for any axis. Compared to the full AQUA test delivery (Gantry Runout), the shorter AQUA test version (Gantry Runout short) gave a root mean square MV isocenter (± 1 SD) difference of 0.18 mm ± 0.08. Across 7 machines, root mean square differences between AQUA and PIPS PRO reported MV isocenter (for 6 MV and 6 MV FFF beams) was 0.1 mm ± 0.1 mm, with most of the difference observed in the gun-target (Y-axis). AQUA 6 MV isocentre position was offset to gantry relative to Elekta XVI customer acceptance test (CAT) workflow by 0.25 mm to 0.52 mm. For 6 MV FFF beams, AQUA reported an MV isocenter of between 0.27 and 0.39 mm offset to target relative to Elekta CAT.

Positioning accuracy of a single-isocenter multiple targets SRS treatment: A comparison between Varian TrueBeam CBCT and Brainlab ExacTrac

PURPOSE

This study compared the positioning accuracy between cone-beam CT (CBCT) and ExacTrac (ETX) for a single-isocenter multiple target stereotactic radiosurgery (SRS) on two TrueBeam STx systems.

METHODS

A single-isocenter treatment plan was simulated on an anthropomorphic head phantom with six spherical steel ball bearings (BBs). One of the BBs was chosen to be the isocenter. The five off-isocenter targets were located at various distances from the isocenter. MV portal images were generated to evaluate the deviations between the expected and the real center of the targets after CBCT and ETX positioning, respectively.

RESULTS

The evaluation of the positioning accuracy for the isocenter target showed that CBCT and ETX positioning provided comparable, sub-millimetric results. Deviations in positioning accuracy were also calculated for all other targets, also showing comparable results for CBCT and ETX. Moreover, our study showed that the deviation between CBCT and ETX positioning were in better agreement for TBSTx1 and deviated slightly higher on TBSTx2 (maximum: 1.23 mm at S/I direction), due to a less perfect alignment between the CBCT coordinate system and the ETX coordinate system on TBSTx2 compared to TBSTx1. This study also showed a correlation between the target positioning accuracy and the distance to the isocenter.

CONCLUSION

The positioning accuracy of ETX and CBCT for targets located at isocenter and off-isocenter locations was compared on two treatment machines and found comparable. Our study highlights the importance of a proper calibration procedure, to ensure correct alignment between the CBCT, ETX and machine coordinate systems.

Independent 6D quality assurance of stereotactic radiotherapy repositioning on linacs

PURPOSE

A high level of accuracy while positioning the patient is mandatory for frameless stereotactic radiotherapy (SRT), as large doses in multiple fractions can be delivered near organs at risk. The objective of this study is to propose an end-to-end quality assurance method to verify that submillimetre alignment can be achieved with stereotactic conventional linacs.

METHODS

We used a TrueBeam® linear accelerator equipped with a 6DOF robotic couch. The "ISO Cube" phantom was used with a homemade stand designed to generate known translational and rotational offsets. A reference CT scan was performed with straight alignment of the phantom. The procedure introduced 1.6° angular offset for the couch pitch and roll, at various gantry angles. The couch base was also moved between 0° and 270°. We compared the results with the daily machine performance check tests (MPC, Varian).

RESULTS

The mean isocentre size, MV and kV imager offsets were found to agree to within 0.1mm, 0.1mm and 0.3mm respectively, and were in close agreement between the methods. For a total four months data collection period, the mean deviation between requested and measured 6DOF couch shifts was 0.6mm and 0.2°. Errors on field size were smaller than 1mm for 97.7% of the 324 data points.

CONCLUSION

Results demonstrate that the linac equipped with a 6DOF robotic positioner and CBCT imaging satisfies requirements for SRT. Our methodology, based on a modified Winston-Lutz quality control, allowed us to quantitatively assess end-to-end accuracy of a linac in order to safely deliver SRT.

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.

Feasibility of virtual starshot analysis providing submillimeter radiation isocenter accuracy: A long-term multi-institutional analysis

PURPOSE

We developed a technique to calculate the offset between room lasers and the radiation isocenter using a digital Winston-Lutz (WL) test with a starshot technique. We have performed isocenter localization quality assurance (QA) with submillimeter accuracy for a long period. Here we evaluated the feasibility and accuracy of this virtual starshot (VS) analysis for isocenter localization QA.

METHODS

A 6-MV photon beam with a square multileaf collimator field was used to irradiate a WL sphere positioned at the intersection of the room lasers. Images were acquired using an electronic portal imaging device. A four-field WL test was performed, and the path of each beam was calculated from the offset between the beam and sphere. Virtual starshot analysis was used to analyze the radiation isocenter, which calculates the center of the beam paths by using a least-squares method, similar to the starshot analysis. Then, eight coplanar and 12 noncoplanar beams were irradiated to evaluate isocenter localization accuracy.

RESULTS

Several VS analyses, using different WL spheres, were performed at three institutions, and the calculated accuracies were within 0.1 mm at all institutions. Long-term analysis showed that the isocenter localization accuracy was appropriately managed with three-dimensional accuracy within ± 0.5 mm for 90 months after the first laser adjustments. The offset between each beam and the room laser was within 0.6 mm and within 1.0 mm for eight coplanar and 12 noncoplanar beams, respectively, for 90 months. Cone-beam computed tomography images, acquired after verification beams, showed that the offset between the radiation isocenter and the imaging center was within 0.66 mm for 90 months. The isocenter localization accuracy within 1 mm was kept for long period at other four institutions.

CONCLUSIONS

Long-term analysis showed the feasibility of VS analysis for isocenter localization QA, including room laser re-alignment, noncoplanar irradiation verification, and image guidance accuracy.

A simulation-based method for evaluating geometric tests of a linac c-arm in quality control in radiotherapy

Purpose

Assessment of the accuracy of geometric tests of a linac used in external beam therapy is crucial for ensuring precise dose delivery. In this paper, a new simulation‐based method for assessing accuracy of such geometric tests is proposed and evaluated on a set of testing procedures.

Methods

Linac geometry testing methods used in this study are based on an established design of a two‐module phantom. Electronic portal imaging device (EPID) images of fiducial balls contained in these modules can be used to automatically reconstruct linac geometry. The projection of the phantom modules fiducial balls onto the EPID detector plane is simulated for assumed nominal geometry of a linac. Then, random errors are added to the coordinates of the projections of the centers of the fiducial balls and the linac geometry is reconstructed from these data.

Results

Reconstruction is performed for a set of geometric test designs and it is shown how the dispersion of the reconstructed values of geometric parameters depends on the design of a geometric test. Assuming realistic accuracy of EPID image analysis, it is shown that for selected testing plans the reconstruction accuracy of geometric parameters can be significantly better than commonly used action thresholds for these parameters.

Conclusions

Proposed solution has the potential to improve geometric testing design and practice. It is an important part of a fully automated geometric testing solution.

Automatic detection of graticule isocenter and scale from kV and MV images

PURPOSE

To automate the detection of isocenter and scale of the mechanical graticule on kilo-voltage (kV) or mega-voltage (MV) films or electronic portal imaging device (EPID) images.

METHODS

We developed a robust image processing approach to automatically detect isocenter and scale of mechanical graticule from digitized kV or MV films and EPID images. After a series of preprocessing steps applied to the digital images, a combination of Hough transform and Radon transform was performed to detect the graticule axes and isocenter. The magnification of the graticule was automatically detected by solving an optimization problem using golden section search and parabolic interpolation algorithm. Tick marks of the graticule were then determined by extending from isocenter along the graticule axes with multiples of the magnification value. This approach was validated using 23 kV films, 26 MV films, and 91 EPID images in different anatomical sites (head-and-neck, thorax, and pelvis). Accuracy was measured by comparing computer detected results with manually selected results.

RESULTS

The proposed approach was robust for kV and MV films of varying image quality. The isocenter was detected within 1 mm for 98% of the images. The exceptions were three kV films where the graticule was not actually visible. Of all images with correct isocenter detection, 99% had a magnification detection error less than 1% and tick mark detection error less than 1 mm, with the exception of 1 kV film (magnification error: 3.17%; tick mark error: 1.29 mm) and 1 MV film (magnification error: 0.45%; tick mark error: 1.11 mm).

CONCLUSION

We developed an approach to robustly and automatically detect graticule isocenter and scale from two-dimensionla (2D) kV and MV films. This is a first step toward automated treatment planning based on 2D x-ray images.

A novel phantom and procedure providing submillimeter accuracy in daily QA tests of accelerators used for stereotactic radiosurgery*

Stereotactic radiosurgery (SRS) places great demands on spatial accuracy. Steel BBs used as markers in quality assurance (QA) phantoms are clearly visible in MV and planar kV images, but artifacts compromise cone‐beam CT (CBCT) isocenter localization. The purpose of this work was to develop a QA phantom for measuring with sub‐mm accuracy isocenter congruence of planar kV, MV, and CBCT imaging systems and to design a practical QA procedure that includes daily Winston‐Lutz (WL) tests and does not require computer aid. The salient feature of the phantom (Universal Alignment Ball (UAB)) is a novel marker for precisely localizing isocenters of CBCT, planar kV, and MV beams. It consists of a 25.4 mm diameter sphere of polymethylmetacrylate (PMMA) containing a concentric 6.35 mm diameter tungsten carbide ball. The large density difference between PMMA and the polystyrene foam in which the PMMA sphere is embedded yields a sharp image of the sphere for accurate CBCT registration. The tungsten carbide ball serves in finding isocenter in planar kV and MV images and in doing WL tests. With the aid of the UAB, CBCT isocenter was located within 0.10±0.05 mm of its true positon, and MV isocenter was pinpointed in planar images to within 0.06±0.04 mm. In clinical morning QA tests extending over an 18 months period the UAB consistently yielded measurements with sub‐mm accuracy. The average distance between isocenter defined by orthogonal kV images and CBCT measured 0.16±0.12 mm. In WL tests the central ray of anterior beams defined by a 1.5×1.5 cm2 MLC field agreed with CBCT isocenter within 0.03±0.14 mm in the lateral direction and within 0.10±0.19 mm in the longitudinal direction. Lateral MV beams approached CBCT isocenter within 0.00±0.11 mm in the vertical direction and within ‐0.14±0.15 mm longitudinally. It took therapists about 10 min to do the tests. The novel QA phantom allows pinpointing CBCT and MV isocenter positions to better than 0.2 mm, using visual image registration. Under CBCT guidance, MLC‐defined beams are deliverable with sub‐mm spatial accuracy. The QA procedure is practical for daily tests by therapists.

PACS number(s): 87.53.Ly, 87.56.Fc

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.

Note: A simple image processing based fiducial auto-alignment method for sample registration

A simple method for the location and auto-alignment of sample fiducials for sample registration using widely available MATLAB/LabVIEW software is demonstrated. The method is robust, easily implemented, and applicable to a wide variety of experiment types for improved reproducibility and increased setup speed. The software uses image processing to locate and measure the diameter and center point of circular fiducials for distance self-calibration and iterative alignment and can be used with most imaging systems. The method is demonstrated to be fast and reliable in locating and aligning sample fiducials, provided here by a nanofabricated array, with accuracy within the optical resolution of the imaging system. The software was further demonstrated to register, load, and sample the dynamically wetted array.

Automated analysis of dynamic behavior of single cells in picoliter droplets

We present a droplet-based microfluidic platform to automatically track and characterize the behavior of single cells over time. This high-throughput assay allows encapsulation of single cells in micro-droplets and traps intact droplets in arrays of miniature wells on a PDMS-glass chip. Automated time-lapse fluorescence imaging and image analysis of the incubated droplets on the chip allows the determination of the viability of individual cells over time. In order to automatically track the droplets containing cells, we developed a simple method based on circular Hough transform to identify droplets in images and quantify the number of live and dead cells in each droplet. Here, we studied the viability of several hundred single isolated HEK293T cells over time and demonstrated a high survival rate of the encapsulated cells for up to 11 hours. The presented platform has a wide range of potential applications for single cell analysis, e.g. monitoring heterogeneity of drug action over time and rapidly assessing the transient behavior of single cells under various conditions and treatments in vitro.

Quantifying the gantry sag on linear accelerators and introducing an MLC-based compensation strategy

PURPOSE

Gantry sag is one of the well-known sources of mechanical imperfections that compromise the spatial accuracy of radiation dose delivery. The objectives of this study were to quantify the gantry sag on multiple linear accelerators (linacs), to investigate a multileaf collimator (MLC)-based strategy to compensate for gantry sag, and to verify the gantry sag and its compensation with film measurements.

METHODS

The authors used the Winston-Lutz method to measure gantry sag on three Varian linacs. A ball bearing phantom was imaged with megavolt radiation fields at 10° gantry angle intervals. The images recorded with an electronic portal imaging device were analyzed to derive the radiation isocenter and the gantry sag, that is, the superior-inferior wobble of the radiation field center, as a function of the gantry angle. The authors then attempted to compensate for the gantry sag by applying a gantry angle-specific correction to the MLC leaf positions. The gantry sag and its compensation were independently verified using film measurements.

RESULTS

Gantry sag was reproducible over a six-month measurement period. The maximum gantry sag was found to vary from 0.7 to 1.0 mm, depending on the linac and the collimator angle. The radiation field center moved inferiorly (i.e., away from the gantry) when the gantry was rotated from 0° to 180°. After the MLC leaf position compensation was applied at 90° collimator angle, the maximum gantry sag was reduced to <0.2 mm. The film measurements at gantry angles of 0° and 180° verified the inferior shift of the radiation fields and the effectiveness of MLC compensation.

CONCLUSIONS

The results indicate that gantry sag on a linac can be quantitatively measured using a simple phantom and an electronic portal imaging device. Reduction of gantry sag is feasible by applying a gantry angle-specific correction to MLC leaf positions at 90° collimator angle.

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.

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.

A quality assurance procedure to evaluate cone-beam CT image center congruence with the radiation isocenter of a linear accelerator

A quality assurance (QA) procedure was developed to evaluate the congruence between the cone‐beam computed tomography (CBCT) image center and the radiation isocenter on a Varian Trilogy linac. In contrast to the published QA procedures, this method did not require a ball bearing (BB) phantom to be placed exactly at the radiation isocenter through precalibrated room lasers or light field crosshairs. The only requirement was that the BB phantom be in a stationary position near the radiation isocenter during the image acquisition process. The radiation isocenter was determined with respect to the center of the BB using a Winston‐Lutz test. The CBCT image center was found to have excellent short‐term positional reproducibility (i.e., less than 0.1 mm of wobble in each of the x (lateral), y (vertical), and z (longitudinal) directions) in 10 consecutive acquisitions. Measured over a seven‐month period, the CBCT image center deviated from the radiation isocenter by 0.40±0.12mm(x),0.43±0.04mm(y), and 0.34±0.14mm(z). The z displacement of the 3D CBCT image center was highly correlated (ρ=0.997) with that of the 2D kV portal image center. The correlation coefficients in the x and y directions were poor (ρ=0.66 and ‐0.35, respectively). Systematic discrepancies were found between the CBCT image center and the 2D MV, kV portal image centers. For the linear accelerator studied, we detected a 0.8 mm discrepancy between the CBCT image center and the MV EPID image center in the anterior‐posterior direction. This discrepancy was demonstrated in a clinical case study where the patient was positioned with CBCT followed by MV portal verification. The results from the new QA procedure are useful for guiding high‐precision patient positioning in stereotactic body radiation therapy.

PACS number: 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.