Radar DEM generation; achievements and benefits


 HOME

  Radar DEM generation; achievements and benefits 

Parviz Tarikhi

Head, Microwave Remote Sensing Research Core

Mahdasht Satellite Receiving Station, ISA

parviz_tarikhi@hotmail.com

  

[Prepared to be presented to the ISNET (Inter-Islamic Network on Space Science and Technology)/GORS (Syrian General Organization of Remote Sensing) Workshop on SAR and Optical Data for Natural Resources and Environmental Management, Damascus, Syria, 2-5 November 2009 (I failed to attend this workshop.)]

  

Abstract

InSAR is a radar technique for combining synthetic aperture radar (SAR) single look complex images to form interferogram and utilizing its phase contribution to land topography, surface movement and target velocity. In recent years considerable applications of Interferometric SAR technique are developed. It is an established technique for generating high quality digital elevation models (DEM) from spaceborne and airborne data, and that it has advantages over other methods for the generation of large area DEM. InSAR is capable of producing DEMs with the precision of a couple of ten meters whereas its movement map results have sub-centimeter precision.

DEMs are used in many applications in the context of earth sciences such as in topographic mapping, environmental modelling, rainfall-runoff studies, landslide hazard zonation, and seismic source modelling.

In this paper the experience gained on the subject in course of the continuous research work since 1994 in the former Iranian Remote Sensing Center and the Iranian Space Agency is given and the case studies on Izmit quake of August 1999 as well as the Bam quake of December 2003 using the related single look complex (SLC) images collected by the European ERS-1&2 and Envisat remote sensing satellites is discussed. Using tandem images of ERS-1&2 with the temporal baseline of one day leads to very good results for DEM generation while the special baselines smaller than 40m may not be suitable to this mean since slopes will be under-sampled.   

________________________________________

Introduction

Digital Elevation Model (DEM) refers to the process of demonstrating terrain elevation characteristics in 3-D space, but very often it specifically means the raster or regular grid of spot heights. DEM is the simplest form of digital representation of topography, while Digital Surface Model (DSM) describes the visible surface of the Earth (Fig. 1).

Fig. 1: DEM concept

Radar is a useful tool for land and planetary surface mapping. It is a good mean for obtaining a general idea of the geological setting of the area before proceeding for field work. Time, incidence angle, resolutions and coverage area all play important role at the outcome.

InSAR is a sophisticated processing of radar data for combining synthetic aperture radar (SAR) single look complex (SLC) images to form interferogram and utilizing its phase contribution to generate DEM, surface deformation and movement maps and target velocity. The interferogram contains phase difference of two images to which the imaging geometry, topography, surface displacement, atmospheric change and noise are the contributing factors.

Considerable applications of InSAR have been developed leaving it an established technique for high-quality DEM generation from spaceborne and airborne data and that it has advantages over other methods for the large-area DEM generation. It is capable of producing DEMs with the precision of a couple of ten meters while its movement map results have sub-centimeter precision over time spans of days to years. Terrestrial use of InSAR for DEM generation was first reported in 1974. It is used for different means particularly in geo-hazards like earthquakes, volcanoes, landslides and land subsidence.

 

Data sources and software

Satellite-based InSAR began in the 1980s using Seasat data, although the technique’s potential was expanded in the 1990s with launch of ERS-1 (1991), JERS-1 (1992), Radarsat-1 and ERS-2 (1995). They provided the stable well-defined orbits and short baselines necessary for InSAR. The 11-day NASA STS-99 mission in February 2000 used two SAR antenna with 60-m separation to collect data for the Shuttle Radar Topography Mission (SRTM). As a successor to ERS, in 2002 ESA launched the Advanced SAR (ASAR) aboard Envisat. Majority of InSAR systems has utilized the C-band sensors, but recent missions like ALOS PALSAR and TerraSAR-X are using L- and X-band. ERS and Radarsat use the frequency of 5.375GHz for instance. Numerous InSAR processing packages are also used commonly. IMAGINE-InSAR, EarthView-InSAR, ROI-PAC, DORIS, SAR-e2, Gamma, SARscape, Pulsar, IDIOT and DIAPASON are common for interferometry and DEM generation.

 

DEM generation steps using InSAR

