Search NASA⌕ Search

SEARCH · Search NASA

Results for “Histogram”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

Tomotherapy dose distribution verification using MAGIC-f polymer gel dosimetry

Purpose: This paper presents the application of MAGIC-f gel in a three-dimensional dose distribution measurement and its ability to accurately measure the dose distribution from a tomotherapy unit. Methods: A prostate intensity-modulated radiation therapy (IMRT) irradiation was simulated in the gel phantom and the treatment was delivered by a TomoTherapy equipment. Dose distribution was evaluated by the R2 distribution measured in magnetic resonance imaging. Results: A high similarity was found by overlapping of isodoses of the dose distribution measured with the gel and expected by the treatment planning system (TPS). Another analysis was done by comparing the relative absorbed dose profiles in the measured and in the expected dose distributions extracted along indicated lines of the volume and the results were also in agreement. The gamma index analysis was also applied to the data and a high pass rate was achieved (88.4% for analysis using 3%/3 mm and of 96.5% using 4%/4 mm). The real three-dimensional analysis compared the dose-volume histograms measured for the planning volumes and expected by the treatment planning, being the results also in good agreement by the overlapping of the curves. Conclusions: These results show that MAGIC-f gel is a promise for tridimensional dose distribution measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗

Evolving treatment plan quality criteria from institution-specific experience

Purpose: The dosimetric aspects of radiation therapy treatment plan quality are usually evaluated and reported with dose volume histogram (DVH) endpoints. For clinical practicality, a small number of representative quantities derived from the DVH are often used as dose endpoints to summarize the plan quality. National guidelines on reference values for such quantities for some standard treatment approaches are often used as acceptance criteria to trigger treatment plan review. On the other hand, treatment prescription and planning approaches specific to each institution warrants the need to report plan quality in terms of practice consistency and with respect to institution-specific experience. The purpose of this study is to investigate and develop a systematic approach to record and characterize the institution-specific plan experience and use such information to guide the design of plan quality criteria. In the clinical setting, this approach will assist in (1) improving overall plan quality and consistency and (2) detecting abnormal plan behavior for retrospective analysis. Methods: The authors propose a self-evolving methodology and have developed an in-house prototype software suite that (1) extracts the dose endpoints from a treatment plan and evaluates them against both national standard and institution-specific criteria and (2) evolves the statistics for the dose endpoints and updates institution-specific criteria. Results: The validity of the proposed methodology was demonstrated with a database of prostate stereotactic body radiotherapy cases. As more data sets are accumulated, the evolving institution-specific criteria can serve as a reliable and stable consistency measure for plan quality and reveals the potential use of the ''tighter'' criteria than national standards or projected criteria, leading to practice that may push to shrink the gap between plans deemed acceptable and the underlying unknown optimality. Conclusions: The authors have developed a rationale to improve plan quality and consistency, by evolving the plan quality criteria from institution-specific experience, complementary to national standards. The validity of the proposed method was demonstrated with a prototype system on prostate stereotactic body radiotherapy (SBRT) cases. The current study uses direct and indirect DVH endpoints for plan quality evaluation, but the infrastructure proposed here applies to general outcome data as well. The authors expect forward evaluation together with intelligent update based on evidence-based learning, which will evolve the clinical practice for improved efficiency, consistency, and ultimately better treatment outcome.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dosimetric and deformation effects of image-guided interventions during stereotactic body radiation therapy of the prostate using an endorectal balloon

Purpose: During stereotactic body radiation therapy (SBRT) for the treatment of prostate cancer, an inflatable endorectal balloon (ERB) may be used to reduce motion of the target and reduce the dose to the posterior rectal wall. This work assessed the dosimetric impact of manual interventions on ERB position in patients receiving prostate SBRT and investigated the impact of ERB interventions on prostate shape. Methods: The data of seven consecutive patients receiving SBRT for the treatment of clinical stage T1cN0M0 prostate cancer enrolled in a multi-institutional, IRB-approved trial were analyzed. The SBRT dose was 50 Gy in five fractions to a planning target volume (PTV) that included the prostate (implanted with three fiducial markers) with a 3-5 mm margin. All plans were based on simulation images that included an ERB inflated with 60 cm{sup 3} of air. Daily kilovoltage cone-beam computed tomography (CBCT) imaging was performed to localize the PTV, and an automated fusion with the planning images yielded displacements required for PTV relocalization. When the ERB volume and/or position were judged to yield inaccurate repositioning, manual adjustment (ERB reinflation and/or repositioning) was performed. Based on all 59 CBCT image sets acquired, a deformable registration algorithm was used to determine the dose received by, displacement of, and deformation of the prostate, bladder (BLA), and anterior rectal wall (ARW). This dose tracking methodology was applied to images taken before and after manual adjustment of the ERB (intervention), and the delivered dose was compared to that which would have been delivered in the absence of intervention. Results: Interventions occurred in 24 out of 35 (69%) of the treated fractions. The direct effect of these interventions was an increase in the prostate radiation dose that included 95% of the PTV (D95) from 9.6 {+-} 1.0 to 10.0 {+-} 0.2 Gy (p = 0.06) and an increase in prostate coverage from 94.0% {+-} 8.5% to 97.8% {+-} 1.9% (p = 0.03). Additionally, ERB interventions reduced prostate deformation in the anterior-posterior (AP) direction, reduced errors in the sagittal rotation of the prostate, and increased the similarity in shape of the prostate to the radiotherapy plan (increased Dice coefficient from 0.76 {+-} 0.06 to 0.80 {+-} 0.04, p = 0.01). Postintervention decreases in prostate volume receiving less than the prescribed dose and decreases in the voxel-wise displacement of the prostate, bladder, and anterior rectal wall were observed, which resulted in improved dose-volume histogram (DVH) characteristics. Conclusions: Image-guided interventions in ERB volume and/or position during prostate SBRT were necessary to ensure the delivery of the dose distribution as planned. ERB interventions resulted in reductions in prostate deformations that would have prevented accurate localization of patient anatomy.

