1
|
# variables derived from the dem following commands in http://www.gdal.org/gdaldem.html#gdaldem_slope
|
2
|
# to check ls tif/Smoothed*.tif | xargs -n 1 -P 30 bash -c $' gdalinfo -mm $1 |grep Max | awk \'{ gsub("[=,]", " ") ; print $3,$4 }\' ' _
|
3
|
|
4
|
|
5
|
# find /mnt/data2/dem_variables/{altitude,aspect,roughness,slope,tpi,tri} -name *.tif | xargs -n 1 -P 10 rm ;
|
6
|
|
7
|
# ls /mnt/data/jetzlab/Data/environ/global/dem/tiles/Smoothed_*.bil | xargs -n 1 -P 12 bash /mnt/data2/dem_variables/scripts/sc1_dem_variables.sh
|
8
|
|
9
|
# scandinavia
|
10
|
# ls /mnt/data/jetzlab/Data/environ/global/dem/tiles/Smoothed_N{55,60,65,70}E0{05,10,15,20,25,30}.bil | xargs -n 1 -P 12 bash /mnt/data2/dem_variables/scripts/sc1_dem_variables.sh
|
11
|
|
12
|
export OUTDIR=/mnt/data2/dem_variables
|
13
|
export file=$1
|
14
|
export filename=`basename $file .bil`
|
15
|
|
16
|
(
|
17
|
|
18
|
# gdal_translate -ot Int16 -co COMPRESS=LZW $file $OUTDIR/tif/$filename.tif
|
19
|
|
20
|
echo altitude variables with file $OUTDIR/tif/$filename.tif
|
21
|
|
22
|
# # median
|
23
|
# pkfilter -m -32768 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/tif/$filename.tif -o $OUTDIR/altitude/median/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
24
|
# # stdev
|
25
|
# pkfilter -m -32768 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/tif/$filename.tif -o $OUTDIR/altitude/stdev/tiles/tmp_$filename.tif -co COMPRESS=LZW -ot Int32 # max 1552385.000 sqrt(1245.947)
|
26
|
# gdal_calc.py -A $OUTDIR/altitude/stdev/tiles/tmp_$filename.tif --calc="sqrt(A)" --type Int16 --overwrite --outfile $OUTDIR/altitude/stdev/tiles/$filename.tif
|
27
|
# rm -f $OUTDIR/altitude/stdev/tiles/tmp_$filename
|
28
|
# # min
|
29
|
# pkfilter -m -32768 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/tif/$filename.tif -o $OUTDIR/altitude/min/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
30
|
# # max
|
31
|
# pkfilter -m -32768 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/tif/$filename.tif -o $OUTDIR/altitude/max/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
32
|
# # mean
|
33
|
# pkfilter -m -32768 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/tif/$filename.tif -o $OUTDIR/altitude/mean/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
34
|
|
35
|
# # starting to use gdaldem to compute variables. Gdaldem use -9999 as no data.
|
36
|
|
37
|
# echo slope with file $OUTDIR/tif/$filename.tif
|
38
|
# gdaldem slope -s 111120 -co COMPRESS=LZW $OUTDIR/tif/$filename.tif $OUTDIR/slope/tiles/$filename.tif # -s to consider xy in degree and z in meters
|
39
|
# # slope median
|
40
|
# pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/slope/tiles/$filename.tif -o $OUTDIR/slope/median/tiles/$filename.tif -co COMPRESS=LZW -ot Byte
|
41
|
# # slope stdev
|
42
|
# pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/slope/tiles/$filename.tif -o $OUTDIR/slope/stdev/tiles/tmp_$filename.tif -co COMPRESS=LZW -ot Int32
|
43
|
# gdal_calc.py -A $OUTDIR/slope/stdev/tiles/tmp_$filename.tif --calc="sqrt(A)" --type Int16 --overwrite --outfile $OUTDIR/slope/stdev/tiles/$filename.tif
|
44
|
# rm -f $OUTDIR/slope/stdev/tiles/tmp_$filename
|
45
|
# # slope min
|
46
|
# pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/slope/tiles/$filename.tif -o $OUTDIR/slope/min/tiles/$filename.tif -co COMPRESS=LZW -ot Byte
|
47
|
# # slope max
|
48
|
# pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/slope/tiles/$filename.tif -o $OUTDIR/slope/max/tiles/$filename.tif -co COMPRESS=LZW -ot Byte
|
49
|
# # slope mean
|
50
|
# pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/slope/tiles/$filename.tif -o $OUTDIR/slope/mean/tiles/$filename.tif -co COMPRESS=LZW -ot Byte
|
51
|
|
52
|
# rm -f $OUTDIR/slope/tiles/$filename.tif
|
53
|
|
54
|
# echo generate a Terrain Ruggedness Index TRI with file $file
|
55
|
# gdaldem TRI -co COMPRESS=LZW $OUTDIR/tif/$filename.tif $OUTDIR/tri/tiles/$filename.tif
|
56
|
# # tri median
|
57
|
# pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/tri/tiles/$filename.tif -o $OUTDIR/tri/median/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
58
|
# # tri stdev
|
59
|
# pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/tri/tiles/$filename.tif -o $OUTDIR/tri/stdev/tiles/tmp_$filename.tif -co COMPRESS=LZW
|
60
|
# gdal_calc.py -A $OUTDIR/tri/stdev/tiles/tmp_$filename.tif --calc="sqrt(A)" --type Int16 --overwrite --outfile $OUTDIR/tri/stdev/tiles/$filename.tif
|
61
|
# rm -f $OUTDIR/tri/stdev/tiles/tmp_$filename
|
62
|
# # tri min
|
63
|
# pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/tri/tiles/$filename.tif -o $OUTDIR/tri/min/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
64
|
# # tri max
|
65
|
# pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/tri/tiles/$filename.tif -o $OUTDIR/tri/max/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
66
|
# # tri mean
|
67
|
# pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/tri/tiles/$filename.tif -o $OUTDIR/tri/mean/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
68
|
|
69
|
# rm -f $OUTDIR/tri/tiles/$filename.tif
|
70
|
|
71
|
|
72
|
echo generate a Topographic Position Index TPI with file $OUTDIR/tif/$filename.tif
|
73
|
# gdaldem TPI -co COMPRESS=LZW $OUTDIR/tif/$filename.tif $OUTDIR/tpi/tiles/$filename.tif # tpi has negative number
|
74
|
|
75
|
oft-calc -ot Float32 $OUTDIR/tpi/tiles/$filename.tif $OUTDIR/tpi/tiles/$filename"_t10.tif" &> /dev/null <<EOF
|
76
|
1
|
77
|
#1 10 *
|
78
|
EOF
|
79
|
|
80
|
# tpi median
|
81
|
pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/tpi/tiles/$filename"_t10.tif" -o $OUTDIR/tpi/median/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
82
|
echo tpi stdev
|
83
|
pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/tpi/tiles/$filename"_t10.tif" -o $OUTDIR/tpi/stdev/tiles/tmp_$filename.tif -co COMPRESS=LZW
|
84
|
gdal_calc.py -A $OUTDIR/tpi/stdev/tiles/tmp_$filename.tif --calc="sqrt(A)" --type Int32 --overwrite --outfile $OUTDIR/tpi/stdev/tiles/$filename.tif
|
85
|
rm -f $OUTDIR/tpi/stdev/tiles/tmp_$filename
|
86
|
# tpi min
|
87
|
pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/tpi/tiles/$filename"_t10.tif" -o $OUTDIR/tpi/min/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
88
|
# tpi max
|
89
|
pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/tpi/tiles/$filename"_t10.tif" -o $OUTDIR/tpi/max/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
90
|
# tpi mean
|
91
|
pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/tpi/tiles/$filename"_t10.tif" -o $OUTDIR/tpi/mean/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
92
|
|
93
|
|
94
|
|
95
|
|
96
|
|
97
|
# echo generate roughness with file $OUTDIR/tif/$filename.tif
|
98
|
|
99
|
# gdaldem roughness -co COMPRESS=LZW $OUTDIR/tif/$filename.tif $OUTDIR/roughness/tiles/$filename.tif
|
100
|
|
101
|
# # roughness median
|
102
|
# pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/roughness/tiles/$filename.tif -o $OUTDIR/roughness/median/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
103
|
# echo roughness stdev
|
104
|
# pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/roughness/tiles/$filename.tif -o $OUTDIR/roughness/stdev/tiles/tmp_$filename.tif -co COMPRESS=LZW
|
105
|
# gdal_calc.py -A $OUTDIR/roughness/stdev/tiles/tmp_$filename.tif --calc="sqrt(A)" --type Int32 --overwrite --outfile $OUTDIR/roughness/stdev/tiles/$filename.tif
|
106
|
|
107
|
# rm -f $OUTDIR/roughness/stdev/tiles/tmp_$filename
|
108
|
|
109
|
# # roughness min
|
110
|
# pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/roughness/tiles/$filename.tif -o $OUTDIR/roughness/min/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
111
|
# # roughness max
|
112
|
# pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/roughness/tiles/$filename.tif -o $OUTDIR/roughness/max/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
113
|
# # roughness mean
|
114
|
# pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/roughness/tiles/$filename.tif -o $OUTDIR/roughness/mean/tiles/$filename.tif -co COMPRESS=LZW -ot Int16
|
115
|
|
116
|
# rm -f $OUTDIR/roughness/tiles/$filename.tif
|
117
|
|
118
|
|
119
|
|
120
|
echo aspect with file $OUTDIR/tif/$filename.tif
|
121
|
|
122
|
gdaldem aspect -zero_for_flat -co COMPRESS=LZW $OUTDIR/tif/$filename.tif $OUTDIR/aspect/tiles/$filename.tif
|
123
|
|
124
|
# r1 aspect , r2 slope
|
125
|
|
126
|
gdal_calc.py --NoDataValue -9999 -A $OUTDIR/aspect/tiles/$filename.tif --calc="(sin(A))" --outfile $OUTDIR/aspect/tiles/$filename"_sin.tif" --overwrite --type Float32
|
127
|
gdal_calc.py --NoDataValue -9999 -A $OUTDIR/aspect/tiles/$filename.tif --calc="(cos(A))" --outfile $OUTDIR/aspect/tiles/$filename"_cos.tif" --overwrite --type Float32
|
128
|
|
129
|
gdal_calc.py --NoDataValue -9999 -A $OUTDIR/slope/tiles/$filename.tif -B $OUTDIR/aspect/tiles/$filename"_sin.tif" --calc="((sin(A))*B)" --outfile $OUTDIR/aspect/tiles/$filename"_Ew.tif" --overwrite --type Float32
|
130
|
gdal_calc.py --NoDataValue -9999 -A $OUTDIR/slope/tiles/$filename.tif -B $OUTDIR/aspect/tiles/$filename"_cos.tif" --calc="((sin(A))*B)" --outfile $OUTDIR/aspect/tiles/$filename"_Nw.tif" --overwrite --type Float32
|
131
|
|
132
|
echo aspect sin cos Ew Nw median
|
133
|
# sin
|
134
|
pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/aspect/tiles/$filename"_sin.tif" -o $OUTDIR/aspect/median/tiles/$filename"_sin_f.tif" -co COMPRESS=LZW -ot Float32
|
135
|
oft-calc -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_sin_f.tif" $OUTDIR/aspect/median/tiles/$filename"_sin_t10k_tmp.tif" &> /dev/null <<EOF
|
136
|
1
|
137
|
#1 10000 *
|
138
|
EOF
|
139
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_sin_t10k_tmp.tif" $OUTDIR/aspect/median/tiles/$filename"_sin_t10k.tif"
|
140
|
rm $OUTDIR/aspect/median/tiles/$filename"_sin_t10k_tmp.tif"
|
141
|
# cos
|
142
|
pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/aspect/tiles/$filename"_cos.tif" -o $OUTDIR/aspect/median/tiles/$filename"_cos_f.tif" -co COMPRESS=LZW -ot Float32
|
143
|
oft-calc -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_cos_f.tif" $OUTDIR/aspect/median/tiles/$filename"_cos_t10k_tmp.tif" &> /dev/null <<EOF
|
144
|
1
|
145
|
#1 10000 *
|
146
|
EOF
|
147
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_cos_t10k_tmp.tif" $OUTDIR/aspect/median/tiles/$filename"_cos_t10k.tif"
|
148
|
rm $OUTDIR/aspect/median/tiles/$filename"_cos_t10k_tmp.tif"
|
149
|
# Ew
|
150
|
pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/aspect/tiles/$filename"_Ew.tif" -o $OUTDIR/aspect/median/tiles/$filename"_Ew_f.tif" -co COMPRESS=LZW -ot Float32
|
151
|
oft-calc -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_Ew_f.tif" $OUTDIR/aspect/median/tiles/$filename"_Ew_t10k_tmp.tif" &> /dev/null <<EOF
|
152
|
1
|
153
|
#1 10000 *
|
154
|
EOF
|
155
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_Ew_t10k_tmp.tif" $OUTDIR/aspect/median/tiles/$filename"_Ew_t10k.tif"
|
156
|
rm $OUTDIR/aspect/median/tiles/$filename"_Ew_t10k_tmp.tif"
|
157
|
# Nw
|
158
|
pkfilter -m -9999 -dx 10 -dy 10 -f median -d 10 -i $OUTDIR/aspect/tiles/$filename"_Nw.tif" -o $OUTDIR/aspect/median/tiles/$filename"_Nw_f.tif" -co COMPRESS=LZW -ot Float32
|
159
|
oft-calc -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_Nw_f.tif" $OUTDIR/aspect/median/tiles/$filename"_Nw_t10k_tmp.tif" &> /dev/null <<EOF
|
160
|
1
|
161
|
#1 10000 *
|
162
|
EOF
|
163
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/median/tiles/$filename"_Nw_t10k_tmp.tif" $OUTDIR/aspect/median/tiles/$filename"_Nw_t10k.tif"
|
164
|
rm $OUTDIR/aspect/median/tiles/$filename"_Nw_t10k_tmp.tif"
|
165
|
|
166
|
|
167
|
echo aspect sin cos Ew Nw mean
|
168
|
# sin
|
169
|
pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/aspect/tiles/$filename"_sin.tif" -o $OUTDIR/aspect/mean/tiles/$filename"_sin_f.tif" -co COMPRESS=LZW -ot Float32
|
170
|
oft-calc -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_sin_f.tif" $OUTDIR/aspect/mean/tiles/$filename"_sin_t10k_tmp.tif" &> /dev/null <<EOF
|
171
|
1
|
172
|
#1 10000 *
|
173
|
EOF
|
174
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_sin_t10k_tmp.tif" $OUTDIR/aspect/mean/tiles/$filename"_sin_t10k.tif"
|
175
|
rm $OUTDIR/aspect/mean/tiles/$filename"_sin_t10k_tmp.tif"
|
176
|
# cos
|
177
|
pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/aspect/tiles/$filename"_cos.tif" -o $OUTDIR/aspect/mean/tiles/$filename"_cos_f.tif" -co COMPRESS=LZW -ot Float32
|
178
|
oft-calc -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_cos_f.tif" $OUTDIR/aspect/mean/tiles/$filename"_cos_t10k_tmp.tif" &> /dev/null <<EOF
|
179
|
1
|
180
|
#1 10000 *
|
181
|
EOF
|
182
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_cos_t10k_tmp.tif" $OUTDIR/aspect/mean/tiles/$filename"_cos_t10k.tif"
|
183
|
rm $OUTDIR/aspect/mean/tiles/$filename"_cos_t10k_tmp.tif"
|
184
|
# Ew
|
185
|
pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/aspect/tiles/$filename"_Ew.tif" -o $OUTDIR/aspect/mean/tiles/$filename"_Ew_f.tif" -co COMPRESS=LZW -ot Float32
|
186
|
oft-calc -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_Ew_f.tif" $OUTDIR/aspect/mean/tiles/$filename"_Ew_t10k_tmp.tif" &> /dev/null <<EOF
|
187
|
1
|
188
|
#1 10000 *
|
189
|
EOF
|
190
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_Ew_t10k_tmp.tif" $OUTDIR/aspect/mean/tiles/$filename"_Ew_t10k.tif"
|
191
|
rm $OUTDIR/aspect/mean/tiles/$filename"_Ew_t10k_tmp.tif"
|
192
|
# Nw
|
193
|
pkfilter -m -9999 -dx 10 -dy 10 -f mean -d 10 -i $OUTDIR/aspect/tiles/$filename"_Nw.tif" -o $OUTDIR/aspect/mean/tiles/$filename"_Nw_f.tif" -co COMPRESS=LZW -ot Float32
|
194
|
oft-calc -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_Nw_f.tif" $OUTDIR/aspect/mean/tiles/$filename"_Nw_t10k_tmp.tif" &> /dev/null <<EOF
|
195
|
1
|
196
|
#1 10000 *
|
197
|
EOF
|
198
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/mean/tiles/$filename"_Nw_t10k_tmp.tif" $OUTDIR/aspect/mean/tiles/$filename"_Nw_t10k.tif"
|
199
|
rm $OUTDIR/aspect/mean/tiles/$filename"_Nw_t10k_tmp.tif"
|
200
|
|
201
|
echo aspect sin cos Ew Nw max
|
202
|
# sin
|
203
|
pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/aspect/tiles/$filename"_sin.tif" -o $OUTDIR/aspect/max/tiles/$filename"_sin_f.tif" -co COMPRESS=LZW -ot Float32
|
204
|
oft-calc -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_sin_f.tif" $OUTDIR/aspect/max/tiles/$filename"_sin_t10k_tmp.tif" &> /dev/null <<EOF
|
205
|
1
|
206
|
#1 10000 *
|
207
|
EOF
|
208
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_sin_t10k_tmp.tif" $OUTDIR/aspect/max/tiles/$filename"_sin_t10k.tif"
|
209
|
rm $OUTDIR/aspect/max/tiles/$filename"_sin_t10k_tmp.tif"
|
210
|
# cos
|
211
|
pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/aspect/tiles/$filename"_cos.tif" -o $OUTDIR/aspect/max/tiles/$filename"_cos_f.tif" -co COMPRESS=LZW -ot Float32
|
212
|
oft-calc -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_cos_f.tif" $OUTDIR/aspect/max/tiles/$filename"_cos_t10k_tmp.tif" &> /dev/null <<EOF
|
213
|
1
|
214
|
#1 10000 *
|
215
|
EOF
|
216
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_cos_t10k_tmp.tif" $OUTDIR/aspect/max/tiles/$filename"_cos_t10k.tif"
|
217
|
rm $OUTDIR/aspect/max/tiles/$filename"_cos_t10k_tmp.tif"
|
218
|
# Ew
|
219
|
pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/aspect/tiles/$filename"_Ew.tif" -o $OUTDIR/aspect/max/tiles/$filename"_Ew_f.tif" -co COMPRESS=LZW -ot Float32
|
220
|
oft-calc -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_Ew_f.tif" $OUTDIR/aspect/max/tiles/$filename"_Ew_t10k_tmp.tif" &> /dev/null <<EOF
|
221
|
1
|
222
|
#1 10000 *
|
223
|
EOF
|
224
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_Ew_t10k_tmp.tif" $OUTDIR/aspect/max/tiles/$filename"_Ew_t10k.tif"
|
225
|
rm $OUTDIR/aspect/max/tiles/$filename"_Ew_t10k_tmp.tif"
|
226
|
# Nw
|
227
|
pkfilter -m -9999 -dx 10 -dy 10 -f max -d 10 -i $OUTDIR/aspect/tiles/$filename"_Nw.tif" -o $OUTDIR/aspect/max/tiles/$filename"_Nw_f.tif" -co COMPRESS=LZW -ot Float32
|
228
|
oft-calc -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_Nw_f.tif" $OUTDIR/aspect/max/tiles/$filename"_Nw_t10k_tmp.tif" &> /dev/null <<EOF
|
229
|
1
|
230
|
#1 10000 *
|
231
|
EOF
|
232
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/max/tiles/$filename"_Nw_t10k_tmp.tif" $OUTDIR/aspect/max/tiles/$filename"_Nw_t10k.tif"
|
233
|
rm $OUTDIR/aspect/max/tiles/$filename"_Nw_t10k_tmp.tif"
|
234
|
|
235
|
echo aspect sin cos Ew Nw min
|
236
|
# sin
|
237
|
pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/aspect/tiles/$filename"_sin.tif" -o $OUTDIR/aspect/min/tiles/$filename"_sin_f.tif" -co COMPRESS=LZW -ot Float32
|
238
|
oft-calc -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_sin_f.tif" $OUTDIR/aspect/min/tiles/$filename"_sin_t10k_tmp.tif" &> /dev/null <<EOF
|
239
|
1
|
240
|
#1 10000 *
|
241
|
EOF
|
242
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_sin_t10k_tmp.tif" $OUTDIR/aspect/min/tiles/$filename"_sin_t10k.tif"
|
243
|
rm $OUTDIR/aspect/min/tiles/$filename"_sin_t10k_tmp.tif"
|
244
|
# cos
|
245
|
pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/aspect/tiles/$filename"_cos.tif" -o $OUTDIR/aspect/min/tiles/$filename"_cos_f.tif" -co COMPRESS=LZW -ot Float32
|
246
|
oft-calc -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_cos_f.tif" $OUTDIR/aspect/min/tiles/$filename"_cos_t10k_tmp.tif" &> /dev/null <<EOF
|
247
|
1
|
248
|
#1 10000 *
|
249
|
EOF
|
250
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_cos_t10k_tmp.tif" $OUTDIR/aspect/min/tiles/$filename"_cos_t10k.tif"
|
251
|
rm $OUTDIR/aspect/min/tiles/$filename"_cos_t10k_tmp.tif"
|
252
|
# Ew
|
253
|
pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/aspect/tiles/$filename"_Ew.tif" -o $OUTDIR/aspect/min/tiles/$filename"_Ew_f.tif" -co COMPRESS=LZW -ot Float32
|
254
|
oft-calc -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_Ew_f.tif" $OUTDIR/aspect/min/tiles/$filename"_Ew_t10k_tmp.tif" &> /dev/null <<EOF
|
255
|
1
|
256
|
#1 10000 *
|
257
|
EOF
|
258
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_Ew_t10k_tmp.tif" $OUTDIR/aspect/min/tiles/$filename"_Ew_t10k.tif"
|
259
|
rm $OUTDIR/aspect/min/tiles/$filename"_Ew_t10k_tmp.tif"
|
260
|
# Nw
|
261
|
pkfilter -m -9999 -dx 10 -dy 10 -f min -d 10 -i $OUTDIR/aspect/tiles/$filename"_Nw.tif" -o $OUTDIR/aspect/min/tiles/$filename"_Nw_f.tif" -co COMPRESS=LZW -ot Float32
|
262
|
oft-calc -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_Nw_f.tif" $OUTDIR/aspect/min/tiles/$filename"_Nw_t10k_tmp.tif" &> /dev/null <<EOF
|
263
|
1
|
264
|
#1 10000 *
|
265
|
EOF
|
266
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/min/tiles/$filename"_Nw_t10k_tmp.tif" $OUTDIR/aspect/min/tiles/$filename"_Nw_t10k.tif"
|
267
|
rm $OUTDIR/aspect/min/tiles/$filename"_Nw_t10k_tmp.tif"
|
268
|
|
269
|
echo aspect sin cos Ew Nw stdev
|
270
|
# stdev
|
271
|
pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/aspect/tiles/$filename"_sin.tif" -o $OUTDIR/aspect/stdev/tiles/$filename"_sin_f.tif" -co COMPRESS=LZW -ot Float32
|
272
|
gdal_calc.py -A $OUTDIR/aspect/stdev/tiles/$filename"_sin_f.tif" --calc="(sqrt(A))*10000" --type Int16 --overwrite --outfile $OUTDIR/aspect/stdev/tiles/$filename"_sin_t10k_tmp.tif"
|
273
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/stdev/tiles/$filename"_sin_t10k_tmp.tif" $OUTDIR/aspect/stdev/tiles/$filename"_sin_t10k.tif"
|
274
|
rm $OUTDIR/aspect/stdev/tiles/$filename"_sin_t10k_tmp.tif"
|
275
|
# cos
|
276
|
pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/aspect/tiles/$filename"_cos.tif" -o $OUTDIR/aspect/stdev/tiles/$filename"_cos_f.tif" -co COMPRESS=LZW -ot Float32
|
277
|
gdal_calc.py -A $OUTDIR/aspect/stdev/tiles/$filename"_cos_f.tif" --calc="(sqrt(A))*10000" --type Int16 --overwrite --outfile $OUTDIR/aspect/stdev/tiles/$filename"_cos_t10k_tmp.tif"
|
278
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/stdev/tiles/$filename"_cos_t10k_tmp.tif" $OUTDIR/aspect/stdev/tiles/$filename"_cos_t10k.tif"
|
279
|
rm $OUTDIR/aspect/stdev/tiles/$filename"_cos_t10k_tmp.tif"
|
280
|
# Ew
|
281
|
pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/aspect/tiles/$filename"_Ew.tif" -o $OUTDIR/aspect/stdev/tiles/$filename"_Ew_f.tif" -co COMPRESS=LZW -ot Float32
|
282
|
gdal_calc.py -A $OUTDIR/aspect/stdev/tiles/$filename"_Ew_f.tif" --calc="(sqrt(A))*10000" --type Int16 --overwrite --outfile $OUTDIR/aspect/stdev/tiles/$filename"_Ew_t10k_tmp.tif"
|
283
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/stdev/tiles/$filename"_Ew_t10k_tmp.tif" $OUTDIR/aspect/stdev/tiles/$filename"_Ew_t10k.tif"
|
284
|
rm $OUTDIR/aspect/stdev/tiles/$filename"_Ew_t10k_tmp.tif"
|
285
|
# Nw
|
286
|
pkfilter -m -9999 -dx 10 -dy 10 -f var -d 10 -i $OUTDIR/aspect/tiles/$filename"_Nw.tif" -o $OUTDIR/aspect/stdev/tiles/$filename"_Nw_f.tif" -co COMPRESS=LZW -ot Float32
|
287
|
gdal_calc.py -A $OUTDIR/aspect/stdev/tiles/$filename"_Nw_f.tif" --calc="(sqrt(A))*10000" --type Int16 --overwrite --outfile $OUTDIR/aspect/stdev/tiles/$filename"_Nw_t10k_tmp.tif"
|
288
|
gdal_translate -co COMPRESS=LZW -ot Int16 $OUTDIR/aspect/stdev/tiles/$filename"_Nw_t10k_tmp.tif" $OUTDIR/aspect/stdev/tiles/$filename"_Nw_t10k.tif"
|
289
|
rm $OUTDIR/aspect/stdev/tiles/$filename"_Nw_t10k_tmp.tif"
|
290
|
|
291
|
) 2>&1 | tee /mnt/data2/dem_variables/log.txt
|