MR-linac MLC positioning QA by digitally stitching dual double-exposed films.

PURPOSE: The picket fence (PF) test is highly recommended for multi-leaf collimator (MLC) quality assurance. However, since the electronic portal imaging device (EPID) on the Elekta Unity only covers a small area, it is not feasible to perform the PF test for the entire MLC. Here, we propose a technique for the PF test by stitching two double-exposed films. METHODS: Two EBT3 films were used to encompass the entire MLC, with each one covering one half of the area. Two fields were employed to apply double exposure: a PF pattern consisting of 11 2 mm wide pickets and a 2.84 cm x 22 cm open field. The edges of the open field defined by the diaphragms were used to correct film rotation as well as align them horizontally. The PF pattern was also measured with the EPID where the pickets were used to align the films vertically. Individual leaf positions were detected on the merged film for quantitative analysis. Various MLC positioning errors were introduced to evaluate the technique's sensitivity. RESULTS: The merged films covered 72 leaf pairs properly (four leaf pairs on both sides were outside the treatment couch). With the EPID, the leaf positioning accuracy was -0.02 ± 0.07 mm (maximum: 0.29 mm) and the picket width variation was 0.00 ± 0.03 mm (maximum: 0.11 mm); with the films, the position accuracy and width variation were -0.03 ± 0.13 mm (maximum: 0.80 mm) and 0.00 ± 0.13 mm (maximum: 0.74 mm), respectively. The EPID was able to detect errors of 0.5 mm or above with submillimeter accuracy; the films were only able to detect errors > 1.0 mm. CONCLUSION: We developed a quantitative technique for the PF test on the Elekta Unity. The merged films covered nearly the entire MLC leaf banks. The technique exhibited clinically acceptable accuracy and sensitivity to MLC positioning errors.

Reconstruction of volume averaging effect-free continuous photon beam profiles from discrete ionization chamber array measurements using a machine learning technique.

PURPOSE: The use of the ionization chamber array ICProfiler (ICP) is limited by its relatively poor detector spatial resolution and the inherent volume averaging effect (VAE). The purpose of this work is to study the feasibility of reconstructing VAE-free continuous photon beam profiles from ICP measurements with a machine learning technique. METHODS: In- and cross-plane photon beam profiles of a 6 MV beam from an Elekta linear accelerator, ranging from 2 × 2 to 10 × 10 cm2 at 1.5 cm, 5 cm, and 10 cm depth, were measured with an ICP. The discrete measurements were interpolated with a Makima method to obtain continuous beam profiles. Artificial neural networks (ANNs) were trained to restore the penumbra of the beam profiles. Plane-specific (in- and cr-plane) ANNs and a combined ANN were separately trained. The performance of the ANNs was evaluated using the penumbra width difference (PWD, the difference between the penumbra widths of the reconstructed and the reference profile). The plane-specific and the combined ANNs were compared to study the feasibility of using a single ANN for both in- and cross-plane. RESULTS: The profiles reconstructed with all the ANNs had excellent agreement with the reference. For in-plane, the ANNs reduced the PWD from 1.6 ± 0.7 mm at 1.5 cm depth to 0.1 ± 0.1 mm, from 1.8 ± 0.6 mm at 5.0 cm depth to 0.1 ± 0.1 mm, and from 2.4 ± 0.1 mm at 10.0 cm depth to 0.0 ± 0.0 mm; for cross-plane, the ANNs reduced the PWD from 1.2 ± 0.4 mm at 1.5 cm depth, 1.2 ± 0.3 mm at 5.0 cm depth, and 1.6 ± 0.1 mm at 10.0 cm depth, to 0.1 ± 0.1 mm. CONCLUSIONS: This study demonstrated the feasibility of using simple ANNs to reconstruct VAE-free continuous photon beam profiles from discrete ICP measurements. A combined ANN can restore the penumbra of in- and cross-plane beam profiles of various fields at different depths.

Evaluation of a neural network-based photon beam profile deconvolution method.

