AAPM Medical physics practice guideline 15.A: Peer review in clinical physics.

The American Association of Physicists in Medicine (AAPM) is a nonprofit professional society whose primary purposes are to advance the science, education, and professional practice of medical physics. The AAPM has more than 8000 members and is the principal organization of medical physicists in the United States. The AAPM will periodically define new practice guidelines for medical physics practice to help advance the science of medical physics and to improve the quality of service to patients throughout the United States. Existing medical physics practice guidelines will be reviewed for the purpose of revision or renewal, as appropriate, on their fifth anniversary or sooner. Each medical physics practice guideline represents a policy statement by the AAPM, has undergone a thorough consensus process in which it has been subjected to extensive review, and requires the approval of the Professional Council. The medical physics practice guidelines recognize that the safe and effective use of diagnostic and therapeutic radiology requires specific training, skills, and techniques, as described in each document. Reproduction or modification of the published practice guidelines and technical standards by those entities not providing these services is not authorized. The following terms are used in the AAPM practice guidelines: Must and Must Not: Used to indicate that adherence to the recommendation is considered necessary to conform to this practice guideline. While must is the term to be used in the guidelines, if an entity that adopts the guideline has shall as the preferred term, the AAPM considers that must and shall have the same meaning. Should and Should Not: Used to indicate a prudent practice to which exceptions may occasionally be made in appropriate circumstances.

Operational consistency of medical linear accelerators manufactured and commissioned in series.

The purpose of this study was to determine if medical linear accelerators (linac) produced by the same manufacturer exhibit operational consistency within their subsystems and components. Two linacs that were commissioned together and installed at the same facility were monitored. Each machine delivered a daily robust quality assurance (QA) irradiation. Linacs and their components operate consistently, but have different operational parameter levels even when produced by the same manufacturer and commissioned in series. These findings have implications on the feasibility of true clinical beam matching.

A model for preemptive maintenance of medical linear accelerators-predictive maintenance.

BACKGROUND: Unscheduled accelerator downtime can negatively impact the quality of life of patients during their struggle against cancer. Currently digital data accumulated in the accelerator system is not being exploited in a systematic manner to assist in more efficient deployment of service engineering resources. The purpose of this study is to develop an effective process for detecting unexpected deviations in accelerator system operating parameters and/or performance that predicts component failure or system dysfunction and allows maintenance to be performed prior to the actuation of interlocks. METHODS: The proposed predictive maintenance (PdM) model is as follows: 1) deliver a daily quality assurance (QA) treatment; 2) automatically transfer and interrogate the resulting log files; 3) once baselines are established, subject daily operating and performance values to statistical process control (SPC) analysis; 4) determine if any alarms have been triggered; and 5) alert facility and system service engineers. A robust volumetric modulated arc QA treatment is delivered to establish mean operating values and perform continuous sampling and monitoring using SPC methodology. Chart limits are calculated using a hybrid technique that includes the use of the standard SPC 3σ limits and an empirical factor based on the parameter/system specification. RESULTS: There are 7 accelerators currently under active surveillance. Currently 45 parameters plus each MLC leaf (120) are analyzed using Individual and Moving Range (I/MR) charts. The initial warning and alarm rule is as follows: warning (2 out of 3 consecutive values ≥ 2σ hybrid) and alarm (2 out of 3 consecutive values or 3 out of 5 consecutive values ≥ 3σ hybrid). A customized graphical user interface provides a means to review the SPC charts for each parameter and a visual color code to alert the reviewer of parameter status. Forty-five synthetic errors/changes were introduced to test the effectiveness of our initial chart limits. Forty-three of the forty-five errors (95.6 %) were detected in either the I or MR chart for each of the subsystems monitored. CONCLUSION: Our PdM model shows promise in providing a means for reducing unscheduled downtime. Long term monitoring will be required to establish the effectiveness of the model.

Initial investigation using statistical process control for quality control of accelerator beam steering.

