1 |
73e9d9c6
|
Giuseppe
|
# 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
|