PURPOSE: The authors have previously shown the feasibility of using an artificial neural network (ANN) to eliminate the volume average effect (VAE) of scanning ionization chambers (ICs). The purpose of this work was to evaluate the method when applied to beams of different energies (6 and 10 MV) and modalities [flattened (FF) vs unflattened (FFF)], measured with ICs of various sizes. METHODS: The three-layer ANN extracted data from transverse photon beam profiles using a sliding window, and output deconvolved value corresponding to the location at the center of the window. Beam profiles of seven fields ranging from 2 × 2 to 10 × 10 cm2 at four depths (1.5, 5, 10 and 20 cm) were measured with three ICs (CC04, CC13, and FC65-P) and an EDGE diode detector for 6 MV FF and FFF. Similar data for the 10 MV FF beam was also collected with CC13 and EDGE. The EDGE-measured profiles were used as reference data to train and test the ANNs. Separate ANNs were trained by using the data of each beam energy and modality. Combined ANNs were also trained by combining data of different beam energies and/or modalities. The ANN's performance was quantified and compared by evaluating the penumbra width difference (PWD) between the deconvolved and reference profiles. RESULTS: Excellent agreement between the deconvolved and reference profiles was achieved with both separate and combined ANNs for all studied ICs, beam energies, beam modalities, and geometries. After deconvolution, the average PWD decreased from 1-3 mm to under 0.15 mm with separate ANNs and to under 0.20 mm with combined ANN. CONCLUSIONS: The ANN-based deconvolution method can be effectively applied to beams of different energies and modalities measured with ICs of various sizes. Separate ANNs yielded marginally better results than combined ANNs. An IC-specific, combined ANN can provide clinically acceptable results as long as the training data includes data of each beam energy and modality.

A self-checking treatment couch coordinate calculation system in radiotherapy.

PURPOSE: Traditionally, the treatment couch coordinates (TCCs) for patients undergoing radiotherapy can only be determined at the time of treatment, placing pressure on the treating therapists and leaving several pathways for errors such as wrong-site treatment or wrong treatment table shift from a reference point. The purpose of this work is to propose an accurate, robust, and streamlined system that calculates TCC in advance. METHODS: The proposed system combines the advantages of two different calculation methods that use an indexed immobilization device. The first method uses an array of reference ball bearings (BBs) embedded in the CT scanner's couch-top. To obtain the patient-specific TCC, the spatial offset of the treatment planning isocenter from the reference BB is used. The second method performs a calculation using the one-to-one mapping relationship between the CT scanner's DICOM (Digital Imaging and Communications in Medicine) coordinate system and the TCC system. Both methods use a reference point in the CT coordinate system to correlate a point in the TCC system to perform the coordinate transfer between the two systems. Both methods were used to calculate the TCC and the results were checked against each other, creating an integrated workflow via automated self-checking. The accuracy of the calculation system was retrospectively evaluated with 275 patients, where the actual treatment position determined with cone-beam CT was used as a reference. RESULTS: An efficient workflow transparent to the therapists at both CT simulation and treatment was created. It works with any indexed immobilization device and can be universally applied to all treatment sites. The two methods had comparable accuracy, with 95% of the calculations within 3 mm. The inter-fraction variation was within ± 1.0 cm for 95% of the coordinates across all the treatment sites. CONCLUSIONS: A robust, accurate, and streamlined system was implemented to calculate TCCs in advance. It eases the pressure on the treating therapists, reduces patient setup time, and enhances the patient safety by preventing setup errors.

Beam flatness modulation for a flattening filter free photon beam utilizing a novel direct leaf trajectory optimization model.

