Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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↗
DOE OSTI2013

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.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Strategies for automatic online treatment plan reoptimization using clinical treatment planning system: A planning parameters study

Purpose: Adaptive radiation therapy for prostate cancer using online reoptimization provides an improved control of interfractional anatomy variations. However, the clinical implementation of online reoptimization is currently limited by the low efficiency of current strategies and the difficulties associated with integration into the current treatment planning system. This study investigates the strategies for performing fast (∼2 min) automatic online reoptimization with a clinical fluence-map-based treatment planning system; and explores the performance with different input parameters settings: dose-volume histogram (DVH) objective settings, starting stage, and iteration number (in the context of real time planning).Methods: Simulated treatments of 10 patients were reoptimized daily for the first week of treatment (5 fractions) using 12 different combinations of optimization strategies. Options for objective settings included guideline-based RTOG objectives, patient-specific objectives based on anatomy on the planning CT, and daily-CBCT anatomy-based objectives adapted from planning CT objectives. Options for starting stages involved starting reoptimization with and without the original plan's fluence map. Options for iteration numbers were 50 and 100. The adapted plans were then analyzed by statistical modeling, and compared both in terms of dosimetry and delivery efficiency.Results: All online reoptimized plans were finished within ∼2 min with excellent coverage and conformity to the daily target. The three input parameters, i.e., DVH objectives, starting stage, and iteration number, contributed to the outcome of optimization nearly independently. Patient-specific objectives generally provided better OAR sparing compared to guideline-based objectives. The benefit in high-dose sparing from incorporating daily anatomy into objective settings was positively correlated with the relative change in OAR volumes from planning CT to daily CBCT. The use of the original plan fluence map as the starting stage reduced OAR dose at the mid-dose region, but increased the monitor units by 17%. Differences of only 2cc or less in OAR V50%/V70Gy/V76Gy were observed between 100 and 50 iterations.Conclusions: It is feasible to perform automatic online reoptimization in ∼2 min using a clinical treatment planning system. Selecting optimal sets of input parameters is the key to achieving high quality reoptimized plans, and should be based on the individual patient's daily anatomy, delivery efficiency, and time allowed for plan adaptation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Breast Imaging Reporting and Data System (BI-RADS) breast composition descriptors: Automated measurement development for full field digital mammography

Purpose: The Breast Imaging Reporting and Data System (BI-RADS) breast composition descriptors are used for standardized mammographic reporting and are assessed visually. This reporting is clinically relevant because breast composition can impact mammographic sensitivity and is a breast cancer risk factor. New techniques are presented and evaluated for generating automated BI-RADS breast composition descriptors using both raw and calibrated full field digital mammography (FFDM) image data.Methods: A matched case-control dataset with FFDM images was used to develop three automated measures for the BI-RADS breast composition descriptors. Histograms of each calibrated mammogram in the percent glandular (pg) representation were processed to create the new BR{sub pg} measure. Two previously validated measures of breast density derived from calibrated and raw mammograms were converted to the new BR{sub vc} and BR{sub vr} measures, respectively. These three measures were compared with the radiologist-reported BI-RADS compositions assessments from the patient records. The authors used two optimization strategies with differential evolution to create these measures: method-1 used breast cancer status; and method-2 matched the reported BI-RADS descriptors. Weighted kappa (κ) analysis was used to assess the agreement between the new measures and the reported measures. Each measure's association with breast cancer was evaluated with odds ratios (ORs) adjusted for body mass index, breast area, and menopausal status. ORs were estimated as per unit increase with 95% confidence intervals.Results: The three BI-RADS measures generated by method-1 had κ between 0.25–0.34. These measures were significantly associated with breast cancer status in the adjusted models: (a) OR = 1.87 (1.34, 2.59) for BR{sub pg}; (b) OR = 1.93 (1.36, 2.74) for BR{sub vc}; and (c) OR = 1.37 (1.05, 1.80) for BR{sub vr}. The measures generated by method-2 had κ between 0.42–0.45. Two of these measures were significantly associated with breast cancer status in the adjusted models: (a) OR = 1.95 (1.24, 3.09) for BR{sub pg}; (b) OR = 1.42 (0.87, 2.32) for BR{sub vc}; and (c) OR = 2.13 (1.22, 3.72) for BR{sub vr}. The radiologist-reported measures from the patient records showed a similar association, OR = 1.49 (0.99, 2.24), although only borderline statistically significant.Conclusions: A general framework was developed and validated for converting calibrated mammograms and continuous measures of breast density to fully automated approximations for the BI-RADS breast composition descriptors. The techniques are general and suitable for a broad range of clinical and research applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Toward adaptive radiotherapy for head and neck patients: Feasibility study on using CT-to-CBCT deformable registration for “dose of the day” calculations

