Space-based and airplane-based synthetic aperture RADAR (SAR) can monitor ground height using interferometric SAR (InSAR) collections. However, fielding the airplane-based SAR is expensive and coordinating the frequency and timing of ground experiments with space-based SAR is challenging. This research explored the possibility of using a small, mobile unmanned aerial vehicle- base (UAV) SAR to see if it could provide a quick and inexpensive InSAR option for the Source Physics Experiment (SPE) Phase III project. Firstly, a local feasibility collection using a UAV-based SAR showed that InSAR products and height measurements were possible, but that in-scene fiducials were needed to assist in digital elevation model (DEM) construction. Secondly, an InSAR collection was planned and executed over the SPE Phase III site using the same platform configuration. We found that the image formation by the SAR manufacturer creates discontinuities, and that noise impacted the generation and accuracy of height maps. These processing artifacts need to be overcome to generate an accurate height map.
Seismic signals from 1993 show a series of magnitude (Mw) 3.7 or less seismic events in Rock Valley on the Nevada National Security Site (NNSS). Historic synthetic aperture radar images of that location were found that could provide interferometric synthetic aperture radar (InSAR) measures of the ground height during 1993. Given this historic SAR imagery, we explore answering the question if ground movement from the 1993 Rock Valley earthquake activity could be sensed by remote sensing means. Finding earthquake surface movement would assist in locating the Rock Valley fault and the 1993 earthquake hypocenter where the Source Physics Experiment Phase III series of experiments will be conducted. In this report, we show that InSAR can sense very small height differences, and for the European Radar Satellite-1 InSAR collections during 1992 and 1993 over Rock Valley earth surface movements were measured with 8 mm uplift and 12.5 mm subsidence over isolated areas. One of these earth movement areas coincides with an InSAR image pair coherence drop between March 5, 1993 and June 18, 1993. The coherence drop is over an approximately 13 square km area south southeast of Skull Mountain centered at 36° 43' 30" N latitude and 116° 05' 00" W longitude. Measured small surface movement and a loss of InSAR coherence may be caused by the series of earthquakes. The location of these InSAR detections may also coincide with water drainage or erosion displacement. There are no records to disambiguate the earthquake and erosion earth surface motion possibilities. Therefore, the InSAR findings of earth surface movement by InSAR are inconclusive.
Performing terrain classification with data from heterogeneous imaging modalities is a very challenging problem. The challenge is further compounded by very high spatial resolution. (In this paper we consider very high spatial resolution to be much less than a meter.) At very high resolution many additional complications arise, such as geometric differences in imaging modalities and heightened pixel-by-pixel variability due to inhomogeneity within terrain classes. In this paper we consider the fusion of very high resolution hyperspectral imaging (HSI) and polarimetric synthetic aperture radar (PolSAR) data. We introduce a framework that utilizes the probabilistic feature fusion (PFF) one-class classifier for data fusion and demonstrate the effect of making pixelwise, superpixel, and pixelwise voting (within a superpixel) terrain classification decisions. We show that fusing imaging modality data sets, combined with pixelwise voting within the spatial extent of superpixels, gives a robust terrain classification framework that gives a good balance between quantitative and qualitative results.
The Source Physics Experiment (SPE) Phase I conducted six underground chemical explosions at the same experimental pad with the goal of characterizing underground explosions to enhance the United States (U.S.) ability to detect and discriminate underground nuclear explosions (UNEs). A fully polarimetric synthetic aperture RADAR (PolSAR) collected imagery in VideoSAR mode during the fifth and sixth explosions in the series (SPE-5 and SPE-6). Previously, we reported the prompt PolSAR surface changes cause by SPE-5 and SPE-6 explosions within seconds or minutes of the underground chemical explosions, including a drop of spatial coherence and polarimetric scattering changes. Therein it was hypothesized that surface changes occurred when surface particles experienced upward acceleration greater than 1 g. Because the SPE site was instrumented with surface accelerometers, we explore that hypothesis and report our findings in this article. We equate explosion-caused prompt surface expressions measured by PolSAR to the prompt surface movement measured by accelerometers. We tie these findings to UNE detection by comparing the PolSAR and accelerometer results to empirical ground motion predictions derived from accelerometer recordings of UNEs collected prior to cessation of U.S. nuclear testing. We find the single threshold greater than 1 g hypothesis is not correct for it does not explain the PolSAR results. Our findings show PolSAR surface coherence spatial extent is highly correlated with surface velocity, both measured and predicted, and the resulting surface deformation extent is corroborated by accelerometer records and the predicted lateral spall extent. PolSAR scattering changes measured during SPE-6 are created by the prompt surface displacement being larger than the spall gap.
Phase I of the Source Physics Experiment (SPE) series involved six underground chemical explosions, all of which were conducted at the same experimental pad. Research from the sixth explosion of the series (SPE-6) demonstrated that polarimetric synthetic aperture radar (PolSAR) is a viable technology for monitoring an underground chemical explosion when the geologic structure is Cretaceous granitic intrusive. It was shown that a durable signal is measurable by the H/A/alpha polarimetric decomposition parameters. After the SPE-6 experiment, the SPE program moved to the Phase II location, which is composed of dry alluvium geology (DAG). The loss of wavefront energy is greater through dry alluvium than through granite. In this article, we compare the SPE-6 analysis to the second DAG (DAG-2) experiment. We hypothesize that despite the geology at the DAG site being more challenging than at the Phase I location, combined with the DAG-2 experiment having a 3.37 times deeper scaled depth of burial than the SPE-6, a durable nonprompt signal is still measurable by a PolSAR sensor. We compare the PolSAR time-series measures from videoSAR frames, from the SPE-6 and DAG-2 experiments, with accelerometer data. We show which PolSAR measures are invariant to the two types of geology and which are geology dependent. We compare a coherent change detection (CCD) map from the DAG-2 experiment with the data from a fiber-optic distributed acoustic sensor to show the connection between the spatial extent of coherence loss in CCD maps and spallation caused by the explosion. Finally, we also analyze the spatial extent of the PolSAR measures from both explosions.
Deciding on an imaging modality for terrain classification can be a challenging problem. For some terrain classes a given sensing modality may discriminate well, but may not have the same performance on other classes that a different sensor may be able to easily separate. The most effective terrain classification will utilize the abilities of multiple sensing modalities. The challenge of utilizing multiple sensing modalities is then determining how to combine the information in a meaningful and useful way. In this paper, we introduce a framework for effectively combining data from optical and polarimetric synthetic aperture radar sensing modalities. We demonstrate the fusion framework for two vegetation classes and two ground classes and show that fusing data from both imaging modalities has the potential to improve terrain classification from either modality, alone.
There are several factors that should be considered for robust terrain classification. We address the issue of high pixel-wise variability within terrain classes from remote sensing modalities, when the spatial resolution is less than one meter. Our proposed method segments an image into superpixels, makes terrain classification decisions on the pixels within each superpixel using the probabilistic feature fusion (PFF) classifier, then makes a superpixel-level terrain classification decision by the majority vote of the pixels within the superpixel. We show that this method leads to improved terrain classification decisions. We demonstrate our method on optical, hyperspectral, and polarimetric synthetic aperture radar data.
Sandia National Laboratories flew its Facility for Advanced RF and Algorithm Development X-Band (9.6-GHz center frequency), fully polarimetric synthetic aperture radar (PolSAR) in VideoSAR mode to collect complex-valued SAR imagery before, during, and after the sixth Source Physics Experiment's (SPE-6) underground explosion. The VideoSAR products generated from the data sets include 'movies' of single-and quad-polarization coherence maps, magnitude imagery, and polarimetric decompositions. Residual defocus, due to platform motion during data acquisition, was corrected with a digital elevation model-based autofocus algorithm. We generated and exploited the VideoSAR image products to characterize the surface movement effects caused by the underground explosion. Unlike seismic sensors, which measure local area seismic waves using sparse spacing and subterranean positioning, these VideoSAR products captured high-spatial resolution, 2-D, time-varying surface movement. The results from the fifth SPE (SPE-5) used single-polarimetric VideoSAR data. In this paper, we present single-polarimetric and fully polarimetric VideoSAR results while monitoring the SPE-6 underground chemical explosion. We show that fully polarimetric VideoSAR imaging provides a unique, coherent, time-varying measure of the surface expression of the SPE-6 underground chemical explosion. We include new surface characterization results from the measured PolSAR SPE-6 data via H/A/α polarimetric decomposition.
Sandia National Laboratories (SNL) flew its Facility for Advanced RF and Algorithm Development (FARAD) X-Band (9.6 GHz center frequency), fully-polarimetric synthetic aperture radar (PolSAR) in VideoSAR-mode to collect complex-valued SAR imagery before, during, and after the fifth and sixth Source Physics Experiment's (SPE-5 and SPE-6) underground explosion. The results from the fifth Source Physics Experiment (SPE-5) used single-polarimetric VideoSAR data while SPE-6 used single and fully-polarimetric VideoSAR data. We show that SAR can provide surface change products indicative of disturbances caused by the underground chemical explosions. These are surface coherence measures, Po1SAR change signatures, and differential interferometric SAR (InSAR) height change.
Fully-polarimetric X-band (9.6 GHz center frequency) VideoSAR with 0.125-meter ground resolution flew collections before, during, and after the fifth Source Physics Experiment (SPE-5) underground chemical explosion. We generate and exploit synthetic aperture RADAR (SAR) and VideoSAR products to characterize surface effects caused by the underground explosion. To our knowledge, this has never been done. Exploited VideoSAR products are "movies" of coherence maps, phase-difference maps, and magnitude imagery. These movies show two-dimensional, time-varying surface movement. However, objects located on the SPE pad created unwanted, vibrating signatures during the event which made registration and coherent processing more difficult. Nevertheless, there is evidence that dynamic changes are captured by VideoSAR during the event. VideoSAR provides a unique, coherent, time-varying measure of surface expression of an underground chemical explosion.
In past research, two-pass repeat-geometry synthetic aperture radar (SAR) coherent change detection (CCD) predominantly utilized the sample degree of coherence as a measure of the temporal change occurring between two complex-valued image collects. Previous coherence-based CCD approaches tend to show temporal change when there is none in areas of the image that have a low clutter-to-noise power ratio. Instead of employing the sample coherence magnitude as a change metric, in this paper, we derive a new maximum-likelihood (ML) temporal change estimate-the complex reflectance change detection (CRCD) metric to be used for SAR coherent temporal change detection. The new CRCD estimator is a surprisingly simple expression, easy to implement, and optimal in the ML sense. This new estimate produces improved results in the coherent pair collects that we have tested.
In previous research, two-pass repeat-geometry synthetic aperture radar (SAR) coherent change detection (CCD) predominantly utilized the sample degree of coherence as a measure of the temporal change occurring between two complex-valued image collects. Previous coherence-based CCD approaches tend to show temporal change when there is none in areas of the image that have a low clutter-to-noise power ratio. Instead of employing the sample coherence magnitude as a change metric, in this paper, we derive a new maximum-likelihood (ML) temporal change estimate—the complex reflectance change detection (CRCD) metric to be used for SAR coherent temporal change detection. The new CRCD estimator is a surprisingly simple expression, easy to implement, and optimal in the ML sense. As a result, this new estimate produces improved results in the coherent pair collects that we have tested.
Typical synthetic aperture RADAR (SAR) imaging employs a co-located RADAR transmitter and receiver. Bistatic SAR imaging separates the transmitter and receiver locations. A bistatic SAR configuration allows for the transmitter and receiver(s) to be in a variety of geometric alignments. Sandia National Laboratories (SNL) / New Mexico proposed the deployment of a ground-based RADAR receiver. This RADAR receiver was coupled with the capability of digitizing and recording the signal collected. SNL proposed the possibility of creating an image of targets the illuminating SAR observes. This document describes the developed hardware, software, bistatic SAR configuration, and its deployment to test the concept of a ground-based bistatic SAR. In the proof-of-concept experiments herein, the RADAR transmitter will be a commercial SAR satellite and the RADAR receiver will be deployed at ground level, observing and capturing RADAR ground/targets illuminated by the satellite system.
While typical SAR imaging employs a co-located (monostatic) RADAR transmitter and receiver, bistatic SAR imaging separates the transmitter and receiver locations. The transmitter and receiver geometry determines if the scattered signal is back scatter, forward scatter, or side scatter. The monostatic SAR image is backscatter. Therefore, depending on the transmitter/receiver collection geometry, the captured imagery may be quite different that that sensed at the monostatic SAR. This document presents imagery and image products formed from captured signals during the validation stage of the bistatic SAR research. Image quality and image characteristics are discussed first. Then image products such as two-color multi-view (2CMV) and coherent change detection (CCD) are presented.
This report describes the significant processing steps that were used to take the raw recorded digitized signals from the bistatic synthetic aperture RADAR (SAR) hardware built for the NCNS Bistatic SAR project to a final bistatic SAR image. In general, the process steps herein are applicable to bistatic SAR signals that include the direct-path signal and the reflected signal. The steps include preprocessing steps, data extraction to for a phase history, and finally, image format. Various plots and values will be shown at most steps to illustrate the processing for a bistatic COSMO SkyMed collection gathered on June 10, 2013 on Kirtland Air Force Base, New Mexico.
In this paper, we derive a new optimal change metric to be used in synthetic aperture RADAR (SAR) coherent change detection (CCD). Previous CCD methods tend to produce false alarm states (showing change when there is none) in areas of the image that have a low clutter-to-noise power ratio (CNR). The new estimator does not suffer from this shortcoming. It is a surprisingly simple expression, easy to implement, and is optimal in the maximum-likelihood (ML) sense. The estimator produces very impressive results on the CCD collects that we have tested.
In this paper we show that the technique for spotlight-mode SAR image formation generally known as "backprojection" or "time- domain" is most easily derived and described in terms of the well-known methods of phased-array beamforming. By contrast, backprojection has been typically developed via analogy to tomographic imaging [1], which restricts this technique to the case of planar wavefronts. We demonstrate how the very simple notion of delay-and-sum beamforming leads directly to the backprojection algorithm for SAR, including the case for curved wavefronts. We further explain why backprojection offers a certain elegant simplicity for SAR imaging, and allows direct one-step computation of several useful SAR products, including an orthographically correct image free of any geometric or defocus effects from wavefront curvature and also free of the effects of terrain-elevation-induced defocus. (This product requires as an input a pre-existing digital elevation map (DEM) of the scene to be imaged.) In addition, we'll demonstrate why beamforming yields a mode-independent SAR image formation algorithm, i.e. one that can just as easily accommodate strip-map or spotlight-mode phase histories collected on an arbitrary flight path.
Coherent stereo pairs from cross-track synthetic aperture radar (SAR) collects allow fully automated correlation matching using magnitude and phase data. Yet, automated feature matching (correspondence) becomes more difficult when imaging rugged terrain utilizing large stereo crossing angle geometries because high-relief features can undergo significant spatial distortions. These distortions sometimes cause traditional, shift-only correlation matching to fail. This paper presents a possible solution addressing this difficulty. Changing the complex correlation maximization search from shift-only to shift-and-scaling using the downhill simplex method results in higher correlation. This is shown on eight coherent spotlight-mode cross-track stereo pairs with stereo crossing angles averaging 93.7° collected over terrain with slopes greater than 20°. The resulting digital elevation maps (DEMs) are compared to ground truth. Using the shift-scaling correlation approach to calculate disparity, height errors decrease and the number of reliable DEM posts increase.
Useful products generated from interferometric synthetic aperture radar (IFSAR) complex data include height measurement, coherent change detection, and classification. The IFSAR coherence is a spatial measure of complex correlation between two collects, a product of IFSAR signal processing. A tacit assumption in such IFSAR signal processing is that the terrain height is constant across an averaging box used in the process of correlating the two images. This paper presents simulations of IFSAR coherence if two targets with different heights exist in a given correlation cell, a condition in IFSAR collections produced by layover. It also includes airborne IFSAR data confirming the simulation results. The paper concludes by exploring the implications of the results on IFSAR height measurements and classification.