Я использую {растр}, чтобы обрезать (или обрезать) растр на основе неправильного шейп-файла (биома Amazon), но выход всегда имеет прямоугольный экстент. Однако мне нужен вывод в точно такой же геометрии шейп-файла. Какие-нибудь советы? Приветствия.

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