library(raster)
library(leaflet)

ph <- getData("GADM", country = "PHL", level = 1)

leaflet(data = ph) %>%
  addProviderTiles("OpenWeatherMap.Precipitation") %>%
  addPolygons(fillColor = "red",
              fillOpacity = 1,
              color = "#FFFFFF",
              weight = 1.3)

library(mapview)
m <- mapview(ph, map.types = "OpenWeatherMap.Precipitation")
m

library(rgdal)
library(RCurl)
library(httr)
library(jsonlite)
library(leaflet)
library(mapview)
library(owmr)
library(raster)
library(RColorBrewer)
setwd("D:/Works/Fernando/CRS/R/shapes/")
getwd()
file.exists(".","contorno.shp")
ph <- raster::getData('worldclim',var='tmin',res=0.5, lon=-86.921641,lat=14.761664)
pais<-rgdal::readOGR(dsn=getwd(),layer="contorno")
zona<-rgdal::readOGR(dsn=getwd(),layer="La_Paz_Valle")
#newproj <- "+proj=longlat +datum=WGS84"
crs(ph$tmin1_23)
crs(pais)
prj<-"+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
pais2<-spTransform(pais,crs(prj))
zona2<-spTransform(zona,crs(prj))
crs(pais2)
cp<-crop(ph$tmin1_23,pais2)
s<-stack(ph$tmin1_23)
coors<-read.csv("sequia.csv")
coords<-data.frame(lon=coors$LONGITUDE, lat=coors$LATITUDE)
coordinates(coords)<-c("lon","lat")
val<-extract(x=ph$tmin1_23, y=coords)
#proj4string(xy) <- CRS(prj)  ## for example
#sfc<-st_sfc(xy,prj)
#xy<-data.frame(x=coors$LONGITUDE,Y=coors$LATITUDE)
#require(sf)
#res<-SpatialPoints(cbind(coors$LONGITUDE, -coors$LATITUDE), proj4string=CRS("+proj=longlat"))
#as(res,"SpatialPoints")
#rasValue<-lapply(ph$tmin1_23,extract, res)
colores <- c('red', 'green', 'blue', 'chocolate', 'deeppink', 'grey')
at <- seq(0, 4800, 800)
cb <- colorBin(palette = colores, bins = at, domain = at)
m<-leaflet() %>%
  addTiles()%>%
  addRasterImage(cp,opacity = 0.8)%>%
  addLegend(pal = cb, values = at)%>%
  addPolygons(data = zona2,fillColor = "#ff0000",opacity = 0.1)
m<-addMarkers(m,lng=coors$LONGITUDE,lat=coors$LATITUDE,popup=paste("Region:",coors$Municipality,"<br>",
                                                                   "Precipitacion:","<br>",
                                                                   val,"mm/mes"))
#m<-addMarkers(m,lng=-87.8333,lat = 14.1167,popup=paste("Guajiquiro","prec:",val,sep = " "))
#m<-addMarkers(m,lng = -87.7,lat = 13.9667,popup="San Juan")
#m<-addMarkers(m,lng = -87.95,lat = 14.0667,popup = "Santa Ana")
#m<-addMarkers(m,lng = -87.5333,lat = 13.65,popup = "San Francisco de Coray")
#m<-addMarkers(m,lng = -87.6525,lat = 13.62083,popup = "Langue")
#m<-addMarkers(m,lng = -87.7167,lat = 13.7333,popup = "Aramecina")
#m %>% mapview(ph, map.types = "OpenWeatherMap.Precipitation")
m