Use of local noise power spectrum and wavelet analysis in quantitative image quality assurance for EPIDs.

PURPOSE: To investigate the use of local noise power spectrum (NPS) to characterize image noise and wavelet analysis to isolate defective pixels and inter-subpanel flat-fielding artifacts for quantitative quality assurance (QA) of electronic portal imaging devices (EPIDs). METHODS: A total of 93 image sets including custom-made bar-pattern images and open exposure images were collected from four iViewGT a-Si EPID systems over three years. Global quantitative metrics such as modulation transform function (MTF), NPS, and detective quantum efficiency (DQE) were computed for each image set. Local NPS was also calculated for individual subpanels by sampling region of interests within each subpanel of the EPID. The 1D NPS, obtained by radially averaging the 2D NPS, was fitted to a power-law function. The r-square value of the linear regression analysis was used as a singular metric to characterize the noise properties of individual subpanels of the EPID. The sensitivity of the local NPS was first compared with the global quantitative metrics using historical image sets. It was then compared with two commonly used commercial QA systems with images collected after applying two different EPID calibration methods (single-level gain and multilevel gain). To detect isolated defective pixels and inter-subpanel flat-fielding artifacts, Haar wavelet transform was applied on the images. RESULTS: Global quantitative metrics including MTF, NPS, and DQE showed little change over the period of data collection. On the contrary, a strong correlation between the local NPS (r-square values) and the variation of the EPID noise condition was observed. The local NPS analysis indicated image quality improvement with the r-square values increased from 0.80 ± 0.03 (before calibration) to 0.85 ± 0.03 (after single-level gain calibration) and to 0.96 ± 0.03 (after multilevel gain calibration), while the commercial QA systems failed to distinguish the image quality improvement between the two calibration methods. With wavelet analysis, defective pixels and inter-subpanel flat-fielding artifacts were clearly identified as spikes after thresholding the inversely transformed images. CONCLUSIONS: The proposed local NPS (r-square values) showed superior sensitivity to the noise level variations of individual subpanels compared with global quantitative metrics such as MTF, NPS, and DQE. Wavelet analysis was effective in detecting isolated defective pixels and inter-subpanel flat-fielding artifacts. The proposed methods are promising for the early detection of imaging artifacts of EPIDs.

Monitoring ABC-assisted deep inspiration breath hold for left-sided breast radiotherapy with an optical tracking system.

PURPOSE: Recent knowledge on the effects of cardiac toxicity warrants greater precision for left-sided breast radiotherapy. Different breath-hold (BH) maneuvers (abdominal vs thoracic breathing) can lead to chest wall positional variations, even though the patient's tidal volume remains consistent. This study aims to investigate the feasibility of using optical tracking for real-time quality control of active breathing coordinator (ABC)-assisted deep inspiration BH (DIBH). METHODS: An in-house optical tracking system (OTS) was used to monitor ABC-assisted DIBH. The stability and localization accuracy of the OTS were assessed with a ball-bearing phantom. Seven patients with left-sided breast cancer were included. A free-breathing (FB) computed tomography (CT) scan and an ABC-assisted BH CT scan were acquired for each patient. The OTS tracked an infrared (IR) marker affixed over the patient's xiphoid process to measure the positional variation of each individual BH. Using the BH within which the CT scan was performed as the reference, the authors quantified intra- and interfraction BH variations for each patient. To estimate the dosimetric impact of BH variations, the authors studied the positional correlation between the marker and the left breast using the FB CT and BH CT scans. The positional variations of 860 BHs as measured by the OTS were retrospectively incorporated into the original treatment plans to evaluate their dosimetric impact on breast and cardiac organs [heart and left anterior descending (LAD) artery]. RESULTS: The stability and localization accuracy of the OTS was within 0.2 mm along each direction. The mean intrafraction variation among treatment BHs was less than 2.8 mm in all directions. Up to 12.6 mm anteroposterior undershoot, where the patient's chest wall displacement of a BH is less than that of a reference BH, was observed with averages of 4.4, 3.6, and 0.1 mm in the anteroposterior, craniocaudal, and mediolateral directions, respectively. A high positional correlation between the marker and the breast was found in the anteroposterior and craniocaudal directions with respective Pearson correlation values of 0.95 and 0.93, but no mediolateral correlation was found. Dosimetric impact of BH variations on breast coverage was negligible. However, the mean heart dose, mean LAD dose, and max LAD dose were estimated to increase from 1.4/7.4/18.6 Gy (planned) to 2.1/15.7/31.0 Gy (delivered), respectively. CONCLUSIONS: In ABC-assisted DIBH, large positional variation can occur in some patients, due to their different BH maneuvers. The authors' study has shown that OTS can be a valuable tool for real-time quality control of ABC-assisted DIBH.

