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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2014

Assessment of individual organ doses in a realistic human phantom from neutron and gamma stimulated spectroscopy of the breast and liver

Purpose: Understanding the radiation dose to a patient is essential when considering the use of an ionizing diagnostic imaging test for clinical diagnosis and screening. Using Monte Carlo simulations, the authors estimated the three-dimensional organ-dose distribution from neutron and gamma irradiation of the male liver, female liver, and female breasts for neutron- and gamma-stimulated spectroscopic imaging. Methods: Monte Carlo simulations were developed using the Geant4 GATE application and a voxelized XCAT human phantom. A male and a female whole body XCAT phantom was voxelized into 256 × 256 × 600 voxels (3.125 × 3.125 × 3.125 mm{sup 3}). A monoenergetic rectangular beam of 5.0 MeV neutrons or 7.0 MeV photons was made incident on a 2 cm thick slice of the phantom. The beam was rotated at eight different angles around the phantom ranging from 0° to 180°. Absorbed dose was calculated for each individual organ in the body and dose volume histograms were computed to analyze the absolute and relative doses in each organ. Results: The neutron irradiations of the liver showed the highest organ dose absorption in the liver, with appreciably lower doses in other proximal organs. The dose distribution within the irradiated slice exhibited substantial attenuation with increasing depth along the beam path, attenuating to ∼15% of the maximum value at the beam exit side. The gamma irradiation of the liver imparted the highest organ dose to the stomach wall. The dose distribution from the gammas showed a region of dose buildup at the beam entrance, followed by a relatively uniform dose distribution to all of the deep tissue structures, attenuating to ∼75% of the maximum value at the beam exit side. For the breast scans, both the neutron and gamma irradiation registered maximum organ doses in the breasts, with all other organs receiving less than 1% of the breast dose. Effective doses ranged from 0.22 to 0.37 mSv for the neutron scans and 41 to 66 mSv for the gamma scans. Conclusions: Neutron and gamma irradiation of a primary target organ was found to impart the majority of the total dose to the primary target organ (and other large organs) within the beam plane and considerably lower dose to proximal organs outside of the beam. These results also indicate that despite the use of a highly scattering particle such as a neutron, the dose from neutron stimulated emission computed tomography scans is on par with other clinical imaging techniques such as x-ray computed tomography (x-ray CT). Given the high nonuniformity in the dose across an organ during the neutron scan, care must be taken when computing average doses from neutron irradiations. The effective doses from neutron scanning were found to be comparable to x-ray CT. Further technique modifications are needed to reduce the effective dose levels from the gamma scans.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

A critical evaluation of worst case optimization methods for robust intensity-modulated proton therapy planning

Purpose: To critically evaluate and compare three worst case optimization methods that have been previously employed to generate intensity-modulated proton therapy treatment plans that are robust against systematic errors. The goal of the evaluation is to identify circumstances when the methods behave differently and to describe the mechanism behind the differences when they occur. Methods: The worst case methods optimize plans to perform as well as possible under the worst case scenario that can physically occur (composite worst case), the combination of the worst case scenarios for each objective constituent considered independently (objectivewise worst case), and the combination of the worst case scenarios for each voxel considered independently (voxelwise worst case). These three methods were assessed with respect to treatment planning for prostate under systematic setup uncertainty. An equivalence with probabilistic optimization was used to identify the scenarios that determine the outcome of the optimization. Results: If the conflict between target coverage and normal tissue sparing is small and no dose-volume histogram (DVH) constraints are present, then all three methods yield robust plans. Otherwise, they all have their shortcomings: Composite worst case led to unnecessarily low plan quality in boundary scenarios that were less difficult than the worst case ones. Objectivewise worst case generally led to nonrobust plans. Voxelwise worst case led to overly conservative plans with respect to DVH constraints, which resulted in excessive dose to normal tissue, and less sharp dose fall-off than the other two methods. Conclusions: The three worst case methods have clearly different behaviors. These behaviors can be understood from which scenarios that are active in the optimization. No particular method is superior to the others under all circumstances: composite worst case is suitable if the conflicts are not very severe or there are DVH constraints whereas voxelwise worst case is advantageous if there are severe conflicts but no DVH constraints. The advantages of composite and voxelwise worst case outweigh those of objectivewise worst case.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

ARCHER{sub RT} – A GPU-based and photon-electron coupled Monte Carlo dose computing engine for radiation therapy: Software development and application to helical tomotherapy

Purpose: Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHER{sub RT} is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy® cases: the prostate, lung, and head and neck. Methods: To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHER{sub RT}. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHER{sub RT} and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU. Results: For the water phantom, the depth dose curve and dose profiles from ARCHER{sub RT} agree well with DOSXYZnrc. For clinical cases, results from ARCHER{sub RT} are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head and neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHER{sub RT} achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7–1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9–13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500–800 s. Conclusions: ARCHER{sub RT} was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Independent calculation of dose distributions for helical tomotherapy using a conventional treatment planning system