60 APPLIED LIFE SCIENCES↗

Quantifying the margin sharpness of lesions on radiological images for content-based image retrieval

Purpose: To develop a method to quantify the margin sharpness of lesions on CT and to evaluate it in simulations and CT scans of liver and lung lesions. Methods: The authors computed two attributes of margin sharpness: the intensity difference between a lesion and its surroundings, and the sharpness of the intensity transition across the lesion boundary. These two attributes were extracted from sigmoid curves fitted along lines automatically drawn orthogonal to the lesion margin. The authors then represented the margin characteristics for each lesion by a feature vector containing histograms of these parameters. The authors created 100 simulated CT scans of lesions over a range of intensity difference and margin sharpness, and used the concordance correlation between the known parameter and the corresponding computed feature as a measure of performance. The authors also evaluated their method in 79 liver lesions (44 patients: 23 M, 21 F, mean age 61) and 58 lung nodules (57 patients: 24 M, 33 F, mean age 66). The methodology presented takes into consideration the boundary of the liver and lung during feature extraction in clinical images to ensure that the margin feature do not get contaminated by anatomy other than the normal organ surrounding the lesions. For evaluation in these clinical images, the authors created subjective independent reference standards for pairwise margin sharpness similarity in the liver and lung cohorts, and compared rank orderings of similarity used using our sharpness feature to that expected from the reference standards using mean normalized discounted cumulative gain (NDCG) over all query images. In addition, the authors compared their proposed feature with two existing techniques for lesion margin characterization using the simulated and clinical datasets. The authors also evaluated the robustness of their features against variations in delineation of the lesion margin by simulating five types of deformations of the lesion margin. Equivalence across deformations was assessed using Schuirmann's paired two one-sided tests. Results: In simulated images, the concordance correlation between measured gradient and actual gradient was 0.994. The mean (s.d.) and standard deviation NDCG score for the retrieval of K images, K = 5, 10, and 15, were 84% (8%), 85% (7%), and 85% (7%) for CT images containing liver lesions, and 82% (7%), 84% (6%), and 85% (4%) for CT images containing lung nodules, respectively. The authors' proposed method outperformed the two existing margin characterization methods in average NDCG scores over all K, by 1.5% and 3% in datasets containing liver lesion, and 4.5% and 5% in datasets containing lung nodules. Equivalence testing showed that the authors' feature is more robust across all margin deformations (p < 0.05) than the two existing methods for margin sharpness characterization in both simulated and clinical datasets. Conclusions: The authors have described a new image feature to quantify the margin sharpness of lesions. It has strong correlation with known margin sharpness in simulated images and in clinical CT images containing liver lesions and lung nodules. This image feature has excellent performance for retrieving images with similar margin characteristics, suggesting potential utility, in conjunction with other lesion features, for content-based image retrieval applications.

60 APPLIED LIFE SCIENCES↗

Reduction of prostate intrafraction motion using gas-release rectal balloons

Purpose: To analyze prostate intrafraction motion using both non-gas-release (NGR) and gas-release (GR) rectal balloons and to evaluate the ability of GR rectal balloons to reduce prostate intrafraction motion. Methods: Twenty-nine patients with NGR rectal balloons and 29 patients with GR balloons were randomly selected from prostate patients treated with proton therapy at University of Florida Proton Therapy Institute (Jacksonville, FL). Their pretreatment and post-treatment orthogonal radiographs were analyzed, and both pretreatment setup residual error and intrafraction-motion data were obtained. Population histograms of intrafraction motion were plotted for both types of balloons. Population planning target-volume (PTV) margins were calculated with the van Herk formula of 2.5{Sigma}+ 0.7{sigma} to account for setup residual errors and intrafraction motion errors. Results: Pretreatment and post-treatment radiographs indicated that the use of gas-release rectal balloons reduced prostate intrafraction motion along superior-inferior (SI) and anterior-posterior (AP) directions. Similar patient setup residual errors were exhibited for both types of balloons. Gas-release rectal balloons resulted in PTV margin reductions from 3.9 to 2.8 mm in the SI direction, 3.1 to 1.8 mm in the AP direction, and an increase from 1.9 to 2.1 mm in the left-right direction. Conclusions: Prostate intrafraction motion is an important uncertainty source in radiotherapy after image-guided patient setup with online corrections. Compared to non-gas-release rectal balloons, gas-release balloons can reduce prostate intrafraction motion in the SI and AP directions caused by gas buildup.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Three-dimensional radiobiological dosimetry of kidneys for treatment planning in peptide receptor radionuclide therapy

