Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An optical to IR sky brightness model for the LSST

To optimize the observing strategy of a large survey such as the LSST, one needs an accurate model of the night sky emission spectrum across a range of atmospheric conditions and from the near-UV to the near-IR. We have used the ESO SkyCalc Sky Model Calculator to construct a library of template spectra for the Chilean night sky. The ESO model includes emission from the upper and lower atmosphere, scattered starlight, scattered moonlight, and zodiacal light. We have then extended the ESO templates with an empirical fit to the twilight sky emission as measured by a Canon all-sky camera installed at the LSST site. With the ESO templates and our twilight model we can quickly interpolate to any arbitrary sky position and date and return the full sky spectrum or surface brightness magnitudes in the LSST filter system. Comparing our model to all-sky observations, we find typical residual RMS values of ±0.2-0.3 magnitudes per square arcsecond.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Performance of 12 DIR algorithms in low-contrast regions for mass and density conserving deformation

Purpose: Deformable image registration (DIR) has become a key tool for adaptive radiotherapy to account for inter- and intrafraction organ deformation. Of contemporary interest, the application to deformable dose accumulation requires accurate deformation even in low contrast regions where dose gradients may exist within near-uniform tissues. One expects high-contrast features to generally be deformed more accurately by DIR algorithms. The authors systematically assess the accuracy of 12 DIR algorithms and quantitatively examine, in particular, low-contrast regions, where accuracy has not previously been established.Methods: This work investigates DIR algorithms in three dimensions using deformable gel (DEFGEL) [U. J. Yeo, M. L. Taylor, L. Dunn, R. L. Smith, T. Kron, and R. D. Franich, “A novel methodology for 3D deformable dosimetry,” Med. Phys. 39, 2203–2213 (2012)], for application to mass- and density-conserving deformations. CT images of DEFGEL phantoms with 16 fiducial markers (FMs) implanted were acquired in deformed and undeformed states for three different representative deformation geometries. Nonrigid image registration was performed using 12 common algorithms in the public domain. The optimum parameter setup was identified for each algorithm and each was tested for deformation accuracy in three scenarios: (I) original images of the DEFGEL with 16 FMs; (II) images with eight of the FMs mathematically erased; and (III) images with all FMs mathematically erased. The deformation vector fields obtained for scenarios II and III were then applied to the original images containing all 16 FMs. The locations of the FMs estimated by the algorithms were compared to actual locations determined by CT imaging. The accuracy of the algorithms was assessed by evaluation of three-dimensional vectors between true marker locations and predicted marker locations.Results: The mean magnitude of 16 error vectors per sample ranged from 0.3 to 3.7, 1.0 to 6.3, and 1.3 to 7.5 mm across algorithms for scenarios I to III, respectively. The greatest accuracy was exhibited by the original Horn and Schunck optical flow algorithm. In this case, for scenario III (erased FMs not contributing to driving the DIR calculation), the mean error was half that of the modified demons algorithm (which exhibited the greatest error), across all deformations. Some algorithms failed to reproduce the geometry at all, while others accurately deformed high contrast features but not low-contrast regions—indicating poor interpolation between landmarks.Conclusions: The accuracy of DIR algorithms was quantitatively evaluated using a tissue equivalent, mass, and density conserving DEFGEL phantom. For the model studied, optical flow algorithms performed better than demons algorithms, with the original Horn and Schunck performing best. The degree of error is influenced more by the magnitude of displacement than the geometric complexity of the deformation. As might be expected, deformation is estimated less accurately for low-contrast regions than for high-contrast features, and the method presented here allows quantitative analysis of the differences. The evaluation of registration accuracy through observation of the same high contrast features that drive the DIR calculation is shown to be circular and hence misleading.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

An automated voxelized dosimetry tool for radionuclide therapy based on serial quantitative SPECT/CT imaging

