Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial interpolation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

bifacialVF - Bifacial PV View Factor model [SWR 17-58]

bifacialvf is a self-contained view factor (or configuration factor) model which replicates a 5-row PV system of infinite extent perpendicular to the module rows. Single-axis tracking is supported, and hourly output files based on TMY inputs are saved. Spatial nonuniformity is reported, with multiple rear-facing irradiances collected on the back of each module row. Bilinear interpolation code add-on to bifacialvf (description below) to pre-generate IV arrays and bifacial coefficients, and to examine the energy production with back side irradiance mismatch for either a portrait or landscape module. Included are IV curves and bifacial info for a Yingli (standard) module. Based on the publication: "A Practical Irradiance Model for Bifacial PV Modules" B. Marion, S. MacAlpine, C. Deline, A. Asgharzadeh, F. Toor, D. Riley, J. Stein, C. Hansen 2017 IEEE Photovoltaic Specialists Conference, Washington DC, 2017 URL: https://www.nrel.gov/docs/fy17osti/67847.pdf

Deline, Chris↗

bifacialVF - Bifacial PV View Factor model 2022 [SWR 17-58]

bifacialvf is a self-contained view factor (or configuration factor) model which replicates a 5-row PV system of infinite extent perpendicular to the module rows. Single-axis tracking is supported, and hourly output files based on TMY inputs are saved. Spatial nonuniformity is reported, with multiple rear-facing irradiances collected on the back of each module row. Bilinear interpolation code add-on to bifacialvf (description below) to pre-generate IV arrays and bifacial coefficients, and to examine the energy production with back side irradiance mismatch for either a portrait or landscape module. Included are IV curves and bifacial info for a Yingli (standard) module. Based on the publication: "A Practical Irradiance Model for Bifacial PV Modules" B. Marion, S. MacAlpine, C. Deline, A. Asgharzadeh, F. Toor, D. Riley, J. Stein, C. Hansen 2017 IEEE Photovoltaic Specialists Conference, Washington DC, 2017 URL: https://www.nrel.gov/docs/fy17osti/67847.pdf

Deline, Chris↗

Deep Learning-based Non-Stationary Bias Correction (NSBC)

This work develops the NSBC (non-stationary bias correction) methodology to correct temperature projection bias from E3SM. The NSBC deep learning framework consists of a three-part architecture: an auto-encoder for compressing the spatial information, an LSTM for predicting annual temperature mean, and a U-Net for capturing the residual bias in temperature. The non-stationary bias correction (NSBC) framework can correct the non-stationarity of the biases of the climate models, which significantly improves the accuracy of future temperature prediction and improves the overestimation of extreme high temperatures that many existing bias correction methods suffer from. Getting started 1. Obtain the historical climate simulation and observation data. The E3SM simulation data are available through https://aims2.llnl.gov/search/cmip6/. The pseudo observations, the Geophysical Fluid Dynamics Laboratory (GFDL)-ESM4 model (Krasting et al., 2018) are available through https://aims2.llnl.gov/search/cmip6/. The spatial resolution of E3SM and pseudo observation datasets are both regridded to a common 1° resolution grid using conservative interpolation. The regridded E3SM and pseudo observation with 1° resolution can be found throught ./data/. 2. Train the Auto-encoder model. Python 0-autoencoder.py 3. Train the LSTM Python 1-LSTM.py 4. Generate the annual mean temperature based on trained LSTM Python 2-generate_annual_mean_LSTM.py 5. Train the U-Net. Python 3-unet.py 6. Evaluation and compared with the baseline Python 4_evaluation.py Is there a deadline approaching that requires the release of yo

Lucas, Donald↗

Enhancing control of grid distribution in algebraic grid generation

Three techniques are presented to enhance the control of grid-point distribution for a class of algebraic grid generation methods known as the two-, four- and six-boundary methods. First, multidimensional stretching functions are presented, and a technique is devised to construct them based on the desired distribution of grid points along certain boundaries. Second, a normalization procedure is proposed which allows more effective control over orthogonality of grid lines at boundaries and curvature of grid lines near boundaries. And third, interpolating functions based on tension splines are introduced to control curvature of grid lines in the interior of the spatial domain. In addition to these three techniques, consistency conditions are derived which must be satisfied by all user-specified data employed in the grid generation process to control grid-point distribution. The usefulness of the techniques developed in this study was demonstrated by using them in conjunction with the two- and four-boundary methods to generate several grid systems, including a three-dimensional grid system in the coolant passage of a radial turbine blade with serpentine channels and pin fins.