Purpose: Peptide receptor radionuclide therapy (PRRT) delivers high absorbed doses to kidneys and may lead to permanent nephropathy. Reliable dosimetry of kidneys is thus critical for safe and effective PRRT. The aim of this work was to assess the feasibility of planning PRRT based on 3D radiobiological dosimetry (3D-RD) in order to optimize both the amount of activity to administer and the fractionation scheme, while limiting the absorbed dose and the biological effective dose (BED) to the renal cortex. Methods: Planar and SPECT data were available for a patient examined with {sup 111}In-DTPA-octreotide at 0.5 (planar only), 4, 24, and 48 h post-injection. Absorbed dose and BED distributions were calculated for common therapeutic radionuclides, i.e., {sup 111}In, {sup 90}Y and {sup 177}Lu, using the 3D-RD methodology. Dose-volume histograms were computed and mean absorbed doses to kidneys, renal cortices, and medullae were compared with results obtained using the MIRD schema (S-values) with the multiregion kidney dosimetry model. Two different treatment planning approaches based on (1) the fixed absorbed dose to the cortex and (2) the fixed BED to the cortex were then considered to optimize the activity to administer by varying the number of fractions. Results: Mean absorbed doses calculated with 3D-RD were in good agreement with those obtained with S-value-based SPECT dosimetry for {sup 90}Y and {sup 177}Lu. Nevertheless, for {sup 111}In, differences of 14% and 22% were found for the whole kidneys and the cortex, respectively. Moreover, the authors found that planar-based dosimetry systematically underestimates the absorbed dose in comparison with SPECT-based methods, up to 32%. Regarding the 3D-RD-based treatment planning using a fixed BED constraint to the renal cortex, the optimal number of fractions was found to be 3 or 4, depending on the radionuclide administered and the value of the fixed BED. Cumulative activities obtained using the proposed simulated treatment planning are compatible with real activities administered to patients in PRRT. Conclusions: The 3D-RD treatment planning approach based on the fixed BED was found to be the method of choice for clinical implementation in PRRT by providing realistic activity to administer and number of cycles. While dividing the activity in several cycles is important to reduce renal toxicity, the clinical outcome of fractionated PRRT should be investigated in the future.

60 APPLIED LIFE SCIENCES↗

A digitally reconstructed radiograph algorithm calculated from first principles

Purpose: To develop an algorithm for computing realistic digitally reconstructed radiographs (DRRs) that match real cone-beam CT (CBCT) projections with no artificial adjustments. Methods: The authors used measured attenuation data from cone-beam CT projection radiographs of different materials to obtain a function to convert CT number to linear attenuation coefficient (LAC). The effects of scatter, beam hardening, and veiling glare were first removed from the attenuation data. Using this conversion function the authors calculated the line integral of LAC through a CT along rays connecting the radiation source and detector pixels with a ray-tracing algorithm, producing raw DRRs. The effects of scatter, beam hardening, and veiling glare were then included in the DRRs through postprocessing. Results: The authors compared actual CBCT projections to DRRs produced with all corrections (scatter, beam hardening, and veiling glare) and to uncorrected DRRs. Algorithm accuracy was assessed through visual comparison of projections and DRRs, pixel intensity comparisons, intensity histogram comparisons, and correlation plots of DRR-to-projection pixel intensities. In general, the fully corrected algorithm provided a small but nontrivial improvement in accuracy over the uncorrected algorithm. The authors also investigated both measurement- and computation-based methods for determining the beam hardening correction, and found the computation-based method to be superior, as it accounted for nonuniform bowtie filter thickness. The authors benchmarked the algorithm for speed and found that it produced DRRs in about 0.35 s for full detector and CT resolution at a ray step-size of 0.5 mm. Conclusions: The authors have demonstrated a DRR algorithm calculated from first principles that accounts for scatter, beam hardening, and veiling glare in order to produce accurate DRRs. The algorithm is computationally efficient, making it a good candidate for iterative CT reconstruction techniques that require a data fidelity term based on the matching of DRRs and projections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dosimetric accuracy of a deterministic radiation transport based {sup 192}Ir brachytherapy treatment planning system. Part III. Comparison to Monte Carlo simulation in voxelized anatomical computational models