Purpose: To create an accurate map of the distribution of radiation dose deposition in healthy and target tissues during radionuclide therapy.Methods: Serial quantitative SPECT/CT images were acquired at 4, 24, and 72 h for 28 {sup 177}Lu-octreotate peptide receptor radionuclide therapy (PRRT) administrations in 17 patients with advanced neuroendocrine tumors. Deformable image registration was combined with an in-house programming algorithm to interpolate pharmacokinetic uptake and clearance at a voxel level. The resultant cumulated activity image series are comprised of values representing the total number of decays within each voxel's volume. For PRRT, cumulated activity was translated to absorbed dose based on Monte Carlo-determined voxel S-values at a combination of long and short ranges. These dosimetric image sets were compared for mean radiation absorbed dose to at-risk organs using a conventional MIRD protocol (OLINDA 1.1).Results: Absorbed dose values to solid organs (liver, kidneys, and spleen) were within 10% using both techniques. Dose estimates to marrow were greater using the voxelized protocol, attributed to the software incorporating crossfire effect from nearby tumor volumes.Conclusions: The technique presented offers an efficient, automated tool for PRRT dosimetry based on serial post-therapy imaging. Following retrospective analysis, this method of high-resolution dosimetry may allow physicians to prescribe activity based on required dose to tumor volume or radiation limits to healthy tissue in individual patients.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Resolution enhancement of lung 4D-CT via group-sparsity

Purpose: 4D-CT typically delivers more accurate information about anatomical structures in the lung, over 3D-CT, due to its ability to capture visual information of the lung motion across different respiratory phases. This helps to better determine the dose during radiation therapy for lung cancer. However, a critical concern with 4D-CT that substantially compromises this advantage is the low superior-inferior resolution due to less number of acquired slices, in order to control the CT radiation dose. To address this limitation, the authors propose an approach to reconstruct missing intermediate slices, so as to improve the superior-inferior resolution.Methods: In this method the authors exploit the observation that sampling information across respiratory phases in 4D-CT can be complimentary due to lung motion. The authors’ approach uses this locally complimentary information across phases in a patch-based sparse-representation framework. Moreover, unlike some recent approaches that treat local patches independently, the authors’ approach employs the group-sparsity framework that imposes neighborhood and similarity constraints between patches. This helps in mitigating the trade-off between noise robustness and structure preservation, which is an important consideration in resolution enhancement. The authors discuss the regularizing ability of group-sparsity, which helps in reducing the effect of noise and enables better structural localization and enhancement.Results: The authors perform extensive experiments on the publicly available DIR-Lab Lung 4D-CT dataset [R. Castillo, E. Castillo, R. Guerra, V. Johnson, T. McPhail, A. Garg, and T. Guerrero, “A framework for evaluation of deformable image registration spatial accuracy using large landmark point sets,” Phys. Med. Biol. 54, 1849–1870 (2009)]. First, the authors carry out empirical parametric analysis of some important parameters in their approach. The authors then demonstrate, qualitatively as well as quantitatively, the ability of their approach to achieve more accurate and better localized results over bicubic interpolation as well as a related state-of-the-art approach. The authors also show results on some datasets with tumor, to further emphasize the clinical importance of their method.Conclusions: The authors have proposed to improve the superior-inferior resolution of 4D-CT by estimating intermediate slices. The authors’ approach exploits neighboring constraints in the group-sparsity framework, toward the goal of achieving better localization and noise robustness. The authors’ results are encouraging, and positively demonstrate the role of group-sparsity for 4D-CT resolution enhancement.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

A cascaded model of spectral distortions due to spectral response effects and pulse pileup effects in a photon-counting x-ray detector for CT