Purpose: The aim of this study was to evaluate the appropriateness of using computed tomography (CT) to cone-beam CT (CBCT) deformable image registration (DIR) for the application of calculating the “dose of the day” received by a head and neck patient. Methods: NiftyReg is an open-source registration package implemented in our institution. The affine registration uses a Block Matching-based approach, while the deformable registration is a GPU implementation of the popular B-spline Free Form Deformation algorithm. Two independent tests were performed to assess the suitability of our registrations methodology for “dose of the day” calculations in a deformed CT. A geometric evaluation was performed to assess the ability of the DIR method to map identical structures between the CT and CBCT datasets. Features delineated in the planning CT were warped and compared with features manually drawn on the CBCT. The authors computed the dice similarity coefficient (DSC), distance transformation, and centre of mass distance between features. A dosimetric evaluation was performed to evaluate the clinical significance of the registrations errors in the application proposed and to identify the limitations of the approximations used. Dose calculations for the same intensity-modulated radiation therapy plan on the deformed CT and replan CT were compared. Dose distributions were compared in terms of dose differences (DD), gamma analysis, target coverage, and dose volume histograms (DVHs). Doses calculated in a rigidly aligned CT and directly in an extended CBCT were also evaluated. Results: A mean value of 0.850 in DSC was achieved in overlap between manually delineated and warped features, with the distance between surfaces being less than 2 mm on over 90% of the pixels. Deformable registration was clearly superior to rigid registration in mapping identical structures between the two datasets. The dose recalculated in the deformed CT is a good match to the dose calculated on a replan CT. The DD is smaller than 2% of the prescribed dose on 90% of the body's voxels and it passes a 2% and 2 mm gamma-test on over 95% of the voxels. Target coverage similarity was assessed in terms of the 95%-isodose volumes. A mean value of 0.962 was obtained for the DSC, while the distance between surfaces is less than 2 mm in 95.4% of the pixels. The method proposed provided adequate dose estimation, closer to the gold standard than the other two approaches. Differences in DVH curves were mainly due to differences in the OARs definition (manual vs warped) and not due to differences in dose estimation (dose calculated in replan CT vs dose calculated in deformed CT). Conclusions: Deforming a planning CT to match a daily CBCT provides the tools needed for the calculation of the “dose of the day” without the need to acquire a new CT. The initial clinical application of our method will be weekly offline calculations of the “dose of the day,” and use this information to inform adaptive radiotherapy (ART). The work here presented is a first step into a full implementation of a “dose-driven” online ART.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Dosimetric comparison of flattened and unflattened beams for stereotactic ablative radiotherapy of stage I non-small cell lung cancer

