# calculate altitude area percent above a specific treshold

# e.g  75% of class200 means that 75% of the pixels is >= 200 

# for this case the treshold is set every 100 meter

# geting the integer of the altitued / 100

#  for file in /mnt/data2/dem_variables/GMTED2010/tiles/mi75_grd_tif/*.tif ; do echo $file mi ; done  | xargs -n 2 -P 25  bash /mnt/data2/dem_variables/GMTED2010/scripts/sc2_class_treshold_percent.sh 

# rm  /tmp/ramdisk/* ; for file in /mnt/data2/dem_variables/GMTED2010/tiles/mx75_grd_tif/*.tif ; do echo $file mx ; done  | xargs -n 2 -P 25  bash /mnt/data2/dem_variables/GMTED2010/scripts/sc2_class_treshold_percent.sh 

# to check 

# ny=45 ; for n in `seq 1 9` ;do echo `gdallocationinfo -valonly altitude/class/5_1_class_tmp3.tif 2$n $ny` `gdallocationinfo  -valonly  tiles/mn75_grd_tif/5_1.tif 2$n $ny` ; done 

export file=$1

export mm=$2

export filename=`basename $file .tif`

(

INDIR=/mnt/data2/dem_variables/GMTED2010/tiles/$mm"75_grd_tif"

OUTDIR=/mnt/data2/dem_variables/GMTED2010/altitude/class_$mm/class

echo  calcualte altitude  every 100 meter  using the the Int16 factors for $file 

oft-calc  -X -ot Int16 $INDIR/$filename.tif $OUTDIR/../tmp/$filename"_class_tmp1.tif"  &> /dev/null   <<EOF

1

#1 49 - 100 /

EOF

oft-calc  -X -ot Int16  $OUTDIR/../tmp/$filename"_class_tmp1.tif"  $OUTDIR/../tmp/$filename"_class_tmp2.tif"    &> /dev/null   <<EOF

1

#1 100 *

EOF

rm $OUTDIR/../tmp/$filename"_class_tmp1.tif" 

gdal_translate -co COMPRESS=LZW -co ZLEVEL=9  $OUTDIR/../tmp/$filename"_class_tmp2.tif" $OUTDIR/$filename"_class.tif"

rm $OUTDIR/../tmp/$filename"_class_tmp2.tif"

gdalinfo -mm $OUTDIR/$filename"_class.tif"  | grep Computed | awk '{ gsub ("[=,]"," "); print int($(NF-1)), int($(NF))}' > $OUTDIR/../tmp/$filename"_min_max_"$mm.txt

echo min and max for $OUTDIR/../tmp/$filename"_min_max_"$mm.txt

for class in $(seq `awk '{ print $1}'  $OUTDIR/../tmp/${filename}_min_max_${mm}.txt` 100  `awk '{ print $2}' $OUTDIR/../tmp/${filename}_min_max_${mm}.txt`) ; do 

    class_start=$class

    class_end=$(awk '{ print $2}'  $OUTDIR/../tmp/${filename}_min_max_${mm}.txt)

    class_string=$(for class_unique in $(seq $class_start 100  $class_end  ) ; do echo -n "-class $class_unique "  ; done) 

    echo applaid the filter for the classes  $class_string

    pkfilter -dx 4 -dy 4 $class_string -f density -d 4 -i $OUTDIR/$filename"_class.tif" -o $OUTDIR/../class$class/$filename"_C"$class"Perc.tif" -co COMPRESS=LZW -co ZLEVEL=9 -ot Byte 

done 

) 2>&1 | tee  /mnt/data2/dem_variables/GMTED2010/log/log_of_class_$mm"_"$filename".txt"