1
|
# 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
|
|
76
|
|