A further investigation of the centroid-to-centroid method for stereotactic lung radiotherapy: a phantom study.

PURPOSE: Our previous study [B. Lu et al., "A patient alignment solution for lung SBRT setups based on a deformable registration technique," Med. Phys. 39(12), 7379-7389 (2012)] proposed a deformable-registration-based patient setup strategy called the centroid-to-centroid (CTC) method, which can perform an accurate alignment of internal-target-volume (ITV) centroids between averaged four-dimensional computed tomography and cone-beam computed tomography (CBCT) images. Scenarios with variations between CBCT and simulation CT caused by irregular breathing and/or tumor change were not specifically considered in the patient study [B. Lu et al., "A patient alignment solution for lung SBRT setups based on a deformable registration technique," Med. Phys. 39(12), 7379-7389 (2012)] due to the lack of both a sufficiently large patient data sample and a method of tumor tracking. The aim of this study is to thoroughly investigate and compare the impacts of breathing pattern and tumor change on both the CTC and the translation-only (T-only) gray-value mode strategies by employing a four-dimensional (4D) lung phantom. METHODS: A sophisticated anthropomorphic 4D phantom (CIRS Dynamic Thorax Phantom model 008) was employed to simulate all desired respiratory variations. The variation scenarios were classified into four groups: inspiration to expiration ratio (IE ratio) change, tumor trajectory change, tumor position change, tumor size change, and the combination of these changes. For each category the authors designed several scenarios to demonstrate the effects of different levels of breathing variation on both of the T-only and the CTC methods. Each scenario utilized 4DCT and CBCT scans. The ITV centroid alignment discrepancies for CTC and T-only were evaluated. The dose-volume-histograms (DVHs) of ITVs for two extreme cases were analyzed. RESULTS: Except for some extreme cases in the combined group, the accuracy of the CTC registration was about 2 mm for all cases for both the single and the combined scenarios. The performance of the CTC method was insensitive to region-of-registration (ROR) size selections, as suggested by the comparable accuracy between 1 and 2 cm expansions of the ROR selections for the method. The T-only method was suitable for some single scenarios, such as trajectory variation, position variation, and size variation. However, for combined scenarios and/or a large variation in the IE ratio, the T-only method failed to produce reasonable registration results (within 3 mm). The discrepancy was close to, or even greater than, 1 cm. In addition, unlike the CTC method, the T-only method was sensitive to the ROR size selection. The DVH analysis suggested that a large ITV to PTV margin should be considered if a breathing pattern variation is observed. CONCLUSIONS: The phantom study demonstrated that the CTC method was reliable for scenarios in which breathing pattern variation was involved. The T-only gray value method worked for some scenarios, but not for scenarios that involved an IE ratio variation. For scenarios involving position variation, the T-only method worked only with a careful selection of the ROR, whereas the CTC method was independent of ROR size as long as the ITVs were included in the ROR. One indication of the dose consequence analysis was that a large ITV to PTV margin should be considered if a breathing pattern variation is observed.

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.

Feasibility of using two-dimensional array dosimeter for in vivo dose reconstruction via transit dosimetry.

Two-dimensional array dosimeters are commonly used to perform pretreatment quality assurance procedures, which makes them highly desirable for measuring transit fluences for in vivo dose reconstruction. The purpose of this study was to determine if an in vivo dose reconstruction via transit dosimetry using a 2D array dosimeter was possible. To test the accuracy of measuring transit dose distribution using a 2D array dosimeter, we evaluated it against the measurements made using ionization chamber and radiochromic film (RCF) profiles for various air gap distances (distance from the exit side of the solid water slabs to the detector distance; 0 cm, 30 cm, 40 cm, 50 cm, and 60 cm) and solid water slab thicknesses (10 cm and 20 cm). The backprojection dose reconstruction algorithm was described and evaluated. The agreement between the ionization chamber and RCF profiles for the transit dose distribution measurements ranged from -0.2% ~ 4.0% (average 1.79%). Using the backprojection dose reconstruction algorithm, we found that, of the six conformal fields, four had a 100% gamma index passing rate (3%/3 mm gamma index criteria), and two had gamma index passing rates of 99.4% and 99.6%. Of the five IMRT fields, three had a 100% gamma index passing rate, and two had gamma index passing rates of 99.6% and 98.8%. It was found that a 2D array dosimeter could be used for backprojection dose reconstruction for in vivo dosimetry.

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.