Purpose: To compare contribution and accuracy of delivery for two flattening filter free (FFF) beams of the nominal energy 6 and 10 MV and a 6 MV flattened beam for early stage lung cancer. Methods: For each of 11 patients with stage I nonsmall cell lung cancer three volumetric modulated arc therapy plans were prepared utilizing a 6 MV flattened photon beam (X6FF) and two nonflattened beams of nominal energy 6 and 10 MV (X6FFF, X10FFF). Optimization constraints were set to produce dose distributions that meet the criteria of the RTOG-0915 protocol. The radiation schedule used for plan comparison in all patients was 50 Gy in five fractions. Dosimetric parameters of planning target volume (PTV) and organs-at-risk and delivery times were assessed and compared. All plans were subject to verification using Delta{sup 4} unit (Scandidos, Sweden) and absolutely calibrated gafchromic films in a thorax phantom. Results: All plans had a qualitatively comparable outcome. Obtained dose distributions were conformal (CI < 1.17) and exhibited a steep dose fall-off outside the PTV. The ratio of monitor units for FFF versus FF plans in the authors' study ranged from 0.95 to 1.21 and from 0.93 to 1.25 for X6FFF/X6FF and X10FFF/X6FF comparisons, respectively. The ratio systematically increased with increasing size of the PTV (up to +25% for 150 cm{sup 3} PTV). Yet the integral dose to healthy tissue did not follow this trend. Comparison of cumulative dose volume histograms for a patient's body showed that X6FFF plans exhibit improved conformity and reduced the volume of tissue that received more than 50% of the prescription dose. Parameters related to dose gradient showed statistically significant improvement. CI{sub 50%}, CI{sub 60%}, CI{sub 80%}, and CI{sub 100%} were on average reduced by 4.6% (p < 0.001), 4.6% (p = 0.002), 3.1% (p = 0.002), and 1.2% (p = 0.039), respectively. Gradient measure was on average reduced by 4.2% (p < 0.001). Due to dose reduction in the surrounding lung tissue, the V{sub 20} {sub Gy} and V{sub 12.5} {sub Gy} were reduced by 5.5% (p = 0.002) and 4.5% (p < 0.001). These dosimetric improvements in the fall-off were not observed for the X10FFF plans. Differences in sparing of normal tissues were not found to be statistically significant for either of the two FFF beams. Mean beam-on times were 111 s (2SD = 11 s) for X10FFF, 128 s (2SD = 19 s) for X6FFF, and X6FF plans required on average 269 s (2SD = 71 s). While the mean dose rate was 1555 ± 264 and 1368 ± 63 MU/min, for X10FFF and X6FFF, plans using the conventional X6FF were delivered with the constant maximum dose rate of 600 MU/min. Verification of all plans showed acceptable and comparable results for all plans in homogeneous as well as heterogeneous phantoms. Mean GS (3%, 2 mm) using the Delta{sup 4} phantom were 98.9% (2SD = 3.2%), 99.2% (2SD = 2.3%), and 99.2% (2SD = 2.3%) for X6FFF, X6FF, and X10FFF modalities. Verification using a thorax phantom showed GS > 98% in all cases. Conclusions: The use of FFF beams for stereotactic radiation therapy of nonsmall cell lung cancer patients yielded dose distributions qualitatively comparable to flattened beams and significantly reduced treatment delivery time. Utilizing the X6FFF beam improved conformity of dose distribution. On the other hand, X10FFF beam offered a slight improvement in treatment efficiency, and lower skin and peripheral dose. All effects were relatively small.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Beam geometry selection using sequential beam addition

