Search NASA⌕ Search

SEARCH · Search NASA

Results for “Topographic correction”

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 55 records · Page 3

Preprocessing: Geocoding of AVIRIS data using navigation, engineering, DEM, and radar tracking system data

Remotely sensed data have geometric characteristics and representation which depend on the type of the acquisition system used. To correlate such data over large regions with other real world representation tools like conventional maps or Geographic Information System (GIS) for verification purposes, or for further treatment within different data sets, a coregistration has to be performed. In addition to the geometric characteristics of the sensor there are two other dominating factors which affect the geometry: the stability of the platform and the topography. There are two basic approaches for a geometric correction on a pixel-by-pixel basis: (1) A parametric approach using the location of the airplane and inertial navigation system data to simulate the observation geometry; and (2) a non-parametric approach using tie points or ground control points. It is well known that the non-parametric approach is not reliable enough for the unstable flight conditions of airborne systems, and is not satisfying in areas with significant topography, e.g. mountains and hills. The present work describes a parametric preprocessing procedure which corrects effects of flight line and attitude variation as well as topographic influences and is described in more detail by Meyer.

Meyer, Peter↗

Characteristics of Forests in Western Sayani Mountains, Siberia from SAR Data

This paper investigated the possibility of using spaceborne radar data to map forest types and logging in the mountainous Western Sayani area in Siberia. L and C band HH, HV, and VV polarized images from the Shuttle Imaging Radar-C instrument were used in the study. Techniques to reduce topographic effects in the radar images were investigated. These included radiometric correction using illumination angle inferred from a digital elevation model, and reducing apparent effects of topography through band ratios. Forest classification was performed after terrain correction utilizing typical supervised techniques and principal component analyses. An ancillary data set of local elevations was also used to improve the forest classification. Map accuracy for each technique was estimated for training sites based on Russian forestry maps, satellite imagery and field measurements. The results indicate that it is necessary to correct for topography when attempting to classify forests in mountainous terrain. Radiometric correction based on a DEM (Digital Elevation Model) improved classification results but required reducing the SAR (Synthetic Aperture Radar) resolution to match the DEM. Using ratios of SAR channels that include cross-polarization improved classification and

Ranson, K. Jon↗

LiDAR-Inertial Based Navigation and Mapping for Precision Landing

Future lander missions will travel to ambitious, scientifically interesting locations near rough and dangerous terrain. They will need to operate with limited prior information about the terrain, and under varying lighting conditions. Landing safely and precisely in the face of these challenges is difficult for existing vision-based landing systems, which require detailed orbital reconnaissance, a priori hazard maps, and impose time-of-day restrictions on landing to ensure similar lighting conditions in orbital and descent imagery. Advanced 3D imaging LiDAR systems currently under development, and originally intended for single-scan hazard detection, have the potential to be operated continuously from altitudes of up to 5 km. Used together with existing inertial measurement units (IMUs), these sensors open a path-to-flight for a full navigation and mapping system, which could replace or augment a traditional landing sensor suite. A landing system based around these sensors can perform accurate altimetry, map-relative localization (MRL), LiDAR-inertial odometry, and map refinement in an illumination-insensitive manner, over unknown or partially known terrain. This paper outlines preliminary work on a LiDAR-inertial landing system that: estimates the spacecraft trajectory during entry, descent, and landing (EDL); and maps the topography of the terrain below, for future use in hazard detection and avoidance. An incremental, factor graph based, smoothing approach is used to solve for the maximum a posteriori trajectory of spacecraft states. Integrated IMU measurements and features tracked in adjacent range and intensity images are used to estimate motion (LiDAR-inertial odometry). LiDAR scans are binned into motion-corrected digital elevation models (DEMs), which are matched to an existing orbital topographic map to provide absolute position information (MRL). The estimated trajectory is then used to project the LiDAR scans into the map frame, creating a variable-resolution quadtree topographic map suitable for hazard detection and avoidance. Existing topographic maps from throughout the solar system (i.e., Earth, the Moon, Mars, Ceres, Vesta, Europa, Enceladus, and Eros) are upsampled for use in EDL simulations. The Mars 2020 Lander Vision System Simulator (LVSS) is extended to simulate LiDAR-inertial data for realistic EDL trajectories. Results of the algorithm operating on the simulated data are presented. Estimated spacecraft trajectory and refined map are compared to ground truth to assess estimation accuracy.

Katake, Anup↗

Topographic mapping using a monopulse SAR system