Characterization of a real-time surface image-guided stereotactic positioning system.

PURPOSE: The AlignRT3C system is an image-guided stereotactic positioning system (IGSPS) that provides real-time target localization. This study involves the first use of this system with three camera pods. The authors have evaluated its localization accuracy and tracking ability using a cone-beam computed tomography (CBCT) system and an optical tracking system in a clinical setting. METHODS: A modified Rando head-and-neck phantom and five patients receiving intracranial stereotactic radiotherapy (SRT) were used to evaluate the calibration, registration, and position-tracking accuracies of the AlignRT3C system and to study surface reconstruction uncertainties, including the effects due to interfractional and intrafractional motion, skin tone, room light level, camera temperature, and image registration region of interest selection. System accuracy was validated through comparison with the Elekta kV CBCT system (XVI) and the Varian frameless SonArray (FSA) optical tracking system. Surface-image data sets were acquired with the AlignRT3C daily for the evaluation of pretreatment and interfractional and intrafractional motion for each patient. Results for two different reference image sets, planning CT surface contours (CTS) and previously recorded AlignRT3C optical surface images (ARTS), are reported. RESULTS: The system origin displacements for the AlignRT3C and XVI systems agreed to within 1.3 mm and 0.7 degrees. Similar results were seen for AlignRT3C vs FSA. For the phantom displacements having couch angles of 0 degrees, those that utilized ART_S references resulted in a mean difference of 0.9 mm/0.4 degrees with respect to XVI and 0.3 mm/0.2 degrees with respect to FSA. For phantom displacements of more than +/- 10 mm and +/- 3 degrees, the maximum discrepancies between AlignRT and the XVI and FSA systems were 3.0 and 0.4 mm, respectively. For couch angles up to +/- 90 degrees, the mean (max.) difference between the AlignRT3C and FSA was 1.2 (2.3) mm/0.7 degrees (1.2 degrees). For all tests, the mean registration errors obtained using the CT_S references were approximately 1.3 mm/1.0 degrees larger than those obtained using the ART_S references. For the patient study, the mean differences in the pretreatment displacements were 0.3 mm/0.2 degrees between the AlignRT3C and XVI systems and 1.3 mm/1 degrees between the FSA and XVI systems. For noncoplanar treatments, interfractional motion displacements obtained using the ART_S and CT_S references resulted in 90th percentile differences within 2.1 mm/0.8 degrees and 3.3 mm/0.3 degrees, respectively, compared to the FSA system. Intrafractional displacements that were tracked for a maximum of 14 min were within 1 mm/1 degrees of those obtained with the FSA system. Uncertainties introduced by the bite-tray were as high as 3 mm/2 degrees for one patient. The combination of gantry, aSi detector panel, and x-ray tube blockage effects during the CBCT acquisition resulted in a registration error of approximately 3 mm. No skin-tone or surface deformation effects were seen with the limited patient sample. CONCLUSIONS: AlignRT3C can be used as a nonionizing IGSPS with accuracy comparable to current image/marker-based systems. IGSPS and CBCT can be combined for high-precision positioning without the need for patient-attached localization devices.

High-precision GAFCHROMIC EBT film-based absolute clinical dosimetry using a standard flatbed scanner without the use of a scanner non-uniformity correction.

To report a study of the use of GAFCHROMIC EBT radiochromic film (RCF) digitized with a commercially available flatbed document scanner for accurate and reliable all-purpose two-dimensional (2D) absolute dosimetry within a clinical environment. We used a simplified methodology that yields high-precision dosimetry measurements without significant postirradiation correction. The Epson Expression 1680 Professional scanner and the Epson Expression 10000XL scanner were used to digitize the films. Both scanners were retrofitted with light-diffusing glass to minimize the effects of Newton rings. Known doses were delivered to calibration films. Flat and wedge fields were irradiated with variable depth of solid water and 5 cm back scatter solid water. No particular scanner nonuniformity effect corrections or significant post-scan image processing were carried out. The profiles were compared with CC04 ionization chamber profiles. The depth dose distribution was measured at a source-to-surface distance (SSD) of 100 cm with a field size of 10 x 10 cm2. Additionally, 22 IMRT fields were measured and evaluated using gamma index analysis. The overall accuracy of RCF with respect to CC04 was found to be 2%-4%. The overall accuracy of RCF was determined using the absolute mean of difference for all flat and wedge field profiles. For clinical IMRT fields, both scanners showed an overall gamma index passing rate greater than 90%. This work demonstrated that EBT films, in conjunction with a commercially available flatbed scanner, can be used as an accurate and precise absolute dosimeter. Both scanners showed that no significant scanner nonuniformity correction is necessary for accurate absolute dosimetry using the EBT films for field sizes smaller than or equal to 15 x 15 cm2.

