Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Radiological Data”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-580: On the Significance of Model Based Dosimetry for Breast and Head and Neck 192Ir HDR Brachytherapy

Purpose: To assess the dosimetric and radiobiological differences between TG43-based and model-based dosimetry in the treatment planning of {sup 192}Ir HDR brachytherapy for breast and head and neck cancer. Methods: Two cohorts of 57 Accelerated Partial Breast Irradiation (APBI) and 22 head and neck (H&N) patients with oral cavity carcinoma were studied. Dosimetry for the treatment plans was performed using the TG43 algorithm of the Oncentra Brachy v4.4 treatment planning system (TPS). Corresponding Monte Carlo (MC) simulations were performed using MCNP6 with input files automatically prepared by the BrachyGuide software tool from DICOM RT plan data. TG43 and MC data were compared in terms of % dose differences, Dose Volume Histograms (DVHs) and related indices of clinical interest for the Planning Target Volume (PTV) and the Organs-At-Risk (OARs). A radiobiological analysis was also performed using the Equivalent Uniform Dose (EUD), mean survival fraction (S) and Tumor Control Probability (TCP) for the PTV, and the Normal Tissue Control Probability (N TCP) and the generalized EUD (gEUD) for the OARs. Significance testing of the observed differences performed using the Wilcoxon paired sample test. Results: Differences between TG43 and MC DVH indices, associated with the increased corresponding local % dose differences observed, were statistically significant. This is mainly attributed to their consistency however, since TG43 agrees closely with MC for the majority of DVH and radiobiological parameters in both patient cohorts. Differences varied considerably among patients only for the ipsilateral lung and ribs in the APBI cohort, with a strong correlation to target location. Conclusion: While the consistency and magnitude of differences in the majority of clinically relevant DVH indices imply that no change is needed in the treatment planning practice, individualized dosimetry improves accuracy and addresses instances of inter-patient variability observed. Research co-financed by the ESF and Greek funds through the Operational Program Education and Lifelong Learning Investing in Knowledge Society of the NSRF. Research Funding Program Aristeia. Nucletron, an Elekta company (Veenendaal, The Netherlands) is gratefully acknowledged for providing Oncentra Brachy v4.4 for research purposes.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Using compressive sensing to recover images from PET scanners with partial detector rings

Purpose: Most positron emission tomography/computed tomography (PET/CT) scanners consist of tightly packed discrete detector rings to improve scanner efficiency. The authors’ aim was to use compressive sensing (CS) techniques in PET imaging to investigate the possibility of decreasing the number of detector elements per ring (introducing gaps) while maintaining image quality. Methods: A CS model based on a combination of gradient magnitude and wavelet domains (wavelet-TV) was developed to recover missing observations in PET data acquisition. The model was designed to minimize the total variation (TV) and L1-norm of wavelet coefficients while constrained by the partially observed data. The CS model also incorporated a Poisson noise term that modeled the observed noise while suppressing its contribution by penalizing the Poisson log likelihood function. Three experiments were performed to evaluate the proposed CS recovery algorithm: a simulation study, a phantom study, and six patient studies. The simulation dataset comprised six disks of various sizes in a uniform background with an activity concentration of 5:1. The simulated image was multiplied by the system matrix to obtain the corresponding sinogram and then Poisson noise was added. The resultant sinogram was masked to create the effect of partial detector removal and then the proposed CS algorithm was applied to recover the missing PET data. In addition, different levels of noise were simulated to assess the performance of the proposed algorithm. For the phantom study, an IEC phantom with six internal spheres each filled with F-18 at an activity-to-background ratio of 10:1 was used. The phantom was imaged twice on a RX PET/CT scanner: once with all detectors operational (baseline) and once with four detector blocks (11%) turned off at each of 0 °, 90 °, 180 °, and 270° (partially sampled). The partially acquired sinograms were then recovered using the proposed algorithm. For the third test, PET images from six patient studies were investigated using the same strategy of the phantom study. The recovered images using WTV and TV as well as the partially sampled images from all three experiments were then compared with the fully sampled images (the baseline). Comparisons were done by calculating the mean error (%bias), root mean square error (RMSE), contrast recovery (CR), and SNR of activity concentration in regions of interest drawn in the background as well as the disks, spheres, and lesions. Results: For the simulation study, the mean error, RMSE, and CR for the WTV (TV) recovered images were 0.26% (0.48%), 2.6% (2.9%), 97% (96%), respectively, when compared to baseline. For the partially sampled images, these results were 22.5%, 45.9%, and 64%, respectively. For the simulation study, the average SNR for the baseline was 41.7 while for WTV (TV), recovered image was 44.2 (44.0). The phantom study showed similar trends with 5.4% (18.2%), 15.6% (18.8%), and 78% (60%), respectively, for the WTV (TV) images and 33%, 34.3%, and 69% for the partially sampled images. For the phantom study, the average SNR for the baseline was 14.7 while for WTV (TV) recovered image was 13.7 (11.9). Finally, the average of these values for the six patient studies for the WTV-recovered, TV, and partially sampled images was 1%, 7.2%, 92% and 1.3%, 15.1%, 87%, and 27%, 25.8%, 45%, respectively. Conclusions: CS with WTV is capable of recovering PET images with good quantitative accuracy from partially sampled data. Such an approach can be used to potentially reduce the cost of scanners while maintaining good image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