Steinthorsson, E.↗

Data Reduction Pipeline for the CHARIS Integral-Field Spectrograph I: Detector Readout Calibration and Data Cube Extraction

We present the data reduction pipeline for CHARIS, a high-contrast integral-field spectrograph for the Subaru Telescope. The pipeline constructs a ramp from the raw reads using the measured nonlinear pixel response and reconstructs the data cube using one of three extraction algorithms: aperture photometry, optimal extraction, or chi-squared fitting. We measure and apply both a detector flatfield and a lenslet flatfield and reconstruct the wavelength- and position-dependent lenslet point-spread function (PSF) from images taken with a tunable laser. We use these measured PSFs to implement a chi-squared-based extraction of the data cube, with typical residuals of approximately 5 percent due to imperfect models of the under-sampled lenslet PSFs. The full two-dimensional residual of the chi-squared extraction allows us to model and remove correlated read noise, dramatically improving CHARIS's performance. The chi-squared extraction produces a data cube that has been deconvolved with the line-spread function and never performs any interpolations of either the data or the individual lenslet spectra. The extracted data cube also includes uncertainties for each spatial and spectral measurement. CHARIS's software is parallelized, written in Python and Cython, and freely available on github with a separate documentation page. Astrometric and spectrophotometric calibrations of the data cubes and PSF subtraction will be treated in a forthcoming paper.

Brandt, Timothy D.↗

Spatial characteristics of observed precipitation fields: A catalog of summer storms in Arizona, Volume 1

Eight years of summer raingage observations are analyzed for a dense, 93 gage, network operated by the U. S. Department of Agriculture, Agricultural Research Service, in their 150 sq km Walnut Gulch catchment near Tucson, Arizona. Storms are defined by the total depths collected at each raingage during the noon to noon period for which there was depth recorded at any of the gages. For each of the resulting 428 storms, the 93 gage depths are interpolated onto a dense grid and the resulting random field is anlyzed. Presented are: storm depth isohyets at 2 mm contour intervals, first three moments of point storm depth, spatial correlation function, spatial variance function, and the spatial distribution of total rainstorm depth.

Fennessey, N. M.↗

Edge detection - Image-plane versus digital processing

To optimize edge detection with the familiar Laplacian-of-Gaussian operator, it has become common to implement this operator with a large digital convolution mask followed by some interpolation of the processed data to determine the zero crossings that locate edges. It is generally recognized that this large mask causes substantial blurring of fine detail. It is shown that the spatial detail can be improved by a factor of about four with either the Wiener-Laplacian-of-Gaussian filter or an image-plane processor. The Wiener-Laplacian-of-Gaussian filter minimizes the image-gathering degradations if the scene statistics are at least approximately known and also serves as an interpolator to determine the desired zero crossings directly. The image-plane processor forms the Laplacian-of-Gaussian response by properly combining the optical design of the image-gathering system with a minimal three-by-three lateral-inhibitory processing mask. This approach, which is suggested by Marr's model of early processing in human vision, also reduces data processing by about two orders of magnitude and data transmission by up to an order of magnitude.

Huck, Friedrich O.↗

Numerical Method for Fractional Advection-Diffusion Equation with Heredity

We propose a method of construction of difference schemes for fractional partial differential equations with delay in time. For the fractional equation with two-sided diffusion, fractional transfer in time, and a functional aftereffect, we construct an implicit difference scheme. We use the shifted Grünwald–Letnikov formulas for the approximation of fractional derivatives with respect to spatial variables and the L1-algorithm for the approximation of fractional derivatives in time. Also we use piecewise constant interpolation and extrapolation by extending the discrete prehistory of the model in time. The algorithm is a fractional analog of a purely implicit method; on each time step, it is reduced to the solution of linear algebraic systems. We prove the stability of the method and find its order of convergence.

97 MATHEMATICS AND COMPUTING↗

Concept for the subresolution measurement of earthquake strain fields using SPOT panchromatic imagery

The use of satellite-derived imagery for measuring subresolution horizontal terrain displacments associated with present-day earthquakes is discussed with reference to data from the French SPOT satellite whose sensor array provides 10-m panchromatic imagery. The measured terrain displacements can be up to several meters, but usually no more than 6-8 m even for major earthquakes. The general approach is to spatially match the after image to the before image at each point on a half-kilometer grid by iteratively interpolating one and testing its correlation with the other. The discussion covers the basic algorithm, results of initial tests, error types and limitations, and future work.

Crippen, Robert E.↗