Quality assessment of frameless fractionated stereotactic radiotherapy using cone beam computed tomography.

PURPOSE: A quality assessment of intracranial stereotactic radiotherapy was performed using cone beam computed tomography (CBCT). Setup errors were analyzed for two groups of patients: (1) those who were positioned using a frameless SonArray (FSA) system and immobilized with a bite plate and thermoplastic (TP) mask (the bFSA group); and (2) those who were positioned by room laser and immobilized using a TP mask (the mLAS group). METHODS AND MATERIALS: A quality assurance phantom was used to study the system differences between FSA and CBCT. The quality assessment was performed using an Elekta Synergy imager (XVI) (Elekta Oncology Systems, Norcross, GA) and an On-Board Imager (OBI) (Varian Medical Systems, Palo Alto, CA) for 25 patients. For the first three fractions, and weekly thereafter, the FSA system was used for patient positioning, after which CBCT was performed to obtain setup errors. RESULTS: (1) Phantom tests: The mean differences in the isocenter displacements for the two systems was 1.2 ± 0.7 mm. No significant variances were seen between the XVI and OBI units (p~0.208). (2)Patient tests: The mean of the displacements between FSA and CBCT were independent of the CBCT system used; mean setup errors for the bFSA group were smaller (1.2 mm) than those of the mLAS group (3.2 mm) (p < 0.005). For the mLAS patients, the 90th percentile and the maximum rotational displacements were 3° and 5°, respectively. A 4-mm drift in setup accuracy occurred over the treatment course for 1 bFSA patient. CONCLUSIONS: System differences of less than 1 mm between CBCT and FSA were seen. Error regression was observed for the bFSA patients, using CBCT (up to 4 mm) during the treatment course. For the mLAS group, daily CBCT imaging was needed to obtain acceptable setup accuracies.

Optimal image-guidance scenario with cone-beam computed tomography in conventionally fractionated radiotherapy for lung tumors.

PURPOSE: To determine the residual setup errors of several image guidance scenarios, using cone-beam computed tomography (CBCT) in conventionally fractionated radiotherapy for lung tumors. METHODS: Thirteen lung cancer patients were treated with conventionally fractionated radiotherapy, using daily image guidance with CBCT, resulting in 389 CBCT scans which were registered to the planning scan using automated soft-tissue registration. Using the resulting daily alignment data, 4 imaging frequency scenarios were analyzed: (A) no imaging; (B) weekly imaging with a 3-mm threshold; (C) first 5 fractions imaged, then weekly imaging with a patient-specific threshold; and (D) imaging every other day. RESULTS: The systematic setup error (Sigma) was reduced with increasing frequency of imaging from 3.4 mm for no imaging to 1.0 mm for imaging every other day. Random setup error (sigma), however, varied little regardless of the frequency of imaging: 2.9, 3.0, 3.4, and 3.2 mm for scenarios A, B, C, and D, respectively. The setup margins required to account for the residual error of each imaging scenario were 1 to 1.6 cm for scenario A, 4 to 6 mm for scenarios B and C, and 4 to 5 mm for scenario D. As the residual error of daily CBCT was not included in this analysis, these margins compare with a margin of zero for daily CBCT. CONCLUSIONS: Daily image guidance is ideal as the setup margin can be reduced by about 5 mm versus a nondaily imaging scenario. However, if daily image guidance is not possible, there is little benefit in imaging more often than once a week.

Use of a line-pair resolution phantom for comprehensive quality assurance of electronic portal imaging devices based on fundamental imaging metrics.