A Study of the Effects of Different Urban Wind Models on Dispersion Patterns Using Joint Urban 2003 Data

The Quick Urban & Industrial Complex (QUIC) Dispersion Modeling System has been developed to rapidly compute the transport and dispersion of toxic agent releases in the vicinity of buildings. It is composed of a wind solver, an 'urbanized' Lagrangian random-walk model, and a graphical user interface. QUIC has two different wind models: (a) The QUIC-URB wind solver, an empirically-based diagnostic wind model and (b) The QUIC-CFD (RANS) solver, based on the 3D Reynolds-Averaged Navier-Stokes (RANS) equations. In this paper, we discuss the effect of different wind models on dispersion patterns in dense built-up areas. The model-computed wind from the two urban wind models- QUIC-URB and QUIC-CFD are used to drive the dispersion model. The concentration fields are then compared to measurements from the Oklahoma City Joint Urban 2003 field experiment. QUIC produces high-resolution 3-D mean wind and concentration fields around buildings, in addition to deposition on the ground and building surfaces. It has options for different release types, including point, moving point, line, area, and volumetric sources, as well as dense gas, explosive buoyant rise, multi-particle size, bioslurry, and two-phase releases. Other features include indoor infiltration, a pressure solver, outer grid simulations, vegetative canopies, and population exposure calculations. It has been used for biological agent sensor siting in cities, vulnerability assessments for heavier-than-air chemical releases at industrial facilities, and clean-up assessments for radiological dispersal device (RDD) releases in cities (e.g., see Linger et al., 2005; Brown, 2006a, b). QUIC has also been used for dust transport studies (Bowker et al., 2007a) and for the impact of highway sound barriers on the transport and dispersion of vehicle emissions (Bowker et al., 2007b).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

SU-E-T-558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry

Purpose: Monte Carlo simulation on GPU has experienced rapid advancements over the past a few years and tremendous accelerations have been achieved. Yet existing packages were developed only in voxelized geometry. In some applications, e.g. radioactive seed modeling, simulations in more complicated geometry are needed. This abstract reports our initial efforts towards developing a quadric geometry module aiming at expanding the application scope of GPU-based MC simulations. Methods: We defined the simulation geometry consisting of a number of homogeneous bodies, each specified by its material composition and limiting surfaces characterized by quadric functions. A tree data structure was utilized to define geometric relationship between different bodies. We modified our GPU-based photon MC transport package to incorporate this geometry. Specifically, geometry parameters were loaded into GPU’s shared memory for fast access. Geometry functions were rewritten to enable the identification of the body that contains the current particle location via a fast searching algorithm based on the tree data structure. Results: We tested our package in an example problem of HDR-brachytherapy dose calculation for shielded cylinder. The dose under the quadric geometry and that under the voxelized geometry agreed in 94.2% of total voxels within 20% isodose line based on a statistical t-test (95% confidence level), where the reference dose was defined to be the one at 0.5cm away from the cylinder surface. It took 243sec to transport 100million source photons under this quadric geometry on an NVidia Titan GPU card. Compared with simulation time of 99.6sec in the voxelized geometry, including quadric geometry reduced efficiency due to the complicated geometry-related computations. Conclusion: Our GPU-based MC package has been extended to support photon transport simulation in quadric geometry. Satisfactory accuracy was observed with a reduced efficiency. Developments for charged particle transport in this geometry are currently in progress.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Comparison of flat cleaved and cylindrical diffusing fibers as treatment sources for interstitial photodynamic therapy

Purpose: For interstitial photodynamic therapy (iPDT) of bulky tumors, careful treatment planning is required in order to ensure that a therapeutic dose is delivered to the tumor, while minimizing damage to surrounding normal tissue. In clinical contexts, iPDT has typically been performed with either flat cleaved or cylindrical diffusing optical fibers as light sources. Here, the authors directly compare these two source geometries in terms of the number of fibers and duration of treatment required to deliver a prescribed light dose to a tumor volume. Methods: Treatment planning software for iPDT was developed based on graphics processing unit enhanced Monte Carlo simulations. This software was used to optimize the number of fibers, total energy delivered by each fiber, and the position of individual fibers in order to deliver a target light dose (D{sub 90}) to 90% of the tumor volume. Treatment plans were developed using both flat cleaved and cylindrical diffusing fibers, based on tissue volumes derived from CT data from a head and neck cancer patient. Plans were created for four cases: fixed energy per fiber, fixed number of fibers, and in cases where both or neither of these factors were fixed. Results: When the number of source fibers was fixed at eight, treatment plans based on flat cleaved fibers required each to deliver 7180–8080 J in order to deposit 90 J/cm{sup 2} in 90% of the tumor volume. For diffusers, each fiber was required to deliver 2270–2350 J (333–1178 J/cm) in order to achieve this same result. For the case of fibers delivering a fixed 900 J, 13 diffusers or 19 flat cleaved fibers at a spacing of 1 cm were required to deliver the desired dose. With energy per fiber fixed at 2400 J and the number of fibers fixed at eight, diffuser fibers delivered the desired dose to 93% of the tumor volume, while flat cleaved fibers delivered this dose to 79%. With both energy and number of fibers allowed to vary, six diffusers delivering 3485–3600 J were required, compared to ten flat cleaved fibers delivering 2780–3600 J. Conclusions: For the same number of fibers, cylindrical diffusers allow for a shorter treatment duration compared to flat cleaved fibers. For the same energy delivered per fiber, diffusers allow for the insertion of fewer fibers in order to deliver the same light dose to a target volume.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Monte Carlo simulations of photoneutron activation in a medical LINAC

The first step of the present study was to simulate the γ-radioactivity produced by photoneutrons in the treatment head of a decommissioned medical LINAC. These results, obtained with MCNPX and CINDER'90, are then compared with γ-ray spectrometric measurements of parts of the treatment head. Our ultimate objective is to simulate the associated production by photonuclear reactions of pure β-emitters, which are more difficult to access experimentally than the γ-ray emitters. (authors)

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Use of Monte Carlo simulations with a realistic rat phantom for examining the correlation between hematopoietic system response and red marrow absorbed dose in Brown Norway rats undergoing radionuclide therapy with {sup 177}Lu- and {sup 90}Y-BR96 mAbs

