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,351 records · Page 75

Design of teleoperation system with a force-reflecting real-time simulator

We developed a force-reflecting teleoperation system that uses a real-time graphic simulator. This system eliminates the effects of communication time delays in remote robot manipulation. The simulator provides the operator with predictive display and feedback of computed contact forces through a six-degree of freedom (6-DOF) master arm on a real-time basis. With this system, peg-in-hole tasks involving round-trip communication time delays of up to a few seconds were performed at three support levels: a real image alone, a predictive display with a real image, and a real-time graphic simulator with computed-contact-force reflection and a predictive display. The experimental results indicate the best teleoperation efficiency was achieved by using the force-reflecting simulator with two images. The shortest work time, lowest sensor maximum, and a 100 percent success rate were obtained. These results demonstrate the effectiveness of simulated-force-reflecting teleoperation efficiency.

Hirata, Mitsunori↗

Benefits of adaptive radiation therapy in lung cancer as a function of replanning frequency

Purpose: To quantify the potential benefit associated with daily replanning in lung cancer in terms of normal tissue dose sparing and to characterize the tradeoff between adaptive benefit and replanning frequency. Methods: A set of synthetic images and contours, derived from weekly active breathing control images of 12 patients who underwent radiation therapy treatment for nonsmall cell lung cancer, is generated for each fraction of treatment using principal component analysis in a way that preserves temporal anatomical trends (e.g., tumor regression). Daily synthetic images and contours are used to simulate four different treatment scenarios: (1) a “no-adapt” scenario that simulates delivery of an initial plan throughout treatment, (2) a “midadapt” scenario that implements a single replan for fraction 18, (3) a “weekly adapt” scenario that simulates weekly adaptations, and (4) a “full-adapt” scenario that simulates daily replanning. An initial intensity modulated radiation therapy plan is created for each patient and replanning is carried out in an automated fashion by reoptimizing beam apertures and weights. Dose is calculated on each image and accumulated to the first in the series using deformable mappings utilized in synthetic image creation for comparison between simulated treatments. Results: Target coverage was maintained and cord tolerance was not exceeded for any of the adaptive simulations. Average reductions in mean lung dose (MLD) and volume of lung receiving 20 Gy or more (V20{sub lung}) were 65 ± 49 cGy (p = 0.000 01) and 1.1% ± 1.2% (p = 0.0006), respectively, for all patients. The largest reduction in MLD for a single patient was 162 cGy, which allowed an isotoxic escalation of the target dose of 1668 cGy. Average reductions in cord max dose, mean esophageal dose (MED), dose received by 66% of the heart (D66{sub heart}), and dose received by 33% of the heart (D33{sub heart}), were 158 ± 280, 117 ± 121, 37 ± 77, and 99 ± 120 cGy, respectively. Average incremental reductions in MLD for the midadapt, weekly adapt, and full-adapt treatments were 38, 18, and 8 cGy, respectively. Incremental reductions in MED for the same treatments were 57, 37, and 23 cGy. Reductions in MLD and MED for the full-adapt treatment were correlated with the absolute decrease in the planning target volume (r = 0.34 and r = 0.26). Conclusions: Adaptive radiation therapy for lung cancer yields clinically relevant reductions in normal tissue doses for frequencies of adaptation ranging from a single replan up to daily replanning. Increased frequencies of adaptation result in additional benefit while magnitude of benefit decreases.

60 APPLIED LIFE SCIENCES↗

SU-G-TeP4-04: An Automated Monte Carlo Based QA Framework for Pencil Beam Scanning Treatments

Purpose: Prior to treating new PBS field, multiple (three) patient-field-specific QA measurements are performed: two 2D dose distributions at shallow depth (M1) and at the tumor depth (M2) with treatment hardware at zero gantry angle; one 2D dose distribution at iso-center (M3) without patient specific devices at the planned gantry angle. This patient-specific QA could be simplified by the use of MC model. The results of MC model commissioning for a spot-scanning system and the fully automated TOPAS/MC-based QA framework will be presented. Methods: We have developed in-house MC interface to access a TPS (Astroid) database from a computer cluster remotely. Once a plan is identified, the interface downloads information for the MC simulations, such as patient images, apertures points, and fluence maps and initiates calculations in both the patient and QA geometries. The resulting calculations are further analyzed to evaluate the TPS dose accuracy and the PBS delivery. Results: The Monte Carlo model of our system was validated within 2.0 % accuracy over the whole range of the dose distribution (proximal/shallow part, as well as target dose part) due to the location of the measurements. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. Conclusion: As M1 depths range typically from 1 cm to 4 cm from the phantom surface, The Monte Carlo model of our system was validated within +− 2.0 % in absolute dose level over a whole treatment range. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. This work was supported by NIH/NCI under CA U19 21239.