Flattening filter free (FFF) linear accelerators produce a fluence distribution that is forward peaked. Various dosimetric benefits, such as increased dose rate, reduced leakage and out of field dose has led to the growth of FFF technology in the clinic. The literature has suggested the idea of vendors offering dedicated FFF units where the flattening filter (FF) is removed completely and manipulating the beam to deliver conventional flat radiotherapy treatments. This work aims to develop an effective way to deliver modulated flat beam treatments, rather than utilizing a physical FF. This novel optimization model is an extension of the direct leaf trajectory optimization (DLTO) previously developed for volumetric modulated radiation therapy (VMAT) and is capable of accounting for all machine and multileaf collimator (MLC) dynamic delivery constraints, using a combination of linear constraints and a convex objective function. Furthermore, the tongue and groove (T&G) effect was also incorporated directly into our model without introducing nonlinearity to the constraints, nor nonconvexity to the objective function. The overall beam flatness, machine deliverability, and treatment time efficiency were assessed. Regular square fields, including field sizes of 10 × 10 cm2 to 40 × 40 cm2 were analyzed, as well as three clinical fields, and three arbitrary contours with "concave" features. Quantitative flatness was measured for all modulated FFF fields, and the results were comparable or better than their open FF counterparts, with the majority having a quantitative flatness of less than 3.0%. The modulated FFF beams, due to the included efficiency constraint, were able to achieve acceptable delivery time compared to their open FF counterpart. The results indicated that the dose uniformity and flatness for the modulated FFF beams optimized with the DLTO model can successfully match the uniformity and flatness of their conventional FF counterparts, and may even provide further benefit by taking advantage of the unique FFF beam characteristics.

Efficient independent planar dose calculation for FFF IMRT QA with a bivariate Gaussian source model.

The aim of this study is to perform a direct comparison of the source model for photon beams with and without flattening filter (FF) and to develop an efficient independent algorithm for planar dose calculation for FF-free (FFF) intensity-modulated radiotherapy (IMRT) quality assurance (QA). The source model consisted of a point source modeling the primary photons and extrafocal bivariate Gaussian functions modeling the head scatter, monitor chamber backscatter, and collimator exchange effect. The model parameters were obtained by minimizing the difference between the calculated and measured in-air output factors (Sc ). The fluence of IMRT beams was calculated from the source model using a backprojection and integration method. The off-axis ratio in FFF beams were modeled with a fourth degree polynomial. An analytical kernel consisting of the sum of three Gaussian functions was used to describe the dose deposition process. A convolution-based method was used to account for the ionization chamber volume averaging effect when commissioning the algorithm. The algorithm was validated by comparing the calculated planar dose distributions of FFF head-and-neck IMRT plans with measurements performed with a 2D diode array. Good agreement between the measured and calculated Sc was achieved for both FF beams (<0.25%) and FFF beams (<0.10%). The relative contribution of the head-scattered photons reduced by 34.7% for 6 MV and 49.3% for 10 MV due to the removal of the FF. Superior agreement between the calculated and measured dose distribution was also achieved for FFF IMRT. In the gamma comparison with a 2%/2 mm criterion, the average passing rate was 96.2 ± 1.9% for 6 MV FFF and 95.5 ± 2.6% for 10 MV FFF. The efficient independent planar dose calculation algorithm is easy to implement and can be valuable in FFF IMRT QA.

Image-domain shading correction for cone-beam CT without prior patient information.

In the era of high-precision radiotherapy, cone-beam CT (CBCT) is frequently utilized for on-board treatment guidance. However, CBCT images usually contain severe shading artifacts due to strong photon scatter from illumination of a large volume and non-optimized patient-specific data measurements, limiting the full clinical applications of CBCT. Many algorithms have been proposed to alleviate this problem by data correction on projections. Sophisticated methods have also been designed when prior patient information is available. Nevertheless, a standard, efficient, and effective approach with large applicability remains elusive for current clinical practice. In this work, we develop a novel algorithm for shading correction directly on CBCT images. Distinct from other image-domain correction methods, our approach does not rely on prior patient information or prior assumption of patient data. In CBCT, projection errors (mostly from scatter and non-ideal usage of bowtie filter) result in dominant low-frequency shading artifacts in image domain. In circular scan geometry, these artifacts often show global or local radial patterns. Hence, the raw CBCT images are first preprocessed into the polar coordinate system. Median filtering and polynomial fitting are applied on the transformed image to estimate the low-frequency shading artifacts (referred to as the bias field) angle-by-angle and slice-by-slice. The low-pass filtering process is done firstly along the angular direction and then the radial direction to preserve image contrast. The estimated bias field is then converted back to the Cartesian coordinate system, followed by 3D low-pass filtering to eliminate possible high-frequency components. The shading-corrected image is finally obtained as the uncorrected volume divided by the bias field. The proposed algorithm was evaluated on CBCT images of a pelvis patient and a head patient. Mean CT number values and spatial non-uniformity on the reconstructed images were used as image quality metrics. Within selected regions of interest, the average CT number error was reduced from around 300 HU to 42 and 38 HU, and the spatial nonuniformity error was reduced from above 17.5% to 2.1% and 1.7% for the pelvis and the head patients, respectively. As our method suppresses only low-frequency shading artifacts, patient anatomy and contrast were retained in the corrected images for both cases. Our shading correction algorithm on CBCT images offers several advantages. It has a high efficiency, since it is deterministic and directly operates on the reconstructed images. It requires no prior information or assumptions, which not only achieves the merits of CBCT-based treatment monitoring by retaining the patient anatomy, but also facilitates its clinical use as an efficient image-correction solution.

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.