Terrain height variations in mountainous areas cause two problems in the radiometric correction of SAR images: the first being that the wrong elevation angle may be used in correcting for the radiometric variation of the antenna pattern; the second that the local incidence angle used in correcting the projection of the pixel area from slant range to ground range coordinates may vary from that given by the flat earth assumption. We propose a novel design of a SAR system which exploits the monopulse principle to determine the elevation angle and thus the height at the different parts of the image. The key element of such a phase monopulse system is an antenna, which can be divided into a lower and upper half in elevation using a monopulse comparator. In addition to the usual sum pattern, the elevation difference pattern can be generated by a -pi phase shift on one half of the antenna. From the ratios of images radiometrically modulated by the difference and sum antenna pattern in cross-track direction, we can derive the appropriate elevation angle at any point in the image. Together with the slant range we can calculate the height of the platform above this point using information on the antenna pointing and the platform attitude. This operation, repeated at many locations throughout the image, allows us to build up a topographic map of the height of the aircraft above each location. Inversion of this map, using the precisely determined aircraft altitude and the accurate flight path, leads to the actual topography of the imaged surface. The precise elevation of one point in the image could also be used to convert the height map to a topographic map. In this paper, we present design considerations for a corresponding airborne SAR system in X-Band and give estimates of the error due to system noise and azimuth ambiguities as well as the expected performance and precision in topographic mapping.

Zink, M.↗

The regression intersection method of adjusting image data for band ratioing

Estimation of the combined path radiance and sensor offset terms is essential in adjusting multispectral image radiance measurements for band ratioing. Commonly applied techniques for making this estimate have required assumptions or ancillary information regarding the reflectance properties of surface materials. This paper presents a technique that is unique in that it provides absolute (not just relative) statistically-derived estimates without the use of ancillary information. It is termed the regression intersection method (RIM). RIM is based on contrasts between the spectral properties of various homogeneous areas in rugged terrain. These areas are selected by examination of the image data alone. Bispectral regression (first principal component) lines are determined for each area and are projected, in pairs, to intersection points. Ideally, the coordinates of these points must equal the measurements for zero ground radiance since that is the only condition under which spectrally different materials can have the same radiance values. The median result from several site-pair and band-pair comparisons is used in order to statistically mitigate noise and minor variations due to natural variability. Tests show that the method is successful in determining correction values for the image data that result in maximum removal of the topographic effect in ratio images.

Crippen, Robert E.↗

Topographic effect on tilt, strain, and displacement measurements

An approximate analytical method is used to give simple expressions for the effect of arbitrarily shaped topography on vertical, horizontal and surface tilt measurements. Corrections to tilt fields are found to be directly proportional to the local slope of the topography. Comparison is made between expressions for tilt and those for strain and displacement. Although the topographic effects on secular tilt, strain, and vertical displacement are not yet measurable, the topographic effect on precision tidal tilt and strain measurements can be significant.

Meertens, Charles M.↗

A field measure of the shade fraction

'Shade' has a technical definition peculiar to linear spectral mixture analysis of imaging spectrometer data: it is the reduction in radiance from a surface due to lighting conditions and geometry, and includes topographic shading described by photometric functions as well as shadowing at all scales. 'Shade' is an important constituent of nearly all remotely sensed images, and is one endmember resolved in spectral mixture analysis, where it is represented as a fraction of the measured radiance and a characteristic spectrum. This spectrum is typically the null vector, provided the data have been corrected for atmospheric and instrument effects: i.e., 'shade' is the radiance from an ideal black surface. In topographic shading, irradiance is reduced - typically in proportion to cos(i), where i (incidence angle) is the angle between the sun and the local surface normal vectors. Therefore, the radiance is lowered by a multiplicative factor. Shadowing occurs when i is greater than 90 deg, or when sunlight is blocked by adjacent high terrain; the only irradiance is down-welling skylight and bounce light from adjacent terrain. In spectral mixture analysis, 'shade' is regarded as an additive term. In this regard, it is an accurate description of the proportion of a scene that consists of ideal shadows ('checkerboard mixing'); however, 'shade' represents the multiplicative cos(i) factor as well, as here it should be interpreted as the proportion of shadow that would darken the scene an equivalent amount. In either case, the 'shade' fraction is lessened by adjacency effects, because the scene has a non-zero reflectivity instead of the ideal black surface generally assumed.

Gillespie, Alan R.↗

The reduction correction in North America

