|
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
|
|