BACKGROUND: This study seeks to increase clinical operational efficiency and accelerator beam consistency by retrospectively investigating the application of statistical process control (SPC) to linear accelerator beam steering parameters to determine the utility of such a methodology in detecting changes prior to equipment failure (interlocks actuated). METHODS: Steering coil currents (SCC) for the transverse and radial planes are set such that a reproducibly useful photon or electron beam is available. SCC are sampled and stored in the control console computer each day during the morning warm-up. The transverse and radial - positioning and angle SCC for photon beam energies were evaluated using average and range (Xbar-R) process control charts (PCC). The weekly average and range values (subgroup n = 5) for each steering coil were used to develop the PCC. SCC from September 2009 (annual calibration) until two weeks following a beam steering failure in June 2010 were evaluated. PCC limits were calculated using the first twenty subgroups. Appropriate action limits were developed using conventional SPC guidelines. RESULTS: PCC high-alarm action limit was set at 6 standard deviations from the mean. A value exceeding this limit would require beam scanning and evaluation by the physicist and engineer. Two low alarms were used to indicate negative trends. Alarms received following establishment of limits (week 20) are indicative of a non-random cause for deviation (Xbar chart) and/or an uncontrolled process (R chart). Transverse angle SCC for 6 MV and 15 MV indicated a high-alarm 90 and 108 days prior to equipment failure respectively. A downward trend in this parameter continued, with high-alarm, until failure. Transverse position and radial angle SCC for 6 and 15 MV indicated low-alarms starting as early as 124 and 116 days prior to failure, respectively. CONCLUSION: Radiotherapy clinical efficiency and accelerator beam consistency may be improved by instituting SPC methods to monitor the beam steering process and detect abnormal changes prior to equipment failure.PACS numbers: 87.55n, 87.55qr, 87.56bd.

Improved volumetric imaging in tomosynthesis using combined multiaxial sweeps.

This study explores the volumetric reconstruction fidelity attainable using tomosynthesis with a kV imaging system which has a unique ability to rotate isocentrically and with multiple degrees of mechanical freedom. More specifically, we seek to investigate volumetric reconstructions by combining multiple limited-angle rotational image acquisition sweeps. By comparing these reconstructed images with those of a CBCT reconstruction, we can gauge the volumetric fidelity of the reconstructions. In surgical situations, the described tomosynthesis-based system could provide high-quality volumetric imaging without requiring patient motion, even with rotational limitations present. Projections were acquired using the Digital Integrated Brachytherapy Unit, or IBU-D. A phantom was used which contained several spherical objects of varying contrast. Using image projections acquired during isocentric sweeps around the phantom, reconstructions were performed by filtered backprojection. For each image acquisition sweep configuration, a contrasting sphere is analyzed using two metrics and compared to a gold standard CBCT reconstruction. Since the intersection of a reconstructed sphere and an imaging plane is ideally a circle with an eccentricity of zero, the first metric presented compares the effective eccentricity of intersections of reconstructed volumes and imaging planes. As another metric of volumetric reconstruction fidelity, the volume of one of the contrasting spheres was determined using manual contouring. By comparing these manually delineated volumes with a CBCT reconstruction, we can gauge the volumetric fidelity of reconstructions. The configuration which yielded the highest overall volumetric reconstruction fidelity, as determined by effective eccentricities and volumetric contouring, consisted of two orthogonally-offset 60° L-arm sweeps and a single C-arm sweep which shared a pivot point with one the L-arm sweeps. When compared to a similar configuration that lacked the C-arm component, it is shown that the C-arm improves the delineation of volumes along the transverse axis. The results described herein suggest that volumetric reconstruction using multiple, unconstrained orthogonal sweeps can provide an improvement compared with traditional cone beam CT using standard axial rotations.

Incorporation of a Laguerre-Gauss channelized Hotelling observer for false-positive reduction in a mammographic mass CAD system.

Previously, we developed a simple Laguerre-Gauss (LG) channelized Hotelling observer (CHO) for incorporation into our mass computer-aided detection (CAD) system. This LG-CHO was trained using initial detection suspicious region data and was empirically optimized for free parameters. For the study presented in this paper, we wish to create a more optimal mass detection observer based on a novel combination of LG channels. A large set of LG channels with differing free parameters was created. Each of these channels was applied to the suspicious regions, and an output test statistic was determined. A stepwise feature selection algorithm was used to determine which LG channels would combine best to detect masses. These channels were combined using a HO to create a single template for the mass CAD system. Results from free-response receiver operating characteristic curves demonstrated that the incorporation of the novel LG-CHO into the CAD system slightly improved performance in high-sensitivity regions.

