library(raster)
## Loading required package: sp
setwd("~/Downloads/") et0 <- raster("eto.201906.tif") opad <- raster("rgrs.201906.tif") et0_2 <- projectRaster(from=et0, to=opad) plot(opad-et0_2)