60 APPLIED LIFE SCIENCES↗

EFFECT OF BACKGROUND MAGNETIC FIELD ON TURBULENCE DRIVEN BY MAGNETOROTATIONAL INSTABILITY IN ACCRETION DISKS

We investigate the background magnetic field dependence of the saturated state of a magnetorotational instability (MRI) in an accretion disk by performing three-dimensional magnetohydrodynamic simulations. We assume an unstratified disk by employing the local shearing box approximation. Three different uniform background magnetic field configurations are treated for a wide range of field intensities. These simulations indicate that the time variations of the turbulent stress and the magnetic energy are altered by the presence of a poloidal component of the background field. We find that the saturation amplitude of the turbulent stress and the magnetic energy are determined by both the poloidal and azimuthal components of the field. In particular, when the poloidal component has the same intensity, the obtained turbulent stress for {beta}{sub y0} Almost-Equal-To 200 becomes smaller than those for a purely poloidal field case. Despite the fact that the background field affects the MRI turbulence, the correlation between the obtained turbulent stress and the magnetic energy in the nonlinear stage is independent of the field topology. Our results indicate that the saturated turbulent stress has a stronger correlation with the power of the perturbed component of the magnetic field than with the power of the total magnetic field. These results suggest that both the intensity and the direction of the background magnetic field significantly affect the turbulent motion of the MRI in accretion disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of dual energy quantitative CT for determining the spatial distributions of red marrow and bone for dosimetry in internal emitter radiation therapy

Purpose: To evaluate a three-equation three-unknown dual-energy quantitative CT (DEQCT) technique for determining region specific variations in bone spongiosa composition for improved red marrow dose estimation in radionuclide therapy. Methods: The DEQCT method was applied to 80/140 kVp images of patient-simulating lumbar sectional body phantoms of three sizes (small, medium, and large). External calibration rods of bone, red marrow, and fat-simulating materials were placed beneath the body phantoms. Similar internal calibration inserts were placed at vertebral locations within the body phantoms. Six test inserts of known volume fractions of bone, fat, and red marrow were also scanned. External-to-internal calibration correction factors were derived. The effects of body phantom size, radiation dose, spongiosa region segmentation granularity [single (∼17 × 17 mm) region of interest (ROI), 2 × 2, and 3 × 3 segmentation of that single ROI], and calibration method on the accuracy of the calculated volume fractions of red marrow (cellularity) and trabecular bone were evaluated. Results: For standard low dose DEQCT x-ray technique factors and the internal calibration method, the RMS errors of the estimated volume fractions of red marrow of the test inserts were 1.2–1.3 times greater in the medium body than in the small body phantom and 1.3–1.5 times greater in the large body than in the small body phantom. RMS errors of the calculated volume fractions of red marrow within 2 × 2 segmented subregions of the ROIs were 1.6–1.9 times greater than for no segmentation, and RMS errors for 3 × 3 segmented subregions were 2.3–2.7 times greater than those for no segmentation. Increasing the dose by a factor of 2 reduced the RMS errors of all constituent volume fractions by an average factor of 1.40 ± 0.29 for all segmentation schemes and body phantom sizes; increasing the dose by a factor of 4 reduced those RMS errors by an average factor of 1.71 ± 0.25. Results for external calibrations exhibited much larger RMS errors than size matched internal calibration. Use of an average body size external-to-internal calibration correction factor reduced the errors to closer to those for internal calibration. RMS errors of less than 30% or about 0.01 for the bone and 0.1 for the red marrow volume fractions would likely be satisfactory for human studies. Such accuracies were achieved for 3 × 3 segmentation of 5 mm slice images for: (a) internal calibration with 4 times dose for all size body phantoms, (b) internal calibration with 2 times dose for the small and medium size body phantoms, and (c) corrected external calibration with 4 times dose and all size body phantoms. Conclusions: Phantom studies are promising and demonstrate the potential to use dual energy quantitative CT to estimate the spatial distributions of red marrow and bone within the vertebral spongiosa.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A finite mixtures algorithm for finding classes in images

