1 |
7526fb1c
|
Jim Regetz
|
# Script for creating Slope and Aspect dem (tif) for SRTM aster straddle, and both uncorrected
|
2 |
|
|
# fused layer (straddle) and blend gaussian fused layer (blendgau) for each nine region.
|
3 |
|
|
#
|
4 |
|
|
# 9-Dec-2011
|
5 |
|
|
# Yuina Nunokawa
|
6 |
|
|
|
7 |
|
|
|
8 |
|
|
|
9 |
|
|
|
10 |
|
|
|
11 |
|
|
|
12 |
|
|
------------------------------------------------------------------------
|
13 |
|
|
|
14 |
|
|
# SRTM
|
15 |
|
|
|
16 |
|
|
export SRTMDIR="/data/project/organisms/DEM/Yuni/Data/srtm"
|
17 |
|
|
|
18 |
|
|
|
19 |
|
|
# SRTM Slope
|
20 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_e020e059_below.tif $SRTMDIR/srtm_e020e059_below_below_s.tif
|
21 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_e060e099_below.tif $SRTMDIR/srtm_e060e099_below_below_s.tif
|
22 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_e100e139_below.tif $SRTMDIR/srtm_e100e139_below_below_s.tif
|
23 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_e140e179_below.tif $SRTMDIR/srtm_e140e179_below_below_s.tif
|
24 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_w020e019_below.tif $SRTMDIR/srtm_w020e019_below_below_s.tif
|
25 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_w060w021_below.tif $SRTMDIR/srtm_w060w021_below_below_s.tif
|
26 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_w100w061_below.tif $SRTMDIR/srtm_w100w061_below_below_s.tif
|
27 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_w140w101_below.tif $SRTMDIR/srtm_w140w101_below_below_s.tif
|
28 |
|
|
gdaldem slope -s 111120 $SRTMDIR/srtm_w180w141_below.tif $SRTMDIR/srtm_w180w141_below_below_s.tif
|
29 |
|
|
|
30 |
|
|
|
31 |
|
|
# SRTM Aspect
|
32 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_e020e059_below.tif $SRTMDIR/srtm_e020e059_below_a.tif
|
33 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_e060e099_below.tif $SRTMDIR/srtm_e060e099_below_a.tif
|
34 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_e100e139_below.tif $SRTMDIR/srtm_e100e139_below_a.tif
|
35 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_e140e179_below.tif $SRTMDIR/srtm_e140e179_below_a.tif
|
36 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_w020e019_below.tif $SRTMDIR/srtm_w020e019_below_a.tif
|
37 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_w060w021_below.tif $SRTMDIR/srtm_w060w021_below_a.tif
|
38 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_w100w061_below.tif $SRTMDIR/srtm_w100w061_below_a.tif
|
39 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_w140w101_below.tif $SRTMDIR/srtm_w140w101_below_a.tif
|
40 |
|
|
gdaldem aspect -s 111120 $SRTMDIR/srtm_w180w141_below.tif $SRTMDIR/srtm_w180w141_below_a.tif
|
41 |
|
|
|
42 |
|
|
|
43 |
|
|
------------------------------------------------------------------------
|
44 |
|
|
|
45 |
|
|
# ASTER
|
46 |
|
|
|
47 |
|
|
export ASTDIR="/data/project/organisms/DEM/Yuni/Data/aster2"
|
48 |
|
|
|
49 |
|
|
# Create slop files of aster straddle layer (straddle) for each nine region
|
50 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_e020e059_straddle.tif $ASTDIR/aster2_e020e059_straddle_s.tif
|
51 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_e060e099_straddle.tif $ASTDIR/aster2_e060e099_straddle_s.tif
|
52 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_e100e139_straddle.tif $ASTDIR/aster2_e100e139_straddle_s.tif
|
53 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_e140e179_straddle.tif $ASTDIR/aster2_e140e179_straddle_s.tif
|
54 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_w020e019_straddle.tif $ASTDIR/aster2_w020e019_straddle_s.tif
|
55 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_w060w021_straddle.tif $ASTDIR/aster2_w060w021_straddle_s.tif
|
56 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_w100w061_straddle.tif $ASTDIR/aster2_w100w061_straddle_s.tif
|
57 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_w140w101_straddle.tif $ASTDIR/aster2_w140w101_straddle_s.tif
|
58 |
|
|
gdaldem slope -s 111120 $ASTDIR/aster2_w180w141_straddle.tif $ASTDIR/aster2_w180w141_straddle_s.tif
|
59 |
|
|
|
60 |
|
|
|
61 |
|
|
|
62 |
|
|
# Create aspect files of aster straddle layer (straddle) for each nine region
|
63 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_e020e059_straddle.tif $ASTDIR/aster2_e020e059_straddle_a.tif
|
64 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_e060e099_straddle.tif $ASTDIR/aster2_e060e099_straddle_a.tif
|
65 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_e100e139_straddle.tif $ASTDIR/aster2_e100e139_straddle_a.tif
|
66 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_e140e179_straddle.tif $ASTDIR/aster2_e140e179_straddle_a.tif
|
67 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_w020e019_straddle.tif $ASTDIR/aster2_w020e019_straddle_a.tif
|
68 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_w060w021_straddle.tif $ASTDIR/aster2_w060w021_straddle_a.tif
|
69 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_w100w061_straddle.tif $ASTDIR/aster2_w100w061_straddle_a.tif
|
70 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_w140w101_straddle.tif $ASTDIR/aster2_w140w101_straddle_a.tif
|
71 |
|
|
gdaldem aspect -s 111120 $ASTDIR/aster2_w180w141_straddle.tif $ASTDIR/aster2_w180w141_straddle_a.tif
|
72 |
|
|
|
73 |
|
|
|
74 |
|
|
|
75 |
|
|
|
76 |
|
|
# Create slop files of uncorrected fused layer (straddle) for each nine region
|
77 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e020e059_straddle.tif $ASTDIR/fused_e020e059_straddle_s.tif
|
78 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e060e099_straddle.tif $ASTDIR/fused_e060e099_straddle_s.tif
|
79 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e100e139_straddle.tif $ASTDIR/fused_e100e139_straddle_s.tif
|
80 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e140e179_straddle.tif $ASTDIR/fused_e140e179_straddle_s.tif
|
81 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w020e019_straddle.tif $ASTDIR/fused_w020e019_straddle_s.tif
|
82 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w060w021_straddle.tif $ASTDIR/fused_w060w021_straddle_s.tif
|
83 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w100w061_straddle.tif $ASTDIR/fused_w100w061_straddle_s.tif
|
84 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w140w101_straddle.tif $ASTDIR/fused_w140w101_straddle_s.tif
|
85 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w180w141_straddle.tif $ASTDIR/fused_w180w141_straddle_s.tif
|
86 |
|
|
|
87 |
|
|
|
88 |
|
|
|
89 |
|
|
|
90 |
|
|
# Create aspect files of uncorrected fused layer (straddle) for each nine region
|
91 |
|
|
|
92 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e020e059_straddle.tif $ASTDIR/fused_e020e059_straddle_a.tif
|
93 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e060e099_straddle.tif $ASTDIR/fused_e060e099_straddle_a.tif
|
94 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e100e139_straddle.tif $ASTDIR/fused_e100e139_straddle_a.tif
|
95 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e140e179_straddle.tif $ASTDIR/fused_e140e179_straddle_a.tif
|
96 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w020e019_straddle.tif $ASTDIR/fused_w020e019_straddle_a.tif
|
97 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w060w021_straddle.tif $ASTDIR/fused_w060w021_straddle_a.tif
|
98 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w100w061_straddle.tif $ASTDIR/fused_w100w061_straddle_a.tif
|
99 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w140w101_straddle.tif $ASTDIR/fused_w140w101_straddle_a.tif
|
100 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w180w141_straddle.tif $ASTDIR/fused_w180w141_straddle_a.tif
|
101 |
|
|
|
102 |
|
|
|
103 |
|
|
|
104 |
|
|
|
105 |
|
|
# Create slope files of blend gaussian fused layer (blendgau) for nine each region
|
106 |
|
|
|
107 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e020e059_blendgau.tif $ASTDIR/fused_e020e059_blendgau_s.tif
|
108 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e060e099_blendgau.tif $ASTDIR/fused_e060e099_blendgau_s.tif
|
109 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e100e139_blendgau.tif $ASTDIR/fused_e100e139_blendgau_s.tif
|
110 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_e140e179_blendgau.tif $ASTDIR/fused_e140e179_blendgau_s.tif
|
111 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w020e019_blendgau.tif $ASTDIR/fused_w020e019_blendgau_s.tif
|
112 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w060w021_blendgau.tif $ASTDIR/fused_w060w021_blendgau_s.tif
|
113 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w100w061_blendgau.tif $ASTDIR/fused_w100w061_blendgau_s.tif
|
114 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w140w101_blendgau.tif $ASTDIR/fused_w140w101_blendgau_s.tif
|
115 |
|
|
gdaldem slope -s 111120 $ASTDIR/fused_w180w141_blendgau.tif $ASTDIR/fused_w180w141_blendgau_s.tif
|
116 |
|
|
|
117 |
|
|
|
118 |
|
|
|
119 |
|
|
|
120 |
|
|
# Create aspect files of blend gaussian fused layer (blendgau) for each nine region
|
121 |
|
|
|
122 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e020e059_blendgau.tif $ASTDIR/fused_e020e059_blendgau_a.tif
|
123 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e060e099_blendgau.tif $ASTDIR/fused_e060e099_blendgau_a.tif
|
124 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e100e139_blendgau.tif $ASTDIR/fused_e100e139_blendgau_a.tif
|
125 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_e140e179_blendgau.tif $ASTDIR/fused_e140e179_blendgau_a.tif
|
126 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w020e019_blendgau.tif $ASTDIR/fused_w020e019_blendgau_a.tif
|
127 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w060w021_blendgau.tif $ASTDIR/fused_w060w021_blendgau_a.tif
|
128 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w100w061_blendgau.tif $ASTDIR/fused_w100w061_blendgau_a.tif
|
129 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w140w101_blendgau.tif $ASTDIR/fused_w140w101_blendgau_a.tif
|
130 |
|
|
gdaldem aspect -s 111120 $ASTDIR/fused_w180w141_blendgau.tif $ASTDIR/fused_w180w141_blendgau_a.tif
|
131 |
|
|
|