Purpose: To compare TG43-based and Acuros deterministic radiation transport-based calculations of the BrachyVision treatment planning system (TPS) with corresponding Monte Carlo (MC) simulation results in heterogeneous patient geometries, in order to validate Acuros and quantify the accuracy improvement it marks relative to TG43. Methods: Dosimetric comparisons in the form of isodose lines, percentage dose difference maps, and dose volume histogram results were performed for two voxelized mathematical models resembling an esophageal and a breast brachytherapy patient, as well as an actual breast brachytherapy patient model. The mathematical models were converted to digital imaging and communications in medicine (DICOM) image series for input to the TPS. The MCNP5 v.1.40 general-purpose simulation code input files for each model were prepared using information derived from the corresponding DICOM RT exports from the TPS. Results: Comparisons of MC and TG43 results in all models showed significant differences, as reported previously in the literature and expected from the inability of the TG43 based algorithm to account for heterogeneities and model specific scatter conditions. A close agreement was observed between MC and Acuros results in all models except for a limited number of points that lay in the penumbra of perfectly shaped structures in the esophageal model, or at distances very close to the catheters in all models. Conclusions: Acuros marks a significant dosimetry improvement relative to TG43. The assessment of the clinical significance of this accuracy improvement requires further work. Mathematical patient equivalent models and models prepared from actual patient CT series are useful complementary tools in the methodology outlined in this series of works for the benchmarking of any advanced dose calculation algorithm beyond TG43.

60 APPLIED LIFE SCIENCES↗

PTV-based IMPT optimization incorporating planning risk volumes vs robust optimization

Purpose: Robust optimization leads to intensity-modulated proton therapy (IMPT) plans that are less sensitive to uncertainties and superior in terms of organs-at-risk (OARs) sparing, target dose coverage, and homogeneity compared to planning target volume (PTV)-based optimized plans. Robust optimization incorporates setup and range uncertainties, which implicitly adds margins to both targets and OARs and is also able to compensate for perturbations in dose distributions within targets and OARs caused by uncertainties. In contrast, the traditional PTV-based optimization considers only setup uncertainties and adds a margin only to targets but no margins to the OARs. It also ignores range uncertainty. The purpose of this work is to determine if robustly optimized plans are superior to PTV-based plans simply because the latter do not assign margins to OARs during optimization. Methods: The authors retrospectively selected from their institutional database five patients with head and neck (H and N) cancer and one with prostate cancer for this analysis. Using their original images and prescriptions, the authors created new IMPT plans using three methods: PTV-based optimization, optimization based on the PTV and planning risk volumes (PRVs) (i.e., 'PTV+PRV-based optimization'), and robust optimization using the 'worst-case' dose distribution. The PRVs were generated by uniformly expanding OARs by 3 mm for the H and N cases and 5 mm for the prostate case. The dose-volume histograms (DVHs) from the worst-case dose distributions were used to assess and compare plan quality. Families of DVHs for each uncertainty for all structures of interest were plotted along with the nominal DVHs. The width of the 'bands' of DVHs was used to quantify the plan sensitivity to uncertainty. Results: Compared with conventional PTV-based and PTV+PRV-based planning, robust optimization led to a smaller bandwidth for the targets in the face of uncertainties {l_brace}clinical target volume [CTV] bandwidth: 0.59 [robust], 3.53 [PTV+PRV], and 3.53 [PTV] Gy (RBE){r_brace}. It also resulted in higher doses to 95% of the CTV {l_brace}D{sub 95%}: 60.8 [robust] vs 59.3 [PTV+PRV] vs 59.6 [PTV] Gy (RBE){r_brace}, smaller D{sub 5%} (doses to 5% of the CTV) minus D{sub 95%} {l_brace}D{sub 5%}- D{sub 95%}: 13.2 [robust] vs 17.5 [PTV+PRV] vs 17.2 [PTV] Gy (RBE){r_brace}. At the same time, the robust optimization method irradiated OARs less {l_brace}maximum dose to 1 cm{sup 3} of the brainstem: 48.3 [robust] vs 48.8 [PTV+PRV] vs 51.2 [PTV] Gy (RBE); mean dose to the oral cavity: 22.3 [robust] vs 22.9 [PTV+PRV] vs 26.1 [PTV] Gy (RBE); maximum dose to 1% of the normal brain: 66.0 [robust] vs 68.0 [PTV+PRV] vs 69.3 [PTV] Gy (RBE){r_brace}. Conclusions: For H and N cases studied, OAR sparing in PTV+PRV-based optimization was inferior compared to robust optimization but was superior compared to PTV-based optimization; however, target dose robustness and homogeneity were comparable in the PTV+PRV-based and PTV-based optimizations. The same pattern held for the prostate case. The authors' data suggest that the superiority of robust optimization is not due simply to its inclusion of margins for OARs, but that this is due mainly to the ability of robust optimization to compensate for perturbations in dose distributions within target volumes and normal tissues caused by uncertainties.