Image guided radiation therapy solutions based on megavoltage computed tomography (MVCT) involve the extension of electronic portal imaging devices (EPIDs) from their traditional role of weekly localization imaging and planar dose mapping to volumetric imaging for 3D setup and dose verification. To sustain the potential advantages of MVCT, EPIDs are required to provide improved levels of portal image quality. Therefore, it is vital that the performance of EPIDs in clinical use is maintained at an optimal level through regular and rigorous quality assurance (QA). Traditionally, portal imaging QA has been carried out by imaging calibrated line-pair and contrast resolution phantoms and obtaining arbitrarily defined QA indices that are usually dependent on imaging conditions and merely indicate relative trends in imaging performance. They are not adequately sensitive to all aspects of image quality unlike fundamental imaging metrics such as the modulation transfer function (MTF), noise power spectrum (NPS), and detective quantum efficiency (DQE) that are widely used to characterize detector performance in radiographic imaging and would be ideal for QA purposes. However, due to the difficulty of performing conventional MTF measurements, they have not been used for routine clinical QA. The authors present a simple and quick QA methodology based on obtaining the MTF, NPS, and DQE of a megavoltage imager by imaging standard open fields and a bar-pattern QA phantom containing 2 mm thick tungsten line-pair bar resolution targets. Our bar-pattern based MTF measurement features a novel zero-frequency normalization scheme that eliminates normalization errors typically associated with traditional bar-pattern measurements at megavoltage x-ray energies. The bar-pattern QA phantom and open-field images are used in conjunction with an automated image analysis algorithm that quickly computes the MTF, NPS, and DQE of an EPID system. Our approach combines the fundamental advantages of linear systems metrics such as robustness, sensitivity across the full spatial frequency range of interest, and normalization to imaging conditions (magnification, system gain settings, and exposure), with the simplicity, ease, and speed of traditional phantom imaging. The algorithm was analyzed for accuracy and sensitivity by comparing with a commercial portal imaging QA method (PIPSPRO, Standard Imaging, Middleton, WI) on both first-generation lens-coupled and modern a-Si flat-panel based clinical EPID systems. The bar-pattern based QA measurements were found to be far more sensitive to even small levels of degradation in spatial resolution and noise. The bar-pattern based QA methodology offers a comprehensive image quality assessment tool suitable for both commissioning and routine EPID QA.

Evaluation of similarity measures for use in the intensity-based rigid 2D-3D registration for patient positioning in radiotherapy.

PURPOSE: Rigid 2D-3D registration is an alternative to 3D-3D registration for cases where largely bony anatomy can be used for patient positioning in external beam radiation therapy. In this article, the authors evaluated seven similarity measures for use in the intensity-based rigid 2D-3D registration using a variation in Skerl's similarity measure evaluation protocol. METHODS: The seven similarity measures are partitioned intensity uniformity, normalized mutual information (NMI), normalized cross correlation (NCC), entropy of the difference image, pattern intensity (PI), gradient correlation (GC), and gradient difference (GD). In contrast to traditional evaluation methods that rely on visual inspection or registration outcomes, the similarity measure evaluation protocol probes the transform parameter space and computes a number of similarity measure properties, which is objective and optimization method independent. The variation in protocol offers an improved property in the quantification of the capture range. The authors used this protocol to investigate the effects of the downsampling ratio, the region of interest, and the method of the digitally reconstructed radiograph (DRR) calculation [i.e., the incremental ray-tracing method implemented on a central processing unit (CPU) or the 3D texture rendering method implemented on a graphics processing unit (GPU)] on the performance of the similarity measures. The studies were carried out using both the kilovoltage (kV) and the megavoltage (MV) images of an anthropomorphic cranial phantom and the MV images of a head-and-neck cancer patient. RESULTS: Both the phantom and the patient studies showed the 2D-3D registration using the GPU-based DRR calculation yielded better robustness, while providing similar accuracy compared to the CPU-based calculation. The phantom study using kV imaging suggested that NCC has the best accuracy and robustness, but its slow function value change near the global maximum requires a stricter termination condition for an optimization method. The phantom study using MV imaging indicated that PI, GD, and GC have the best accuracy, while NCC and NMI have the best robustness. The clinical study using MV imaging showed that NCC and NMI have the best robustness. CONCLUSIONS: The authors evaluated the performance of seven similarity measures for use in 2D-3D image registration using the variation in Skerl's similarity measure evaluation protocol. The generalized methodology can be used to select the best similarity measures, determine the optimal or near optimal choice of parameter, and choose the appropriate registration strategy for the end user in his specific registration applications in medical imaging.

