1 |
73e9d9c6
|
Giuseppe
|
# calculate different variables for the dem
|
2 |
|
|
# for dir1 in altitude ; do for dir2 in stdev ; do echo $dir1/$dir2 ; done ; done | xargs -n 1 -P 5 bash /mnt/data2/dem_variables/scripts/sc2_dem_merge.sh
|
3 |
|
|
# for dir1 in altitude slope tri tpi roughness ; do for dir2 in max mean median min stdev; do echo $dir1/$dir2 ; done ; done | xargs -n 1 -P 12 bash /mnt/data2/dem_variables/scripts/sc2_dem_merge.sh
|
4 |
|
|
|
5 |
|
|
|
6 |
|
|
OUTDIR=/mnt/data2/dem_variables/resol_x10
|
7 |
|
|
|
8 |
|
|
DIR=$1
|
9 |
|
|
dir1=$(echo ${DIR%/*}) # cancel the part after /
|
10 |
|
|
dir2=$(echo ${DIR#*/}) # cancel the part before /
|
11 |
|
|
|
12 |
|
|
|
13 |
|
|
if [ $dir1 = altitude ] ; then type=Int16 ; fi
|
14 |
|
|
if [ $dir1 = aspect ] ; then type=Int16 ; fi
|
15 |
|
|
if [ $dir1 = slope ] ; then type=Byte ; fi
|
16 |
|
|
if [ $dir1 = tri ] ; then type=Int16 ; fi
|
17 |
|
|
if [ $dir1 = tpi ] ; then type=Int16 ; fi
|
18 |
|
|
if [ $dir1 = roughness ] ; then type=Int16 ; fi
|
19 |
|
|
|
20 |
|
|
if [ $dir1 != aspect ]; then
|
21 |
|
|
|
22 |
|
|
echo processing merging tiles in $dir1
|
23 |
|
|
|
24 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2.tif
|
25 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NW".tif
|
26 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NE".tif
|
27 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SW".tif
|
28 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SE".tif
|
29 |
|
|
|
30 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_NW".tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_N*W*.tif
|
31 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_NE".tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_N*E*.tif
|
32 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_SW".tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_S*W*.tif
|
33 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_SE".tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_S*E*.tif
|
34 |
|
|
|
35 |
|
|
exit
|
36 |
|
|
|
37 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2.tif
|
38 |
|
|
|
39 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/$dir1"_"$dir2"_NW".tif $OUTDIR/$DIR/$dir1"_"$dir2"_NE".tif $OUTDIR/$DIR/$dir1"_"$dir2"_SW".tif $OUTDIR/$DIR/$dir1"_"$dir2"_SE".tif
|
40 |
|
|
|
41 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NW".tif
|
42 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NE".tif
|
43 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SW".tif
|
44 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SE".tif
|
45 |
|
|
|
46 |
|
|
else # just applaied to the aspect variables
|
47 |
|
|
|
48 |
|
|
echo processing merging tiles for the different aspect variables in $dir1
|
49 |
|
|
|
50 |
|
|
for aspect_var in "_sin_t10k" "_cos_t10k" "_Ew_t10k" "_Nw_t10k" ; do
|
51 |
|
|
|
52 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NW"$aspect_var.tif
|
53 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NE"$aspect_var.tif
|
54 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SW"$aspect_var.tif
|
55 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SE"$aspect_var.tif
|
56 |
|
|
|
57 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_NW"$aspect_var.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_N*W*$aspect_var.tif
|
58 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_NE"$aspect_var.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_N*E*$aspect_var.tif
|
59 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_SW"$aspect_var.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_S*W*$aspect_var.tif
|
60 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2"_SE"$aspect_var.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/tiles/Smoothed_S*E*$aspect_var.tif
|
61 |
|
|
|
62 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2$aspect_var.tif
|
63 |
|
|
|
64 |
|
|
gdal_merge.py -o $OUTDIR/$DIR/$dir1"_"$dir2$aspect_var.tif -ot $type -co COMPRESS=LZW $OUTDIR/$DIR/$dir1"_"$dir2"_NW"$aspect_var.tif $OUTDIR/$DIR/$dir1"_"$dir2"_NE"$aspect_var.tif $OUTDIR/$DIR/$dir1"_"$dir2"_SW"$aspect_var.tif $OUTDIR/$DIR/$dir1"_"$dir2"_SE"$aspect_var.tif
|
65 |
|
|
|
66 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NW"$aspect_var.tif
|
67 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_NE"$aspect_var.tif
|
68 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SW"$aspect_var.tif
|
69 |
|
|
rm -f $OUTDIR/$DIR/$dir1"_"$dir2"_SE"$aspect_var.tif
|
70 |
|
|
|
71 |
|
|
done
|
72 |
|
|
|
73 |
|
|
fi
|
74 |
|
|
|
75 |
|
|
|