Issues in Data Fusion for Satellite Aerosol Measurements for Applications with GIOVANNI System at NASA GES DISC

We look at issues, barriers and approaches for Data Fusion of satellite aerosol data as available from the GES DISC GIOVANNI Web Service. Daily Global Maps of AOT from a single satellite sensor alone contain gaps that arise due to various sources (sun glint regions, clouds, orbital swath gaps at low latitudes, bright underlying surfaces etc.). The goal is to develop a fast, accurate and efficient method to improve the spatial coverage of the Daily AOT data to facilitate comparisons with Global Models. Data Fusion may be supplemented by Optimal Interpolation (OI) as needed.

Gopalan, Arun↗

Development of Multi-Sensor Global Cloud and Radiance Composites for DSCOVR EPIC Imager with Subpixel Definition

The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). EPIC delivers adequate spatial resolution imagery but only in shortwave bands (317-780 nm), while NISTAR measures the top-of-atmosphere (TOA) whole-disk radiance in shortwave and longwave broadband windows. Accurate calculation of albedo and outgoing longwave flux requires a high-resolution scene identification such as the radiance observations and cloud properties retrievals from low earth orbit (LEO, including NASA Terra and Aqua MODIS, Suomi-NPP VIIRS, and NOAA AVHRR) and geosynchronous (GEO, including GOES east and west, METEOSAT, INSAT-3D, MTSAT-2, and Himawari-8) satellite imagers. The cloud properties are derived using the Clouds and the Earth's Radiant Energy System (CERES) mission Cloud Subsystem group algorithms. These properties have to be co-located with EPIC pixels to provide the scene identification and to select anisotropic directional models (ADMs), which are then used to adjust the NISTAR-measured radiance and subsequently obtain the global daytime shortwave and longwave fluxes. This work presents an algorithm for optimal merging of selected radiance and cloud property parameters derived from multiple satellite imagers to obtain seamless global hourly composites at 5-km resolution. Selection of satellite data for each 5-km pixel is based on an aggregated rating that incorporates five parameters: nominal satellite resolution, pixel time relative to the EPIC time, viewing zenith angle, distance from day/night terminator, and probability of sun glint. To provide a smoother transition in the merged output, in regions where candidate pixel data from two satellite sources have comparable aggregated rating, the selection decision is defined by the cumulative function of the normal distribution so that abrupt changes in the visual appearance of the composite data are avoided. Higher spatial accuracy in the composite product is achieved by using the inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling.

Khlopenkov, Konstantin V.↗

Influence of strain and dislocations on GaSb/GaAs quantum dots: From nested to staggered band alignment

We investigate the influence of strain and dislocations on band alignment in GaSb/GaAs quantum dot systems. Composition profiles from cross-sectional scanning tunneling microscopy images are interpolated onto a finite element mesh in order to calculate the distribution of local elastic strain, which is converted to a spatially varying band alignment using deformation potential theory. Our calculations predict that dislocation-induced strain relaxation and charging lead to significant local variations in band alignment. Furthermore, misfit strain induces a transition from a nested (type I) to a staggered (type II) band alignment. Although dislocation-induced strain relaxation prevents the type I to type II transition, electrostatic charging at dislocations induces the staggered band alignment once again.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.

Lemieux, Jean-François (ORCID:0000000320845759)↗

Software for Generating Troposphere Corrections for InSAR Using GPS and Weather Model Data

Atmospheric errors due to the troposphere are a limiting error source for spaceborne interferometric synthetic aperture radar (InSAR) imaging. This software generates tropospheric delay maps that can be used to correct atmospheric artifacts in InSAR data. The software automatically acquires all needed GPS (Global Positioning System), weather, and Digital Elevation Map data, and generates a tropospheric correction map using a novel algorithm for combining GPS and weather information while accounting for terrain. Existing JPL software was prototypical in nature, required a MATLAB license, required additional steps to acquire and ingest needed GPS and weather data, and did not account for topography in interpolation. Previous software did not achieve a level of automation suitable for integration in a Web portal. This software overcomes these issues. GPS estimates of tropospheric delay are a source of corrections that can be used to form correction maps to be applied to InSAR data, but the spacing of GPS stations is insufficient to remove short-wavelength tropospheric artifacts. This software combines interpolated GPS delay with weather model precipitable water vapor (PWV) and a digital elevation model to account for terrain, increasing the spatial resolution of the tropospheric correction maps and thus removing short wavelength tropospheric artifacts to a greater extent. It will be integrated into a Web portal request system, allowing use in a future L-band SAR Earth radar mission data system. This will be a significant contribution to its technology readiness, building on existing investments in in situ space geodetic networks, and improving timeliness, quality, and science value of the collected data

