A weighted least squares solution for space intersection of spaceborne stereo SAR data

The use of stereoscopic SAR images offers an alternative to interferometric SAR for the generation of digital elevation models (DEMs). The stereo radargrammetric method is robust and can generate DEMs of sufficient accuracy to geocode SAR images. Previous work has shown that ground coordinates with accuracy of four times the resolution cell can be obtained from ERS data without using any ground control points (GCPs), where the high accuracy of the orbit and satellite position of the order of metres introduce insignificant errors into the intersection procedure. The orbit data for RADARSAT is not as accurate as that for ERS, and the perpendicular relationship between the resultant velocity vector and the resultant range vector is uncertain in terms of image geometry. Hence, it is necessary to refine the method to allow for possible errors. This paper introduces a weighted space intersection algorithm based on an analysis of the predicted errors. A radargrammetric error model for observation errors is also formulated to predict the accuracy of the algorithm. The revised method can be used without any GCPs, but this can lead to systematic errors due to less accurate orbit data, and it has been found that the use of two GCPs provides a reasonable solution. The method is insensitive to the spatial distribution of GCPs, which is often critical in traditional methods. The error statistics of the results generated from 32 independent check points, distributed through the entire SAR image, approach the predicted errors and give positional accuracy of 38 m in three dimensions

Matching stereoscopic SAR images for radargrammetric applications

International audienceThe aim of this paper is to present our studies about extraction of 3D information from radar images. Several radargrammetric methods allow DEM (Digital Elevation Model) generation from SAR images and we take a special interest to stereoscopic method. The main idea is to match image stereo pairs, to create a disparity map from one image to the other and to compute elevation thanks to the incidences angle

Height Accuracy of Radargrammetric Generated Digital Elevation Model

