Kirk, R. L.; Callahan, P.; Seu, R.; Lorenz, R. D.; Paganelli, F.; Lopes, R.; Elachi, C.
The Cassini Titan RADAR Mapper is a K(sub u)-band (13.78 GHz, lambda = 2.17 cm) linear polarized RADAR instrument capable of operating in synthetic aperture (SAR), scatterometer, altimeter and radiometer modes. During the first targeted flyby of Titan on 26 October, 2004 (referred to as Ta) observations were made in all modes. Evidence for topographic relief based on the Ta altimetry and SAR data are presented here. Additional SAR and altimetry observations are planned for the T3 encounter on 15 February, 2005, but have not been carried out at this writing. Results from the T3 encounter relevant to topography will be included in our presentation. Data obtained in the Ta encounter include a SAR image swath
Farr, T. G.; Kobrick, M.
The Shuttle Radar Topography Mission (SRTM), which flew successfully aboard Endeavour in February 2000, is a cooperative project between NASA and the National Imagery and Mapping Agency (NIMA). The mission was designed to use a single-pass radar interferometer to produce a digital elevation model of the Earth's land surface between about 60 degrees north and 56 degrees south latitude. The DEM will have 30 m horizontal resolution and about 15 m vertical errors. Two ortho-rectified C-band image mosaics are also planned. SRTM used a modification of the radar instrument that comprised the Spaceborne Radar Laboratory that flew twice on the Shuttle Endeavour in 1994. To collect the interferometric data, a 60 m mast, additional C-band antenna, and improved tracking and navigation devices were added. A second X-band antenna was also added by the German Space Agency, and produced higher resolution topographic measurements in strips nested within the full, C-band coverage. First results indicate that the radars and ancillary instruments worked very well. Data played back to the ground during the flight were processed to DEMs and products released hours after acquisition. An extensive program for calibration and verification of the SRTM data is now underway. When complete later this year, systematic processing of the data will begin, with final products emerging a continent at a time. Data processing will be completed by the end of 2002. Products will be transferred to the US Geological Survey's EROS Data Center for civilian archive and distribution. NIMA will handle Department of Defense distribution. * Work performed under contract to NASA.
Farr, T. G.; Kobrick, M.
The Shuttle Radar Topography Mission (SRTM), which flew successfully aboard Endeavour in February 2000, is a cooperative project between NASA, the National Imagery and Mapping Agency, and the German and Italian Space Agencies. The mission was designed to use a single-pass radar interferometer to produce a digital elevation model of the Earth's land surface between about 60 degrees north and 56 degrees south latitude. The DEM will have 30 m horizontal resolution and better than 15 m vertical errors. Two ortho-rectified C-band image mosaics are also planned. Data processing will be completed by the end of 2002. SRTM used a modification of the radar instrument that comprised the Spaceborne Radar Laboratory that flew twice on the Shuttle Endeavour in 1994. To collect the interferometric data, a 60 m mast, additional C-band antenna, and improved tracking and navigation devices were added. A second X-band antenna was also added by the German Space Agency, and produced higher resolution topographic measurements in strips nested within the full, C-band coverage. First results indicate that the radars and ancillary instruments worked very well. Data played back to the ground during the flight were processed to DEMs and products released hours after acquisition. An extensive program for calibration and verification of the SRTM data is now underway. When complete later this year, systematic processing of the data will begin, with final products emerging a continent at a time. Products will be transferred to the US Geological Survey's EROS Data Center for civilian archive and distribution. NIMA will handle Department of Defense distribution. * Work performed under contract to NASA.
U.S. Geological Survey
Under an agreement with the National Aeronautics and Space Administration (NASA) and the Department of Defense's National Geospatial-Intelligence Agency (NGA), the U.S. Geological Survey (USGS) is distributing elevation data from the Shuttle Radar Topography Mission (SRTM). The SRTM is a joint project of NASA and NGA to map the Earth's land surface in three dimensions at an unprecedented level of detail. As part of space shuttle Endeavour's flight during February 11-22, 2000, the SRTM successfully collected data over 80 percent of the Earth's land surface for most of the area between latitudes 60 degrees north and 56 degrees south. The SRTM hardware included the Spaceborne Imaging Radar-C (SIR-C) and X-band Synthetic Aperture Radar (X-SAR) systems that had flown twice previously on other space shuttle missions. The SRTM data were collected with a technique known as interferometry that allows image data from dual radar antennas to be processed for the extraction of ground heights.
U.S. Geological Survey
Under an agreement with the National Aeronautics and Space Administration (NASA) and the Department of Defense's National Imagery and Mapping Agency (NIMA), the U.S. Geological Survey (USGS) is now distributing elevation data from the Shuttle Radar Topography Mission (SRTM). The SRTM is a joint project between NASA and NIMA to map the Earth's land surface in three dimensions at a level of detail unprecedented for such a large area. Flown aboard the NASA Space Shuttle Endeavour February 11-22, 2000, the SRTM successfully collected data over 80 percent of the Earth's land surface, for most of the area between 60? N. and 56? S. latitude. The SRTM hardware included the Spaceborne Imaging Radar-C (SIR-C) and X-band Synthetic Aperture Radar (X-SAR) systems that had flown twice previously on other space shuttle missions. The SRTM data were collected specifically with a technique known as interferometry that allows image data from dual radar antennas to be processed for the extraction of ground heights.
Podest, E.; McDonald, K.; Kimball, J.; Randerson, J. T.
The annual freeze/thaw cycle drives the length of the growing season in the boreal forest, and is a major factor determining annual productivity and associated exchange of CO2 with the atmosphere. Variations in freeze/thaw processes are spatially and temporally complex in boreal environments, particularly in areas of complex topography and in fire disturbance regimes. We investigate the spatial and temporal characteristics of seasonal freeze/thaw dynamics in complex boreal landscapes, as derived from radar backscatter measured with ERS (C-band, VV polarization, 200m resolution) and JERS-1 (L-band, HH polarization, 100m resolution) Synthetic Aperture Radars (SARs), and with the SeaWinds scatterometer (Ku-band, 25km resolution). C- and L-band backscatter are applied to characterize freeze/thaw transitions for a chronosequence of recovering burn sites near Delta Junction, Alaska, and for a region of complex topography on the Kenai Peninsula, Alaska. We characterize differences in radar-derived freeze/thaw state, examining transitions over complex terrain and landscape disturbance regimes. In areas of complex terrain, we explore freeze/thaw dynamics related to elevation, slope aspect and varying landcover. In the burned regions, we explore the timing of seasonal freeze/thaw transition as related to the recovering landscape, relative to that of a nearby control site. We apply in situ biophysical measurements, including flux tower measurements to validate and interpret the remotely sensed parameters. A multi-scale analysis is performed relating high-resolution SAR backscatter and moderate resolution scatterometer measurements to assess trade-offs in spatial and temporal resolution in the remotely sensed fields. A temporal change discriminator is applied to classify time series radar imagery to classify the landscape freeze-thaw state. We apply a 30m-resolution digital elevation model (DEM) derived from Shuttle Radar Topography Mission (SRTM) data to orthorectify the time
U.S. Geological Survey
In February 2000, the Shuttle Radar Topography Mission (SRTM) successfully collected Interferometric C-Band Synthetic Aperture Radar data over 80 percent of the Earth's land surface, for most of the area between 60?N and 56?S latitude. NASA and the National Geospatial-Intelligence Agency (NGA), formerly known as the National Imagery and Mapping Agency (NIMA), co-sponsored the mission. NASA's Jet Propulsion Laboratory (JPL) performed preliminary processing of SRTM data and forwarded partially finished data directly to NGA for finishing by NGA contractors and subsequent monthly deliveries to the NGA Digital Products Data Warehouse (DPDW). All data products delivered by the contractors conform to NGA SRTM Data Products and NGA Digital Terrain Elevation Data? (DTED?) specifications. The DPDW ingests the SRTM data products, checks them for formatting errors, loads the public SRTM DTED? into the NGA data distribution system, and ships them to the U.S. Geological Survey (USGS) Center for Earth Resources Observation and Science (EROS). In addition to NGA's SRTM DTED? format, USGS EROS has reformatted the data into a non-proprietary, generic raster binary SRTM format that is readable by most remote sensing software packages. The SRTM format is also publicly available from USGS EROS.
Lin, Qian; Vesecky, John F.; Zebker, Howard A.
Methods are presented for using Synthetic Aperture Radar (SAR) interferometry data to estimate surface topography. An expression is given to relate the elevation of a ground point to the phase difference of SAR images received from two spatially separated antennas. An iterative algorithm which solves for the position and elevation of each point in the image simultaneously is developed. One of the critical issues that determines the accuracy of the terrain mapping is the phase unwrapping. An approach to the problem by fringe line detection is proposed. The algorithms are tested with two Seasat SAR images of terrain near Yellowstone National Park. The resultant elevation map is compared with a USGS terrain elevation model. The error of the SAR elevation with respect to the digital terrain map is about 8.2 percent of the total terrain variation.
Farr, T. G.; Werner, M.; Kobrick, M.
The Shuttle Radar Topography Mission (SRTM), which flew successfully aboard Endeavour in February 2000, is a cooperative project between NASA, the National Imagery and Mapping Agency, and the German and Italian Space Agencies. The mission was designed to use a single-pass radar interferometer to produce a digital elevation model of the Earth's land surface between about 60^o north and 56^o south latitude. The DEM has 30 m horizontal resolution and better than 15 m vertical errors. Two ortho-rectified C-band image mosaics are also produced. SRTM used a modification of the radar instrument that comprised the Spaceborne Radar Laboratory that flew twice on the Shuttle Endeavour in 1994. To collect the interferometric data, a 60 m mast, additional C-band antenna, and improved tracking and navigation devices were added. A second X-band antenna was also added by the German Space Agency, and produced higher resolution topographic measurements in strips nested within the full, C-band coverage. First results indicate that the radars and ancillary instruments worked very well. Data played back to the ground during the flight were processed to DEMs and preliminary products released hours after acquisition. Precision processing of the C-band data was completed at the end of 2002. An extensive program for calibration and verification of the SRTM data is now underway. Data have been released so far for the US and a few test areas for scientific analysis. Public release of the data will occur in stages throughout 2003. Products are being transferred to the US Geological Survey's EROS Data Center for civilian archive and distribution. NIMA will handle Department of Defense distribution. X-band data are being processed at the German and Italian Space Agencies. As of late 2002, Europe and Africa had been completed and the remaining continents were on schedule to be completed by the end of 2003. This special session will highlight applications of this new high-resolution view of the
Molchanov, Pavlo A.
All digital radar architecture requires exclude mechanical scan system. The phase antenna array is necessarily large because the array elements must be co-located with very precise dimensions and will need high accuracy phase processing system for aggregate and distribute T/R modules data to/from antenna elements. Even phase array cannot provide wide field of view. New nature inspired all digital radar architecture proposed. The fly's eye consists of multiple angularly spaced sensors giving the fly simultaneously thee wide-area visual coverage it needs to detect and avoid the threats around him. Fly eye radar antenna array consist multiple directional antennas loose distributed along perimeter of ground vehicle or aircraft and coupled with receiving/transmitting front end modules connected by digital interface to central processor. Non-steering antenna array allows creating all-digital radar with extreme flexible architecture. Fly eye radar architecture provides wide possibility of digital modulation and different waveform generation. Simultaneous correlation and integration of thousands signals per second from each point of surveillance area allows not only detecting of low level signals ((low profile targets), but help to recognize and classify signals (targets) by using diversity signals, polarization modulation and intelligent processing. Proposed all digital radar architecture with distributed directional antenna array can provide a 3D space vector to the jammer by verification direction of arrival for signals sources and as result jam/spoof protection not only for radar systems, but for communication systems and any navigation constellation system, for both encrypted or unencrypted signals, for not limited number or close positioned jammers.
Bürgmann, Roland; Rosen, Paul A.; Fielding, Eric J.
Synthetic aperture radar interferometry (InSAR) from Earth-orbiting spacecraft provides a new tool to map global topography and deformation of the Earth's surface. Radar images taken from slightly different viewing directions allow the construction of digital elevation models of meter-scale accuracy. These data sets aid in the analysis and interpretation of tectonic and volcanic landscapes. If the Earth's surface deformed between two radar image acquisitions, a map of the surface displacement with tens-of-meters resolution and subcentimeter accuracy can be constructed. This review gives a basic overview of InSAR for Earth scientists and presents a selection of geologic applications that demonstrate the unique capabilities of InSAR for mapping the topography and deformation of the Earth.
Amezquita, R.; Rincon, O. J.; Torres, Y. M.; Amezquita, S.
The optical characteristics of Diffractive Optical Elements are determined by the properties of the photosensitive film on which they are produced. When working with photoresist plates, the most important property is the change in the plate's topography for different exposures. In this case, the required characterization involves a topographic measurement that can be made using digital holography. This work presents a digital holography system in which a hologram's phase map is obtained from a single recorded image. The phase map is calculated by applying a phase-shifting algorithm to a set of images that are created using a digital phase-shifting/tilteliminating procedure. Also, the curvatures, introduced by the imaging elements used in the experimental setup, are digitally compensated for using a polynomial fitting-method. The object's topography is then obtained from this modified phase map. To demonstrate the proposed procedure, the topography of patches exposed on a Shipley 1818 photoresist plate by microlithography equipment-which is currently under construction-is shown.
Moller, Delwyn; Heavey, Brandon; Sadowy, Gregory
Compact, highly customizable digital receivers are being developed for the system described in 'Radar Interferometer for Topographic Mapping of Glaciers and Ice Sheets' (NPO-43962), NASA Tech Briefs, Vol. 31, No. 7 (August 2007), page 72. The receivers are required to operate in unison, sampling radar returns received by the antenna elements in a digital beam-forming (DBF) mode. The design of these receivers could also be adapted to commercial radar systems. At the time of reporting the information for this article, there were no commercially available digital receivers capable of satisfying all of the operational requirements and compact enough to be mounted directly on the antenna elements. A provided figure depicts the overall system of which the digital receivers are parts. Each digital receiver includes an analog-to-digital converter (ADC), a demultiplexer (DMUX), and a field-programmable gate array (FPGA). The ADC effects 10-bit band-pass sampling of input signals having frequencies up to 3.5 GHz. The input samples are demultiplexed at a user-selectable rate of 1:2 or 1:4, then buffered in part of the FPGA that functions as a first-in/first-out (FIFO) memory. Another part of the FPGA serves as a controller for the ADC, DMUX, and FIFO memory and as an interface between (1) the rest of the receiver and (2) a front-panel data port (FPDP) bus, which is an industry-standard parallel data bus that has a high data-rate capability and multichannel configuration suitable for DBF. Still other parts of the FPGA in each receiver perform signal-processing functions. The digital receivers can be configured to operate in a stand-alone mode, or in a multichannel mode as needed for DBF. The customizability of the receiver makes it applicable to a broad range of system architectures. The capability for operation of receivers in either a stand-alone or a DBF mode enables the use of the receivers in an unprecedentedly wide variety of radar systems.
A theoretical model which explains basic properties of radar imaging of underwater bottom topography in tidal channels is presented. The surface roughness modulation is described by weak hydrodynamic interaction theory in the relaxation time approximation. In contrast to previous theories on short wave modulation by long ocean waves, a different approximation is used to describe short wave modulation by tidal flow over underwater bottom topography. The modulation depth is proportional to the relaxation time of the Bragg waves. The large modulation of radar reflectivity observed in SEASAT-SAR imagery of sand banks in the Southern Bight of the North Sea are explained by assuming that the relaxation time of 34 cm Bragg waves is of the order of 30-40 seconds.
Podest, Erika; McDonald, Kyle; Kimball, John; Randerson, James
We characterize differences in radar-derived freeze/thaw state, examining transitions over complex terrain and landscape disturbance regimes. In areas of complex terrain, we explore freezekhaw dynamics related to elevation, slope aspect and varying landcover. In the burned regions, we explore the timing of seasonal freeze/thaw transition as related to the recovering landscape, relative to that of a nearby control site. We apply in situ biophysical measurements, including flux tower measurements to validate and interpret the remotely sensed parameters. A multi-scale analysis is performed relating high-resolution SAR backscatter and moderate resolution scatterometer measurements to assess trade-offs in spatial and temporal resolution in the remotely sensed fields.
Fault growth and propagation during incipient continental rifting: Insights from a combined aeromagnetic and Shuttle Radar Topography Mission digital elevation model investigation of the Okavango Rift Zone, northwest Botswana
Kinabo, B. D.; Hogan, J. P.; Atekwana, E. A.; Abdelsalam, M. G.; Modisi, M. P.
Digital Elevation Models (DEM) extracted from the Shuttle Radar Topography Mission (SRTM) data and high-resolution aeromagnetic data are used to characterize the growth and propagation of faults associated with the early stages of continental extension in the Okavango Rift Zone (ORZ), northwest Botswana. Significant differences in the height of fault scarps and the throws across the faults in the basement indicate extended fault histories accompanied by sediment accumulation within the rift graben. Faults in the center of the rift either lack topographic expressions or are interpreted to have become inactive, or have large throws and small scarp heights indicating waning activity. Faults on the outer margins of the rift exhibit either (1) large throws or significant scarp heights and are considered older and active or (2) throws and scarp heights that are in closer agreement and are considered young and active. Fault linkages between major fault systems through a process of "fault piracy" have combined to establish an immature border fault for the ORZ. Thus, in addition to growing in length (by along-axis linkage of segments), the rift is also growing in width (by transferring motion to younger faults along the outer margins while abandoning older faults in the middle). Finally, utilization of preexisting zones of weakness allowed the development of very long faults (>100 km) at a very early stage of continental rifting, explaining the apparent paradox between the fault length versus throw for this young rift. This study clearly demonstrates that the integration of the SRTM DEM and aeromagnetic data provides a 3-D view of the faults and fault systems, providing new insight into fault growth and propagation during the nascent stages of continental rifting.
[figure removed for brevity, see original site] Click on the image for the animationAbout the animation: This simulated view of the potential effects of storm surge flooding on Lake Pontchartrain and the New Orleans area was generated with data from the Shuttle Radar Topography Mission. Although it is protected by levees and sea walls against storm surges of 18 to 20 feet, much of the city is below sea level, and flooding due to storm surges caused by major hurricanes is a concern. The animation shows regions that, if unprotected, would be inundated with water. The animation depicts flooding in one-meter increments. About the image: The city of New Orleans, situated on the southern shore of Lake Pontchartrain, is shown in this radar image from the Shuttle Radar Topography Mission (SRTM). In this image bright areas show regions of high radar reflectivity, such as from urban areas, and elevations have been coded in color using height data also from the SRTM mission. Dark green colors indicate low elevations, rising through yellow and tan, to white at the highest elevations. New Orleans is near the center of this scene, between the lake and the Mississippi River. The line spanning the lake is the Lake Pontchartrain Causeway, the world's longest overwater highway bridge. Major portions of the city of New Orleans are actually below sea level, and although it is protected by levees and sea walls that are designed to protect against storm surges of 18 to 20 feet, flooding during storm surges associated with major hurricanes is a significant concern. Data used in this image were acquired by the Shuttle Radar Topography Mission aboard the Space Shuttle Endeavour, launched on Feb. 11, 2000. SRTM used the same radar instrument that comprised the Spaceborne Imaging Radar-C/X-Band Synthetic Aperture Radar (SIR-C/X-SAR) that flew twice on the Space Shuttle Endeavour in 1994. SRTM was designed to collect 3-D measurements of the Earth's surface
Lovelace, J.; Bellamy, H.; Snell, E. H.; Borgstahl, G.
A low-cost, real-time digital topography system is under development which will replace x-ray film and nuclear emulsion plates. The imaging system is based on an inexpensive surveillance camera that offers a 1000x1000 array of 8 im square pixels, anti-blooming circuitry, and very quick read out. Currently, the system directly converts x-rays to an image with no phosphor. The system is small and light and can be easily adapted to work with other crystallographic equipment. Preliminary images have been acquired of cubic insulin at the NSLS x26c beam line. NSLS x26c was configured for unfocused monochromatic radiation. Six reflections were collected with stills spaced from 0.002 to 0.001 degrees apart across the entire oscillation range that the reflections were in diffracting condition. All of the reflections were rotated to the vertical to reduce Lorentz and beam related effects. This particular CCD is designed for short exposure applications (much less than 1 sec) and so has a relatively high dark current leading to noisy raw images. The images are processed to remove background and other system noise with a multi-step approach including the use of wavelets, histogram, and mean window filtering. After processing, animations were constructed with the corresponding reflection profile to show the diffraction of the crystal volume vs. the oscillation angle as well as composite images showing the parts of the crystal with the strongest diffraction for each reflection. The final goal is to correlate features seen in reflection profiles captured with fine phi slicing to those seen in the topography images. With this development macromolecular topography finally comes into the digital age.
Lorenz, R. D.; Cassini RADAR Team
Cassini RADAR data are used to construct a global, albeit sparsely-sampled, topography map, and to generate a hypsometric profile to compare with other planetary bodies. Titan’s hypsogram is unimodal and strikingly narrow compared with the terrestrial planets. To investigate topographic extremes, a novel variant on the classic hypsogram is introduced, with a logarithmic abscissa to highlight mountainous terrain. In such a plot, the top of the terrestrial hypsogram is quite distinct from those of Mars and Venus due to the ‘glacial buzz-saw’ that clips terrestrial topography above the snowline. In contrast to the positive skew seen in other hypsograms, with a long tail of positive relief due to mountains, there is an indication (weak, given the limited data for Titan so far) that the Titan hypsogram appears slightly negatively skewed, suggesting a significant population of unfilled depressions. Limited data permit only a simplistic comparison of Titan topography with other icy satellites but we find that the standard deviation of terrain height (albeit at different scales) is similar to those of Ganymede and Europa. The topography of terrestrial planets is sampled with the same coverage that we have for Titan to gauge what as-yet-undiscovered topographic surprises may yet be hidden by Titan’s haze.
An architecture for a Direct RF-digitization Type Digital Mode RADAR was developed at GSFC in 2008. Two variations of a basic architecture were developed for use on RADAR imaging missions using aircraft and spacecraft. Both systems can operate with a pulse repetition rate up to 10 MHz with 8 received RF samples per pulse repetition interval, or at up to 19 kHz with 4K received RF samples per pulse repetition interval. The first design describes a computer architecture for a Continuous Mode RADAR transceiver with a real-time signal processing and display architecture. The architecture can operate at a high pulse repetition rate without interruption for an infinite amount of time. The second design describes a smaller and less costly burst mode RADAR that can transceive high pulse repetition rate RF signals without interruption for up to 37 seconds. The burst-mode RADAR was designed to operate on an off-line signal processing paradigm. The temporal distribution of RF samples acquired and reported to the RADAR processor remains uniform and free of distortion in both proposed architectures. The majority of the RADAR's electronics is implemented in digital CMOS (complementary metal oxide semiconductor), and analog circuits are restricted to signal amplification operations and analog to digital conversion. An implementation of the proposed systems will create a 1-GHz, Direct RF-digitization Type, L-Band Digital RADAR--the highest band achievable for Nyquist Rate, Direct RF-digitization Systems that do not implement an electronic IF downsample stage (after the receiver signal amplification stage), using commercially available off-the-shelf integrated circuits.
Geudtner, D.; Zink, M.; Gierull, C.; Shaffer, S.
The on-orbit alignment of the antenna beams of both the X-band and C-band radar systems during operations of the shuttle radar topography mission/X-band synthetic aperture radar (SRTM/X-SAR)was a key requirement for achieving best interferometric performance.
Hamelin, Jennifer L.; Jackson, Mark C.; Kirchwey, Christopher B.; Pileggi, Roberto A.
The Shuttle Radar Topography Mission (SRTM) flew aboard Space Shuttle Endeavor February 2000 and used interferometry to map 80% of the Earth's landmass. SRTM employed a 200-foot deployable mast structure to extend a second antenna away from the main antenna located in the Shuttle payload bay. Mapping requirements demanded precision pointing and orbital trajectories from the Shuttle on-orbit Flight Control System (PCS). Mast structural dynamics interaction with the FCS impacted stability and performance of the autopilot for attitude maneuvers and pointing during mapping operations. A damper system added to ensure that mast tip motion remained with in the limits of the outboard antenna tracking system while mapping also helped to mitigate structural dynamic interaction with the FCS autopilot. Late changes made to the payload damper system, which actually failed on-orbit, required a redesign and verification of the FCS autopilot filtering schemes necessary to ensure rotational control stability. In-flight measurements using three sensors were used to validate models and gauge the accuracy and robustness of the pre-mission notch filter design.
Makineci, H. B.; Karabörk, H.
Digital elevation model, showing the physical and topographical situation of the earth, is defined a tree-dimensional digital model obtained from the elevation of the surface by using of selected an appropriate interpolation method. DEMs are used in many areas such as management of natural resources, engineering and infrastructure projects, disaster and risk analysis, archaeology, security, aviation, forestry, energy, topographic mapping, landslide and flood analysis, Geographic Information Systems (GIS). Digital elevation models, which are the fundamental components of cartography, is calculated by many methods. Digital elevation models can be obtained terrestrial methods or data obtained by digitization of maps by processing the digital platform in general. Today, Digital elevation model data is generated by the processing of stereo optical satellite images, radar images (radargrammetry, interferometry) and lidar data using remote sensing and photogrammetric techniques with the help of improving technology. One of the fundamental components of remote sensing radar technology is very advanced nowadays. In response to this progress it began to be used more frequently in various fields. Determining the shape of topography and creating digital elevation model comes the beginning topics of these areas. It is aimed in this work , the differences of evaluation of quality between Sentinel-1A SAR image ,which is sent by European Space Agency ESA and Interferometry Wide Swath imaging mode and C band type , and DTED-2 (Digital Terrain Elevation Data) and application between them. The application includes RMS static method for detecting precision of data. Results show us to variance of points make a high decrease from mountain area to plane area.
Sadr, Ramin; Satorius, Edgar; Robinett, J. Loris, Jr.; Olson, Erlend
Report discusses conceptual digital frequency synthesizer part of programmable local oscillator in radar-astronomy system. Phase must remain continuous during adjustments of frequency, phase noise must be low, and spectral purity must be high. Discusses theory of operation in some mathematical detail and presents new analysis of spectral purity of output.
The primary objective of the STS-99 mission was to complete high resolution mapping of large sections of the Earth's surface using the Shuttle Radar Topography Mission (SRTM), a specially modified radar system. This videotape shows technicians in clean room suits working on the SRTM in the Multi-Payload Processing Facility (MPPF).
Topography adjacent to Signal Corps Radar (S.C.R.) 296 Station 5, showing conditions before construction, May 28, 1943, this drawing shows the Bonita Ridge access road retaining wall and general conditions at Bonita Ridge before the construction of Signal Corps Radar (S.C.R.) 296 Station 5 - Fort Barry, Signal Corps Radar 296, Station 5, Transmitter Building Foundation, Point Bonita, Marin Headlands, Sausalito, Marin County, CA
Wagner, H. L.; Shuchman, R. A.
A digital processing and analysis scheme for use with digitized synthetic aperture radar data was developed. Using data from a four channel system, the imagery is preprocessed using specially designed software and then analyzed using preexisting facilities originally intended for use with MSS type data. Geometric and radiometric correction may be performed if desired, as well as classification analysis, Fast Fourier transform, filtering and level slice and display functions. The system provides low cost output in real time, permitting interactive imagery analysis. System information flow diagrams as well as sample output products are shown.
Kim, J.; Wan, W.; Lee, S.; Choi, Y.
Topographic reconstruction is a high priority task for the solid planet and satellite exploration missions. Laser/radar altimetry and stereo analyses have been widely used for this purpose and achieve high quality 3D topographic data over various planetary surfaces such as Venus, Mercury, Moon and Mars. However, in contrast with inner planet and satellite, the base data sets to compose digital topography over outer planets and satellites are very limited. Titan, the largest satellite of Saturn has also too limited data inventory to achieve sufficient spatial resolution in topographic data, in spite of increasing interests about the detailed topography owing to the recent interesting discoveries on methane fluvial system, aeolian geomorphologies and possible tectonic activity. Therefore the endeavours to increase the coverage of digital topography employing radargrammetry (Kirk et al. 2009), radar altimetry (Elachi, et al. 2005) and SARtopo (Stiles et al. 2009) have been actively conducted. Although these efforts result in the construction of a global topographic map, the consequent spatial resolutions of global topography are still poor (Lorentz et al. 2013). In this study, we tried to improve the coverage and the quality of Titan digital terrain model employing approaches as follows; 1) A semi-automated stereo matching scheme manipulating low signal-to-noise SAR image pair incorporating adaptive filtering and base topography, 2) the geodetic control improvement of stereo SAR pair based on altimetric measurements, 3) introduction of radarclinometry to refine the topography from stereo analyses. Especially together with the technical improvements to exploit SAR stereo pair, the possibility to mine height information from Visual Infrared Mapping Spectrometer (VIMS) was actively explored by the means of hybrid stereogrammetry between VIMS and SAR image pairs and photoclinometry. The developed scheme was applied for a few testing areas especially over Xanadu which is
There is some confusion within the glaciological community as to the accuracy of the basal topography derived from radar measurements. A number of texts and papers state that basal topography cannot be determined to better than one quarter of the wavelength of the radar system. On the other hand King et al (Nature Geoscience, 2009) claimed that features of the bed topography beneath Rutford Ice Stream, Antarctica can be distinguished to +/- 3m using a 3 MHz radar system (which has a quarter wavelength of 14m in ice). These statements of accuracy are mutually exclusive. I will show in this presentation that the measurement of ice thickness is a radar range determination to a single strongly-reflective target. This measurement has much higher accuracy than the resolution of two targets of similar reflection strength, which is governed by the quarter-wave criterion. The rise time of the source signal and the sensitivity and digitisation interval of the recording system are the controlling criteria on radar range accuracy. A dataset from Pine Island Glacier, West Antarctica will be used to illustrate these points, as well as the repeatability or precision of radar range measurements, and the influence of gridding parameters and positioning accuracy on the final DEM product.
Sadr, R.; Satorius, E.; Robinett, L.; Olson, E.
The digital frequency synthesizer (DFS) is an integral part of the programmable local oscillator (PLO) which is being developed for the NASA's Deep Space Network (DSN) and radar astronomy. Here, the theory of operation and the design of the DFS are discussed, and the design parameters in application for the Goldstone Solar System Radar (GSSR) are specified. The spectral purity of the DFS is evaluated by analytically evaluating the output spectrum of the DFS. A novel architecture is proposed for the design of the DFS with a frequency resolution of 1/2(exp 48) of the clock frequency (0.35 mu Hz at 100 MHz), a phase resolution of 0.0056 degrees (16 bits), and a frequency spur attenuation of -96 dBc.
Lorenz, Ralph; Kirk, R.; Stofan, E.; Lunine, J.; Hayes, A.; Stiles, B.; Mitchell, K.; Le Gall, A.; Zebker, H.; Wye, L.; Encrenaz, P.; Aharonson, O.; Lucas, A.; Janssen, M.; Notarnicola, C.; Casarano, D.; Ventura, B.; Cassini RADAR Team
Ligeia Mare is the best-mapped of Titan’s three seas, and has attracted particular interest as the target of the proposed TiME (Titan Mare Explorer) mission. Here we summarize radar observations of this 400km wide feature and its environs from Cassini flybys T25, T28, T29 and T65. As noted in studies of Ontario Lacus (Hayes ref), radar reflectivity can be used with assumptions to assess liquid depth in shallow areas. Most of Ligeia is well below the noise floor of our observations (which varies across the scene - we use the most sensitive central beam where available to pose the tightest sigma-0 constraint) indicating depths likely > 10m, although we delineate some possibly shallow margins to aid in future modeling of tidal currents. In addition, the brightness temperature measured by passive radiometry (Janssen et al., 2009) places a joint constraint on the surface temperature and the emissivity, suggesting an upper limit of 10% on suspended solid material. Combination of SAR imaging from the flybys permits construction of a stereo Digital Elevation Model. This stereo topography is compared with SARtopo measurements and shows a number of 1km high mountains in the surrounding terrain: the peaks of these mountains would be above the horizon as seen from much of Ligeia. The model also places constraints on the watershed of Ligeia and thus on the hydrological balance of precipitation and evaporation. We will also report on further observations of Ligeia planned in the T86 flyby, shortly before the DPS meeting.
Weighting is employed in synthetic aperture radar (SAR) processing to reduce the sidelobe response at the expense of peak center response height and mainlobe resolution. The weighting effectiveness in digital processing depends not only on the choice of weighting function, but on the fineness of sampling and quantization, on the time bandwidth product, on the quadratic phase error, and on the azimuth antenna pattern. The results of simulations conducted to uncover the effect of these parameters on azimuth weighting effectiveness are presented. In particular, it is shown that multilook capabilities of future SAR systems may obviate the need for consideration of the antenna pattern, and that azimuth time-bandwidth products of over 200 are probably required before the digital results begin to approach the ideal results.
Griffin, C. R.; Estes, J. M.
A modified APQ-102 sidelooking radar collected synthetic aperture radar (SAR) data which was digitized and recorded on wideband magnetic tape. These tapes were then ground processed into computer compatible tapes (CCT's). The CCT's may then be processed into high resolution radar images by software on the CYBER computer.
Griffin, C. R.; Estes, J. M.
A modified APQ-102 sidelooking radar collected synthetic aperture radar (SAR) data which was digitized and recorded on wideband magnetic tape. These tapes were then ground processed into computer compatible tapes (CCT's). The CCT's may then be processed into high resolution radar images by software on the CYBER computer.
Ramsey, Elijah W., III; Nelson, G.A.; Laine, S.C.; Kirkman, R.G.; Topham, W.
A topographic surface of a low lying coastal marsh was created by using three flood extent vectors digitized from ERS-1 SAR images and two elevation contours from U.S. Geological Survey topographic quadrangles. Point measurement of water depth at the times of the SAR collections allowed conversion of the radar measured flood extent vectors to topographic contours. Generation of the topographic surface was accomplished with a surface gridding algorithm. SAR and on-site measures. Errors in the generated topography were mainly associated with the lack of input contours covering narrow to broad plateaus and topographic highs and lows. The misplacement of SAR derived flood extent vectors also caused errors in sparsely vegetated high marsh at convoluted marsh-forest boundaries, and at topographic depressions. Overall, the standard deviation of differences between measured and predicted elevations at 747 points was 19 cm. Excluding the above mentioned abrupt boundaries and topographic highs and lows outside the range of available contours, standard deviation differences averaged about 14 cm, but most often averaged about 8 cm. This suggested a 5 to 9 factor improvement over the 150 cm topographic resolution currently available for this area.
Ramsey, Elijah W., III; Nelson, G.A.; Laine, S.C.; Kirkman, R.G.; Topham, W.
A topographic surface of a low lying coastal marsh was created by using three flood extent vectors digitized from ERS-1 SAR images and two elevation contours from U.S. Geological Survey topographic quadrangles. Point measurement of water depth at the times of the SAR collections allowed conversion of the radar measured flood extent vectors to topographic contours. Generation of the topographic surface was accomplished with a surface gridding algorithm, SAR and on-site measures. Errors in the generated topography were mainly associated with the lack of input contours covering narrow to broad plateaus and topographic highs and lows. The misplacement of SAR derived flood extent vectors also caused errors in sparsely vegetated high marsh, at convoluted marsh-forest boundaries, and at topographic depressions. Overall, the standard deviation of differences between measured and predicted elevations at 747 points was 19 cm. Excluding the above mentioned abrupt boundaries and topographic highs and lows outside the range of available contours, standard deviation differences averaged about 14 cm, but most often averaged about 8 cm. This suggested a 5 to 9 factor improvement over the 150 cm topographic resolution currently available for this area.
These four images of the Long Valley region of east-central California illustrate the steps required to produced three dimensional data and topographics maps from radar interferometry. All data displayed in these images were acquired by the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) aboard the space shuttle Endeavour during its two flights in April and October, 1994. The image in the upper left shows L-band (horizontally transmitted and received) SIR-C radar image data for an area 34 by 59 kilometers (21 by 37 miles). North is toward the upper right; the radar illumination is from the top of the image. The bright areas are hilly regions that contain exposed bedrock and pine forest. The darker gray areas are the relatively smooth, sparsely vegetated valley floors. The dark irregular patch near the lower left is Lake Crowley. The curving ridge that runs across the center of the image from top to bottom is the northeast rim of the Long Valley Caldera, a remnant crater from a massive volcanic eruption that occurred about 750,000 years ago. The image in the upper right is an interferogram of the same area, made by combining SIR-C L-band data from the April and October flights. The colors in this image represent the difference in the phase of the radar echoes obtained on the two flights. Variations in the phase difference are caused by elevation differences. Formation of continuous bands of phase differences, known as interferometric 'fringes', is only possible if the two observations were acquired from nearly the same position in space. For these April and October data takes, the shuttle tracks were less than 100 meters (328 feet) apart. The image in the lower left shows a topographic map derived from the interferometric data. The colors represent increments of elevation, as do the thin black contour lines, which are spaced at 50-meter (164-foot) elevation intervals. Heavy contour lines show 250-meter intervals (820-foot). Total relief in
Treuhaft, Robert N.
This paper first gives a heuristic description of the sensitivity of Interferometric Synthetic Aperture Radar to vertical vegetation distributions and underlying surface topography. A parameter estimation scenario is then described in which the Interferometric Synthetic Aperture Radar cross-correlation amplitude and phase are the observations from which vegetation and surface topographic parameters are estimated. It is shown that, even in the homogeneous-layer model of the vegetation, the number of parameters needed to describe the vegetation and underlying topography exceeds the number of Interferometric Synthetic Aperture Radar observations for single-baseline, single-frequency, single-incidence-angle, single-polarization Interferometric Synthetic Aperture Radar. Using ancillary ground-truth data to compensate for the underdetermination of the parameters, forest depths are estimated from the INSAR data. A recently-analyzed multibaseline data set is also discussed and the potential for stand-alone Interferometric Synthetic Aperture Radar parameter estimation is assessed. The potential of combining the information content of Interferometric Synthetic Aperture Radar with that of infrared/optical remote sensing data is briefly discussed.
Rincon, Rafael F.; Perrine, Martin; McLinden, Matthew; Valett, Susan
The Digital Beamforming Synthetic Aperture Radar (DBSAR) is a state-of-the-art radar system developed at NASA/Goddard Space Flight Center for the development and implementation of digital beamforming radar techniques. DBSAR was recently upgraded to polarimetric operation in order to enhance its capability as a science instrument. Two polarimetric approaches were carried out which will be demonstrated in upcoming flight campaigns.
Wu, S. S. C.
The application of digital processing techniques to spacecraft television pictures and radar images is discussed. The use of digital rectification to produce contour maps from spacecraft pictures is described; images with azimuth and elevation angles are converted into point-perspective frame pictures. The digital correction of the slant angle of radar images to ground scale is examined. The development of orthophoto and stereoscopic shaded relief maps from digital terrain and digital image data is analyzed. Digital image transformations and rectifications are utilized on Viking Orbiter and Lander pictures of Mars.
Wu, S. S. C.
The application of digital processing techniques to spacecraft television pictures and radar images is discussed. The use of digital rectification to produce contour maps from spacecraft pictures is described; images with azimuth and elevation angles are converted into point-perspective frame pictures. The digital correction of the slant angle of radar images to ground scale is examined. The development of orthophoto and stereoscopic shaded relief maps from digital terrain and digital image data is analyzed. Digital image transformations and rectifications are utilized on Viking Orbiter and Lander pictures of Mars.
Johnson, Jeffrey R.; Gaddis, Lisa
The FLOW computer model of McEwen and Malin (1989) modified for application to the study of Venus fluidized ejecta blankets (FEBs) demonstrates that relatively low viscosities, yield strengths, and initial velocities are required to duplicate the observed flow paths of the outflow materials. The model calculates the velocities and simulated flow paths of gravity flows over Magellan topography. The model is formulated to determine flow movements from initial conditions, gravitational acceleration, and resistance to motion as described by Coulomb, viscous, and turbulent resistance forces. Successful duplication of observed FEB flow paths has been achieved for the FEB craters Addams, Isabella, and Cochran. When used as a simple energy-line model, the model requires low coefficients of friction to extend FEBs to near their observed termini in the synthetic aperture radar (SAR) imagery, although the resulting straight flow lines do not follow the observed flow paths well. For Bingham flow, the model requires low values of viscosity and yield strength which are more similar to pyroclastic or debris flows than basaltic lavas. Flows of 100-m depth require 1 to 2 orders of magnitude higher values of both viscosity and yield strength than 10-m-deep flows. The complicated nature of the flow lines for the low velocity model suggests that FEBs were probably emplaced under variably laminar and turbulent flow conditions, where underlying topography influenced both the direction and energy of flow materials.
Domik, G.; Leberl, F.; Raggam, J.
The refinement of radar image analysis methods has led to a need for a systems approach to radar image processing software. Developments stimulated through satellite radar are combined with standard image processing techniques to create a user environment to manipulate and analyze airborne and satellite radar images. One aim is to create radar products for the user from the original data to enhance the ease of understanding the contents. The results are called secondary image products and derive from the original digital images. Another aim is to support interactive SAR image analysis. Software methods permit use of a digital height model to create ortho images, synthetic images, stereo-ortho images, radar maps or color combinations of different component products. Efforts are ongoing to integrate individual tools into a combined hardware/software environment for interactive radar image analysis.
Digital image transformation and rectification can be described in three categories: (1) digital rectification of spacecraft pictures on workable stereoplotters; (2) digital correction of radar image geometry; and (3) digital reconstruction of shaded relief maps and perspective views including stereograms. Digital rectification can make high-oblique pictures workable on stereoplotters that would otherwise not accommodate such extreme tilt angles. It also enables panoramic line-scan geometry to be used to compile contour maps with photogrammetric plotters. Rectifications were digitally processed on both Viking Orbiter and Lander pictures of Mars as well as radar images taken by various radar systems. By merging digital terrain data with image data, perspective and three-dimensional views of Olympus Mons and Tithonium Chasma, also of Mars, are reconstructed through digital image processing. ?? 1985.
Lerche, H. D.; Tumbreagel, F.
The discussed navigation update system was designed for an unmanned platform with fire and forget capability. It meets the requirement due to fully autonomous operation. The system concept will be characterized by complementary use of the radar seeker for target identification as well as for navigation function. The system works in the navigation mode during preprogrammable phases where the primary target identification function is not active or in parallel processing. The dual function radar seeker system navigates the drone during the midcourse and terminal phases of the mission. Its high resolution due to range measurement and doppler beam sharpening in context with its radar reflectivity sensing capability are the basis for topography referenced navigation computation. The detected height jumps (coming from terrain elevation and cultural objects) and radar reflectivity features will be matched together with topography referenced features. The database comprises elevation data and selected radar reflectivity features that are robust against seasonal influences. The operational benefits of the discussed system are as follows: (1) the improved navigation performance with high probability of position fixing, even over flat terrain; (2) the operation within higher altitudes; and (3) bad weather capability. The developed software modules were verified with captive flight test data running in a hardware-in-the-loop simulation.
Wilson, Kelce; Patrick, Dale; Blair, James
The statistic results for a digital terrain model are presented that closely match measurements for 77% of the 189 possible combinations of 7 radar bands, 3 polarizations, and 9 terrain types. The model produces realistic backscatter coefficient values for the scenarios over all incidence angles from normal to grazing. The generator was created using measured data sets reported in the Handbook of Radar Scattering Statistics for Terrain covering L, C, S, X, Ka, Ku, and W frequency bands; HH, HV, and VV polarizations; and soil and rock, shrub, tree, short vegetation, grass, dry snow, wet snow, road surface, and urban area terrain types. The first two statistical moments match published values precisely, and a Chi-Square histogram test failed to reject the generator at a 95% confidence level for the 146 terrain models implemented. A Sea State model provides the grazing angle extension for predictions beyond the available measurements. This work will contain a comprehensive set of plots of mean and standard deviation versus incidence angle.
Fu, Lee-Lueng; Rodriguez, Ernesto
We propose to apply the technique of synthetic aperture radar interferometry to the measurement of ocean surface topography at spatial resolution approaching 1 km. The measurement will have wide ranging applications in oceanography, hydrology. and marine geophysics. The oceanographic and related societal applications are briefly discussed in the paper. To meet the requirements for oceanographic applications, the instrument must be flown in an orbit with proper sampling of ocean tides.
An earth-based radar topography (ERT) map has been constructed of the Mare Crisium area. Systematic and random sources of error are discussed. A comparison between the ERT map and Lunar Topographic Orthophotomaps shows a random mean discrepancy of less than 100 m between the two maps, except for small-scale (20 km or less in diameter) features, where systematic smoothing reduces the ERT elevation contrast
Hou, Qingkai; Liu, Yang; Chen, Zengping; Su, Shaoying
Digital orthogonal receiver is one of the key techniques in digital receiver of soft radar, and compressed sensing is attracting more and more attention in radar signal processing. In this paper, we propose a CS digital orthogonal receiver for wideband radar which utilizes compressed sampling in the acquisition of radar raw data. In order to reconstruct complex signal from sub-sampled raw data, a novel sparse dictionary is proposed to represent the real-valued radar raw signal sparsely. Using our dictionary and CS algorithm, we can reconstruct the complex-valued radar signal from sub-sampled echoes. Compared with conventional digital orthogonal radar receiver, the architecture of receiver in this paper is more simplified and the sampling frequency of ADC is reduced sharply. At the same time, the range profile can be obtained during the reconstruction, so the matched filtering can be eliminated in the receiver. Some experiments on ISAR imaging based on simulated data prove that the phase information of radar echoes is well reserved in our orthogonal receiver and the whole design is effective for wideband radar.
The Space-Shuttle Radar Topography Mission provided geologists with a detailed digital elevation model of most of Earth's land surface. This new database is used here for structural analysis of grooved surfaces interpreted to be the exhumed footwalls of three active or recently active extensional detachment faults. Exhumed fault footwalls, each with an areal extent of one hundred to several hundred square kilometers, make up much of Dayman dome in eastern Papua New Guinea, the western Gurla Mandhata massif in the central Himalaya, and the northern Tokorondo Mountains in central Sulawesi, Indonesia. Footwall curvature in profile varies from planar to slightly convex upward at Gurla Mandhata to strongly convex upward at northwestern Dayman dome. Fault curvature decreases away from the trace of the bounding detachment fault in western Dayman dome and in the Tokorondo massif, suggesting footwall flattening (reduction in curvature) following exhumation. Grooves of highly variable wavelength and amplitude reveal extension direction, although structural processes of groove genesis may be diverse.
Rincon, Rafael F.; Fatoyinbo, Temilola; Carter, Lynn; Ranson, K. Jon; Vega, Manuel; Osmanoglu, Batuhan; Lee, SeungKuk; Sun, Guoqing
The Digital Beamforming Synthetic Aperture radar (DBSAR) is a state-of-the-art airborne radar developed at NASA/Goddard for the implementation, and testing of digital beamforming techniques applicable to Earth and planetary sciences. The DBSAR measurements have been employed to study: The estimation of vegetation biomass and structure - critical parameters in the study of the carbon cycle; The measurement of geological features - to explore its applicability to planetary science by measuring planetary analogue targets. The instrument flew two test campaigns over the East coast of the United States in 2011, and 2012. During the campaigns the instrument operated in full polarimetric mode collecting data from vegetation and topography features.
Moller, D. K.; Aaron, K.; Gim, Y.; Heavey, B.; Hodges, R.; Nicolson, A.; Rengarajan, S.; Rignot, E.; Rogez, F.; Sadowy, G.; Simard, M.; Zawadzki, M.
The estimation of the mass balance of ice sheets and glaciers on Earth is a problem of considerable scientific and societal importance. The Greenland and Antarctic ice sheets together hold enough ice to raise global sea level by 80 m. The annual exchange of mass on the ice sheets is equivalent to 8mm/yr sea level, so that any fluctuation in that level of exchange is significant on the global scale. A key measurement to understanding, monitoring and forecasting these changes is ice-surface topography, both for ice-sheet and glacial regions. As such NASA identified "ice topographic mapping instruments capable of providing precise elevation and detailed imagery data for measurements on glacial scales for detailed monitoring of ice sheet, and glacier changes" as a science priority for the most recent ESTO- Instrument Incubator Program (IIP) opportunities. Funded under this opportunity is the technological development for a Ka-Band (35GHz) single-pass digitally beamformed interferometric synthetic aperture radar (InSAR). Unique to this concept is the ability to map a significant swath impervious of cloud cover with measurement accuracies comparable to lidar altimeters but with variable resolution as appropriate to the differing scales-of-interest over ice-sheets and glaciers. By diverging from the more traditional profiling measurements employed to date (ie radar altimeters and lidars) we are able to offer the potential to significantly advance the spaciotemporal observational capabilities of both ice sheets and glaciers. Dubbed the Glacier and Land Ice Surface Topography Interferometer (GLISTIN), the instrument and mission presents several significant challenges. In particular, under the IIP program we are designing, building and demonstrating a large Ka-band antenna array with integrated digital receivers and utilizing digital beamforming to preserve both antenna gain and swath. These technology items will ultimately be integrated into a complete interferometric
Stepinski, T. F.; Collier, M. L.
We have developed a novel method for delineating valley networks on Mars. The valleys are inferred from digital topography by an autonomous computer algorithm as drainage networks, instead of being manually mapped from images. Individual drainage basins are precisely defined and reconstructed to restore flow continuity disrupted by craters. Drainage networks are extracted from their underlying basins using the contributing area threshold method. We demonstrate that such drainage networks coincide with mapped valley networks verifying that valley networks are indeed drainage systems. Our procedure is capable of delineating and analyzing valley networks with unparalleled speed and consistency. We have applied this method to 28 Noachian locations on Mars exhibiting prominent valley networks. All extracted networks have a planar morphology similar to that of terrestrial river networks. They are characterized by a drainage density of approx.0.1/km, low in comparison to the drainage density of terrestrial river networks. Slopes of "streams" in Martian valley networks decrease downstream at a slower rate than slopes of streams in terrestrial river networks. This analysis, based on a sizable data set of valley networks, reveals that although valley networks have some features pointing to their origin by precipitation-fed runoff erosion, their quantitative characteristics suggest that precipitation intensity and/or longevity of past pluvial climate were inadequate to develop mature drainage basins on Mars.
Young, R.A.; Sun, Jingsheng
The topography of a clay aquitard is defined by 3D Ground Penetrating Radar (GPR) data at Hill Air Force Base, Utah. Conventional processing augmented by multichannel domain filtering shows a strong reflection from a depth of 20-30 ft despite attenuation by an artificial clay cap approximately 2 ft thick. This reflection correlates very closely with the top of the aquitard as seen in lithology logs at 3 wells crossed by common offset radar profiles from the 3D dataset. Lateral and vertical resolution along the boundary are approximately 2 ft and 1 ft, respectively. The boundary shows abrupt topographic variation of 5 ft over horizontal distances of 20 ft or less and is probably due to vigorous erosion by streams during lowstands of ancient Lake Bonneville. This irregular topography may provide depressions for accumulation of hydrocarbons and chlorinated organic pollutants. A ridge running the length of the survey area may channel movement of ground water and of hydrocarbons trapped at the surface of the water table. Depth slices through a 3D volume, and picked points along the aquitard displayed in depth and relative elevation perspectives provide much more useful visualization than several 2D lines by themselves. The three-dimensional CPR image provides far more detailed definition of geologic boundaries than does projection of soil boring logs into two-dimensional profiles.
Mercuri, Pablo Alberto
Digital Elevation Models (DEMs) are increasingly used even in low relief landscapes for multiple mapping applications and modeling approaches such as surface hydrology, flood risk mapping, agricultural suitability, and generation of topographic attributes. The National Aeronautics and Space Administration (NASA) has produced a nearly global database of highly accurate elevation data, the Shuttle Radar Topography Mission (SRTM) DEM. The main goals of this thesis were to investigate quality issues of SRTM, provide measures of vertical accuracy with emphasis on low relief areas, and to analyze the performance for the generation of physical boundaries and streams for watershed modeling and characterization. The absolute and relative accuracy of the two SRTM resolutions, at 1 and 3 arc-seconds, were investigated to generate information that can be used as a reference in areas with similar characteristics in other regions of the world. The absolute accuracy was obtained from accurate point estimates using the best available federal geodetic network in Indiana. The SRTM root mean square error for this area of the Midwest US surpassed data specifications. It was on the order of 2 meters for the 1 arc-second resolution in flat areas of the Midwest US. Estimates of error were smaller for the global coverage 3 arc-second data with very similar results obtained in the flat plains in Argentina. In addition to calculating the vertical accuracy, the impacts of physiography and terrain attributes, like slope, on the error magnitude were studied. The assessment also included analysis of the effects of land cover on vertical accuracy. Measures of local variability were described to identify the adjacency effects produced by surface features in the SRTM DEM, like forests and manmade features near the geodetic point. Spatial relationships among the bare-earth National Elevation Data and SRTM were also analyzed to assess the relative accuracy that was 2.33 meters in terms of the total
Jassar, H. K. Al; Rao, K. S.
Using different combinations of 29 Advanced Synthetic Aperture Radar (ASAR) images, 43 Digital Elevations Models (DEM) were generated adopting SAR Interferometry (InSAR) technique. Due to sand movement in desert terrain, there is a poor phase correlation between different SAR images. Therefore, suitable methodology for generating DEMs of Kuwait desert terrain using InSAR technique were worked out. Time series analysis was adopted to derive the best DEM out of 43 DEMs. The problems related to phase de-correlation over desert terrain are discussed. Various errors associated with the DEM generation are discussed which include atmospheric effects, penetration into soil medium, sand movement. The DEM of Shuttle Radar Topography Mission (SRTM) is used as a reference. The noise levels of DEM of SRTM are presented.
Huang, J.; Turcotte, D. L.
The concept of fractal mapping is introduced and applied to digitized topography of Arizona. It is shown that the fractal statistics satisfy the topography of the state to a good approximation. The fractal dimensions and roughness amplitudes from subregions are used to construct maps of these quantities. It is found that the fractal dimension of actual two-dimensional topography is not affected by the adding unity to the fractal dimension of one-dimensional topographic tracks. In addition, consideration is given to the production of fractal maps from synthetically derived topography.
74. Transmitter building no. 102, view of radar digital test and maintenance cabinet area control panel and date storage system showing ampex tape storage devices. - Clear Air Force Station, Ballistic Missile Early Warning System Site II, One mile west of mile marker 293.5 on Parks Highway, 5 miles southwest of Anderson, Anderson, Denali Borough, AK
Treuhaft, Robert N.
Drawing from recently submitted work, this paper first gives a heuristic description of the sensitivity of interferometric synthetic aperture radar (INSAR) to vertical vegetation distribution and under laying surface topography. A parameter estimation scenario is then described in which the INSAR cross correlation amplitude and phase are the observations from which vegetation and surface topographic parameters are estimated. It is shown that, even in the homogeneous layer model of the vegetation, the number of parameters needed to describe the vegetation and underlying topography exceeds the number of INSAR observations for single baseline, single frequency, single incidence-angle, single polarization INSAR. Using ancillary ground truth data to compensate for the under determination of the parameters, forest depths are estimated from the INSAR data. A recently analyzed multi-baseline data set is also discussed and the potential for stand alone INSAR parameter estimation is assessed. The potential of combining the information content of INSAR with that of infrared/optical remote sensing data is briefly discussed.
Morlighem, M.; Rignot, E. J.; Mouginot, J.; Seroussi, H. L.
Bed topography, together with ice thickness, is an essential characteristic of glaciers and ice sheets for many glaciological applications. Despite significant technical advances, it remains challenging to measure ice thickness remotely, especially in deep troughs occupied by outlet glaciers. The method of mass conservation, that combines radar-derived ice thickness data with high-resolution InSAR-derived ice velocity vectors, provides an effective method for generating a high-resolution bed from sparse radar sounding profiles, and has been successfully applied along the coast of the Greenland Ice Sheet. Applying the same technique to the coast of the Antarctic Ice Sheet presents a number of challenges. The coverage of ice thickness data collected in Antarctica, for example, is much less comprehensive compared to Greenland, especially in the wake of NASA's Operation IceBridge (OIB) Mission in 2010-2015. Here, we combine radar sounder data collected by various centers (OIB/Center for Remote Sensing of Ice Sheets, the British Antarctic Survey and University of Texas) acquired between 1998 and 2011, with high-resolution ice motion data from interferometric SAR (ALOS PALSAR, RADARSAT-2 and Envisat ASAR) to reconstruct bed topography beneath major Antarctic outlet glaciers at an unprecedented level of detail. The results reveal some important features not known previously at that level of detail and shed light on the vulnerability of these glaciers in a warming climate. We find for example that Recovery glacier is deeper than in previous mappings and has long grooves parallel to the flow direction. Denman Glacier, East Antarctica, flow along a deep, narrow trough more than 2,000 m below sea level that extends more than 100 km inland. We find ridges and bumps in the vicinity of the grounding line of Thwaites Glacier, in the Amundsen Sea sector, that are consistent with the pattern of grounding line retreat. We have also a new mapping of the trough upstream of David
Adams, John W.; Nelson, Jeffrey E.; Banh, N. D.; Moncada, John J.; Bayma, Robert W.
Novel weighted-least-squares approaches to the design of digital filters for SAR applications are presented. The filters belong to three different categories according to their combinations of minimax passband, least-squares stopband, minimax stopband, and maximally-flat passband. For real-time applications, it is important to design the sets of digital filter coefficient tables in an offline environment; the appropriate precomputed filter is then selected for each SAR signal-processing function, as a function of both mode and mapping geometry during real-time processing.
Moller, Delwyn K.; Heavey, Brandon; Hodges, Richard; Rengarajan, Sembiam; Rignot, Eric; Rogez, Francois; Sadowy, Gregory; Simard, Marc; Zawadzki, Mark
The estimation of the mass balance of ice sheets and glaciers on Earth is a problem of considerable scientific and societal importance. A key measurement to understanding, monitoring and forecasting these changes is ice-surface topography, both for ice-sheet and glacial regions. As such NASA identified 'ice topographic mapping instruments capable of providing precise elevation and detailed imagery data for measurements on glacial scales for detailed monitoring of ice sheet, and glacier changes' as a science priority for the most recent Instrument Incubator Program (IIP) opportunities. Funded under this opportunity is the technological development for a Ka-Band (35GHz) single-pass digitally beamformed interferometric synthetic aperture radar (InSAR). Unique to this concept is the ability to map a significant swath impervious of cloud cover with measurement accuracies comparable to laser altimeters but with variable resolution as appropriate to the differing scales-of-interest over ice-sheets and glaciers.
A tutorial on synthetic aperture radar (SAR) is presented with emphasis on digital data collection and processing. Background information on waveform frequency and phase notation, mixing, Q conversion, sampling and cross correlation operations is included for clarity. The fate of a SAR signal from transmission to processed image is traced in detail, using the model of a single bright point target against a dark background. Some of the principal problems connected with SAR processing are also discussed.
Arnold, Stephen; Hsu, Charles C.; Zaghloul, Mona E.; Szu, Harold H.; Karangelen, Nicholas E.; Buss, James R.
A high performance, fully digital Foliage Penetrating Synthetic Aperture Radar (FOPEN SAR) system is described. The FOPEN SAR algorithm is illustrated using Matlab. Digital implementation is derived and simulated using VHDL. The complex mathematical functions required by the algorithm have been demonstrated. Simulations have achieved an SNR equals 290 dB when compared to the baseline results from Matlab. The accuracy of the simulation was limited by the resolution of certain trigonometric and exponential functions implemented using VHDL, and thus can be improved upon. This would allow greater flexibility between speed/area considerations without degradation of the target resolution (100dB-signal accuracy).
Shirzaei, M.; Bürgmann, R.
Atmospheric delay is one of the major sources of error in repeat pass interferometry. We propose a new approach for correcting the topography-correlated components of this artifact. To this aim we use multiresolution wavelet analysis to identify the components of the unwrapped interferogram that correlate with topography. By using a forward wavelet transform we break down the digital elevation model and the unwrapped interferogram into their building blocks based on their frequency properties. We apply a cross-correlation analysis to identify correlated coefficients that represent the effect of the atmospheric delay. Thus, the correction to the unwrapped interferogram is obtained by down-weighting the correlated coefficients during inverse wavelet transform. We test this approach on real and synthetic data sets that are generated over the San Francisco Bay Area. We find that even in the presence of tectonic signals, this method is able to reduce the correlated component of the atmospheric delay by up to 75% and improves the signal in areas of high relief. The remaining part is most likely due to 3D heterogeneities of the atmosphere and can be reduced by integrating temporal information or using complementary observations or models of atmospheric delay.
Horwitz, Dennis N.
The STS-99 Shuttle Radar Topography Mission (SRTM) employed radar interferometry to gather high resolution imagery used to generate the most detailed 3D map of the earth's surface ever produced. Such a map has a broad range of both military and commercial uses. This 11-day mission of the Space Shuttle Endeavour took place from February 11 to 22, 2000, and covered 80% of the earth's surface. The SRTM project gathered 12.3 Terabytes of imaging data, which is equivalent to more than 20,418 compact disks, and approximately equal to the entire contents of the Library of Congress.
Imhoff, Marc Lee; Gesch, Dean B.
Synthetic aperture radar data from the Shuttle Imaging Radar-B Mission were combined with the tide surface information to create a digital terrain model for a 70-km by 40-km section of the Mouths of the Ganges forests in southern Bangladesh. The dominance of the interaction phenomenon (canopy to surface or surface to canopy reflection) in flooded forests was exploited to create sub-canopy flood boundary maps for two different tide times. The boundary maps were digitally combined in x, y, z space with tide elevation models created from tide gauge data gridding the survey site and used as input to interpolation routines to create a terrain model. The end product represents a significant step in our ability to characterize the topography and hydrology of wetland ecosystems. The model derived here can be used for simulating tidal flow and nutrient transport from the forest to the marine habitat.
Lovelace, J.; Soares, A. S.; Bellamy, H.; Sweet, R. M.; Snell, E. H.; Borgstahl, G.
An inexpensive digital CCD camera was used to record X-ray topographs directly from large imperfect crystals of cubic insulin. The topographs recorded were not as detailed as those which can be measured with film or emulsion plates but do show great promise. Six reflections were recorded using a set of finely spaced stills encompassing the rocking curve of each reflection. A complete topographic reflection profile could be digitally imaged in minutes. Interesting and complex internal structure was observed by this technique.The CCD chip used in the camera has anti-blooming circuitry and produced good data quality even when pixels became overloaded.
Lovelace J. J.; Soares A.; Bellamy, H. D.; Sweet, R. M.; Snell, E. H.; Borgstahl, G. E. O.
An inexpensive digital CCD camera was used to record X-ray topographs directly from large imperfect crystals of cubic insulin. The topographs recorded were not as detailed as those which can be measured with film or emulsion plates, but do show great promise. Six reflections were recorded using a set of finely spaced stills encompassing the rocking curve of each reflection. A complete topographic reflection profile could be digitally imaged in minutes. Interesting and complex internal structure was observed by this technique. The CCD chip used in the camera has anti-blooming circuitry and produced good data quality, even when pixels became overloaded.
Fischman, Mark; Berkun, Andrew; Chu, Anhua; Freedman, Adam; Jourdan, Michael; McWatters, Dalia; Paller, Mimi
A digital receiver in a 1.26-GHz spaceborne radar scatterometer now undergoing development includes a module for detecting radio-frequency interference (RFI) that could contaminate scientific data intended to be acquired by the scatterometer. The role of the RFI-detection module is to identify time intervals during which the received signal is likely to be contaminated by RFI and thereby to enable exclusion, from further scientific data processing, of signal data acquired during those intervals. The underlying concepts of detection of RFI and rejection of RFI-contaminated signal data are also potentially applicable in advanced terrestrial radio receivers, including software-defined radio receivers in general, receivers in cellular telephones and other wireless consumer electronic devices, and receivers in automotive collision-avoidance radar systems.
Walsh, Edward J.
Work continues on estimating tilt modulation distortions in the wave topography measured by a scanning radar altimeter. To quantify this effect, a two-dimensional simulation has been performed in the cross-track plane only which assumed that sinusoidal waves of constant wavelength propagate in the cross-track direction (and are infinitely long-crested in the along-track direction). The initial results reported earlier for a Gaussian surface scattering model indicated that when the highest reasonable value of mss is used in the simulation (the Plant limit of 0.08), the nadir values of the ratio of the apparent to actual wave height are the same as for the omnidirectional scattering case. But as the off-nadir angle increased, the apparent wave height increased and became larger than the actual wave height by about 10 degrees off-nadir. And the shorter the wavelength, the larger the apparent wave height increase. This represented a systematic over-estimate of the wave amplitude for waves propagating in the cross-track direction. For lower values of mss the situation worsened. The simulation has been improved by incorporating actual variations of backscattered power with incidence angle measured by the SRA instead of the Gaussian model. The resulting distortion was about half that originally reported. The 3-dimensional simulation to model waves propagating at various azimuthal angles relative to the cross-track plane is still in progress. The results of this model will be verified by comparing them with Scanning Radar Altimeter (SRA) data and optimum correction procedures will be developed. An assessment of the improved 2-dimensional model indicates that the distortion will generally be small in the data taken during the Southern Ocean Waves Experiment (SOWEX) presently be analyzed. But tilt modulation effects are of great concern for SRA data collected during the 1998 hurricane season since the minimum aircraft altitude was 1.5 km and it frequently flew higher. For a
Liu, Wen-Cheng; Huang, Wei-Che
Landslide monitoring is a crucial tool for the prevention of hazards. It is often the only solution for the survey and the early-warning of large landslides cannot be stabilized. The objective of present study is to use a low-cost image system to monitor the active landslides. We adopted the direct linear transformation (DLT) method in close range digital photogrammetry to measure terrain of landslide at the Huoyen Shan, Miaoli of central Taiwan and to compare measured results with e-GPS. The results revealed that the relative error in surface area was approximately 1.7% as comparing the photogrammetry with DLT method and e-GPS measurement. It showed that the close range digital photogrammetry with DLT method had the availability and capability to measure the landslides. The same methodology was then applied to measure the terrain before landslide and after landslide in the study area. The digital terrain model (DTM) was established and then was used to calculate the volume of the terrain before landslide and after landslide. The volume difference before and after landslides was 994.16 m3.
Yamamoto, Masayuki K.; Fujita, Toshiyuki; Abdul Aziz, Noor Hafizah Binti; Gan, Tong; Hashiguchi, Hiroyuki; Yu, Tian-You; Yamamoto, Mamoru
In this paper, we describe a new digital receiver developed for a 1.3-GHz range imaging atmospheric radar. The digital receiver comprises a general-purpose software-defined radio receiver referred to as the Universal Software Radio Peripheral 2 (USRP2) and a commercial personal computer (PC). The receiver is designed to collect received signals at an intermediate frequency (IF) of 130 MHz with a sample rate of 10 MS s-1. The USRP2 digitizes IF received signals, produces IQ time series, and then transfers the IQ time series to the PC through Gigabit Ethernet. The PC receives the IQ time series, performs range sampling, carries out filtering in the range direction, decodes the phase-modulated received signals, integrates the received signals in time, and finally saves the processed data to the hard disk drive (HDD). Because only sequential data transfer from the USRP2 to the PC is available, the range sampling is triggered by transmitted pulses leaked to the receiver. For range imaging, the digital receiver performs real-time signal processing for each of the time series collected at different frequencies. Further, the receiver is able to decode phase-modulated oversampled signals. Because the program code for real-time signal processing is written in a popular programming language (C++) and widely used libraries, the signal processing is easy to implement, reconfigure, and reuse. From radar experiments using a 1-μs subpulse width and 1-MHz frequency span (i.e., 2-MHz frequency bandwidth), we demonstrate that range imaging in combination with oversampling, which was implemented for the first time by the digital receiver, is able to resolve the fine-scale structure of turbulence with a vertical scale as small as 100 m or finer.
Smith Charles M.
This work was performed under NASA's Verification and Validation (V&V) Program as an independent check of data supplied by EarthWatch, Incorporated, through the Earth Science Enterprise Scientific Data Purchase (SDP) Program. This document serves as the basis of reporting results associated with validation of orthorectified interferometric interferometric radar imagery and digital elevation models (DEM). This validation covers all datasets provided under the first campaign (Central America & Virginia Beach) plus three earlier missions (Indonesia, Red River: and Denver) for a total of 13 missions.
King, E. C.; Pritchard, H. D.; Smith, A. M.
We present a digital elevation model of the bed of Rutford Ice Stream, Antarctica derived from radio-echo sounding data. The data cover an 18 km × 40 km area immediately upstream of the grounding line of the ice stream. This area is of particular interest because repeated seismic surveys have shown that rapid erosion and deposition of subglacial sediments has taken place. The bed topography shows a range of different subglacial landforms including mega-scale glacial lineations, drumlins and hummocks. This dataset will form a baseline survey which, when compared to future surveys, should reveal how active subglacial landscapes change over time. These data also allow comparison between subglacial landforms in an active system with those observed in deglaciated areas in both polar regions. The dataset comprises observed ice thickness data, an interpolated bed elevation grid, observed surface elevation data and a surface elevation grid. The dataset is available at http://doi.org/269.
King, Edward C.; Pritchard, Hamish D.; Smith, Andrew M.
We present a digital elevation model of the bed of Rutford Ice Stream, Antarctica, derived from radio-echo sounding data. The data cover an 18 × 40 km area immediately upstream of the grounding line of the ice stream. This area is of particular interest because repeated seismic surveys have shown that rapid erosion and deposition of subglacial sediments has taken place. The bed topography shows a range of different subglacial landforms including mega-scale glacial lineations, drumlins and hummocks. This data set will form a baseline survey which, when compared to future surveys, should reveal how active subglacial landscapes change over time. These data also allow comparison between subglacial landforms in an active system with those observed in deglaciated areas in both polar regions. The data set comprises observed ice thickness data, an interpolated bed elevation grid, observed surface elevation data and a surface elevation grid. The data set is available at http://doi.org/269.
Altman, F. J.
A model for spatial distributions of reflectivity in storm cells was fitted to digital radar data. The data were taken with a modified WSR-57 weather radar with 2.6-km resolution. The data consisted of modified B-scan records on magnetic tape of storm cells tracked at 0 deg elevation for several hours. The MIT L-band radar with 0.8-km resolution produced cross-section data on several cells at 1/2 deg elevation intervals. The model developed uses ellipses for contours of constant effective-reflectivity factor Z with constant orientation and eccentricity within a horizontal cell cross section at a given time and elevation. The centers of the ellipses are assumed to be uniformly spaced on a straight line, with areas linearly related to log Z. All cross sections are similar at different heights (except for cell tops, bottoms, and splitting cells), especially for the highest reflectivities; wind shear causes some translation and rotation between levels. Goodness-of-fit measures and parameters of interest for 204 ellipses are considered.
Berkun, Andrew; Andraka, Ray
High-performance digital electronic circuits for onboard processing of return signals in an airborne precipitation- measuring radar system have been implemented in commercially available field-programmable gate arrays (FPGAs). Previously, it was standard practice to downlink the radar-return data to a ground station for postprocessing a costly practice that prevents the nearly-real-time use of the data for automated targeting. In principle, the onboard processing could be performed by a system of about 20 personal- computer-type microprocessors; relative to such a system, the present FPGA-based processor is much smaller and consumes much less power. Alternatively, the onboard processing could be performed by an application-specific integrated circuit (ASIC), but in comparison with an ASIC implementation, the present FPGA implementation offers the advantages of (1) greater flexibility for research applications like the present one and (2) lower cost in the small production volumes typical of research applications. The generation and processing of signals in the airborne precipitation measuring radar system in question involves the following especially notable steps: The system utilizes a total of four channels two carrier frequencies and two polarizations at each frequency. The system uses pulse compression: that is, the transmitted pulse is spread out in time and the received echo of the pulse is processed with a matched filter to despread it. The return signal is band-limited and digitally demodulated to a complex baseband signal that, for each pulse, comprises a large number of samples. Each complex pair of samples (denoted a range gate in radar terminology) is associated with a numerical index that corresponds to a specific time offset from the beginning of the radar pulse, so that each such pair represents the energy reflected from a specific range. This energy and the average echo power are computed. The phase of each range bin is compared to the previous echo
Lambot, S.; Minet, J.; Jadoon, K. Z.; Slob, E.; Vereecken, H.
Sustainable and optimal agricultural and environmental management of water and land resources particularly relies on the description and understanding of soil water distribution and dynamics at different scales. We present an advanced ground penetrating radar (GPR) method for mapping the shallow soil water content and unsaturated hydraulic properties at the field scale. The radar system is based on vector network analyzer technology, for which calibration is simple and constitutes an international standard. A directive horn antenna is used as both transmitter and receiver and operates off the ground. A full-waveform model describes accurately the radar signal, and is based on a linear system of complex transfer functions for efficiently describing electromagnetic phenomena within the antenna and its interaction with soil, and a specific solution of the three-dimensional Maxwell's equations for wave propagation in multilayered media. The soil electromagnetic properties and their vertical distribution are estimated by resorting to full-waveform inverse modeling using iterative global optimization methods. The proposed methodology has been validated for a series of model configurations of increasing complexity. The method is now routinely used for real-time mapping of soil surface water content and reconstruct a few number of shallow soil layers. For more complex configurations, it is necessary to regularize the inverse problem. We have shown that constraining radar data inversion using soil hydrodynamic modeling has the potential to reconstruct time-lapse, continuously variable, vertical soil water content profiles and identify the shallow unsaturated hydraulic properties. The proposed approach shows great promise for quantitative imaging of the soil properties at the field scale. The technique will be combined with electromagnetic induction in a mechanistic data fusion framework to further extend its capabilities in a digital soil mapping context.
Torres, R.; Mouginis-Mark, P.; Garbeil, H.; Bautista, L.; Ramos, E.
Central Luzon Island (13-16°N, 120-122°E), which is bounded to the east by Philippine Trench, to the west by Manila Trench, to the north by Digdig-Dingalan Fault (DDF) and to the south by Verde Island Passage Fault (VIPF), is one of the most seismically and volcanologically active regions in the Philippines. Active seismicity and violent earthquakes in the region are evidently related to the activities along the subduction zones and branches of the Philippine Fault system. Volcanic eruptions and periodic swarms of volcanic earthquakes were also observed in three active volcanoes, i.e., Pinatubo, Taal Volcano Island and Banahaw, while young calderas of Taal and Laguna de Bay are demonstrably fault-bounded. We use the Shuttle Radar Topography Mission (SRTM) data with 90 m spatial resolution to conduct regional mapping of the faults and volcanic structures in this region. Of particular interests are the NE-SW set of normal faults within the Macolod Corridor, the right-lateral Marikina Valley Fault System (MVFS), the prevalence of N-S trending structures and the series of NW-SE structures that parallel to sub-parallel the active branches of the Philippine Fault. Using ENVI software package, we processed the SRTM data into shaded relief images and examined the lineament features from different azimuth directions and angles of artificial illumination. The prominent NW-SE structures in this area revealed by SRTM data were formed as sinistral shears that parallel the seismically active DDF and VIPF. The N-S trending structures, including some segments of MVFS and N-S oriented fold axes, were apparently generated by an earlier E-W compression, but recently displayed dextral movement with localized vertical component and pull-apart zones. The overprinting of recent fault kinematics on previously formed structures suggest a dramatic shift of regional stress distribution in Central Luzon. The dextral movement along MVFS and the extensional NE-SW faults within the Macolod
Esteban-Fernandez, Daniel; Peral, Eva; McWatters, Dalia; Pollard, Brian; Rodriguez, Ernesto; Hughes, Richard
Over the last two decades, several nadir profiling radar altimeters have provided our first global look at the ocean basin-scale circulation and the ocean mesoscale at wavelengths longer than 100 km. Due to sampling limitations, nadir altimetry is unable to resolve the small wavelength ocean mesoscale and sub-mesoscale that are responsible for the vertical mixing of ocean heat and gases and the dissipation of kinetic energy from large to small scales. The proposed Surface Water and Ocean Topography (SWOT) mission would be a partnership between NASA, CNES (Centre National d'Etudes Spaciales) and the Canadian Space Agency, and would have as one of its main goals the measurement of ocean topography with kilometer-scale spatial resolution and centimeter scale accuracy. In this paper, we provide an overview of all ocean error sources that would contribute to the SWOT mission.
Herzfeld, Ute C.
The central objective of this project has been the development of geostatistical methods fro mapping elevation and ice surface characteristics from satellite radar altimeter (RA) and Syntheitc Aperture Radar (SAR) data. The main results are an Atlas of elevation maps of Antarctica, from GEOSAT RA data and an Atlas from ERS-1 RA data, including a total of about 200 maps with 3 km grid resolution. Maps and digital terrain models are applied to monitor and study changes in Antarctic ice streams and glaciers, including Lambert Glacier/Amery Ice Shelf, Mertz and Ninnis Glaciers, Jutulstraumen Glacier, Fimbul Ice Shelf, Slessor Glacier, Williamson Glacier and others.
This paper describes the results of digital image analysis and techniques applied to acoustic sounder data and topographic relief in the Geyser's region. The two dimensional fast Fourier transform (2DFFT) represents the spacial variability of a photographic image. The spacial variability of topography in complex terrain can be represented in this way and insight into degree of complexity and dominating spacial wavelengths can be gained. This was performed for a 16 km square digitized topographic map of the Geyser's region with 63.5 m resolution. It was also of interest to compare facsimile recordings of acoustic sounder data to optical turbulence measurements.
Modern digital ionosondes, with both direction finding and doppler capabilities can provide large scale pictures of the Spread-F irregularity regions. A morphological framework has been developed that allows interpretation of the hf radar data. A large scale irregularity structure is found to be nightward of the dusk terminator, stationary in the solar reference frame. As the plasma moves through this foehn-wall-like structure it descends, and irregularities may be generated. Localized upwellings, or bubbles, may be produced, and they drift with the background plasma. The spread-F irregularity region is found to be best characterized as a partly cloudy sky, due to the patchiness of the substructures. 13 references, 16 figures.
Buckley, S.; Agram, P. S.; Belz, J. E.; Crippen, R. E.; Gurrola, E. M.; Hensley, S.; Kobrick, M.; Lavalle, M.; Martin, J. M.; Neumann, M.; Nguyen, Q.; Rosen, P. A.; Shimada, J.; Simard, M.; Tung, W.
NASADEM is a significant modernization of SRTM digital elevation model (DEM) data supported by the NASA MEaSUREs program. We are reprocessing the raw radar signal data using improved algorithms and incorporating ICESat and ASTER-derived DEM data unavailable during the original processing. The NASADEM products will be freely-available through the Land Processes Distributed Active Archive Center (LPDAAC) at 1-arcsecond spacing. The most significant processing improvements involve void reduction through improved phase unwrapping and using ICESat data for control. The updated unwrapping strategy now includes the use of SNAPHU for data processing patches where the unwrapped coverage from the original residue-based unwrapper falls below a coverage threshold. In North America continental processing, first experiments show the strip void area is reduced by more than 50% and the number of strip void patches is reduced by 40%. Patch boundary voids are mitigated by reprocessing with a different starting burst and merging the unwrapping results. We also updated a low-resolution elevation database to aid with unwrapping bootstrapping, retaining isolated component of unwrapped phase, and assessing the quality of the strip DEMs. We introduce a height ripple error correction to reduce artifacts in the strip elevation data. These ripples are a few meters in size with along-track spatial scales of tens of kilometers and are due to uncompensated mast motion most pronounced after Shuttle roll angle adjustment maneuvers. We developed an along-track filter utilizing differences between the SRTM heights and ICESat lidar elevation data. For a test using all data over North America, the algorithm reduced the ICESat-SRTM bias from 80 cm to 3 cm and the RMS from 5m to 4m. After merging and regridding the SRTM strip DEMs into 1x1-degree tiles, remaining voids are primarily filled with the ASTER-derived Global DEM. We use a Delta Surface Fill method to rubbersheet fill data across the void for
Abeywickrema, Ujitha; Banerjee, Partha; Kota, Akash; Lakhtakia, Akhlesh; Swiontek, Stephen E.
The analysis of fingerprints is important for biometric identification. Two-wavelength digital holographic interferometry is used to study the topography of various types of fingerprints. This topography depends on several conditions such as the temperature, time of the day, and the proportions of eccrine and sebaceous sweat. With two-wavelength holographic interferometry, surface information can be measured with a better accuracy compared to single-wavelength phase-retrieving techniques. Latent fingerprints on transparent glass, a forensically relevant substrate are first developed by the deposition of 50-1000-nm-thick columnar thin films, and then analyzed using the transmission-mode two-wavelength digital holographic technique. In this technique, a tunable Argon-ion laser (457.9 nm to 514.5 nm) is used and holograms are recorded on a CCD camera sequentially for several sets of two wavelengths. Then the phase is reconstructed for each wavelength, and the phase difference which corresponds to the synthetic wavelength (4 μm to 48 μm) is calculated. Finally, the topography is obtained by applying proper phase-unwrapping techniques to the phase difference. Interferometric setups that utilize light reflected from the surface of interest have several disadvantages such as the effect of multiple reflections as well as the effects of the tilt of the object and its shadow (for the Mach-Zehnder configuration). To overcome these drawbacks, digital holograms of fingerprints in a transmission geometry are used. An approximately in-line geometry employing a slightly tilted reference beam to facilitate separation of various diffraction orders during holographic reconstruction is employed.
Pawul, Rudolf A.
This thesis is a reference for the Advanced Application Flight Experiment (AAFE) altimeter. The transmitter and receiver subsections are described and measurements of their current state is provided. During the 1994 NASA Greenland Experiment, the altimeter experienced several hardware malfunctions. The process of returning the radar to its fully operational state is presented in detail and necessary design modifications are explained. An updated radar user's manual is included along with various circuit designs which need to be implemented. The thesis is intended to provide an incoming graduate student with a solid foundation of the fundamentals of AAFE altimeter operation.
38. Perimeter acquisition radar building room #414, digital/electrical repair shop; showing work areas available for maintenance and equipment repair - Stanley R. Mickelsen Safeguard Complex, Perimeter Acquisition Radar Building, Limited Access Area, between Limited Access Patrol Road & Service Road A, Nekoma, Cavalier County, ND
Simard, M.; Riel, Bryan; Hensley, S.; Lavalle, Marco
Radar backscatter data contain both geometric and radiometric distortions due to underlying topography and the radar viewing geometry. Our objective is to develop a radiometric correction algorithm specific to the UAVSAR system configuration that would improve retrieval of forest structure parameters. UAVSAR is an airborne Lband radar capable of repeat?pass interferometry producing images with a spatial resolution of 5m. It is characterized by an electronically steerable antenna to compensate for aircraft attitude. Thus, the computation of viewing angles (i.e. look, incidence and projection) must include aircraft attitude angles (i.e. yaw, pitch and roll) in addition to the antenna steering angle. In this presentation, we address two components of radiometric correction: area projection and vegetation reflectivity. The first correction is applied by normalization of the radar backscatter by the local ground area illuminated by the radar beam. The second is a correction due to changes in vegetation reflectivity with viewing geometry.
Rincon, Rafael; Fatoyinbo, Temilola; Osmanoglu, Batuhan; Lee, Seung-Kuk; Ranson, K. Jon; Marrero, Victor; Yeary, Mark
NASA's Next generation Digital Beamforming SAR (DBSAR-2) is a state-of-the-art airborne L-band radar developed at the NASA Goddard Space Flight Center (GSFC). The instrument builds upon the advanced architectures in NASA's DBSAR-1 and EcoSAR instruments. The new instrument employs a 16-channel radar architecture characterized by multi-mode operation, software defined waveform generation, digital beamforming, and configurable radar parameters. The instrument has been design to support several disciplines in Earth and Planetary sciences. The instrument was recently completed, and tested and calibrated in a anechoic chamber.
Fredenslund Levinsen, Joanna; Smith, Ben; Sørensen, Louise S.; Forsberg, René
When estimating elevation changes of ice-covered surfaces from radar altimetry, it is important to correct for slope-induced errors. They cause the reflecting point of the pulse to move up-slope and thus return estimates in the wrong coordinates. Slope-induced errors can be corrected for by introducing a Digital Elevation Model (DEM). In this work, such a DEM is developed for the Greenland Ice Sheet using a combination of Envisat radar and ICESat laser altimetry. If time permits, CryoSat radar altimetry will be included as well. The reference year is 2010 and the spatial resolution 2.5 x 2.5 km. This is in accordance with the results obtained in the ESA Ice Sheets CCI project showing that a 5 x 5 km grid spacing is reasonable for ice sheet-wide change detection (Levinsen et al., 2013). Separate DEMs will be created for the given data sets, and the geostatistical spatial interpolation method collocation will be used to merge them, thus adjusting for potential inter-satellite biases. The final DEM is validated with temporally and spatially agreeing airborne lidar data acquired in the NASA IceBridge and ESA CryoVex campaigns. The motivation for developing a new DEM is based on 1) large surface changes presently being observed, and mainly in margin regions, hence necessitating updated topography maps for accurately deriving and correcting surface elevation changes, and 2) although radar altimetry is subject to surface penetration of the signal into the snowpack, data is acquired continuously in time. This is not the case with e.g. ICESat, where laser altimetry data were obtained in periods of active lasers, i.e. three times a year with a 35-day repeat track. Previous DEMs e.g. have 2007 as the nominal reference year, or they are built merely from ICESat data. These have elevation errors as small as 10 cm, which is lower than for Envisat and CryoSat. The advantage of an updated DEM consisting of combined radar and laser altimetry therefore is the possibility of
Griffin, C. R.; Estes, J. M.
A modified APQ-102 sidelooking array radar (SLAR) in a B-57 aircraft test bed is used, with other optical and infrared sensors, in remote sensing of Earth surface features for various users at NASA Johnson Space Center. The video from the radar is normally recorded on photographic film and subsequently processed photographically into high resolution radar images. Using a high speed sampling (digitizing) system, the two receiver channels of cross-and co-polarized video are recorded on wideband magnetic tape along with radar and platform parameters. These data are subsequently reformatted and processed into digital synthetic aperture radar images with the image data available on magnetic tape for subsequent analysis by investigators. The system design and results obtained are described.
Huyghebaert, D. R.; Hussey, G. C.; McWilliams, K. A.; St-Maurice, J. P.
A new 50 MHz ionospheric E-region radar is currently being developed and will be operational for the summer of 2016. The radar group in the Institute of Space and Atmospheric Studies (ISAS) at the University of Saskatchewan is designing and building the radar which will be located near the university in Saskatoon, SK, Canada and will have a field of view over Wollaston Lake in northern Saskatchewan. This novel radar will simultaneously obtain high spatial and temporal resolution through the use of a bistatic setup and pulse modulation techniques. The bistatic setup allows the radar to transmit and receive continuously, while pulse modulation techniques allow for enhanced spatial resolution, only constrained by the radio bandwidth licensing available. A ten antenna array will be used on both the transmitter and receiver sides, with each antenna having an independent radio path. This enables complete digital control of the transmitted 1 kW signal at each antenna, allowing for digital beam steering and multimode broadcasting. On the receiver side the raw digitized signal will be recorded from each antenna, allowing for complete digital post-processing to be performed on the data. From the measurements provided using these modern digital radar capabilities, further insights into the physics of E-region phenomena, such as Alfvén waves propagating from the magnetosphere above and ionospheric irregularities, may be investigated.
Baru, C.; Arrowsmith, R.; Crosby, C.; Nandigam, V.; Phan, M.; Cowart, C.
OpenTopography is a cyberinfrastructure-based facility for online access to high-resolution topography and tools. The project is an outcome of the Geosciences Network (GEON) project, which was a research project funded several years ago in the US to investigate the use of cyberinfrastructure to support research and education in the geosciences. OpenTopography provides online access to large LiDAR point cloud datasets along with services for processing these data. Users are able to generate custom DEMs by invoking DEM services provided by OpenTopography with custom parameter values. Users can track the progress of their jobs, and a private myOpenTopo area retains job information and job outputs. Data available at OpenTopography are provided by a variety of data acquisition groups under joint agreements and memoranda of understanding (MoU). These include national facilities such as the National Center for Airborne Lidar Mapping, as well as local, state, and federal agencies. OpenTopography is also being designed as a hub for high-resolution topography resources. Datasets and services available at other locations can also be registered here, providing a "one-stop shop" for such information. We will describe the OpenTopography system architecture and its current set of features, including the service-oriented architecture, a job-tracking database, and social networking features. We will also describe several design and development activities underway to archive and publish datasets using digital object identifiers (DOIs); create a more flexible and scalable high-performance environment for processing of large datasets; extend support for satellite-based and terrestrial lidar as well as synthetic aperture radar (SAR) data; and create a "pluggable" infrastructure for third-party services. OpenTopography has successfully created a facility for sharing lidar data. In the next phase, we are developing a facility that will also enable equally easy and successful sharing of
Earth-based albedo maps of Mars were compared with Mariner 9 television data and ground-based radar profiles to investigate the nature of the bright and dark albedo features. Little correlation was found except at the boundaries of classical albedo features, where some topographic control is indicated. Wind-blown dust models for seasonal and secular albedo variations are supported, but it is not clear whether the fines are derived from bright or dark parent rock. Mars, like the Earth and Moon, has probably generated two distinct types of crustal material.
Ford, R. A.
This report describes the use of an array processor for real time radar signal processing. Pulse compression, range marking, and monopulse error computation are some of the functions that will be performed in the array processor for the millimeter wave ALCOR radar augmentation. Real time software design, processor architecture, and system interfaces are discussed in the report.
Kang, Edward; Jeong, Gi Seok; Choi, Yoon Young; Lee, Kwang Ho; Khademhosseini, Ali; Lee, Sang-Hoon
Heterotypic functional materials with compositional and topographical properties that vary spatiotemporally on the micro- or nanoscale are common in nature. However, fabricating such complex materials in the laboratory remains challenging. Here we describe a method to continuously create microfibres with tunable morphological, structural and chemical features using a microfluidic system consisting of a digital, programmable flow control that mimics the silk-spinning process of spiders. With this method we fabricated hydrogel microfibres coded with varying chemical composition and topography along the fibre, including gas micro-bubbles as well as nanoporous spindle-knots and joints that enabled directional water collection. We also explored the potential use of the coded microfibres for tissue engineering applications by creating multifunctional microfibres with a spatially controlled co-culture of encapsulated cells.
Bindschadler, Robert A.; Zwally, H. Jay; Major, Judith A.; Brenner, Anita C.
Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying.
Leitao, C. D.; Huang, N. E.; Parra, C. G.
Both quasi-stationary and dynamic departures from the marine geoid were successfully detected using altitude measurements from the GEOS-3 radar altimeter. The quasi-stationary departures are observed either as elevation changes in single pass profiles across the Gulf Stream or at the crowding of contour lines at the western and northern areas of topographic maps generated using altimeter data spanning one month or longer. Dynamic features such as current meandering and spawned eddies can be monitored by comparing monthly mean maps. Comparison of altimeter inferred eddies with IR detected thermal rings indicates agreement of the two techniques. Estimates of current velocity are made using derived slope estimates in conjunction with the geostrophic equation.
Terblanche, Deon Etienne
This thesis describes the development, testing and implementation of a new method to process the output from a weather radar's logarithmic receiver. The processing method, called DISPLACE, has proven to have many applications, and is computationally efficient and accurate. Its applications include the processing of digitized logarithmic receiver output in order to simulate different receiver transfer functions, the processing of multi-parameter radar measurements and the filtering of ground clutter. It facilitates the computation of CAPPI's and radar-rainfall accumulation. The thesis also deals with the upgrading of South African weather radars since about 1990 through the in-house developed radar data acquisition system and the procedures established to ensure accurate calibrations. In addition, the hydrometeorological infrastructure deployed in the Bethlehem research are is used in an integrated manner to verify data obtained using the new method. This work is well timed to address the needs that are now emerging in South Africa and clearly illustrate the role the NPRP played in reviving radar meteorology. The DISPLACE method is proving once again that the potential of conventional weather radar has not been fully exploited. It has also stimulated the interest of young technicians and scientists in the field of radar meteorology. This augurs well for the future use of weather radar in South Africa, both in the field of rainfall stimulation and as an integral part of systems designed to forecast and to help manage the effects of severe weather conditions.
Heinemann, Paul H.; Martsolf, J. David; Gerber, John F.; Smith, Daniel L.
Manually Digitized Radar (MDR) and Geostationary Operational Environmental Satellite (GOES) thermal infrared (IR) data were merged to form a higher-resolution radar/IR product than that represented by the MDR. The combination MDR/IR maps were processed into a color coded map form and disseminated on a real-time basis through a computer network to users in the Florida agricultural community.
Doerry, A. W.; Dubbert, D. F.; Tise, B. L.
High-performance radar operation, particularly Ground Moving Target Indicator (GMTI) radar modes, are very sensitive to anomalous effects of system nonlinearities. System nonlinearities generate harmonic spurs that at best degrade, and at worst generate false target detections. One significant source of nonlinear behavior is the Analog to Digital Converter (ADC). One measure of its undesired nonlinearity is its Integral Nonlinearity (INL) specification. We examine in this paper the relationship of INL to radar performance; in particular its manifestation in a range-Doppler map or image.
Naraghi, M.; Stromberg, W.; Daily, M.
Geologic analysis of radar imagery requires accurate spatial rectification to allow rock type discrimination and meaningful exploitation of multisensor data files. A procedure is described which removes distortions produced by most sources including the heretofore elusive problem of terrain induced effects. Rectified imagery is presented which displays geologic features not apparent in the distorted data.
Hussey, G. C.; Huyghebaert, D. R.; St-Maurice, J. P.; McWilliams, K. A.
A new fully digital bistatic 50-MHz VHF radar is currently being developed by the radar group in the Institute of Space and Atmospheric Studies (ISAS) at the University of Saskatchewan. This paper presents the scientific motivation for the new radar. Traditionally bistatic radars have had excellent time resolution, but were significantly lacking in range resolution. With the now available accurate timing abilities and advanced pulse modulation techniques, bistatic radar configurations with both excellent temporal and spatial resolution are able to map or 'image' the E-region. The E-region portion of the ionosphere being the base of the magnetosphere has both global (ionosphere-magnetosphere system) and local phenomena of interest. The currents in the magnetosphere close in the E-region. Field-aligned currents (FACs) and Alfven waves are phenomena with origins in the magnetosphere which present their 'signatures' in the E-region. For example, Alfven waves (produced by the Alfven wave resonator) have different time scales, from less than a Hertz to periods of tens of minutes --- and the high temporal and spatial resolution of this new digital E-region radar will be able to detect them all. The E-region is also a dynamic plasma medium with the two-steam and gradient drift instabilities present and the improved measurement abilities will give fresh physical insight.
Richardson, J.; Graves, K.; Bowling, T.
Previous studies of the combined effects of asteroid shape, spin, and self-gravity have focused primarily upon the failure limits for bodies with a variety of standard shapes, friction, and cohesion values [1,2,3]. In this study, we look in the opposite direction and utilize 22 asteroid shape-models derived from radar inversion  and 7 small body shape-models derived from spacecraft observations  to investigate the region in shape/spin space [1,2] wherein self-gravity and rotation combine to produce a stable minimum state with respect to surface potential differences, dynamic topography, slope magnitudes, and erosion rates. This erosional minimum state is self-correcting, such that changes in the body's rotation rate, either up or down, will increase slope magnitudes across the body, thereby driving up erosion rates non-linearly until the body has once again reached a stable, minimized surface state . We investigated this phenomenon in a systematic fashion using a series of synthesized, increasingly prolate spheroid shape models. Adjusting the rotation rate of each synthetic shape to minimize surface potential differences, dynamic topography, and slope magnitudes results in the magenta curve of the figure (right side), defining the zone of maximum surface stability (MSS). This MSS zone is invariant both with respect to body size (gravitational potential and rotational potential scale together with radius), and density when the scaled-spin of  is used. Within our sample of observationally derived small-body shape models, slow rotators (Group A: blue points), that are not in the maximum surface stability (MSS) zone and where gravity dominates the slopes, will generally experience moderate erosion rates (left plot) and will tend to move up and to the right in shape/spin space as the body evolves (right plot). Fast rotators (Group C: red points), that are not in the MSS zone and where spin dominates the slopes, will generally experience high erosion rates
Gesch, D.; Williams, J.; Miller, W.
Elevation models produced from Shuttle Radar Topography Mission (SRTM) data will be the most comprehensive, consistently processed, highest resolution topographic dataset ever produced for the Earth's land surface. Many applications that currently use elevation data will benefit from the increased availability of data with higher accuracy, quality, and resolution, especially in poorly mapped areas of the globe. SRTM data will be produced as seamless data, thereby avoiding many of the problems inherent in existing multi-source topographic databases. Serving as precursors to SRTM datasets, the U.S. Geological Survey (USGS) has produced and is distributing seamless elevation datasets that facilitate scientific use of elevation data over large areas. GTOPO30 is a global elevation model with a 30 arc-second resolution (approximately 1-kilometer). The National Elevation Dataset (NED) covers the United States at a resolution of 1 arc-second (approximately 30-meters). Due to their seamless format and broad area coverage, both GTOPO30 and NED represent an advance in the usability of elevation data, but each still includes artifacts from the highly variable source data used to produce them. The consistent source data and processing approach for SRTM data will result in elevation products that will be a significant addition to the current availability of seamless datasets, specifically for many areas outside the U.S. One application that demonstrates some advantages that may be realized with SRTM data is delineation of land surface drainage features (watersheds and stream channels). Seamless distribution of elevation data in which a user interactively specifies the area of interest and order parameters via a map server is already being successfully demonstrated with existing USGS datasets. Such an approach for distributing SRTM data is ideal for a dataset that undoubtedly will be of very high interest to the spatial data user community.
Maurer, H. E.; Oderman, W.; Crosswell, W. F.
A data set is described which consists of digitized synthetic aperture radar (SAR) imagery plus correlative data and some preliminary analysis results. This data set should be of value to experimenters who are interested in the SAR instrument and its application to the detection and monitoring of oil on water and other distributed targets.
Heidelbach, Robert; Bolus, R.; Chadwick, J.
Digital Terrain Elevations (DTE) that can be rapidly generated, and that have better fidelity and accuracy than Digital Terrain Elevation Data (DTED) Levels 1 or 2, would be extremely beneficial to Department of Defense (DOD) military operations, civil works programs, and various commercial applications. As a result, the Advanced Research Projects Agency (ARPA), along with the U.S. Army Topographic Engineering Center (TEC), are developing an Interferometric Synthetic Aperture Radar (IFSAR) elevation mapping capability. This system, the Interferometric Synthetic Aperture Radar for Digital Radar Elevations (IFSARE), is capable of collecting and providing data in all weather (reasonable), in day or night scenarios, and where obscurants are present. The IFSARE, which is currently undergoing Integration and Test, will allow for rapid on-line automatic processing of the collected digital radar data into DTE and high quality imagery. The prime contractor is the Environmental Research Institute of Michigan (ERIM). This paper addresses the proof of concept for civil works applications by analyzing a data set taken by the Wright Labs/ERIM Data Collection System (DCS). The objective was to demonstrate the capability of an IFSAR system to provide high fidelity, fine resolution DTE that can be employed in hydraulic models of the Mississippi River watershed. The demonstration was sponsored by ARPA and TEC.
A digital signal processing system was studied for the determination of the spectral frequency distribution of echo signals from a teleoperator radar system. The system consisted of a sample and hold circuit, an analog to digital converter, a digital filter, and a Fast Fourier Transform. The system is interfaced to a 16 bit microprocessor. The microprocessor is programmed to control the complete digital signal processing. The digital filtering and Fast Fourier Transform functions are implemented by a S2815 digital filter/utility peripheral chip and a S2814A Fast Fourier Transform chip. The S2815 initially simulates a low-pass Butterworth filter with later expansion to complete filter circuit (bandpass and highpass) synthesizing.
Yasodha, Polisetti; Jayaraman, Achuthan; Thriveni, A.
Digital receiver extracts the received echo signal information, and is a potential subsystem for atmospheric radar, also referred to as wind profiling radar (WPR), which provides the vertical profiles of 3-dimensional wind vector in the atmosphere. This paper presents the development of digital receiver using COTS board based Software Defined Radio technique, which can be used for atmospheric radars. The developmental work is being carried out at National Atmospheric Research Laboratory (NARL), Gadanki. The digital receiver consists of a commercially available software defined radio (SDR) board called as universal software radio peripheral B210 (USRP B210) and a personal computer. USRP B210 operates over a wider frequency range from 70 MHz to 6 GHz and hence can be used for variety of radars like Doppler weather radars operating in S/C bands, in addition to wind profiling radars operating in VHF, UHF and L bands. Due to the flexibility and re-configurability of SDR, where the component functionalities are implemented in software, it is easy to modify the software to receive the echoes and process them as per the requirement suitable for the type of the radar intended. Hence, USRP B210 board along with the computer forms a versatile digital receiver from 70 MHz to 6 GHz. It has an inbuilt direct conversion transceiver with two transmit and two receive channels, which can be operated in fully coherent 2x2 MIMO fashion and thus it can be used as a two channel receiver. Multiple USRP B210 boards can be synchronized using the pulse per second (PPS) input provided on the board, to configure multi-channel digital receiver system. RF gain of the transceiver can be varied from 0 to 70 dB. The board can be controlled from the computer via USB 3.0 interface through USRP hardware driver (UHD), which is an open source cross platform driver. The USRP B210 board is connected to the personal computer through USB 3.0. Reference (10 MHz) clock signal from the radar master oscillator
Desai, Nilesh; Vachhani, J. G.; Soin, Sumit; Agrawal, Rinku; Rao, C. V. N.; Gujraty, Virendra; Rana, Surindersingh
Technology development related to digital, antenna and RF subsystems for Microwave Radar Sensors like Synthetic Aperture Radar, Scatterometer, Altimeter and Radiometer is one of the major activities under ISRO's microwave remote sensing programme, since 1980s. These technologies are now being gainfully utilized for building ISRO's operational Earth Observation missions involving microwave sensors like Radar Imaging Satellite, RISAT SAR, Oceansat-2 Scatterometer, Megha-Tropiques, MADRAS and Airborne SAR for Disaster Management, DMSAR. Concurrently, advanced technology developments in these fields are underway to meet the major technological challenges of building ISRO's proposed advanced microwave missions like ultra-high resolution SAR's, Synthetic Aperture Radiometer (SARAD), Milli-meter and sub-millimeter wave sounders and SAR Constellations for Disaster management as well as Interferometric, Polarmetric and polarmetric interferometry applications. Also, these hardware are being designed with core radar electronics concept, in which the same RF and digital hardware sub-units / modules will be utilized to build different microwave radar sensors. One of the major and common requirements for all these active and passive microwave sensors is the moderate to highspeed data acquisition and signal processing system. Traditionally, the Data acquisition units for all these radar sensors are implemented as stand-alone units, following the radar receivers. For ISRO's C-band airborne SAR (ASAR) and RISAT high resolution SAR, we have designed and developed High Speed 8-bit ADC based I/Q Digitisers, operating at 30.814 MHz and 250 MHz sampling rates, respectively. With the increasing demand of wide bandwidth and ultra-high resolution in imaging and non-imaging radar systems, the technology trend worldwide is towards a digital receiver, involving bandpass or IF sampling, thus eliminating the need for RF down converters and analog IQ demodulators. In order to evolve a generic
Carter, L. M.; Rincon, R. F.
Many important questions in planetary science depends on our ability to detect and map surface and subsurface layers of planetary bodies. We are developing a P-band (435 MHz, 70 cm wavelength) digital beamforming radar, called Space Exploration SAR (SESAR), capable of providing the measurement flexibility needed to address multiple types of science goals. SESAR will provide high spatial resolution imaging, full polarimetry, multibeam scatterometry and altimetry of planetary targets such as the Moon and Mars by using beamforming technology that can adjust the radar experiment to meet the specific science goals of each target.
Foley, D. . Dept. of Earth Sciences); McEwen, A.; Duffield, W. ); Heiken, G. )
The authors propose, based on reconnaissance geology studies and interpretation of landforms as depicted by Landsat Thematic Mapper (TM) images combined with digitized topography, that the Quezaltenango basin of Guatemala is part of a caldera. The Quezaltenango basin is an elliptical depression, about 12 by 25 km and about 500 m deep. The proposed Xela Caldera extends beyond the basin more than 10 km to the north. The geomorphological features of the area that are typical of a geologically young large-scale caldera include bounding walls that have steep interior and gentle exterior slopes; broad flat areas at the base of the walls; at least one large block, about 3 by 12 km, that only partly floundered as the caldera collapsed; resurgence of a younger volcanic dome, flow and small-scale caldera complex (last active in 1818); younger volcanoes located along the structural margin of the major caldera (one of which is currently active) lobate features on the caldera margins that may indicate a multiple sequence of eruptions; and an active, high-temperature geothermal system. The valley is coincident with a gravity low. Extensive ash-flow tuff sheets that have no identified source are located north of the caldera, and may be the outflow deposits. The Xela caldera is similar in size to the Atitlan caldera, which lies about 50 km southeast of Quezaltenango. The Xela Caldera, if confirmed by future studies, may contain undiscovered geothermal resources, may present a significant geologic hazard to the more than 400,000 people who occupy the Quezaltenango valley, and may be a new member of the list of magmatic systems that have the capability to change global climate for several years.
Oberg, J. M.; Ulaby, F. T.
The design of the MAS 2-8 (2 to 8 GHz microwave-active spectrometer), a ground-based sensor system, is presented. A major modification in 1974 to the MAS 2-8, that of a control subsystem to automate the data-taking operation, is the prime focus. The digital control unit automatically changes all system parameters except FM rate and records the return signal on paper tape. The overall system operation and a detailed discussion of the design and operation of the digital control unit are presented.
Oh, T. J.
A new topography generation tool utilizing spectral transformation technique for both structured and unstructured grids is presented. For the source global digital elevation data, the NASA Shuttle Radar Topography Mission (SRTM) 15 arc-second dataset (gap-filling by Jonathan de Ferranti) is used and for land/water mask source, the NASA Moderate Resolution Imaging Spectroradiometer (MODIS) 30 arc-second land water mask dataset v5 is used. The original source data is coarsened to a intermediate global 2 minute lat-lon mesh. Then, spectral transformation to the wave space and inverse transformation with wavenumber truncation is performed for isotropic topography smoothness control. Target grid topography mapping is done by bivariate cubic spline interpolation from the truncated 2 minute lat-lon topography. Gibbs phenomenon in the water region can be removed by overwriting ocean masked target coordinate grids with interpolated values from the intermediate 2 minute grid. Finally, a weak smoothing operator is applied on the target grid to minimize the land/water surface height discontinuity that might have been introduced by the Gibbs oscillation removal procedure. Overall, the new topography generation approach provides spectrally-derived, smooth topography with isotropic resolution and minimum damping, enabling realistic topography forcing in the numerical model. Topography is generated for the cubed-sphere grid and tested on the KIAPS Integrated Model (KIM).
Anderson, K. A.
Papers are presented which were published as a result of a project involving the preparation of a topographical elevation contour map of Mars from all data sources available through 1969, as well as the observation of Mars by spectroscopic methods in 1971 to provide additional pressure data for topographic information. Topics of the papers include: the analysis of large-scale Martian topography variations - data preparation from earth based radar, earth based CO2 spectroscopy, and Mariners 6 and 7 CO2 spectroscopy; the analysis of water content in observed Martian white clouds; and Martian, lunar, and terrestrial crusts - a three-dimensional exercise in comparative geophysics.
Frankel, K. L.; Pazzaglia, F. J.
Mountain ranges in the southern Rocky Mountains have departed on unique landscape evolutionary pathways in the late Cenozoic that are directly dependent upon the degree of post-orogenic tectonic activity they have experienced. The topography of Sierra Nacimiento, a Laramide uplift in west-central New Mexico lacking an active range-front fault, is shaped primarily by erosional exhumation that is continuous, but not steady, being driven by distal base level fall from Rio Grande incision and resultant south to north knickpoint migration. In contrast, the topography of the Taos Range, a rift flank uplift in north-central New Mexico is shaped by contrasting active stream incision and aggradation astride an active range front normal fault. The distinction between exhumation-dominated and tectonically-dominated mountain fronts is best quantified by analyses of a new metric we call the drainage basin volume to drainage basin area ratio (V-A ratio) as well as the gradients of first-order streams. Drainage basin volume and area are calculated by constructing topographic envelope maps from 10 m resolution digital elevation models (DEM). The envelope maps are pinned by the watershed divide and cover the maximum elevations in each drainage basin. Subtracting the original DEM from the maximum elevation envelope map produces a topographic residual map from which area and volume data can be obtained. The erosionally exhumed Sierra Nacimiento has a mean V-A ratio of 88 m while the tectonically active Taos Range has a mean V-A ratio of 140 m. Similarly, there are systematic differences in the gradients of first order streams measured both in the range block and approximately 5 km of adjacent piedmont. Streams were defined and subsequently Strahler ordered by a flow accumulation threshold of 250 water-equivalent grid cell units. First order stream channel long profiles were extracted from the DEM at 30 meter increments and gradients were calculated by a FORTRAN program. Gradients of
This paper describes a digital processing algorithm and its associated system design for producing images from Synthetic Aperture Radar (SAR) data. The proposed system uses the Fast Fourier Transform (FFT) approach to perform the two-dimensional correlation process. The range migration problem, which is often a major obstacle to efficient processing, can be alleviated by approximating the locus of echoes from a point target by several linear segments. SAR data corresponding to each segment is correlated separately, and the results are coherently summed to produce full-resolution images. This processing approach exhibits greatly improved computation efficiency relative to conventional digital processing methods.
Hall, John K.
The Dead Sea Depression occupies the central part of that segment of the Syrian-African rift extending from the mountains of southern Lebanon to the Arava Valley north of the Gulf of Aqaba on the Red Sea. The lowest continental feature in the world, it encloses the Dead Sea, whose surface now lies at about -409 m. This shrunken sea consists of a shallow southern basin occupied by evaporation pans, and a northern basin, 50 km in length by 13-17 km in width, with a flat floor some -731 m below MSL. The depression is asymmetric, with the mountains to the east being up to 500 m higher. The closed depression has two outlets, one at 60.5 m altitude into the Mediterranean via the Jezreel Valley near Afula, and a second into the Red Sea at about 200 m altitude in the Arava Valley. A digital terrain model (DTM), based upon the local 1:50,000 scale topographic maps and modern bathymetric surveys, gives elevations to decimeter resolution on a grid with 25 m spacing. When visualized as shaded relief or in coloured hypsometric format, the DTM graphically shows the overall morphology and its many features. These include the local tectonics, from the major boundary faults to the associated subsidiary faults and lineaments, and their possible relationships with the incised canyons offshore. Also clearly visible are the volcanic cones on the Golan heights, the coast-paralleling kurkar ridges, and the halokinetic disruptions to the deep offshore topography. The DTM was used to compute a hypsometric curve for the closed depression from -731 m to +60 m. The maximum area and volume of a lake extending to the Afula 'spillway' would be 5985 km 2 and 1602 km 3, respectively. Inflection points were observed at -715 m, where the Dead Sea's flat bottom meets its steep flanks, at -385 m where the northern basin meets the flat floor of the southern basin and surrounding plain, at -228 with the step up to the basin of the Sea of Galilee, and at Mediterranean sea level.
Parsons, Reid; Holt, John
Lobate debris aprons (LDAs) are midlatitude deposits of debris-covered ice formed during one or more periods of glaciation during the Amazonian period. However, little is known about the climate conditions that led to LDA formation. We explore a hypothesis in which a single, extended period of precipitation of ice on the steep slopes of Euripus Mons (45°S, 105°E—east of the Hellas Basin) produced a flowing ice deposit which was protected from subsequent ablation to produce the LDA found at this location. We test this hypothesis with a numerical ice flow model using an ice rheology based on low-temperature ice deformation experiments. The model simulates ice accumulation and flow for the northern and southern lobes of the Euripus Mons LDA using basal topography constrained by data from the Shallow Radar (SHARAD) and a range of ice viscosities (determined by ice temperature and ice grain size). Simulations for the northern lobe of the Euripus LDA produce good fits to the surface topography. Assuming an LDA age of ˜60 Myr and an expected temperature range of 200 to 204 K (for various obliquities) gives an ice grain size of ≈2 mm. Simulations of the southern section produce poor fits to surface topography and result in much faster flow timescales unless multiple ice deposition events or higher ice viscosities are considered.
Olson, Richard F.; Braselton, William J.; Mohlere, Richard D.
Matched filter processing for pulse compression of phase coded waveforms is a classic method for increasing radar range measurement resolution. A generic approach for simulating high resolution range extended radar scenes in a Hardware in the Loop (HWIL) test environment is to pass the phase coded radar transmit pulse through an RF tapped delay line comprised of individually amplitude- and phase-weighted output taps. In the generic approach, the taps are closely spaced relative to time intervals equivalent to the range resolution of the compressed radar pulse. For a range-extended high resolution clutter scene, the increased number of these taps can make an analog implementation of an RF tapped delay system impractical. Engineers at the U.S. Army Aviation and Missile Research, Development and Engineering Center (AMRDEC) have addressed this problem by transferring RF tapped delay line signal operations to the digital domain. New digital tapped delay line (DTDL) systems have been designed and demonstrated which are physically compact compared to analog RF TDLs, leverage low cost FPGA and data converter technology, and may be readily expanded using open slots in a VME card cage. In initial HWIL applications, the new DTDLs have been shown to produce better dynamic range in pulse compressed range profiles than their analog TDL predecessors. This paper describes the signal requirements and system architecture for digital tapped delay lines. Implementation, performance, and HWIL simulation integration issues for AMRDEC's first generation DTDLs are addressed. The paper concludes with future requirements and plans for ongoing DTDL technology development at AMRDEC.
Levinsen, J. F.; Smith, B. E.; Sandberg Sørensen, L.; Khvorostovsky, K.; Forsberg, R.
With the launch of the first radar altimeter by ESA in 1992, more than two decades of radar altimetry data are now available. Therefore, one goal of ESA's Ice Sheet Climate Change Initiative is the estimation of surface elevation changes of the Greenland Ice Sheet (GrIS) based on ERS-1, -2, Envisat, CryoSat-2, and, in the longer term, Sentinel-3 data. This will create a data record from 1992 until present date. In addition to elevation-change records, such data can be processed to produce digital elevation models, or DEMs, of the ice sheets. The DEMs can be used to correct radar altimetry data for slope-induced errors resulting from the large footprint (e.g. 2-10 km for Envisat vs. 60 m for ICESat laser altimetry) or to correct for the underlying surface topography when applying the repeat-track method. DEMs also provide key information in e.g. SAR remote sensing of ice velocities to remove the interferograms' topographic signal or in regional climate modeling. This work focuses on the development of a GrIS DEM from Envisat and CryoSat-2 altimetry, corrected with temporally and spatially coincident NASA ICESat, ATM, and LVIS laser data. The spatial resolution is 2 x 2 km and the reference year 2010. It is based on 2009 and 2010 data, the 2009 data adjusted to 2010 by accounting for the intermediate elevation changes. This increases the spatial data coverage and reduces data errors. The GIMP DEM has been corrected for negative elevations and errors in the north, and used to constrain the final DEM. The recently acquired observations and increased data coverage give a strong advantage to this DEM relative to previous models, based on lower-resolution, more temporally scattered data (e.g. a decade of observations or only ICESat data, limited to three annual 35-day acquisition periods). Furthermore, as surface changes occur continuously, an up-to-date DEM is necessary to correctly constrain the observations, thereby ensuring an accurate change detection or modeling
Salamunićcar, Goran; Lončarić, Sven
In our previous work, in order to extend the GT-57633 catalogue [PSS, 56 (15), 1992-2008] with still uncatalogued impact-craters, the following has been done [GRS, 48 (5), in press, doi:10.1109/TGRS.2009.2037750]: (1) the crater detection algorithm (CDA) based on digital elevation model (DEM) was developed; (2) using 1/128° MOLA data, this CDA proposed 414631 crater-candidates; (3) each crater-candidate was analyzed manually; and (4) 57592 were confirmed as correct detections. The resulting GT-115225 catalog is the significant result of this effort. However, to check such a large number of crater-candidates manually was a demanding task. This was the main motivation for work on improvement of the CDA in order to provide better classification of craters as true and false detections. To achieve this, we extended the CDA with the machine learning capability, using support vector machines (SVM). In the first step, the CDA (re)calculates numerous terrain morphometric attributes from DEM. For this purpose, already existing modules of the CDA from our previous work were reused in order to be capable to prepare these attributes. In addition, new attributes were introduced such as ellipse eccentricity and tilt. For machine learning purpose, the CDA is additionally extended to provide 2-D topography-profile and 3-D shape for each crater-candidate. The latter two are a performance problem because of the large number of crater-candidates in combination with the large number of attributes. As a solution, we developed a CDA architecture wherein it is possible to combine the SVM with a radial basis function (RBF) or any other kernel (for initial set of attributes), with the SVM with linear kernel (for the cases when 2-D and 3-D data are included as well). Another challenge is that, in addition to diversity of possible crater types, there are numerous morphological differences between the smallest (mostly very circular bowl-shaped craters) and the largest (multi-ring) impact
Schmidt, G.; Ruster, R.; Czechowsky, P.
The SOUSY-VHF-Radar operates at a frequency of 53.5 MHz in a valley in the Harz mountains, Germany, 90 km from Hanover. The radar controller, which is programmed by a 16-bit computer holds 1024 program steps in core and controls, via 8 channels, the whole radar system: in particular the master oscillator, the transmitter, the transmit-receive-switch, the receiver, the analog to digital converter, and the hardware adder. The high-sensitivity receiver has a dynamic range of 70 dB and a video bandwidth of 1 MHz. Phase coding schemes are applied, in particular for investigations at mesospheric heights, in order to carry out measurements with the maximum duty cycle and the maximum height resolution. The computer takes the data from the adder to store it in magnetic tape or disc. The radar controller is programmed by the computer using simple FORTRAN IV statements. After the program has been loaded and the computer has started the radar controller, it runs automatically, stopping at the program end. In case of errors or failures occurring during the radar operation, the radar controller is shut off caused either by a safety circuit or by a power failure circuit or by a parity check system.
Doerry, Armin Walter; Dubbert, Dale F.; Tise, Bertice L.
Radar operation, particularly Ground Moving Target Indicator (GMTI) radar modes, are very sensitive to anomalous effects of system nonlinearities. These throw off harmonic spurs that are sometimes detected as false alarms. One significant source of nonlinear behavior is the Analog to Digital Converter (ADC). One measure of its undesired nonlinearity is its Integral Nonlinearity (INL) specification. We examine in this report the rela tionship of INL to GMTI performance. - 4 - Acknowledgements This report is the result of a n unfunded Research and Development effort . Sandia National Laboratories is a multi - program laboratory managed and operated by Sandia Corporation, a wholly owned subsidia ry of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE - AC04 - 94AL85000.
Chirico, Peter G.; Warner, Michael B.
EXPLANATION The digital elevation model of Ecuador represented in this data set was produced from over 40 individual tiles of elevation data from the Shuttle Radar Topography Mission (SRTM). Each tile was downloaded, converted from its native Height file format (.hgt), and imported into a geographic information system (GIS) for additional processing. Processing of the data included data gap filling, mosaicking, and re-projection of the tiles to form one single seamless digital elevation model. For 11 days in February of 2000, NASA, the National Geospatial-Intelligence Agency (NGA), the German Aerospace Center (DLR), and the Italian Space Agency (ASI) flew X-band and C-band radar interferometry onboard the Space Shuttle Endeavor. The mission covered the Earth between 60?N and 57?S and will provide interferometric digital elevation models (DEMs) of approximately 80% of the Earth's land mass when processing is complete. The radar-pointing angle was approximately 55? at scene center. Ascending and descending orbital passes generated multiple interferometric data scenes for nearly all areas. Up to eight passes of data were merged to form the final processed SRTM DEMs. The effect of merging scenes averages elevation values recorded in coincident scenes and reduces, but does not completely eliminate, the amount of area with layover and terrain shadow effects. The most significant form of data processing for the Ecuador DEM was gap-filling areas where the SRTM data contained a data void. These void areas are a result of radar shadow, layover, standing water, and other effects of terrain, as well as technical radar interferometry phase unwrapping issues. To fill these gaps, topographic contours were digitized from 1:50,000 - scale topographic maps which date from the mid-late 1980's (Souris, 2001). Digital contours were gridded to form elevation models for void areas and subsequently were merged with the SRTM data through GIS and remote sensing image-processing techniques
Kütter, Sissy; Franke-Börner, Antje; Börner, Ralph-Uwe; Spitzer, Klaus
Marine volcanoes are particularly demanding when it comes to applying electric or electromagnetic methods to investigate their interiors. First, the surrounding highly conductive sea water represents a significant difference in conductivity with respect to the volcanic edifice, second, the volcano's topography has great impact on the electromagnetic response, and, third, the surrounding sea bed topography heavily distorts electromagnetic fields in frequency bands that interfere with a certain spatial wavelength and amplitude of the bathymetry. By neglecting these issues severe misinterpretations are the inevitable consequence. We present different approaches to 3D vector finite element simulation on unstructured grids which are able to compute plain-wave magnetotelluric fields for models including arbitrary surface and sea bed topography. As an example, we consider Stromboli volcano. One major issue is the incorporation of the Stromboli topography using a digital terrain model so that nearly all geometric features affecting the electromagnetic response are considered and an electromagnetic view on Stromboli's interior becomes possible. By carrying out a number of different synthetic experiments it has become obvious that not only the topography of Stromboli island itself is influencing the behavior of the fields but, even stronger, the topography of the surrounding sea bed within a radius of several tens of kilometers. The experiment therefore comprises three steps which gradually approach the complex setting of the target and map the entire volcanic environment with increasing accuracy. The first step outlines the volcano as a resistive geometric frustum surrounded by conductive sea water and underlain by a resistive substratum. This model already gives fundamental answers concerning the principal frequency-dependent current flow pattern within the edifice and the surrounding sea. For this purpose, the MT response was calculated at the earth/sea and the earth
Clark, W. L.; Vanzandt, T. E.; Gage, K. S.; Einaudi, F. E.; Rottman, J. W.; Hollinger, S. E.
The paper describes a six-station digital-barometer network centered on the Flatland ST radar to support observational studies of gravity waves and other mesoscale features at the Flatland Atmospheric Observatory in central Illinois. The network's current mode of operation is examined, and a preliminary example of an apparent group of waves evident throughout the network as well as throughout the troposphere is presented. Preliminary results demonstrate the capabilities of the current operational system to study wave convection, wave-front, and other coherent mesoscale interactions and processes throughout the troposphere. Unfiltered traces for the pressure and horizontal zonal wind, for days 351 to 353 UT, 1990, are illustrated.
Masuoka, Penny M.; Harris, Jeff; Lowman, Paul D., Jr.; Blodget, Herbert W.
Various digital enhancement techniques for SAR are compared using SIR-B and Seasat images of the Canadian Shield. The three best methods for enhancing geological structure were found to be: (1) a simple linear contrast stretch; (2) a mean or median low-pass filter to reduce speckle prior to edge enhancement or a K nearest-neighbor average to cosmetically reduce speckle; and (3) a modification of the Moore-Waltz (1983) technique. Three look directions were coregistered and several means of data display were investigated as means of compensating for radar azimuth biasing.
Space-shuttle radar topography data from central Sulawesi, Indonesia, reveal two corrugated, domal landforms, covering hundreds to thousands of square kilometers, that are bounded to the north by an abrupt transition to typical hilly to mountainous topography. These domal landforms are readily interpreted as metamorphic core complexes, an interpretation consistent with a single previous field study, and the abrupt northward transition in topographic style is interpreted as marking the trace of two extensional detachment faults that are active or were recently active. Fault dip, as determined by the slope of exhumed fault footwalls, ranges from 4?? to 18??. Application of critical-taper theory to fault dip and hanging-wall surface slope, and to similar data from several other active or recently active core complexes, suggests a theoretical limit of three degrees for detachment-fault dip. This result appears to conflict with the dearth of seismological evidence for slip on faults dipping less than ~. 30??. The convex-upward form of the gently dipping fault footwalls, however, allows for greater fault dip at depths of earthquake initiation and dominant energy release. Thus, there may be no conflict between seismological and mapping studies for this class of faults. ?? 2011 Elsevier B.V.
Brown Jr., C G; Sarabandi, K; Pierce, L E
In this paper, accurate tree stand height retrieval is demonstrated using C-band Shuttle Radar Topography Mission (SRTM) height and ancillary data. The tree height retrieval algorithm is based on modeling uniform tree stands with a single layer of randomly oriented vegetation particles. For such scattering media, the scattering phase center height, as measured by SRTM, is a function of tree height, incidence angle, and the extinction coefficient of the medium. The extinction coefficient for uniform tree stands is calculated as a function of tree height and density using allometric equations and a fractal tree model. The accuracy of the proposed algorithm is demonstrated using SRTM and TOPSAR data for 15 red pine and Austrian pine stands (TOPSAR is an airborne interferometric synthetic aperture radar). The algorithm yields root-mean-square (rms) errors of 2.5-3.6 m, which is a substantial improvement over the 6.8-8.3-m rms errors from the raw SRTM minus National Elevation Dataset Heights.
Rial, F. I.; Mendez-Rial, Roi; Lawadka, Lukasz; Gonzalez-Huici, Maria A.
In this paper we show how radio frequency interference (RFI) generated by digital video broadcasting terrestrial and digital audio broadcasting transmitters can be an important noise source for forward-looking ground penetrating radar (FLGPR) systems. Even in remote locations the average interference power sometimes exceeds ultra-wideband signals by many dB, becoming the limiting factor in the system sensitivity. The overall problem of RFI and its impact in GPR systems is briefly described and several signal processing approaches to removal of RFI are discussed. These include spectral estimation and coherent subtraction algorithms and various filter approaches which have been developed and applied by the research community in similar contexts. We evaluate the performance of these methods by simulating two different scenarios submitted to real RFI acquired with a FLGPR system developed at the Fraunhofer Institute for High Frequency Physics and Radar Techniques (FHR), (GER). The effectiveness of these algorithms in removing RFI is presented using some performance indices after suppression.
Cook, A. C.; Spudis, P. D.; Robinson, M. S.; Watters, T. R.
Planet-wide (1 km/pixel and 5 km/pixel) Digital Elevation Models (DEM) of the Moon have been produced using Clementine UVVIS (Ultraviolet-Visible) stereo. Six new basins have been discovered, two suspected basins have been confirmed, and the dimensions of existing basins better defined. Additional information is contained in the original extended abstract.
Delikaraoglou, D.; Mintourakis, I.; Kallianou, F.
With the realization of the Shuttle Radar Topographic Mission (SRTM) and the free distribution of its global elevation dataset with 3 arcsec (90 m) resolution and less than 16 m vertical accuracy, together with the availability of the higher resolution (30 m) and accuracy (10 m) Digital Terrain Models (DTM) from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), these two valuable sources of uniform DEM data represent a revolution in the world of terrain modelling. DEMs are an important source of data for the generation of high resolution geoids since they provide the high-frequency content of the gravity field spectrum and are suitable for the computation of terrain effects to gravity and indirect effects to the geoid, thus allowing the combination of global geopotential models, local gravity anomalies and information about the earth's topography (represented by a given DEM). However, although such models are available for land, there are no readily accessible Digital Bathymetry Models (DBMs) of equivalent quality for the coastal and oceanic regions. Most of the global DBM's (e.g. ETOPO1, SRTM30, and GEBCO global bathymetric grid) are compilations of heterogeneous data with medium resolution and accuracy. This prevents to exploit the potential of the recent high resolution (1 arcmin) marine free-air gravity anomalies datasets derived from satellite altimetry (such as the DNSC08, and the Sandwell & Smith v18.1 (S&Sv18.1) global solutions) in conjunction with such global DBM's. Fortunately, for some regions, recently have become available DBM's of much better accuracy and resolution, such as the DBM of 1 km resolution for many regions of the Mediterranean Sea which is distributed by IFREMER, the French Research Institute for Exploitation of the Sea. The scope of this study is to use this latest regional DBM in combination with the newly available DNSC08 and SSV18.1 global marine free-air gravity anomalies datasets for marine and near shore
Sadowy, Gregory; Ghaemi, Hirad; Heavy, Brandon; Perkovic, Dragana; Quddus, Momin; Zawadzki, Mark; Moller, Delwyn
GLISTIN is an instrument concept for a single-pass interferometric SAR operating at 35.6 GHz. To achieve large swath widths using practical levels of transmitter power, a digitally-beamformed planar waveguide array is used. This paper describes results from a ground-based demonstration of a 16-receiver prototype. Furthermore, SweepSAR is emerging as promising technique for achieving very wide swaths for surface change detection. NASA and DLR are studying this approach for the DESDynI and Tandem-L missions. SweepSAR employs a reflector with a digitally-beamformed array feed. We will describe development of an airborne demonstration of SweepSAR using the GLISTIN receiver array and a reflector.
Corbett, Francis J.; Groden, Michael; Dryden, Gordon L.; Pfeiffer, George; Boos, Robert; Youmans, Douglas G.
Textron has designed and built a high-powered CO2 laser radar for long range targeting and remote sensing. This is a coherent, multi-wavelength system with a 2D, wide-band image processing capability. The digital processor produces several output products from the transmitter return signals including range, velocity, angle, and 2D range-Doppler images of hard-body targets (LADAR mode). In addition, the processor sorts and reports on data acquired from gaseous targets by wavelength and integrated path absorption (LIDAR mode). The digital processor has been developed from commercial components with a SUN SPARC 20 serving as the operator workstation and display. The digital output products are produced in real time and stored off-line for post-mission analysis and further target enhancements. This LADAR is distinguished from other designs primarily by the waveforms produced by the laser for target interrogation. The digital processing algorithms are designed to extract certain features through operation on each of the two waveforms. The waveforms are a pulse-tone and a pulse-burst designed for target acquisition and track, and 2D imaging respectively. The algorithms are categorized by function as acquisition/track, 2D imaging, integrated absorption for gaseous targets, and post mission enhancements such as tomographic reconstruction for multiple looks at targets from different perspectives. Field tests are now in process and results acquired from Feb.-June '96 will be reported on. The digital imaging system, its architecture, algorithms, simulations, and products will be described.
Williams, David R.; Gaddis, Lisa
The tectonics of the Tellus Region highland on Venus is examined using the altimetry and gravity data collected by Pioneer Venus, which were incorporated into a thin elastic shell model to calculate both the global (long-wavelength) and the regional (short-wavelength) stresses for various assumed values of crust, lithosphere, and mantle thickness and modes of compensation. The resultant stress fields were compared to the surface morphology observed in the Venera 15/16 radar images and interpreted in terms of stress history of Tellus Regio. The best fitting parameters were found to be consistent with minor amounts of lithospheric flexure being necessary to produce the observed surface features of this region.
Zelenka, Richard E.
Avionic systems that depend on digitized terrain elevation data for guidance generation or navigational reference require accurate absolute and relative distance measurements to the terrain, especially as they approach lower altitudes. This is particularly exacting in low-altitude helicopter missions, where aggressive terrain hugging maneuvers create minimal horizontal and vertical clearances and demand precise terrain positioning. Sole reliance on airborne precision navigation and stored terrain elevation data for above-ground-level (AGL) positioning severely limits the operational altitude of such systems. A Kalman filter is presented which blends radar altimeter returns, precision navigation, and stored terrain elevation data for AGL positioning. The filter is evaluated using low-altitude helicopter flight test data acquired over moderately rugged terrain. The proposed Kalman filter is found to remove large disparities in predicted AGL altitude (i.e., from airborne navigation and terrain elevation data) in the presence of measurement anomalies and dropouts. Previous work suggested a minimum clearance altitude of 220 ft AGL for a near-terrain guidance system; integration of a radar altimeter allows for operation of that system below 50 ft, subject to obstacle-avoidance limitations.
Robinson, S.H.; Morrison, R.E.
A continuous wave doppler radar system has been designed which is portable, easily deployable and can be remotely controlled. The system is immune to ground clutter and is used for wind speed detection and direction determination. Nearly real time digital signal processing is performed by an Analog Devices ADSP-21020, a 40-bit floating point Digital Signal Processing (DSP) microprocessor. This paper provides an overview of the design of the system including the radio frequency (RF) to digital interface. The various DSP detection algorithms are discussed and compared to system performance and sensitivity. Finally, DSP performance is compared to the performance of an earlier system using Analog Device's ADSP-2100. 6 refs.
Liskovich, Diana; Simard, Marc
Using radar and lidar data, the aim is to improve 3D rendering of terrain, including digital elevation models (DEM) and estimates of vegetation height and biomass in a variety of forest types and terrains. The 3D mapping of vegetation structure and the analysis are useful to determine the role of forest in climate change (carbon cycle), in providing habitat and as a provider of socio-economic services. This in turn will lead to potential for development of more effective land-use management. The first part of the project was to characterize the Shuttle Radar Topography Mission DEM error with respect to ICESat/GLAS point estimates of elevation. We investigated potential trends with latitude, canopy height, signal to noise ratio (SNR), number of LiDAR waveform peaks, and maximum peak width. Scatter plots were produced for each variable and were fitted with 1st and 2nd degree polynomials. Higher order trends were visually inspected through filtering with a mean and median filter. We also assessed trends in the DEM error variance. Finally, a map showing how DEM error was geographically distributed globally was created.
Neish, Catherine D.; Carter, Lynn M.
This chapter describes the principles of planetary radar, and the primary scientific discoveries that have been made using this technique. The chapter starts by describing the different types of radar systems and how they are used to acquire images and accurate topography of planetary surfaces and probe their subsurface structure. It then explains how these products can be used to understand the properties of the target being investigated. Several examples of discoveries made with planetary radar are then summarized, covering solar system objects from Mercury to Saturn. Finally, opportunities for future discoveries in planetary radar are outlined and discussed.
Salcher, Bernhard; Baumann, Sebastian; Kober, Florian; Robl, Jörg; Heiniger, Lukas
The analysis of the slope-area relationship in bedrock streams is a common way for discriminating the channel from the hillslope domain and associated landscape processes. Spatial variations of these domains are important indicators of landscape change. In fluvial catchments, this relationship is a function of contributing drainage area, channel slope and the threshold drainage area for fluvial erosion. The resulting pattern is related to climate, tectonic and underlying bedrock. These factors may become secondary in catchments affected by glacial erosion, as it is the case in many mid- to high-latitude mountain belts. The perturbation (i.e. the destruction) of an initial steady state fluvial bedrock morphology (where uplift is balanced by surface lowering rates) will tend to become successively larger if the repeated action of glacial processes exceeds the potential of fluvial readjustment during deglaciated periods. Topographic change is associated with a decrease and fragmentation of the channel network and an extension of the hillslope domain. In case of glacially conditioned catchments discrimination of the two domains remains problematic and a discrimination inconsistent. A definition is therefore highly needed considering that (i) a spatial shift in the domains affect the process and rate of erosion and (ii) topographic classifications of alpine catchments often base on channel and hillslope parameters (i.e.channel or hillslope relief). Here we propose a novel numerical approach to topographically define channel heads from digital topography in glacially conditioned mountain range catchments in order to discriminate the channel from the hillslope domain. We analyzed the topography of the southern European Central Alps, a region which (i) has been glaciated multiple times during the Quaternary, shows (ii) little lithological variations, is (iii) home of very low erodible rocks and is (iv) known as a region were tectonic processes have largely ceased. The
Lahaie, I. J.; Dias, A. R.; Darling, G. D.
The digital processing requirements of several algorithms for extracting the spectrum of a detected synthetic aperture radar (SAR) image from the raw SAR data are described and compared. The most efficient algorithms for image spectrum extraction from raw SAR data appear to be those containing an intermediate image formation step. It is shown that a recently developed compact formulation of the image spectrum in terms of the raw data is computationally inefficient when evaluated directly, in comparison with the classical method where matched-filter image formation is an intermediate result. It is also shown that a proposed indirect procedure for digitally implementing the same compact formulation is somewhat more efficient than the classical matched-filtering approach. However, this indirect procedure includes the image formation process as part of the total algorithm. Indeed, the computational savings afforded by the indirect implementation are identical to those obtained in SAR image formation processing when the matched-filtering algorithm is replaced by the well-known 'dechirp-Fourier transform' technique. Furthermore, corrections to account for slant-to-ground range conversion, spherical earth, etc., are often best implemented in the image domain, making intermediate image formation a valuable processing feature.
Salamunićcar, Goran; Vinković, Dejan; Lončarić, Sven; Vučina, Damir; Pehnec, Igor; Vojković, Marin; Gomerčić, Mladen; Hercigonja, Tomislav
In our previous work the following has been done: (1) the crater detection algorithm (CDA) based on digital elevation model (DEM) has been developed and the GT-115225 catalog has been assembled [GRS, 48 (5), in press, doi:10.1109/TGRS.2009.2037750]; and (2) the results of comparison between explosion-induced laboratory craters in stone powder surfaces and GT-115225 have been presented using depth/diameter measurements [41stLPSC, Abstract #1428]. The next step achievable using the available technology is to create 3D scans of such labo-ratory craters, in order to compare different properties with simple Martian craters. In this work, we propose a formal method for evaluation of laboratory craters, in order to provide objective, measurable and reproducible estimation of the level of achieved similarity between these laboratory and real impact craters. In the first step, the section of MOLA data for Mars (or SELENE LALT for Moon) is replaced with one or several 3D-scans of laboratory craters. Once embedment was done, the CDA can be used to find out whether this laboratory crater is similar enough to real craters, as to be recognized as a crater by the CDA. The CDA evaluation using ROC' curve represents how true detection rate (TDR=TP/(TP+FN)=TP/GT) depends on the false detection rate (FDR=FP/(TP+FP)). Using this curve, it is now possible to define the measure of similarity between laboratory and real impact craters, as TDR or FDR value, or as a distance from the bottom-right origin of the ROC' curve. With such an approach, the reproducible (formally described) method for evaluation of laboratory craters is provided.
Mukul, Manas; Srivastava, Vinee; Mukul, Malay
The Shuttle Radar Topography Mission (SRTM) carried out in February 2000 has provided near global topographic data that has been widely used in many fields of earth sciences. The mission goal of an absolute vertical accuracy within 16 m (with 90% confidence)/RMSE ˜10 m was achieved based on ground validation of SRTM data through various studies using global positioning system (GPS). We present a new and independent assessment of the vertical accuracy of both the X- and C-band SRTM datasets using data from the International GNSS Service (IGS) network of high-precision static GPS stations. These stations exist worldwide, have better spatial distribution than previous studies, have a vertical accuracy of 6 mm and constitute the most accurate ground control points (GCPs) possible on earth; these stations are used as fiducial stations to define the International Terrestrial Reference Frame (ITRF). Globally, for outlier-filtered data (135 X-band stations and 290 C-band stations), the error or difference between IGS and SRTM heights exhibits a non-normal distribution with a mean and standard error of 8.2 ± 0.7 and 6.9 ± 0.5 m for X- and C-band data, respectively. Continent-wise, Africa, Australia and North America comply with the SRTM mission absolute vertical accuracy of 16 m (with 90% confidence)/RMSE ˜10 m. However, Asia, Europe and South America have vertical errors higher than the SRTM mission goal. At stations where both the X- and C-band SRTM data were present, the root mean square error (RMSE) of both the X- and C-bands was identical at 11.5 m, indicating similar quality of both the X- and C-band SRTM data.
Morrison, R.E.; Robinson, S.H.
A continuous wave Doppler radar system has been designed which is portable, easily deployed, and remotely controlled. The heart of this system is a DSP/control board using Analog Devices ADSP-21020 40-bit floating point digital signal processor (DSP) microprocessor. Two 18-bit audio A/D converters provide digital input to the DSP/controller board for near real time target detection. Program memory for the DSP is dual ported with an Intel 87C51 microcontroller allowing DSP code to be up-loaded or down-loaded from a central controlling computer. The 87C51 provides overall system control for the remote radar and includes a time-of-day/day-of-year real time clock, system identification (ID) switches, and input/output (I/O) expansion by an Intel 82C55 I/O expander. 5 refs., 8 figs., 2 tabs.
Hoffman, James Patrick; Peral, Eva; Veilluex, Louise; Perkovic, Dragana; Shaffer, Scott
Real-time digital beamforming, combined with lightweight, large aperture reflectors, enable SweepSAR architectures such as that of the proposed DESDynI [Deformation, Ecosystem Structure, and Dynamics of Ice] SAR [Synthetic Aperture Radar] Instrument (or DSI). SweepSAR promises significant increases in instrument capability for solid earth and biomass remote sensing, while reducing mission mass and cost. This new instrument concept requires new methods for calibrating the multiple channels, which must be combined on-board, in real-time. We are developing new methods for digitally calibrating digital beamforming arrays to reduce development time, risk and cost of precision calibrated TR modules for array architectures by accurately tracking modules' characteristics through closed-loop Digital Calibration, thus tracking systematic changes regardless of temperature
Gruber, Fabian E.; Baruck, Jasmin; Hastik, Richard; Geitner, Clemens
All major soil description and classification systems, including the World Reference Base (WRB) and the German Soil description guidelines (KA5), require the characterization of landform and topography for soil profile sites. This is commonly done at more than one scale, for instance at macro-, meso- and micro scale. However, inherent when humans perform such a task, different surveyors will reach different conclusions due to their subjective perception of landscape structure, based on their individual mind-model of soil-landscape structure, emphasizing different aspects and scales of the landscape. In this study we apply a work-flow using the GRASS GIS extension module r.geomorphon to make use of high resolution digital elevation models (DEMs) to characterize the landform elements and topography of soil profile sites at different scales, and compare the results with a large number of soil profile site descriptions performed during the course of forestry surveys in South and North Tyrol (Italy and Austria, respectively). The r.geomorphon extension module for the open source geographic information system GRASS GIS applies a pattern recognition algorithm to delineate landform elements based on an input DEM. For each raster cell it computes and characterizes the visible neighborhood using line-of-sight calculations and then applies a lookup-table to classify the raster cell into one of ten landform elements (flat, peak, ridge, shoulder, slope, spur, hollow, footslope, valley and pit). The input parameter search radius (L) represents the maximum number of pixels for line-of-sight calculation, resulting in landforms larger than L to be split into landform components. The use of these visibility calculations makes this landform delineation approach suitable for comparison with the landform descriptions of soil surveyors, as their spatial perception of the landscape surrounding a soil profile site certainly influences their classification of the landform on which the
Wells, R. A.
Using direct radar ranging of surface heights on Mars and spectrophotometric observations of absorptions produced by carbon dioxide molecules in the Martian atmosphere, data have been obtained on Martian topographical variations at spatial resolutions ranging from about 100 to 1000 km. These data have been studied and analyzed. As a result, a surface height contour map has been produced which clearly reveals a structural complex of blocks and basins whose distribution enhances the magnitude of low-degree surface harmonics. It is emphasized that Mars possesses unexpectedly pronounced topography which can have important geophysical consequences.
Nakamura, Yoshitaka; Yoshikawa, Eiichi; Akita, Manabu; Morimoto, Takeshi; Ushio, Tomoo; Kawasaki, Zen-Ichiro; Saito, Toshiya; Nishida, Takashi; Sakazume, Norio
We propose a high-resolution precipitation and lightning monitoring for meteorological application. This monitoring is mainly utilized the Ku-band broadband radar (BBR) and the VHF broadband digital interferometer (DITF). The BBR can accurately measure the radar reflectivity factor and the mean Doppler velocity with 5 m resolution over a range from 40 m to several kilometers for 10 W power using a pulse compression technique. The two or more DITFs make us visualize lightning channel propagations in three dimensions. Moreover, we add new functions that integrate these observation data and disclose integration analyses results with the quasi real-time information disclosure system. Initial observations for severe storms with lightning during summer and winter thunderstorm season by these monitoring instruments indicate that we obtain detailed precipitation distribution and detect active convective cells with lightning discharges.
Blair, J. Bryan; Rabine, David L.; Hofton, Michelle A.
The Laser Vegetation Imaging Sensor (LVIS) is an airborne, scanning laser altimeter designed and developed at NASA's Goddard Space Flight Center. LVIS operates at altitudes up to 10 km above ground, and is capable of producing a data swath up to 1000 m wide nominally with 25 m wide footprints. The entire time history of the outgoing and return pulses is digitized, allowing unambiguous determination of range and return pulse structure. Combined with aircraft position and attitude knowledge, this instrument produces topographic maps with decimeter accuracy and vertical height and structure measurements of vegetation. The laser transmitter is a diode-pumped Nd:YAG oscillator producing 1064 nm, 10 nsec, 5 mJ pulses at repetition rates up to 500 Hz. LVIS has recently demonstrated its ability to determine topography (including sub-canopy) and vegetation height and structure on flight missions to various forested regions in the U.S. and Central America. The LVIS system is the airborne simulator for the Vegetation Canopy Lidar (VCL) mission (a NASA Earth remote sensing satellite due for launch in 2000), providing simulated data sets and a platform for instrument proof-of-concept studies. The topography maps and return waveforms produced by LVIS provide Earth scientists with a unique data set allowing studies of topography, hydrology, and vegetation with unmatched accuracy and coverage.
Harding, D. J.; Carabajal, C. C.
Interferometric Synthetic Aperture Radar (InSAR) and laser altimeter measurements of topography provide complimentary approaches to characterize landforms. Results from the Shuttle Radar Topography Mission (SRTM) will provide an unprecedented, near-global, public-domain topography data set at 90 m resolution using a single pass C-band (5.6 cm wavelength) radar interferometer. In vegetated terrains, the C-band radar energy will penetrate part way into vegetation cover. The elevation of the resulting radar phase center, somewhere between the canopy top and underlying ground, will depend on the vegetation density, structure, and presence or absence of foliage. The high vertical accuracy and spatial resolution achieved by laser altimeters, and their capability to directly measure vegetation height and ground topography beneath vegetation cover, provides a method to evaluate InSAR representations of topography. Here a preliminary C-band SRTM digital elevation model (DEM) for a portion of western Washington State is evaluated using laser altimeter data to assess its elevation accuracy and the extent of vegetation penetration. The SRTM DEM extends from the Cascades Range westward to the Olympic Peninsula. The laser altimeter data includes two profiles acquired by the second flight of the Shuttle Laser Altimeter (SLA-02) in August, 1997, numerous transects acquired by the airborne Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) in September, 1995, and comprehensive mapping in the Puget Lowland region acquired by Terrapoint, LLC for the Puget Sound Lidar Consortium in the winters of 2000 and 2001. SLA-02 and SLICER acquired waveforms that record the height distribution of illuminated surfaces within 120 m and 10 m diameter footprints, respectively. The Terrapoint elevations consist of up to four discrete returns from 1 m footprints spaced 1.5 apart, with all areas mapped twice. Methods for comparing laser altimeter and SRTM topography developed here will be
Held, D.N.; Gasparovic, R.F.; Mansfield, A.W.; Melville, W.K.; Mollo-Christensen, E.L.; Zebker, H.A.
Radar non-acoustic anti-submarine warfare (NAASW) became the subject of considerable scientific investigation and controversy in the West subsequent to the discovery by the Seasat satellite in 1978 that manifestations of underwater topography, thought to be hidden from the radar, were visible in synthetic aperture radar (SAR) images of the ocean. In addition, the Seasat radar produced images of ship wakes where the observed angle between the wake arms was much smaller than expected from classical Kelvin wake theory. These observations cast doubt on the radar oceanography community's ability to adequately explain these phenomena, and by extension on the ability of existing hydrodynamic and radar scattering models to accurately predict the observability of submarine-induced signatures. If one is of the opinion that radar NAASW is indeed a potentially significant tool in detecting submerged operational submarines, then the Soviet capability, as evidenced throughout this report, will be somewhat daunting. It will be shown that the Soviets have extremely fine capabilities in both theoretical and experimental hydrodynamics, that Soviet researchers have been conducting at-sea radar remote sensing experiments on a scale comparable to those of the United States for several years longer than we have, and that they have both an airborne and spaceborne SAR capability. The only discipline that the Soviet Union appears to be lacking is in the area of digital radar signal processing. If one is of the opinion that radar NAASW can have at most a minimal impact on the detection of submerged submarines, then the Soviet effort is of little consequence and poses not threat. 280 refs., 31 figs., 12 tabs.
Yasodha, Polisetti; Jayaraman, Achuthan; Kamaraj, Pandian; Durga rao, Meka; Thriveni, A.
Modern phased array radars depend highly on digital signal processing (DSP) to extract the echo signal information and to accomplish reliability along with programmability and flexibility. The advent of ASIC technology has made various digital signal processing steps to be realized in one DSP chip, which can be programmed as per the application and can handle high data rates, to be used in the radar receiver to process the received signal. Further, recent days field programmable gate array (FPGA) chips, which can be re-programmed, also present an opportunity to utilize them to process the radar signal. A multi-channel direct IF/RF digital receiver (MCDRx) is developed at NARL, taking the advantage of high speed ADCs and high performance DSP chips/FPGAs, to be used for atmospheric radars working in HF/VHF bands. Multiple channels facilitate the radar t be operated in multi-receiver modes and also to obtain the wind vector with improved time resolution, without switching the antenna beam. MCDRx has six channels, implemented on a custom built digital board, which is realized using six numbers of ADCs for simultaneous processing of the six input signals, Xilinx vertex5 FPGA and Spartan6 FPGA, and two ADSPTS201 DSP chips, each of which performs one phase of processing. MCDRx unit interfaces with the data storage/display computer via two gigabit ethernet (GbE) links. One of the six channels is used for Doppler beam swinging (DBS) mode and the other five channels are used for multi-receiver mode operations, dedicatedly. Each channel has (i) ADC block, to digitize RF/IF signal, (ii) DDC block for digital down conversion of the digitized signal, (iii) decoding block to decode the phase coded signal, and (iv) coherent integration block for integrating the data preserving phase intact. ADC block consists of Analog devices make AD9467 16-bit ADCs, to digitize the input signal at 80 MSPS. The output of ADC is centered around (80 MHz - input frequency). The digitized data is fed
This document describes the design and operation of a low resolution radar data recording system for precipitation measurements. This system records a full azimuth scan on seven track magnetic tapes every five minutes. It is designed to operate on a continuous basis with operator intervention required only for changing tape reels and calibration.
Young, Steven D.; deHaag, Maatren Uijt
Flight in Instrument Meteorological Conditions requires pilots to manipulate flight controls while referring to a Primary Flight Display. The Primary Flight Display indicates aircraft attitude along with, in some cases, many other state variables such as altitude, speed, and guidance cues. Synthetic Vision Systems have been proposed that overlay the traditional information provided on Primary Flight Displays onto a scene depicting the location of terrain and other geo-spatial features.Terrain models used by these displays must have sufficient quality to avoid providing misleading information. This paper describes how X-band radar measurements can be used as part of a monitor, and/or maintenance system, to quantify the integrity of terrain models that are used by systems such as Synthetic Vision. Terrain shadowing effects, as seen by the radar, are compared in a statistical manner against estimated shadow feature elements extracted from the stored terrain model from the perspective of the airborne observer. A test statistic is defined that enables detection of errors as small as the range resolution of the radar. Experimental results obtained from two aircraft platforms hosting certified commercial-off-the-shelf X-band radars test the premise and illustrate its potential.
Danielson, E. F.; Hipskind, R. S.; Gaines, S. E.
Results are presented from computer processing and digital filtering of radiosonde and radar tracking data obtained during the ITCZ experiment when coordinated measurements were taken daily over a 16 day period across the Panama Canal Zone. The temperature relative humidity and wind velocity profiles are discussed.
Young, Steve; UijtdeHaag, Maarten; Sayre, Jonathon
Synthetic Vision Systems (SVS) provide pilots with displays of stored geo-spatial data representing terrain, obstacles, and cultural features. As comprehensive validation is impractical, these databases typically have no quantifiable level of integrity. Further, updates to the databases may not be provided as changes occur. These issues limit the certification level and constrain the operational context of SVS for civil aviation. Previous work demonstrated the feasibility of using a realtime monitor to bound the integrity of Digital Elevation Models (DEMs) by using radar altimeter measurements during flight. This paper describes an extension of this concept to include X-band Weather Radar (WxR) measurements. This enables the monitor to detect additional classes of DEM errors and to reduce the exposure time associated with integrity threats. Feature extraction techniques are used along with a statistical assessment of similarity measures between the sensed and stored features that are detected. Recent flight-testing in the area around the Juneau, Alaska Airport (JNU) has resulted in a comprehensive set of sensor data that is being used to assess the feasibility of the proposed monitor technology. Initial results of this assessment are presented.
Wilson, G. S.
An analysis is conducted of a small 4-year climatological data base of Manually Digitized Radar (MDR) data to infer the monthly and seasonal distributions of the relative frequency of occurrence of potential flash flood-producing rains over the Central and Eastern U.S. Some possible meteorological mechanisms for producing potential flash flooding rains are discussed in terms of the relative maxima and minima in the monthly and seasonal frequency distributions over the MDR network. Frequencies were found to be generally higher in more southern locations and lower farther north in all months and seasons. However, most locations experienced an annual cycle in the frequency of occurrence with maxima in summer and minima in winter. In given seasons and months, local areas of maximum and minimum occurrences may be related to quasi-stationary meteorological processes that trigger and organize intense convection over a common area.
An area near Long Valley, California, was mapped by the Spaceborne Imaging Radar-C and X-band Synthetic Aperture Radar aboard the space shuttle Endeavor on April 13, 1994, during the first flight of the radar instrument, and on October 4, 1994, during the second flight of the radar instrument. The orbital configurations of the two data sets were ideal for interferometric combination -- that is overlaying the data from one image onto a second image of the same area to create an elevation map and obtain estimates of topography. Once the topography is known, any radar-induced distortions can be removed and the radar data can be geometrically projected directly onto a standard map grid for use in a geographical information system. The 50 kilometer by 50 kilometer (31 miles by 31 miles) map shown here is entirely derived from SIR-C L-band radar (horizontally transmitted and received) results. The color shown in this image is produced from the interferometrically determined elevations, while the brightness is determined by the radar backscatter. The map is in Universal Transverse Mercator (UTM) coordinates. Elevation contour lines are shown every 50 meters (164 feet). Crowley Lake is the dark feature near the south edge of the map. The Adobe Valley in the north and the Long Valley in the south are separated by the Glass Mountain Ridge, which runs through the center of the image. The height accuracy of the interferometrically derived digital elevation model is estimated to be 20 meters (66 feet) in this image. Spaceborne Imaging Radar-C and X-band Synthetic Aperture Radar (SIR-C/X-SAR) is part of NASA's Mission to Planet Earth. The radars illuminate Earth with microwaves, allowing detailed observations at any time, regardless of weather or sunlight conditions. SIR-C/X-SAR uses three microwave wavelengths: L-band (24 cm), C-band (6 cm) and X-band (3 cm). The multi-frequency data will be used by the international scientific community to better understand the global
Li, Weixing; Zhang, Yue; Lin, Jianzhi; Chen, Zengping
Amplitude-phase errors and mutual coupling errors among multi-channels in digital array radar (DAR) will seriously deteriorate the performance of signal processing such as digital beam-forming (DBF) and high resolution direction finding. In this paper, a combined algorithm for error calibration in DAR has been demonstrated. The algorithm firstly estimates the amplitude-phase errors of each channel using interior calibration sources with the help of the calibration network. Then the signals from far field are received and the amplitude-phase errors are compensated. According to the subspace theories, the relationship between the principle eigenvectors and distorted steering vectors is expressed, and the cost function containing the mutual coupling matrix (MCM) and incident directions is established. Making use of the properties of MCM of uniform linear array, Gauss-Newton method is implied to iteratively compute the MCM and the direction of arrival (DOA). Simulation results have shown the effectiveness and performance of proposed algorithm. Based on an 8-elements DAR test-bed, experiments are carried out in anechoic chamber. The results illustrate that the algorithm is feasible in actual systems.
Adams, John W.; Hudson, Ralph E.; Bayma, Robert W.; Nelson, Jeffrey E.
The polar format approach to SAR image formation requires data to be interpolated from a warped grid onto a Cartesian lattice. In general, this requires that data be interpolated between varying sampling rates. In this paper, frequency-domain optimality criteria for polar format interpolators are defined and justified, and an approach to designing the corresponding digital filters is described.
Salamunićcar, Goran; Lončarić, Sven
Crater detection algorithms (CDAs) are an important subject of recent scientific research, as evident from the numerous recent publications in the field [ASR, 42 (1), 6-19]. In our previous work: (1) all the craters from the major currently available manually assembled catalogues have been merged into the catalogue with 57633 known Martian impact-craters [PSS, 56 (15), 1992-2008]; and (2) the CDA (developed to search for still uncatalogued impact-craters using 1/128° MOLA data) has been used to extend GT-57633 catalogue with 57592 additional craters resulting in GT-115225 catalog [GRS, 48 (5), in press, doi:10.1109/TGRS.2009.2037750]. On the other hand, the most complete catalog for Moon is the Morphological catalog of Lunar craters [edited by V. V. Shevchenko], which includes information on 14923 craters larger than 10km, visible on the lunar nearside and farside. This was the main motivation for application of our CDA to newly available Lunar SELENE LALT data. However, one of the main differences between MOLA and LALT data is the highest available resolution, wherein MOLA is available in 1/128° and LALT in 1/16° . The consequence is that only the largest craters can be detected using LALT dataset. However, this is still an excellent opportunity for further work on CDA in order to prepare it for forthcoming LRO LOLA data (which is expected to be in even better resolution than MOLA). The importance is in the fact that morphologically Martian and Lunar craters are not the same. Therefore, it is important to use the dataset for Moon in order to work on the CDA which is meant for detection of Lunar craters as well. In order to overcome the problem of currently available topography data in low resolution only, we particularly concentrated our work on the CDA's capability to detect very small craters relative to available dataset (up to the extreme case wherein the radius is as small as only two pixels). For this purpose, we improved the previous CDA with a new
Rodriguez, Ernesto; Esteban-Fernandez, Daniel
The Surface Water and Ocean Topography (SWOT) mission will study ocean mesoscale and submesoscale phenomena and provide an inventory of storage change and discharge for fresh water bodies and rivers. In this paper, we examine the combination of measurements that will be used by SWOT to achieve a globally consistent data set. We introduce a new channel in the SWOT measurement that combines data transmitted by the interferometer antennas and received by the radiometer antenna allows the closing of the SWOT nadir coverage gap. This new mode also allows for improved calibration between the nadir altimeter and the interferometer, resulting in consistent range measurements. Consistency in the phase measurements is achieved using a mixture of cross-over calibration combined with optimal estimation of system error drift.
Gubbels, T.L.; Isacks, B.L. ); Ellis, J.M. )
The central Andean plateau is one of the Earth's most remote and poorly mapped regions. The plateau has an average elevation of 3.7 km, and extends from central Peru to at least 30[degrees]S latitude. The plateau and flanking Subandean foldthrust belt (FTB) reach their greatest width near 20[degrees]S, and at this latitude both the FTB and the basin within the plateau (Altiplano basin) are areas of active hydrocarbon exploration. We have used Landsat TM imagery, stereoscopic SPOT imagery, and digital topography to construct a crustal-scale transect across the central Andes in order to better understand Andean tectonics at this latitude. Beginning at the Peru-Chile trench and continuing to the east, the transect crosses the Coastal Cordillera, Longitudinal Valley, Active Magmatic Arc, Altiplano basin, Eastern Cordillera, Subandean fold-thrust belt, and Subandean foreland basin. A digital elevation model across the entire region illustrates that the magmatic arc, Altiplano basin, and Eastern cordillera all lie within the plateau region. Satellite imagery across the transect illustrates the characteristic geology, structure, and geomorphology of each of the major morphotectonic regions, as well as the nature of their boundaries. The transect has led us to a number of new insights on Andean tectonics at this latitude. Most importantly, it supports a two-stage model of Andean Cenozoic growth in which a widespread Oligocene to mid-Miocene compressional deformation in the Altiplano and Eastern Cordillera is followed in the late Miocene and Pliocene by thrusting localized east of the Eastern Cordillera, forming the Subandean fold-thrust belt.
Ford, J. P.
Digitally correlated Seasat synthetic-aperture radar (SAR) images of the Alaska Range, Alaska, and the drumlin-drift belt in Ireland are analyzed for the perception and identification of geomorphic features. The two terrains display strongly contrasted types of glacial topography whose identification in each case is related to the geometry of the Seasat imaging radar. Identification of terrain shape and form is important within the caveats imposed by the intrinsic distortions on the radar images. Image texture serves coarsely to distinguish topography. Image tones are scene-dependent and do not uniquely identify specific targets. Extensive alignments of linear and curvilinear features provide some of the more important image information from which to make geologic interpretations in each case.
Johansson, G.; Borg, J.; Johansson, J.; Lundberg Nordenvaad, M.; Wannberg, G.
This paper presents simulations and methods developed to investigate the feasibility of using a Fractional-Sample-Delay (FSD) system in the planned EISCAT_3D incoherent scatter radar. Key requirements include a frequency-independent beam direction over a 30 MHz band centered around 220 MHz, with correct reconstruction of pulse lengths down to 200 ns. The clock jitter from sample to sample must be extremely low for the integer sample delays. The FSD must also be able to delay the 30 MHz wide signal band by 1/1024th of a sample without introducing phase shifts, and it must operate entirely in baseband. An extensive simulation system based on mathematical models has been developed, with inclusion of performance-degrading aspects such as noise, timing error, and bandwidth. Finite Impulse Response (FIR) filters in the baseband of a band-pass-sampled signal have been used to apply true time delay beamforming. It has been confirmed that such use is both possible and well behaved. The target beam-pointing accuracy of 0.06° is achievable using optimized FIR filters with lengths of 36 taps and an 18 bit coefficient resolution. Even though the minimum fractional delay step necessary for beamforming is ˜13.1 ps, the maximum sampling timing error allowed in the array is found to be σ ≤ 120 ps if the errors are close to statistically independent.
LeDantec, Fernando A.
This thesis is concerned with the development of a model to analyze a Digital Image Synthesizer (DIS) integrated circuit designed to create false target images to deceive Inverse Synthetic Aperture Radar (ISAR). The DIS is able to recreate the scattering effect of a moving target by using appropriate phase and gain modulations on an intercepted ISAR chirp signal before retransmitting it with the proper time delay. The DIS signal processing and the ISAR compression of the modulated return are modeled to examine the range-Doppler profile of a synthesized false target image. The image is representative of the image that would appear on an ISAR display. ISAR image quality is used to evaluate different DIS architectures and bit formats. Evaluation of the image quality is based on the deviation from an infinite resolution false target image. The results obtained from evaluating different DIS architectures indicate that the design is tolerant of significant quantization errors. The model is used to validate the architecture of the integrated circuit being fabricated. Finally, various different ISAR integration times and pulse repetition frequencies are used to confirm the integrity of the model.
This is a three-dimensional perspective view of Saline Valley, about 30 km (19 miles) east of the town of Independence, California created by combining two spaceborne radar images using a technique known as interferometry. Visualizations like this one are helpful to scientists because they clarify the relationships of the different types of surfaces detected by the radar and the shapes of the topographic features such as mountains and valleys. The view is looking southwest across Saline Valley. The high peaks in the background are the Inyo Mountains, which rise more than 3,000 meters (10,000 feet) above the valley floor. The dark blue patch near the center of the image is an area of sand dunes. The brighter patches to the left of the dunes are the dry, salty lake beds of Saline Valley. The brown and orange areas are deposits of boulders, gravel and sand known as alluvial fans. The image was constructed by overlaying a color composite radar image on top of a digital elevation map. The radar image was taken by the Spaceborne Imaging Radar-C/X-bandSynthetic Aperture Radar (SIR-C/X-SAR) on board the space shuttleEndeavour in October 1994. The digital elevation map was producedusing radar interferometry, a process in which radar data are acquired on different passes of the space shuttle. The two data passes are compared to obtain elevation information. The elevation data were derived from a 1,500-km-long (930-mile) digital topographic map processed at JPL. Radar image data are draped over the topography to provide the color with the following assignments: red is L-band vertically transmitted, vertically received; green is C-band vertically transmitted, vetically received; and blue is the ratio of C-band vertically transmitted, vertically received to L-band vertically transmitted, vertically received. This image is centered near 36.8 degrees north latitude and 117.7 degrees west longitude. No vertical exaggeration factor has been applied to the data. SIR-C/X-SAR, a joint
Child, S.; Stearns, L. A.; Purdon, K.; Li, J.; Rodriguez-Morales, F.; Crowe, R.; Gomez-Garcia, D.
Bed topography controls the pattern and magnitude of ice velocity, far into the catchment basin of many Antarctic outlet glaciers. Predictive models of glacier dynamics and ice sheet mass balance rely on well-prescribed bed topography, but in many regions bed topography is largely unknown. This particular study investigates the relationship between bed and surface topography and ice velocity in the catchment basin and trunk of Byrd Glacier. Byrd Glacier drains ~19 % of the area of the East Antarctic Ice Sheet (1,070,400 km2), and has the potential to play a significant role in the ice sheet's total mass balance. In 2011/2012, the Center for Remote Sensing of Ice Sheets (CReSIS) collected airborne radar data over Byrd Glacier. These new measurements of bed topography, along with updated surface digital elevation models (DEMs) and basin-wide ice velocity maps, are used to investigate the flow dynamics of Byrd Glacier with improved accuracy. Surface DEMs are derived from new high-resolution WorldView imagery; ice velocity is derived from repeat visible imagery, coupled with InSAR results (Rignot et al., 2011). Results exhibit relatively smooth depressions surrounding the inferred subglacial lakes, ~200 km upflow of the grounding line on Byrd Glacier (Stearns et al., 2008). Downflow of the subglacial lakes is a complex pattern of hills and valleys as ice enters the glacier trunk. At the mouth of the trunk is a large overdeepening (~2500 m) that coincides with faster ice flow. We use along- and across-flow radar profiles to perform detailed comparisons of ice velocity, bed topography and surface topography throughout the Byrd Glacier region. Gridded products are used to complete an updated force balance assessment. These results provide us with a better understanding of Byrd Glacier's flow dynamics and sensitivity to external perturbations.
The coherent radar technique is reviewed with special emphasis to mesosphere-stratosphere-troposphere (MST) radars operating in the VHF band. Some basic introduction to Doppler radar measurements and the radar equation is followed by an outline of the characteristics of atmospheric turbulence, viewed from the scattering and reflection processes of radar signals. Radar signal acquisition and preprocessing, namely coherent detection, digital sampling, pre-integration and coding, is briefly discussed. The data analysis is represented in terms of the correlation and spectrum analysis, yielding the essential parameters: power, signal-to-noise ratio, average and fluctuating velocity and persistency. The techniques to measure wind velocities, viz. the different modes of the Doppler method as well as the space antenna method are surveyed and the feasibilities of the MST radar interferometer technique are elucidated. A general view on the criteria to design phased array antennas is given. An outline of the hardware of a typical MST radar system is presented.
Widyatmanti, Wirastuti; Wicaksono, Ikhsan; Dinta Rahma Syam, Prima
Dense vegetation that covers most landscapes in Indonesia becomes a common limitation in mapping the landforms in tropical region. This paper aims to examine the use of radar interferometry for landform mapping in tropical region; to examine the application of segmentation method to develop landform type boundaries; and to identify the topographic elements composition for each type of landform. Using Idrisi® and “eCognition ®” softwares, toposhape analysis, segmentation and multi-spectral classification were applied to identify the composition of topographic elements i.e. the types of land-cover from Landsat 8, elevation, slope, relief intensity and curvatures from SRTM (DEM). Visual interpretation on DEM and land-cover fusion imagery was conducted to derive basic control maps of landform and land-cover. The result shows that in segmentation method, shape and compactness levels are essential in obtaining land-cover, elevation, and slope class units to determine the most accurate class borders of each element. Despite a complex procedure applied in determining landform classification, the combination of topographic elements segmentation result presents a distinct border of each landform class. The comparison between landform maps derived from segmentation process and visual interpretation method demonstrates slight dissimilarities, meaning that multi-stage segmentation approach can improve and provide more effective digital landform mapping method in tropical region. Topographic elements on each type of landforms show distinctive composition key containing the percentage of each curvature elements per area unit. Supported by GIS programming and modeling in the future, this finding is significant in reducing effort in landform mapping using visual interpretation method for a very large coverage but in detail scale level.
Mukul, Manas; Srivastava, Vinee; Mukul, Malay
Global Shuttle Radar Topography Mission (SRTM) data products have been widely used in Earth Sciences without an estimation of their accuracy and reliability even though large outliers exist in them. The global 1 arc-sec, 30 m resolution, SRTM C-Band (C-30) data collected in February 2000 has been recently released (2014-2015) outside North America. We present the first global assessment of the vertical accuracy of C-30 data using Ground Control Points (GCPs) from the International GNSS Service (IGS) Network of high-precision static fiducial stations that define the International Terrestrial Reference Frame (ITRF). Large outliers (height error ranging from -1285 to 2306 m) were present in the C-30 dataset and 14% of the data were removed to reduce the root mean square error (RMSE) of the dataset from ˜187 to 10.3 m which is close to the SRTM goal of an absolute vertical accuracy of RMSE ˜10 m. Globally, for outlier-filtered data from 287 GCPs, the error or difference between IGS and SRTM heights exhibited a non-normal distribution with a mean and standard error of 6.5 ± 0.5 m. Continent-wise, only Australia, North and South America complied with the SRTM goal. At stations where all the X- and C-Band SRTM data were present, the RMSE of the outlier-filtered C-30 data was 11.7 m. However, the RMSE of outlier-included dataset where C- and X-Band data were present was ˜233 m. The results suggest that the SRTM data must only be used after regional accuracy analysis and removal of outliers. If used raw, they may produce results that are statistically insignificant with RMSE in 100s of meters.
Liu, C.; Wen, H.; Liu, J.; Ko, M.; Yan, H.; Chang, L.
translative offset between two similar images to be rapidly estimated. To meet the requirements in remote sensing and biomedical imaging, the technology of phase correlation has been extended to the sub-pixel level. Liu and Yan (2008) developed a robust phase correlation model using the based feature matching for image co-registration and DEM generation. Considering the fact that the Formosat-2 consecutive images are intrinsically stereo pairs with very narrow baselines, this innovative stereo-matching algorithm based on SPPC technique is employed to process Formosat-2 daily revisit stereo pairs with very narrow baselines. The detailed accuracy and efficiency analysis is investigated for the study area, Namasha, Kaohsiung, using the 50cm resolution aerial photo and the 2m resolution DEM derived from airborne LiDAR data. The archive of Formosat-2 images in Taiwan area collected from 2005 to 2012 was screened out, with the intention to select the consecutive pairs of those areas where major slope disasters occurred in the past eight years. This research encourages the repeated topography surveys of geomorphic changes using digital surface models deriving from Formosat-2 daily revisit stereo pair with very narrow baseline.
Woodman, R. F.
Decoding techniques and equipment used by MST radars are described and some recommendations for new systems are presented. Decoding can be done either by software in special-purpose (array processors, etc.) or general-purpose computers or in specially designed digital decoders. Both software and hardware decoders are discussed and the special case of decoding for bistatic radars is examined.
Rosen, P. A.; Hensley, S.; Joughin, I. R.; Li, F.; Madsen, S. N.; Rodriguez, E.; Goldstein, R. M.
Synthetic aperture radar interferometry is an imaging technique for measuring the topography of a surface, its changes over time, and other changes in the detailed characteristics of the surface. This paper reviews the techniques of interferometry, systems and limitations, and applications in a rapidly growing area of science and engineering.
This image of Fort Irwin in California's Mojave Desert compares interferometric radar signatures topography -- data that were obtained by multiple imaging of the same region to produce three-dimensional elevation maps -- as it was obtained on October 7-8, 1994 by the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar aboard the space shuttle Endeavour. Data were acquired using the L-band (24 centimeter wavelength) and C-band (6 centimeter wavelength). The image covers an area about 25 kilometers by 70 kilometers (15.5 miles by 43 miles). North is to the lower right of the image. The color contours shown are proportional to the topographic elevation. With a wavelength one-fourth that of the L-band, the results from the C-band cycle through the color contours four times faster for a given elevation change. Detailed comparisons of these multiple frequency data over different terrain types will provide insights in the future into wavelength-dependent effects of penetration and scattering on the topography measurement accuracy. Fort Irwin is an ideal site for such detailed digital elevation model comparisons because a number of high precision digital models of the area already exist from conventional measurements as well as from airborne interferometric SAR data. Spaceborne Imaging Radar-C and X-band Synthetic Aperture Radar (SIR-C/X-SAR) is part of NASA's Mission to Planet Earth. The radars illuminate Earth with microwaves, allowing detailed observations at any time, regardless of weather or sunlight conditions. SIR-C/X-SAR uses three microwave wavelengths: L-band (24 cm), C-band (6 cm) and X-band (3 cm). The multi-frequency data will be used by the international scientific community to better understand the global environment and how it is changing. The SIR-C/X-SAR data, complemented by aircraft and ground studies, will give scientists clearer insights into those environmental changes which are caused by nature and those changes which are induced by human
Ostro, S. J.
The present investigation is concerned with planetary radar research reported during the time from 1979 to 1982. A brief synopsis of radar definitions and technical terminology is also provided. In connection with the proximity of the moon to earth, lunar radar studies have been performed over a wider range of wavelengths than radar investigations of other planetary targets. The most recent study of lunar quasispecular scattering is due to Simpson and Tyler (1982). The latest efforts to interpret the lunar radar maps focus on maria-highlands regolith differences and models of crater ejecta evolution. The highly successful Pioneer Venus Radar Mapper experiment has provided a first look at Venus' global distributions of topography, lambda 17-cm radar reflectivity, and rms surface slopes. Attention is given to recent comparisons of Viking Orbiter images of Mars to groundbased radar altimetry of the planet, the icy Galilean satellites, radar observations of asteroids and comets, and lambda 4-cm and lambda 13-cm observations of Saturn's rings.
Raymond, C. A.; Jaumann, R.; Nathues, A.; Sierks, H.; Roatsch, T.; Preusker, E; Scholten, F.; Gaskell, R. W.; Jorda, L.; Keller, H.-U.; Zuber, M. T.; Smith, D. E.; Mastrodemos, N.; Mottola, S.
The objective of the Dawn topography investigation is to derive the detailed shapes of 4 Vesta and 1 Ceres in order to create orthorectified image mosaics for geologic interpretation, as well as to study the asteroids' landforms, interior structure, and the processes that have modified their surfaces over geologic time. In this paper we describe our approaches for producing shape models, plans for acquiring the needed image data for Vesta, and the results of a numerical simulation of the Vesta mapping campaign that quantify the expected accuracy of our results. Multi-angle images obtained by Dawn's framing camera will be used to create topographic models with 100 m/pixel horizontal resolution and 10 m height accuracy at Vesta, and 200 m/pixel horizontal resolution and 20 m height accuracy at Ceres. Two different techniques, stereophotogrammetry and stereophotoclinometry, are employed to model the shape; these models will be merged with the asteroidal gravity fields obtained by Dawn to produce geodetically controlled topographic models for each body. The resulting digital topography models, together with the gravity data, will reveal the tectonic, volcanic and impact history of Vesta, and enable co-registration of data sets to determine Vesta's geologic history. At Ceres, the topography will likely reveal much about processes of surface modification as well as the internal structure and evolution of this dwarf planet.
vanZyl, Jakob J.
The field of synthetic aperture radar changed dramatically over the past decade with the operational introduction of advance radar techniques such as polarimetry and interferometry. Radar polarimetry became an operational research tool with the introduction of the NASA/JPL AIRSAR system in the early 1980's, and reached a climax with the two SIR-C/X-SAR flights on board the space shuttle Endeavour in April and October 1994. Radar interferometry received a tremendous boost when the airborne TOPSAR system was introduced in 1991 by NASA/JPL, and further when data from the European Space Agency ERS-1 radar satellite became routinely available in 1991. Several airborne interferometric SAR systems are either currently operational, or are about to be introduced. Radar interferometry is a technique that allows one to map the topography of an area automatically under all weather conditions, day or night. The real power of radar interferometry is that the images and digital elevation models are automatically geometrically resampled, and could be imported into GIS systems directly after suitable reformatting. When combined with polarimetry, a technique that uses polarization diversity to gather more information about the geophysical properties of the terrain, a very rich multi-layer data set is available to the remote sensing scientist. This talk will discuss the principles of radar interferometry and polarimetry with specific application to the automatic categorization of land cover. Examples will include images acquired with the NASA/JPL AIRSAR/TOPSAR system in Australia and elsewhere.
Ausherman, D.; Larson, R.; Liskow, C.
Digitized four-channel radar images corresponding to particular areas from the Phoenix and Huntington test sites were generated in conjunction with prior experiments performed to collect X- and L-band synthetic aperture radar imagery of these two areas. The methods for generating this imagery are documented. A secondary objective was the investigation of digital processing techniques for extraction of information from the multiband radar image data. Following the digitization, the remaining resources permitted a preliminary machine analysis to be performed on portions of the radar image data. The results, although necessarily limited, are reported.
Shafique, Muhammad; Anggraeni, Dita; Bakker, Wim; van der Meijde, Mark
The intensity of earthquake triggered ground shaking is influenced by the characteristics of earthquake source, medium and site effects. These site effects are often not included in the regional ground shaking models, especially the local topography. It is being experimentally proved and noticed during many previous earthquakes, that topography has significant impact on variation of ground shaking and subsequent building damages. Majority of the previous studies investigating the topographic impact on seismic response are limited to synthetic environments or isolated hills. This study deals with exploring the impact of topography on variation of ground shaking caused by the 2005 Kashmir earthquake, at a regional scale. With the proliferation of remote sensing technologies, digital elevation models (DEMs) are freely and readily available at medium resolution, and with global cover. DEMs derived from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), with 30m resolution, and Shuttle Radar Topography Mission (SRTM), with 90m resolution, can therefore be utilized to model and predict the impact of topography on seismic response, also quickly after a seismic event. The topography of the 2005 Kashmir earthquake affected area is derived from ASTER and SRTM DEMs and analyzed using a 3D spectral finite element code (SPECFEM3D). SPECFEM3D takes into account the seismic source parameters, medium and topography to generate shake maps and earthquake simulations. The ground shaking simulations and peak ground acceleration maps were generated initially assuming the homogenous ground surface and later by including the topography to assess the role of topography in seismic amplification. Topography derived from ASTER and SRTM DEMs were simulated separately to predict the impact of DEM resolution on computed ground shaking simulations and maps. The preliminary result from the model simulations shows that seismic waves were dispersed at topographic
This is a three-dimensional perspective of the eastern front range of the Rocky Mountains, about 120 kilometers (75 miles) west of Great Falls, Montana. The image was created by combining two spaceborne radar images using a technique known as interferometry. Visualizations like this are useful to scientists because they show the shapes of the topographic features such as mountains and valleys. This technique helps to clarify the relationships of the different types of materials on the surface detected by the radar. The view is looking south-southeast. Along the right edge of the image is the valley of the north fork of the Sun River. The western edge of the Great Plains appears on the left side. The valleys in the lower center, running off into the plains on the left, are branches of the Teton River. The highest mountains are at elevations of 2,860 meters (9,390 feet), and the plains are about 1,400 meters (4,500 feet) above sea level. The dark brown areas are grasslands, bright green areas are farms, light brown, orange and purple areas are scrub and forest, and bright white and blue areas are steep rocky slopes. The two radar images were taken on successive days by the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) on board the space shuttle Endeavour in October 1994. The digital elevation map was produced using radar interferometry, a process in which radar data are acquired on different passes of the space shuttle. The two data passes are compared to obtain elevation information. Radar image data are draped over the topography to provide the color with the following assignments: red is L-band vertically transmitted, vertically received; green is C-band vertically transmitted, vertically received; and blue are the differences seen in the L-band data between the two days. This image is centered near 47.7 degrees north latitude and 112.7 degrees west longitude. No vertical exaggeration factor has been applied to the data. SIR-C/X-SAR, a
Ostro, S. J.
Radar investigations of asteroids, including observations during 1984 to 1985 of at least 8 potential targets and continued analyses of radar data obtained during 1980 to 1984 for 30 other asteroids is proposed. The primary scientific objectives include estimation of echo strength, polarization, spectral shape, spectral bandwidth, and Doppler shift. These measurements yield estimates of target size, shape, and spin vector; place constraints on topography, morphology, density, and composition of the planetary surface; yield refined estimates of target orbital parameters; and reveals the presence of asteroidal satellites.
Jara-Olivares, V. A.; Player, K.; Rodriguez-Morales, F.; Gogineni, P.
Glaciers in several parts of the world are reported to be retreating and thinning rapidly over the last decade. Radar instruments can be used to provide a wealth of information regarding the internal and basal conditions of large and small ice masses. These instruments typically operate in the VHF and UHF regions of the electromagnetic spectrum. For temperate-ice sounding, however, the high water content produces scattering and attenuation in propagating radar waves at VHF and UHF frequencies, which significantly reduce the penetration depths. Radars operating in the HF band are better suited for systematic surveys of the thickness and sub-glacial topography of temperate-ice regions. We are developing a dual-frequency Temperate-Ice-Depth Sounding Radar (TIDSoR) that can penetrate through water pockets, thus providing more accurate measurements of temperate ice properties such as thickness and basal conditions. The radar is a light-weight, low power consumption portable system for surface-based observations in mountainous terrain or aerial surveys. TIDSoR operates at two different center frequencies: 7.7 MHz and 14 MHz, with a maximum output peak power of 20 W. The transmit waveform is a digitally generated linear frequency-modulated chirp with 1 MHz bandwidth. The radar can be installed on aircrafts such as the CReSIS UAV , DCH-6 (Twin Otter), or P-3 Orion for aerial surveys, where it could be supported by the airplane power system. For surface based experiments, TIDSoR can operate in a backpack configuration powered by a compact battery system. The system can also be installed on a sled towed by a motorized vehicle, in which case the power supply can be replaced by a diesel generator. The radar consists of three functional blocks: the digital section, the radio-frequency (RF) section, and the antenna, and is designed to weigh less than 2 kg, excluding the power supply. The digital section generates the transmit waveforms as well as timing and control signals
Landais, Francois; Schmidt, Frédéric; Lovejoy, Shaun
In the last decades, a huge amount of topographic data has been obtained by several techniques (laser and radar altimetry, DTM…) for different bodies in the solar system. In each case, topographic fields exhibit an extremely high variability with details at each scale, from millimeters to thousands of kilometers. In our study, we investigate the statistical properties of the topography. Our statistical approach is motivated by the well known scaling behavior of topography that has been widely studied in the past. Indeed, scaling laws are strongly present in geophysical field and can be studied using fractal formalism. More precisely, we expect multifractal behavior in global topographic fields. This behavior reflects the high variability and intermittency observed in topographic fields that can not be generated by simple scaling models. In the multifractal formalism, each statistical moment exhibits a different scaling law characterized by a function called the moment scaling function. Previous studies were conducted at regional scale to demonstrate that topography present multifractal statistics (Gagnon et al., 2006, NPG). We have obtained similar results on Mars (Landais et al. 2015) and more recently on different body in the the solar system including the Moon, Venus and Mercury. We present the result of different multifractal approaches performed on global and regional basis and compare the fractal parameters from a body to another.
Environmental monitoring, earth-resource mapping, and military systems require broad-area imaging at high resolutions. Many times the imagery must be acquired in inclement weather or during night as well as day. Synthetic aperture radar (SAR) provides such a capability. SAR systems take advantage of the long-range propagation characteristics of radar signals and the complex information processing capability of modern digital electronics to provide high resolution imagery. SAR complements photographic and other optical imaging capabilities because of the minimum constrains on time-of-day and atmospheric conditions and because of the unique responses of terrain and cultural targets to radar frequencies. Interferometry is a method for generating a three-dimensional image of terrain. The height projection is obtained by acquiring two SAR images from two slightly differing locations. It is different from the common method of stereoscopic imaging for topography. The latter relies on differing geometric projections for triangulation to define the surface geometry whereas interferometry relies on differences in radar propagation times between the two SAR locations. This paper presents the capabilities of SAR, explains how SAR works, describes a few SAR applications, provides an overview of SAR development at Sandia, and briefly describes the motion compensation subsystem.
Watters, T.R.; Robinson, M.S.
Earth-based radar altimetry and image derived photoclinometric profiles were analyzed to examine both the long- and short-wavelength topography associated with wrinkle ridges on Mars. Photoclinometrically derived elevation data across wrinkle ridges were evaluated to determine the sensitivity of profiles to two empirical photoclinometric parameters, the horizontal digital number (HDN) and the scattered light value (SLV). The photoclinometric profiles are extremely sensitive to small variations in HDN. The sense of slope of a profile can be completely reversed over a range in HDN of as little as ?? 1. Comparably small variations in the SLV have relatively minor effects on the photoclinometrically derived elevations. The existence of elevation offsets from one side of the ridge to the other, reported in previous photoclinometric studies of martian wrinkles, were not confirmed through photoclinometry. In addition, no evidence of elevation offsets were found in Earth-based radar altimetry profiles across wrinkle ridges. In order to more accurately model wrinkle ridge topography, we controlled photoclinometrically derived elevations with long-wavelength topography obtained from the radar altimetry. The results of this study do not support kinematic models for the origin of planetary wrinkle ridges that involve deeply rooted thrust faults which separate crustal blocks at different elevations. A kinematic model involving buckling of shallow crustal layers into concentric folds that close, leading to the development of thrust faults, is consistent with wrinkle ridge morphology and terrestrial analogs. Recent geophysical studies of terrestrial analogs and the influence of shallow subsurface structures, particularly buried craters, on the localization of many wrinkle ridges on Mars suggest that thrust faults associated with the ridges are confined to the ridged plains material and do not extend into the lithosphere. Copyright 1997 by the American Geophysical Union.
Le, Charles; Spencer, Michael W.; Veilleux, Louise; Chan, Samuel; He, Yutao; Zheng, Jason; Nguyen, Kayla
An approach for algorithm specifications and development is described for SMAP's radar onboard processor with multi-stage demodulation and decimation bandpass digital filter. Point target simulation is used to verify and validate the filter design with the usual radar performance parameters. Preliminary FPGA implementation is also discussed.
Haldemann, A. F. C.; Jurgens, R. F.; Slade, M. A.
The entire fourteen-year database of Goldstone Solar System Radar Mars near-nadir radar scattering model fits is being revised using the latest topography from the Mars Global Surveyor Mars Orbiter Laser Altimeter instrument.
Brenner, Anita; Frey, Herb; DiMarzio, John; Tsaoussi, Lucia
The results of the surface topography mapping of South America during the ERS-1 geodetic mission are presented. The altimeter waveforms, the range measurement, and the internal and Doppler range corrections were obtained. The atmospheric corrections and solid tides were calculated. Comparisons between Shuttle laser altimetry and ERS-1 altimetry grid showed good agreement. Satellite radar altimetry data can be used to improve the topographic knowledge of regions for which only poor elevation data currently exist.
Bills, B. G.; Ferrari, A. J.
The paper describes an improved model of Martian global topography which has been obtained by fitting a sixteenth-degree harmonic series to occultation, radar, spectral, and photogrammetric measurements. Empirical elevation data based on photographic data are used to supplement the observations in areas without data. Values for the mean radius, the mean density, and the displacement of the center of the figure from the center of mass are presented. The reported geometric flattening is too great and the reported dynamic flattening is too small for Mars to be homogeneous and hydrostatic. Maps of the data distribution, global topography, and Bouguer gravity anomaly are interpreted in terms of a crustal thickness map which is consistent with gravity, topography, and recent preliminary Viking seismic results.
A brief history and development of the Newcastle Meteor Radar system is given. Also described are its geographical coordinates and its method of operation. The initial objective when the project was commenced was to develop an entirely digital analyzer capable of recognizing meteor echo signals and recording as many of their parameters as possible. This objective was achieved.
Slade, Martin A.; Jurgens, R. F.
Goldstar radar can provide topography 'profiles', statistical surface roughness, and radar images within a few degrees of the sub-Earth point. Goldstone/Very Large Array (VLA) bistatic radar observations can image the whole disk of Mars with integration times on the order of ten min before pixel smearing occurs. Data from all these radar techniques can be useful for observing the local surface conditions relating to landing safety issues for Mars Pathfinder.
Neish, C.; Kirk, R.; Lorenz, R.; Bray, V.; Schenk, P.; Stiles, B.; Turtle, E.; Cassini Radar Team
Unique among the icy satellites, Titan’s surface shows evidence for extensive modification by fluvial and aeolian erosion, which act to change the topography of its surface over time. Quantifying the extent of this landscape evolution is difficult, since the original, ‘non-eroded’ surface topography is generally unknown. However, fresh craters on icy satellites have a well-known shape and morphology, which has been determined from extensive studies on the airless worlds of the outer solar system (Schenk et al., 2004). By comparing the topography of craters on Titan to similarly sized, pristine analogues on airless bodies, we can obtain one of the few direct measures of the amount of erosion that has occurred on Titan. Cassini RADAR has imaged >30% of the surface of Titan, and more than 60 potential craters have been identified in this data set (Wood et al., 2010; Neish and Lorenz, 2012). Topographic information for these craters can be obtained from a technique known as ‘SARTopo’, which estimates surface heights by comparing the calibration of overlapping synthetic aperture radar (SAR) beams (Stiles et al., 2009). We present topography data for several craters on Titan, and compare the data to similarly sized craters on Ganymede, for which topography has been extracted from stereo-derived digital elevation models (Bray et al., 2012). We find that the depths of craters on Titan are generally within the range of depths observed on Ganymede, but several hundreds of meters shallower than the average (Fig. 1). A statistical comparison between the two data sets suggests that it is extremely unlikely that Titan’s craters were selected from the depth distribution of fresh craters on Ganymede, and that is it much more probable that the relative depths of Titan are uniformly distributed between ‘fresh’ and ‘completely infilled’. This is consistent with an infilling process that varies linearly with time, such as aeolian infilling. Figure 1: Depth of
Saturn's moon, Titan is a geomorphologically active planetary object, and its surface is influenced by multiple processes like impact cratering, fluvial and aeolian erosion, lacustrine processes, tectonics, cryovolcanism and mantling. Disentangling the processes that compete to shape Titan's landscape is difficult in the absence of global topography data. In this thesis, I utilize techniques in topographic statistics, fractal theory, study of terrestrial analogs and landscape evolution modeling to characterize Titan's topography and surface roughness and investigate the relative roles of surface processes in sculpting its landscape. I mapped the shorelines of 290 North Polar Titanian lakes using the Cassini Synthetic Aperture Radar dataset. The fractal dimensions of the shorelines were calculated via the divider/ruler method and box-counting method, at length scales of (1--10) km and found to average 1.27 and 1.32, respectively. The inferred power-spectral exponent of Titan's topography was found to be ≤ 2, which is lower than the values obtained from the global topography of the Earth or Venus. In order to interpret fractal dimensions of Titan's shorelines in terms of the surficial processes at work, I repeated a similar statistical analysis with 114 terrestrial analogous lakes formed by different processes, using C-band radar backscatter data from the Shuttle Radar Topography Mission (SRTM). I found different lake generation mechanisms on Earth produce 'statistically different' shorelines; however, no specific set of processes could be identified for forming Titanian lake basins. Using the Cassini RADAR altimetry data, I investigated Titan's global surface roughness and calculated median absolute slopes, average relief and Hurst exponent (H) for the surface of Titan. I detected a clear trend with latitude in these roughness parameters. Equatorial regions had the smallest slopes, lowest values of H and smallest intra-footprint relief, compared to the mid
Yoder, Paul R.; Macri, Timothy F.; Telfair, William B.; Bennett, Peter S.; Martin, Clifford A.; Warner, John W.
We describe a new electro-optical device being developed to provide precise measurements of the three-dimensional topography of the human cornea. This device, called a digital keratoscope, is intended primarily for use in preparing for and determining the effect of corneal surgery procedures such as laser refractive keratectomy, radial keratotomy or corneal transplant on the refractive power of the cornea. It also may serve as an aid in prescribing contact lenses. The basic design features of the hardware and of the associated computer software are discussed, the means for alignment and calibration are described and typical results are given.
Fountain, A.G.; Jacobel, R.W.
South Cascade Glacier, Washington, U.S.A., is one of the most extensively studied glaciers in the Western Hemisphere. In addition to mass-balance measurements, which date to 1958, numerous hydrological investigations have been carried out during the last three decades, and repeated ice-thickness determinations have been made using a variety of techniques. In the late 1960s, the basal topography was initially determined by gravitimetric methods. In the mid-1970s some of the first depth measurements using radar on temperate ice were made. The basal topography was remapped soon after from a series of point radar measurements and boreholes drilled to the glacier bottom. During the 1990s, the ice thickness was remapped using digital recording of continuous profiles that obtained over 5000 ice-thickness measurements. Profiles have been corrected for the finite beamwidth of the antenna radiation pattern and reflections in steep terrain, resulting in a significantly improved depiction of the basal surface and internal structures. The map based on our recent radar profiles confirms the large-scale features of the basal topography previously depicted and reveals more structural detail. A bright reflector was detected at the base of the glacier and could be traced in adjacent profiles. Comparison with results from water-level measurements in boreholes drilled to the bed indicates that the reflector is a subglacial conduit.
Application of ground-penetrating radar, digital optical borehole images, and cores for characterization of porosity hydraulic conductivity and paleokarst in the Biscayne aquifer, southeastern Florida, USA
This paper presents examples of ground-penetrating radar (GPR) data from two study sites in southeastern Florida where karstic Pleistocene platform carbonates that comprise the unconfined Biscayne aquifer were imaged. Important features shown on resultant GPR profiles include: (1) upward and lateral qualitative interpretative distribution of porosity and hydraulic conductivity; (2) paleotopographic relief on karstic subaerial exposure surfaces; and (3) vertical stacking of chronostratigraphic high-frequency cycles (HFCs). These characteristics were verified by comparison to rock properties observed and measured in core samples, and identified in digital optical borehole images. Results demonstrate that an empirical relation exists between measured whole-core porosity and hydraulic conductivity, observed porosity on digital optical borehole images, formation conductivity, and GPR reflection amplitudes-as porosity and hydraulic conductivity determined from core and borehole images increases, formation conductivity increases, and GPR reflection amplitude decreases. This relation allows for qualitative interpretation of the vertical and lateral distribution of porosity and hydraulic conductivity within HFCs. Two subtidal HFCs in the uppermost Biscayne aquifer have significantly unique populations of whole-core porosity values and vertical hydraulic conductivity values. Porosity measurements from one cycle has a median value about two to three times greater than the values from the other HFC, and median values of vertical hydraulic-conductivity about three orders of magnitude higher than the other HFC. The HFC with the higher porosity and hydraulic conductivity values is shown as a discrete package of relatively low-amplitude reflections, whereas the HFC characterized by lower porosity and hydraulic-conductivity measurements is expressed by higher amplitude reflections. Porosity and hydraulic-conductivity values measured from whole-core samples, and vuggy porosity
During the fifty years since its initial development as a means of providing early warning of airborne attacks against allied countries during World War II, radar systems have developed to the point of being highly mobile and versatile systems capable of supporting a wide variety of remote sensing applications. Instead of being tied to stationary land-based sites, radar systems have found their way into highly mobile land vehicles as well as into aircraft, missiles, and ships of all sizes. Of all these applications, however, the most exciting revolution has occurred in the airborne platform arena where advanced technology radars can be found in all shapes and sizes...ranging from the large AWACS and Joint STARS long range surveillance and targeting systems to small millimeter wave multi-spectral sensors on smart weapons that can detect and identify their targets through the use of highly sophisticated digital signal processing hardware and software. This paper presents an overview of these radar applications with the emphasis on modern airborne sensors that span the RF spectrum. It will identify and describe the factors that influence the parameters of low frequency and ultra wide band radars designed to penetrate ground and dense foliage environments and locate within them buried mines, enemy armor, and other concealed or camouflaged weapons of war. It will similarly examine the factors that lead to the development of airborne radar systems that support long range extended endurance airborne surveillance platforms designed to detect and precision-located both small high speed airborne threats as well as highly mobile time critical moving and stationary surface vehicles. The mission needs and associated radar design impacts will be contrasted with those of radar systems designed for high maneuverability rapid acquisition tactical strike warfare platforms, and shorter range cued air-to-surface weapons with integral smart radar sensors.
A method is presented for the determination of radar frequency radiation power densities that the PAVE PAWS radar system could produce in its air and ground environment. The effort was prompted by the concern of the people in the vicinity of OTIS AFB MA and BEALE AFB CA about the possible radar frequency radiation hazard of the PAVE PAWS radar. The method is based on the following main assumptions that: (a) the total field can be computed as the vector summation of the individual fields due to each antenna element; (b) the individual field can be calculated using distances for which the field point is in the far field of the antenna element. An RFR computer program was coded for the RADC HE 6180 digital computer and exercised to calculate the radiation levels in the air and ground space for the present baseline and the possible Six DB and 10 DB growth systems of the PAVE PAWS radar system at OTIS AFB MA. The average radiation levels due to the surveillance fence were computed for three regions: in the air space in front of the radar, at the radar hazard fence at OTIS AFB MA and at representative ground points in the OTIS AFB vicinity. It was concluded that the radar frequency radiation of PAVE PAWS does not present a hazard to personnel provided there is no entry to the air hazard zone or to the area within the hazard fence. The method developed offers a cost effective way to determine radiation levels from a phased array radar especially in the near field and transition regions.
Tuccari, G.; Bezrukovs, Vl.; Nechaeva, M.
A digital base band converter (DBBC) system has been developed by the Istituto di Radioastronomia (Noto, Italy) for increasing the sensitivity of European VLBI Network (EVN) by expanding the full observed bandwidth using numerical methods. The output data rate of this VLBI-backend is raised from 1 to 4 Gbps for each radiotelescope. All operations related to the signal processing (frequency translation, amplification, frequency generation with local oscillators, etc.) are transferred to the digital domain, which allows - in addition to well-known advantages coming from digital technologies - achieving better repeatability, precision, simplicity, etc. The maximum input band of DBBC system is 3.5 GHz, and the instantaneous bandwidth is up to 1 GHz for each radio frequency/intermediate frequency (RF/IF) out of the eight possible. This backend is a highly powerful platform for other radioastronomy applications, and a number of additional so-called personalities have been developed and used. This includes PFB (polyphase filter bank) receivers and Spectra for high resolution spectroscopy. An additional new development with the same aim - to use the DBBC system as a multi-purpose backend - is related to the bi-static radar observations including Radar VLBI. In such observations it is possible to study the population of space debris, with detection of even centimetre class fragments. A powerful transmitter is used to illuminate the sky region to be analyzed, and the echoes coming from known or unknown objects are reflected to one or more groundbased telescopes thus producing a single-dish or interferometric detection. The DBBC Radar VLBI personality is able to realize a high-resolution spectrum analysis, maintaining in the central area the echo signal at the expected frequency including the Doppler shift of frequency. For extremely weak signals a very large integration time is needed, so for this personality different input parameters are provided. The realtime information
Rao, Kota S.; Al Jassar, Hala K.
The aim of this paper is to analyze the errors in the Digital Elevation Models (DEMs) derived through repeat pass SAR interferometry (InSAR). Out of 29 ASAR images available to us, 8 are selected for this study which has unique data set forming 7 InSAR pairs with single master image. The perpendicular component of baseline (B highmod) varies between 200 to 400 m to generate good quality DEMs. The Temporal baseline (T) varies from 35 days to 525 days to see the effect of temporal decorrelation. It is expected that all the DEMs be similar to each other spatially with in the noise limits. However, they differ very much with one another. The 7 DEMs are compared with the DEM of SRTM for the estimation of errors. The spatial and temporal distribution of errors in the DEM is analyzed by considering several case studies. Spatial and temporal variability of precipitable water vapour is analysed. Precipitable water vapour (PWV) corrections to the DEMs are implemented and found to have no significant effect. The reasons are explained. Temporal decorrelation of phases and soil moisture variations seem to have influence on the accuracy of the derived DEM. It is suggested that installing a number of corner reflectors (CRs) and the use of Permanent Scatter approach may improve the accuracy of the results in desert test sites.
Fu, Lee-Lueng; Leconte, Denis; Pihos, Greg; Davidson, Roger; Kruizinga, Gerhard; Tapley, Byron
Three video loops showing various aspects of the dynamic ocean topography obtained from the TOPEX/POSEIDON radar altimetry data will be presented. The first shows the temporal change of the global ocean topography during the first year of the mission. The time-averaged mean is removed to reveal the temporal variabilities. Temporal interpolation is performed to create daily maps for the animation. A spatial smoothing is also performed to retain only the large-sale features. Gyre-scale seasonal changes are the main features. The second shows the temporal evolution of the Gulf Stream. The high resolution gravimetric geoid of Rapp is used to obtain the absolute ocean topography. Simulated drifters are used to visualize the flow pattern of the current. Meanders and rings of the current are the main features. The third is an animation of the global ocean topography on a spherical earth. The JGM-2 geoid is used to obtain the ocean topography...
Erenskjold Moeslund, Jesper; Klith Bøcher, Peder; Svenning, Jens-Christian; Mølhave, Thomas; Arge, Lars
This study examines the potential impact of 21st century sea-level rise on Aarhus, the second largest city in Denmark, emphasizing the economic risk to the city's real estate. Furthermore, it assesses which possible adaptation measures that can be taken to prevent flooding in areas particularly at risk from flooding. We combine a new national Digital Elevation Model in very fine resolution (~2 meter), a new highly computationally efficient flooding algorithm that accurately models the influence of barriers, and geospatial data on real-estate values to assess the economic real-estate risk posed by future sea-level rise to Aarhus. Under the A2 and A1FI (IPCC) climate scenarios we show that relatively large residential areas in the northern part of the city as well as areas around the river running through the city are likely to become flooded in the event of extreme, but realistic weather events. In addition, most of the large Aarhus harbour would also risk flooding. As much of the area at risk represent high-value real estate, it seems clear that proactive measures other than simple abandonment should be taken in order to avoid heavy economic losses. Among the different possibilities for dealing with an increased sea level, the strategic placement of flood-gates at key potential water-inflow routes and the construction or elevation of existing dikes seems to be the most convenient, most socially acceptable, and maybe also the cheapest solution. Finally, we suggest that high-detail flooding models similar to those produced in this study will become an important tool for a climate-change-integrated planning of future city development as well as for the development of evacuation plans.
Michelson, M.; Shrader, W. W.; Wieler, J. G.
The terminal Doppler weather radar (TDWR) system, now under development, will provide automatic detection of microbursts and low-level wind shear. This paper discusses the TDWR performance parameters and describes its structural elements, including the antenna subsystem, the transmitter, the receiver/exciter, the digital signal processor, and the radar product generator/remote monitoring subsystem. Attention is also given to the processes of the base data formation, point target removal, signal-to-noise thresholding, and velocity de-aliasing and to the TDWR algorithms and displays. A schematic diagram of the TDWR system is presented.
Esmaeilzade, M.; Amini, J.; Zakeri, S.
Due to the SAR1 geometry imaging, SAR images include geometric distortions that would be erroneous image information and the images should be geometrically calibrated. As the radar systems are side looking, geometric distortion such as shadow, foreshortening and layover are occurred. To compensate these geometric distortions, information about sensor position, imaging geometry and target altitude from ellipsoid should be available. In this paper, a method for geometric calibration of SAR images is proposed. The method uses Range-Doppler equations. In this method, for the image georeferencing, the DEM2 of SRTM with 30m pixel size is used and also exact ephemeris data of the sensor is required. In the algorithm proposed in this paper, first digital elevation model transmit to range and azimuth direction. By applying this process, errors caused by topography such as foreshortening and layover are removed in the transferred DEM. Then, the position of the corners on original image is found base on the transferred DEM. Next, original image registered to transfer DEM by 8 parameters projective transformation. The output is the georeferenced image that its geometric distortions are removed. The advantage of the method described in this article is that it does not require any control point as well as the need to attitude and rotational parameters of the sensor. Since the ground range resolution of used images are about 30m, the geocoded images using the method described in this paper have an accuracy about 20m (subpixel) in planimetry and about 30m in altimetry. 1 Synthetic Aperture Radar 2 Digital Elevation Model
Mason, David; Garcia-Pintado, Javier; Cloke, Hannah; Dance, Sarah
A basic data requirement of a river flood inundation model is a Digital Terrain Model (DTM) of the reach being studied. The scale at which modeling is required determines the accuracy required of the DTM. For modeling floods in urban areas, a high resolution DTM such as that produced by airborne LiDAR (Light Detection And Ranging) is most useful, and large parts of many developed countries have now been mapped using LiDAR. In remoter areas, it is possible to model flooding on a larger scale using a lower resolution DTM, and in the near future the DTM of choice is likely to be that derived from the TanDEM-X Digital Elevation Model (DEM). A variable-resolution global DTM obtained by combining existing high and low resolution data sets would be useful for modeling flood water dynamics globally, at high resolution wherever possible and at lower resolution over larger rivers in remote areas. A further important data resource used in flood modeling is the flood extent, commonly derived from Synthetic Aperture Radar (SAR) images. Flood extents become more useful if they are intersected with the DTM, when water level observations (WLOs) at the flood boundary can be estimated at various points along the river reach. To illustrate the utility of such a global DTM, two examples of recent research involving WLOs at opposite ends of the spatial scale are discussed. The first requires high resolution spatial data, and involves the assimilation of WLOs from a real sequence of high resolution SAR images into a flood model to update the model state with observations over time, and to estimate river discharge and model parameters, including river bathymetry and friction. The results indicate the feasibility of such an Earth Observation-based flood forecasting system. The second example is at a larger scale, and uses SAR-derived WLOs to improve the lower-resolution TanDEM-X DEM in the area covered by the flood extents. The resulting reduction in random height error is significant.
This is a three-dimensional perspective view of Missoula, Montana, created by combining two spaceborne radar images using a technique known as interferometry. Visualizations like this are useful because they show scientists the shapes of the topographic features such as mountains and valleys. This technique helps to clarify the relationships of the different types of materials on the surface detected by the radar. The view is looking north-northeast. The blue circular area at the lower left corner is a bend of the Bitterroot River just before it joins the Clark Fork, which runs through the city. Crossing the Bitterroot River is the bridge of U.S. Highway 93. Highest mountains in this image are at elevations of 2,200 meters (7,200 feet). The city is about 975 meters (3,200 feet) above sea level. The bright yellow areas are urban and suburban zones, dark brown and blue-green areas are grasslands, bright green areas are farms, light brown and purple areas are scrub and forest, and bright white and blue areas are steep rocky slopes. The two radar images were taken on successive days by the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) onboard the space shuttle Endeavour in October 1994. The digital elevation map was produced using radar interferometry, a process in which radar data are acquired on different passes of the space shuttle. The two data passes are compared to obtain elevation information. Radar image data are draped over the topography to provide the color with the following assignments: red is L-band vertically transmitted, vertically received; green is C-band vertically transmitted, vertically received; and blue are differences seen in the L-band data between the two days. This image is centered near 46.9 degrees north latitude and 114.1 degrees west longitude. No vertical exaggeration factor has been applied to the data. SIR-C/X-SAR, a joint mission of the German, Italian and United States space agencies, is part of NASA
Matsuoka, Kenichi; Saito, Ryoji; Naruse, Renji
We have developed a novel ice-penetrating radar system that can be carried on a backpack. Including batteries for a 3 hour continuous measurement, the total weight is 13 kg. In addition, it operates reliably down to -25°C, has a low power consumption of 24 W, and is semi-waterproof. The system has a built-in-one controller with a high-brightness display for reading data quickly, a receiver with 12-bit digitizing, and a 1 kV pulse transmitter in which the pulse amplitude varies by <0.2%. Optical communications between components provides low-noise data acquisition and allows synchronizing of the pulse transmission with sampling. Measurements with the system revealed the 300 m deep bed topography of a temperate valley glacier in the late ablation season.
Parra, C. G.; Forsythe, R. G.; Parsons, C. L.
The dynamic topography of the sea surface was measured. The radar altimeter measurements yield average ocean topographic data which are mapped. Seasonal deviations from a 3 year mean topography are presented. The altimeters are also instrumented with sample and hold gates which provide information about the shape and amplitude of the return waveform. Parameters including ocean surface wind speed and the significant wave height are determined. One hundred eighty six wind speed and significant wave height histograms are presented.
Mancus, E. F.; Baker, L. H.
The baseline performance characteristics of the moving target detector (MTD) and radar data acquisition system (RDAS) as an integral part of the Mode S sensor, were determined. The MTD and RDAS were separately evaluated to determine their capability to provide radar data suitable for utilization by the Mode S sensor and automated radar terminal system (ARTS). The design modifications made to the Mode S sensor to provide the capability of interfacing to either an MTD or RDAS were evaluated to determine if they were in compliance with the Federal Aviation Administration engineering requirement, FAA-ER-240-26. Radar baseline technical performance data was provided to characterize the MTD, RDAS, Mode S, and ARTS. The minimum radar tracking requirements are studied to determine if they are adequate to provide reliable radar track data to an air traffic control facility. It was concluded that the Mode S sensor, when integrated with an MTD-2 radar digitizer, can provide reliable primary radar track data to the ARTS III system for automated radar track acquisition.
Dynamic Topography Revisited Dynamic topography is usually considered to be one of the trinity of contributing causes to the Earth's non-hydrostatic topography along with the long-term elastic strength of the lithosphere and isostatic responses to density anomalies within the lithosphere. Dynamic topography, thought of this way, is what is left over when other sources of support have been eliminated. An alternate and explicit definition of dynamic topography is that deflection of the surface which is attributable to creeping viscous flow. The problem with the first definition of dynamic topography is 1) that the lithosphere is almost certainly a visco-elastic / brittle layer with no absolute boundary between flowing and static regions, and 2) the lithosphere is, a thermal / compositional boundary layer in which some buoyancy is attributable to immutable, intrinsic density variations and some is due to thermal anomalies which are coupled to the flow. In each case, it is difficult to draw a sharp line between each contribution to the overall topography. The second definition of dynamic topography does seem cleaner / more precise but it suffers from the problem that it is not measurable in practice. On the other hand, this approach has resulted in a rich literature concerning the analysis of large scale geoid and topography and the relation to buoyancy and mechanical properties of the Earth [e.g. refs 1,2,3] In convection models with viscous, elastic, brittle rheology and compositional buoyancy, however, it is possible to examine how the surface topography (and geoid) are supported and how different ways of interpreting the "observable" fields introduce different biases. This is what we will do. References (a.k.a. homework)  Hager, B. H., R. W. Clayton, M. A. Richards, R. P. Comer, and A. M. Dziewonski (1985), Lower mantle heterogeneity, dynamic topography and the geoid, Nature, 313(6003), 541-545, doi:10.1038/313541a0.  Parsons, B., and S. Daly (1983), The
Shapiro, Irwin I.
We aid in a study of the solar system by means of ground-based radar. We have concentrated on (1) developing the ephemerides needed to acquire radar data at Arecibo Observatory and (2) analyzing the resultant data to: test fundamental laws of gravitation; determine the size, shape, topography, and spin vectors of the targets; and study the surface properties of these objects, through their scattering law and polarization characteristics.
We consider a mathematical model of synthetic aperture radar with a known, possibly non-flat, topography. In this context we consider the problem of recovering the wavefront set of the ground reflectivity, given radar data measured along a curved flight path. We show that if singularities are located at ‘mirror points,’ then the resulting data may be smooth; in effect, the singularities ‘cancel.’ With a flat topography, these mirror points are always discrete, but we show that in a non-flat topography there may be infinite families of mirror points.
Balchin, A. A.
Discusses some aspects in X-ray topography, including formation of dislocations, characteristics of stacking faults, x-ray contrast in defect inspection, Berg-Barrett technique, and Lang traversing crystal and Borrmann's methods. (CC)
Moreno Yeras, A.
For several decades measurement optical techniques have been used in different branches of Science and Technology and in medicine. One of these techniques is the so-called Moire topography that allows the accurate measurement of different parts of the human body topography. This investigation presents the measurement of topographies of teeth and gums using an automated system of shadow moire, with which precision can be reached up to the order of the microns by the phase shift instrumentation in an original way. Advantages and disadvantages of using the Moire topography and its comparison with other techniques used in the optical metrology are presented. Also, some positive and negative aspects of the implementation of this technique are shown in dentistry.
Peters, M. E.; Blankenship, D. D.; Morse, D. L.; Holt, J. W.; Kempf, S. D.; Richter, T. G.; Falola, B.; Oliason, S.
Lake Vostok was discovered using airborne ice-sounding radar in East Antarctica during the mid 1970's, but interest in this largest known subglacial lake has increased in recent years. Frozen microbial discoveries from ice cores taken just above Lake Vostok suggest its potential for being an isolated biological ecosystem. Also, the lake's unique combination of glaciologic, hydrologic and geological processes make it a possible terrestrial analogue for sub-ice water on other planetary bodies. Satellite radar has mapped the spatial extent of the lake from surface topography, and Russian ground traverses have gathered radar and seismic data along select profiles, but the full subglacial environment has remained uncharted. In response to a proposal by R.E. Bell and M. Studinger at Lamont Doherty Earth Observatory, the University of Texas Institute for Geophysics (UTIG) conducted an airborne geophysical survey over Lake Vostok and its surroundings during the 2000/01 field season. The survey included 21,000 line-km of geophysical observations with a line spacing of 7.5 km and a tie-line spacing of 11.25 or 22.5 km. The instrument suite included incoherent ice-sounding radar, laser altimetry, and precise GPS positioning and navigation, as well as airborne gravity and magnetics measurements. The radar system consisted of a 60 MHz, 8000 watt peak power transmitter operating in pulsed continuous-wave mode at 12.5 kHz (with 250 ns pulse width), a log-detection incoherent receiver (with 80 dB dynamic range), and a signal digitizer with a unique capability to average signals rapidly. Incoherent radar observations constructed from 2048 averaged transmissions occurred roughly every 12 m along-track. Ice thicknesses in excess of 4000 m were routinely sounded over Lake Vostok using this system. In addition to the incoherent radar, a new acquisition system was developed on an experimental basis to coherently integrate radar signals utilizing synthetic aperture radar techniques
Chapman, B.; Freeman, A.
A perspective of NASA's Jet Propulsion Laboratory as a center of planetary exploration, and its involvement in studying the earth from space is given. Remote sensing, radar maps, land topography, snow cover properties, vegetation type, biomass content, moisture levels, and ocean data are items discussed related to earth orbiting satellite imaging radar. World Wide Web viewing of this content is discussed.
This radar image acquired by SRTM shows an area south of the Sao Francisco River in Brazil. The area is predominantly scrub forest. Areas such as these are difficult to map by traditional methods because of frequent cloud cover and local inaccessibility. Image brightness differences in this image are caused by differences in vegetation type and density. Tributaries of the Sao Francisco are visible in the upper right. The Sao Francisco River is a major source of water for irrigation and hydroelectric power. Mapping such regions will allow scientists to better understand the relationships between flooding cycles, forestation and human influences on ecosystems.This radar image was obtained by the Shuttle Radar Topography Mission as part of its mission to map the Earth's topography. The image was acquired by just one of SRTM's two antennas, and consequently does not show topographic data but only the strength of the radar signal reflected from the ground. This signal, known as radar backscatter, provides insight into the nature of the surface, including its roughness, vegetation cover, and urbanization.The Shuttle Radar Topography Mission (SRTM), launched on February 11, 2000, uses the same radar instrument that comprised the Spaceborne Imaging Radar-C/X-Band Synthetic Aperture Radar (SIR-C/X-SAR) that flew twice on the Space Shuttle Endeavour in 1994. The mission is designed to collect three-dimensional measurements of the Earth's surface. To collect the 3-D data, engineers added a 60-meter-long (200-foot) mast, an additional C-band imaging antenna and improved tracking and navigation devices. The mission is a cooperative project between the National Aeronautics and Space Administration (NASA), the National Imagery and Mapping Agency (NIMA) and the German and Italian space agencies. It is managed by NASA's Jet Propulsion Laboratory, Pasadena, CA, for NASA's Earth Science Enterprise, Washington, DC.
This is a three-dimensional view of Isabela, one of the Galapagos Islands located off the western coast of Ecuador, South America. This view was constructed by overlaying a Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) image on a digital elevation map produced by TOPSAR, a prototype airborne interferometric radar which produces simultaneous image and elevation data. The vertical scale in this image is exaggerated by a factor of 1.87. The SIR-C/X-SAR image was taken on the 40th orbit of space shuttle Endeavour. The image is centered at about 0.5 degree south latitude and 91 degrees west longitude and covers an area of 75 by 60 kilometers (47 by 37 miles). The radar incidence angle at the center of the image is about 20 degrees. The western Galapagos Islands, which lie about 1,200 kilometers (750 miles)west of Ecuador in the eastern Pacific, have six active volcanoes similar to the volcanoes found in Hawaii and reflect the volcanic processes that occur where the ocean floor is created. Since the time of Charles Darwin's visit to the area in 1835, there have been more than 60 recorded eruptions on these volcanoes. This SIR-C/X-SAR image of Alcedo and Sierra Negra volcanoes shows the rougher lava flows as bright features, while ash deposits and smooth pahoehoe lava flows appear dark. Vertical exaggeration of relief is a common tool scientists use to detect relationships between structure (for example, faults, and fractures) and topography. Spaceborne Imaging Radar-C and X-Synthetic Aperture Radar (SIR-C/X-SAR) is part of NASA's Mission to Planet Earth. The radars illuminate Earth with microwaves allowing detailed observations at any time, regardless of weather or sunlight conditions. SIR-C/X-SAR uses three microwave wavelengths: L-band (24 cm), C-band (6 cm) and X-band (3 cm). The multi-frequency data will be used by the international scientific community to better understand the global environment and how it is changing. The SIR-C/X-SAR data
Zebker, Howard A.; Villasensor, John
A radar interferometric technique for topographic mapping of surfaces promises a high resolution, globally consistent approach to generation of digital elevation models. One implementation approach, that of utilizing a single SAR system in a nearly repeating orbit, is attractive not only for cost and complexity reasons but also in that it permits inference of changes in the surface over the orbit repeat cycle from the correlation properties of the radar echoes. The various sources contributing to the echo correlation statistics are characterized, and the term which most closely describes surficial change is isolated. There is decorrelation increasing with time, but digital terrain model generation remains feasible.
Esteban-Fernandez, Daniel; Rodriquez, Ernesto; Peral, Eva; Clark, Duane I.; Wu, Xiaoqing
An interferometric synthetic aperture radar (SAR) onboard processor concept and algorithm has been developed for the Ka-band radar interferometer (KaRIn) instrument on the Surface and Ocean Topography (SWOT) mission. This is a mission- critical subsystem that will perform interferometric SAR processing and multi-look averaging over the oceans to decrease the data rate by three orders of magnitude, and therefore enable the downlink of the radar data to the ground. The onboard processor performs demodulation, range compression, coregistration, and re-sampling, and forms nine azimuth squinted beams. For each of them, an interferogram is generated, including common-band spectral filtering to improve correlation, followed by averaging to the final 1 1-km ground resolution pixel. The onboard processor has been prototyped on a custom FPGA-based cPCI board, which will be part of the radar s digital subsystem. The level of complexity of this technology, dictated by the implementation of interferometric SAR processing at high resolution, the extremely tight level of accuracy required, and its implementation on FPGAs are unprecedented at the time of this reporting for an onboard processor for flight applications.
Frost, E.; Downing, J.
Tests were conducted at several sites over the coastal lowlands of New Jersey and over a region of high plains and low mountains in Oklahoma. In one area, a salt marsh in New Jersey, extensive ground tests were combined with laboratory data on expected insect backscatter to arrive at an extremely convincing model of the insect origin of most Dot Angels. A great deal of insight was studied from radar on the buildup and dispersal of insect swarms, since radar can follow where other means of trapping and observation cannot. Data on large-scale behavior as a function of wind and topography are presented. Displayed techniques which show individual or small swarm motion within some larger cloud or mass, or which can show the overall motion over great distances were developed. The influence of wind and terrain on insect motion and dispersal is determined from radar data.
Thesenga, David; Town, James
In February 2000, the Space Shuttle Endeavour flew a specially modified radar system during an 11-day mission. The purpose of the multinational Shuttle Radar Topography Mission (SRTM) was to "obtain elevation data on a near-global scale to generate the most complete high-resolution digital topographic database of Earth" by using radar interferometry. The data and resulting products are now publicly available for download and give a view of the landscape removed of vegetation, buildings, and other structures. This new view of the Earth's topography allows us to see previously unmapped or poorly mapped regions of the Earth as well as providing a level of detail that was previously unknown using traditional topographic mapping techniques. Understanding and appreciating the geographic terrain is a complex but necessary requirement for middle school aged (11-14yo) students. Abstract in nature, topographic maps and other 2D renderings of the Earth's surface and features do not address the inherent spatial challenges of a concrete-learner and traditional methods of teaching can at times exacerbate the problem. Technological solutions such as 3D-imaging in programs like Google Earth are effective but lack the tactile realness that can make a large difference in learning comprehension and retention for these young students. First developed in the 1980's, 3D printers were not commercial reality until recently and the rapid rise in interest has driven down the cost. With the advent of sub US1500 3D printers, this technology has moved out of the high-end marketplace and into the local office supply store. Schools across the US and elsewhere in the world are adding 3D printers to their technological workspaces and students have begun rapid-prototyping and manufacturing a variety of projects. This project attempted to streamline the process of transforming SRTM data from a GeoTIFF format by way of Python code. The resulting data was then inputted into a CAD-based program for
Moore, R. K.; Eckerman, J.; Meneghini, R.; Atlas, D.; Boerner, W. M.; Cherry, S.; Clark, J. F.; Doviak, R. J.; Goldhirsh, J.; Lhermitte, R. M.
The spaceborne radar panel considered how radar could be used to measure precipitation from satellites. The emphasis was on how radar could be used with radiometry (at microwave, visible (VIS), and infrared (IR) wavelengths) to reduce the uncertainties of measuring precipitation with radiometry alone. In addition, the fundamental electromagnetic interactions involved in the measurements were discussed to determine the key work areas for research and development to produce effective instruments. Various approaches to implementing radar systems on satellites were considered for both shared and dedicated instruments. Finally, a research and development strategy was proposed for establishing the parametric relations and retrieval algorithms required for extracting precipitation information from the radar and associated radiometric data.
Holt, J. W.; Peters, M. E.; Kempf, S. D.; Morse, D. L.; Blankenship, D. D.
radar sounding data. The first technique simulates radar data using a digital elevation model (DEM) of surface topography to predict the location and shape of surface echoes in the radar data. This is complemented by the cross-track migration of radar echoes onto the surface. These migrated echoes are superimposed on imagery in order to correlate them with potential surface sources. Using these techniques enabled us to identify a number of echoes in a 24-km segment of the Dry Valleys flight path as arising from the surface and to identify subsurface echoes under the main trunk of Taylor Glacier and possibly multiple reflectors beneath the toe of Taylor Glacier. Surface-based radar confirms the thickness of the glacier at three crossing points. In the ice-free section of the test segment no real subsurface reflectors were found, indicating that the electromagnetic properties of the ground there do not allow significant radar penetration at 60 MHz and/or no radar-significant subsurface interfaces exist. These results illustrate the importance of using complementary techniques, the usefulness of a DEM, and the limitations of single-pass radar sounding data. Advanced processing techniques utilizing radar phase information show promise for achieving better clutter removal for single-pass data. Multi-pass data that we recently collected in the Dry Valleys should allow for the development of techniques to reduce or eliminate the need for a surface elevation model.
Kim, Jungrack; Wan, Wanhui
altimetry and stereo analyses have been widely used for this purpose and achieved high quality 3D topographic data over various planetary surfaces such as Venus, Mercury, Moon and Mars. However, in contrast with inner plane and satellite, the base data sets to compose digital topography over outer planets and satellites are very limited. Titan, the largest satellite of Saturn has also too limited data inventory to achieve sufficient spatial resolution in topographic data, in spite of increasing interests about the detailed topography owing to the recent interesting discoveries on methane fluvial system, aeolian geomorphologies and possible tectonic activity. Therefore the endeavours to increase the coverage of digital topography employing radargrammetry (Kirk et al. 2009), radar altimetry (Elachi, et al. 2005) and SARtopo (Stiles et al. 2009) have been actively conducted. Although these efforts result in the construction of a global topographic map, the consequent spatial resolutions of global topography is still poor and cover the resolution ranges from 520m to 1700m (Lorentz et al. 2013). In this study, we tried to improve the coverage and the quality of Titan digital terrain model employing approaches as follows; 1) A semi-automated stereo matching scheme manipulating low signal-to-noise SAR image pair incorporating adaptive filtering and base topography, 2) the geodetic control improvement of stereo SAR pair based on generic sensor model and tie points, 3) the introduction of radarclinometry to refine the topography from stereo analyses. The developed scheme was applied for a few testing areas especially over the fluvial channels and the lakes which are only the acting hydrological system in solar system except terrestrial one and well covered by SAR images. Considering geodetic controllability over the SAR images is better than 3-4 pixel when the images were projected onto zero height plane, it was evaluated that the stereo processed using the generic sensor
Zebker, Howard A. (Inventor); Held, Daniel N. (Inventor); Vanzyl, Jakob J. (Inventor); Dubois, Pascale C. (Inventor); Norikane, Lynne (Inventor)
Two alternative methods are presented for digital reduction of synthetic aperture multipolarized radar data using scattering matrices, or using Stokes matrices, of four consecutive along-track pixels to produce averaged data for generating a synthetic polarization image.
Zargli, Eleni; Liodakis, Stelios; Kyriakidis, Phaedon; Savvaidis, Alexandros
Continuous topography from Digital Elevation Model (DEM) data is frequently segmented into terrain classes based on local morphological characteristics of terrain elevation, e.g., local slope gradient and convexity. The resulting classes are often used as proxies for the average shear wave velocity up to 30 m, and the determination of ground types as required by the Eurocode (EC8) for computing elastic design spectra. In this work, we investigate the links between terrain related variables, particularly slope gradient, extracted for the area of Greece from the Shuttle Radar Topography Mission (SRTM) 30 arc second global topographic data available from the United States Geological Survey (USGS), with: (a) the global terrain classification product of Iwahashi and Pike (2007) in which 16 terrain types are identified for the same spatial resolution, and (b) information on geological units extracted at the same resolution from the geological map of Greece at a scale of 1/500000 as published from the Institute of Geology and Mineral Exploration (IGME). An interpretation of these links is presented within the context of understanding the reliability of using geology, slope and terrain classes for site characterizations of earthquake risk in a high seismicity area like Greece. Our results indicate that slope is a somewhat biased proxy for solid rocks, whereas in Alluvial deposits the distance to and type of the nearest geological formation appears to provide qualitative information on the size of the sedimentary deposit.
These two images of the eastern part of the island of Oahu, Hawaii provide information on regional topography and show the relationship between urban development and sensitive ecosystems. On the left is a topographic radar image collected by the Shuttle Radar Topography Mission (SRTM.) On the right is an optical image acquired by a digital camera on the Space Shuttle Endeavour, which carried SRTM. Features of interest in this scene include Diamond Head (an extinct volcano at the lower center), Waikiki Beach (just left of Diamond Head), the Punchbowl National Cemetery (another extinct volcano, at the foot of the Koolau Mountains), downtown Honolulu and Honolulu airport (lower left of center), and Pearl Harbor (at the left edge.)The topography shows the steep, high central part of the island surrounded by flatter coastal areas. The optical image shows the urban areas and a darker, forested region on the mountain slopes. The clouds in the optical image and the black areas on the topographic image are both a result of the steep topography. In this tropical region, high mountain peaks are usually covered in clouds. These steep peaks also cause shadows in the radar data, resulting in missing data 'holes.' A second pass over the island was obtained by SRTM and will be used to fill in the holes.The left image combines two types of SRTM data. Brightness corresponds to the strength of the radar signal reflected from the ground, while colors show the elevation. Each color cycle (from pink through blue and back to pink) represents 400 meters (1,300 feet) of elevation difference, like the contour lines on a topographic map. This image contains about 2,400 meters (8,000 feet) of total relief. The optical image was acquired by the Shuttle Electronic Still Camera with a lens focal length of 64 millimeters (2.5 inches) for the Earth Knowledge Acquired by Middle school students (EarthKAM) project. EarthKAM has flown on five space shuttle missions since 1996. Additional
Moller, Delwyn K.; Sadowy, Gregory A.; Rignot, Eric J.; Madsen, Soren N.
A report discusses Ka-band (35-GHz) radar for mapping the surface topography of glaciers and ice sheets at high spatial resolution and high vertical accuracy, independent of cloud cover, with a swath-width of 70 km. The system is a single- pass, single-platform interferometric synthetic aperture radar (InSAR) with an 8-mm wavelength, which minimizes snow penetration while remaining relatively impervious to atmospheric attenuation. As exhibited by the lower frequency SRTM (Shuttle Radar Topography Mission) AirSAR and GeoSAR systems, an InSAR measures topography using two antennas separated by a baseline in the cross-track direction, to view the same region on the ground. The interferometric combination of data received allows the system to resolve the pathlength difference from the illuminated area to the antennas to a fraction of a wavelength. From the interferometric phase, the height of the target area can be estimated. This means an InSAR system is capable of providing not only the position of each image point in along-track and slant range as with a traditional SAR but also the height of that point through interferometry. Although the evolution of InSAR to a millimeter-wave center frequency maximizes the interferometric accuracy from a given baseline length, the high frequency also creates a fundamental problem of swath coverage versus signal-to-noise ratio. While the length of SAR antennas is typically fixed by mass and stowage or deployment constraints, the width is constrained by the desired illuminated swath width. As the across-track beam width which sets the swath size is proportional to the wavelength, a fixed swath size equates to a smaller antenna as the frequency is increased. This loss of antenna size reduces the two-way antenna gain to the second power, drastically reducing the signal-to-noise ratio of the SAR system. This fundamental constraint of high-frequency SAR systems is addressed by applying digital beam-forming (DBF) techniques to
Landais, F.; Schmidt, F.; Lovejoy, S.
In the present study, we investigate the scaling properties of the topography of Mars. Planetary topographic fields are well known to roughly exhibit (mono)fractal behavior. Indeed, the fractal formalism reproduces much of the variability observed in topography. Still, a single fractal dimension is not enough to explain the huge variability and intermittency. Previous studies have claimed that fractal dimensions might be different from one region to another, excluding a general description at the planetary scale. In this article, we analyze the Martian topographic data with a multifractal formalism to study the scaling intermittency. In the multifractal paradigm, the apparent local variation of the fractal dimension is interpreted as a statistical property of multifractal fields. We analyze the topography measured with the Mars Orbiter Laser altimeter (MOLA) at 300 m horizontal resolution, 1 m vertical resolution. We adapted the Haar fluctuation method to the irregularly sampled signal. The results suggest a multifractal behavior from the planetary scale down to 10 km. From 10 to 300 m, the topography seems to be simple monofractal. This transition indicates a significant change in the geological processes governing the Red Planet's surface.
Landais, F.; Schmidt, F.; Lovejoy, S.
In the present study, we investigate the scaling properties of the topography of Mars. Planetary topographic fields are well known to roughly exhibit (mono)fractal behavior. Indeed, the fractal formalism is reproduces much of the variability observed in topography. Still, a single fractal dimension is not enough to explain the huge variability and intermittency. Previous studies have claimed that fractal dimensions might be different from one region to an other, excluding a general description at the planetary scale. In this article, we are analyzing the Martian topographic data with a multifractal formalism to study the scaling intermittency. In the multifractal paradigm, the apparent local variation of the fractal dimension is interpreted as a statistical property of multifractal fields. We analyze the topography measured with the laser altimeter MOLA at 300 m horizontal resolution, 1 m vertical resolution. We adapted the Haar fluctuation method to the the irregularly sampled signal. The results suggest a multifractal behavior from planetary scale down to 10 km. From 10 km to 300 m, the topography seems to be simple monofractal. This transition indicates a significant change in the geological processes governing the Red Planet's surface.
Radar networks for automtovie short-range applications (up to 30m) based on powerful but inexpensive 24GHz high range resolution pulse or FMCW radar systems have been developed at the Technical University of Hamburg-Harburg. The described system has been integrated in to an experimental vehicle and tested in real street environment. This paper considers the general network design, the individual pulse or FMCW radar sensors, the network signal processing scheme, the tracking procedure and possible automotive applications, respectively. Object position estimation is accomplished by the very precise range measurement of each individual sensor and additional trilateration procedures. The paper concludes with some results obtained in realistic traffic conditions with multiple target situations using 24 GHz radar network.
The invention of radar, as mentioned in Chris Lavers' article on warship stealth technology (March pp21-25), continues to be a subject of discussion. Here in Malvern we have just unveiled a blue plaque to commemorate the physicist Albert Percival Rowe, who arrived in 1942 as the head of the Telecommunications Research Establishment (TRE), which was the Air Ministry research facility responsible for the first British radar systems.
Kirk, John C.; Lin, Kai; Gray, Andrew; Hseih, Chung; Darden, Scott; Kwong, Winston; Majumder, Uttam; Scarborough, Steven
A small and lightweight dual-channel radar has been developed for SAR data collections. Using standard Displaced Phase Center Antenna (DPCA) radar digital signal processing, SAR GMTI images have been obtained. The prototype radar weighs 5-lbs and has demonstrated the extraction of ground moving targets (GMTs) embedded in high-resolution SAR imagery data. Heretofore this type of capability has been reserved for much larger systems such as the JSTARS. Previously, small lightweight SARs featured only a single channel and only displayed SAR imagery. Now, with the advent of this new capability, SAR GMTI performance is now possible for small UAV class radars.
Duxbury, T. C.
Surface topography is an important constraint when investigating the evolution of solar system bodies. Topography is typically obtained from stereo photogrammetric or photometric (shape from shading) analyses of overlapping / stereo images and from laser / radar altimetry data. The ESA Mars Express Mission  carries a Super Resolution Channel (SRC) as part of the High Resolution Stereo Camera (HRSC) . The SRC can build up overlapping / stereo coverage of Mars, Phobos and Deimos by viewing the surfaces from different orbits. The derivation of high precision topography data from the SRC raw images is degraded because the camera is out of focus. The point spread function (PSF) is multi-peaked, covering tens of pixels. After registering and co-adding hundreds of star images, an accurate SRC PSF was reconstructed and is being used to restore the SRC images to near blur free quality. The restored images offer a factor of about 3 in improved geometric accuracy as well as identifying the smallest of features to significantly improve the stereo photogrammetric accuracy in producing digital elevation models. The difference between blurred and restored images provides a new derived image product that can provide improved feature recognition to increase spatial resolution and topographic accuracy of derived elevation models. Acknowledgements: This research was funded by the NASA Mars Express Participating Scientist Program.  Chicarro, et al., ESA SP 1291(2009)  Neukum, et al., ESA SP 1291 (2009). A raw SRC image (h4235.003) of a Martian crater within Gale crater (the MSL landing site) is shown in the upper left and the restored image is shown in the lower left. A raw image (h0715.004) of Phobos is shown in the upper right and the difference between the raw and restored images, a new derived image data product, is shown in the lower right. The lower images, resulting from an image restoration process, significantly improve feature recognition for improved derived
Crosby, Christopher; Nandigam, Viswanath; Baru, Chaitan; Arrowsmith, J. Ramon
High-resolution topography data acquired with lidar (light detection and ranging) technology are revolutionizing the way we study the Earth's surface and overlying vegetation. These data, collected from airborne, tripod, or mobile-mounted scanners have emerged as a fundamental tool for research on topics ranging from earthquake hazards to hillslope processes. Lidar data provide a digital representation of the earth's surface at a resolution sufficient to appropriately capture the processes that contribute to landscape evolution. The U.S. National Science Foundation-funded OpenTopography Facility (http://www.opentopography.org) is a web-based system designed to democratize access to earth science-oriented lidar topography data. OpenTopography provides free, online access to lidar data in a number of forms, including the raw point cloud and associated geospatial-processing tools for customized analysis. The point cloud data are co-located with on-demand processing tools to generate digital elevation models, and derived products and visualizations which allow users to quickly access data in a format appropriate for their scientific application. The OpenTopography system is built using a service-oriented architecture (SOA) that leverages cyberinfrastructure resources at the San Diego Supercomputer Center at the University of California San Diego to allow users, regardless of expertise level, to access these massive lidar datasets and derived products for use in research and teaching. OpenTopography hosts over 500 billion lidar returns covering 85,000 km2. These data are all in the public domain and are provided by a variety of partners under joint agreements and memoranda of understanding with OpenTopography. Partners include national facilities such as the NSF-funded National Center for Airborne Lidar Mapping (NCALM), as well as non-governmental organizations and local, state, and federal agencies. OpenTopography has become a hub for high-resolution topography
Holt, J. W.; Blankenship, D. D.; Peters, M. E.; Kempf, S. D.; Williams, B. J.
seasons. Post-processing of the positioning data yields accuracies of ~ 0.10 m for samples at ~ 15 m intervals. Precise positioning was accomplished through the use of two carrier-phase GPS receivers on the aircraft and two at McMurdo Station. Surface and shallow subsurface properties are being supplied by glacial geomorphologists conducting ground-based studies in Taylor and Beacon Valleys. Two techniques are being used in parallel to discriminate subsurface echoes from surface echoes due to surrounding topography. In the first method, surface returns are simulated using aircraft position data, the modeled radar antenna pattern, and surface topography from a digital elevation model (DEM) recently acquired by the USGS and NASA in the Dry Valleys with 2-meter postings. These predicted surface returns are then compared with the actual data to reveal side echoes. The second method identifies all echoes in the radar data and maps them into possible correlative surface features to the sides of the aircraft through range estimation. This uses the measured time delay of the echo and known surface topography. We map the echoes onto the DEM (and optical imagery) at the appropriate range in order to identify candidate surface return sources. The two methods should identify all echoes that are not from the subsurface. The comparison of different radar configurations and parallel tracks where they are available will also be utilized to identify the source of any ambiguous echoes.
Fu, L. L.; Rodriguez, E.
We propose to apply the technique of synthetic aperture radar interferometry to the measurement of ocean surface topography at spatial resolu tion approaching 1 km . The measurement w ill have wide ranging applications in oceanography , hydrology , and marine geophysics. The oceanographic and related societal applications are briefly discussed in the paper. To meet the requirements for oceanographic application s, the in strument must be flown in an orbit w ith proper samp ling of ocean tides.
Fu, Lee-Lueng; Rodriquez, Ernesto
We propose to apply the technique of synthetic aperture radar interferometry to the measurement of ocean surface topography at spatial resolution approaching 1 km. The measurement will have wide ranging applications in oceanography, hydrology, and marine geophysics. The oceanographic and related societal applications are briefly discussed in the paper. To meet the requirements for oceanographic applications, the instrument must be flown in an orbit with proper sampling of ocean tides.
Williams, David R.; Greeley, Ronald
The Tellus Regio area of Venus represents a subset of a narrow latitude band where Pioneer Venus Orbiter (PVO) altimetry data, line-of-sight (LOS) gravity data, and Venera 15/16 radar images have all been obtained with good resolution. Tellus Regio also has a wide variety of surface morphologic features, elevations ranging up to 2.5 km, and a relatively low LOS gravity anomaly. This area was therefore chosen in order to examine the theoretical stress distributions resulting from various models of compensation of the observed topography. These surface stress distributions are then compared with the surface morphology revealed in the Venera 15/16 radar images. Conclusions drawn from these comparisons will enable constraints to be put on various tectonic parameters relevant to Tellus Regio. The stress distribution is calculated as a function of the topography, the equipotential anomaly, and the assumed model parameters. The topography data is obtained from the PVO altimetry. The equipotential anomaly is estimated from the PVO LOS gravity data. The PVO LOS gravity represents the spacecraft accelerations due to mass anomalies within the planet. These accelerations are measured at various altitudes and angles to the local vertical and therefore do not lend themselves to a straightforward conversion. A minimum variance estimator of the LOS gravity data is calculated, taking into account the various spacecraft altitudes and LOS angles and using the measured PVO topography as an a priori constraint. This results in an estimated equivalent surface mass distribution, from which the equipotential anomaly is determined.
Bertani, Thiago de Castilho; Rossetti, Dilce de Fátima; Albuquerque, Paulo Cesar Gurgel
Reconstructing the evolution of large tropical fluvial systems over the geological time is challenging, particularly in areas such as the Amazonian lowlands where basic geological and geomorphological data are still scarce relatively to the large dimension of the region. In such areas, remote sensing data are useful for detecting ancient morphological features that may reveal past fluvial dynamics. In this study, we explored object-based image analysis (OBIA) in the Madeira-Purus interfluve, Southwestern Brazilian Amazonia, integrating geospatial data including Landsat satellite multispectral images, the digital elevation model (DEM) acquired during the Shuttle Radar Topography Mission (SRTM), and stream channels digitized from topographic maps. This approach provided the basis to categorize automatically classes with contrasting vegetation and/or topographic characteristics within the dense tropical forest over an extensive and relatively flat forested area. The main goal was to use these classes as a surrogate for the recognition of ancient geomorphic features consisting mainly of paleochannels that may help reconstructing fluvial history in space and time. Landsat optical images with stream vector were appropriate to classify open vegetation areas that grow over paleochannels, but failed to identify these objects when they were located over forested areas. However, the digital elevation model (DEM) derived from the Shuttle Radar Topography Mission (SRTM) was successful to detect these objects even in forested areas. Topographic survey undertaken in the field increased the classification reliability by demonstrating true terrain variations along transects measured across the paleochannels. Based on this technique, networks of dendritic paleochannels were mapped and related to ancient tributaries of the Madeira River that had their courses flowing opposite to main modern streams. This denotes a significant change in fluvial dynamics over time, most likely
Margot, J. L.; Campbell, D. B.; Jurgens, R. F.; Slade, M. A.
Earth-based radar interferometry  has been used to map the lunar polar regions and Tycho Crater at high spatial ( ~ 100 m) and height ( ~ 50 m) resolutions. Compared to existing topographic data sets, the radar observations offer digital elevation models with dense horizontal spacing and improved height resolution. Earth-based radars can also provide measurements of the largely unknown topography in the polar regions. Elevation data and radar imagery obtained with the Goldstone X-band system (lambda = 3.5 cm) are presented for the Tycho Crater area, with a spatial resolution of 200 m and a height resolution of 30 m. A careful comparison of the radar-derived topography with Clementine altimetry points  reveals a very good agreement between the two techniques. Rms deviations between the radar-derived heights and 87 Clementine points available over the 200 x 200 km scene are ~ 100 m. The digital elevation model allows detailed morphometry of the 85 km diameter crater: the floor of Tycho lies 3970 m below a 1738 km radius sphere, and the crater's central peak rises 2400 m above the floor. The average rim crest elevation is 730 m above the 1738 km datum, giving a mean rim to floor depth of 4700 m. The floor has two distinct units with the western section being higher in elevation by ~ 200 m. This dichotomy is consistent with an asymmetry in the crater shape which reveals that maximum wall slumping occured in the western and southwestern regions of the crater. Digital elevation models of the polar regions are being used to estimate the location of permanently shadowed areas which may harbor ice deposits . The range of illumination conditions over the lunar polar regions could be sampled by an imaging instrument in a polar orbit during a full terrestrial year. Alternatively, topographic maps obtained with Earth-based radar can be used to model the illumination conditions over the entire solar illumination cycle.  I. I. Shapiro et al. (1972). Science, 178, 939
Moller, D.; Hensley, S.; Sadowy, G.; Wu, X.; Carswell, J.; Fisher, C.; Michel, T.; Lou, Y.
In May 2009 a new radar technique for mapping ice surface topography was demonstrated in a Greenland campaign as part of the NASA International Polar Year (IPY) activities. This was achieved with the airborne Glacier and Ice Surface Topography Interferometer (GLISTIN-A): a 35.6 GHz single-pass interferometer. Although the technique of using radar interferometry for mapping terrain has been demonstrated before, this is the first such application at millimeter-wave frequencies. The proof-of-concept demonstration was achieved by interfacing Ka-band RF and antenna hardware with the Uninhabited Airborne Vehicle Synthetic Aperture Radar (UAVSAR). The GLISTIN-A was implemented as a custom installation of the NASA Dryden Flight Research Center Gulfstream III. Instrument performance indicates swath widths over the ice between 5-7km, with height precisions ranging from 30cm-3m at a posting of 3m x 3m. Processing challenges were encountered in achieving the accuracy requirements on several fronts including, aircraft motion sensitivity, multipath and systematic drifts. However, through a combination of processor optimization, a modified phase-screen and motion-compensation implementations were able to minimize the impact of these systematic error sources. We will present results from the IPY data collections including system performance evaluations and imagery. This includes a large area digital elevation model (DEM) collected over Jakobshavn glacier as an illustrative science data product. Further, by intercomparison with the NASA Wallops Airborne Topographic Mapper (ATM) and calibration targets we quantify the interferometric penetration bias of the Ka-band returns into the snow cover. Following the success of the IPY campaign, we are funded under the Earth Science Techonology Office (ESTO) Airborne Innovative Technology Transition (AITT) program to transition GLISTIN-A to a permanently-available pod-only system compatible with unpressurized operation. In addition
Willner, K.; Shi, X.; Oberst, J.
The global shape and the dynamic environment are fundamental properties of a body. Other properties such as volume, bulk density, and models for the dynamic environment can subsequently be computed based on such models. Stereo-photogrammetric methods were applied to derive a global digital terrain model (DTM) with 100 m/pixel resolution using High Resolution Stereo Camera images of the Mars Express mission and Viking Orbiter images. In a subsequent least-squares fit, coefficients of the spherical harmonic function to degree and order 45 are computed. The dynamic models for Phobos were derived from a polyhedron representation of the DTM. The DTM, spherical harmonic function model, and dynamic models, have been refined and represent Phobos' dynamic and geometric topography with much more detail when compared to Shi et al. (2012) and Willner et al. (2010) models, respectively. The volume of Phobos has been re-determined to be in the order of 5741 km3 with an uncertainty of only 0.6% of the total volume. This reduces the bulk density to 1.86±0.013 g/cm3 in comparison to previous results. Assuming a homogeneous mass distribution a forced libration amplitude for Phobos of 1.14° is computed that is in better agreement with observations by Willner et al. (2010) than previous estimates.
Sayegh, Samir; Jiang, Yanshui
Commercial OCT systems provide pachymetry measurements. Full corneal topographic information of anterior and posterior corneal surfaces for use in cataract surgery and refractive procedures is a desirable goal and would add to the usefulness of anterior and posterior segment evaluation. While substantial progress has been made towards obtaining "average" central corneal power (D Huang), power in different meridians and topography are still missing. This is usually reported to be due to eye movement. We analyze the role of centration, eye movements and develop a model that allows for the formulation of criteria for obtaining reliable topographic data within ¼ diopter.
Taylor, R. M.
The radar astronomy activities supported by the Deep Space Network during June, July, and August 1980 are reported. The planetary bodies observed were Venus, Mercury, and the asteroid Toro. Data were obtained at both S and X band, and the observations were considered successful.
This three-dimensional perspective of the remote Karakax Valley in the northern Tibetan Plateau of western China was created by combining two spaceborne radar images using a technique known as interferometry. Visualizations like this are helpful to scientists because they reveal where the slopes of the valley are cut by erosion, as well as the accumulations of gravel deposits at the base of the mountains. These gravel deposits, called alluvial fans, are a common landform in desert regions that scientists are mapping in order to learn more about Earth's past climate changes. Higher up the valley side is a clear break in the slope, running straight, just below the ridge line. This is the trace of the Altyn Tagh fault, which is much longer than California's San Andreas fault. Geophysicists are studying this fault for clues it may be able to give them about large faults. Elevations range from 4000 m (13,100 ft) in the valley to over 6000 m (19,700 ft) at the peaks of the glaciated Kun Lun mountains running from the front right towards the back. Scale varies in this perspective view, but the area is about 20 km (12 miles) wide in the middle of the image, and there is no vertical exaggeration. The two radar images were acquired on separate days during the second flight of the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) aboard the space shuttle Endeavour in October 1994. The interferometry technique provides elevation measurements of all points in the scene. The resulting digital topographic map was used to create this view, looking northwest from high over the valley. Variations in the colors can be related to gravel, sand and rock outcrops. This image is centered at 36.1 degrees north latitude, 79.2 degrees east longitude. Radar image data are draped over the topography to provide the color with the following assignments: Red is L-band vertically transmitted, vertically received; green is the average of L-band vertically transmitted