The procedure for producing interferograms involves applied software and needed precision, but it usually includes some basic steps. InSAR operational steps for DEM generation from SLC-SAR data are performed in three major stages including data search, processing and product validation, while data processing stage comprises of five steps of data pre-processing, co-registration, interferogram generation, phase unwrapping, and geo-coding.

Two coherent SAR images are required to produce an interferogram. The images are first co-registered for finding the offset and difference in geometry between two amplitude images. Normally the baselines of 80-300m are suitable for DEM generation. One SAR image is then re-sampled to match the geometry of the other, meaning each pixel represents the same ground area in both images. The interferogram is generated by multiplication of the first image (master) to the complex conjugate of the second image (slave), and the interferometric phase due to the reference ellipsoid is removed; the process is referred to as flattening.

When produced, the interferogram is filtered using an adaptive power-spectrum filter to amplify the phase signal. For most quantitative applications the consecutive fringes present in the interferogram have to be unwrapped, which involves interpolating over the 0-2π phase jumps to produce a continuous deformation field. At some point, before or after unwrapping, incoherent areas of the image may be masked out. The final processing stage requires geo-coding, which involves resampling the interferogram from the acquisition geometry related to direction of satellite path into the desired geographic projection. By applying phase information in the interferogram, extraction of DEM with meter accuracy, height change information, and fine scale temporal change measurements will be possible.

 

DEM generation methods

DEM is important for surveying and other applications in engineering. Its accuracy is paramount; for some applications high accuracy does not matter but for some others it does. Numerous DEM generation techniques with different accuracies for various means are used. DEMs can be generated through different methods which are classified in three groups that are DEM generation by (i) geodesic measurements, (ii) photogrammetry and (iii) remote sensing.

In DEM generation by geodesic measurements, the planimetric coordinates and height values of each point of the feature are summed point-by-point and using the acquired data the topographic maps are generated with contour lines. The 1:25000-scale topographic maps are common example. The method uses contour-grid transfer to turn the vector data from the maps into digital data. For DEM generation by photogrammetry, the photographs are taken from an aircraft or spacecraft and evaluated as stereo-pairs and consequently 3-D height information is obtained.

DEM generation by remote sensing can be made in some ways, including stereo-pairs, laser scanning (LIDAR) and InSAR. There are three types of InSAR technique that is single-pass, double-pass and three-pass (Fig.2). In double-pass InSAR, a single SAR instrument passes over the same area two times while through the differences between these observations, height can be extracted. In three-pass interferometry (or DInSAR) the obtained interferogram of a double-pass InSAR for the commonly tandem image pairs is subtracted from the third image with wider temporal baseline respective to the two other images. In single-pass InSAR, space-craft has two SAR instrument aboard which acquire data for same area from different view angles at the same time. With single-pass, third dimension can be extracted and the phase difference between the first and second radar imaging instruments give the height value of the point of interest with some mathematical method. SRTM used the single-pass interferometry technique in C- and X-band. Earth’s height model generated by InSAR-SRTM with 90-m horizontal resolution is available while the DEM with 4-to-4.5-m relative accuracy is also available for restricted areas around the world.

Fig. 2: Differential interferometric SAR data collection scheme

 

 

Comparison

InSAR ability to generate topographic and displacement maps in wide applications like earthquakes, mining, landslide, volcanoes has been proven. Although other facilities like GPS, total stations, laser altimeters are also used, comparison between InSAR and these tools reveals its reliability. Laser altimeters can generate high resolution DEM and low resolution displacement maps in contrary to InSAR with the spatial resolution of 25m. However, most laser altimeters record narrow swaths. Therefore, for constructing a DEM by laser altimeter, more overlapping images are required. Displacement map precision obtained by terrestrial surveying using GPS and total stations is similar or better than InSAR. GPS generally provides better estimation of horizontal displacement and with permanent benchmarks slow deformations is monitored for years without being concerned about surface de-correlation. The most important advantage of InSAR over GPS and total stations are wide continuous coverage with no need for fieldwork. Therefore, wide and continuous coverage, high precision, cost effectiveness and feasibility of recording data in all weather conditions are its main privileges. However, it is important that the InSAR displacement result is in the line-of-the-sight direction and to decompose this vector to parallel and normal components the terrestrial data or extra interferograms with different imaging geometry are required. It is shown that DEM generated by photogrammetric method is more accurate than the others. It has approximately 5.5m accuracy for open and 6.5m for forest areas. SRTM X-band DSM is 4m less accurate for open and 4.5m less accurate for forest areas.