High performance computing for deformable image registration: towards a new paradigm in adaptive radiotherapy.

The advent of readily available temporal imaging or time series volumetric (4D) imaging has become an indispensable component of treatment planning and adaptive radiotherapy (ART) at many radiotherapy centers. Deformable image registration (DIR) is also used in other areas of medical imaging, including motion corrected image reconstruction. Due to long computation time, clinical applications of DIR in radiation therapy and elsewhere have been limited and consequently relegated to offline analysis. With the recent advances in hardware and software, graphics processing unit (GPU) based computing is an emerging technology for general purpose computation, including DIR, and is suitable for highly parallelized computing. However, traditional general purpose computation on the GPU is limited because the constraints of the available programming platforms. As well, compared to CPU programming, the GPU currently has reduced dedicated processor memory, which can limit the useful working data set for parallelized processing. We present an implementation of the demons algorithm using the NVIDIA 8800 GTX GPU and the new CUDA programming language. The GPU performance will be compared with single threading and multithreading CPU implementations on an Intel dual core 2.4 GHz CPU using the C programming language. CUDA provides a C-like language programming interface, and allows for direct access to the highly parallel compute units in the GPU. Comparisons for volumetric clinical lung images acquired using 4DCT were carried out. Computation time for 100 iterations in the range of 1.8-13.5 s was observed for the GPU with image size ranging from 2.0 x 10(6) to 14.2 x 10(6) pixels. The GPU registration was 55-61 times faster than the CPU for the single threading implementation, and 34-39 times faster for the multithreading implementation. For CPU based computing, the computational time generally has a linear dependence on image size for medical imaging data. Computational efficiency is characterized in terms of time per megapixels per iteration (TPMI) with units of seconds per megapixels per iteration (or spmi). For the demons algorithm, our CPU implementation yielded largely invariant values of TPMI. The mean TPMIs were 0.527 spmi and 0.335 spmi for the single threading and multithreading cases, respectively, with <2% variation over the considered image data range. For GPU computing, we achieved TPMI =0.00916 spmi with 3.7% variation, indicating optimized memory handling under CUDA. The paradigm of GPU based real-time DIR opens up a host of clinical applications for medical imaging.

Novel image registration quality evaluator (RQE) with an implementation for automated patient positioning in cranial radiation therapy.

In external beam radiation therapy, digitally reconstructed radiographs (DRRs) and portal images are used to verify patient setup based either on a visual comparison or, less frequently, with automated registration algorithms. A registration algorithm can be trapped in local optima due to irregularity of patient anatomy, image noise and artifacts, and/or out-of-plane shifts, resulting in an incorrect solution. Thus, human observation, which is subjective, is still required to check the registration result. We propose to use a novel image registration quality evaluator (RQE) to automatically identify misregistrations as part of an algorithm-based decision-making process for verification of patient positioning. A RQE, based on an adaptive pattern classifier, is generated from a pair of reference and target images to determine the acceptability of a registration solution given an optimization process. Here we applied our RQE to patient positioning for cranial radiation therapy. We constructed two RQEs-one for the evaluation of intramodal registrations (i.e., portal-portal); the other for intermodal registrations (i.e., portal-DRR). Mutual information, because of its high discriminatory ability compared with other measures (i.e., correlation coefficient and partitioned intensity uniformity), was chosen as the test function for both RQEs. We adopted 1 mm translation and 1 degree rotation as the maximal acceptable registration errors, reflecting desirable clinical setup tolerances for cranial radiation therapy. Receiver operating characteristic analysis was used to evaluate the performance of the RQE, including computations of sensitivity and specificity. The RQEs showed very good performance for both intramodal and intermodal registrations using simulated and phantom data. The sensitivity and the specificity were 0.973 and 0.936, respectively, for the intramodal RQE using phantom data. Whereas the sensitivity and the specificity were 0.961 and 0.758, respectively, for the intermodal RQE using phantom data. Phantom experiments also indicated our RQEs detected out-of-plane deviations exceeding 2.5 mm and 2.50. A preliminary retrospective clinical study of the RQE on cranial portal imaging also yielded good sensitivity > or = 0.857) and specificity (> or = 0.987). Clinical implementation of a RQE could potentially reduce the involvement of the human observer for routine patient positioning verification, while increasing setup accuracy and reducing setup verification time.

