Я использую {raster} для клипа (или обрезки) растра на основе нерегулярного шейп файла (биома 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")