The stereoscopic concept used in photogrammetry was successfully extended to the optical sensor of SPOT (System Probatoire d'Observation de la Terre) to generate the digital elevation model (DEM) and to extract the planimetric features. Following this success, scientists are impressively encouraged to explore the feasibility of using stereoscopy on synthetic aperture radar (SAR) images, applying the radargrammetric technique. In the regions where cloud cover or darkness prevails, active microwave remote sensing data such as SAR can be fully utilized to procure information about land surface and forest canopy. The objective of this study is to analyze the accuracy of the elevations extracted from the DEM generated using the radargrammetric technique as compared with the elevations generated from the photogrammetric technique. The capability of the radargrammetric technique and its potential in extracting the altimetric information were subsequently assessed. The stereo RADARSAT images of Klang Valley with coverage of 100 Ian by 100 Ian were acquired. A total of 199 Ground Control Points (GCPs) were selected on relatively low terrains as backscattering radar data is very sensitive to the slope and high terrains. The primary input data was the coordinates of GCPs which were extracted from the topographical maps. When the errors produced in the GCP collection report were acceptable, the next process was creating epipolar images and generating DEM. It was then followed by generating the geocoded DEM

Extraction of spatial information from sterioscopic SAR images

Synthetic Aperture Radar (SAR) is now widely used for generating Digital Elevation Models (DEMs) and has advantages over optical data in terms of availability as it allows all-day and all-weather operations. The stereoscopic SAR method, which allows direct extraction of spatial information in three-dimensional space, has been established for decades. However, the traditional stereoscopic methods developed for SAR data depend on many human operations and need ground control points (GCPs), to set up geometric models. The aims of the thesis are not only to propose a refined rigorous stereoscopic SAR method and a new error model to predict theoretic errors, but also to achieve a higher level of automation and accuracy. By using a weighting matrix, which is derived by considering different observations in the space intersection algorithm, the minimal number of the GCPs required for the refined algorithm is only two. To achieve a high degree of automation, an optimized strategy of parameter selection for the pyramidal image correlation scheme employing a region-growing technique has been proposed. This avoids a trial-and-error approach to produce digital parallax data from the same-side SAR image pairs. A new method to derive GCPs automatically has been developed using a SAR image simulation technique, under the condition that a known DEM chip is available, to minimize human interventions and operator error. The proposed method for providing GCPs and the DEMs generated from space intersection have been incorporated into the procedures for geocoding SAR images to validate the proposed algorithms. The results derived show that the stereoscopic SAR data can be applied to geometric rectification in flat-to-moderate areas, and other applications of extraction of spatial information are promising

Effects of Orbit and Pointing Geometry of a Spaceborne Formation for Monostatic-Bistatic Radargrammetry on Terrain Elevation Measurement Accuracy

During the last decade a methodology for the reconstruction of surface relief by Synthetic Aperture Radar (SAR) measurements – SAR interferometry – has become a standard. Different techniques developed before, such as stereo-radargrammetry, have been experienced from space only in very limiting geometries and time series, and, hence, branded as less accurate. However, novel formation flying configurations achievable by modern spacecraft allow fulfillment of SAR missions able to produce pairs of monostatic-bistatic images gathered simultaneously, with programmed looking angles. Hence it is possible to achieve large antenna separations, adequate for exploiting to the utmost the stereoscopic effect, and to make negligible time decorrelation, a strong liming factor for repeat-pass stereo-radargrammetric techniques. This paper reports on design of a monostatic-bistatic mission, in terms of orbit and pointing geometry, and taking into account present generation SAR and technology for accurate relative navigation. Performances of different methods for monostatic-bistatic stereo-radargrammetry are then evaluated, showing the possibility to determine the local surface relief with a metric accuracy over a wide range of Earth latitudes

High resolution radargrammetry with COSMO-SkyMed, TerraSAR-X and RADARSAT-2 imagery: development and implementation of an image orientation model for Digital Surface Model generation

Digital Surface and Terrain Models (DSM/DTM) have large relevance in several territorial applications, such as topographic mapping, monitoring engineering, geology, security, land planning and management of Earth's resources. The satellite remote sensing data offer the opportunity to have continuous observation of Earth's surface for territorial application, with short acquisition and revisit times. Meeting these requirements, the SAR (Synthetic Aperture Radar) high resolution satellite imagery could offer night-and-day and all-weather functionality (clouds, haze and rain penetration). Two different methods may be used in order to generate DSMs from SAR data: the interferometric and the radargrammetric approaches. The radargrammetry uses only the intensity information of the SAR images and reconstructs the 3D information starting from a couple of images similarly to photogrammetry. Radargrammetric DSM extraction procedure consists of two basic steps: the stereo pair orientation and the image matching for the automatic detection of homologous points. The goal of this work is the definition and the implementation of a geometric model in order to orientate SAR imagery in zero Doppler geometry. The radargrammetric model implemented in SISAR (Software per Immagini Satellitari ad Alta Risoluzione - developed at the Geodesy and Geomatic Division - University of Rome "La Sapienza") is based on the equation of radar target acquisition and zero Doppler focalization Moreover a tool for the SAR Rational Polynomial Coefficients (RPCs) generation has been implemented in SISAR software, similarly to the one already developed for the optical sensors. The possibility to generate SAR RPCs starting from a radargrammetric model sounds of particular interest since, at present, the most part of SAR imagery is not supplied with RPCs, although the RPFs model is available in several commercial software. Only RADARSAT-2 data are supplied with vendors RPCs. To test the effectiveness of the implemented RPCs generation tool and the SISAR radargrammetric orientation model the reference results were computed: the stereo pairs were orientated with the two model. The tests were carried out on several test site using COSMO-SkyMed, TerraSAR-X and RADARSAT-2 data. Moreover, to evaluate the advantages and the different accuracy between the orientation models computed without GCPs and the orientation model with GCPs a Monte Carlo test was computed. At last, to define the real effectiveness of radargrammetric technique for DSM extraction and to compare the radrgrammetric tool implemented in a commercial software PCI-Geomatica v. 2012 and SISAR software, the images acquired on Beauport test site were used for DSM extraction. It is important underline that several test were computed. Part of this tests were carried out under the supervision of Prof. Thierry Toutin at CCRS (Canada Centre of Remote Sensing) where the PCI-Geomatica orientation model was developed, in order to check the better parameters solution to extract radargrammetric DSMs. In conclusion, the results obtained are representative of the geometric potentialities of SAR stereo pairs as regards 3D surface reconstruction

Disaster debris estimation using high-resolution polarimetric stereo-SAR

AbstractThis paper addresses the problem of debris estimation which is one of the most important initial challenges in the wake of a disaster like the Great East Japan Earthquake and Tsunami. Reasonable estimates of the debris have to be made available to decision makers as quickly as possible. Current approaches to obtain this information are far from being optimal as they usually rely on manual interpretation of optical imagery. We have developed a novel approach for the estimation of tsunami debris pile heights and volumes for improved emergency response. The method is based on a stereo-synthetic aperture radar (stereo-SAR) approach for very high-resolution polarimetric SAR. An advanced gradient-based optical-flow estimation technique is applied for optimal image coregistration of the low-coherence non-interferometric data resulting from the illumination from opposite directions and in different polarizations. By applying model based decomposition of the coherency matrix, only the odd bounce scattering contributions are used to optimize echo time computation. The method exclusively considers the relative height differences from the top of the piles to their base to achieve a very fine resolution in height estimation. To define the base, a reference point on non-debris-covered ground surface is located adjacent to the debris pile targets by exploiting the polarimetric scattering information. The proposed technique is validated using in situ data of real tsunami debris taken on a temporary debris management site in the tsunami affected area near Sendai city, Japan. The estimated height error is smaller than 0.6m RMSE. The good quality of derived pile heights allows for a voxel-based estimation of debris volumes with a RMSE of 1099m3. Advantages of the proposed method are fast computation time, and robust height and volume estimation of debris piles without the need for pre-event data or auxiliary information like DEM, topographic maps or GCPs

Elevation modelling and palaeo-environmental interpretation in the Siwa area (Egypt): Application of SAR interferometry and radargrammetry to COSMO-SkyMed imagery

PublishedJournal Article© 2015 Elsevier B.V. Digital elevation models produced from COSMO-SkyMed imagery were used to delineate palaeo-drainage in a wide area surrounding the Siwa and Al-Jaghbub oases of the western Sahara Desert (Egypt and Libya).This new generation of synthetic aperture radar imagery is suitable for this purpose because of its high spatial resolution and capacity to penetrate dry surface sediments. Different techniques such as radar interferometry and radargrammetry were used to produce digital elevation models. These were assessed for accuracy and then combined to produce a single elevation model of the area. The resulting elevation model was used to support the geological study and palaeo-environmental interpretation of the area. It revealed buried features of the landscape, including inactive palaeo-drainage systems. Drainage features were extracted from the elevation model using geographical information systems; results were combined and assessed with respect to geological field data, as well as data from the literature. Previous studies in the area suggest that a wide river, probably the old Nile River, flowed into the Libyan palaeo-Sirte before the Late Messinian drawdown of the Mediterranean Sea. During the Late Messinian lowering of the sea the fluvial system changed shape and carved deep canyons throughout north-eastern Africa. The reported findings on the key Siwa area were used to precisely delineate the physiography of the modern drainage network and to confirm findings from our previous geological research in the area.We gratefully acknowledge Ruggero Matteucci and Johannes Pignatti (La Sapienza, University of Rome), Francesco Checchi (ENI S.p.A., IOEC), Filippo Bonciani and Debora Graziosi (University of Siena) for their collaboration. Research was supported by the ASI (Id 2262) in the framework of the COSMO-SkyMed Announcement of Opportunity project “Application of COSMO-SkyMed data for geological researches in Egypt and Libya”

Interferometric Synthetic Aperture RADAR and Radargrammetry towards the Categorization of Building Changes

The purpose of this work is the investigation of SAR techniques relying on multi image acquisition for fully automatic and rapid change detection analysis at building level. In particular, the benefits and limitations of a complementary use of two specific SAR techniques, InSAR and radargrammetry, in an emergency context are examined in term of quickness, globality and accuracy. The analysis is performed using spaceborne SAR data