Estoy usando {ráster} para recortar (o recortar) un ráster basado en un archivo de forma irregular (el bioma de Amazon) pero la salida siempre tiene una extensión rectangular. Sin embargo, necesito la salida en la misma geometría exacta del archivo de forma. ¿Algun consejo? Salud.

library(raster)
library(rgdal)

myshp <- readOGR("Amazon.shp", layer="Amazon")
e <- extent(myshp)
myraster <- raster("Temperature.tif")
myraster.crop <- crop(myraster, e, snap="out", filename="myoutput.tif")