##########################################################################################
##readme zum Bearbeiten von Rasterdaten im NetCDF (.nc)-Format des Referenzdatensatz 3.0##
################################Erstellt am 26.08.2024####################################
#########################Erstellt von Philipp Koerner, iamk GmbH###########################
##########################################################################################

##Herunterladen der Daten##
1. Moeglichkeit: durch Linksklick auf die Datei, oder Rechtsklick, dann 'verlinkten Inhalt speichern unter'
2. Moeglichkeit, bei mehreren Dateien: Browserplugin "DownloadThemAll"
3. Moeglichkeit: in der shell ueber 'wget', oder in R ueber die Funktion 'download.file'. Unter Windows zusätzlich die Option <mode="wb"> verwenden.

##Einlesen in QGIS##
#QGIS oeffnen#
#Klick auf Layer --> Layer hinzufuegen --> Rasterlayer hinzufuegen#
#Klick bei "Quelle" auf die drei Punkte, Datensatz auswaehlen --> Datensatz erscheint im Fenster "Layer" (in den Standardeinstellungen unten links)#
#Rechsklick auf den Datensatz im Fenster "Layer" --> Einstellungen#
#Bei "Symbolisierung" klick auf die gewuenschte Darstellungsart#
#Bei "Kanal" (je nach Darstellungsart auch grauer, roter, gruener, blauer Kanal) laesst sich der gewuenschte Tag des Jahres auswaehlen#


##Einlesen in R##
#R oeffnen#
#Packet terra laden mit 'library(terra)', wenn noch nicht installiert, vorher mit 'install.packages("terra")' installieren (immer ohne'-Zeichen)#
#Pfad wechseln mit 'setwd("...Pfad wo die .nc-Daten lokal abgelegt sind ...")'
#Daten einlesen mit 'meine_Daten<-rast("1961.nc")'

##Bearbeiten in R##
#Bearbeitung fuer das komplette Jahr gleichzeitig: 'sum(meine_Daten)' ergibt beispielsweise die Jahressumme der Groesse, 'mean' ergaebe den Mittelwert#
#Bearbeitung einzelner Tage: 'plot(meine_Daten[[111]])' plottet die raeumliche Verteilung am 111. Tag des Jahres. 'time(meine_Daten[[111]])' gibt das Datum des Datensatzes an#
#Erstellen von Zeitreihen: Mittlere Werte pro Tag: 'global(meine_Daten,mean,na.rm=T)'
#Zuschneiden auf Einzugsgebiet: 'crop(meine_Daten,mein_EZG,mask=T)'


##ESRI ArcGis##
#hier wird auf die Hilfeseiten von ESRI verwiesen#