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