Purpose: Energy discriminating, photon-counting detectors (PCDs) are an emerging technology for computed tomography (CT) with various potential benefits for clinical CT. The photon energies measured by PCDs can be distorted due to the interactions of a photon with the detector and the interaction of multiple coincident photons. These effects result in distorted recorded x-ray spectra which may lead to artifacts in reconstructed CT images and inaccuracies in tissue identification. Model-based compensation techniques have the potential to account for the distortion effects. This approach requires only a small number of parameters and is applicable to a wide range of spectra and count rates, but it needs an accurate model of the spectral distortions occurring in PCDs. The purpose of this study was to develop a model of those spectral distortions and to evaluate the model using a PCD (model DXMCT-1; DxRay, Inc., Northridge, CA) and various x-ray spectra in a wide range of count rates. Methods: The authors hypothesize that the complex phenomena of spectral distortions can be modeled by: (1) separating them into count-rate independent factors that we call the spectral response effects (SRE), and count-rate dependent factors that we call the pulse pileup effects (PPE), (2) developing separate models for SRE and PPE, and (3) cascading the SRE and PPE models into a combined SRE+PPE model that describes PCD distortions at both low and high count rates. The SRE model describes the probability distribution of the recorded spectrum, with a photo peak and a continuum tail, given the incident photon energy. Model parameters were obtained from calibration measurements with three radioisotopes and then interpolated linearly for other energies. The PPE model used was developed in the authors’ previous work [K. Taguchi et al., “Modeling the performance of a photon counting x-ray detector for CT: Energy response and pulse pileup effects,” Med. Phys. 38(2), 1089–1102 (2011)]. The agreement between the x-ray spectra calculated by the cascaded SRE+PPE model and the measured spectra was evaluated for various levels of deadtime loss ratios (DLR) and incident spectral shapes, realized using different attenuators, in terms of the weighted coefficient of variation (COV{sub W}), i.e., the root mean square difference weighted by the statistical errors of the data and divided by the mean. Results: At low count rates, when DLR < 10%, the distorted spectra measured by the DXMCT-1 were in agreement with those calculated by SRE only, with COV{sub W}'s less than 4%. At higher count rates, the measured spectra were also in agreement with the ones calculated by the cascaded SRE+PPE model; with PMMA as attenuator, COV{sub W} was 5.6% at a DLR of 22% and as small as 6.7% for a DLR as high as 55%. Conclusions: The x-ray spectra calculated by the proposed model agreed with the measured spectra over a wide range of count rates and spectral shapes. The SRE model predicted the distorted, recorded spectra with low count rates over various types and thicknesses of attenuators. The study also validated the hypothesis that the complex spectral distortions in a PCD can be adequately modeled by cascading the count-rate independent SRE and the count-rate dependent PPE.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-D-19A-02: Electron and Photon Absorbed Fractions for Tumors of Varying Sizes and Compositions

Purpose: To calculate absorbed fractions for mono-energetic photons and electrons in tumors of varying compositions using Monte Carlo simulations in MCNPX. Although tumor dosimetry has been previously investigated, these studies are very limited as they only consider absorbed fractions for soft-tissue tumors. Methods: The tumors were modeled as spheres with radii ranging from 0.10 cm to 6.0 cm and with compositions varying from 100% soft tissue to 100% bone. The energies of both the photons and electrons were varied from 10 keV to 10 MeV and were homogenously distributed throughout the tumor volume. Furthermore, this investigation addresses the issue of spherical versus elliptical tumors. Both prolate and oblate spheroid tumors of different compositions were modeled, and absorbed fractions were calculated for various electron and photon energies. Results: The data clearly shows an absorbed fraction dependence on tumor composition. For example, a soft-tissue model for a 3 MeV electron emitted in a 1 cm diameter bone tumor would have an 83% error, and this same assumption for a 500 keV photon would yield a 74% error. Ultimately, empirical fits were created for each of the five material compositions in order to facilitate the absorbed fraction calculation, requiring only the tumor size and particle energy. Furthermore, the data shows that absorbed fractions for moderate spheroids can be well approximated by spherical tumors of equal mass to within 8%, but in the extreme cases where the spheroid resembles more of a disk, the errors can be as high as 30%. Conclusion: This comprehensive data set is most valuable for nuclear medicine dosimetry because it incorporates particle type, particle energy, tumor size, and tumor composition. Although mono-energetic particles were modeled, absorbed fractions and S-values may be calculated for any radionuclide via linear interpolation, as long as the particle energies or spectra are known.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-95: Personalized Radiography Technical Parameters for Each Patient and Exam