An inverse Poisson integral technique was used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1 x 1 deg, 5 x 5', 5/8 x 5/8') from a decreasing area of integration. Reduction corrections are computed at 5/8 x 5/8' granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5' intervals. The 1 x 1 deg averages of the reduction corrections thus obtained for 172 1 x 1 deg squares in western North America are discussed. The 1 x 1 deg mean reduction corrections are predominantly positive, varying from -3 to +15 mgal, with values in excess of 5 mgal for 26 squares. Their mean and rms values are +2.4 and 3.6 mgal respectively and they correlate well with the mean terrain corrections. The mean and rms contributions from the three stages of computation are: 1 x 1 deg stage +0.15 and 0.7 mgal; 5 x 5' stage + 1.0 and 1.6 mgal; and 5/8 x 5/8' stage +1.3 and 1.8 mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.

Martzen, P. D.↗

The Collection 6 'dark-target' MODIS Aerosol Products

Aerosol retrieval algorithms are applied to Moderate resolution Imaging Spectroradiometer (MODIS) sensors on both Terra and Aqua, creating two streams of decade-plus aerosol information. Products of aerosol optical depth (AOD) and aerosol size are used for many applications, but the primary concern is that these global products are comprehensive and consistent enough for use in climate studies. One of our major customers is the international modeling comparison study known as AEROCOM, which relies on the MODIS data as a benchmark. In order to keep up with the needs of AEROCOM and other MODIS data users, while utilizing new science and tools, we have improved the algorithms and products. The code, and the associated products, will be known as Collection 6 (C6). While not a major overhaul from the previous Collection 5 (C5) version, there are enough changes that there are significant impacts to the products and their interpretation. In its entirety, the C6 algorithm is comprised of three sub-algorithms for retrieving aerosol properties over different surfaces: These include the dark-target DT algorithms to retrieve over (1) ocean and (2) vegetated-dark-soiled land, plus the (3) Deep Blue (DB) algorithm, originally developed to retrieve over desert-arid land. Focusing on the two DT algorithms, we have updated assumptions for central wavelengths, Rayleigh optical depths and gas (H2O, O3, CO2, etc.) absorption corrections, while relaxing the solar zenith angle limit (up to 84) to increase pole-ward coverage. For DT-land, we have updated the cloud mask to allow heavy smoke retrievals, fine-tuned the assignments for aerosol type as function of season location, corrected bugs in the Quality Assurance (QA) logic, and added diagnostic parameters such as topographic altitude. For DT-ocean, improvements include a revised cloud mask for thin-cirrus detection, inclusion of wind speed dependence in the retrieval, updates to logic of QA Confidence flag (QAC) assignment, and additions of important diagnostic information. At the same time as we have introduced algorithm changes, we have also accounted for upstream changes including: new instrument calibration, revised land-sea masking, and changed cloud masking. Upstream changes also impact the coverage and global statistics of the retrieved AOD. Although our responsibility is to the DT code and products, we have also added a product that merges DT and DB product over semi-arid land surfaces to provide a more gap-free dataset, primarily for visualization purposes. Preliminary validation shows that compared to surface-based sunphotometer data, the C6, Level 2 (along swath) DT-products compare at least as well as those from C5. C6 will include new diagnostic information about clouds in the aerosol field, including an aerosol cloud mask at 500 m resolution, and calculations of the distance to the nearest cloud from clear pixels. Finally, we have revised the strategy for aggregating and averaging the Level 2 (swath) data to become Level 3 (gridded) data. All together, the changes to the DT algorithms will result in reduced global AOD (by 0.02) over ocean and increased AOD (by 0.02) over land, along with changes in spatial coverage. Changes in calibration will have more impact to Terras time series, especially over land. This will result in a significant reduction in artificial differences in the Terra and Aqua datasets, and will stabilize the MODIS data as a target for AEROCOM studie

Aerosol retrieval algorithms↗

Measuring Directional Wave Spectra and Wind Speed with a Scanning Radar Altimeter

The geometry for the NASA Scanning Radar Altimeter (SRA) is shown. It transmits a 8-ns duration pulse at Ka-band (8.3 mm) and measures time of flight as it scans a 1 degree (two-way) beam from left to right across the aircraft ground track. The most recent configuration determines the surface elevation at 64 points spaced at uniform angular intervals of about 0.7 across a swath whose width is about 0.8 times the aircraft altitude. The system generates these raster lines of the surface topography beneath the aircraft at about a 10 Hz rate. In postflight processing the SRA wave topographic data are transformed with a two-dimensional Fast Fourier Transformation (FFT) and Doppler corrected to produce directional wave spectra. The SRA is not absolutely calibrated in power, but by measuring the relative fall-off of backscatter with increasing incidence angle, the SRA can also determine the mean square slope (mss) of the sea surface, a surrogate for wind speed. For the slope-dependent specular point model of radar sea surface scattering, an expression approximated by a geometric optics form, for the relative variation with incidence angle of the normalized backscatter radar cross section would be sigma (sup 0) (sub rel) = sec (exp 4) theta exp (-tan squared theta/mss) where theta is the off-nadir incidence angle.