Moore, Angelyn W.↗

Evaluation of Time-Series Gap-Filling Methods for Solar Irradiance Applications

A complete solar resource data set is essential for any stage of a solar energy project - from feasibility studies to daily operations. But measured or modeled solar resource data are prone to data gaps and data quality issues. To mitigate these issues, a data imputation process should be implemented to obtain a complete and reliable temporal and spatial data series. This study focused on imputing temporal scales by applying random and artificial data gaps and then implementing eight imputation methods, including the Kalman filtering and smoothing and stine interpolations. These methods were implemented on 1-minute to half hourly irradiance data for 1 year using a few locations from the National Solar Radiation Database (NSRDB) and ground measurement data set. The results demonstrated that some of the simpler methods, such as the stine and linear interpolation methods, were the relatively best models based on the statistical metrics for imputing NSRDB and ground measurement data, respectively.

14 SOLAR ENERGY↗

Estimates of Long Term Surface Soil Moisture in the Midwestern U.S. Derived from Satellite Microwave Observations

Soil moisture is a key component of the water and energy balances of the Earth's surface, and has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, soil moisture is often somewhat difficult to measure accurately in both space and time, especially at large spatial scales. Soil moisture is highly variable, and while point measurements are typically quite accurate, subsequent areal averaging of these measurements often leads to large errors. Since remotely sensed land surface observations are already a spatially averaged or areally integrated value, they are a logical input parameter to regional or larger scale land process models. A database of long-term soil moisture was compared to satellite microwave observations over test sites in the Midwestern United States. Ground measurements of average volumetric surface soil moisture in the top ten cm were made bimonthly at 19 locations throughout the state of Illinois. Nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and a temporal frequency over the test area of about 3 nighttime orbits per week. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Because the time of satellite coverage does not always coincide with the ground measurements of soil moisture, the existing ground data were used to calibrate a water balance for the top IO cm surface layer in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area average soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture, Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, M.↗

Estimating biotic integrity to capture existence value of freshwater ecosystems

The US Environmental Protection Agency (EPA) uses a water quality index (WQI) to estimate benefits of proposed Clean Water Act regulations. The WQI is relevant to human use value, such as recreation, but may not fully capture aspects of nonuse value, such as existence value. Here, we identify an index of biological integrity to supplement the WQI in a forthcoming national stated preference survey that seeks to capture existence value of streams and lakes more accurately within the conterminous United States (CONUS). We used literature and focus group research to evaluate aquatic indices regularly reported by the EPA’s National Aquatic Resource Surveys. We chose an index that quantifies loss in biodiversity as the observed-to-expected (O/E) ratio of taxonomic composition because focus group participants easily understood its meaning and the environmental changes that would result in incremental improvements. However, available datasets of this index do not provide the spatial coverage to account for how conditions near survey respondents affect their willingness to pay for its improvement. Therefore, we modeled and interpolated the values of this index from sampled sites to 1.1 million stream segments and 297,071 lakes across the CONUS to provide the required coverage. The models explained 13 to 36% of the variation in O/E scores and demonstrate how modeling can provide data at the required density for benefits estimation. We close by discussing future work to improve performance of the models and to link biological condition with water quality and habitat models that will allow us to forecast changes resulting from regulatory options.

59 BASIC BIOLOGICAL SCIENCES↗

NASTRAN supplemental documentation for modal forced vibration analysis of aerodynamically excited turbosystems

This report is a supplemental NASTRAN document for a new capability to determine the vibratory response of turbosystems subjected to aerodynamic excitation. Supplements of NASTRAN Theoretical, User's, Programmer's, and Demonstration Manuals are included. Turbosystems such as advanced turbopropellers with highly swept blades, and axial-flow compressors and turbines can be analyzed using this capability, which has been developed and implemented in the April 1984 release of the general purpose finite element program NASTRAN. The dynamic response problem is addressed in terms of the normal modal coordinates of these tuned rotating cyclic structures. Both rigid and flexible hubs/disks are considered. Coriolis and centripetal accelerations, as well as differential stiffness effects are included. Generally nonuniform steady inflow fields and uniform flow fields arbitrarily inclined at small angles with respect to the axis of rotation of the turbosystem are considered as the sources of aerodynamic excitation. The spatial nonuniformities are considered to be small deviations from a principally uniform inflow. Subsonic relative inflows are addressed, with provision for linearly interpolating transonic airloads.

Elchuri, V.↗