60 APPLIED LIFE SCIENCES↗

Objective models of compressed breast shapes undergoing mammography

Purpose: To develop models of compressed breasts undergoing mammography based on objective analysis, that are capable of accurately representing breast shapes in acquired clinical images and generating new, clinically realistic shapes. Methods: An automated edge detection algorithm was used to catalogue the breast shapes of clinically acquired cranio-caudal (CC) and medio-lateral oblique (MLO) view mammograms from a large database of digital mammography images. Principal component analysis (PCA) was performed on these shapes to reduce the information contained within the shapes to a small number of linearly independent variables. The breast shape models, one of each view, were developed from the identified principal components, and their ability to reproduce the shape of breasts from an independent set of mammograms not used in the PCA, was assessed both visually and quantitatively by calculating the average distance error (ADE). Results: The PCA breast shape models of the CC and MLO mammographic views based on six principal components, in which 99.2% and 98.0%, respectively, of the total variance of the dataset is contained, were found to be able to reproduce breast shapes with strong fidelity (CC view mean ADE = 0.90 mm, MLO view mean ADE = 1.43 mm) and to generate new clinically realistic shapes. The PCA models based on fewer principal components were also successful, but to a lesser degree, as the two-component model exhibited a mean ADE = 2.99 mm for the CC view, and a mean ADE = 4.63 mm for the MLO view. The four-component models exhibited a mean ADE = 1.47 mm for the CC view and a mean ADE = 2.14 mm for the MLO view. Paired t-tests of the ADE values of each image between models showed that these differences were statistically significant (max p-value = 0.0247). Visual examination of modeled breast shapes confirmed these results. Histograms of the PCA parameters associated with the six principal components were fitted with Gaussian distributions. The six-component model was also used to generate CC and MLO view mammogram breast shapes, using the mean PCA parameter values of these distributions and randomly generated values based on the fitted Gaussian distributions, which resemble clinically encountered breasts. A spreadsheet with the data necessary to apply this model is provided as the supplementary material. Conclusions: Our PCA models of breast shapes in both mammographic views successfully reproduce analyzed breast shapes and generate new clinically relevant shapes. This work can aid in research applications which incorporate breast shape modeling, such as x-ray scatter correction, dosimetry, and image registration.

60 APPLIED LIFE SCIENCES↗

A gEUD-based inverse planning technique for HDR prostate brachytherapy: Feasibility study

Purpose: The purpose of this work was to study the feasibility of a new inverse planning technique based on the generalized equivalent uniform dose for image-guided high dose rate (HDR) prostate cancer brachytherapy in comparison to conventional dose-volume based optimization. Methods: The quality of 12 clinical HDR brachytherapy implants for prostate utilizing HIPO (Hybrid Inverse Planning Optimization) is compared with alternative plans, which were produced through inverse planning using the generalized equivalent uniform dose (gEUD). All the common dose-volume indices for the prostate and the organs at risk were considered together with radiobiological measures. The clinical effectiveness of the different dose distributions was investigated by comparing dose volume histogram and gEUD evaluators. Results: Our results demonstrate the feasibility of gEUD-based inverse planning in HDR brachytherapy implants for prostate. A statistically significant decrease in D{sub 10} or/and final gEUD values for the organs at risk (urethra, bladder, and rectum) was found while improving dose homogeneity or dose conformity of the target volume. Conclusions: Following the promising results of gEUD-based optimization in intensity modulated radiation therapy treatment optimization, as reported in the literature, the implementation of a similar model in HDR brachytherapy treatment plan optimization is suggested by this study. The potential of improved sparing of organs at risk was shown for various gEUD-based optimization parameter protocols, which indicates the ability of this method to adapt to the user's preferences.

60 APPLIED LIFE SCIENCES↗

Needle segmentation using 3D Hough transform in 3D TRUS guided prostate transperineal therapy

Purpose: Prostate adenocarcinoma is the most common noncutaneous malignancy in American men with over 200 000 new cases diagnosed each year. Prostate interventional therapy, such as cryotherapy and brachytherapy, is an effective treatment for prostate cancer. Its success relies on the correct needle implant position. This paper proposes a robust and efficient needle segmentation method, which acts as an aid to localize the needle in three-dimensional (3D) transrectal ultrasound (TRUS) guided prostate therapy. Methods: The procedure of locating the needle in a 3D TRUS image is a three-step process. First, the original 3D ultrasound image containing a needle is cropped; the cropped image is then converted to a binary format based on its histogram. Second, a 3D Hough transform based needle segmentation method is applied to the 3D binary image in order to locate the needle axis. The position of the needle endpoint is finally determined by an optimal threshold based analysis of the intensity probability distribution. The overall efficiency is improved through implementing a coarse-fine searching strategy. The proposed method was validated in tissue-mimicking agar phantoms, chicken breast phantoms, and 3D TRUS patient images from prostate brachytherapy and cryotherapy procedures by comparison to the manual segmentation. The robustness of the proposed approach was tested by means of varying parameters such as needle insertion angle, needle insertion length, binarization threshold level, and cropping size. Results: The validation results indicate that the proposed Hough transform based method is accurate and robust, with an achieved endpoint localization accuracy of 0.5 mm for agar phantom images, 0.7 mm for chicken breast phantom images, and 1 mm for in vivo patient cryotherapy and brachytherapy images. The mean execution time of needle segmentation algorithm was 2 s for a 3D TRUS image with size of 264 Multiplication-Sign 376 Multiplication-Sign 630 voxels. Conclusions: The proposed needle segmentation algorithm is accurate, robust, and suitable for 3D TRUS guided prostate transperineal therapy.