Characterization of difference of Gaussian filters in the detection of mammographic regions.

In this article, we present a characterization of the effect of difference of Gaussians (DoG) filters in the detection of mammographic regions. DoG filters have been used previously in mammographic mass computer-aided detection (CAD) systems. As DoG filters are constructed from the subtraction of two bivariate Gaussian distributions, they require the specification of three parameters: the size of the filter template and the standard deviations of the constituent Gaussians. The influence of these three parameters in the detection of mammographic masses has not been characterized. In this work, we aim to determine how the parameters affect (1) the physical descriptors of the detected regions, (2) the true and false positive rates, and (3) the classification performance of the individual descriptors. To this end, 30 DoG filters are created from the combination of three template sizes and four values for each of the Gaussians' standard deviations. The filters are used to detect regions in a study database of 181 craniocaudal-view mammograms extracted from the Digital Database for Screening Mammography. To describe the physical characteristics of the identified regions, morphological and textural features are extracted from each of the detected regions. Differences in the mean values of the features caused by altering the DoG parameters are examined through statistical and empirical comparisons. The parameters' effects on the true and false positive rate are determined by examining the mean malignant sensitivities and false positives per image (FPpI). Finally, the effect on the classification performance is described by examining the variation in FPpI at the point where 81% of the malignant masses in the study database are detected. Overall, the findings of the study indicate that increasing the standard deviations of the Gaussians used to construct a DoG filter results in a dramatic decrease in the number of regions identified at the expense of missing a small number of malignancies. The sharp reduction in the number of identified regions allowed the identification of textural differences between large and small mammographic regions. We find that the classification performances of the features that achieve the lowest average FPpI are influenced by all three of the parameters.

Predicting radiotherapy-induced cardiac perfusion defects.

The purpose of this work is to compare the efficacy of mathematical models in predicting the occurrence of radiotherapy-induced left ventricular perfusion defects assessed using single-photon emission computed tomography (SPECT). The basis of this study is data from 73 left-sided breast/ chestwall patients treated with tangential photon fields. The mathematical models compared were three commonly used parametric models [Lyman normal tissue complication probability (LNTCP), relative serialty (RS), generalized equivalent uniform dose (gEUD)] and a nonparametric model (Linear discriminant analysis--LDA). Data used by the models were the left ventricular dose--volume histograms, or SPECT-based dose-function histograms, and the presence/absence of SPECT perfusion defects 6 months postradiation therapy (21 patients developed defects). For the parametric models, maximum likelihood estimation and F-tests were used to fit the model parameters. The nonparametric LDA model step-wise selected features (volumes/function above dose levels) using a method based on receiver operating characteristics (ROC) analysis to best separate the groups with and without defects. Optimistic (upper bound) and pessimistic (lower bound) estimates of each model's predictive capability were generated using ROC curves. A higher area under the ROC curve indicates a more accurate model (a model that is always accurate has area = 1). The areas under these curves for different models were used to statistically test for differences between them. Pessimistic estimates of areas under the ROC curve using dose-volume histogram/ dose-function histogram inputs, in order of increasing prediction accuracy, were LNTCP (0.79/0.75), RS (0.80/0.77), gEUD (0.81/0.78), and LDA (0.84/0.86). Only the LDA model benefited from SPECT-based regional functional information. In general, the LDA model was statistically superior to the parametric models. The LDA model selected as features the left ventricular volumes above approximately 23 Gy (V23), essentially volume in field, and 33 Gy (V33), as best separating the groups with and without defects. In conclusion, the nonparametric LDA model appears to be a more accurate predictor of radiotherapy-induced left ventricular perfusion defects than commonly used parametric models.

Penalized likelihood fluence optimization with evolutionary components for intensity modulated radiation therapy treatment planning.