Data availability and atmospheric effects limit using InSAR, however processing of its data is challenging. For each selected image pair, several processing steps have to be performed. One of the current challenges is to bring the techniques to a level where DEM generation can be performed on an operational basis. This is important not only for commercial exploitation of InSAR data, but also for many government and scientific applications. Multi pass interferometry is affected by the atmospheric effects. Spatial and temporal changes due to the 20% of relative humidity produce an error of 10cm in deformation. Moreover, for the image pairs with inappropriate baseline the error introduced to the topographic maps is almost 100m. In topographic mapping this error can be reduced by choosing interferometric pairs with relatively long baselines, while in the displacement case the solution is to average independent interferograms.

 

Why InSAR DEMs are better?

Distinction between SAR imaging and the optical systems are more profound than the ability of SAR to operate in conditions that would cause optical instruments to fail. There are basic differences in the physical principles dominating the two approaches. Optical sensors record the intensity of radiation beamed from the sun and reflected from the features. The intensity of the detected light characterizes each element of the resulting image or pixel. SAR antenna illuminates its target with coherent radiation. Since the crests and troughs of the emitted electromagnetic wave follow a regular sinusoidal pattern, both the intensity and the phase of returned waves can be measured.

InSAR has some similarities to stereo-optical imaging in that two images of the common area, viewed from different angles, are appropriately combined to extract the topographic information. The main difference between interferometry and stereo imaging is the way to obtain topography from stereo-optical images. Distance information is inherent in SAR data that enables the automatic generation of topography through interferometry. In other words DEMs can be generated by SAR interferometry with greater automation and less errors than optical techniques. Moreover, using DInSAR surface deformations can be measured accurately.

Different DEM generation methods of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) stereoscopy, ERS tandem InSAR, and SRTM-InSAR are used. Both the ERS-InSAR and SRTM DEMs are free of weather conditions, but ASTER DEM quality may be affected by cloud coverage in some local areas. InSAR has the potential of providing DEMs with 1-10cm accuracy, which can be improved to millimeter level by DInSAR. Its developments are rapid however it is our requirements that say which one is better for use.   

Examples of our practical studies and achievements

In course of the years of studies and verifications by our Microwave Remote Sensing Group good and valuable achievements gained on InSAR technology applications. In this paper the results and outcomes related to the generation of DEMs using Interferometric SAR are given as example. 

(Image credit: Parviz Tarikhi)

Fig. 3: Satellite view of Bam area with the pre and post quake photos of the Citadel Bam as insets  

(I) Advent of the awful earthquake of December 26, 2003 in Bam, Iran has drawn the attention of the many scientific and humanitarian organizations to study the phenomenon and its causes and origins as well as its impacts and developments. Fig.3 shows the position of the area on a satellite image with the photos of the pre and post quake views of the Citadel Bam as insets. In Fig.4, the left image is the topo-DInSAR product acquired from the Envisat-ASAR data of 11 June and 3 December 2003, while the image at right is the topo-DInSAR product of the 3 December 2003 and 7 January 2004. The team conducted by the author used the data provided by ESRIN, and the DORIS and IDIOT softwares to generate the products. The middle image obtained by NASA scientists is the 3-D perspective view of vertical displacement of the land surface south of Bam during the 3.5 years after the 6.6 earthquake of December 26, 2003 that is derived from analysis of radar images. Blue and magenta tones show where the ground surface moved downward; yellow and red tones show upward motion (particularly in south of Bam). Displacements are superimposed on a false-color Landsat Thematic Mapper image taken on October 1, 1999. In the right image that is obtained from the ASAR data pre and post earthquake the curl-shape pattern south of Bam is distinguishable where such the torsion is not visible in the left image that obtained from pre-earthquake data. For the left image the normal baseline is 476.9 m and parallel baseline is 141.6 m, while for image at right the normal baseline is 521.9 m and the parallel baseline is 268.3 m. The left image demonstrates that the related interferogram includes four lobes. Since the displacement in the east is greater than that in the west, the related lobes are larger.  The displacements measured along the radar line-of-sight direction are 30cm and 16cm at south-east and north-east lobes of the interferogram, respectively. However, the displacement related to the western part of the area is about 5cm along the radar line-of-sight direction.    

 