The authors describe an algorithm for estimating the proportions of classes in a SAR (synthetic aperture radar) image without any user interaction. The method assumes that the image is a mixture of a known number of different pixel types. A maximum likelihood estimate of the parameters of the resulting mixture distribution is then used to find the proportions for the various classes. The technique was successfully applied to aircraft SAR images of sea ice. Computer simulations were used to determine the relative errors of the technique. The technique performs well even with extremely noisy images.

Samadani, Ramin↗

Evaluation of Techniques for Vegetation Removal from UAV-Based Photogrammetric DSM near WIPP Land Withdrawal Act Boundary, NM

The Waste Isolation Pilot Plant is located within New Mexico's karst landscape and is the United States' only deep geologic radioactive waste repository which isolates transuranic waste from defense activities underground in a bedded salt formation. Characteristic karst topography in this region is made up of material (such as limestone) that can be dissolved by water over a long period of time, and includes features such as springs, surface streams, sinking streams, caves, and sinkholes. Karst features can potentially impact the integrity of the waste repository in the future and compromise WIPP's performance, particularly when coupled with incompatible land-use activities within the Land Withdrawal Act (LWA) boundary as well as increased water withdrawals outside the LWA. In order to better capture the ground surface topography to support hydrology modeling efforts, the U.S. DOE has enlisted FIU to develop a high-resolution digital elevation model (DEM). In a previous field study, a high-resolution digital surface model (DSM) of an adjacent representative site, Basin 6, was developed by photogrammetry using aerial images captured with unmanned aerial vehicles. The DSM generated was based on point clouds representing the entire landscapes, including points of terrain, vegetation and infrastructure. These above-ground features need to be removed to create a DSM/DEM with good representation of the bare ground. In recent years, many techniques including both pre- and post-processing of DSMs, were developed for vegetation removal with varied success. This research focuses on applying and comparing various vegetation removal methods using point clouds generated of the Basin 6 pilot study area. The 3D RGB-based point cloud classification using Python scripts was a modification of the procedure by Themistocleous (2019). The method using Liblas classifies the point cloud three-dimensionally, whereas the method by Themistocleous (2019) is two-dimensional. The preliminary result of the method using Liblas, however, still had a large amount of vegetation remaining compared to the Pix4D Machine Learning method. The vegetation layer remains highly visible in the DEMs processed using the vegetation indices as opposed to the Pix4D Machine Learning method, which produces a DEM where mostly bare ground is seen. Further modifications to the input parameters are required to improve the vegetation vs. bare ground classification using this approach. The next step will be to improve the method using Liblas and to also test other vegetation removal methods found in the literature review to compare them against the two described here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Four-Dimensional Dose Reconstruction for Scanned Proton Therapy Using Liver 4DCT-MRI

Purpose: Four-dimensional computed tomography-magnetic resonance imaging (4DCT-MRI) is an image-processing technique for simulating many 4DCT data sets from a static reference CT and motions extracted from 4DMRI studies performed using either volunteers or patients. In this work, different motion extraction approaches were tested using 6 liver cases, and a detailed comparison between 4DCT-MRI and 4DCT was performed. Methods and Materials: 4DCT-MRI has been generated using 2 approaches. The first approach used motion extracted from 4DMRI as being “most similar” to that of 4DCT from the same patient (subject-specific), and the second approach used the most similar motion obtained from a motion library derived from 4DMRI liver studies of 13 healthy volunteers (population-based). The resulting 4DCT-MRI and 4DCTs were compared using scanned proton 4D dose calculations (4DDC). Results: Dosimetric analysis showed that 93% ± 8% of points inside the clinical target volume (CTV) agreed between 4DCT and subject-specific 4DCT-MRI (gamma analysis: 3%/3 mm). The population-based approach however showed lower dosimetric agreement with only 79% ± 14% points in the CTV reaching the 3%/3 mm criteria. Conclusions: 4D CT-MRI extends the capabilities of motion modeling for dose calculations by accounting for realistic and variable motion patterns, which can be directly employed in clinical research studies. We have found that the subject-specific liver modeling appears more accurate than the population-based approach. The former is particularly interesting for clinical applications, such as improved target delineation and 4D dose reconstruction for patient-specific QA to allow for inter- and/or intra-fractional plan corrections.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Imaging surface topography with coherent x-ray reflectivity: Theory, kinematics, and simulations