Purpose: Biokinetic and dosimetry studies in laboratory animals often precede clinical radionuclide therapies in humans. A reliable evaluation of therapeutic efficacy is essential and should be based on accurate dosimetry data from a realistic dosimetry model. The aim of this study was to develop an anatomically realistic dosimetry model for Brown Norway rats to calculate S factors for use in evaluating correlations between absorbed dose and biological effects in a preclinical therapy study. Methods: A realistic rat phantom (Roby) was used, which has some flexibility that allows for a redefinition of organ sizes. The phantom was modified to represent the anatomic geometry of a Brown Norway rat, which was used for Monte Carlo calculations of S factors. Kinetic data for radiolabeled BR96 monoclonal antibodies were used to calculate the absorbed dose. Biological data were gathered from an activity escalation study with {sup 90}Y- and {sup 177}Lu-labeled BR96 monoclonal antibodies, in which blood cell counts and bodyweight were examined up to 2 months follow-up after injection. Reductions in white blood cell and platelet counts and declines in bodyweight were quantified by four methods and compared to the calculated absorbed dose to the bone marrow or the total body. Results: A red marrow absorbed dose-dependent effect on hematological parameters was observed, which could be evaluated by a decrease in blood cell counts. The absorbed dose to the bone marrow, corresponding to the maximal tolerable activity that could safely be administered, was determined to 8.3 Gy for {sup 177}Lu and 12.5 Gy for {sup 90}Y. Conclusions: There was a clear correlation between the hematological effects, quantified with some of the studied parameters, and the calculated red marrow absorbed doses. The decline in body weight was stronger correlated to the total body absorbed dose, rather than the red marrow absorbed dose. Finally, when considering a constant activity concentration, the phantom weight, ranging from 225 g to 300 g, appeared to have no substantial effect for the estimated absorbed dose.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Addendum to the AAPM's TG-51 protocol for clinical reference dosimetry of high-energy photon beams

An addendum to the AAPM's TG-51 protocol for the determination of absorbed dose to water in megavoltage photon beams is presented. This addendum continues the procedure laid out in TG-51 but new k{sub Q} data for photon beams, based on Monte Carlo simulations, are presented and recommendations are given to improve the accuracy and consistency of the protocol's implementation. The components of the uncertainty budget in determining absorbed dose to water at the reference point are introduced and the magnitude of each component discussed. Finally, the consistency of experimental determination of N{sub D,w} coefficients is discussed. It is expected that the implementation of this addendum will be straightforward, assuming that the user is already familiar with TG-51. The changes introduced by this report are generally minor, although new recommendations could result in procedural changes for individual users. It is expected that the effort on the medical physicist's part to implement this addendum will not be significant and could be done as part of the annual linac calibration.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-I-20: Dead Time Count Loss Compensation in SPECT/CT: Projection Versus Global Correction

Purpose: To compare projection-based versus global correction that compensate for deadtime count loss in SPECT/CT images. Methods: SPECT/CT images of an IEC phantom (2.3GBq 99mTc) with ∼10% deadtime loss containing the 37mm (uptake 3), 28 and 22mm (uptake 6) spheres were acquired using a 2 detector SPECT/CT system with 64 projections/detector and 15 s/projection. The deadtime, Ti and the true count rate, Ni at each projection, i was calculated using the monitor-source method. Deadtime corrected SPECT were reconstructed twice: (1) with projections that were individually-corrected for deadtime-losses; and (2) with original projections with losses and then correcting the reconstructed SPECT images using a scaling factor equal to the inverse of the average fractional loss for 5 projections/detector. For both cases, the SPECT images were reconstructed using OSEM with attenuation and scatter corrections. The two SPECT datasets were assessed by comparing line profiles in xyplane and z-axis, evaluating the count recoveries, and comparing ROI statistics. Higher deadtime losses (up to 50%) were also simulated to the individually corrected projections by multiplying each projection i by exp(-a*Ni*Ti), where a is a scalar. Additionally, deadtime corrections in phantoms with different geometries and deadtime losses were also explored. The same two correction methods were carried for all these data sets. Results: Averaging the deadtime losses in 5 projections/detector suffices to recover >99% of the loss counts in most clinical cases. The line profiles (xyplane and z-axis) and the statistics in the ROIs drawn in the SPECT images corrected using both methods showed agreement within the statistical noise. The count-loss recoveries in the two methods also agree within >99%. Conclusion: The projection-based and the global correction yield visually indistinguishable SPECT images. The global correction based on sparse sampling of projections losses allows for accurate SPECT deadtime loss correction while keeping the study duration reasonable.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-253: The Radiomics Toolbox in the Computational Environment for Radiological Research (CERR)

