library(lubridate)
library(ggplot2)
library(StreamMetabolism)
library(xts)
library(reshape)
library(scales)
19300_Grabow <- sunrise.set(53.2760517,11.5636518, "2024/01/01", timezone="MET", num.days=370)
sunrise <- 19300_Grabow$sunrise
sunset <- 19300_Grabow$sunset
sunrise <- strftime(sunrise, format="%R", tz="MET")
sunset <- strftime(sunset, format="%R", tz="MET")
19300_Grabow["sr"] <- as.POSIXct(sunrise, format = "%H:%M")
19300_Grabow["ss"] <- as.POSIXct(sunset, format = "%H:%M")
19300_Grabow["timestamp"] <- align.time(19300_Grabow$sunrise, 60*10)
19300_Grabow <- 19300_Grabow[c("timestamp", "sr", "ss")]
locsrss <- ggplot(19300_Grabow, aes(x=19300_Grabow$timestamp)) + geom_line(aes(y=19300_Grabow$sr)) + geom_line(aes(y=19300_Grabow$ss)) + labs(title = " Sonnenauf-/Sonnenuntergang - 19300_Grabow 2024", x = "Datum", y = "Zeit")
pdf("19300_Grabow_SA_SU.pdf", paper="a4r", width=11)
locsrss
dev.off()
png(filename="19300_Grabow_SA_SU.png", width = 1400, height = 800, units = "px")
locsrss
dev.off()
19300_Grabow["Sonnenaufgang"] <- strftime(19300_Grabow$sr, format="%H:%M")
19300_Grabow["Sonnenuntergang"] <- strftime(19300_Grabow$ss, format="%H:%M")
write.table(19300_Grabow, file="19300_Grabow_SaSu.csv", dec=',', sep=';', row.names=FALSE)