Extracted Pages 1

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

How to make an InSAR time series from Sentinel-1 TOPS data: Kilauea

Introduction

This document presents a step-by-step recipe for creating InSAR time series from Sentinel-1
TOPS data. We aim to guide you, the user, through system set-up, downloading necessary data
and orbit files, coregistering chosen images, creating interferograms, unwrapping those
interferograms, and running the Small BAseline Subset Analysis (SBAS) to create the final time
series. In the example below, we chose a region of interest on the big island of Hawaii to look at
the deformation surrounding the eruption of Kilauea and the ensuing earthquake that occurred
there in 2018. This recipe can be adapted to any region of interest, and we include tips and
tricks that apply to various types of areas (both areas of good coherence and poor coherence).

To begin this tutorial, we wish to prepare the user for the path ahead with three tips:
(1) This document assumes you have a working knowledge of UNIX or LINUX systems. While
we provide quick commands to accomplish the required steps in most cases, if you are
new to UNIX/LINUX, we suggest you run through a tutorial of these systems so you are
prepared (http://www.ee.surrey.ac.uk/Teaching/Unix/ ).
(2) All updated software and help pages are located at Github under GMTSAR
(https://github.com/gmtsar/gmtsar ). If you find yourself troubleshooting, or having any
problems, check out the issues tab at this site to see if someone else has encountered
your issue (https://github.com/gmtsar/gmtsar/issues ).
(3) We also suggest that you make an account with the Alaska Satellite Facility (ASF) and
with Copernicus (for Sentinel-1 orbit files) ahead of time for downloading data (you will
need your username and password in the steps below).

If you haven’t already done so, visit the ASF EarthData page and click on their VERTEX
tool to find the sign in/register page:
https://asf.alaska.edu/

For European Space Agency’s Copernicus System, visit:


https://scihub.copernicus.eu/dhus/#/home

Lastly, remember to enjoy the fact that we can visualize minute displacements on the Earth’s
surface through the use of orbiting satellites that measure radar range pulses like a bat from
outer space. There is a remarkable amount of high-quality data these days that we can gain
access to, and it is amazing how we can use these data to study our planet’s various changes.
Creating an InSAR time series is just one of those ways we can study our planet from space, so
go forth and experiment and discover!

This document was updated and prepared by Katherine Guns and Julie Gevorgian in October 2020, with
assistance from Xiaohua Xu and David Sandwell. [Last updated 02/2022.]
Table of Contents

1. General Setup
2. Preparing the Topography Grid Files
3. Data Selection and Download for Region of Interest
a. Data Selection
b. Download Data
c. Adjusting Region of Interest
4. Downloading Precise (and/or Restituted) Orbits
a. Navigating the ESA Website
b. Downloading Orbit Files through Command Line
c. Downloading Restituted Orbits (if desired)
5. Choosing Master/Reference Image and Aligning
a. Selecting a Reference Image from a Baseline Plot
b. Aligning Secondary Images to Reference Image
c. Creating Baseline Table
6. Running Interferograms
a. Preparing List of Interferograms
b. Run One Interferogram to Test Settings
c. Run All Interferograms
7. Merging Across Subswaths
8. Unwrapping Interferograms
a. Unwrapping in Regions of Good Correlation
b. Unwrapping in Region of Poor Correlation
i. Stacking Coherence Grids
ii. Creating a Mask
9. Running Small BAseline Subset (SBAS) Analysis
a. Prepare the Input Files
b. Run SBAS
10. How to Compare with GNSS Data
a. Project GNSS Velocities/Displacements into LOS
b. Extract Point Velocities/Displacements from InSAR LOS Grids
11. Correcting InSAR Interferograms with GNSS Displacements
12. Note on Incorporating Ascending and Descending Time Series
13. References

2
1. General Setup
We need to start by setting up our directories in the necessary pattern. Make a top level
directory structure for ascending or descending data all the following directory names:

[Note: Depending on how GMTSAR is installed, a copy of batch_tops.config can be found in:
/usr/local/GMTSAR/gmtsar/csh, though the path to the GMTSAR directory can vary. If you
cannot locate it, visit the Github repository for GMTSAR to download a copy to place in your
F1,F2,F3 directories.]

2. Preparing the Topography DEM Grid Files

Prepare a topography grid (dem.grd)


a. Using the DEM generator web interface (http://topex.ucsd.edu/gmtsar/demgen/),
prepare a topography grid for your area of interest (lon -157.0 to -154.2, lat 18.0
to 20.4 for Kilauea). This will download as dem.grd, along with a kml companion
file. View the kml file in Google Earth and verify that it completely covers the
area (as shown in the image following 3c).

[Note: The DEM generator can only generate 4 degree wide topography areas at a
time. In addition, if your area is near the equator, you cannot have a dem area
straddling the equator--you need to download above and below the equator areas
separately. To stitch multiple dem grids together, we recommend using gmt grdpaste
(see examples here: http://gmt.soest.hawaii.edu/doc/latest/grdpaste.html )]

b. Move the dem.grd file to your upper level topo directory (inside asc or des
directories).

3
c. Next link the dem.grd file to a topo directory inside the F1/F2/F3 directories..
Create a symbolic link in each:
cd asc/F1/
mkdir topo
cd topo
ln -s ../../topo/dem.grd .

d. Also link the dem.grd to the merge directory (doesn’t need to be inside a topo
folder here).
cd asc/merge/
ln -s ../topo/dem.grd .

3. Data Selection and Download for Region of Interest

a. Data Selection: Plan your overall processing strategy


i. Go to the Sentinel open Copernicus data hub GUI, login or make an
account, and search for data:
(https://scihub.copernicus.eu/dhus/#/home)
ii. Use the Pan and Box tools (top right corner) to zoom in on the area of
interest, for example, The Big Island, Hawaii (lon -157.0 to -154.2, lat 18.0
to 20.4, an area that spans the entire frame for making the topography).
Draw a box over this area.
iii. Using the drop down search menu (three horizontal lines icon), set the
search parameters for the sensing dates, product type (SLC), polarisation
(VV) and the sensor mode (IW). As you can see in the example below, we
had ordered the search results by ascending frames (this can be changed
to descending). Ascending frames are oriented NW/SE while descending
frames are oriented NE/SW.

You might also like