A theoretical formalism is described for understanding coherent x-ray reflectivity (CXR) from the surface of a semi-infinite crystal having a variable surface topography, described by the height profile ℎ(𝑥,𝑦). The surface topography is imaged as a complex “effective density,” obtained from the phasing and inversion of the coherent x-ray reflectivity data, measured through a rocking scan centered at a vertical momentum transfer 𝑄$^{0}_{𝑧}$ and a vertical range Δ⁢𝑄 𝑧 . The formalism predicts that the effective density has an amplitude with a maximum located at the surface height for each position within the surface plane. The phase of the effective density has a lateral variation that is controlled by the surface height and a vertical variation that reflects a combination of the interfacial structure and specific choice of measurement conditions. This understanding enables direct observation of nanometer-scale interfacial topography, i.e., ℎ⁡(𝑥,𝑦)⁢𝑐 𝑠 (where 𝑐 𝑠 is the vertical substrate lattice parameter) with Å-scale sensitivity to surface height. Numerical simulations illustrate and confirm the theoretical results. These results show how the interpretation of the interfacial density phase obtained by CXR data inversion (i.e., surface topography with respect to a flat surface) is conceptually similar to that previously known for Bragg coherent diffraction imaging (BCDI) measurements of isolated nanoparticles (i.e., lattice displacements with respect to an ideal crystal lattice). This suggests that CXR can be thought of as a form of dark field imaging with respect to the bright field BCDI approach. An implication of these results is that interfacial imaging may bypass some of the significant challenges associated with BCDI imaging of multiple particles having different orientations.

X-ray imaging↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

PISCES High Contrast Integral Field Spectrograph Simulations and Data Reduction Pipeline

The PISCES (Prototype Imaging Spectrograph for Coronagraphic Exoplanet Studies) is a lenslet array based integral field spectrograph (IFS) designed to advance the technology readiness of the WFIRST (Wide Field Infrared Survey Telescope)-AFTA (Astrophysics Focused Telescope Assets) high contrast Coronagraph Instrument. We present the end to end optical simulator and plans for the data reduction pipeline (DRP). The optical simulator was created with a combination of the IDL (Interactive Data Language)-based PROPER (optical propagation) library and Zemax (a MatLab script), while the data reduction pipeline is a modified version of the Gemini Planet Imager's (GPI) IDL pipeline. The simulations of the propagation of light through the instrument are based on Fourier transform algorithms. The DRP enables transformation of the PISCES IFS data to calibrated spectral data cubes.