Purpose: To present an open source and free platform to facilitate radiomics research — The “Radiomics toolbox” in CERR. Method: There is scarcity of open source tools that support end-to-end modeling of image features to predict patient outcomes. The “Radiomics toolbox” strives to fill the need for such a software platform. The platform supports (1) import of various kinds of image modalities like CT, PET, MR, SPECT, US. (2) Contouring tools to delineate structures of interest. (3) Extraction and storage of image based features like 1st order statistics, gray-scale co-occurrence and zonesize matrix based texture features and shape features and (4) Statistical Analysis. Statistical analysis of the extracted features is supported with basic functionality that includes univariate correlations, Kaplan-Meir curves and advanced functionality that includes feature reduction and multivariate modeling. The graphical user interface and the data management are performed with Matlab for the ease of development and readability of code and features for wide audience. Open-source software developed with other programming languages is integrated to enhance various components of this toolbox. For example: Java-based DCM4CHE for import of DICOM, R for statistical analysis. Results: The Radiomics toolbox will be distributed as an open source, GNU copyrighted software. The toolbox was prototyped for modeling Oropharyngeal PET dataset at MSKCC. The analysis will be presented in a separate paper. Conclusion: The Radiomics Toolbox provides an extensible platform for extracting and modeling image features. To emphasize new uses of CERR for radiomics and image-based research, we have changed the name from the “Computational Environment for Radiotherapy Research” to the “Computational Environment for Radiological Research”.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Monte Carlo simulation of breast imaging using synchrotron radiation

Purpose: Synchrotron radiation (SR), being the brightest artificial source of x-rays with a very promising geometry, has raised the scientific expectations that it could be used for breast imaging with optimized results. The ''in situ'' evaluation of this technique is difficult to perform, mostly due to the limited available SR facilities worldwide. In this study, a simulation model for SR breast imaging was developed, based on Monte Carlo simulation techniques, and validated using data acquired in the SYRMEP beamline of the Elettra facility in Trieste, Italy. Furthermore, primary results concerning the performance of SR were derived. Methods: The developed model includes the exact setup of the SR beamline, considering that the x-ray source is located at almost 23 m from the slit, while the photon energy was considered to originate from a very narrow Gaussian spectrum. Breast phantoms, made of Perspex and filled with air cavities, were irradiated with energies in the range of 16-28 keV. The model included a Gd{sub 2}O{sub 2}S detector with the same characteristics as the one available in the SYRMEP beamline. Following the development and validation of the model, experiments were performed in order to evaluate the contrast resolution of SR. A phantom made of adipose tissue and filled with inhomogeneities of several compositions and sizes was designed and utilized to simulate the irradiation under conventional mammography and SR conditions. Results: The validation results of the model showed an excellent agreement with the experimental data, with the correlation for contrast being 0.996. Significant differences only appeared at the edges of the phantom, where phase effects occur. The initial evaluation experiments revealed that SR shows very good performance in terms of the image quality indices utilized, namely subject contrast and contrast to noise ratio. The response of subject contrast to energy is monotonic; however, this does not stand for contrast to noise ratio, since there is a range of optimal performance for SR (18-21 keV). In comparison to conventional mammography, SR shows improved subject contrast for energies lower than the mean energy of each spectrum. Conclusions: The comparison of the results of the two models, conventional and SR, proved that SR exhibits better performance in the majority of cases. The proposed simulation model offers the possibility to perform exhaustive search to evaluate the performance of SR in clinical applications such as breast imaging.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-I-21: Deformation Mapping and Shape Prediction with 3D Tumor Volume Morphing