Purpose: The dosimetric verification of treatment plans in helical tomotherapy usually is carried out via verification measurements. In this study, a method for independent dose calculation of tomotherapy treatment plans is presented, that uses a conventional treatment planning system with a pencil kernel dose calculation algorithm for generation of verification dose distributions based on patient CT data. Methods: A pencil beam algorithm that directly uses measured beam data was configured for dose calculation for a tomotherapy machine. Tomotherapy treatment plans were converted into a format readable by an in-house treatment planning system by assigning each projection to one static treatment field and shifting the calculation isocenter for each field in order to account for the couch movement. The modulation of the fluence for each projection is read out of the delivery sinogram, and with the kernel-based dose calculation, this information can directly be used for dose calculation without the need for decomposition of the sinogram. The sinogram values are only corrected for leaf output and leaf latency. Using the converted treatment plans, dose was recalculated with the independent treatment planning system. Multiple treatment plans ranging from simple static fields to real patient treatment plans were calculated using the new approach and either compared to actual measurements or the 3D dose distribution calculated by the tomotherapy treatment planning system. In addition, dose–volume histograms were calculated for the patient plans. Results: Except for minor deviations at the maximum field size, the pencil beam dose calculation for static beams agreed with measurements in a water tank within 2%/2 mm. A mean deviation to point dose measurements in the cheese phantom of 0.89% ± 0.81% was found for unmodulated helical plans. A mean voxel-based deviation of −0.67% ± 1.11% for all voxels in the respective high dose region (dose values >80%), and a mean local voxel-based deviation of −2.41% ± 0.75% for all voxels with dose values >20% were found for 11 modulated plans in the cheese phantom. Averaged over nine patient plans, the deviations amounted to −0.14% ± 1.97% (voxels >80%) and −0.95% ± 2.27% (>20%, local deviations). For a lung case, mean voxel-based deviations of more than 4% were found, while for all other patient plans, all mean voxel-based deviations were within ±2.4%. Conclusions: The presented method is suitable for independent dose calculation for helical tomotherapy within the known limitations of the pencil beam algorithm. It can serve as verification of the primary dose calculation and thereby reduce the need for time-consuming measurements. By using the patient anatomy and generating full 3D dose data, and combined with measurements of additional machine parameters, it can substantially contribute to overall patient safety.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-28: Evaluating the Organ Dose From Computed Tomography Using Monte Carlo Calculations

Purpose: To evaluate organ doses from computed tomography (CT) using Monte Carlo (MC) calculations. Methods: A Philips Brilliance CT scanner (64 slice) was simulated using the GMctdospp (IMPS, Germany) based on the EGSnrc user code. The X-ray spectra and a bowtie filter for MC simulations were determined to coincide with measurements of half-value layer (HVL) and off-center ratio (OCR) profile in air. The MC dose was calibrated from absorbed dose measurements using a Farmer chamber and a cylindrical water phantom. The dose distribution from CT was calculated using patient CT images and organ doses were evaluated from dose volume histograms. Results: The HVLs of Al at 80, 100, and 120 kV were 6.3, 7.7, and 8.7 mm, respectively. The calculated HVLs agreed with measurements within 0.3%. The calculated and measured OCR profiles agreed within 3%. For adult head scans (CTDIvol) =51.4 mGy), mean doses for brain stem, eye, and eye lens were 23.2, 34.2, and 37.6 mGy, respectively. For pediatric head scans (CTDIvol =35.6 mGy), mean doses for brain stem, eye, and eye lens were 19.3, 24.5, and 26.8 mGy, respectively. For adult chest scans (CTDIvol=19.0 mGy), mean doses for lung, heart, and spinal cord were 21.1, 22.0, and 15.5 mGy, respectively. For adult abdominal scans (CTDIvol=14.4 mGy), the mean doses for kidney, liver, pancreas, spleen, and spinal cord were 17.4, 16.5, 16.8, 16.8, and 13.1 mGy, respectively. For pediatric abdominal scans (CTDIvol=6.76 mGy), mean doses for kidney, liver, pancreas, spleen, and spinal cord were 8.24, 8.90, 8.17, 8.31, and 6.73 mGy, respectively. In head scan, organ doses were considerably different from CTDIvol values. Conclusion: MC dose distributions calculated by using patient CT images are useful to evaluate organ doses absorbed to individual patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-63: Quantitative Evaluation of the Effects of Orthopedic Metal Artifact Reduction (OMAR) Software On CT Images for Radiotherapy Simulation