Prevention of gross setup errors in radiotherapy with an efficient automatic patient safety system.

Treatment of the wrong body part due to incorrect setup is among the leading types of errors in radiotherapy. The purpose of this paper is to report an efficient automatic patient safety system (PSS) to prevent gross setup errors. The system consists of a pair of charge-coupled device (CCD) cameras mounted in treatment room, a single infrared reflective marker (IRRM) affixed on patient or immobilization device, and a set of in-house developed software. Patients are CT scanned with a CT BB placed over their surface close to intended treatment site. Coordinates of the CT BB relative to treatment isocenter are used as reference for tracking. The CT BB is replaced with an IRRM before treatment starts. PSS evaluates setup accuracy by comparing real-time IRRM position with reference position. To automate system workflow, PSS synchronizes with the record-and-verify (R&V) system in real time and automatically loads in reference data for patient under treatment. Special IRRMs, which can permanently stick to patient face mask or body mold throughout the course of treatment, were designed to minimize therapist's workload. Accuracy of the system was examined on an anthropomorphic phantom with a designed end-to-end test. Its performance was also evaluated on head and neck as well as abdominalpelvic patients using cone-beam CT (CBCT) as standard. The PSS system achieved a seamless clinic workflow by synchronizing with the R&V system. By permanently mounting specially designed IRRMs on patient immobilization devices, therapist intervention is eliminated or minimized. Overall results showed that the PSS system has sufficient accuracy to catch gross setup errors greater than 1 cm in real time. An efficient automatic PSS with sufficient accuracy has been developed to prevent gross setup errors in radiotherapy. The system can be applied to all treatment sites for independent positioning verification. It can be an ideal complement to complex image-guidance systems due to its advantages of continuous tracking ability, no radiation dose, and fully automated clinic workflow.

Performance assessment of variant UNet-based deep-learning dose engines for MR-Linac-based prostate IMRT plans.

Objective. UNet-based deep-learning (DL) architectures are promising dose engines for traditional linear accelerator (Linac) models. Current UNet-based engines, however, were designed differently with various strategies, making it challenging to fairly compare the results from different studies. The objective of this study is to thoroughly evaluate the performance of UNet-based models on magnetic-resonance (MR)-Linac-based intensity-modulated radiation therapy (IMRT) dose calculations.Approach. The UNet-based models, including the standard-UNet, cascaded-UNet, dense-dilated-UNet, residual-UNet, HD-UNet, and attention-aware-UNet, were implemented. The model input is patient CT and IMRT field dose in water, and the output is patient dose calculated by DL model. The reference dose was calculated by the Monaco Monte Carlo module. Twenty training and ten test cases of prostate patients were included. The accuracy of the DL-calculated doses was measured using gamma analysis, and the calculation efficiency was evaluated by inference time.Results. All the studied models effectively corrected low-accuracy doses in water to high-accuracy patient doses in a magnetic field. The gamma passing rates between reference and DL-calculated doses were over 86% (1%/1 mm), 98% (2%/2 mm), and 99% (3%/3 mm) for all the models. The inference times ranged from 0.03 (graphics processing unit) to 7.5 (central processing unit) seconds. Each model demonstrated different strengths in calculation accuracy and efficiency; Res-UNet achieved the highest accuracy, HD-UNet offered high accuracy with the fewest parameters but the longest inference, dense-dilated-UNet was consistently accurate regardless of model levels, standard-UNet had the shortest inference but relatively lower accuracy, and the others showed average performance. Therefore, the best-performing model would depend on the specific clinical needs and available computational resources.Significance. The feasibility of using common UNet-based models for MR-Linac-based dose calculations has been explored in this study. By using the same model input type, patient training data, and computing environment, a fair assessment of the models' performance was present.

