Software
Zur Navigation springen
Zur Suche springen
Das Python-Script
#!/usr/bin/env python
# -*- coding: utf-8 -*-
#--------------------------------
import os, io, sys, time
import csv
import picamera
import Adafruit_DHT
import Adafruit_BMP.BMP085 as BMP085
import decimal
from decimal import *
import MySQLdb
import Image
import ImageFont
import ImageDraw
#-----------------------------------------------------------
# Offsets festlegen
#-----------------------------------------------------------
# Sonnenaufgang und -untergang (YY.csv)
OffsetAuf = 25
OffsetUnter = 30
#-----------------------------------------------------------
# Logfile definieren
#-----------------------------------------------------------
logfile = '/home/pi/timelapse.log'
#-----------------------------------------------------------
# Funktion Logfile schreiben
# log("string"+"\n")
#-----------------------------------------------------------
def log(msg):
file = open(logfile,"a")
file.write(msg)
file.close
#-----------------------------------------------------------
# Funktion Temperatur in Farbe umrechnen
#-----------------------------------------------------------
def tempcolor(temp):
gruppe = int(temp/10) + 2
mod = temp - int(temp/10)*10
if gruppe < 0:
gruppe = 0
if gruppe > 6:
gruppe = 6
if gruppe == 6:
r = 255
g = 0
b = 0
return (r, g, b)
elif gruppe == 5:
r = 200 + int(Decimal(5.5) * mod)
g = 0
b = 0
return (r, g, b)
elif gruppe == 4:
r = 200
g = 200 - int(20 * mod)
b = 0
return (r, g, b)
elif gruppe == 3:
r = 0 + int(20 * mod)
g = 200
b = 0
return (r, g, b)
elif gruppe == 2:
if temp >= 0:
r = 0
g = 200
b = 200 - int(20 * mod)
return (r, g, b)
else:
r = 0
g = 200 + int(20 * mod)
b = 200
return (r, g, b)
elif gruppe == 1:
r = 0 - int(20 * mod)
g = 0
b = 200
return (r, g, b)
else:
r = 200
g = 0
b = 100
return (r, g, b)
#-----------------------------------------------------------
# Funktion Humidity in Farbe umrechnen
#-----------------------------------------------------------
def humcolor(hum):
r = 255 - int(hum*255/100)
g = 255 - int(hum*255/100)
b = 255
return (r, g, b)
#-----------------------------------------------------------
# Function DS18B20 lesen
#-----------------------------------------------------------
def Temp(sensorfile):
# 1-wire Slave Datei lesen
file = open('/sys/bus/w1/devices/' + sensorfile + '/w1_slave')
filecontent = file.read()
file.close()
# Temperaturwerte auslesen und konvertieren
stringvalue = filecontent.split("\n")[1].split(" ")[9]
temperature = float(stringvalue[2:]) / 1000
# Temperatur zurückgeben
return(temperature)
#-----------------------------------------------------------
# Sensoren auswerten
#-----------------------------------------------------------
# DS18B20 auswerten ------------
temp1 = Temp('28-0516851390ff')
temp = Decimal(temp1 -0.05).quantize(Decimal('.1'), rounding=ROUND_HALF_UP)
# DHT-22 auswerten ------------------
sensor = Adafruit_DHT.DHT22
# GPIO-Nummer festlegen -------------
pin = 17
humidity, temperature = Adafruit_DHT.read_retry(sensor, pin)
# humkorr = humidity * 1.13
humkorr = humidity * 1.15
# humkorr = humidity * 1.22
if humkorr > 100:
humkorr = 100
hum = Decimal(humkorr).quantize(Decimal('.1'), rounding=ROUND_HALF_UP)
humraw = Decimal(humidity).quantize(Decimal('.1'), rounding=ROUND_HALF_UP)
# Luftdruck bestimmen ---------------
bmp = BMP085.BMP085()
BMP_temp = bmp.read_temperature() - 0.5
tempint = Decimal(BMP_temp).quantize(Decimal('.1'), rounding=ROUND_HALF_UP)
BMP_hpa = bmp.read_pressure()
# Meereshoehe Standort in m = 315 + Kalibrierwert
Kalib = -9.0
altitude = 315.0 + Kalib
BMP_psea = BMP_hpa / pow(1.0 - altitude/44330.0, 5.255) / 100
hpa = Decimal(BMP_psea).quantize(Decimal('.1'), rounding=ROUND_HALF_UP)
#-----------------------------------------------------------
# Datenbank schreiben
#-----------------------------------------------------------
# Aktuelles Datum und Zeit bestimmen
# Zeit = time.strftime("%y%m%d-%H%M%S")
Zeit = time.strftime("%y%m%d-%H%M")
UTC = int(time.time()/60)*60
YY = Zeit[0:2]
MM = Zeit[2:4]
DD = Zeit[4:6]
Stunde = Zeit[7:9]
Minute = Zeit[9:11]
Datum = Datum = DD + "." + MM + "." + YY
TTMM = DD + MM
JJMMTT = YY + MM + DD
TagMin = int(Stunde) * 60 + int(Minute)
ZeitPIC = DD + "." + MM + "." + "20" + YY + " " + Stunde + ":" + Minute
timestamp = YY + MM + DD + Stunde + Minute
# Vor dem Schreiben Plausibilitätsprüfung von Zeit und Temperatur durchführen
# Datenbank definieren
db = MySQLdb.connect("localhost","pi","pedlemic01","sensor" )
cursor = db.cursor()
# Letzten Datensatz aus werte lesen
sql = "SELECT * FROM werte ORDER BY id DESC LIMIT 1"
cursor.execute(sql)
checktime = cursor.fetchone()
if checktime is not None:
check_utc = checktime[1]
check_temp = checktime[3]
if (UTC - float(check_utc)) > 301:
log("E8: " + timestamp + " " + str((UTC - float(check_utc))/60) + " Minuten Zeitdifferenz!" + "\n")
if temp > 42:
log("E9: " + timestamp + " Temperaturwert > 42: " + str(temp))
temp = check_temp
# 5-Minutenwerte schreiben Tabelle werte --------------------
sql = "INSERT INTO werte (utc, timestamp, temp, hum, humraw, hpa, tempint) \
VALUES ('%s', '%s', '%f', '%f', '%f', '%f', '%f')" % \
(UTC, timestamp, temp, hum, humraw, hpa, tempint)
try:
cursor.execute(sql)
db.commit()
except:
log("E1: " + timestamp + " SQL write to werte not ok!" +"\n")
db.rollback()
# 3-Stundenwert hpa aus Tabelle werte lesen ----------------
UTC_3h = UTC - 10800
sql = "SELECT hpa FROM werte WHERE UTC = %s"
cursor.execute(sql, (UTC_3h,))
UTC_3hDB = cursor.fetchone()
# Ausmaß der Luftdruckänderung feststellen -----------------
hpaSign = ""
hpadiff = 0
if UTC_3hDB is not None:
hpadiff = hpa - UTC_3hDB[0]
if (hpadiff > -1) and (hpadiff < 1):
hpaSign = "="
elif (hpadiff >= 1) and (hpadiff < 3):
hpaSign = "+"
elif (hpadiff >= 3) and (hpadiff < 5):
hpaSign = "++"
elif (hpadiff >= 5):
hpaSign = "+++"
elif (hpadiff <= -1) and (hpadiff > -3):
hpaSign = "-"
elif (hpadiff <= -3) and (hpadiff > -5):
hpaSign = "--"
elif (hpadiff <= -5):
hpaSign = "---"
# Extremwerte behandeln Tabelle hochtief -------------------
# Tabelle hochtief Tageshochtief lesen ---------------------
timestamp_ht = YY + MM + DD
sql = "SELECT temph, tempt, humh, humt, hpah, hpat, id FROM hochtief WHERE timestamp = %s"
cursor.execute(sql, (timestamp_ht,))
ergebnis = cursor.fetchone()
if not ergebnis:
sql = "INSERT INTO hochtief (timestamp, temph, tempt, humh, humt, hpah, hpat) \
VALUES ('%s', '%f', '%f', '%f', '%f', '%f', '%f')" % \
(timestamp_ht, temp, temp, hum, hum, hpa, hpa)
try:
cursor.execute(sql)
db.commit()
except:
log("E2: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
tempstrTh = "%.1f" % temp
tempstrTt = "%.1f" % temp
tempfloath = float(temp)
tempfloatt = float(temp)
humfloath = float(hum)
humfloatt = float(hum)
hpafloath = float(hpa)
hpafloatt = float(hpa)
hpastrTh = "%.1f" % hpa
hpastrTt = "%.1f" % hpa
else:
speichern = False
temph = ergebnis[0]
tempt = ergebnis[1]
tempstrTh = "%.1f" % temph
tempstrTt = "%.1f" % tempt
tempfloath = float(temph)
tempfloatt = float(tempt)
humh = ergebnis[2]
humt = ergebnis[3]
humfloath = float(humh)
humfloatt = float(humt)
hpah = ergebnis[4]
hpat = ergebnis[5]
hpafloath = float(hpah)
hpafloatt = float(hpat)
hpastrTh = "%.1f" % hpah
hpastrTt = "%.1f" % hpat
id = ergebnis[6]
if temp > temph:
temph = temp; speichern = True
tempstrTh = "%.1f" % temph
tempfloath = float(temph)
if temp < tempt:
tempt = temp; speichern = True
tempstrTt = "%.1f" % tempt
tempfloatt = float(tempt)
if hum > humh:
humh = hum; speichern = True
humfloath = float(humh)
if hum < humt:
humt = hum; speichern = True
humfloatt = float(humt)
if hpa > hpah:
hpah = hpa; speichern = True
hpastrTh = "%.1f" % hpah
hpafloath = float(hpah)
if hpa < hpat:
hpat = hpa; speichern = True
hpastrTt = "%.1f" % hpat
hpafloatt = float(hpat)
if speichern:
try:
cursor.execute("""
UPDATE hochtief
SET temph=%s, tempt=%s, humh=%s, humt=%s, hpah=%s, hpat=%s
WHERE id=%s
""", (temph, tempt, humh, humt, hpah, hpat, id))
db.commit()
except:
log("E3: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
# Tabelle hochtief Monatshochtief lesen --------------------
timestamp_ht = YY + MM + "00"
sql = "SELECT temph, tempt, humh, humt, hpah, hpat, id FROM hochtief WHERE timestamp = %s"
cursor.execute(sql, (timestamp_ht,))
ergebnis = cursor.fetchone()
if not ergebnis:
sql = "INSERT INTO hochtief (timestamp, temph, tempt, humh, humt, hpah, hpat) \
VALUES ('%s', '%f', '%f', '%f', '%f', '%f', '%f')" % \
(timestamp_ht, temp, temp, hum, hum, hpa, hpa)
try:
cursor.execute(sql)
db.commit()
except:
log("E4: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
tempstrMh = "%.1f" % temp
tempstrMt = "%.1f" % temp
else:
speichern = False
temph = ergebnis[0]
tempt = ergebnis[1]
tempstrMh = "%.1f" % temph
tempstrMt = "%.1f" % tempt
humh = ergebnis[2]
humt = ergebnis[3]
hpah = ergebnis[4]
hpat = ergebnis[5]
id = ergebnis[6]
if temp > temph:
temph = temp; speichern = True
tempstrMh = "%.1f" % temph
if temp < tempt:
tempt = temp; speichern = True
tempstrMt = "%.1f" % tempt
if hum > humh:
humh = hum; speichern = True
if hum < humt:
humt = hum; speichern = True
if hpa > hpah:
hpah = hpa; speichern = True
if hpa < hpat:
hpat = hpa; speichern = True
if speichern:
try:
cursor.execute("""
UPDATE hochtief
SET temph=%s, tempt=%s, humh=%s, humt=%s, hpah=%s, hpat=%s
WHERE id=%s
""", (temph, tempt, humh, humt, hpah, hpat, id))
db.commit()
except:
log("E5: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
# Tabelle hochtief Jahreshochtief lesen --------------------
timestamp_ht = YY + "0000"
sql = "SELECT temph, tempt, humh, humt, hpah, hpat, id FROM hochtief WHERE timestamp = %s"
cursor.execute(sql, (timestamp_ht,))
ergebnis = cursor.fetchone()
if not ergebnis:
sql = "INSERT INTO hochtief (timestamp, temph, tempt, humh, humt, hpah, hpat) \
VALUES ('%s', '%f', '%f', '%f', '%f', '%f', '%f')" % \
(timestamp_ht, temp, temp, hum, hum, hpa, hpa)
try:
cursor.execute(sql)
db.commit()
except:
log("E6: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
tempstrJh = "%.1f" % temp
tempstrJt = "%.1f" % temp
else:
speichern = False
temph = ergebnis[0]
tempt = ergebnis[1]
tempstrJh = "%.1f" % temph
tempstrJt = "%.1f" % tempt
humh = ergebnis[2]
humt = ergebnis[3]
hpah = ergebnis[4]
hpat = ergebnis[5]
id = ergebnis[6]
if temp > temph:
temph = temp; speichern = True
tempstrJh = "%.1f" % temph
if temp < tempt:
tempt = temp; speichern = True
tempstrJt = "%.1f" % tempt
if hum > humh:
humh = hum; speichern = True
if hum < humt:
humt = hum; speichern = True
if hpa > hpah:
hpah = hpa; speichern = True
if hpa < hpat:
hpat = hpa; speichern = True
if speichern:
try:
cursor.execute("""
UPDATE hochtief
SET temph=%s, tempt=%s, humh=%s, humt=%s, hpah=%s, hpat=%s
WHERE id=%s
""", (temph, tempt, humh, humt, hpah, hpat, id))
db.commit()
except:
log("E7: " + timestamp_ht + " SQL write to hochtief not ok!" +"\n")
db.rollback()
# disconnect from server
db.close()
# -------------------------------------------------------------------------
# PiCam
#--------------------------------------------------------------------------
# Directories festlegen, checken (und anlegen)
Jahr = time.strftime("%Y")
dirJahr = "/media/usbhdd/" + Jahr + "/"
try: os.stat(dirJahr)
except: os.mkdir(dirJahr)
# Ordner für Nachtbilder nicht mehr anlegen
# dirJahrN = "/media/usbhdd/" + Jahr + "N/"
# try: os.stat(dirJahrN)
# except: os.mkdir(dirJahrN)
# Foto machen --------------------------------
# Create the in-memory stream
stream = io.BytesIO()
with picamera.PiCamera() as camera:
camera.resolution = (1920, 1080)
camera.start_preview()
camera.meter_mode = 'matrix'
# camera.sharpness = 10
# camera.saturation = 5
# Camera warm-up time
time.sleep(2)
camera.capture(stream, format='jpeg')
# "Rewind" the stream to the beginning so we can read its content
camera.stop_preview()
stream.seek(0)
bild = Image.open(stream)
camera.close()
# Bild beschriften-------------------------------------------------------
#------------------------------------------------------------------------
draw = ImageDraw.Draw(bild, 'RGBA')
font = ImageFont.truetype("/usr/share/fonts/truetype/dejavu/DejaVuSans.ttf",35)
fontk = ImageFont.truetype("/usr/share/fonts/truetype/dejavu/DejaVuSans.ttf",30)
fontg = ImageFont.truetype("/usr/share/fonts/truetype/dejavu/DejaVuSans.ttf",40)
draw.rectangle((0,991,1919,1079),fill=(80,80,80,144))
draw.rectangle((0,0,80,990),fill=(80,80,80,144))
draw.rectangle((1840,0,1919,990),fill=(80,80,80,144))
# Zeit schreiben --------------------------------------------------------
draw.text((800, 1035),ZeitPIC,(255,255,255),font=font)
# Temperatur schreiben --------------------------------------------------
# -----------------------------------------------------------------------
draw.text((3,994),chr(176)+"C",(255,255,255),font=font)
tempstr = "%.1f" % temp
draw.text((3,1035),tempstr,(255,255,255),font=font)
# Temperaturbalken schreiben --------------------
# oben=0, unten=990, Null Grad=660, 1 Grad=16,5
rot, gruen, blau = tempcolor(temp)
farbe = "rgb("+str(rot)+","+str(gruen)+","+str(blau)+")"
tempfloat = float(temp)
y = int(660 + 0.5 - tempfloat * 16.5)
if y < 0:
y = 0
if y > 990:
y = 990
coor = (0,y,30,990)
draw.rectangle(coor, fill=farbe)
coor = (0,0,30,990)
draw.rectangle(coor, outline="rgb(255,255,255)")
# Hoch- und Tiefmarkierung
y = int(660 + 0.5 - tempfloath * 16.5)
if y < 1: y = 1
if y > 989: y = 989
coor = (0,y-1,30,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
y = int(660 + 0.5 - tempfloatt * 16.5)
if y < 1: y = 1
if y > 989: y = 989
coor = (0,y-1,30,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
# Null-Grad Markierung
coor = (0,659,39,661)
draw.rectangle(coor, fill="rgb(255,255,255)")
draw.text((40,639),"0"+chr(176),(255,255,255),font=fontg)
# max/min schreiben -------------------------------
draw.text((110,1003),"max",(255,255,255),font=fontk)
draw.text((110,1039),"min",(255,255,255),font=fontk)
draw.text((195,1020),"D",(255,255,255),font=font)
draw.text((235,1003),tempstrTh,(255,255,255),font=fontk)
draw.text((235,1039),tempstrTt,(255,255,255),font=fontk)
draw.text((335,1020),"M",(255,255,255),font=font)
draw.text((375,1003),tempstrMh,(255,255,255),font=fontk)
draw.text((375,1039),tempstrMt,(255,255,255),font=fontk)
draw.text((475,1020),"Y",(255,255,255),font=font)
draw.text((515,1003),tempstrJh,(255,255,255),font=fontk)
draw.text((515,1039),tempstrJt,(255,255,255),font=fontk)
# Luftfeuchtigkeit ------------------------------------------------------
# -----------------------------------------------------------------------
humstr = str(int(float(hum) + 0.5))
draw.text ((1850,994),"rF%",(255,255,255),font=font)
draw.text ((1850,1035),humstr,(255,255,255), font=font)
# Luftfeuchtigkeitbalken schreiben ---------------
# oben=0, unten=990, 1 Prozent=9,9
rot, gruen, blau = humcolor(hum)
farbe = "rgb("+str(rot)+","+str(gruen)+","+str(blau)+")"
humfloat = float(hum)
y = int(990 + 0.5 - humfloat * 9.9)
coor = (1889,y,1919,990)
draw.rectangle(coor, fill=farbe)
coor = (1889,0,1919,990)
draw.rectangle(coor, outline="rgb(255,255,255)")
# Hoch- und Tiefmarkierung
y = int(990 + 0.5 - humfloath * 9.9)
if y < 1: y = 1
if y > 989: y = 989
coor = (1889,y-1,1919,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
y = int(990 + 0.5 - humfloatt * 9.9)
if y < 1: y = 1
if y > 989: y = 989
coor = (1889,y-1,1919,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
# Luftdruck -------------------------------------------------------------
#------------------------------------------------------------------------
draw.text((1700,994),"hPa",(255,255,255),font=font)
hpastr = "%.1f" % hpa
draw.text((1700, 1035),hpastr,(255,255,255),font=font)
# Luftdruckbalken schreiben ----------------------
# oben=0, unten=990, 1 hPa=16,5
# Farbe festlegen
farbe = "rgb(170,170,170)"
hpafloat = float(hpa)
hpafloat = hpafloat - 980.0
y = int(990 + 0.5 - hpafloat * 16.5)
if y < 0: y = 0
if y > 990: y = 990
coor = (1839,y,1869,990)
draw.rectangle(coor, fill=farbe)
coor = (1839,0,1869,990)
draw.rectangle(coor, outline="rgb(255,255,255)")
# Hoch- und Tiefmarkierung
y = int(990 + 0.5 - (hpafloath - 980.0) * 16.5)
if y < 1: y = 1
if y > 989: y = 989
coor = (1839,y-1,1869,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
y = int(990 + 0.5 - (hpafloatt - 980.0) * 16.5)
if y < 1: y = 1
if y > 989: y = 989
coor = (1839,y-1,1869,y+1)
draw.rectangle(coor, fill="rgb(255,255,255)")
# Luftdruck steigend oder fallend markieren
if hpaSign == "=":
coor = (1772,1013,1816,1017)
draw.rectangle(coor, fill="rgb(255,225,0)")
elif hpaSign == "+":
draw.polygon(((1772,1025),(1792,1025),(1782,1005)),fill="rgb(0,255,0)")
elif hpaSign == "++":
draw.polygon(((1772,1025),(1792,1025),(1782,1005)),fill="rgb(0,255,0)")
draw.polygon(((1792,1025),(1812,1025),(1802,1005)),fill="rgb(0,255,0)")
elif hpaSign == "+++":
draw.polygon(((1772,1025),(1792,1025),(1782,1005)),fill="rgb(0,255,0)")
draw.polygon(((1792,1025),(1812,1025),(1802,1005)),fill="rgb(0,255,0)")
draw.polygon(((1812,1025),(1832,1025),(1822,1005)),fill="rgb(0,255,0)")
elif hpaSign == "-":
draw.polygon(((1772,1005),(1792,1005),(1782,1025)),fill="rgb(255,0,0)")
elif hpaSign == "--":
draw.polygon(((1772,1005),(1792,1005),(1782,1025)),fill="rgb(255,0,0)")
draw.polygon(((1792,1005),(1812,1005),(1802,1025)),fill="rgb(255,0,0)")
elif hpaSign == "---":
draw.polygon(((1772,1005),(1792,1005),(1782,1025)),fill="rgb(255,0,0)")
draw.polygon(((1792,1005),(1812,1005),(1802,1025)),fill="rgb(255,0,0)")
draw.polygon(((1812,1005),(1832,1005),(1822,1025)),fill="rgb(255,0,0)")
# max/min schreiben -------------------------------
draw.text((1530,1020),"D",(255,255,255),font=font)
draw.text((1570,1003),hpastrTh,(255,255,255),font=fontk)
draw.text((1570,1039),hpastrTt,(255,255,255),font=fontk)
# Eventtext schreiben -----------------------------
EventArray = [];
EventText = ""
# Datei Event.csv in Array EventArray einlesen
csvReader = csv.reader(open('/home/pi/Event.csv', 'r'), delimiter=';');
for row in csvReader:
EventArray.append(row);
# JJMMTT in Array EventArray suchen - mit 1 beginnen wegen Überschrift
for i in range(1, len(EventArray)):
if EventArray[i][0] == JJMMTT:
EventText = EventArray[i][1]
break
if EventText <> "":
draw.rectangle((81,0,1838,55),fill=(80,80,80,144))
draw.text((38,3),EventText,(255,255,255),font=fontg)
# -----------------------------------------------------------------------
# Array mit Sonnenauf- und untergangszeiten (Minuten)
SonneAU = [];
# Datei YY.csv in Array SonneAU einlesen
csvReader = csv.reader(open('/home/pi/' + YY + '.csv', 'r'), delimiter=';');
for row in csvReader:
SonneAU.append(row);
# TTMM und Tagesminute TagMin in Array SonneAU suchen
for i in range(0, len(SonneAU)):
if SonneAU[i][0] == TTMM:
SonneAuf = int(SonneAU[i][3]) - OffsetAuf
SonneUnter = int (SonneAU[i][4]) + OffsetUnter
break
# Bild speichern
if (SonneAuf <= int(TagMin)) and (SonneUnter >= int(TagMin)):
# bild.save(dirJahr+'P'+Zeit+'.jpg', format='JPEG', subsampling=0, quality=100)
bild.save(dirJahr+'P'+Zeit+'.jpg')
# Nachtbild nicht mehr speichern
# else:
# bild.save(dirJahrN+'P'+Zeit+'.jpg')
bild.save('/media/usbhdd/webpic/weather.jpg')
# Programmende----------------------------------------------------------------
Die Datei Event.csv hat folgenden Aufbau:
JJMMTT;Eventtext (mit Blanks und ohne Umlaute) 181201;Dezember 181202;1. Adventsonntag 181209;2. Adventsonntag 181216;3. Adventsonntag 181223;4. Adventsonntag ...
Die Datei JJ.csv (z.B. 17.csv) enthält die Zeiten der Sonnenauf- und untergänge und sieht wie folgt aus:
TTMM;Auf;Unter;Auf;Unter 0101;07:46;16:11;466;971 0201;07:45;16:12;465;972 0301;07:45;16:13;465;973 0401;07:45;16:14;465;974 0501;07:45;16:15;465;975 0601;07:45;16:16;465;976 ...
Die Daten können z.B. hier https://www.zamg.ac.at/cms/de/klima/klimauebersichten/ephemeriden bezogen werden und müssen (mit Excel) entsprechen aufbereitet werden. Die letzten beiden Parameter sind die Anzahl der Minuten ab 00:00 Uhr.