Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling algorithm”

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 559 records · Page 31

Landsat D Thematic Mapper image resampling for scan geometry correction

The Landsat D system is described, and the geometric correction processing for the Thematic Mapper (TM) is reviewed. The resampling procedure (the generation of a TM output image) is analyzed with emphasis on the effect of sampling geometry on output image. Effects of scan gaps and spacecraft jitter on output image are studied by means of a simulation of the sampling and resampling processes for three sampling geometries. Visual inspection of the resampling results shows that the resampling algorithm works excellently under all conditions and that distortion is visible only in the rare case of large gaps between scans.

Prakash, A.↗

Estimation of coefficients and boundary parameters in hyperbolic systems

Semi-discrete Galerkin approximation schemes are considered in connection with inverse problems for the estimation of spatially varying coefficients and boundary condition parameters in second order hyperbolic systems typical of those arising in 1-D surface seismic problems. Spline based algorithms are proposed for which theoretical convergence results along with a representative sample of numerical findings are given.

Banks, H. T.↗

Design of the Galileo scan platform control

Design of the scan platform control for the Galileo spacecraft is presented. Emphasis is placed on the inertial pointing of the scan platform when the spacecraft is in the dual-spin configuration. The various methods of operation used for Galileo scan pointing are outlined. Major design considerations, such as spacecraft flexibility and separation of actuator and sensor by a flexible structure, are discussed. The pointing requirements imposed on the scan platform control are explained. A high-level description of the relevant scan pointing algorithms is included. The performance of the design is demonstrated by a sample slew test case.

Chodas, J. L.↗

Nonlinear filtering for spacecraft attitude estimation

Nonlinear filtering techniques are applied to spacecraft attitude estimation using quaternion parameterization for the attitude kinematics. By replacing the angular velocity vector by the gyro output vector, a state dependent noise vector is introduced in the seven-dimensional system equations. The resulting conditional probability density function from the Ito differential rule is governed by the Fokker Planck partial differential equation which is approximated by the second order mean and covariance differential equations. In order to minimize computer loading, the covariance propagation is carried out in six-dimensional state space using a matrix transformation. The star tracker data is used to update the covariance matrix in the seven-dimensional space. The algorithm is simulated for an earth pointing spacecraft mission, using Monte Carlo samples of gyro and star measurements. The performance of the second order filter is compared with the extended Kalman Filter through several simulation runs and drift rates have been identified.

Vathsal, S.↗

Implicit finite-difference methods for the Euler equations

The present paper is concerned with two-dimensional Euler equations and with schemes which are in use of the time of this writing. Most of the development presented carries over directly to three dimensions. The characteristics of the two-dimensional Euler equations in Cartesian coordinates are considered along with generalized curvilinear coordinate transformations, metric relations, invariants of the transformation, flux Jacobian matrices and eigensystems, numerical algorithms, flux split algorithms, implicit and explicit nonlinear control (smoothing), upwind differencing in supersonic regions, unsteady and steady-state computation, the diagonal form of implicit algorithm, metric differencing and invariants, boundary conditions, geometry and mesh generation, and sample solutions.

Pulliam, T. H.↗

Estimation of coefficients and boundary parameters in hyperbolic systems

Semi-discrete Galerkin approximation schemes in connection with inverse problems for the estimation of spatially varying coefficients and boundary condition parameters in second order hyperbolic systems typical of those arising in 1-D surface seismic problems are considered. Spline based algorithms are proposed for which theoretical convergence results along with a representative sample of numerical findings are given.

Banks, H. T.↗

Adaptive numerical solutions of the Euler equations in 3D using finite elements

The development of an adaptive mesh solution for a flow involving shock interaction on a swept cylinder and an initial solution for a flow past a complex fighter configuration is reported. The finite element solution algorithm, the mesh generation, and the adaptivity of the solution are described. Sample results for the flow past an F-18 configuration at Mach 0.9 and alpha of 3 deg and for shock interaction on a swept cylinder at Mach 8.04 are summarized.

Peraire, J.↗

Finding fixed satellite service orbital allotments with a k-permutation algorithm

