dd

library(sf)

## Linking to GEOS 3.6.1, GDAL 2.2.3, PROJ 4.9.3

library(tidyverse)

## -- Attaching packages ---------------------------------------------- tidyverse 1.2.1 --

## v ggplot2 3.1.0       v purrr   0.3.2  
## v tibble  2.1.1       v dplyr   0.8.0.1
## v tidyr   0.8.3       v stringr 1.4.0  
## v readr   1.3.1       v forcats 0.4.0

## -- Conflicts ------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

library(dplyr)
library(rgdal)

## Loading required package: sp

## rgdal: version: 1.4-3, (SVN revision 828)
##  Geospatial Data Abstraction Library extensions to R successfully loaded
##  Loaded GDAL runtime: GDAL 2.2.3, released 2017/11/20
##  Path to GDAL shared files: C:/Users/TOSHIBA/Documents/R/win-library/3.5/rgdal/gdal
##  GDAL binary built with GEOS: TRUE 
##  Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
##  Path to PROJ.4 shared files: C:/Users/TOSHIBA/Documents/R/win-library/3.5/rgdal/proj
##  Linking to sp version: 1.3-1

library(RColorBrewer)
library(raster)

## 
## Attaching package: 'raster'

## The following object is masked from 'package:dplyr':
## 
##     select

## The following object is masked from 'package:tidyr':
## 
##     extract

a = raster("D:/UC/CURSOS UC - 2019 I/Seminario de investigacion Departamental/hfp-s-america-geo-grid/hfp-s-america-geo-grid/hfp_s_america_grid/hfp_s_amer/w001001.adf")

library(raster)
library(rgdal)
Junin <- read_sf("D:/UC/CURSOS UC - 2019 I/Seminario de investigacion Departamental/Datos para Junin/Junin.shp")

HFPJunin <- crop(a, Junin)
HFPJunin <- HFPJunin <- mask(HFPJunin, Junin)
plot(HFPJunin)

APJu <- read_sf("D:/UC/CURSOS UC - 2019 I/Seminario de investigacion Departamental/Datos para Junin/APJunin.shp") %>% as_Spatial()

plot(HFPJunin)
plot(APJu, add = T)

Dif = read_rds("D:/UC/CURSOS UC - 2019 I/Seminario de investigacion Departamental/Datos para Junin/Dif.rds") 
plot(Dif)
plot(APJu, add = T)