(Image credit: Parviz Tarikhi)

Fig.4: Bam area’s topo DInSAR images and DEM 

(II) In Fig.5, the top images at the left and right show the tandem amplitude data of 12 and 13 August 1999 (4 and 5 days pre-quake) of Izmit, Turkey as master and slave images respectively. The research team conducted by the author used the ERS-1&2 data provided by ESRIN, and the Earth-view and SAR Toolbox softwares to generate the variety of related products. The normal baseline for the image pair is 224.2 m and parallel baseline is 91.1 m. The image at the bottom left is the coherence image while the right image at the bottom depicts the DEM image where the interferogram is overlaid on it. For each product the relevant histogram is seen as inset. The similarity of the histograms of master and slave images is considerable due to the high correlation of the images that is clearly seen in coherence image. It is important to note that lowest coherence values (darkest values) correspond both to steep slopes or vegetated areas (especially visible in the lower part of the image) and to the lakes (image center and left). DEMs generated from the tandem images are accurate because of the high correlation between master and slave images. Although both the master and slave images are pre-earthquake data of the 7.8 earthquake of August 17, 1999, the strain in the disaster area is visible a week before the quake. It could be a useful precursor for the advent of a disaster like the earthquake in Izmit.

(Image credit: Parviz Tarikhi)

Fig.5: InSAR products of Izmit area in western Turkey

Acknowledgement

The SLC images for both the Izmit and Bam areas for our study was generously provided by the European Space Research Institute (ESRIN) of the European Space Agency (ESA) in 1999 and 2004.

Presenter’s CV 

Parviz Tarikhi specializes in radar remote sensing since 1994 and heads the Microwave Remote Sensing Research Core at the Mahdasht Satellite Receiving Station in Karaj, Iran. He holds a PhD in physics focusing on microwave remote sensing. He has been involved with the United Nations Committee on the Peaceful Uses of Outer Space (UN-COPUOS) since 2000, including as Second Vice-Chair and Rapporteur in 2004-06 of the committee bureau. Since 2001 he has co-chaired Action Team number 1 of UNISPACE-III with the mission ‘to develop a comprehensive worldwide environmental monitoring strategy’. From 2004-07 he led the Office for Specialized International Cooperation of the Iranian Space Agency. He has made in the mean time years of research and study on the developments and status of space science and technology with a particular focus on Iran. For more information please see:

 https://parviztarikhi.wordpress.com

   

  

 

 

 

 

 

Note: Apart from any fair dealing for the purposes of research or private study, or criticism or review, this publication may only be reproduced, stored or transmitted in any form or by any means with the prior permission in writing of the author.

BACK                                                        HOME

 

———————————————————————-

Dear visitor,

Thanks for visiting this page!

I would appreciate your comments and suggestions on the content and set-up.

Parviz Tarikhi

 

 

Advertisements

