David Sandwell, Fwebryary 16, 2012

This directory contains code to establish zones of sampling rate for the LOS InSAR data.
Xiaopeng takes the file sample.grd and uses it to subsample the LOS Data.  Here is the approachs.

Hi Xiaopeng,

I put a file on scratch called sample.tar.  Inside there is a file called sample.grd.  This is a grid of number 1-7 based on the second invariant of the strain rate.

1 - sample every point
2 - sample every other point
.
4-7 - sample every 5th point
.
.

You can look at the ps file to see that there will be a 1-sample along the creeping faults.
So here is is a way to use this.  Run blockmedian on all your LOS tiles so there will be 4 sets of tiles.
A
1- downsample to -I.01/.008  (800 m resolution)
2- downsample to -I.02/.016  (1600 m resolution)
3- downsample to -I.04/.032  (3200 m resolution)
4-7 downsample  to -I.08/.064 (6400 m resolution)

I'm guessing at the numbers but you get the idea.  We will be running Bridget's model at 500 m.
Then your script will read the sample.grd and select the appropriate subsampled LOS grid.
I think this will keep all the LOS information.