Purpose: To determine exact electrical parameters (kV, mAs) a radiological technologist shall use taking account the exam and patient's structure, with guarantee of minimum dose and adequate quality image. Methods: A patient's absorbed dose equation was developed by means of Entrance Skin Dose (ESD), irradiated area and patient width for specific anatomy. ESD is calculated from a developed equation, where entrance surface air-KERMA and backscatter factor are included, with air-to-skin coefficient conversion. We developed specific Lambert-Beer attenuation equations derived from mass energy-absorption coefficients data for skin, fat, and muscle and bone as one tissue. Anatomy tissue thickness distribution at central X-ray location in anteroposterior incidence for hand and chest, was estimate by discounting constant skin and bone thickness from patient measured width, assuming the result as muscle and fat. A clinical research at a big hospital were executed when real parameters (kV, mAs, filtration, ripple) used by technologists were combined with the image quality and patient's data: anatomy width, height and weight. A correlation among the best images acquired and electrical parameters used were confronted with patient's data and dose estimation. The best combinations were used as gold standards. Results: For each anatomy, two equations were developed to calculate voltage (kV) and exposure (mAs) to reproduce and interpolate the gold standards. Patient is measured and data are input into equations, giving radiological technologists the right set of electrical parameters for that specific exam. Conclusion: This work indicates that radiological technologist can personalize the exact electrical parameters for each patient exam, instead of using standard values. It also guarantee that patients under or over-sized measures will receive the right dose for the best image. It will stop wrong empiric adjusts technologists do when examining a non-standard patient and reduce probability of radiography retaken because of over or under exposition.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-08: A Hybrid Three Dimensional Registration Framework for Image-Guided Accurate Radiotherapy System ARTS-IGRT

Purpose: The purpose of this work was to develop a registration framework and method based on the software platform of ARTS-IGRT and implement in C++ based on ITK libraries to register CT images and CBCT images. ARTS-IGRT was a part of our self-developed accurate radiation planning system ARTS. Methods: Mutual information (MI) registration treated each voxel equally. Actually, different voxels even having same intensity should be treated differently in the registration procedure. According to their importance values calculated from self-information, a similarity measure was proposed which combined the spatial importance of a voxel with MI (S-MI). For lung registration, Firstly, a global alignment method was adopted to minimize the margin error and achieve the alignment of these two images on the whole. The result obtained at the low resolution level was then interpolated to become the initial conditions for the higher resolution computation. Secondly, a new similarity measurement S-MI was established to quantify how close the two input image volumes were to each other. Finally, Demons model was applied to compute the deformable map. Results: Registration tools were tested for head-neck and lung images and the average region was 128*128*49. The rigid registration took approximately 2 min and converged 10% faster than traditional MI algorithm, the accuracy reached 1mm for head-neck images. For lung images, the improved symmetric Demons registration process was completed in an average of 5 min using a 2.4GHz dual core CPU. Conclusion: A registration framework was developed to correct patient's setup according to register the planning CT volume data and the daily reconstructed 3D CBCT data. The experiments showed that the spatial MI algorithm can be adopted for head-neck images. The improved Demons deformable registration was more suitable to lung images, and rigid alignment should be applied before deformable registration to get more accurate result. Supported by National Natural Science Foundation of China (NO.81101132) and Natural Science Foundation of Anhui Province (NO.11040606Q55)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-10: A Moving-Blocker-Based Strategy for Simultaneous Megavoltage and Kilovoltage Scatter Correction in Cone-Beam Computed Tomography Image Acquired During Volumetric Modulated Arc Therapy

Purpose: To evaluate a moving-blocker-based approach in estimating and correcting megavoltage (MV) and kilovoltage (kV) scatter contamination in kV cone-beam computed tomography (CBCT) acquired during volumetric modulated arc therapy (VMAT). Methods: XML code was generated to enable concurrent CBCT acquisition and VMAT delivery in Varian TrueBeam developer mode. A physical attenuator (i.e., “blocker”) consisting of equal spaced lead strips (3.2mm strip width and 3.2mm gap in between) was mounted between the x-ray source and patient at a source to blocker distance of 232mm. The blocker was simulated to be moving back and forth along the gantry rotation axis during the CBCT acquisition. Both MV and kV scatter signal were estimated simultaneously from the blocked regions of the imaging panel, and interpolated into the un-blocked regions. Scatter corrected CBCT was then reconstructed from un-blocked projections after scatter subtraction using an iterative image reconstruction algorithm based on constraint optimization. Experimental studies were performed on a Catphan 600 phantom and an anthropomorphic pelvis phantom to demonstrate the feasibility of using moving blocker for MV-kV scatter correction. Results: MV scatter greatly degrades the CBCT image quality by increasing the CT number inaccuracy and decreasing the image contrast, in addition to the shading artifacts caused by kV scatter. The artifacts were substantially reduced in the moving blocker corrected CBCT images in both Catphan and pelvis phantoms. Quantitatively, CT number error in selected regions of interest reduced from 377 in the kV-MV contaminated CBCT image to 38 for the Catphan phantom. Conclusions: The moving-blockerbased strategy can successfully correct MV and kV scatter simultaneously in CBCT projection data acquired with concurrent VMAT delivery. This work was supported in part by a grant from the Cancer Prevention and Research Institute of Texas (RP130109) and a grant from the American Cancer Society (RSG-13-326-01-CCE)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-14: A Novel Approach to Evaluate the Dosimetric Effect of Rectal Variation During Image Guided Prostate Radiotherapy

