# Sestavi mapy z gridovych bodu tlaku vzduchu (MSLP) pro sezony s vyskytem sucha
# data: Luterbacher et al. (2002)
# https://www.ncdc.noaa.gov/paleo-search/study/6366

setwd("C:/data/R_data/MFGcviceni")
rm(list = ls()) # vymazani promennych

# potřebné knihovny
if(!require(raster)){install.packages("raster")}
library(raster)
if(!require(maptools)){install.packages("maptools")}
library(maptools)


# načítá SLP pro vybrane suche sezony a normal 61-90
sp.sucho <- read.table("au_sucho.csv", header=TRUE, sep=";")
head(sp.sucho)

# obsah souboru 
# lon, lat - zeměpisná délka a šířka
# au.mean - průměrný tlak v suchých sezónách
# au.normal - průměrný tlak 1961-1990
# au.dif - odchylka od průměru
# au.sign - p-hodnota t-testu

# Prevede data z textoveho formatu x,y,z na raster, prida projekci
ra.sp <- rasterFromXYZ(sp.sucho)
projection(ra.sp)="+init=epsg:4326"

par(mfrow=c(2,2))

plot(ra.sp$au.mean, useRaster=TRUE, interpolate=TRUE, main="Suche sezony")
data(wrld_simpl)
plot(wrld_simpl, add=TRUE)
contour(ra.sp$au.mean, add=TRUE, lwd=2, col="blue")

plot(ra.sp$au.normal, useRaster=TRUE, interpolate=TRUE, main="Normal 1961-1990")
plot(wrld_simpl, add=TRUE)
contour(ra.sp$au.normal, add=TRUE, lwd=2, col="blue")

plot(ra.sp$au.dif, useRaster=TRUE, interpolate=TRUE, main="Vyznamne odchylky od normalu")
plot(wrld_simpl, add=TRUE)
contour(ra.sp$au.dif, add=TRUE, lwd=2, col="blue")

# vybere jen vyznamne diference
sig.au.sucho <-sp.sucho[which(sp.sucho$au.sig<0.051), ]
sig.points <- cbind(sig.au.sucho$lon, sig.au.sucho$lat)
points(sig.points, pch=20)