32 Responses to “Radar DEM generation; achievements and benefits”

  1. khaled Says:

    Many thanks ,great job, once again Excellent.

  2. diablo 3 beta download Says:

    Hey! Do you know if they make any plugins to assist with Search Engine Optimization? I’m trying to get my blog to rank for some targeted keywords but I’m not seeing very good results. If you know of any please share. Appreciate it!

  3. Janella Skeels Says:

    Some really prize articles on this site, saved to bookmarks .

  4. Cheap Underwater Digital Camera Says:

    Thank you I just Linked To You…

    I am entirely perplexed about this, but are you sure you happen to be to the suitable monitor? Please publish examples…

  5. First Grade Reading Books Says:

    Trackback Priority…

    I’ll right away grasp your rss feed as I can’t to find your email subscription link or newsletter service. Do you’ve any? Please allow me know in order that I may subscribe. Thanks….

  6. FZRaymond Says:

    Если понравится буду выкладывать еще.

    http://opera.timedom.org/E-Ten.html

    скачать оперу мини для Hi Tech HT 210i

  7. bicentennial silver proof set Says:

    Nice Focus….

    I really like the direction you’ve decided to take this blog….

  8. camisetas de barcelona Says:

    Thank you for indicating your thinking within this web site. In the same way, a myth in connection with banking institutions wishes at any time sharing foreclosures would be that the standard bank won’t have my payments. We have a fair bit of your energy this loan company usually requires payments from time to time. For anyone who is also deep inside pit, they should often drive for you to spend the money for precise transaction altogether. Having said that, i’m not really proclaiming that they’re going to have any kind of installments in any respect. For those who along with the standard bank may find a way to operate anything away, the exact real estate foreclosure treatment may perhaps prevent. On the other hand, if you ever keep miss repayments in the new strategy, the house or property house foreclosures practice can make up in which this still left out of.

  9. Credit Cards For Bad Credit Says:

    Wow this site rocks…

    Many thanks for writing this in depth publish, it was so valuable for me!…

  10. Chicken Coops Says:

    Great Job…

    keep bringing more please I love this website…

  11. car accident reports Says:

    Informative and precise…

    Its difficult to find informative and precise information but here I found…

  12. inoplelsontix Says:

    del plz erorr

    • Parviz Tarikhi Says:

      Would you please comment?!

  13. salazartwin Says:

    Hack again?!

    • Parviz Tarikhi Says:

      What for?

  14. Dagmarpcpmnet Says:

    шайтанама

    • Parviz Tarikhi Says:

      Could you please comment?!

  15. billot Says:

    Thank you a lot for giving everyone remarkably wonderful chance to read in detail from this website. It is often very awesome plus packed with a great time for me and my office friends to visit your web site at the least thrice a week to learn the newest issues you have got. Of course, I’m so always astounded with the amazing secrets you serve. Selected 1 ideas in this article are completely the most beneficial I have ever had.

    • Parviz Tarikhi Says:

      Thanks a lot for your support and kind words.

  16. paloma119Ronaldo Says:

    Шалки Валки

    • Parviz Tarikhi Says:

      – Would you please comment?!

  17. orsak Says:

    Virtually all of the things you claim is astonishingly appropriate and it makes me wonder the reason why I hadn’t looked at this with this light before. This particular piece really did turn the light on for me as far as this specific subject goes. Nevertheless there is just one factor I am not really too cozy with so whilst I make an effort to reconcile that with the actual core idea of your point, let me observe what the rest of your visitors have to say.Very well done.

  18. Tehaniosechkin2000 Says:

    привет всем

    • Parviz Tarikhi Says:

      What does this mean; I do not know Russian language.

  19. chalcographic Says:

    It’s not afterward time and energy to terminate or simply fault the particular e-book.

  20. Joshalinokulow Says:

    ОПАПУЛИЧКИ

    • Parviz Tarikhi Says:

      I do not know Russian language. Could you please translate your comment.
      Regards, Parviz

  21. macbook pro memory upgrade Says:

    WONDERFUL Post.thanks for share..more wait ..

  22. dell streak 7 Says:

    WONDERFUL Post.thanks for share..more wait ..

  23. PetrDucknoicign Says:

    I am sorry that there can not write, but when uploading photos to Proville writes error java, who can help do that ?

  24. otdrrяemontkv Says:

    ремонт в коммунальной квартире Реставрация и отделка офисов, квартир или загородных домов – это насущная в настоящее время проблема большинства людей, решившихся, наконец, привести в порядок свое жилище. Мало того, что улучшение – это занятие трудоемкое и напряженное, еще и найти заслуживающую доверия компанию, занимающуюся улучшением комнат или загородных домов сегодня почти нельзя. Дело в том, что все меньше и меньше в России осталось организаций, делающих свою работу качественно, эффективно и недорого. собственник ремонт квартиры
    Однако наша организация является приятным исключением из этого неприятного правила. Мы вот уже много лет занимается улучшением помещений, комнат, загородных домов и еще ни разу нам не приходило от наших клиентов ни неприятных отзывов, ни каких-либо жалоб. кредит на ремонт квартиры сбербанк Большой опыт наших специалистов, а также высокий профессионализм всех без исключения наших сотрудников, поможет сделать хорошо, качественно и – главное – доступно как косметический, так и капитальный ремонт. ремонт квартир интерьер Современное оборудование, новые технологии планировки и дизайна, а также индивидуальный подход к каждому клиенту и гибкая система скидок сделали нашу организацию самой успешной на сегодня фирмой, оказывающей услуги ремонту и отделке помещений, комнат или домов. ремонт квартиры хрущевки
    ремонт квартир ремонт квартиры

    • Parviz Tarikhi Says:

      I do not know Russian language. Could you please translate your comment.
      Regards,

      Parviz


Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: