|
1 |
May 5, 2011 / Rick Reeves
|
|
2 |
ReadMe file: Shell scripts to create the mosaic (3 arcsecond resolution)
|
|
3 |
of CGIAR/SRTM and ASTER GDEM imagery spanning the 60 Degree North Latitude
|
|
4 |
boundary within western Canada
|
|
5 |
|
|
6 |
Three shell scripts used:
|
|
7 |
|
|
8 |
1) "mosaicCgiarSrtmBdySRTM.sh" : creates mosaic from CGIAR/SRTM 5 degree tiles
|
|
9 |
(3 arcsecond resolution) converted to .tif format from original ArcMap ASCII
|
|
10 |
format.
|
|
11 |
|
|
12 |
2) "mosaicCgiarSrtmMBdyAster.sh" : creates mosaic from ASTER GDEM 1 degree tiles
|
|
13 |
(1 arcsecond resolution) converted to .tif format from original ArcMap ASCII
|
|
14 |
format. Note: incoming image tiles are resampled to 3 arcsecond resolution
|
|
15 |
to create outgoing image mosaic,
|
|
16 |
|
|
17 |
3) "mosaicSrtmAsterPartsBdy.sh" : creates mosaic from the outputs from scripts
|
|
18 |
1) and 2) (3 arcsecond resolution):
|
|
19 |
mergeCgiarAsterBdySRTM_BL.tif
|
|
20 |
mergeCgiarAsterBdyASTER_BL.tif
|
|
21 |
|
|
22 |
To create the ASTER / CGIAR boundary mosaic:
|
|
23 |
|
|
24 |
1) make current directory '/data/project/organisms/rcr/AsterCgiarMerge/'
|
|
25 |
2) run the script "mosaicCgiarSrtmBdySRTM.sh"
|
|
26 |
3) run the script "mosaicSrtmMBdyAster.sh"
|
|
27 |
4) run the script "mosaicSrtmAsterPartsBdy.sh"
|
|
28 |
|
|
29 |
Note: These initial scripts do NOT specify a new image extent for any of the mosaic components.
|
|
30 |
(gdalwarp parameter '-te xmin ymin xmax ymax' used to specify the extent)
|
|
31 |
So for these initial *test* runs, the output image extent is derived from the input files.
|
|
32 |
|
|
33 |
In 'production' runs, the extent of output mosaic tile compoents will be specified using -te;
|
|
34 |
we may also set other 'gdalwarp' parameters to properly establish other output image parameters.
|
|
35 |
|
ReadMe file describes how to use shell scripts (also checked in) to create example Srtm/Aster DEM image mosaic of Western Canada region spanning 60 deg N Latitude boundary