Purpose: To compare projection-based versus global correction that compensate for deadtime count loss in SPECT/CT images. Methods: SPECT/CT images of an IEC phantom (2.3GBq 99mTc) with ∼10% deadtime loss containing the 37mm (uptake 3), 28 and 22mm (uptake 6) spheres were acquired using a 2 detector SPECT/CT system with 64 projections/detector and 15 s/projection. The deadtime, Ti and the true count rate, Ni at each projection, i was calculated using the monitor-source method. Deadtime corrected SPECT were reconstructed twice: (1) with projections that were individually-corrected for deadtime-losses; and (2) with original projections with losses and then correcting the reconstructed SPECT images using a scaling factor equal to the inverse of the average fractional loss for 5 projections/detector. For both cases, the SPECT images were reconstructed using OSEM with attenuation and scatter corrections. The two SPECT datasets were assessed by comparing line profiles in xyplane and z-axis, evaluating the count recoveries, and comparing ROI statistics. Higher deadtime losses (up to 50%) were also simulated to the individually corrected projections by multiplying each projection i by exp(-a*Ni*Ti), where a is a scalar. Additionally, deadtime corrections in phantoms with different geometries and deadtime losses were also explored. The same two correction methods were carried for all these data sets. Results: Averaging the deadtime losses in 5 projections/detector suffices to recover >99% of the loss counts in most clinical cases. The line profiles (xyplane and z-axis) and the statistics in the ROIs drawn in the SPECT images corrected using both methods showed agreement within the statistical noise. The count-loss recoveries in the two methods also agree within >99%. Conclusion: The projection-based and the global correction yield visually indistinguishable SPECT images. The global correction based on sparse sampling of projections losses allows for accurate SPECT deadtime loss correction while keeping the study duration reasonable.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Use of dMLC for implementation of dynamic respiratory-gated radiation therapy

Purpose: To simulate and evaluate the use of dynamic multileaf collimators (dMLC) in respiratory gating to compensate for baseline drift.Methods: Tumor motion tracking data from 30 lung tumors over 322 treatment fractions was analyzed with the finite state model. A dynamic respiratory gating window was established in real-time by determining the average positions during the previous two end-of-expiration breathing phases and centering the dMLC aperture on a weighted average of these positions. A simulated dMLC with physical motion constraints was used in dynamic gating treatment simulations. Fluence maps were created to provide a statistical description of radiation delivery for each fraction. Duty cycle was also calculated for each fraction.Results: The average duty cycle was 2.3% greater under dynamic gating conditions. Dynamic gating also showed higher fluences and less tumor obstruction. Additionally, dynamic gating required fewer beam toggles and each delivery period was longer on average than with static gating.Conclusions: The use of dynamic gating showed better performance than static gating and the physical constraints of a dMLC were shown to not be an impediment to dynamic gating.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Technical Note: Development of a 3D printed subresolution sandwich phantom for validation of brain SPECT analysis

