|
1
|
/* merge 9 tiles in Oregon around Portland for testing methods
|
|
2
|
|
|
3
|
|
|
4
|
&if %:program%_ eq ARC_ &then grid
|
|
5
|
&if %:program%_ ne GRID_ &then
|
|
6
|
&return This program must be run from ARC or GRID.
|
|
7
|
|
|
8
|
&sv tiledir \\jupiter\\organisms\topo\tiles
|
|
9
|
&sv outdir \\jupiter\\organisms\topo\experimental
|
|
10
|
&sv name srtmv41
|
|
11
|
&sv vernum 1
|
|
12
|
|
|
13
|
/* This can be extended to larger areas (beyond 3x3 used here) but at some stage the
|
|
14
|
/* tilelist variable is going to get too long - AML has a limited line length
|
|
15
|
/* At this point, it would be better to tile by longitude band i.e. do a merge for
|
|
16
|
/* each value of w, then merge the resulting strips as a separate step
|
|
17
|
&sv tilelist
|
|
18
|
&do w = 122 &to 124
|
|
19
|
&do n = 44 &to 46
|
|
20
|
&sv tilelist %tilelist%, %tiledir%\w%w%\n%n%\%vernum%\%name%\w%w%n%n%
|
|
21
|
&end
|
|
22
|
&end
|
|
23
|
|
|
24
|
%outdir%\srtmv41 = merge(%tilelist%)
|
|
25
|
:q
|