Purpose: The selection of optimal beam geometry has been of interest since the inception of conformal radiotherapy. The authors report on sequential beam addition, a simple beam geometry selection method, for intensity modulated radiation therapy. Methods: The sequential beam addition algorithm (SBA) requires definition of an objective function (score) and a set of candidate beam geometries (pool). In the first iteration, the optimal score is determined for each beam in the pool and the beam with the best score selected. In the next iteration, the optimal score is calculated for each beam remaining in the pool combined with the beam selected in the first iteration, and the best scoring beam is selected. The process is repeated until the desired number of beams is reached. The authors selected three treatment sites, breast, lung, and brain, and determined beam arrangements for up to 11 beams from a pool comprised of 25 equiangular transverse beams. For the brain, arrangements were additionally selected from a pool of 22 noncoplanar beams. Scores were determined for geometries comprised equiangular transverse beams (EQA), as well as two tangential beams for the breast case. Results: In all cases, SBA resulted in scores superior to EQA. The breast case had the strongest dependence on beam geometry, for which only the 7-beam EQA geometry had a score better than the two tangential beams, whereas all SBA geometries with more than two beams were superior. In the lung case, EQA and SBA scores monotonically improved with increasing number of beams; however, SBA required fewer beams to achieve scores equivalent to EQA. For the brain case, SBA with a coplanar pool was equivalent to EQA, while the noncoplanar pool resulted in slightly better scores; however, the dose-volume histograms demonstrated that the differences were not clinically significant. Conclusions: For situations in which beam geometry has a significant effect on the objective function, SBA can identify arrangements equivalent to equiangular geometries but using fewer beams. Furthermore, SBA provides the value of the objective function as the number of beams is increased, allowing the planner to select the minimal beam number that achieves the clinical goals. The method is simple to implement and could readily be incorporated into an existing optimization system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Helical tomotherapy with dynamic running-start-stop delivery compared to conventional tomotherapy delivery

Purpose: Despite superior target dose uniformity, helical tomotherapy{sup ®} (HT) may involve a trade-off between longitudinal dose conformity and beam-on time (BOT), due to the limitation of only three available jaw sizes with the conventional HT (1.0, 2.5, and 5.0 cm). The recently introduced dynamic running-start-stop (RSS) delivery allows smaller jaw opening at the superior and inferior ends of the target when a sharp penumbra is needed. This study compared the dosimetric performance of RSS delivery with the fixed jaw HT delivery. Methods: Twenty patient cases were selected and deidentified prior to treatment planning, including 16 common clinical cases (brain, head and neck (HN), lung, and prostate) and four special cases of whole brain with hippocampus avoidance (WBHA) that require a high degree of dose modulation. HT plans were generated for common clinical cases using the fixed 2.5 cm jaw width (HT2.5) and WBHA cases using 1.0 cm (HT1.0). The jaw widths for RSS were preset with a larger size (RSS5.0 vs HT2.5 and RSS2.5 vs HT1.0). Both delivery techniques were planned based on identical contours, prescriptions, and planning objectives. Dose indices for targets and critical organs were compared using dose-volume histograms, BOT, and monitor units. Results: The average BOT was reduced from 4.8 min with HT2.5 to 2.5 min with RSS5.0. Target dose homogeneity with RSS5.0 was shown comparable to HT2.5 for common clinical sites. Superior normal tissue sparing was observed in RSS5.0 for optic nerves and optic chiasm in brain and HN cases. RSS5.0 demonstrated improved dose sparing for cord and esophagus in lung cases, as well as penile bulb in prostate cases. The mean body dose was comparable for both techniques. For the WBHA cases, the target homogeneity was significantly degraded in RSS2.5 without distinct dose sparing for hippocampus, compared to HT1.0. Conclusions: Compared to the fixed jaw HT delivery, RSS combined with a larger jaw width provides faster treatment delivery and improved cranial-caudal target dose conformity. The target coverage achieved by RSS with a large jaw width is comparable to the fixed jaw HT delivery for common cancer sites, but may deteriorate for cases where complex geometry is present in the middle part of the target.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

MRI-based treatment planning with pseudo CT generated through atlas registration