A patient alignment solution for lung SBRT setups based on a deformable registration technique.

PURPOSE: In this work, the authors propose a novel registration strategy for translation-only correction scenarios of lung stereotactic body radiation therapy setups, which can achieve optimal dose coverage for tumors as well as preserve the consistency of registrations with minimal human interference. METHODS: The proposed solution (centroid-to-centroidor CTC solution) uses the average four-dimensional CT (A4DCT) as the reference CT. The cone-beam CT (CBCT) is deformed to acquire a new centroid for the internal target volume (ITV) on the CBCT. The registration is then accomplished by simply aligning the centroids of the ITVs between the A4DCT and the CBCT. Sixty-seven cases using 64 patients (each case is associated with separate isocenters) have been investigated with the CTC method and compared with the conventional gray-value (G) mode and bone (B) mode registration methods. Dosimetric effects among the tree methods were demonstrated by 18 selected cases. The uncertainty of the CTC method has also been studied. RESULTS: The registration results demonstrate the superiority of the CTC method over the other two methods. The differences in the D99 and D95 ITV dose coverage between the CTC method and the original plan is small (within 5%) for all of the selected cases except for one for which the tumor presented significant growth during the period between the CT scan and the treatment. Meanwhile, the dose coverage differences between the original plan and the registration results using either the B or G method are significant, as tumor positions varied dramatically, relative to the rib cage, from their positions on the original CT. The largest differences between the D99 and D95 dose coverage of the ITV using the B or G method versus the original plan are as high as 50%. The D20 differences between any of the methods versus the original plan are all less than 2%. CONCLUSIONS: The CTC method can generate optimal dose coverage to tumors with much better consistency compared with either the G or B method, and it is especially useful when the tumor position varies greatly from its position on the original CT, relative to the rib cage.

An ultra-fast deep-learning-based dose engine for prostate VMAT via knowledge distillation framework with limited patient data.

Objective. Deep-learning (DL)-based dose engines have been developed to alleviate the intrinsic compromise between the calculation accuracy and efficiency of the traditional dose calculation algorithms. However, current DL-based engines typically possess high computational complexity and require powerful computing devices. Therefore, to mitigate their computational burdens and broaden their applicability to a clinical setting where resource-limited devices are available, we proposed a compact dose engine via knowledge distillation (KD) framework that offers an ultra-fast calculation speed with high accuracy for prostate Volumetric Modulated Arc Therapy (VMAT).Approach. The KD framework contains two sub-models: a large pre-trained teacher and a small to-be-trained student. The student receives knowledge transferred from the teacher for better generalization. The trained student serves as the final engine for dose calculation. The model input is patient computed tomography and VMAT dose in water, and the output is DL-calculated patient dose. The ground-truth \dose was computed by the Monte Carlo module of the Monaco treatment planning system. Twenty and ten prostate cases were included for model training and assessment, respectively. The model's performance (teacher/student/student-only) was evaluated by Gamma analysis and inference efficiency.Main results. The dosimetric comparisons (input/DL-calculated/ground-truth doses) suggest that the proposed engine can effectively convert low-accuracy doses in water to high-accuracy patient doses. The Gamma passing rate (2%/2 mm, 10% threshold) between the DL-calculated and ground-truth doses was 98.64 ± 0.62% (teacher), 98.13 ± 0.76% (student), and 96.95 ± 1.02% (student-only). The inference time was 16 milliseconds (teacher) and 11 milliseconds (student/student-only) using a graphics processing unit device, while it was 936 milliseconds (teacher) and 374 milliseconds (student/student-only) using a central processing unit device.Significance. With the KD framework, a compact dose engine can achieve comparable accuracy to that of a larger one. Its compact size reduces the computational burdens and computing device requirements, and thus such an engine can be more clinically applicable.