Purpose: To make an adaptable, head shaped radionuclide phantom to simulate molecular imaging of the brain using clinical acquisition and reconstruction protocols. This will allow the characterization and correction of scanner characteristics, and improve the accuracy of clinical image analysis, including the application of databases of normal subjects. Methods: A fused deposition modeling 3D printer was used to create a head shaped phantom made up of transaxial slabs, derived from a simulated MRI dataset. The attenuation of the printed polylactide (PLA), measured by means of the Hounsfield unit on CT scanning, was set to match that of the brain by adjusting the proportion of plastic filament and air (fill ratio). Transmission measurements were made to verify the attenuation of the printed slabs. The radionuclide distribution within the phantom was created by adding {sup 99m}Tc pertechnetate to the ink cartridge of a paper printer and printing images of gray and white matter anatomy, segmented from the same MRI data. The complete subresolution sandwich phantom was assembled from alternate 3D printed slabs and radioactive paper sheets, and then imaged on a dual headed gamma camera to simulate an HMPAO SPECT scan. Results: Reconstructions of phantom scans successfully used automated ellipse fitting to apply attenuation correction. This removed the variability inherent in manual application of attenuation correction and registration inherent in existing cylindrical phantom designs. The resulting images were assessed visually and by count profiles and found to be similar to those from an existing elliptical PMMA phantom. Conclusions: The authors have demonstrated the ability to create physically realistic HMPAO SPECT simulations using a novel head-shaped 3D printed subresolution sandwich method phantom. The phantom can be used to validate all neurological SPECT imaging applications. A simple modification of the phantom design to use thinner slabs would make it suitable for use in PET.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Sparing Healthy Tissue and Increasing Tumor Dose Using Bayesian Modeling of Geometric Uncertainties for Planning Target Volume Personalization

Objective: To develop a mathematical tool that can update a patient's planning target volume (PTV) partway through a course of radiation therapy to more precisely target the tumor for the remainder of treatment and reduce dose to surrounding healthy tissue. Methods and Materials: Daily on-board imaging was used to collect large datasets of displacements for patients undergoing external beam radiation therapy for solid tumors. Bayesian statistical modeling of these geometric uncertainties was used to optimally trade off between displacement data collected from previously treated patients and the progressively accumulating data from a patient currently partway through treatment, to optimally predict future displacements for that patient. These predictions were used to update the PTV position and margin width for the remainder of treatment, such that the clinical target volume (CTV) was more precisely targeted. Results: Software simulation of dose to CTV and normal tissue for 2 real prostate displacement datasets consisting of 146 and 290 patients treated with a minimum of 30 fractions each showed that re-evaluating the PTV position and margin width after 8 treatment fractions reduced healthy tissue dose by 19% and 17%, respectively, while maintaining CTV dose. Conclusion: Incorporating patient-specific displacement patterns from early in a course of treatment allows PTV adaptation for the remainder of treatment. This substantially reduces the dose to healthy tissues and thus can reduce radiation therapy–induced toxicities, improving patient outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Automated Calibration for Rapid Optical Spectroscopy Sensor Development for Online Monitoring

An automated platform has been developed to assist researchers in the rapid development of optical spectroscopy sensors to quantify species from spectral data. This platform performs calibration and validation measurements simultaneously. Real-time, in situ monitoring of complex systems through optical spectroscopy has been shown to be a useful tool; however, building calibration models requires development time, which can be a limiting factor in the case of radiological or otherwise hazardous systems. While calibration time can be reduced through optimized design of experiments, this study approached the challenge differently through automation. The ATLAS (Automated Transient Learning for Applied Sensors) platform used pneumatic control of stock solutions to cycle flow profiles through desired calibration concentrations for multivariate model construction. Additionally, the transients between desired concentrations based on flow calculations were used as validation measurements to understand model predictive capabilities. This automated approach yielded an incredible 76% reduction in model development time and a 60% reduction in sample volume versus estimated manual sample preparation and static measurements. The ATLAS system was demonstrated on two systems: a three-lanthanide system with Pr/Nd/Ho representing a use case with significant overlap or interference between analyte signatures and an alternate system containing Pr/Nd/Ni to demonstrate a use case in which broad-band corrosion species signatures interfered with more distinct lanthanide absorbance profiles. Both systems resulted in strong model prediction performance (RMSEP < 9%). Lastly, ATLAS was demonstrated as a tool to simulate process monitoring scenarios (e.g., column separation) in which models can be further optimized to account for day-to-day changes as necessary (e.g., baseline correction). Ultimately, ATLAS offers a vital tool to rapidly screen monitoring methods, investigate sensor fusion, and explore more complex systems (i.e., larger numbers of species).

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Robustness of the Voluntary Breath-Hold Approach for the Treatment of Peripheral Lung Tumors Using Hypofractionated Pencil Beam Scanning Proton Therapy

