Search NASA⌕ Search

SEARCH · Search NASA

Results for “image simulation”

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,621 records · Page 90

A low dose simulation tool for CT systems with energy integrating detectors

Purpose: This paper introduces a new strategy for simulating low-dose computed tomography (CT) scans using real scans of a higher dose as an input. The tool is verified against simulations and real scans and compared to other approaches found in the literature. Methods: The conditional variance identity is used to properly account for the variance of the input high-dose data, and a formula is derived for generating a new Poisson noise realization which has the same mean and variance as the true low-dose data. The authors also derive a formula for the inclusion of real samples of detector noise, properly scaled according to the level of the simulated x-ray signals. Results: The proposed method is shown to match real scans in number of experiments. Noise standard deviation measurements in simulated low-dose reconstructions of a 35 cm water phantom match real scans in a range from 500 to 10 mA with less than 5% error. Mean and variance of individual detector channels are shown to match closely across the detector array. Finally, the visual appearance of noise and streak artifacts is shown to match in real scans even under conditions of photon-starvation (with tube currents as low as 10 and 80 mA). Additionally, the proposed method is shown to be more accurate than previous approaches (1) in achieving the correct mean and variance in reconstructed images from pure-Poisson noise simulations (with no detector noise) under photon-starvation conditions, and (2) in simulating the correct noise level and detector noise artifacts in real low-dose scans. Conclusions: The proposed method can accurately simulate low-dose CT data starting from high-dose data, including effects from photon starvation and detector noise. This is potentially a very useful tool in helping to determine minimum dose requirements for a wide range of clinical protocols and advanced reconstruction algorithms.

60 APPLIED LIFE SCIENCES↗

Physically-based parameterization of spatially variable soil and vegetation using satellite multispectral data

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

Jasinski, Michael F.↗

Simulation studies of a full-ring, CZT SPECT system for whole-body imaging of 99m Tc and 177 Lu

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. With the goal of filling this technological gap, we have designed a cadmium zinc telluride (CZT) full-ring SPECT scanner instrumented with a broad-energy tungsten collimator. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT does not rely on the dual- and the triple-head rotational SPECT standard paradigm, enabling a larger effective area in each scan to increase the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized phantoms filled with 99m Tc and 177 Lu. Our SPECT design can resolve 7.9 mm rods for 99m Tc (140 keV) and 9.5 mm for 177 Lu (208 keV) in a hot-rod Derenzo phantom with a 3-min exposure and reach an image contrast of 78% for 99m Tc and 57% for 177 Lu using the NEMA IQ phantom with a 6-min exposure. Our modified scatter correction shows an improved contrast-recovery ratio compared to a standard correction. In this paper, we demonstrate the good performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99m Tc imaging of brain perfusion and 123 I imaging of dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design provides similar IQ and contrast to the commercial full-ring SPECT VERITON for 99m Tc. Regarding 177 Lu imaging of the 208 keV emissions, our design provides similar contrast to that of other state-of-the-art SPECTs with a significant reduction in exposure. In conclusion, the high sensitivity and extended energy range up to 250 keV makes our SPECT design a promising alternative for clinical imaging and theranostics of emerging radionuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Imaging an event horizon: mitigation of scattering toward Sagittarius A*

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ∼50 μas. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry. However, strong-field GR features of interest will be blurred at λ ≥ 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-G-IeP1-08: MR Geometric Distortion Dependency On Imaging Sequence, Acquisition Orientation and Receiver Bandwidth of a Dedicated 1.5T MR-Simulator