60 APPLIED LIFE SCIENCES↗

A support vector machine classifier reduces interscanner variation in the HRCT classification of regional disease pattern in diffuse lung disease: Comparison to a Bayesian classifier

Purpose: To investigate the effect of using different computed tomography (CT) scanners on the accuracy of high-resolution CT (HRCT) images in classifying regional disease patterns in patients with diffuse lung disease, support vector machine (SVM) and Bayesian classifiers were applied to multicenter data. Methods: Two experienced radiologists marked sets of 600 rectangular 20 Multiplication-Sign 20 pixel regions of interest (ROIs) on HRCT images obtained from two scanners (GE and Siemens), including 100 ROIs for each of local patterns of lungs-normal lung and five of regional pulmonary disease patterns (ground-glass opacity, reticular opacity, honeycombing, emphysema, and consolidation). Each ROI was assessed using 22 quantitative features belonging to one of the following descriptors: histogram, gradient, run-length, gray level co-occurrence matrix, low-attenuation area cluster, and top-hat transform. For automatic classification, a Bayesian classifier and a SVM classifier were compared under three different conditions. First, classification accuracies were estimated using data from each scanner. Next, data from the GE and Siemens scanners were used for training and testing, respectively, and vice versa. Finally, all ROI data were integrated regardless of the scanner type and were then trained and tested together. All experiments were performed based on forward feature selection and fivefold cross-validation with 20 repetitions. Results: For each scanner, better classification accuracies were achieved with the SVM classifier than the Bayesian classifier (92% and 82%, respectively, for the GE scanner; and 92% and 86%, respectively, for the Siemens scanner). The classification accuracies were 82%/72% for training with GE data and testing with Siemens data, and 79%/72% for the reverse. The use of training and test data obtained from the HRCT images of different scanners lowered the classification accuracy compared to the use of HRCT images from the same scanner. For integrated ROI data obtained from both scanners, the classification accuracies with the SVM and Bayesian classifiers were 92% and 77%, respectively. The selected features resulting from the classification process differed by scanner, with more features included for the classification of the integrated HRCT data than for the classification of the HRCT data from each scanner. For the integrated data, consisting of HRCT images of both scanners, the classification accuracy based on the SVM was statistically similar to the accuracy of the data obtained from each scanner. However, the classification accuracy of the integrated data using the Bayesian classifier was significantly lower than the classification accuracy of the ROI data of each scanner. Conclusions: The use of an integrated dataset along with a SVM classifier rather than a Bayesian classifier has benefits in terms of the classification accuracy of HRCT images acquired with more than one scanner. This finding is of relevance in studies involving large number of images, as is the case in a multicenter trial with different scanners.

60 APPLIED LIFE SCIENCES↗

Dosimetric impact of Acuros XB deterministic radiation transport algorithm for heterogeneous dose calculation in lung cancer

Purpose: The novel deterministic radiation transport algorithm, Acuros XB (AXB), has shown great potential for accurate heterogeneous dose calculation. However, the clinical impact between AXB and other currently used algorithms still needs to be elucidated for translation between these algorithms. The purpose of this study was to investigate the impact of AXB for heterogeneous dose calculation in lung cancer for intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT). Methods: The thorax phantom from the Radiological Physics Center (RPC) was used for this study. IMRT and VMAT plans were created for the phantom in the Eclipse 11.0 treatment planning system. Each plan was delivered to the phantom three times using a Varian Clinac iX linear accelerator to ensure reproducibility. Thermoluminescent dosimeters (TLDs) and Gafchromic EBT2 film were placed inside the phantom to measure delivered doses. The measurements were compared with dose calculations from AXB 11.0.21 and the anisotropic analytical algorithm (AAA) 11.0.21. Two dose reporting modes of AXB, dose-to-medium in medium (D{sub m,m}) and dose-to-water in medium (D{sub w,m}), were studied. Point doses, dose profiles, and gamma analysis were used to quantify the agreement between measurements and calculations from both AXB and AAA. The computation times for AAA and AXB were also evaluated. Results: For the RPC lung phantom, AAA and AXB dose predictions were found in good agreement to TLD and film measurements for both IMRT and VMAT plans. TLD dose predictions were within 0.4%-4.4% to AXB doses (both D{sub m,m} and D{sub w,m}); and within 2.5%-6.4% to AAA doses, respectively. For the film comparisons, the gamma indexes ({+-}3%/3 mm criteria) were 94%, 97%, and 98% for AAA, AXB{sub Dm,m}, and AXB{sub Dw,m}, respectively. The differences between AXB and AAA in dose-volume histogram mean doses were within 2% in the planning target volume, lung, heart, and within 5% in the spinal cord. However, differences up to 8% between AXB and AAA were found at lung/soft tissue interface regions for individual IMRT fields. AAA was found to be 5-6 times faster than AXB for IMRT, while AXB was 4-5 times faster than AAA for VMAT plan. Conclusions: AXB is satisfactorily accurate for the dose calculation in lung cancer for both IMRT and VMAT plans. The differences between AXB and AAA are generally small except in heterogeneous interface regions. AXB D{sub w,m} and D{sub m,m} calculations are similar inside the soft tissue and lung regions. AXB can benefit lung VMAT plans by both improving accuracy and reducing computation time.