Purpose: The safe clinical implementation of pencil beam scanning (PBS) proton therapy for lung tumors is complicated by the delivery uncertainties caused by breathing motion. The purpose of this feasibility study was to investigate whether a voluntary breath-hold technique could limit the delivery uncertainties resulting from interfractional motion. Methods and Materials: Data from 15 patients with peripheral lung tumors previously treated with stereotactic radiation therapy were included in this study. The patients had 1 computed tomographic (CT) scan in voluntary breath-hold acquired before treatment and 3 scans during the treatment course. PBS proton treatment plans with 2 fields (2F) and 3 fields (3F), respectively, were calculated based on the planning CT scan and subsequently recalculated on the 3 repeated CT scans. Recalculated plans were considered robust if the V{sub 95%} (volume receiving ≥95% of the prescribed dose) of the gross target volume (GTV) was within 5% of what was expected from the planning CT data throughout the simulated treatment. Results: A total of 14/15 simulated treatments for both 2F and 3F met the robustness criteria. Reduced V{sub 95%} was associated with baseline shifts (2F, P=.056; 3F, P=.008) and tumor size (2F, P=.025; 3F, P=.025). Smaller tumors with large baseline shifts were also at risk for reduced V{sub 95%} (interaction term baseline/size: 2F, P=.005; 3F, P=.002). Conclusions: The breath-hold approach is a realistic clinical option for treating lung tumors with PBS proton therapy. Potential risk factors for reduced V{sub 95%} are small targets in combination with large baseline shifts. On the basis of these results, the baseline shift of the tumor should be monitored (eg, through image guided therapy), and appropriate measures should be taken accordingly. The intrafractional motion needs to be investigated to confirm that the breath-hold approach is robust.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

TU-C-17A-08: Improving IMRT Planning and Reducing Inter-Planner Variability Using the Stochastic Frontier Method: Validation Based On Clinical and Simulated Data

Purpose: Intensity modulated radiation therapy always requires compromises between PTV coverage and organs at risk (OAR) sparing. We previously developed metrics that correlate doses to OAR to specific patients’ morphology using stochastic frontier analysis (SFA). Here, we aim to examine the validity of this approach using a large set of realistically simulated dosimetric and geometric data. Methods: SFA describes a set of treatment plans as an asymmetric distribution with respect to a frontier defining optimal plans. Eighty head and neck IMRT plans were used to establish a metric predicting the mean dose to parotids as a function of simple geometric parameters. A database of 140 parotids was used as a basis distribution to simulate physically plausible data of geometry and dose. Distributions comprising between 20 and 5000 were simulated and the SFA was applied to obtain new frontiers, which were compared to the original frontier. Results: It was possible to simulate distributions consistent with the original dataset. Below 160 organs, the SFA could not always describe distributions as asymmetric: a few cases showed a Gaussian or half-Gaussian distribution. In order to converge to a stable solution, the number of organs in a distribution must ideally be above 100, but in many cases stable parameters could be achieved with as low as 60 samples of organ data. Mean RMS value of the error of new frontiers was significantly reduced when additional organs are used. Conclusion: The number of organs in a distribution showed to have an impact on the effectiveness of the model. It is always possible to obtain a frontier, but if the number of organs in the distribution is small (< 160), it may not represent de lowest dose achievable. These results will be used to determine number of cases necessary to adapt the model to other organs.

60 APPLIED LIFE SCIENCES↗