A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling.

PURPOSE: Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect. METHODS AND MATERIALS: 6 MV 1 × 1 m m 2 $1 \times 1{\rm{\;m}}{{\rm{m}}^2}$ pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 c m 2 $40 \times 40{\rm{\;c}}{{\rm{m}}^2}$ fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using fast Fourier transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively. RESULTS: The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for 10 VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm. CONCLUSIONS: We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA.

An approach for online evaluations of dose consequences caused by small rotational setup errors in intracranial stereotactic radiation therapy.

PURPOSE: The purpose of this work is to investigate the impact of small rotational errors on the magnitudes and distributions of spatial dose variations for intracranial stereotactic radiotherapy (SRT) treatment setups, and to assess the feasibility of using the original dose map overlaid with rotated contours (ODMORC) method as a fast, online evaluation tool to estimate dose changes (using DVHs) to clinical target volumes (CTVs) and organs-at-risks (OARs) caused by small rotational setup errors. METHODS: Fifteen intracranial SRT cases treated with either three-dimensional conformal radiation therapy (3DCRT) or intensity-modulated radiation therapy (IMRT) techniques were chosen for the study. Selected cases have a variety of anatomical dimensions and pathologies. Angles of ±3° and ±5° in all directions were selected to simulate the rotational errors. Dose variations in different regions of the brain, CTVs, and OARs were evaluated to illustrate the various spatial effects of dose differences before and after rotations. DVHs accounting for rotations that were recomputed by the treatment planning system (TPS) and those generated by the ODMORC method were compared. A framework of a fast algorithm for multicontour rotation implemented by ODMORC is introduced as well. RESULTS: The average values of relative dose variations between original dose and recomputed dose accounting for rotations were greater than 4.0% and 10.0% in absolute mean and in standard deviation, respectively, at the skull and adjacent regions for all cases. They were less than 1.0% and 2.5% in absolute mean and in standard deviation, respectively, for dose points 3 mm away from the skull. The results indicated that spatial dose to any part of the brain organs or tumors separated from the skull or head surface would be relatively stable before and after rotations. Statistical data of CTVs and OARs indicate the lens and cochleas have the large dose variations before and after rotations, whereas the remaining ROIs have insignificant dose differences. DVH comparisons suggest that the ODMORC method is able to estimate the DVH of CTVs fairly accurately (within 1.5% of relative dose differences for evaluation volumes). The results also show that most of the OARs including the brain stem, spinal cord, chiasm, hippocampuses, optic nerves, and retinas, which were relatively distal from the skull and surface, had good agreement (within 2.0% of relative dose differences for 0.1 cc of the volumes ) between the ODMORC method and the recomputation, whereas OARs more proximate to the bone-tissue interface or surface, such as the lenses and cochlea, had larger dose variations (greater than 5.0%) for some cases due to the incapability of the ODMORC to account for scatter contribution variations proximate to interfaces and intrinsic dose calculation uncertainties for ROIs with small volumes. CONCLUSIONS: The ODMORC method can be implemented as an online evaluation system for rotation-induced dose changes of CTVs and most OARs and for other related dose consequence analyses.

Dosimetric consequences of rotational setup errors with direct simulation in a treatment planning system for fractionated stereotactic radiotherapy.

