Ich benutze {raster}, um ein auf einem unregelmäßigen Shapefile (dem Amazon-Biom) basierendes Raster auszuschneiden (oder zuzuschneiden), aber die Ausgabe hat immer eine rechteckige Ausdehnung. Ich brauche jedoch die Ausgabe in der exakt gleichen Geometrie des Shapefiles. Irgendwelche Tipps? Prost.

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