Purpose: Deformable registration establishes the spatial correspondence back to the reference image in order to accumulate dose. However, in prostate radiotherapy the changing shape and volume of the rectum present a challenge to accurate deformable registration and consequently calculation of delivered dose. We explored an alternative approach to calculating accumulated dose to the rectum, independent of deformable registration. Methods: This study was performed on three patients who received online image-guided radiotherapy (IGRT) with daily CBCT (XVI-system,Elekta) and target localization using intraprostatic fiducials. On each CBCT, the rectum was manually contoured and bulk density assignments were made allowing dose to be calculated for each fraction. Dose-surface maps (DSM) were generated (MATLAB,Mathworks,Natick,MA) by considering the rectum as a cylinder and sampling the dose at 21-equispaced points on each CT slice. The cylinder was “cut” at the posterior-most position on each CT and unfolded to generate a DSM. These were normalised in the longitudinal direction by interpolation creating maps of 21×21 pixels. A DSM was produced for each CBCT and the dose was accumulated. Results: The mean accumulated delivered rectal surface dose was on average 7.5(+/−3.5)% lower than the planned dose. The dose difference maps consistently show that the greatest variation in dose between planned and delivered dose is away from where the rectal surface is adjacent to the prostate. Conclusion: Estimation of dose accumulation using DSM provides an alternative method for determining actual delivered dose to the rectum. The dose difference is greatest in areas away from the region where the rectal surface abuts the prostate, the region where set-up is verified. The change in size and shape of the rectum was shown to resultin a change in the accumulated dose compared to the planned dose and this will have an impact on determining the relationships between dose delivered and toxicity. We acknowledge funding from CRUK and acknowledge NHS funding to the NIHR Biomedical Research Centre for Cancer. Patients were treated within the CHHiP IGRT sub-study (CRUK/06/016, ISRCTN:97182923) funded by CRUK. RayStation was used under an evaluation agreement with RaySearch Laboratories AB.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-246: MU Model Implementation for a Passively Scattered Proton Beam: Inclusion of the Bolus Gap and Nozzle Extension Factors

Purpose: To develop and implement an MU model for a passively scattered proton beam, and eliminate the need for patient specific calibrations for field sizes from 3cm to 15cm. This would enable consistent and timely calibrations for a wide variety of patient portals, streamlining the treatment process. Methods: Measurements were initially made using a Standard Imaging A-16 ion chamber and a modified water tank to determine the bolus gap factors (BGF) for multiple combinations of aperture size, bolus thickness, and air gap. The BGF was then separated into two component factors: the bolus thickness factor (BTF) and the nozzle extension factor (NEF). Polynomial curves were generated using the measured data to produce BTF tables for air gaps from 0cm to 30cm and for bolus thicknesses from 0cm to 10cm, and NEF tables for the full range of clinically used nozzle extensions. Additionally, data tables were created for every factor that affects beam output in the MU model. The MUs were then modeled for 487 patient portals and retrospectively compared to the MUs generated from the physical calibrations previously performed. Results: Of the 487 patient portals tested, 100% of the portals used for the comparison were within 2.5% from the MUs generated using a physical calibration, and 95.9% of the MU model portals tested were within 2%. The patient portals tested had field sizes ranging from 2.1cm to 10.1cm, with air gaps from 2cm to 25cm. Output factors for field sizes below 3cm with irregularly shaped fields demonstrated inconsistent results and will be further studied. Conclusion: The most problematic output factor, the BGF, was modeled accurately and consistently using the lowest order polynomial curve fits and interpolation between measured data. The study results demonstrate the robustness of the MU model and the potential for saving valuable personnel and beam time.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-348: Verification MU Calculation for Conformal Radiotherapy with Multileaf Collimator Using Report AAPM TG 114