The purpose was to determine dose-delivery errors resulting from systematic rotational setup errors for fractionated stereotactic radiotherapy using direct simulation in a treatment planning system. Ten patients with brain tumors who received intensity-modulated radiotherapy had dose distributions re-evaluated to assess the impact of systematic rotational setup errors. The dosimetric effect of rotational setup errors was simulated by rotating images and contours using a 3 by 3 rotational matrix. Combined rotational errors of ± 1°, ± 3°, ± 5° and ± 7° and residual translation errors of 1 mm along each axis were simulated. Dosimetric effects of the rotated images were evaluated by recomputing dose distributions and compared with the original plan. The mean volume of CTV that received the prescription dose decreased from 99.3% ± 0.5% (original) to 98.6% ± 1.6% (± 1°), 97.0% ± 2.0% (± 3°), 93.1% ± 4.6% (± 5°), and 87.8% ± 14.2% (± 7°). Minimal changes in the cold and hot spots were seen in the CTV. In general, the increase in the volumes of the organs at risk (OARs) receiving the tolerance doses was small and did not exceed the tolerance, except for cases where the OARs were in close proximity to the PTV. For intracranial tumors treated with IMRT with a CTV-to-PTV margin of 3 mm, rotational setup errors of 3° or less didn't decrease the CTV coverage to less than 95% in most cases. However, for large targets with irregular or elliptical shapes, the target coverage decreased significantly as rotational errors of 5° or more were present. Our results indicate that setup margins are warranted even in the absence of translational setup errors to account for rotational setup errors. Rotational setup errors should be evaluated carefully for clinical cases involving large tumor sizes and for targets with elliptical or irregular shape, as well as when isocenter is away from the center of the PTV or OARs are in close proximity to the target volumes.

An image quality comparison study between XVI and OBI CBCT systems.

The purpose of this study is to evaluate and compare image quality characteristics for two commonly used and commercially available CBCT systems: the X-ray Volumetric Imager and the On-Board Imager. A commonly used CATPHAN image quality phantom was used to measure various image quality parameters, namely, pixel value stability and accuracy, noise, contrast to noise ratio (CNR), high-contrast resolution, low contrast resolution and image uniformity. For the XVI unit, we evaluated the image quality for four manufacturer-supplied protocols as a function of mAs. For the OBI unit, we did the same for the full-fan and half-fan scanning modes, which were respectively used with the full bow-tie and half bow-tie filters. For XVI, the mean pixel values of regions of interest were found to generally decrease with increasing mAs for all protocols, while they were relatively stable with mAs for OBI. Noise was slightly lower on XVI and was seen to decrease with increasing mAs, while CNR increased with mAs for both systems. For XVI and OBI, the high-contrast resolution was approximately limited by the pixel resolution of the reconstructed image. On OBI images, up to 6 and 5 discs of 1% and 0.5% contrast, respectively, were visible for a high mAs setting using the full-fan mode, while none of the discs were clearly visible on the XVI images for various mAs settings when the medium resolution reconstruction was used. In conclusion, image quality parameters for XVI and OBI have been quantified and compared for clinical protocols under various mAs settings. These results need to be viewed in the context of a recent study that reported the dose-mAs relationship for the two systems and found that OBI generally delivered higher imaging doses than XVI.

Feasibility study of performing IGRT system daily QA using a commercial QA device.

The purpose of this study was to investigate the feasibility of using a single QA device for comprehensive, efficient daily QA of a linear accelerator (Linac) and three image-guided stereotactic positioning systems (IGSPSs). The Sun Nuclear Daily QA 3 (DQA3) device was used to perform daily dosimetry and mechanical accuracy tests for an Elekta Linac, as well as daily image geometric and isocenter coincidence accuracy tests for three IGSPSs: the AlignRT surface imaging system; the frameless SonArray optical tracking System (FSA) and the Elekta kV CBCT. The DQA3 can also be used for couch positioning, repositioning, and rotational tests during the monthly QA. Based on phantom imaging, the Linac coordinate system determined using AlignRT was within 0.7 mm/0.6° of that of the CBCT system. The difference is attributable to the different calibration methods that are utilized for these two systems. The laser alignment was within 0.5 mm of the isocenter location determined with the three IGSPSs. The ODI constancy was ± 0.5 mm. For gantry and table angles of 0°, the mean isocenter displacement vectors determined using the three systems were within 0.7 mm and 0.6° of one another. Isocenter rotational offsets measured with the systems were all ≤ 0.5°. For photon and electron beams tested over a period of eight months, the output was verified to remain within 2%, energy variations were within 2%, and the symmetry and flatness were within 1%. The field size and light-radiation coincidence were within 1mm ± 1 mm. For dosimetry reproducibility, the standard deviation was within 0.2% for all tests and all energies, except for photon energy variation which was 0.6%. The total measurement time for all tasks took less than 15 minutes per QA session compared to 40 minutes with our previous procedure, which utilized an individual QA device for each IGSPS. The DQA3 can be used for accurate and efficient Linac and IGSPS daily QA. It shortens QA device setup time, eliminates errors introduced by changing phantoms to perform different tests, and streamlines the task of performing dosimetric checks.