Walsh, E. J.↗

Application of Core Dynamics Modeling to Core-Mantle Interactions

Observations have demonstrated that length of day (LOD) variation on decadal time scales results from exchange of axial angular momentum between the solid mantle and the core. There are in general four core-mantle interaction mechanisms that couple the core and the mantle. Of which, three have been suggested likely the dominant coupling mechanism for the decadal core-mantle angular momentum exchange, namely, gravitational core-mantle coupling arising from density anomalies in the mantle and in the core (including the inner core), the electromagnetic coupling arising from Lorentz force in the electrically conducting lower mantle (e.g. D-layer), and the topographic coupling arising from non-hydrostatic pressure acting on the core-mantle boundary (CMB) topography. In the past decades, most effort has been on estimating the coupling torques from surface geomagnetic observations (kinematic approach), which has provided insights on the core dynamical processes. In the meantime, it also creates questions and concerns on approximations in the studies that may invalidate the corresponding conclusions. The most serious problem is perhaps the approximations that are inconsistent with dynamical processes in the core, such as inconsistencies between the core surface flow beneath the CMB and the CMB topography, and that between the D-layer electric conductivity and the approximations on toroidal field at the CMB. These inconsistencies can only be addressed with numerical core dynamics modeling. In the past few years, we applied our MoSST (Modular, Scalable, Self-consistent and Three-dimensional) core dynamics model to study core-mantle interactions together with geodynamo simulation, aiming at assessing the effect of the dynamical inconsistencies in the kinematic studies on core-mantle coupling torques. We focus on topographic and electromagnetic core-mantle couplings and find that, for the topographic coupling, the consistency between the core flow and the CMB topography is critical for correct evaluation of the coupling torque.

Kuang, Weijia↗

Planetary Geophysics and Tectonics

We have carried out several studies that explore explanations for the role of chemical density variations in Moon s evolution. Meaningful models for the evolution of the Moon must explain a number of important magmatic characteristics. Volcanic activity subsequent to the formation of its anorthositic crust was dominated by the eruption of mare basalt. 1) The main phase of mare volcanism began approx. 500 Myr after the crystallization of the anorthositic crust and continued for approx. l Gyr. 2) The picitic glasses, considered to be representative of mare basalt least affected by low pressure, near-surface fractionation, were generated by melting, at 400-600 km depth, of a source containing components that, on the basis of the magma ocean hypothesis, should have crystallized at much shallower depth during fractionation of the anorthositic crust. 3) Mare basalts occur primarily in one region of the Moon. Recent topographic data demonstrate that the earlier idea that mare basalt flooded areas of low elevation is not correct. Large areas of very low elevation do not contain mare basalt. The hemispheric asymmetry of mare basalt distribution on the lunar surface must be explained in some other way. 4) A region of the surface roughly correlating with that containing mare basalts also is thought to contain high subsurface concentrations of KREEP which was excavated during the formation of large impact basins. This so-called Procellarum KREEP Terrane (PKT) is responsible for the Imbrium basin-centered thorium anomaly mapped by Lunar Prospector.

Parmentier, E. Marc↗

Modeling the effects of various radiant transfers in mountainous terrain on sensor response

A model is developed for the study of remote sensing spectral data problems encountered in mountainous terrain, where slope orientation with respect to sources of incident radiation must be considered for accurate classification of multispectral data. The model simulates the effects of solar, sky and adjacent slope irradiance fluxes and the bidirectional reflectance properties of the target surface on sensor response for various slope orientations. The errors associated with two practical techniques used to account for topographic effects on sensor response were evaluated in light of the model, and it was found that the simple cosine correction of the direct solar source showed reflectance errors greater than 11% while the technique considering both direct solar and diffuse sky sources was 50% more accurate.

Kimes, D. S.↗

Geocoding of AIRSAR/TOPSAR SAR Data

