Best viewed with Firefox

ROI_PAC processing for Modern Geodetic Methods (Page 5 of 6)

Processing Table of Contents Clean up

Creating Figures:

  • If all goes well, there are only one or two figures that you may want to create: the final interferogram with all modeled orbit and topography removed, and the final unwrapped interferogram with values in LOS deformation.

Creating Interferogram plot:

Unfortunately, the script can work so automatically that it doesn't by default generate an interferogram file that is georeferenced and ready to image with other software such as GMT. Thus, we will need to do a couple of file conversions before we start plotting.
  • Be sure you're in the correct Directory
    % newgrp MGM
    % cd ~/MGM/HectorMine/int_990915_991020/  
  • Convert the interferogram from a complex to hgt format.
  • Now, we must geocode the interferogram (convert it from satellite to an earth reference frame) using the 'geomap_4rlks.trans' transformation grid.
    % geomap_4rlks.trans 

    Be careful to only transform files with the same number of looks (e.g. 4rlks) as the program will assume so, and perform a blind pixel-by-pixel transformation. In order to downsample a raw sized image to one that matches your current transformation grid use ''.

    % 4  # this will create
  • We will also need to convert the HGT format used by ROI_PAC to a grid format read by GMT. This will require matlab, and can be done either by bringing up a matlab window or calling the commands directly from the command line. I show the later here because of its expediency. Matlab simply calls a couple of simple scripts that I wrote called wrap2grd and internally readforwrapped to convert the input file into a NetCDF file for GMT.
    The format for wrap2grd is:
                     AP = 0=Amp, 1=Phase
                     infile = hgt inputfile (must have .rsc file too)
    		 outfile = output NetCDF file
    Looks like the grdwrite component of wrap2grd doens't work on the newest version of matlab. Use matlab2010 on PELE!!
    % matlab  -nosplash -nojvm -nodesktop  << EOF
        addpath '/home/anewman/InSAR/matlab_common'
  • Now we will copy over and run a GMT script, HMwrapped.gmt written specifically for the Hector Mine region. A version of this script can be easily enough modified to plot interferograms for any region.
    % cp ~anewman/MGM/HectorMine/int_990915_991020/HMwrapped.gmt .  # copy my version so you can later modify it
    % HMwrapped.gmt 
  • To create a high resolution png image file suitable for viewing rapidly online you can run pstoimg. Here, I am forcing the image density in dots per inch (DPI) to be 200.
    % pstoimg -antialias -density 200  # this may take a while
    # the program can now be viewed with gqview or xv
    % xv

Creating unwrapped LOS deformation plot:

Because the assumed final product of ROI_PAC is a deformation plot in Line-of-Site (LOS) length, this plot is considerably more straightforward to make.
  • Create an unwrapped deformation grd
    % matlab  -nosplash -nojvm -nodesktop  << EOF
        addpath '/home/anewman/InSAR/matlab_common'
  • Create an unwrapped RADAR amplitude grd
    % matlab  -nosplash -nojvm -nodesktop  << EOF
        addpath '/home/anewman/InSAR/matlab_common'
  • Now, with only minor modifications, you can change HMwrapped.gmt to a couple of new scripts for Amp and Phase plots. You can either do this on your own or copy my versions.
    % cp ~anewman/MGM/HectorMine/int_990915_991020/HMamp.gmt .
    % HMamp.gmt geo_990915-991020.amp
    % cp ~anewman/MGM/HectorMine/int_990915_991020/HMunw.gmt .
    % HMunw.gmt geo_990915-991020.unw

    Example Images:
    1. Radar Amplitude (990915-991020.amp)
    2. Radar Interferogram: no correction (
    3. Radar Interferogram: orbit correction (
    4. Radar Interferogram: SIMulated orbit and topographic corrections (
    5. Radar Interferogram: orbit and topographic corrections (
    6. LOS Displacement (990915-991020.unw)

Processing Clean up

Course Home | | Updated: Mon Apr 9 14:46:47 EDT 2012