May 5, 2011 / Rick Reeves
ReadMe file: Shell scripts to create the mosaic (3 arcsecond resolution) 
of CGIAR/SRTM and ASTER GDEM imagery spanning the 60 Degree North Latitude 
boundary within western Canada

Three shell scripts used:

1) "mosaicCgiarSrtmBdySRTM.sh" : creates mosaic from CGIAR/SRTM 5 degree tiles
    (3 arcsecond resolution) converted to .tif format from original ArcMap ASCII
    format.

2) "mosaicCgiarSrtmMBdyAster.sh" : creates mosaic from ASTER GDEM 1 degree tiles
   (1 arcsecond resolution) converted to .tif format from original ArcMap ASCII
   format. Note: incoming image tiles are resampled to 3 arcsecond resolution
   to create outgoing image mosaic,

3) "mosaicSrtmAsterPartsBdy.sh" : creates mosaic from the outputs from scripts
    1) and 2) (3 arcsecond resolution): 
                                            mergeCgiarAsterBdySRTM_BL.tif 
                                            mergeCgiarAsterBdyASTER_BL.tif

To create the ASTER / CGIAR boundary mosaic: 

1) make current directory '/data/project/organisms/rcr/AsterCgiarMerge/'
2) run the script "mosaicCgiarSrtmBdySRTM.sh"
3) run the script "mosaicSrtmMBdyAster.sh"
4) run the script "mosaicSrtmAsterPartsBdy.sh"

Note: These initial scripts do NOT specify a new image extent for any of the mosaic components.
      (gdalwarp parameter '-te xmin ymin xmax ymax' used to specify the extent)
      So for these initial *test* runs, the output image extent is derived from the input files.
 
      In 'production' runs, the extent of output mosaic tile compoents will be specified using -te;
      we may also set other 'gdalwarp' parameters to properly establish other output image parameters.