Purpose: To evaluate the feasibility and accuracy of magnetic resonance imaging (MRI)-based treatment planning using pseudo CTs generated through atlas registration. Methods: A pseudo CT, providing electron density information for dose calculation, was generated by deforming atlas CT images previously acquired on other patients. The authors tested 4 schemes of synthesizing a pseudo CT from single or multiple deformed atlas images: use of a single arbitrarily selected atlas, arithmetic mean process using 6 atlases, and pattern recognition with Gaussian process (PRGP) using 6 or 12 atlases. The required deformation for atlas CT images was derived from a nonlinear registration of conjugated atlas MR images to that of the patient of interest. The contrasts of atlas MR images were adjusted by histogram matching to reduce the effect of different sets of acquisition parameters. For comparison, the authors also tested a simple scheme assigning the Hounsfield unit of water to the entire patient volume. All pseudo CT generating schemes were applied to 14 patients with common pediatric brain tumors. The image similarity of real patient-specific CT and pseudo CTs constructed by different schemes was compared. Differences in computation times were also calculated. The real CT in the treatment planning system was replaced with the pseudo CT, and the dose distribution was recalculated to determine the difference. Results: The atlas approach generally performed better than assigning a bulk CT number to the entire patient volume. Comparing atlas-based schemes, those using multiple atlases outperformed the single atlas scheme. For multiple atlas schemes, the pseudo CTs were similar to the real CTs (correlation coefficient, 0.787–0.819). The calculated dose distribution was in close agreement with the original dose. Nearly the entire patient volume (98.3%–98.7%) satisfied the criteria of chi-evaluation (<2% maximum dose and 2 mm range). The dose to 95% of the volume and the percentage of volume receiving at least 95% of the prescription dose in the planning target volume differed from the original values by less than 2% of the prescription dose (root-mean-square, RMS < 1%). The PRGP scheme did not perform better than the arithmetic mean process with the same number of atlases. Increasing the number of atlases from 6 to 12 often resulted in improvements, but statistical significance was not always found. Conclusions: MRI-based treatment planning with pseudo CTs generated through atlas registration is feasible for pediatric brain tumor patients. The doses calculated from pseudo CTs agreed well with those from real CTs, showing dosimetric accuracy within 2% for the PTV when multiple atlases were used. The arithmetic mean process may be a reasonable choice over PRGP for the synthesis scheme considering performance and computational costs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

HDRMC, an accelerated Monte Carlo dose calculator for high dose rate brachytherapy with CT-compatible applicators

Purpose: To present a new accelerated Monte Carlo code for CT-based dose calculations in high dose rate (HDR) brachytherapy. The new code (HDRMC) accounts for both tissue and nontissue heterogeneities (applicator and contrast medium). Methods: HDRMC uses a fast ray-tracing technique and detailed physics algorithms to transport photons through a 3D mesh of voxels representing the patient anatomy with applicator and contrast medium included. A precalculated phase space file for the{sup 192}Ir source is used as source term. HDRM is calibrated to calculated absolute dose for real plans. A postprocessing technique is used to include the exact density and composition of nontissue heterogeneities in the 3D phantom. Dwell positions and angular orientations of the source are reconstructed using data from the treatment planning system (TPS). Structure contours are also imported from the TPS to recalculate dose-volume histograms. Results: HDRMC was first benchmarked against the MCNP5 code for a single source in homogenous water and for a loaded gynecologic applicator in water. The accuracy of the voxel-based applicator model used in HDRMC was also verified by comparing 3D dose distributions and dose-volume parameters obtained using 1-mm{sup 3} versus 2-mm{sup 3} phantom resolutions. HDRMC can calculate the 3D dose distribution for a typical HDR cervix case with 2-mm resolution in 5 min on a single CPU. Examples of heterogeneity effects for two clinical cases (cervix and esophagus) were demonstrated using HDRMC. The neglect of tissue heterogeneity for the esophageal case leads to the overestimate of CTV D90, CTV D100, and spinal cord maximum dose by 3.2%, 3.9%, and 3.6%, respectively. Conclusions: A fast Monte Carlo code for CT-based dose calculations which does not require a prebuilt applicator model is developed for those HDR brachytherapy treatments that use CT-compatible applicators. Tissue and nontissue heterogeneities should be taken into account in modern HDR brachytherapy planning.

61 RADIATION PROTECTION AND DOSIMETRY↗