A satellite system synthesis problem, the satellite location problem (SLP), is addressed. In SLP, orbital locations (longitudes) are allotted to geostationary satellites in the fixed satellite service. A linear mixed-integer programming model is presented that views SLP as a combination of two problems: the problem of ordering the satellites and the problem of locating the satellites given some ordering. A special-purpose heuristic procedure, a k-permutation algorithm, has been developed to find solutions to SLPs. Solutions to small sample problems are presented and analyzed on the basis of calculated interferences.

Reilly, Charles H.↗

Faint X-ray source counts and the origin of the X-ray background

A thorough reanalysis of the Einstein Observatory Deep Survey data is presented in order to determine the contribution of detected, discrete sources of X-ray emission to the cosmic X-ray background. Substantial discrepancies with previously published work on this problem are found. A detailed discussion of data editing and source algorithms buttresses a claim of having constructed a complete, flux-limited sample of the faintest sources detectable with the Einstein imaging proportional counter, the most sensitive X-ray instrument yet flown. A total of 33 sources is found in a survey region of about 3.3 sq deg down to a minimum flux threshold of 4 x 10 to the 14th ergs/sq cm/s in the 0.3-3.5 keV band. Roughly 30 percent of the objects are foreground stars, leading to an extragalactic source surface density of 70,000/sr at this threshold. The integrated contribution from discrete sources to the number of cosmic X-ray background photons measured in this same band with the same instrument is 12 percent + or - 3 percent, substantially below previous estimates. Implications of these results for the origin of the background are discussed.

Hamilton, T. T.↗

Using single buffers and data reorganization to implement a multi-megasample fast Fourier transform

Data ordering in large fast Fourier transforms (FFT's) is both conceptually and implementationally difficult. Discribed here is a method of visualizing data orderings as vectors of address bits, which enables the engineer to use more efficient data orderings and reduce double-buffer memory designs. Also detailed are the difficulties and algorithmic solutions involved in FFT lengths up to 4 megasamples (Msamples) and sample rates up to 80 MHz.

Brown, R. D.↗

Parallel Computation of Unsteady Flows on a Network of Workstations

Parallel computation of unsteady flows requires significant computational resources. The utilization of a network of workstations seems an efficient solution to the problem where large problems can be treated at a reasonable cost. This approach requires the solution of several problems: 1) the partitioning and distribution of the problem over a network of workstation, 2) efficient communication tools, 3) managing the system efficiently for a given problem. Of course, there is the question of the efficiency of any given numerical algorithm to such a computing system. NPARC code was chosen as a sample for the application. For the explicit version of the NPARC code both two- and three-dimensional problems were studied. Again both steady and unsteady problems were investigated. The issues studied as a part of the research program were: 1) how to distribute the data between the workstations, 2) how to compute and how to communicate at each node efficiently, 3) how to balance the load distribution. In the following, a summary of these activities is presented. Details of the work have been presented and published as referenced.

Source record↗

In-Flight Edge Response Measurements for High Spatial Resolution Remote Sensing Systems

In-flight measurement of spatial resolution were conducted as part of the ASA Scientific Data Purchase (SDP) Validation and Verification (V&V) process. Characterization included remote sensing systems with ground sample distance (GSD) of 1 meter or less, such as the panchromatic imager on-board the ICONOS satellite and the airborne ADAR System 5500 multispectral instrument. Final image products were used to evaluate the effect of both the image acquisition system (e.g., optics, electronics, motion, jitter, atmosphere) and image post-processing (e.g., resampling, modulation trasfer function (MTF) compensator). Spatial resolution was characterized by full width at half maximum (FWHM) of an edge response-derived line spread function. This was found to be a more robust measure of spatial resolution than the value of NTF at Nyquist frequency The edge responses were analysed using the tilted-edge technique that ovecomes the spatial sampling limitations of the digital imaging systems. As an enhancement to existing algorithms, the slope of the edge response and the orientation of the edge target were determined by a single computational process. Adjacent black and white square panels, either painted on a flat surface or deployed as traps, formed the ground-based edge targets used in the tests. Orientation of the deployable tarps was optimized beforehand, based on simulations of the imaging system. Numerous edge target images were analyzed for each of the tested sensors. The effect of such factors as acquisition geometry, temporal variability, MTF compensation, and GSD on spatial resolution were investigated.

Blonski, Slawomir↗