It has been demonstrated that radar interferometry is a promising method for determination of digital elevation information and terrain slope from synthetic aperture radar (SAR) data. A multipolarimetric radar AIRSAR operates in the P, L, and C bands on board the NASA DC-8 aircraft. The TOPSAR, a special mode of the AIRSAR system, is able to collect single pass interferometric C and/or L band VV polarized data. A possible configuration of the AIRSAR/TOPSAR system is to acquire single pass interferometric data at C-band VV polarization and polarimetric radar data at the two other lower frequencies. The advantage of this configuration is to acquire digital topographic information at the same time the radar data is collected. The digital elevation information can therefore be used to correctly calibrate the SAR data. In this paper, a procedure to geocode the new AIRSAR/TOPSAR data is presented and an earlier AIRSAR/TOPSAR image is geocoded and evaluated in terms of geometric accuracy.

Holecz, Francesco↗

Topography-Dependent Motion Compensation: Application to UAVSAR Data

The UAVSAR L-band synthetic aperture radar system has been designed for repeat track interferometry in support of Earth science applications that require high-precision measurements of small surface deformations over timescales from hours to years. Conventional motion compensation algorithms, which are based upon assumptions of a narrow beam and flat terrain, yield unacceptably large errors in areas with even moderate topographic relief, i.e., in most areas of interest. This often limits the ability to achieve sub-centimeter surface change detection over significant portions of an acquired scene. To reduce this source of error in the interferometric phase, we have implemented an advanced motion compensation algorithm that corrects for the scene topography and radar beam width. Here we discuss the algorithm used, its implementation in the UAVSAR data processor, and the improvement in interferometric phase and correlation achieved in areas with significant topographic relief.

surface deformation↗

Computer-aided analysis of Skylab multispectral scanner data in mountainous terrain for land use, forestry, water resource, and geologic applications

The author has identified the following significant results. One of the most significant results of this Skylab research involved the geometric correction and overlay of the Skylab multispectral scanner data with the LANDSAT multispectral scanner data, and also with a set of topographic data, including elevation, slope, and aspect. The Skylab S192 multispectral scanner data had distinct differences in noise level of the data in the various wavelength bands. Results of the temporal evaluation of the SL-2 and SL-3 photography were found to be particularly important for proper interpretation of the computer-aided analysis of the SL-2 and SL-3 multispectral scanner data. There was a quality problem involving the ringing effect introduced by digital filtering. The modified clustering technique was found valuable when working with multispectral scanner data involving many wavelength bands and covering large geographic areas. Analysis of the SL-2 scanner data involved classification of major cover types and also forest cover types. Comparison of the results obtained wth Skylab MSS data and LANDSAT MSS data indicated that the improved spectral resolution of the Skylab scanner system enabled a higher classification accuracy to be obtained for forest cover types, although the classification performance for major cover types was not significantly different.

Hoffer, R. M.↗

LANDSAT-D investigations in snow hydrology

Progress on the registration of TM data to digital topographic data; on comparison of TM, MSS and NOAA meteorological satellite data for snowcover mapping; and on radiative transfer models for atmospheric correction is reported. Some methods for analyzing spatial contiguity of snow within the snow covered area were selected. The methods are based on a two-channel version of the grey level co-occurence matrix, combined with edge detection derived from an algorithm for computing slopes and exposures from digital terrain data.

Dozier, J.↗

Asteroid approach covariance analysis for the Clementine mission

The Clementine mission is designed to test Strategic Defense Initiative Organization (SDIO) technology, the Brilliant Pebbles and Brilliant Eyes sensors, by mapping the moon surface and flying by the asteroid Geographos. The capability of two of the instruments available on board the spacecraft, the lidar (laser radar) and the UV/Visible camera is used in the covariance analysis to obtain the spacecraft delivery uncertainties at the asteroid. These uncertainties are due primarily to asteroid ephemeris uncertainties. On board optical navigation reduces the uncertainty in the knowledge of the spacecraft position in the direction perpendicular to the incoming asymptote to a one-sigma value of under 1 km, at the closest approach distance of 100 km. The uncertainty in the knowledge of the encounter time is about 0.1 seconds for a flyby velocity of 10.85 km/s. The magnitude of these uncertainties is due largely to Center Finding Errors (CFE). These systematic errors represent the accuracy expected in locating the center of the asteroid in the optical navigation images, in the absence of a topographic model for the asteroid. The direction of the incoming asymptote cannot be estimated accurately until minutes before the asteroid flyby, and correcting for it would require autonomous navigation. Orbit determination errors dominate over maneuver execution errors, and the final delivery accuracy attained is basically the orbit determination uncertainty before the final maneuver.

Ionasescu, Rodica↗