Purpose: The purpose of this work was to develop and validate an open source independent MU dose calculation software for 3D conformal radiotherapy with multileaf high and low resolution according to the report of AAPM TG 11 Methods: Treatment plans were done using Iplan v4.5 BrainLAB TPS. A 6MV photon beam produced by Primus and Novalis linear accelerators equipped with an Optifocus MLC and HDMLC, respectively. TPS dose calculation algorithms were pencil beam and Monte Carlo. 1082 treatments plans were selected for the study. The algorithm was written in free and open source CodeBlocks C++ platform. Treatment plans were imported by the software using RTP format. Equivalent size field is obtained from the positions of the leaves; the effective depth of calculation can be introduced by TPS's dosimetry report or automatically calculated starting from SSD. The inverse square law is calculated by the 3D coordinates of the isocenter and normalization point of the treatment plan. The dosimetric parameters TPR, Sc, Sp and WF are linearly interpolated. Results: 1082 plans of both machines were analyzed. The average uncertainty between the TPS and the independent calculation was −0.43% ± 2.42% [−7.90%, 7.50%]. Specifically for the Primus the variation obtained was −0.85% ± 2.53% and for the Novalis 0.00% ± 2.23%. Data show that 94.8% of the cases the uncertainty was less than or equal to 5%, while 98.9% is less than or equal to 6%. Conclusion: The developed software is appropriate for use in calculation of UM. This software can be obtained upon request.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-397: Include Organ Deformation Into Dose Calculation of Prostate Brachytherapy

Purpose: Prostate brachytherapy is an important curative treatment for patients with localized prostate cancer. In brachytherapy, rectal balloon is generally needed to adjust for unfavorable prostate position for seed placement. However, rectal balloon causes prostate deformation, which is not accounted for in dosimetric planning. Therefore, it is possible that brachytherapy dosimetry deviates significantly from initial plan when prostate returns to its non-deformed state (after procedure). The goal of this study is to develop a method to include prostate deformation into the treatment planning of brachytherapy dosimetry. Methods: We prospectively collected ultrasound images of prostate pre- and post- rectal balloon inflation from thirty five consecutive patients undergoing I-125 brachytherapy. Based on the cylinder coordinate systems, we learned the initial coordinate transformation parameters between the manual segmentations of both deformed and non-deformed prostates of each patient in training set. With the nearest-neighbor interpolation, we searched the best transformation between two coordinate systems to maximum the mutual information of deformed and non-deformed images. We then mapped the implanted seeds of five selected patients from the deformed prostate into non-deformed prostate. The seed position is marked on original pre-inflation US image and it is imported into VariSeed software for dose calculation. Results: The accuracy of image registration is 87.5% as quantified by Dice Index. The prostate coverage V100% dropped from 96.5±0.5% of prostate deformed plan to 91.9±2.6% (p<0.05) of non-deformed plan. The rectum V100% decreased from 0.44±0.26 cc to 0.10±0.18 cc (p<0.05). The dosimetry of the urethra showed mild change but not significant: V150% changed from 0.05±0.10 cc to 0.14±0.15 cc (p>0.05) and D1% changed from 212.9±37.3 Gy to 248.4±42.8 Gy (p>0.05). Conclusion: We have developed a deformable image registration method that allows dosimetry calculation in the non-deformed prostate – which may be used to more accurately reflect “true” dosimetry during planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-458: Determining Threshold-Of-Failure for Dead Pixel Rows in EPID-Based Dosimetry