AAPM Task Group 198 Report: An implementation guide for TG 142 quality assurance of medical accelerators.

The charges on this task group (TG) were as follows: (a) provide specific procedural guidelines for performing the tests recommended in TG 142; (b) provide estimate of the range of time, appropriate personnel, and qualifications necessary to complete the tests in TG 142; and (c) provide sample daily, weekly, monthly, or annual quality assurance (QA) forms. Many of the guidelines in this report are drawn from the literature and are included in the references. When literature was not available, specific test methods reflect the experiences of the TG members (e.g., a test method for door interlock is self-evident with no literature necessary). In other cases, the technology is so new that no literature for test methods was available. Given broad clinical adaptation of volumetric modulated arc therapy (VMAT), which is not a specific topic of TG 142, several tests and criteria specific to VMAT were added.

Calibration of a novel four-dimensional diode array.

PURPOSE: The aim of this work is to develop effective calibration methods for a novel fourdimensional (4D) diode array for pretreatment verification of intensity-modulated radiation therapy (IMRT) and rotational therapy. METHODS: A novel 4D diode array (ArcCHECK, Sun Nuclear, Melbourne, FL) was developed to meet the needs of appropriate and efficient quality assurance for IMRT and especially rotational radiotherapy. The diode array presents a consistent detector image in beam's eye view at arbitrary gantry angles due to isotropic arrangement of diodes in a three-dimensional (3D) cylindrical phantom. The 50 ms simultaneous update of all diodes on the detector array (fourth dimension) makes it capable of time-resolved beam delivery analysis with any rotational delivery techniques. The calibration procedure consisted of delivering and measuring a series of calibration beams with 5.8 degrees angular spacing surrounding the cylindrical diode array. Correction factors for diode intrinsic sensitivity and directional response dependence were derived from these measurements. A real-time algorithm to derive gantry angles based on the detector signal was developed to interpolate and apply the corresponding angular correction factors. RESULTS: The calibration was validated with ion chamber scanned beam profiles in a 3D water tank. Excellent agreement was observed between diode array measurement and treatment planning system calculation. The accuracy of the gantry angle derivation algorithm was within 1 degree which caused a less than 0.2% dosimetric uncertainty. CONCLUSIONS: With the proposed calibration method and the automatic gantry angle derivation algorithm, the 4D diode array achieved isotropic detector response and is suitable for both IMRT and rotational therapy pretreatment verification.

Characterization of a multi-axis ion chamber array.

PURPOSE: The aim of this work was to characterize a multi-axis ion chamber array (IC PROFILER; Sun Nuclear Corporation, Melbourne, FL, USA) that has the potential to simplify the acquisition of LINAC beam data. METHODS: The IC PROFILER (or panel) measurement response was characterized with respect to radiation beam properties, including dose, dose per pulse, pulse rate frequency (PRF), and energy. Panel properties were also studied, including detector-calibration stability, power-on time, backscatter dependence, and the panel's agreement with water tank measurements [profiles, fractional depth dose (FDD), and output factors]. RESULTS: The panel's relative deviation was typically within (+/-) 1% of an independent (or nominal) response for all properties that were tested. Notable results were (a) a detectable relative field shape change of approximately 1% with linear accelerator PRF changes; (b) a large range in backscatter thickness had a minimal effect on the measured dose distribution (typically less than 1%); (c) the error spread in profile comparison between the panel and scanning water tank (Blue Phantom, CC13; IBA Schwarzenbruck, DE) was approximately (+/-) 0.75%. CONCLUSIONS: The ability of the panel to accurately reproduce water tank profiles, FDDs, and output factors is an indication of its abilities as a dosimetry system. The benefits of using the panel versus a scanning water tank are less setup time and less error susceptibility. The same measurements (including device setup and breakdown) for both systems took 180 min with the water tank versus 30 min with the panel. The time-savings increase as the measurement load is increased.