SwampThingPaul

## column names for the different data
col.names.WL <- c("DateTime.UTC","value.m","grade.code","approve.level","qual")
col.names.Q <- c("DateTime.UTC","value.cms","grade.code","approve.level","qual")

## Build the html query
servfull <- "https://ntg.aquaticinformatics.net/Export/DataSet?"
dateMin <- as.Date("2005-01-01 00:00:00")
dateMax <- as.Date("2005-12-01 00:00:00")
location <- "G8140040"

SwampThingPaul

library(zyp)
library(mblm)

# From https://library.virginia.edu/data/articles/theil-sen-regression-programming-and-understanding-an-outlier-resistant-alternative-to-least-squares
theil_sen <- function(x,y){
  n <- length(x)
  max_n_slopes <- (n * (n - 1)) / 2
  slopes <- vector(mode = 'list', length = max_n_slopes) # list of NULL elements
  add_at <- 1
  # Calculate point-to-point slopes if x_i != x_j

SwampThingPaul

library(ggplot2)
library(dplyr)

theme_set(theme_minimal(base_size = 16)+
            theme_bw()+
            theme(panel.border = element_rect("black",fill=NA,linewidth=1)))
# Create sample data
set.seed(123)
n <- 100
groups <- 5

SwampThingPaul

# Libraries
# devtools::install_github("SwampThingPaul/AnalystHelper")
library(AnalystHelper);
library(plyr)
library(reshape2)

# -------------------------------------------------------------------------
# DATA_CSV = data location on PC
caze.ice = read.csv(DATA_CSV,sep="\t",na.strings = "-999")

SwampThingPaul

## Example Recurrence (or return) interval analysis
## functions to calculate return interval for 
## empirical, Gumbel and Log Pearson distributions

library(lmom)

recur.fun=function(x){
  #Sorting data by decreasing order
  sorted.values<-sort(x,decreasing=T)
  

SwampThingPaul

## Example spatial interpolation

# GIS libraries 
library(rgdal)
library(rgeos)
library(raster)

# kriging
library(gstat)

SwampThingPaul

## Objective: To query FDEP's ArcGIS REST server via R
## You must have server link and know a little about the dataset. 
##
## Due to the query limit and area of interest, query is broken into two regions

### Libraries
library(httr)
library(jsonlite)

path <- 'https://services1.arcgis.com/nRHtyn3uE1kyzoYc/arcgis/rest/services/VIEW_FL_Algal_Bloom_Site_Visits_1/FeatureServer/0/query?'

SwampThingPaul

library(sp)
library(rgeos)
library(rworldmap)


# epsg.io
wgs=CRS("+init=epsg:4326")

map_dat=getMap()
map_dat=spTransform(map_dat,wgs)

SwampThingPaul

# This is a very general script to perform gif animations from 
# space-time Generalized Additive Model predictions

# GIS libraries 
library(rgdal)
library(rgeos)
library(raster)
library(tmap)

# GAM libraries

SwampThingPaul

# GIS libraries 
library(sp)
library(rgdal)
library(rgeos)
library(raster)
library(mapmisc)

## defines the coordinate reference system
utm17=CRS("+init=epsg:26917")