Purpose: A pixel correction map is applied to all EPID-based applications on the TrueBeam (Varian Medical Systems, Palo Alto, CA). When dead pixels are detected, an interpolative smoothing algorithm is applied using neighboring-pixel information to supplement missing-pixel information. The vendor suggests that when the number of dead pixels exceeds 70,000, the panel should be replaced. It is common for entire detector rows to be dead, as well as their neighboring rows. Approximately 70 rows can be dead before the panel reaches this threshold. This study determines the number of neighboring dead-pixel rows that would create a large enough deviation in measured fluence to cause failures in portal dosimetry (PD). Methods: Four clinical two-arc VMAT plans were generated using Eclipse's AXB algorithm and PD plans were created using the PDIP algorithm. These plans were chosen to represent those commonly encountered in the clinic: prostate, lung, abdomen, and neck treatments. During each iteration of this study, an increasing number of dead-pixel rows are artificially applied to the correction map and a fluence QA is performed using the EPID (corrected with this map). To provide a worst-case-scenario, the dead-pixel rows are chosen so that they present artifacts in the highfluence region of the field. Results: For all eight arc-fields deemed acceptable via a 3%/3mm gamma analysis (pass rate greater than 99%), VMAT QA yielded identical results with a 5 pixel-width dead zone. When 10 dead lines were present, half of the fields had pass rates below the 99% pass rate. With increasing dead rows, the pass rates were reduced substantially. Conclusion: While the vendor still suggests to request service at the point where 70,000 dead rows are measured (as recommended by the vendor), the authors suggest that service should be requested when there are greater than 5 consecutive dead rows.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-641: Development and Verification of Automatic Reading Dose of Interest From Eclipse's DVH

Purpose: According to clinical and research requirement, we develop a function of automatic reading dose of interest from dose volume histogram(DVH), to replace the traditional method with a mouse one by one point, and it's also verified. Methods: The DVH automatic reading function will be developed in an in-house developed radiotherapy information management system(RTIMS), which is based on Apache+PHP+MySQL. A DVH ASCII file is exported from Varian Eclipse V8.6, which includes the following contents: 1. basic information of patient; 2. dose information of plan; 3. dose information of structures, including basic information and dose volume data of target volume and organ at risk. And the default exported dose volume data also includes relative doses by 1% step and corresponding absolute doses and cumulative relative volumes, and the volumes are 4 decimal fraction. Clinically, we often need read the doses of some integer percent volumes, such as D50 and D30. So it couldn't be directly obtained from the above data, but we can use linear interpolation bye the near volumes and doses: Dx=D2−(V2−Vx)*(D2−D1)/(V2−V1), and program a function to search, read and calculate the corresponding data. And the doses of all preseted volume of interest of all structures can be automatically read one by one patient, and saved as a CSV file. To verify it, we select 24 IMRT plans for prostate cancer, and doses of interest are PTV D98/D95/D5/D2, bladder D30/D50, and rectum D25/D50. Two groups of data, using the automatic reading method(ARM) and pointed dose method(PDM), are analyzed with SPSS 16. The absolute difference=D-ARM-D-PDM, relative difference=absolute difference*100%/prescription dose(7600cGy). Results: The differences are as following: PTV D98/D95/D5/D2: −0.04%/− 0.04%/0.13%/0.19%, bladder D30/D50: −0.02%/0.01%, and rectum D25/D50: 0.03%/0.01%. Conclusion: Using this function, the error is very small, and can be neglected. It could greatly improve the efficiency of clinical work. Project supported by the National Natural Science Foundation of China (Grant No.81101694)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-F-19A-02: Comparison of Absorbed Dose to Water Standards for HDR Ir-192 Brachytherapy Between the LCR, Brazil and NRC, Canada

Purpose: To compare absorbed dose to water standards for HDR brachytherapy dosimetry developed by the Radiological Science Laboratory of Rio de Janeiro State University (LCR) and the National Research Council, Canada (NRC). Methods: The two institutions have separately developed absorbed dose standards based on the Fricke dosimetry system. There are important differences between the two standards, including: preparation and read-out of the Fricke solution, irradiation geometry of the Fricke holder in relation to the Ir-192 source, and determination of the G-value to be used at Ir-192 energies. All measurements for both standards were made directly at the NRC laboratory (i.e., no transfer instrument was used) using a single Ir-192 source (microSelectron v2). In addition, the NRC group has established a self-consistent method to determine the G-value for Ir-192, based on an interpolation between G-values obtained at Co-60 and 250kVp X-rays, and this measurement was repeated using the LCR Fricke solution to investigate possible systematic uncertainties. Results: G-values for Co-60 and 250 kVp x-rays, obtained using the LCR Fricke system, agreed with the NRC values within 0.5 % and 1 % respectively, indicating that the general assumption of universal G-values is appropriate in this case. The standard uncertainty in the determination of G for Ir-192 is estimated to be 0.6 %. For the comparison of absorbed dose measurements at the reference point for Ir-192 (1 cm depth in water, perpendicular to the seed long-axis), the ratio Dw(NRC)/Dw(LCR) was found to be 1.011 with a combined standard uncertainty of 1.7 %, k=1. Conclusion: The agreement in the absorbed dose to water values for the LCR and NRC systems is very encouraging. Combined with the lower uncertainty in this approach compared to the present air-kerma approach, these results reaffirm the use of Fricke solution as a potential primary standard for HDR Ir-192 brachytherapy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