Purpose: CT simulation for patients with metal implants can often be challenging due to artifacts that obscure tumor/target delineation and normal organ definition. Our objective was to evaluate the effectiveness of Orthopedic Metal Artifact Reduction (OMAR), a commercially available software, in reducing metal-induced artifacts and its effect on computed dose during treatment planning. Methods: CT images of water surrounding metallic cylindrical rods made of aluminum, copper and iron were studied in terms of Hounsfield Units (HU) spread. Metal-induced artifacts were characterized in terms of HU/Volume Histogram (HVH) using the Pinnacle treatment planning system. Effects of OMAR on enhancing our ability to delineate organs on CT and subsequent dose computation were examined in nine (9) patients with hip implants and two (2) patients with breast tissue expanders. Results: Our study characterized water at 1000 HU with a standard deviation (SD) of about 20 HU. The HVHs allowed us to evaluate how the presence of metal changed the HU spread. For example, introducing a 2.54 cm diameter copper rod in water increased the SD in HU of the surrounding water from 20 to 209, representing an increase in artifacts. Subsequent use of OMAR brought the SD down to 78. Aluminum produced least artifacts whereas Iron showed largest amount of artifacts. In general, an increase in kVp and mA during CT scanning showed better effectiveness of OMAR in reducing artifacts. Our dose analysis showed that some isodose contours shifted by several mm with OMAR but infrequently and were nonsignificant in planning process. Computed volumes of various dose levels showed <2% change. Conclusions: In our experience, OMAR software greatly reduced the metal-induced CT artifacts for the majority of patients with implants, thereby improving our ability to delineate tumor and surrounding organs. OMAR had a clinically negligible effect on computed dose within tissues. Partially funded by unrestricted educational grant from Philips.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-92: Is Photon Starvation Preventing Metal Artifact Reduction Algorithm From Working in KVCT?

Purpose: High density/high atomic number metallic objects create shading and streaking metal artifacts in the CT image that can cause inaccurate delineation of anatomical structures or inaccurate radiation dose calculation. A modified iterative maximum-likelihood polychromatic algorithm for CT (mIMPACT) that models the energy response of detectors, photon interaction processes and beam polychromaticity has successfully reduced metal artifacts in MVCT. Our extension of mIMPACT in kVCT did not significantly reduce metal artifacts for high density metal like steel. We hypothesize that photon starvation may result in the measured data in a commercial kVCT imaging beam. Methods: We measured attenuation of a range of steel plate thicknesses, sandwiched between two 12cm thick solid water blocks, using a Phillips Big Bore CTTM scanner in scout acquisition mode with 120kVp and 200mAs. The transmitted signal (y) was normalized to the air scan signal (y{sub 0}) to get attenuation [i.e., ln(y/y{sub 0})] data for a detector. Results: Below steel plate thickness of 13.4mm, the variations in measured attenuation as a function of view number are characterized by a quantum noise and show increased attenuation with metal thickness. On or above this thickness the attenuation shows discrete levels in addition to the quantum noise. Some views have saturated attenuation value. The histograms of the measured attenuation for up to 36.7mm of steel show this trend. The detector signal is so small that the quantization levels in the analog to digital (A-to-D) converter are visible, a clear indication of photon starvation. Conclusion: Photons reaching the kVCT detector after passing through a thick metal plate are either so low in number that the signal measured has large quantum noise, or are completely absorbed inside the plate creating photon starvation. This is un-interpretable by the mIMPACT algorithm and cannot reduce metal artifacts in kVCT for certain realistic thicknesses of steel hip implants. Moti Raj Paudel is supported by the Vanier Canada Graduate Scholarship, the Endowed Graduate Scholarship in Oncology and the Dissertation Fellowship at the University of Alberta. The authors acknowledge the CIHR operating grant number MOP 53254.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-J-20: Adaptive Aperture Morphing for Online Correction for Prostate Cancer Radiotherapy

Purpose: Online adaptive aperture morphing is desirable over translational couch shifts to accommodate not only the target position variation but also anatomic changes (rotation, deformation, and relation of target to organ-atrisks). We proposed quick and reliable method for adapting segment aperture leaves for IMRT treatment of prostate. Methods: The proposed method consists of following steps: (1) delineate the contours of prostate, SV, bladder and rectum on kV-CBCT; (2) determine prostate displacement from the rigid body registration of the contoured prostate manifested on the reference CT and the CBCT; (3) adapt the MLC segment apertures obtained from the pre-treatment IMRT planning to accommodate the shifts as well as anatomic changes. The MLC aperture adaptive algorithm involves two steps; first move the whole aperture according to prostate translational/rotational shifts, and secondly fine-tune the aperture shape to maintain the spatial relationship between the planning target contour and the MLC aperture to the daily target contour. Feasibility of this method was evaluated retrospectively on a seven-field IMRT treatment of prostate cancer patient by comparing dose volume histograms of the original plan and the aperture-adjusted plan, with/without additional segments weight optimization (SWO), on two daily treatment CBCTs selected with relative large motion and rotation. Results: For first daily treatment, the prostate rotation was significant (12degree around lateral-axis). With apertureadjusted plan, the D95 to the target was improved 25% and rectum dose (D30, D40) was reduced 20% relative to original plan on daily volumes. For second treatment-fraction, (lateral shift = 6.7mm), after adjustment target D95 improved by 3% and bladder dose (D30, maximum dose) was reduced by 1%. For both cases, extra SWO did not provide significant improvement. Conclusion: The proposed method of adapting segment apertures is promising in treatment position correction, including target translational displacement, rotation and deformation. Additional SWO could improve ROIs dose distribution.

60 APPLIED LIFE SCIENCES↗