60 APPLIED LIFE SCIENCES↗

A study of IMRT planning parameters on planning efficiency, delivery efficiency, and plan quality

Purpose: To improve planning and delivery efficiency of head and neck IMRT without compromising planning quality through the evaluation of inverse planning parameters.Methods: Eleven head and neck patients with pre-existing IMRT treatment plans were selected for this retrospective study. The Pinnacle treatment planning system (TPS) was used to compute new treatment plans for each patient by varying the individual or the combined parameters of dose/fluence grid resolution, minimum MU per segment, and minimum segment area. Forty-five plans per patient were generated with the following variations: 4 dose/fluence grid resolution plans, 12 minimum segment area plans, 9 minimum MU plans, and 20 combined minimum segment area/minimum MU plans. Each plan was evaluated and compared to others based on dose volume histograms (DVHs) (i.e., plan quality), planning time, and delivery time. To evaluate delivery efficiency, a model was developed that estimated the delivery time of a treatment plan, and validated through measurements on an Elekta Synergy linear accelerator. Results: The uncertainty (i.e., variation) of the dose-volume index due to dose calculation grid variation was as high as 8.2% (5.5 Gy in absolute dose) for planning target volumes (PTVs) and 13.3% (2.1 Gy in absolute dose) for planning at risk volumes (PRVs). Comparison results of dose distributions indicated that smaller volumes were more susceptible to uncertainties. The grid resolution of a 4 mm dose grid with a 2 mm fluence grid was recommended, since it can reduce the final dose calculation time by 63% compared to the accepted standard (2 mm dose grid with a 2 mm fluence grid resolution) while maintaining a similar level of dose-volume index variation. Threshold values that maintained adequate plan quality (DVH results of the PTVs and PRVs remained satisfied for their dose objectives) were 5 cm{sup 2} for minimum segment area and 5 MU for minimum MU. As the minimum MU parameter was increased, the number of segments and delivery time were decreased. Increasing the minimum segment area parameter decreased the plan MU, but had less of an effect on the number of segments and delivery time. Our delivery time model predicted delivery time to within 1.8%. Conclusions: Increasing the dose grid while maintaining a small fluence grid allows for improved planning efficiency without compromising plan quality. Delivery efficiency can be improved by increasing the minimum MU, but not the minimum segment area. However, increasing the respective minimum MU and/or the minimum segment area to any value greater than 5 MU and 5 cm{sup 2} is not recommended because it degrades plan quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

3D inpatient dose reconstruction from the PET-CT imaging of {sup 90}Y microspheres for metastatic cancer to the liver: Feasibility study

Purpose: The introduction of radioembolization with microspheres represents a significant step forward in the treatment of patients with metastatic disease to the liver. This technique uses semiempirical formulae based on body surface area or liver and target volumes to calculate the required total activity for a given patient. However, this treatment modality lacks extremely important information, which is the three-dimensional (3D) dose delivered by microspheres to different organs after their administration. The absence of this information dramatically limits the clinical efficacy of this modality, specifically the predictive power of the treatment. Therefore, the aim of this study is to develop a 3D dose calculation technique that is based on the PET imaging of the infused microspheres.Methods: The Fluka Monte Carlo code was used to calculate the voxel dose kernel for {sup 90}Y source with voxel size equal to that of the PET scan. The measured PET activity distribution was converted to total activity distribution for the subsequent convolution with the voxel dose kernel to obtain the 3D dose distribution. In addition, dose-volume histograms were generated to analyze the dose to the tumor and critical structures.Results: The 3D inpatient dose distribution can be reconstructed from the PET data of a patient scanned after the infusion of microspheres. A total of seven patients have been analyzed so far using the proposed reconstruction method. Four patients underwent treatment with SIR-Spheres for liver metastases from colorectal cancer and three patients were treated with Therasphere for hepatocellular cancer. A total of 14 target tumors were contoured on post-treatment PET-CT scans for dosimetric evaluation. Mean prescription activity was 1.7 GBq (range: 0.58–3.8 GBq). The resulting mean maximum measured dose to targets was 167 Gy (range: 71–311 Gy). Mean minimum dose to 70% of target (D70) was 68 Gy (range: 25–155 Gy). Mean minimum dose to 90% of target (D90) was 53 Gy (range: 13–125 Gy).Conclusions: A three-dimensional inpatient dose reconstruction method has been developed that is based on the PET/CT data of a patient treated with {sup 90}Y microspheres. It allows for a complete description of the absorbed dose by the tumor and critical structures. It represents the first step in building predictive models for treatment outcomes for patients receiving this therapeutic modality as well as it allows for better analysis of patients' dose response and will ultimately improve future treatment administration.

