|
weatherstats |
Subversion Repositories: |
Compare with Previous - Blame - Download
#!/usr/bin/perl## Run this periodicaly by cron as root### WeatherStats## version: 0.6## author: Adam Pribyl, covex@lowlevel.cz# url: http://meteolinger.lowlevel.cz/# desc: perl scripts to save weather conditions from TOPCOM 265 NE Weather station into RRD tool chain# deps: wwsr, perl, rrdtools (best >= 1.0.42)# lic: GNU GPL# howto use: call this every 2 minutes from crontab as root#sub help {print "usage: \n";print " weatherstats -a <altitude> reads out the data from the weather station and stores them in round robin database\n";print " weatherstats -a <altitude> -g reads out the data and generates graphs\n";exit(0);}sub callsys {my @cmd = @_;# debug messagesif ( $en_dbg == 1 ) {$out = "";print "---- DBG: WeatherStats: Issuing command:\n @cmd\n";} else {$out = "> /dev/null";}system("@cmd $out");my $err = $?/256;if ( $err > 0 ) { die "ERROR: WeatherStats: System error nr. $err while executing @cmd" }}sub insertdb {my $rrdbn = $_[0];my $value = $_[1];my $scale = $_[2];my $wname = $_[3];my $value_type = $rrdbn;$tm = time;$min = "U";$max = "U"; #unlimitedif ($rrdbn =~ /temperature/) {$min = "-60";$max = "60";$value_type =~ s/temperature//g;$value_type =~ s/_//g;push(@{$temperature{$value_type}}, $wname);push(@{$temperature{$value_type}}, $value);push(@{$temperature{$value_type}}, $scale);$rrdbn = "temperature_".$value_type;} elsif ($rrdbn eq "wind_speed" or $rrdbn eq "wind_gust") {$min = "0";$max = "100"; #m/s > 330km/h$value_type =~ s/wind//g;$value_type =~ s/_//g;push(@{$wind{$value_type}}, $wname);push(@{$wind{$value_type}}, $value);push(@{$wind{$value_type}}, $scale);$rrdbn = "wind_".$value_type;} elsif ($rrdbn eq "wind_direction") {$min = "0";$max = "15";$value_type =~ s/wind//g;$value_type =~ s/_//g;push(@{$wind{$value_type}}, $wname);push(@{$wind{$value_type}}, $value);push(@{$wind{$value_type}}, $scale);$rrdbn = "wind_".$value_type;} elsif ($rrdbn =~ /humidity/) {$min = "0";$max = "100"; #%$value_type =~ s/humidity//g;$value_type =~ s/_//g;push(@{$humidity{$value_type}}, $wname);push(@{$humidity{$value_type}}, $value);push(@{$humidity{$value_type}}, $scale.$scale);$rrdbn = "humidity_".$value_type;} elsif ($rrdbn eq "rain_1h") {$min = "0";$max = "100"; #mm/h$value_type =~ s/rain//g;$value_type =~ s/_//g;push(@{$rain{$value_type}}, $wname);push(@{$rain{$value_type}}, $value);push(@{$rain{$value_type}}, $scale);$rrdbn = "rain_".$value_type;} elsif ($rrdbn eq "rain_total") {$value_type =~ s/rain//g;$value_type =~ s/_//g;push(@{$rain{$value_type}}, $wname);push(@{$rain{$value_type}}, $value);push(@{$rain{$value_type}}, $scale);$rrdbn = "rain_".$value_type;} elsif ($rrdbn =~ /pressure/) {$min = "800";$max = "1200";$value_type =~ s/pressure//g;$value_type =~ s/_//g;push(@{$pressure{$value_type}}, $wname);push(@{$pressure{$value_type}}, $value);push(@{$pressure{$value_type}}, $scale);$rrdbn = "pressure_".$value_type;}if ( !( -e $rrdbp.$rrdbn.".rrd") ) {# every 120s (2min), 600s maximum interval between feeding data# 1. AVERAGE, every 1st value, store 720 (30*24) values (1day)# 2. AVERAGE, every 30th value (each hour), store 720 values (30days)# 3. AVERAGE, every 720th value (each day), store 5*365 values (5 year)# same for MAX values&callsys($rrdtoolp." create ".$rrdbp.$rrdbn.".rrd --start $tm -s 120 \\DS:".$value_type.":GAUGE:600:".$min.":".$max." \\RRA:AVERAGE:0.5:1:720 \\RRA:AVERAGE:0.5:30:720 \\RRA:AVERAGE:0.5:720:1825 \\RRA:MAX:0.5:1:720 \\RRA:MAX:0.5:30:720 \\RRA:MAX:0.5:720:1825 \\");# print "Database $rrdbn created.\n";return(0);} else {&callsys($rrdtoolp." update ".$rrdbp.$rrdbn.".rrd ".$tm.":".$value);};};sub create_graphs_for {my %hasharray = %{$_[0]};my $wname = $_[1];my %gtimes = ("24h" => 86400,"30d" => 2592000,"1y" => 31536000,"5y" => 157680000);my @color = ( "#dc0000", "#dc00dc", "#5800dc", "#009bdc", "#00dc87" );my $coloridx = 0;# N NNE NE ENE E ESE SE SSE S SSW SW WSW W WNW NW NNwmy @direction_cols = ( "#ff0000", "#ff6600", "#ffcc00", "#ccff00", "#66ff00", "#00ff00", "#00ff66", "#00ffcc", "#00ccff", "#0066ff", "#0000ff", "#6600ff", "#cc00ff", "#ff00cc", "#ff0066", "#ff0000");foreach $gtime (keys %gtimes) {$arguments = "";$coloridx = 0;$img = $wwwdir.$wname."_".$gtime.".gif";$cmd_arguments = $rrdtoolp." graph ".$img." --start -".$gtimes{$gtime}." --width 550 --height 150 --vertical-label \"".ucfirst($wname)."\" --watermark \"WeatherStats\" \\"; #jak udelat scale?? --vertical-label".$scale\" \\";$cmd_arguments = $rrdtoolp." graph ".$img." --start -".$gtimes{$gtime}." --width 550 --height 150 --vertical-label \"".ucfirst($wname)."\" \\"; #jak udelat scale?? --vertical-label".$scale\" \\";# compile the graphs argumentsif ($wname eq "temperature") {foreach $atribute (keys %hasharray) {if ($atribute =~ /out/) {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;}}$arguments = $arguments."HRULE:0#474747";$cmd_arguments_tmp = $cmd_arguments;$cmd_arguments_tmp =~ s/$gtime/out_$gtime/g;&callsys($cmd_arguments_tmp.$arguments);$arguments = "";$coloridx = 0;foreach $atribute (keys %hasharray) {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;}$arguments = $arguments."HRULE:0#474747";} elsif ($wname eq "wind") {foreach $atribute (keys %hasharray) {if ($atribute eq "direction") {$scale = ${$hasharray{$atribute}}[2];if ($gtime = "24h") {$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";} else {$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":MAX ";}$arguments = $arguments."CDEF:North=".$atribute.",0,GE,-4,* ";$arguments = $arguments."CDEF:NNE=".$atribute.",1,GE,-4,* ";$arguments = $arguments."CDEF:NE=".$atribute.",2,GE,-4,* ";$arguments = $arguments."CDEF:ENE=".$atribute.",3,GE,-4,* ";$arguments = $arguments."CDEF:East=".$atribute.",4,GE,-4,* ";$arguments = $arguments."CDEF:ESE=".$atribute.",5,GE,-4,* ";$arguments = $arguments."CDEF:SE=".$atribute.",6,GE,-4,* ";$arguments = $arguments."CDEF:SSE=".$atribute.",7,GE,-4,* ";$arguments = $arguments."CDEF:South=".$atribute.",8,GE,-4,* ";$arguments = $arguments."CDEF:SSW=".$atribute.",9,GE,-4,* ";$arguments = $arguments."CDEF:SW=".$atribute.",10,GE,-4,* ";$arguments = $arguments."CDEF:WSW=".$atribute.",11,GE,-4,* ";$arguments = $arguments."CDEF:West=".$atribute.",12,GE,-4,* ";$arguments = $arguments."CDEF:WNW=".$atribute.",13,GE,-4,* ";$arguments = $arguments."CDEF:NW=".$atribute.",14,GE,-4,* ";$arguments = $arguments."CDEF:NNW=".$atribute.",15,GE,-4,* ";$arguments = $arguments."AREA:North".$direction_cols[0].":N ";$arguments = $arguments."AREA:NNE".$direction_cols[1].": ";$arguments = $arguments."AREA:NE".$direction_cols[2].":NE ";$arguments = $arguments."AREA:ENE".$direction_cols[3].": ";$arguments = $arguments."AREA:East".$direction_cols[4].":E ";$arguments = $arguments."AREA:ESE".$direction_cols[5].": ";$arguments = $arguments."AREA:SE".$direction_cols[6].":SE ";$arguments = $arguments."AREA:SSE".$direction_cols[7].": ";$arguments = $arguments."AREA:South".$direction_cols[8].":S ";$arguments = $arguments."AREA:SSW".$direction_cols[9].": ";$arguments = $arguments."AREA:SW".$direction_cols[10].":SW ";$arguments = $arguments."AREA:WSW".$direction_cols[11].": ";$arguments = $arguments."AREA:West".$direction_cols[12].":W ";$arguments = $arguments."AREA:WNW".$direction_cols[13].": ";$arguments = $arguments."AREA:NW".$direction_cols[14].":\"NW\\l\" ";$arguments = $arguments."AREA:NNW".$direction_cols[15].": ";} else {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":MAX ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;}}} elsif ($wname eq "rain") {foreach $atribute (keys %hasharray) {if ($atribute eq "1h") {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;$cmd_arguments_tmp = $cmd_arguments;$cmd_arguments_tmp =~ s/$gtime/1h_$gtime/g;&callsys($cmd_arguments_tmp.$arguments);$arguments = "";} else {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;}}} else {foreach $atribute (keys %hasharray) {$scale = ${$hasharray{$atribute}}[2];$arguments = $arguments."DEF:".$atribute."=".$rrdbp.$wname."_".$atribute.".rrd:".$atribute.":AVERAGE ";$arguments = $arguments."LINE1:".$atribute.$color[$coloridx].":\"".${$hasharray{$atribute}}[0]."\" ";$arguments = $arguments."GPRINT:".$atribute.":MIN:\"Min\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":AVERAGE:\"Avg\\: %4.2lf \" ";$arguments = $arguments."GPRINT:".$atribute.":MAX:\"Max\\: %4.2lf ".$scale."\\j\" ";$coloridx++;}}&callsys($cmd_arguments.$arguments);}}sub create_graphs {&create_graphs_for(\%temperature, "temperature");&create_graphs_for(\%humidity, "humidity");&create_graphs_for(\%wind, "wind");&create_graphs_for(\%rain, "rain");&create_graphs_for(\%pressure, "pressure");}sub print_header {open (INDEX, ">$wwwdir"."index\.html.new") || die "ERROR: WeatherStats: Could not create $wwwdir index.html.new";print INDEX "<html>\n<head><title>Weather statistic from WATSON W-8681</title><link rel=\"StyleSheet\" href=\"styles.css\" type=\"text/css\" media=\"screen\"></head>\n<body>\n";print INDEX "<a href=p3211521.jpg><img src=p3211521_detail.jpg alt=\"WATSON W-8681\" title=\"WATSON W-8681\" align=\"right\" /></a>\n";print INDEX "<h1>Weather statistic</h1>\n";$date = `date`;print INDEX "<p>Generated: ".$date."</p>\n";print INDEX "<p>Station altitude: <span class=\"wvalue\">".$altitude."</span><span class=\"wscale\">mnm</span></p>\n";}sub print_footer {print INDEX "<div id=\"credits\">WeatherStats - <a href=\"http://meteolinger.lowlevel.cz/\">http://meteolinger.lowlevel.cz/</a></div>\n";print INDEX "</body></html>";close INDEX;open (INDEX, ">$wwwdir"."index\.html") || die "ERROR: WeatherStats: Could not create $wwwdir index.html";open (INDEX_NEW, "$wwwdir"."index\.html.new") || die "ERROR: WeatherStats: Could not create $wwwdir index.html";while (<INDEX_NEW>) { print INDEX $_ };close INDEX_NEW;close INDEX;unlink "index.html.new";}sub print_body_line {my @atribute = @_;printf INDEX "<div>\n<div class=\"wname\"><a href=\"$atribute[0].html\">$atribute[1]</a></div><span class=\"wvalue\"><a href=\"$atribute[0].html\">%.1f</a></span><span class=\"wscale\"><a href=\"$atribute[0].html\">%s</a></span>\n</div>\n", $atribute[3], $atribute[4];}sub print_page {my $page = $_[0];open ( PAGE, ">$wwwdir".$page."\.html") || die "ERROR: WeatherStats: Could not create page $wwwdir $rrdbn";print PAGE "<html>\n<head><title>Weather statistic from TOPCOM 265 NE</title></head>\n<body>\n";print PAGE "<h1>".ucfirst($page)."</h1>\n";print PAGE "<img src=\"".$page."_24h.gif\"/><br />\n";print PAGE "<img src=\"".$page."_30d.gif\"/><br />\n";print PAGE "<img src=\"".$page."_1y.gif\"/><br />\n";print PAGE "<img src=\"".$page."_5y.gif\"/><br />\n";print PAGE "</body></html>";close PAGE}sub print_body_begin {my $wtype = $_[0];print INDEX "<div id=\"".$wtype."\">\n";print INDEX "<span class=\"wimg\"><a href=\"".$wtype.".html\"><img src=\"".$wtype."_24h.gif\" width=\"164\" /></a></span><br />\n";print INDEX "<a class=\"mname\" href=\"".$wtype.".html\">".ucfirst($wtype)."</a><br/>\n";}sub print_body {$wtype = "temperature";&print_body_begin($wtype);$wtype = "temperature_out";&print_body_begin($wtype);print INDEX "</div>\n";&print_page($wtype);$wtype = "temperature";foreach $atribute (keys %temperature) {&print_body_line($wtype,$atribute, ${$temperature{$atribute}}[0], ${$temperature{$atribute}}[1], ${$temperature{$atribute}}[2]);}print INDEX "</div>\n";&print_page($wtype);$wtype = "humidity";&print_body_begin($wtype);foreach $atribute (keys %humidity) {&print_body_line($wtype, $atribute, ${$humidity{$atribute}}[0], ${$humidity{$atribute}}[1], ${$humidity{$atribute}}[2]);}print INDEX "</div>\n";&print_page($wtype);$wtype = "wind";&print_body_begin($wtype);foreach $atribute (keys %wind) {&print_body_line($wtype, $atribute, ${$wind{$atribute}}[0], ${$wind{$atribute}}[1], ${$wind{$atribute}}[2]);}print INDEX "</div>\n";&print_page($wtype);$wtype = "rain";&print_body_begin($wtype);$wtype = "rain_1h";&print_body_begin($wtype);print INDEX "</div>\n";&print_page($wtype);$wtype = "rain";foreach $atribute (keys %rain) {&print_body_line($wtype, $atribute, ${$rain{$atribute}}[0], ${$rain{$atribute}}[1], ${$rain{$atribute}}[2]);}print INDEX "</div>\n";&print_page($wtype);$wtype = "pressure";&print_body_begin($wtype);foreach $atribute (keys %pressure) {&print_body_line("$wtype", $atribute, ${$pressure{$atribute}}[0], ${$pressure{$atribute}}[1], ${$pressure{$atribute}}[2]);}print INDEX "</div>\n";&print_page($wtype);}###################### MAIN #######################@wwsr = `/usr/local/bin/wwsr -y`;#@wwsr = `cat wwsr6.out`;#@wwsr = `tail -14 /var/www/html/wwss/log`;$wwwdir = "/home/www/html/wwss/";#$wwwdir = "wwss/";$en_dbg = 0;$rrdtoolp = "/usr/bin/rrdtool";$rrdbp = "/opt/meteolinger/rrd/";#$rrdbp = "wwss/";$graphs = 0;# What options are defined on command line?if ( defined($ARGV[0]) ) {for ($i = 0; $i < @ARGV; $i++) {$argmnt = $ARGV[$i];if ($argmnt eq "-g") {$graphs = 1;}if ($argmnt eq "-a") {if (defined ($ARGV[$i+1] and length($ARGV[$i+1] == 1)) ) {$altitude = $ARGV[$i+1];$i++;} else {die "ERROR: Unknown or missing argument to -a"}}}} else {&help;};open (LOG, ">>$wwwdir"."log") || die "ERROR: no log";print LOG "------- ".`date`;for ($i=0;$i<=($#wwsr-1);$i++) {print LOG $wwsr[$i];}close LOG;for ($i=0;$i<=($#wwsr-1);$i++) {if ($wwsr[$i] =~ /^For postprocessing/) {$wwsr[$i] = "";last;} else {$wwsr[$i] = "";}if ($i==($#wwsr-1)) {print INDEX "ERROR: wwsr output - for postprocessing sequence not found!</br>";for ($ii=0;$ii<=($#wwsr-1);$ii++) { print INDEX $wwsr[$ii]."</br>"; }}}# define main hashes for value groups%temperature = ();%pressure = ();%humidity = ();%wind = ();%rain = ();# store all values that wwsr returnsforeach $line (@wwsr) {chomp $line;if ( $line ne "" ) {@aline = split ("\t", $line);($value, $scale) = split (" ", $aline[1]);# $value = $aline[$#aline]; #posledni prvek je hodnota$rrdbn = $aline[0];$rrdbn =~ s/ /_/g;$rrdbn = lc($rrdbn);if ($rrdbn eq "interval") { next };&insertdb($rrdbn,$value,$scale,$aline[0]);# &print_page($rrdbn,$aline[0]);# &print_body($rrdbn,$value,$aline[0],$scale);}};# additional calculate values$rrdbnt = "temperature_feelout";$wname = "Feels like temperature";$value = ${$wind{'speed'}}[1];$tout = ${$temperature{'outdoor'}}[1];# windchill calculation is only valid for wind speed > 4.5km/h ~ 1.25m/s and $tout < 10Cif ($value > 1.25 and $tout < 10) {#Convert $value=m/s $ws=km/h$ws=$value*3.6;$ws = $ws**0.16;#Wind chill$twc = 13.12+(0.6215*$tout)-(11.37*($ws))+((0.3965*$tout)*($ws));} else {$twc = $tout;}&insertdb($rrdbnt,$twc,"C",$wname);$rrdbnt = "pressure_airsea";$wname = "Air pressure sea";$value = ${$pressure{'air'}}[1];$scale = ${$pressure{'air'}}[2];# aaMADIS calculation$k1 = 0.190284; # discrepency with calculated k1 probably because Smithsonian used less precise gas constant and gravity values$k2 = 8.4184960528E-5; # (standardLapseRate / standardTempK) * (Power(standardSLP, k1)$aps = ((($value-0.3)**$k1)+($k2*$altitude))**(1/$k1); # Power(Power(pressureHPa - 0.3, k1) + (k2 * elevationM), 1/k1);&insertdb($rrdbnt,$aps,$scale,$wname);$rrdbnt = "temperature_dewout";$wname = "Due point temperature";$value = ${$humidity{'outdoor'}}[1];$tdp = $tout-((100-$value)/5); # aproximation&insertdb($rrdbnt,$tdp,"C",$wname);&print_header;&print_body;&print_footer;if ($graphs == 1) {&create_graphs;}