MISR Level 1A CCD Science data, all cameras (MIL1A_V1)

The Level 1A data are raw MISR data that are decommutated, reformatted 12-bit Level 0 data shifted to byte boundaries, i.e., reversal of square-root encoding applied and converted to 16 bit, and annotated (e.g., with time information). These data are used by the Level 1B1 processing algorithm to generate calibrated radiances. The science data output preserves the spatial sampling rate of the Level 0 raw MISR CCD science data. CCD data are collected during routine science observations of the sunlit portion of the Earth. Each product represents one 'granule' of data. A 'granule' is defined to be the smallest unit of data required for MISR processing. Also, included in the Level 1A product are pointers to calibration coefficient files provided for Level 1B processing. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180].

CCD↗

MISR Level 1A CCD Science data, all cameras (MIL1A_V2)

The Level 1A data are raw MISR data that are decommutated, reformatted 12-bit Level 0 data shifted to byte boundaries, i.e., reversal of square-root encoding applied and converted to 16 bit, and annotated (e.g., with time information). These data are used by the Level 1B1 processing algorithm to generate calibrated radiances. The science data output preserves the spatial sampling rate of the Level 0 raw MISR CCD science data. CCD data are collected during routine science observations of the sunlit portion of the Earth. Each product represents one 'granule' of data. A 'granule' is defined to be the smallest unit of data required for MISR processing. Also, included in the Level 1A product are pointers to calibration coefficient files provided for Level 1B processing. [Location=GLOBAL] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180].

AM-1↗

Modeling the Height of Young Forests Regenerating from Recent Disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat-lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, "young" forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS-VCT disturbance products. The height of "young" forest is then mapped based on the derived relationships and the LTSS-VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS-VCT method yielded the lowest cross validation error. The R(exp 2) and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS-VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Li, Ainong↗

Detection of Spectral Evolution in the Bursts Emitted During the 2008-2009 Active Episode of SGR J1550 - 5418

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data of all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations 170ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to less than or equal to 8 degrees. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

von Kienlin, Andreas↗

Uncertainty-aware Continuous Implicit Neural Representations for Remote Sensing Object Counting

Many existing object counting methods rely on density map estimation (DME) of the discrete grid representation by decoding extracted image semantic features from designed convolutional neural networks (CNNs). Relying on discrete density maps not only leads to information loss dependent on the original image resolution, but also has a scalability issue when analyzing high-resolution images with cubically increasing memory complexity. Furthermore, none of the existing methods can offer reliable uncertainty quantification (UQ) for the derived count estimates. To overcome these limitations, we design UNcertainty-aware, hypernetwork-based Implicit neural representations for Counting (UNIC) to assign probabilities and the corresponding counting confidence over continuous spatial coordinates. We derive a sampling-based Bayesian counting loss function and develop the corresponding model training algorithm. UNIC outperforms existing methods on the Remote Sensing Object Counting (RSOC) dataset with reliable UQ and improved interpretability of the derived count estimates. Our code is available at https://github.com/SiyuanXu-tamu/UNIC.

97 MATHEMATICS AND COMPUTING↗

Maximum a posteriori Ly α estimator (MAPLE): band power and covariance estimation of the 3D Ly α forest power spectrum

We present a novel maximum a posteriori estimator to jointly estimate band powers and the covariance of the three-dimensional power spectrum (P3D) of Ly $\alpha$ forest flux fluctuations, called MAPLE. Our Wiener-filter based algorithm reconstructs a window-deconvolved P3D in the presence of complex survey geometries typical for Ly $\alpha$ surveys that are sparsely sampled transverse to and densely sampled along the line of sight. We demonstrate our method on idealized Gaussian random fields with two selection functions: (i) a sparse sampling of 30 background sources per square degree designed to emulate the current Dark Energy Spectroscopic Instrument; (ii) a dense sampling of 900 background sources per square degree emulating the upcoming Prime Focus Spectrograph Galaxy Evolution Survey. Our proof-of-principle shows promise, especially since the algorithm can be extended to marginalize jointly over nuisance parameters and contaminants, i.e. offsets introduced by continuum fitting. Our code is implemented in JAX and is publicly available on GitHub.

79 ASTRONOMY AND ASTROPHYSICS↗