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