MO-F-16A-01: Implementation of MPPG TPS Verification Tests On Various Accelerators

Purpose: To demonstrate the implementation of the Medical Physics Practice Guideline (MPPG) for dose calculation and beam parameters verification of treatment planning systems (TPS). Methods: We implemented the draft TPS MPPG for three linacs: Varian Trilogy, TomoHDA and Elekta Infinity. Static and modulated test plans were created. The static fields are different than used in commissioning. Data was collected using ion chambers and diodes in a scanning water tank, Delta4 phantom and a custom phantom. MatLab and Microsoft Excel were used to create analysis tools to compare reference DICOM dose with scan data. This custom code allowed for the interpolation, registration and gamma analysis of arbitrary dose profiles. It will be provided as open source code. IMRT fields were validated with Delta4 registration and comparison tools. The time for each task was recorded. Results: The tests confirmed the strengths, and revealed some limitations, of our TPS. The agreement between calculated and measured dose was reported for all beams. For static fields, percent depth dose and profiles were analyzed with criteria in the draft MPPG. The results reveal areas of slight mismatch with the model (MLC leaf penumbra, buildup region.) For TomoTherapy, the IMRT plan 2%/2 mm gamma analysis revealed poorest agreement in the low dose regions. For one static test plan for all 10MV Trilogy photon beams, the plan generation, scan queue creation, data collection, data analysis and report took 2 hours, excluding tank setup. Conclusions: We have demonstrated the implementation feasibility of the TPS MPPG. This exercise generated an open source tool for dose comparisons between scan data and DICOM dose data. An easily reproducible and efficient infrastructure with streamlined data collection was created for repeatable robust testing of the TPS. The tests revealed minor discrepancies in our models and areas for improvement that are being investigated.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

TU-A-9A-09: Proton Beam X-Ray Fluorescence CT

Purpose: To evaluate x-ray fluorescence computed tomography induced with proton beams (pXFCT) for imaging of gold contrast agent. Methods: Proton-induced x-ray fluorescence was studied by means of Monte Carlo (MC) simulations using TOPAS, a MC code based on GEANT4. First, proton-induced K-shell and L-shell fluorescence was studied as a function of proton beam energy and 1) depth in water and 2) size of contrast object. Second, pXFCT images of a 2-cm diameter cylindrical phantom with four 5- mm diameter contrast vials and of a 20-cm diameter phantom with 1-cm diameter vials were simulated. Contrast vials were filled with water and water solutions with 1-5% gold per weight. Proton beam energies were varied from 70-250MeV. pXFCT sinograms were generated based on the net number of gold K-shell or L-shell x-rays determined by interpolations from the neighboring 0.5keV energy bins of spectra collected with an idealized 4π detector. pXFCT images were reconstructed with filtered-back projection, and no attenuation correction was applied. Results: Proton induced x-ray fluorescence spectra showed very low background compared to x-ray induced fluorescence. Proton induced L-shell fluorescence had a higher cross-section compared to K-shell fluorescence. Excitation of L-shell fluorescence was most efficient for low-energy protons, i.e. at the Bragg peak. K-shell fluorescence increased with increasing proton beam energy and object size. The 2% and 5% gold contrast vials were accurately reconstructed in K-shell pXFCT images of both the 2-cm and 20-cm diameter phantoms. Small phantom L-shell pXFCT image required attenuation correction and had a higher sensitivity for 70MeV protons compared to 250MeV protons. With attenuation correction, L-shell pXFCT might be a feasible option for imaging of small size (∼2cm) objects. Imaging doses for all simulations were 5-30cGy. Conclusion: Proton induced x-ray fluorescence CT promises to be an alternative quantitative imaging technique to the commonly considered XFCT imaging with x-ray beams.

60 APPLIED LIFE SCIENCES↗