A novel iterative penalized likelihood algorithm with evolutionary components for the optimization of beamlet fluences for intensity modulated radiation therapy (IMRT) is presented. This algorithm is designed to be flexible in terms of the objective function and automatically escalates dose, as long as the objective function increases and all constraints are met. For this study, the objective function employed was the product of target equivalent uniform dose (EUD) and fraction of target tissue within set homogeneity constraints. The likelihood component of the algorithm iteratively attempts to minimize the mean squared error between a homogeneous dose prescription and the actual target dose distribution. The updated beamlet fluences are then adjusted via a quadratic penalty function that is based on the dose-volume histogram (DVH) constraints of the organs at risk. The evolutionary components were included to prevent the algorithm from converging to a local maximum. The algorithm was applied to a prostate cancer dataset, with especially difficult DVH constraints on bladder, rectum, and femoral heads. Dose distributions were generated for manually selected sets of three-, four-, five-, and seven-field treatment plans. Additionally, a global search was performed to find the optimal orientations for an axial three-beam plan. The results from this optimal orientation set were compared to results for manually selected orientation (gantry angle) sets of 3- (0 degrees, 90 degrees, 270 degrees), 4- (0 degrees, 90 degrees, 180 degrees, 270 degrees), 5- (0 degrees, 50 degrees, 130 degrees, 230 degrees, 310 degrees), and 7- (0 degrees, 40 degrees, 90 degrees, 140 degrees, 230 degrees, 270 degrees, 320 degrees) field axial treatment plans. For all the plans generated, all DVH constraints were met and average optimization computation time was approximately 30 seconds. For the manually selected orientations, the algorithm was successful in providing a relatively homogeneous target dose distribution, while simultaneously satisfying dose-volume limits by diverting dose away from proximal critical structures. The global search for an optimal three-beam orientation set yielded gantry angles of 70 degrees, 170 degrees, and 320 degrees. The EUD for this orientation set was 58 Gy, with 96% of the target within the set upper and lower limits. In comparison, optimized EUDs for the manually selected orientation sets of three, four, five and seven beams were 52.3, 52.6, 56.9, and 61.3 Gy, respectively. The orientation optimized three-beam plan yielded higher EUDs than the manually selected three-, four-, and five-beam plans, but lower EUDs than the seven-beam plan. In conclusion, a novel penalized likelihood algorithm with evolutionary components has successfully been implemented to optimize beamlet fluences for IMRT. Initial results are promising for dose conformity and uniformity of dose to target. When combined with optimal beam orientation selection for prostate cancer treatment planning, the results indicate that plans with a small number of optimized beam orientations achieve results comparable to those with a larger number of conventionally oriented beams.

Incorporation of an iterative, linear segmentation routine into a mammographic mass CAD system.

In previous research, we have developed a computer-aided detection (CAD) system designed to detect masses in mammograms. The previous version of our system employed a simple but imprecise method to localize the masses. In this research, we present a more robust segmentation routine for use with mammographic masses. Our hypothesis is that by more accurately describing the morphology of the masses, we can improve the CAD system's ability to distinguish masses from other mammographic structures. To test this hypothesis, we incorporated the new segmentation routine into our CAD system and examined the change in performance. The developed iterative, linear segmentation routine is a gray level-based procedure. Using the identified regions from the previous CAD system as the initial seeds, the new segmentation algorithm refines the suspicious mass borders by making estimates of the interior and exterior pixels. These estimates are then passed to a linear discriminant, which determines the optimal threshold between the interior and exterior pixels. After applying the threshold and identifying the object's outline, two constraints on the border are applied to reduce the influence of background noise. After the border is constrained, the process repeats until a stopping criterion is reached. The segmentation routine was tested on a study database of 183 mammographic images extracted from the Digital Database for Screening Mammography. Eighty-three of the images contained 50 malignant and 50 benign masses; 100 images contained no masses. The previously developed CAD system was used to locate a set of suspicious regions of interest (ROIs) within the images. To assess the performance of the segmentation algorithm, a set of 20 features was measured from the suspicious regions before and after the application of the developed segmentation routine. Receiver operating characteristic (ROC) analysis was employed on the ROIs to examine the discriminatory capabilities of each individual feature before and after the segmentation routine. A statistically significant performance increase was found in many of the individual features, particularly those describing the mass borders. To examine how the incorporation of the segmentation routine affected the performance of the overall CAD system, free-response ROC (FROC) analysis was employed. When considering only malignant masses, the FROC performance of the system with the segmentation routine appeared better than the previous system. When detecting 90% of the malignant masses, the previous system achieved 4.9 false positives per image (FPpI) compared to the post-segmentation system's 4.2 FPpI. At 80% sensitivity, the respective FPpI were 3.5 and 1.6.