Analysis of the kinestatic charge detection system as a high detective quantum efficiency electronic portal imaging device.

Megavoltage x-ray imaging suffers from reduced image quality due to low differential x-ray attenuation and large Compton scatter compared with kilovoltage imaging. Notwithstanding this, electronic portal imaging devices (EPIDs) are now widely used in portal verification in radiotherapy as they offer significant advantages over film, including immediate digital imaging and superior contrast range. However video-camera-based EPIDs (VEPIDs) are limited by problems of low light collection efficiency and significant light scatter, leading to reduced contrast and spatial resolution. Indirect and direct detection-based flat-panel EPIDs have been developed to overcome these limitations. While flat-panel image quality has been reported to exceed that achieved with portal film, these systems have detective quantum efficiency (DQE) limited by the thin detection medium and are sensitive to radiation damage to peripheral read-out electronics. An alternative technology for high-quality portal imaging is presented here: kinesatic charge detection (KCD). The KCD is a scanning tri-electrode ion-chamber containing high-pressure noble gas (xenon at 100 atm) used in conjunction with a strip-collimated photon beam. The chamber is scanned across the patient, and an external electric field is used to regulate the cation drift velocity. By matching the scanning velocity with that of the cation (i.e., ion) drift velocity, the cations remain static in the object frame of reference, allowing temporal integration of the signal. The KCD offers several advantages as a portal imaging system. It has a thick detector geometry with an active detection depth of 6.1 cm, compared to the sub-millimeter thickness of the phosphor layer in conventional phosphor screens, leading to an order of magnitude advantage in quantum efficiency (>0.3). The unique principle of and the use of the scanning strip-collimated x-ray beam provide further integration of charges in time, reduced scatter, and a significantly reduced imaging dose, enhancing the imaging signal-to-noise ratio (SNR) and leading to high DQE. While thick detectors usually suffer from reduced spatial resolution, the KCD provides good spatial resolution due to high gas pressure that limits the spread of scattered electrons, and a strip-collimated beam that significantly reduces the inclusion of scatter in the imaging signal. A 10 cm wide small-field-of-view (SFOV) prototype of the KCD is presented with a complete analysis of its imaging performance. Measurements of modulation transfer function (MTF), noise power spectrum (NPS), and DQE were in good agreement with Monte Carlo simulations. Imaging signal loss from recombination within the KCD chamber was measured at different gas pressures, ion drift velocities, and strip-collimation widths. Image quality for the prototype KCD was also observed with anthropomorphic phantom imaging in comparison with various commercial and research portal imaging systems, including VEPID, flat-panel imager, and conventional and high contrast film systems. KCD-based imaging provided very good contrast and good spatial resolution at very low imaging dose (0.1 cGy per image). For the prototype KCD, measurements yielded DQE(0)=0.19 and DQE(1 cy/mm)=0.004.

Study of a prototype high quantum efficiency thick scintillation crystal video-electronic portal imaging device.

Image quality in portal imaging suffers significantly from the loss in contrast and spatial resolution that results from the excessive Compton scatter associated with megavoltage x rays. In addition, portal image quality is further reduced due to the poor quantum efficiency (QE) of current electronic portal imaging devices (EPIDs). Commercial video-camera-based EPIDs or VEPIDs that utilize a thin phosphor screen in conjunction with a metal buildup plate to convert the incident x rays to light suffer from reduced light production due to low QE (<2% for Eastman Kodak Lanex Fast-B). Flat-panel EPIDs that utilize the same luminescent screen along with an a-Si:H photodiode array provide improved image quality compared to VEPIDs, but they are expensive and can be susceptible to radiation damage to the peripheral electronics. In this article, we present a prototype VEPID system for high quality portal imaging at sub-monitor-unit (subMU) exposures based on a thick scintillation crystal (TSC) that acts as a high QE luminescent screen. The prototype TSC system utilizes a 12 mm thick transparent CsI(Tl) (thallium-activated cesium iodide) scintillator for QE=0.24, resulting in significantly higher light production compared to commercial phosphor screens. The 25 X 25 cm2 CsI(Tl) screen is coupled to a high spatial and contrast resolution Video-Optics plumbicon-tube camera system (1240 X 1024 pixels, 250 microm pixel width at isocenter, 12-bit ADC). As a proof-of-principle prototype, the TSC system with user-controlled camera target integration was adapted for use in an existing clinical gantry (Siemens BEAMVIEW(PLUS)) with the capability for online intratreatment fluoroscopy. Measurements of modulation transfer function (MTF) were conducted to characterize the TSC spatial resolution. The measured MTF along with measurements of the TSC noise power spectrum (NPS) were used to determine the system detective quantum efficiency (DQE). A theoretical expression of DQE(0) was developed to be used as a predictive model to propose improvements in the optics associated with the light detection. The prototype TSC provides DQE(0)=0.02 with its current imaging geometry, which is an order of magnitude greater than that for commercial VEPID systems and comparable to flat-panel imaging systems. Following optimization in the imaging geometry and the use of a high-end, cooled charge-coupled-device (CCD) camera system, the performance of the TSC is expected to improve even further. Based on our theoretical model, the expected DQE(0)=0.12 for the TSC system with the proposed improvements, which exceeds the performance of current flat-panel EPIDs. The prototype TSC provides high quality imaging even at subMU exposures (typical imaging dose is 0.2 MU per image), which offers the potential for daily patient localization imaging without increasing the weekly dose to the patient. Currently, the TSC is capable of limited frame-rate fluoroscopy for intratreatment visualization of patient motion at approximately 3 frames/second, since the achievable frame rate is significantly reduced by the limitations of the camera-control processor. With optimized processor control, the TSC is expected to be capable of intratreatment imaging exceeding 10 frames/second to monitor patient motion.

