Skip to content

CheatSheat.md

Latest commit

 

History

History
143 lines (115 loc) · 3.97 KB

CheatSheat.md

File metadata and controls

143 lines (115 loc) · 3.97 KB

ISCE/Python Cheat Sheet

A collection of common and helpful commands using ISCE, GDAL, GMT, and Python

Command Line

Save unwrapped phase as single band geotiff with control over colorbar

isce2geotiff.py -i filt_topophase.unw.geo -o filt_topophase.unw.geo.tif -b 2 -c -10 10

Create kmz to view in Google Earth

(Note: Google Earth Pro is now free, and can open geotifs directly, so do not really need this...)

gdal_translate -of KMLSUPEROVERLAY filt_topophase.unw.geo.tif filt_topophase.unw.geo.kmz

Convert unwrapped radians to displacements [cm]

(Note: example for C-band Sentinel-1A Wavelength 0.555 cm)

isce2gis.py vrt -i los.rdr.geo
gdal_calc.py -A filt_topophase.unw.geo.vrt --A_band=2 --calc="A*0.05546576/12.5663706" --outfile=filt_topophase.unw_m.geo  --format=ENVI --NoDataValue=-9999 --overwrite

Help on particular ISCE component

iscehelp.py -t Sensor -a sensor=SENTINEL1A

Look down an ISCE file

looks.py -i filt_topophase.flat -r 4 -a 4 -o filt_topophase.flat.4lks

Get quick stats on an image

For example mean incidence and heading angles. Band 1 is incidence (radar los to surface normal), Band 2 is heading (positive clockwise from due east):

isce2gis.py vrt -i los.rdr.geo
gdal_edit.py -a_nodata 0.0 los.rdr.geo.vrt
gdalinfo -stats los.rdr.geo.vrt

Geocode File

isce2gis.py vrt -i filt_topophase.flat --lat lat.rdr --lon lon.rdr

Alternatively, run single step of main processing script with xml input:

topsApp.py —-dostep=geocode topsApp_geocodeonly.xml

where topApp_geocodeonly.xml has:

<?xml version="1.0" encoding="UTF-8"?>
<topsApp>
	<property name="geocode list">
		<value>["filt_topophase.flat.4lks"]</value>
	</property>
</topsApp>

Unwrap specific file with Snaphu

Merge wrapped, geocoded, Sentinel-IW subswaths

gdalwarp -of ENVI -ot CFloat32 -srcnodata 0 -dstnodata 0 20160322_20160415_*/merged/*filt_topophase.flat.geo.vrt filt_topophase_merged.flat.geo.vrt

Put separate geocoded phase files on same grid (with GDAL + GMT)

isce2gis.py envi -i date1_date2_filt_topophase.unw.geo
isce2gis.py envi -i date3_date4_filt_topophase.unw.geo
gdal_translate -of GMT -b 2 date1_date2_filt_topophase.unw.geo raster1.grd
gdal_translate -of GMT -b 2 date3_date4_filt_topophase.unw.geo raster2.grd
grdcut raster2.grd -Rraster1.grd -Graster2_regrid.grd -N

Extract wrapped phase from complex-valued files

isce2gis.py vrt -i filt_topophase.flat.8alks_8rlks.geo

gdal_calc.py --type Float32 -A filt_topophase.flat.8alks_8rlks.geo.vrt --calc="numpy.angle(A)" --outfile=filt_topophase.flat.8alks_8rlks.geo.phs.tif --NoDataValue=0.0 --overwrite

Create RGBA Geotiff for Google Earth

gdaldem color-relief -alpha filt_topophase.flat.8alks_8rlks.geo.phs.tif colors.txt filt_topophase.flat.geo.8alks_8rlks.phs.rgba.tif

where colors.txt is something like this:

-3.14   46 154 88 255
-1.57   251 255 128 255
0       224 108 31 255
1.57    200 55 55 255
3.14    215 244 244 255
nv      0   0   0   0

Python Module

Load Processing information from PICKLE (when running insarApp.py with --steps)

import isce
import pickle
with open('PICKLE/updatepreprocinfo', 'rb') as f:
    insar = pickle.load(f)
# insar is a dictionary in this case:
insar['sensor']

Plot array in map view and profile

Other stuff

Download All Sentinel-1A SLC products for a particular point and orbit from ASF

export LAT=1.22
export LON=-77.37
export ORBIT=120
curl https://api.daac.asf.alaska.edu/services/search/param?intersectsWith=point%28$LON+$LAT%29\&platform=Sentinel-1A\&processingLevel=SLC\&relativeOrbit=$ORBIT\&output=metalink > query$ORBIT.metalink
aria2c --http-auth-challenge=true --conf-path=asf.conf query$ORBIT.metalink

Where asf.conf contains the following:

http-user=<CHANGE>
http-passwd=<CHANGE>
max-concurrent-downloads=5
check-certificate=false
allow-overwrite=false
auto-file-renaming=false
always-resume=true