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 process_2pass.pl 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.
% cpx2rmg.pl filt_990915991020sim_ODR_4rlks.int filt_990915991020sim_ODR_4rlks.int.hgt
 Now, we must geocode the interferogram (convert it from satellite to an earth reference frame) using the 'geomap_4rlks.trans'
transformation grid.
% geocode.pl geomap_4rlks.trans filt_990915991020sim_ODR_4rlks.int.hgt geo_filt_990915991020sim_ODR_4rlks.int.hgt
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 pixelbypixel
transformation. In order to downsample a raw sized image to one that matches your current transformation grid use 'look.pl'.
% look.pl 990915991020.int 4 # this will create 990915991020_4rlks.int
 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: wrap2grd(AP,infile,outfile)
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'
wrap2grd(1,'geo_filt_990915991020sim_ODR_4rlks.int.hgt','geo_filt_990915991020sim_ODR_4rlks.int.hgt.grd');
EOF
 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 geo_filt_990915991020sim_ODR_4rlks.int.hgt
 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 geo_filt_990915991020sim_ODR_4rlks.int.ps # this may take a while
# the program can now be viewed with gqview or xv
% xv geo_filt_990915991020sim_ODR_4rlks.int.png
Creating unwrapped LOS deformation plot:
Because the assumed final product of ROI_PAC is a deformation plot in LineofSite (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'
unw2grd(1,'geo_990915991020.unw','geo_990915991020.unw.grd')
EOF
 Create an unwrapped RADAR amplitude grd
% matlab nosplash nojvm nodesktop << EOF
addpath '/home/anewman/InSAR/matlab_common'
unw2grd(0,'geo_990915991020.unw','geo_990915991020.amp.grd')
EOF
 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_990915991020.amp
% cp ~anewman/MGM/HectorMine/int_990915_991020/HMunw.gmt .
% HMunw.gmt geo_990915991020.unw
Example Images:
 Radar Amplitude (990915991020.amp)
 Radar Interferogram: no correction (990915991020_4rlks.int)
 Radar Interferogram: orbit correction (flat_ODR_990915991020_4rlks.int)
 Radar Interferogram: SIMulated orbit and topographic corrections (990915991020sim_SIM_4rlks.int)
 Radar Interferogram: orbit and topographic corrections (filt_990915991020sim_ODR_4rlks.int)
 LOS Displacement (990915991020.unw)


