Estou usando {raster} para recortar (ou cortar) uma raster com base em um arquivo de forma irregular (o bioma Amazônia), mas a saída sempre tem uma extensão retangular. No entanto, eu preciso da saída exatamente na mesma geometria do shapefile. Alguma dica? Felicidades.

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")