Purpose: To investigate the 3D geometric distortion of four potential MR sequences for radiotheraptic applications, and its dependency on sequence-type, acquisition-orientation and receiver-bandwidth from a dedicated 1.5T 700mm-wide bore MR-simulator (Magnetom-Aera, Sienmens Healthcare, Erlangen, Germany), using a large customized geometric accuracy phantom. Methods: This work studied 3D gradient-echo (VIBE) and spin-echo (SPACE) sequences for anatomical imaging; a specific ultra-short-TE sequence (PETRA) potentially for bone imaging and MR-based dosimetry; and a motion-insensitive sequence (BLADE) for dynamic applications like 4D-MRI. Integrated geometric-correction was employed, three orthogonal acquisition-orientations and up to three receiver-bandwidths were used, yielding 27 acquisitions for testing (Table 1a).A customized geometric accuracy phantom (polyurethane, MR/CT invisible, W×L×H:55×55×32.5cm3) was constructed and filled with 3892 spherical markers (6mm diameter, MR/CT visible) arranged on a 25mm-interval 3D isotropic-grid (Fig.1). The marker positions in MR images were quantitatively calculated and compared against those in the CT-reference using customized MatLab scripts. Results: The average distortion within various diameter-of-spherical-volumes (DSVs) and the usable DSVs under various distortion limits were measured (Tables 1b-c). It was observed that distortions fluctuated when sequence-type, acquisition-orientation or receiver-bandwidth changed (e.g. within 300mm-DSV, the lowest/highest average distortions of VIBE were 0.40mm/0.59mm, a 47.5% difference). According to AAPM-TG66 (<1mm distortion, left-most column of Table 1c), PETRA (Largest-DSV:253.9mm) has the potential on brain treatment, while BLADE (Largest-DSV:207.2mm) may need improvement for thoracic/abdominal applications. The results of VIBE (Largest-DSVs:294.3mm, the best among tested acquisitions) and SPACE (Largest-DSVs:267.7mm) suggests their potentials on head and neck applications. These Largest-DSVs were attained on different acquisition-orientations and receiver-bandwidths. Conclusion: Geometric distortion was shown to be dependent on sequence-type, acquisition-orientation and receiver-bandwidth. In the experiment, no configuration in any one of these factors could consistently reduce distortion while the others were varying. The distortion analysis result is a valuable guideline for sequence selection and optimization for MR-aided radiotherapy applications.

60 APPLIED LIFE SCIENCES↗

Automatic CT simulation optimization for radiation therapy: A general strategy

Purpose: In radiation therapy, x-ray computed tomography (CT) simulation protocol specifications should be driven by the treatment planning requirements in lieu of duplicating diagnostic CT screening protocols. The purpose of this study was to develop a general strategy that allows for automatically, prospectively, and objectively determining the optimal patient-specific CT simulation protocols based on radiation-therapy goals, namely, maintenance of contouring quality and integrity while minimizing patient CT simulation dose. Methods: The authors proposed a general prediction strategy that provides automatic optimal CT simulation protocol selection as a function of patient size and treatment planning task. The optimal protocol is the one that delivers the minimum dose required to provide a CT simulation scan that yields accurate contours. Accurate treatment plans depend on accurate contours in order to conform the dose to actual tumor and normal organ positions. An image quality index, defined to characterize how simulation scan quality affects contour delineation, was developed and used to benchmark the contouring accuracy and treatment plan quality within the predication strategy. A clinical workflow was developed to select the optimal CT simulation protocols incorporating patient size, target delineation, and radiation dose efficiency. An experimental study using an anthropomorphic pelvis phantom with added-bolus layers was used to demonstrate how the proposed prediction strategy could be implemented and how the optimal CT simulation protocols could be selected for prostate cancer patients based on patient size and treatment planning task. Clinical IMRT prostate treatment plans for seven CT scans with varied image quality indices were separately optimized and compared to verify the trace of target and organ dosimetry coverage. Results: Based on the phantom study, the optimal image quality index for accurate manual prostate contouring was 4.4. The optimal tube potentials for patient sizes of 38, 43, 48, 53, and 58 cm were 120, 140, 140, 140, and 140 kVp, respectively, and the corresponding minimum CTDIvol for achieving the optimal image quality index 4.4 were 9.8, 32.2, 100.9, 241.4, and 274.1 mGy, respectively. For patients with lateral sizes of 43–58 cm, 120-kVp scan protocols yielded up to 165% greater radiation dose relative to 140-kVp protocols, and 140-kVp protocols always yielded a greater image quality index compared to the same dose-level 120-kVp protocols. The trace of target and organ dosimetry coverage and the γ passing rates of seven IMRT dose distribution pairs indicated the feasibility of the proposed image quality index for the predication strategy. Conclusions: A general strategy to predict the optimal CT simulation protocols in a flexible and quantitative way was developed that takes into account patient size, treatment planning task, and radiation dose. The experimental study indicated that the optimal CT simulation protocol and the corresponding radiation dose varied significantly for different patient sizes, contouring accuracy, and radiation treatment planning tasks.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow

To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54deg turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJ/kg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilising the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly underpredicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJ/kg condition, but with up to a factor of three overprediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJ/kg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.

expansion tube↗

Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow

To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54 turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJkg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilizing the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly under predicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJkg condition, but with up to a factor of three over prediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJkg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.

Earth entry↗

A Framework Based on 2-D Taylor Expansion for Quantifying the Impacts of Subpixel Reflectance Variance and Covariance on Cloud Optical Thickness and Effective Radius Retrievals Based on the Bispectral Method

The bispectral method retrieves cloud optical thickness (τ) and cloud droplet effective radius (re) simultaneously from a pair of cloud reflectance observations, one in a visible or near-infrared (VIS/NIR) band and the other in a shortwave infrared (SWIR) band. A cloudy pixel is usually assumed to be horizontally homogeneous in the retrieval. Ignoring subpixel variations of cloud reflectances can lead to a significant bias in the retrieved τ and re. In the literature, the retrievals of τ and re are often assumed to be independent and considered separately when investigating the impact of subpixel cloud reflectance variations on the bispectral method. As a result, the impact on τ is contributed only by the subpixel variation of VIS/NIR band reflectance and the impact on re only by the subpixel variation of SWIR band reflectance. In our new framework, we use the Taylor expansion of a two-variable function to understand and quantify the impacts of subpixel variances of VIS/NIR and SWIR cloud reflectances and their covariance on the τ and re retrievals. This framework takes into account the fact that the retrievals are determined by both VIS/NIR and SWIR band observations in a mutually dependent way. In comparison with previous studies, it provides a more comprehensive understanding of how subpixel cloud reflectance variations impact the τ and re retrievals based on the bispectral method. In particular, our framework provides a mathematical explanation of how the subpixel variation in VIS/NIR band influences the re retrieval and why it can sometimes outweigh the influence of variations in the SWIR band and dominate the error in re retrievals, leading to a potential contribution of positive bias to the re retrieval. We test our framework using synthetic cloud fields from a large-eddy simulation and real observations from Moderate Resolution Imaging Spectroradiometer. The predicted results based on our framework agree very well with the numerical simulations. Our framework can be used to estimate the retrieval uncertainty from subpixel reflectance variations in operational satellite cloud products and to help understand the differences in τ and re retrievals between two instruments.

reflectance↗

Plasma Pressure in the Topside Ionosphere

A previous three year NASA-funded project resulted in the first 2-D maps of magnetotail pressure, density and temperature. A proposal to continue the work was declined, but modest funding was provided for one year to ramp down of the work. During the phase-out year, we used a time when 5 DMSP satellites were simultaneously active to produce the first instantaneous partial image of the magnetotail. The results have been submitted to the proceedings of the 1998 Huntsville Meeting on "The New Millennium Magnetosphere: Integrating Imaging, Discrete Observations and Global Simulations". A method of inferring central plasma sheet (CPS) temperature, density, and pressure from ionospheric observations was developed under a previous 3-year grant. These particles properties are calculated from data taken by particle instruments on DMSP satellites. Ion spectra occurring in conjunction with electron acceleration events are excluded. Because of the variability of magnetotail stretching, mapping to the plasma sheet was done using a modified Tsyganenko 1989 magnetic field model adjusted to agree with the actual magnetotail stretch. On May 25, 1997, five DMSP satellites (F10-F14) passed through the southern hemisphere nightside oval within a 19 minute period. Attached is the first magnetotail image, which results from applying our technique to that data set.

Newell, Patrick T.↗