Computer aided detection of masses in mammography using subregion Hotelling observers.

We propose to investigate the use of the subregion Hotelling observer for the basis of a computer aided detection scheme for masses in mammography. A database of 1320 regions of interest (ROIs) was selected from the DDSM database collected by the University of South Florida using the Lumisys scanner cases. The breakdown of the cases was as follows: 656 normal ROIs, 307 benign ROIs, and 357 cancer ROIs. Each ROI was extracted at a size of 1024 x 1024 pixels and sub-sampled to 128 x 128 pixels. For the detection task, cancer and benign cases were considered positive and normal was considered negative. All positive cases had the lesion centered in the ROI. We chose to investigate the subregion Hotelling observer as a classifier to detect masses. The Hotelling observer incorporates information about the signal, the background, and the noise correlation for prediction of positive and negative and is the optimal detector when these are known. For our study, 225 subregion Hotelling observers were set up in a 15 x 15 grid across the center of the ROIs. Each separate observer was designed to "observe," or discriminate, an 8 x 8 pixel area of the image. A leave one out training and testing methodology was used to generate 225 "features," where each feature is the output of the individual observers. The 225 features derived from separate Hotelling observers were then narrowed down by using forward searching linear discriminants (LDs). The reduced set of features was then analyzed using an additional LD with receiver operating characteristic (ROC) analysis. The 225 Hotelling observer features were searched by the forward searching LD, which selected a subset of 37 features. This subset of 37 features was then analyzed using an additional LD, which gave a ROC area under the curve of 0.9412 +/- 0.006 and a partial area of 0.6728. Additionally, at 98% sensitivity the overall classifier had a specificity of 55.9% and a positive predictive value of 69.3%. Preliminary results suggest that using subregion Hotelling observers in combination with LDs can provide a strong backbone for a CAD scheme to help radiologists with detection. Such a system could be used in conjunction with CAD systems for false positive reduction.

Chest radiography: optimization of X-ray spectrum for cesium iodide-amorphous silicon flat-panel detector.

PURPOSE: To ascertain the optimum x-ray spectrum for chest radiography with a cesium iodide-amorphous silicon flat-panel detector. MATERIALS AND METHODS: End points for optimization included the ratio of tissue contrast to bone contrast and a figure of merit (FOM) equal to the square of the signal-to-noise ratio of tissue divided by incident exposure to the patient. Studies were conducted with both computer spectrum modeling and experimental measurement in narrow-beam and full-field exposure conditions for four tissue thicknesses (8-32 cm). Three parameters that affect spectra were considered: the atomic number (Z) of filter material (Z = 13, 26, 29, 42, 50, 56, 64, 74, and 82), kilovoltage (from 50 to 150 kVp), and filter thickness (from 0.25 to 2.00 half-value layer [HVL]). RESULTS: Computer modeling and narrow-beam experimental data showed similar trends for the full range of parameters evaluated. Spectrum model results showed that copper filtration at 120 kVp or more was optimum for FOM. The ratio of contrasts showed a trend to be higher with higher kilovoltage and only a minor variation with filter material. Full-field experimental results, which reflect the added contribution of x-ray scatter, differed in magnitude but not trends from the narrow-beam data in all cases except the ratio of contrasts in the mediastinum. CONCLUSION: The best performance overall, including both FOM and ratio of contrasts, was at 120 kVp with 1-HVL copper filtration (0.2 mm). With this beam spectrum and an increase in tube output (ie, milliampere seconds) of about 50%, a chest radiograph can be obtained with image quality approximately equal to that with a conventional spectrum but with about 25% less patient exposure.