61 RADIATION PROTECTION AND DOSIMETRY↗

Cobalt-60 tomotherapy: Clinical treatment planning and phantom dose delivery studies

Purpose: Investigations have shown that a Cobalt-60 (Co-60) radioactive source has the potential to play a role in intensity modulated radiation therapy (IMRT). In this paper, Co-60 tomotherapy's conformal dose delivery potential is evaluated by delivering conformal dose plans on a cylindrical homogeneous phantom containing clinical structures similar to those found in a typical head and neck (H and N) cancer. Also, the clinical potential of Co-60 tomotherapy is investigated by generating 2D clinical treatment plans for H and N and prostate anatomical regions. These plans are compared with the 6 MV based treatment plans for modalities such as linear accelerator-based tomotherapy and broad beam IMRT, and 15 MV based 3D conformal radiation therapy (3DCRT).Methods: For experimental validation studies, clinical and nonclinical conformal dose patterns were delivered on circular, homogeneous phantoms containing GafChromic film. For clinical planning study, dose calculations were performed with the EGSnrc Monte Carlo program, where a Theratronics 780C Co-60 unit and a 6 MV linear accelerator were modeled with a MIMiC binary multileaf collimator. An inhouse inverse treatment planning system was used to optimize tomotherapy plans using the same optimization parameters for both Co-60 and 6 MV beams. The IMRT and 3DCRT plans for the clinical cases were generated entirely in the Eclipse treatment planning system based on inhouse IMRT and 3DCRT site specific protocols.Results: The doses delivered to the homogeneous phantoms agreed with the calculations, indicating that it is possible to deliver highly conformal doses with the Co-60 unit. The dose distributions for Co-60 tomotherapy clinical plans for both clinical cases were similar to those obtained with 6 MV based tomotherapy and IMRT, and much more conformal compared to 3DCRT plans. The dose area histograms showed that the Co-60 plans achieve the dose objectives for the targets and organs at risk.Conclusions: These results confirm that Co-60 tomotherapy is capable of providing state-of-the-art conformal dose delivery and could be used for the treatment of targets in both small and larger separation anatomical regions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Low-frequency quantitative ultrasound imaging of cell death in vivo

Purpose: Currently, no clinical imaging modality is used routinely to assess tumor response to cancer therapies within hours to days of the delivery of treatment. Here, the authors demonstrate the efficacy of ultrasound at a clinically relevant frequency to quantitatively detect changes in tumors in response to cancer therapies using preclinical mouse models.Methods: Conventional low-frequency and corresponding high-frequency ultrasound (ranging from 4 to 28 MHz) were used along with quantitative spectroscopic and signal envelope statistical analyses on data obtained from xenograft tumors treated with chemotherapy, x-ray radiation, as well as a novel vascular targeting microbubble therapy.Results: Ultrasound-based spectroscopic biomarkers indicated significant changes in cell-death associated parameters in responsive tumors. Specifically changes in the midband fit, spectral slope, and 0-MHz intercept biomarkers were investigated for different types of treatment and demonstrated cell-death related changes. The midband fit and 0-MHz intercept biomarker derived from low-frequency data demonstrated increases ranging approximately from 0 to 6 dBr and 0 to 8 dBr, respectively, depending on treatments administrated. These data paralleled results observed for high-frequency ultrasound data. Statistical analysis of ultrasound signal envelope was performed as an alternative method to obtain histogram-based biomarkers and provided confirmatory results. Histological analysis of tumor specimens indicated up to 61% cell death present in the tumors depending on treatments administered, consistent with quantitative ultrasound findings indicating cell death. Ultrasound-based spectroscopic biomarkers demonstrated a good correlation with histological morphological findings indicative of cell death (r{sup 2}= 0.71, 0.82; p < 0.001).Conclusions: In summary, the results provide preclinical evidence, for the first time, that quantitative ultrasound used at a clinically relevant frequency, in addition to high-frequency ultrasound, can detect tissue changes associated with cell death in vivo in response to cancer treatments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