Quantitative Monte Carlo-based holmium-166 SPECT reconstruction

Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 ({sup 166}Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative {sup 166}Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of {sup 166}Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full {sup 166}Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (A{sup est}) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six {sup 166}Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80% (SPECT-ppMC+DSW) to 76%–103% (SPECT-fMC). Furthermore, SPECT-fMC recovered whole-body activities were most accurate (A{sup est}= 1.06 × A − 5.90 MBq, R{sup 2}= 0.97) and SPECT-fMC tumor absorbed doses were significantly higher than with SPECT-DSW (p = 0.031) and SPECT-ppMC+DSW (p = 0.031).Conclusions: The quantitative accuracy of {sup 166}Ho SPECT is improved by Monte Carlo-based modeling of the image degrading factors. Consequently, the proposed reconstruction method enables accurate estimation of the radiation absorbed dose in clinical practice.

61 RADIATION PROTECTION AND DOSIMETRY↗

Observation simulation experiments with regional prediction models

Research efforts in FY 1990 included studies employing regional scale numerical models as aids in evaluating potential contributions of specific satellite observing systems (current and future) to numerical prediction. One study involves Observing System Simulation Experiments (OSSEs) which mimic operational initialization/forecast cycles but incorporate simulated Advanced Microwave Sounding Unit (AMSU) radiances as input data. The objective of this and related studies is to anticipate the potential value of data from these satellite systems, and develop applications of remotely sensed data for the benefit of short range forecasts. Techniques are also being used that rely on numerical model-based synthetic satellite radiances to interpret the information content of various types of remotely sensed image and sounding products. With this approach, evolution of simulated channel radiance image features can be directly interpreted in terms of the atmospheric dynamical processes depicted by a model. Progress is being made in a study using the internal consistency of a regional prediction model to simplify the assessment of forced diabatic heating and moisture initialization in reducing model spinup times. Techniques for model initialization are being examined, with focus on implications for potential applications of remote microwave observations, including AMSU and Special Sensor Microwave Imager (SSM/I), in shortening model spinup time for regional prediction.

Diak, George↗

Automatic Segmentation and Online virtualCT in Head-and-Neck Adaptive Radiation Therapy

Purpose: The purpose of this work was to develop and validate an efficient and automatic strategy to generate online virtual computed tomography (CT) scans for adaptive radiation therapy (ART) in head-and-neck (HN) cancer treatment. Method: We retrospectively analyzed 20 patients, treated with intensity modulated radiation therapy (IMRT), for an HN malignancy. Different anatomical structures were considered: mandible, parotid glands, and nodal gross tumor volume (nGTV). We generated 28 virtualCT scans by means of nonrigid registration of simulation computed tomography (CTsim) and cone beam CT images (CBCTs), acquired for patient setup. We validated our approach by considering the real replanning CT (CTrepl) as ground truth. We computed the Dice coefficient (DSC), center of mass (COM) distance, and root mean square error (RMSE) between correspondent points located on the automatically segmented structures on CBCT and virtualCT. Results: Residual deformation between CTrepl and CBCT was below one voxel. Median DSC was around 0.8 for mandible and parotid glands, but only 0.55 for nGTV, because of the fairly homogeneous surrounding soft tissues and of its small volume. Median COM distance and RMSE were comparable with image resolution. No significant correlation between RMSE and initial or final deformation was found. Conclusion: The analysis provides evidence that deformable image registration may contribute significantly in reducing the need of full CT-based replanning in HN radiation therapy by supporting swift and objective decision-making in clinical practice. Further work is needed to strengthen algorithm potential in nGTV localization.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Interval From Imaging to Treatment Delivery in the Radiation Surgery Age: How Long Is Too Long?

Purpose: The purpose of this study was to evaluate workflow and patient outcomes related to frameless stereotactic radiation surgery (SRS) for brain metastases. Methods and Materials: We reviewed all treatment demographics, clinical outcomes, and workflow timing, including time from magnetic resonance imaging (MRI), computed tomography (CT) simulation, insurance authorization, and consultation to the start of SRS for brain metastases. Results: A total of 82 patients with 151 brain metastases treated with SRS were evaluated. The median times from consultation, insurance authorization, CT simulation, and MRI for treatment planning were 15, 7, 6, and 11 days to SRS. Local freedom from progression (LFFP) was lower in metastases with MRI ≥14 days before treatment (P=.0003, log rank). The 6- and 12-month LFFP rate were 95% and 75% for metastasis with interval of <14 days from MRI to treatment compared to 56% and 34% for metastases with MRI ≥14 days before treatment. On multivariate analysis, LFFP remained significantly lower for lesions with MRI ≥14 days at SRS (P=.002, Cox proportional hazards; hazard ratio: 3.4, 95% confidence interval: 1.6-7.3). Conclusions: Delay from MRI to SRS treatment delivery for brain metastases appears to reduce local control. Future studies should monitor the timing from imaging acquisition to treatment delivery. Our experience suggests that the time from MRI to treatment should be <14 days.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulations of coherent scattering experiments at storage ring synchrotron radiation sources in the hard x-ray range

Detailed simulations of experiments carried out at modern light sources are directly related to the most efficient and productive use of these facilities for research in multiple branches of science and technology. The “Synchrotron Radiation Workshop” computer code with its Python interface, and Sirepo web-browser-based graphical user interface, currently supports physical optics simulations of coherent X-ray scattering and imaging experiments on user-defined virtual samples. We present examples of simulations of coherent scattering experiments that are typically performed at the Coherent Hard X-ray beamline at Brookhaven National Laboratory’s (BNL) National Synchrotron Light Source II. We also present several comparisons of the simulations with the results of actual coherent X-ray scattering experiments with nano-fabricated test samples produced at BNL’s Center for Functional Nanomaterials.

36 MATERIALS SCIENCE↗

In-Situ Visualization Experiments with ParaView Cinema in RAGE

A previous paper described some numerical experiments performed using the ParaView/Catalyst in-situ visualization infrastructure deployed in the Los Alamos RAGE radiation-hydrodynamics code to produce images from a running large scale 3D ICF simulation. One challenge of the in-situ approach apparent in these experiments was the difficulty of choosing parameters likes isosurface values for the visualizations to be produced from the running simulation without the benefit of prior knowledge of the simulation results and the resultant cost of recomputing in-situ generated images when parameters are chosen suboptimally. A proposed method of addressing this difficulty is to simply render multiple images at runtime with a range of possible parameter values to produce a large database of images and to provide the user with a tool for managing the resulting database of imagery. Recently, ParaView/Catalyst has been extended to include such a capability via the so-called Cinema framework. Here I describe some initial experiments with the first delivery of Cinema and make some recommendations for future extensions of Cinema’s capabilities.

in-situ visualization↗

Domain Adaptation for Measurements of Strong Gravitational Lenses

Upcoming surveys are predicted to discover galaxy-scale strong lenses on the magnitude of 105, making deep learning methods necessary in lensing data analysis. Currently, there is insufficient real lensing data to train deep learning algorithms, but training only on simulated data results in poor performance on real data. Domain adaptation can bridge the gap between simulated and real datasets. We adopt domain adaptation on the estimation of Einstein radius in simulated galaxy-scale gravitational lensing images. We evaluate two domain adaptation techniques - domain adversarial neural networks (DANN) and maximum mean discrepancy (MMD). We train on a source domain of simulated lenses and apply it to a target domain with emulation of DES survey conditions. We show that both domain adaptation techniques can significantly improve the model performance on the more complex target domain datasets. Our results show the potential of using domain adaptation to perform analysis on future survey data with a deep neural network trained on simulated data.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimization of an on-board imaging system for extremely rapid radiation therapy

Purpose: Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system. Methods: Three different CT detectors are proposed: low-resolution (all 4 × 4 mm pixels), medium-resolution (a combination of 4 × 4 mm pixels and 2 × 2 mm pixels), and high-resolution (all 1 × 1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called “realistic” images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a “gold standard” deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields—the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed. Results: By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm. Conclusions: The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

60 APPLIED LIFE SCIENCES↗