Theoretical analysis and experimental evaluation of a Csl(TI) based electronic portal imaging system.

This article discusses the design and analysis of a portal imaging system based on a thick transparent scintillator. A theoretical analysis using Monte Carlo simulation was performed to calculate the x-ray quantum detection efficiency (QDE), signal to noise ratio (SNR) and the zero frequency detective quantum efficiency [DQE(0)] of the system. A prototype electronic portal imaging device (EPID) was built, using a 12.7 mm thick, 20.32 cm diameter, Csl(Tl) scintillator, coupled to a liquid nitrogen cooled CCD TV camera. The system geometry of the prototype EPID was optimized to achieve high spatial resolution. The experimental evaluation of the prototype EPID involved the determination of contrast resolution, depth of focus, light scatter and mirror glare. Images of humanoid and contrast detail phantoms were acquired using the prototype EPID and were compared with those obtained using conventional and high contrast portal film and a commercial EPID. A theoretical analysis was also carried out for a proposed full field of view system using a large area, thinned CCD camera and a 12.7 mm thick CsI(TI) crystal. Results indicate that this proposed design could achieve DQE(0) levels up to 11%, due to its order of magnitude higher QDE compared to phosphor screen-metal plate based EPID designs, as well as significantly higher light collection compared to conventional TV camera based systems.

Verification of multileaf collimator leaf positions using an electronic portal imaging device.

An automated method is presented for determining individual leaf positions of the Siemens dual focus multileaf collimator (MLC) using the Siemens BEAMVIEW(PLUS) electronic portal imaging device (EPID). Leaf positions are computed with an error of 0.6 mm at one standard deviation (sigma) using separate computations of pixel dimensions, image distortion, and radiation center. The pixel dimensions are calculated by superimposing the film image of a graticule with the corresponding EPID image. A spatial correction is used to compensate for the optical distortions of the EPID, reducing the mean distortion from 3.5 pixels (uncorrected) per localized x-ray marker to 2 pixels (1 mm) for a rigid rotation and 1 pixel for a third degree polynomial warp. A correction for a nonuniform dosimetric response across the field of view of the EPID images is not necessary due to the sharp intensity gradients across leaf edges. The radiation center, calculated from the average of the geometric centers of a square field at 0 degrees and 180 degrees collimator angles, is independent of graticule placement error. Its measured location on the EPID image was stable to within 1 pixel based on 3 weeks of repeated extensions/retractions of the EPID. The MLC leaf positions determined from the EPID images agreed to within a pixel of the corresponding values measured using film and ionization chamber. Several edge detection algorithms were tested: contour, Sobel, Roberts, Prewitt, Laplace, morphological, and Canny. These agreed with each other to within < or = 1.2 pixels for the in-air EPID images. Using a test pattern, individual MLC leaves were found to be typically within 1 mm of the corresponding record-and-verify values, with a maximum difference of 1.8 mm, and standard deviations of <0.3 mm in the daily reproducibility. This method presents a fast, automatic, and accurate alternative to using film or a light field for the verification and calibration of the MLC.