#!/bin/bash
printf "Content-type: text/html\n\n"
#QUERY_STRING='AlblasserdamXDordrechtXGorinchemXHardinxveld-GiessendamXHendrik-Ido-AmbachtXHoekscheWaardXMolenlandenXPapendrechtXSliedrechtXZwijndrecht'
#QUERY_STRING='AmsterdamXRotterdamXs-GravenhageXUtrechtXEindhoven'
#QUERY_STRING='gem=VRZuid-Holland-Zuid'
#QUERY_STRING='gem=Alblasserdam'
#echo $QUERY_STRING
STRING=$( echo  "$QUERY_STRING" | sed 's/%28/(/g;s/%29/)/g;s/%3A/:/g' | sed 's/[^a-zA-Z0-9&:=._-]/Q/g' )
#echo $STRING
IFS="&"  read -a ARRAY <<< "$STRING"
str=${ARRAY[0]}
gem=${str#gem=}
#gem=${gem:-Land_van_Cuijk}
gem=${gem:-Papendrecht}
#gem=${gem/%./}
#echo $gem
#exit
if [[ -z "${gem##VR*}" ]]
then
	gem=$(grep $gem < /home/pi/corona/regios | awk '{print $2}')
fi
gem=${gem//X/\\\>\\\|\\\<}
#cat /home/pi/gem$(date +%Y%m%d) | awk '{FS=";";if($3=="'${gem}'"){print $1,$5,$6,$7}' | tr ";" " " | sed 's/ /T/' | awk '{w=$0;str=$1;sub("T"," ",str);  cmd = "date +\"%Y %j\" -d \""str"\""; if (cmd | getline){print $0,w } close(cmd); }' > /tmp/gem$$
#echo $STRING "<br>"$gem
echo "<html> <head> <meta charset=-8> </head>"
if [ -e /home/pi/corona/gem$(date +%Y%m%d) ]
then
	date
	off=0
else
	date -d "-1day" +%Y%m%d
	off=1
fi


echo "<br>"

if [[ "${gem}" == "${gem/|/Q}" ]]
then
cat /home/pi/corona/gem$(date +%Y%m%d -d -${off}day) | sed 's/ [^;]*//' | cut -f1,3,4,6,7,8 -d";"  | sed 's/ /_/g;s/(\|)//g' | grep "[;']${gem};" | tr ";" " " | awk '{w=$0;str=$1;sub("T"," ",str);  cmd = "date +\"%Y %j\" -d \""str"\""; if (cmd | getline){print $1,$2+366*($1-2020),w } close(cmd); }' | sort -k4,4 -k2,2g | awk '(q!=$4){print "\n"; w++;}{print $0,w;q=$4}' > /tmp/gem$$
	tail /tmp/gem$$ | sed 's/$/<br>/'
echo 'set timestamp "%Y-%m-%dT%H:%M:%S"  font "DejaVuSansCondensed,8"; set key left; set term png size 1600,900; set out "/var/www/temp/img/gem.png";  set log y; set grid; set title "'"${gem}"'" ;plot "/tmp/gem'$$'" u 2:5  t "besmettingen", "" u 2:6 t "ziekenhuis", "" u 2:7 t "dood" w l;'  | gnuplot 
cat /home/pi/corgemdens | grep  "[' ]${gem} "   | sort -k2,2 -k1,1g | awk '{a[1000+$1]=$4;$4=$4-a[1000+$1-7];$5=a[1000+$1];$6=a[1000+$1]-a[1000+$1-1];print $0}'  | awk '(q!=$2){print "\n"; w++;}{$4=$4/$7*1000;print $0,w;q=$2}' > /tmp/cgd 
echo '<table>'
tail /tmp/cgd | sed 's/^/<tr><td>/;s/ /<\/td><td>/g;s/$/<\/td><\/tr>/'
echo '</table>'
cb=$(awk '{print $2,$9}' < /tmp/cgd | uniq  | awk '{printf "\"%s\" %d ,", $1,$2}' | sed 's/,$//' | tr "_" " " ); 
if [ $gem == "Alblasserdam" ]
then
echo "set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8'; set palette defined ( 0 0 1 0, 0.5 1 1 1, 1 1 0 0 );set key left; set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [:];set grid;set multi; set title '${gem}'; plot [250:]'/tmp/cgd' u 1:4  t 'actieve gevallen per 1000 (laatste 7 dagen)' w l lw 3 , '' u ((\$1==693)?\$1:NaN):4w p pt 7 ps 5 lc 5 t 'gele helm',  '' u 1:(\$6/\$7*7000):(\$4<\$6/\$7*7000) w i pal  lt -1 lw 2 t'dagelijks*7'; " | gnuplot 2>&1 > /var/www/temp/img/gnuout
else
#echo "set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8'; set palette defined ( 0 0 1 0, 0.5 1 1 1, 1 1 0 0 );set key left; set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [:];set grid; set title '${gem}'; set multiplot; set size 1,1; set origin 0,0; set xrange [250:];plot '/tmp/cgd' u 1:4  t 'actieve gevallen per 1000 (laatste 7 dagen)' w l lw 3 , '' u 1:(\$6/\$7*7000):(\$4<\$6/\$7*7000) w i pal  lt -1 lw 2 t'dagelijks*7'; set size 0.5,0,5; set origin 0.2,0.4; rep; uns multi  " | gnuplot 2>&1 > /var/www/temp/img/gnuout
echo "set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8';set xtics rotate ; set palette defined ( 0 0 1 0, 0.5 1 1 1, 1 1 0 0 ); set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [:];set grid; set title '${gem}';  set cbtics (); set key left; set xtics 14;set yrange [:];set grid; set title ''; set multiplot;set size 1,1; set origin 0,0; set xrange [250:];plot '/tmp/cgd' u 1:4  t 'actieve gevallen per 1000 (laatste 7 dagen)' w l lw 3 lc -1 , '' u 1:(\$6/\$7*7000):(\$4<\$6/\$7*7000) w i pal  lt -1 lw 1 t'dagelijks*7'; set size 0.5,0.5 ; set origin 0.1,0.4 ; set object 1 rectangle from screen 0.1,0.4to screen 0.6,0.9 fillcolor rgb '#ffffff' behind ; uns timestamp; set xtics 7; set xrange [GPVAL_X_MAX-100:GPVAL_X_MAX+2] ; uns key; uns colorbox ; rep ; uns multi ;"  | gnuplot 2>&1 > /var/www/temp/img/gnuout

fi
else
cat /home/pi/corona/gem$(date +%Y%m%d -d -${off}day) | sed 's/ [^;]*//' | cut -f1,3,4,6,7,8 -d";"  | sed 's/ /_/g;s/(\|)//g' | grep "\<${gem}\>" | tr ";" " " | awk '{w=$0;str=$1;sub("T"," ",str);  cmd = "date +\"%Y %j\" -d \""str"\""; if (cmd | getline){print $1,$2+366*($1-2020),w } close(cmd); }' | sort -k4,4 -k2,2g | awk '(q!=$4){print "\n"; w++;}{print $0,w;q=$4}' > /tmp/gem$$
	cb=$(awk '{print $4,$8}' < /tmp/gem$$ | uniq  | awk '{printf "\"%s\" %d ,", $1,$2}' | sed 's/,$//' | tr "_" " " )
echo 'set timestamp "%Y-%m-%dT%H:%M:%S"  font "DejaVuSansCondensed,8";set palette defined ( 0 0 0 0, 0.04 0 0 0 , 0.14 1 0 0 , 0.4 1 1 0 ,0.6 0.5 1 0, 0.85 0 0 1, 1 1 0 1 );set key left; set term png size 1600,900; set xtics 7; set xlabel "2020 (dag)";set cbtics ('$cb'); set out "/var/www/temp/img/gem.png";  set grid; set title "'"${gem}"'" ;plot [:]"/tmp/gem'$$'" u 2:5:8 w l lw 3 pal   t "besmettingen" ; '  | gnuplot 
#cat /home/pi/corgemdens | grep  "\<${gem}\>"   | sort -k2,2 -k1,1g | awk '{a[1000+$1]=$4;$4=$4-a[1000+$1-7];$5=a[1000+$1];$6=a[1000+$1]-a[1000+$1-1];print $0}'  | awk '(q!=$2){print "\n"; w++;}{$4=$4/$7*1000;print $0,w;q=$2}' > /tmp/cgd 
#cb=$(awk '{print $2,$9}' < /tmp/cgd | uniq  | awk '{printf "\"%s\" %d ,", $1,$2}' | sed 's/,$//' | tr "_" " " ); 
#echo "set colorsequence classic;set xlabel 'tijd (dag)';set ylabel 'actieve gevallen per 1000 (laatste 7 dagen)';set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8';set palette defined ( 0 0 0 0, 0.04 0 0 0 , 0.14 1 0 0 , 0.4 1 1 0 ,0.6 0.5 1 0, 0.85 0 0 1, 1 1 0 1 );set key left; set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [0:];set grid; plot [250:]'/tmp/cgd' u 1:4:9  not w l lw 3 pal , 75./200 lt -1 lc 1 lw 2 t '75 per 100000 per 14 dagen' , 2.5 lt -1 lw 2 lc 13  t '500 per 100000 in 14 dagen' " | gnuplot 2>&1 > /var/www/temp/img/gnuout
cat /home/pi/corgemdens | grep  "\<${gem}\>"   | sort -k2,2 -k1,1g | awk '{a[1000+$1]=$4;$4=$4-a[1000+$1-13];$5=a[1000+$1];$6=a[1000+$1]-a[1000+$1-1];print $0}'  | awk '(q!=$2){print "\n"; w++;}{$4=$4/$7*100000;print $0,w;q=$2}' > /tmp/cgd 
cb=$(awk '{print $2,$9}' < /tmp/cgd | uniq  | awk '{printf "\"%s\" %d ,", $1,$2}' | sed 's/,$//' | tr "_" " " ); 
echo "set colorsequence classic;set xlabel 'tijd (dag)';set ylabel 'actieve gevallen per 100000 (laatste 14 dagen)';set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8';set palette defined ( 0 0 0 0, 0.04 0 0 0 , 0.14 1 0 0 , 0.4 1 1 0 ,0.6 0.5 1 0, 0.85 0 0 1, 1 1 0 1 );set key left; set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [0:];set grid; plot [532:]'/tmp/cgd' u 1:4:9  not w l lw 3 pal , 75 lt -1 lc 1 lw 2 t '75 per 100000 per 14 dagen' , 500 lt -1 lw 2 lc 13  t '500 per 100000 in 14 dagen' " | gnuplot 2>&1 > /var/www/temp/img/gnuout
#echo "set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8';set palette defined ( 0 0 0 0, 0.04 0 0 0 , 0.14 1 0 0 , 0.4 1 1 0 ,0.6 0.5 1 0, 0.85 0 0 1, 1 1 0 1 );set key left; set term png size 1600,900; set out '/var/www/temp/img/gemdens.png'; set cbtics ($cb); set key left; set xtics 14;set yrange [0:];set grid;set multi; plot [250:]'/tmp/cgd' u 1:4:9  t 'actieve gevallen per 1000 (laatste 7 dagen)' w l lw 3 pal ; uns timestamp; set size 0.25,0.25; set origin 0.65,0.6;uns colorbox; ; set format y ''; set object 1 rectangle from screen 0.65,0.6 to screen 0.9,0.85 fillcolor rgb '#ffffff' behind ; set xtics rotate ;p '' u 1:4:9 not w l pal  ; uns multi; " | gnuplot 2>&1 > /var/www/temp/img/gnuout
fi

grep "${gem}" < /home/pi/corgemdens  | /home/pi/bin/hisdev -k 1 -c 4 -g| awk '{$1=$5;print $0}' | awk '{a[1000+$1]=$4;print $0,$4-a[1000+$1-7],a[1000+$1-7]-a[1000+$1-13]}' > /tmp/localR

echo "set timestamp '%Y-%m-%dT%H:%M:%S'  font 'DejaVuSansCondensed,8'; set term png size 1600,900; set out '/var/www/temp/img/localR.png'; set xtics 14;set grid; set xrange [210:]; set ytics 0.1;  set yrange [0:2]; plot [250:]'/tmp/localR' u 1:(log(\$6/\$7)+1) w lp t '$gem', '/home/pi/fullR' u (\$1+3.5):(log(\$5/\$4)+1) w l lw 3 lc 7 t 'Landelijke R' ; " | gnuplot 2>&1 > /var/www/temp/img/gnuout

echo "<img src=http://lus.petergeelhoed.nl/gemdens.png?$(date +%s)><br>"
echo "<img src=http://lus.petergeelhoed.nl/localR.png?$(date +%s)><br>"
#rm  /tmp/cgd
echo "<img src=http://lus.petergeelhoed.nl/gem.png?$(date +%s)><br>"
echo "<img src=http://lus.petergeelhoed.nl/gempoly.png?$(date +%s)>"
echo "<img src=http://lus.petergeelhoed.nl/gemrel.png?$(date +%s)>"
echo "<img src=http://lus.petergeelhoed.nl/gempoli.png?$(date +%s)>"
echo "<img src=http://lus.petergeelhoed.nl/gempole.png?$(date +%s)><br>"
echo "<img src=http://lus.petergeelhoed.nl/vrs.png?$(date +%s)><br>"
echo "</html>"
echo  /tmp/gem$$
#rm  /tmp/gem$$ 