The PISCES (Prototype Imaging Spectrograph for Cor↗

A 3D MR-acquisition scheme for nonrigid bulk motion correction in simultaneous PET-MR

Purpose: Positron emission tomography (PET) is a highly sensitive medical imaging technique commonly used to detect and assess tumor lesions. Magnetic resonance imaging (MRI) provides high resolution anatomical images with different contrasts and a range of additional information important for cancer diagnosis. Recently, simultaneous PET-MR systems have been released with the promise to provide complementary information from both modalities in a single examination. Due to long scan times, subject nonrigid bulk motion, i.e., changes of the patient's position on the scanner table leading to nonrigid changes of the patient's anatomy, during data acquisition can negatively impair image quality and tracer uptake quantification. A 3D MR-acquisition scheme is proposed to detect and correct for nonrigid bulk motion in simultaneously acquired PET-MR data. Methods: A respiratory navigated three dimensional (3D) MR-acquisition with Radial Phase Encoding (RPE) is used to obtain T1- and T2-weighted data with an isotropic resolution of 1.5 mm. Healthy volunteers are asked to move the abdomen two to three times during data acquisition resulting in overall 19 movements at arbitrary time points. The acquisition scheme is used to retrospectively reconstruct dynamic 3D MR images with different temporal resolutions. Nonrigid bulk motion is detected and corrected in this image data. A simultaneous PET acquisition is simulated and the effect of motion correction is assessed on image quality and standardized uptake values (SUV) for lesions with different diameters. Results: Six respiratory gated 3D data sets with T1- and T2-weighted contrast have been obtained in healthy volunteers. All bulk motion shifts have successfully been detected and motion fields describing the transformation between the different motion states could be obtained with an accuracy of 1.71 ± 0.29 mm. The PET simulation showed errors of up to 67% in measured SUV due to bulk motion which could be reduced to less than 10% with the proposed motion compensation approach. Conclusions: A MR acquisition scheme which yields both high resolution 3D anatomical data and highly accurate nonrigid motion information without an increase in scan time is presented. The proposed method leads to a strong improvement in both MR and PET image quality and ensures an accurate assessment of tracer uptake.

60 APPLIED LIFE SCIENCES↗

Simultaneous reduction of radiation dose and scatter for CBCT by using collimators

Purpose: On-board cone-beam CT (CBCT) imaging has been widely available in radiotherapy clinic for target localization. However, the extra radiation dose from CBCT is always a concern for its frequent use. Additionally, the relatively large scatter in CBCT often degrades the image quality. By using collimators, some of the X-rays can be stopped from reaching the patient and the detectors, hence both the scatter and the patient doses are simultaneously reduced. The authors show in this work that the collimated CBCT data can be reconstructed without any noticeable artifacts for certain collimator blocking ratios and blocking patterns, and the focus of this work is to study the relationship between the image quality and these two collimator factors.Methods: A CBCT system with collimators was simulated following the typical geometry used in clinic. Different collimator designs were tested by varying the size and the number of the collimator slits, and at the same time, the ratio of transmitted beams to total beams was varied from 100% to 10%, resulting in hundreds of different simulation scenarios. Lung and pelvis phantoms created from patients CT images were used in the simulations, and an iterative reconstruction algorithm using the compressed sensing technique was adopted. The image quality was examined by root mean square errors (RMSEs) and compared with the conventional CBCT images.Results: The CBCT image quality increases as the amount of beams passing through the collimators increases, and decreases as the size of the collimator slits increases. With ultra-high resolution collimators, the RMSEs were comparable to the conventional CBCT image quality until the beam transmission rate is reduced below 25%.Conclusions: Collimators can reduce the scatters and radiation dose, however, the collimated CBCT image quality is strongly dependent on both the collimator blocking ratio and the blocking pattern. To achieve image quality comparable to the conventional CBCT, the amount of information and data format must be adequate.

61 RADIATION PROTECTION AND DOSIMETRY↗

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗

Atomic Fingerprinting of Heteroatoms Using Noncontact Atomic Force Microscopy

We report immense strides have been made in increasing the resolution of scanning probe microscopy. Noncontact atomic force microscopy (n-AFM) now offers one the ability to characterize and visualize single molecules with subatomic resolution. Specifically, nc-AFM with a carbon monoxide (CO) functionalized tip has the ability to discriminate functional groups (–C≡C–, –CH 2 , –C$=$O, …), although discriminating atomic species often remains as an ongoing challenge. Here, real-space pseudopotentials constructed within density functional theory are employed to accurately simulate nc-AFM images of molecules containing heteroatoms (S, I, and N) within dibenzothiophene (DBT), 2-iodotriphenylene (ITP), acridine (ACR) and ferrous phthalocyanine (FePc). It is found that S and I atoms can be easily identified from C based on their unique features. For N atoms, a use of tip functionalization is proposed to effectively discriminate them from C atoms.

36 MATERIALS SCIENCE↗

Numerical simulations of early-stage dynamics of electron bunches emitted from plasmonic photocathodes

High-brightness electron sources are a key ingredient to the development of compact accelerator-based light sources. The electron sources are commonly based on (linear) photoemission process where a laser pulse with proper wavelength impinges on the surface of a metallic or semiconductor photocathode. Very recently, the use of plasmonic cathodes - cathodes with a nano-patterned surface - have demonstrated great enhancement in quantum efficiencies (Li et al., 2013 [1]). Alternatively, this type of photocathodes could support the formation of structured beams composed of transversely-separated beamlets. In this paper we discuss numerical simulations of the early-stage beam dynamics of the emission process from plasmonic cathodes carried out using the Warp (Friedman et al., 2014 [2]) framework. Furthermore, the model is used to investigate the properties of beams emitted from these photocathode and subsequently combined with particle-in-cell simulations to explore the imaging of cathode pattern after acceleration in a radiofrequency gun.

43 PARTICLE ACCELERATORS↗