L'écrêtage raster à l'aide de shapefile dans R, mais en gardant la géométrie de la shapefile

Je suis à l'aide de {raster} à clip (ou de cultures) un raster basé sur un fichier de forme irrégulière (le biome Amazonien), mais la sortie a toujours un rectangulaire mesure. Cependant, j'ai besoin de la sortie dans la même géométrie de la shapefile. Des conseils à donner? Des acclamations.

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

OriginalL'auteur ALS.Meyer | 2014-04-15

  1. 24

    Une option est d'utiliser raster::mask()

    library(maptools)  ## For wrld_simpl
library(raster)

## Example SpatialPolygonsDataFrame
data(wrld_simpl)
SPDF <- subset(wrld_simpl, NAME=="Brazil")

## Example RasterLayer
r <- raster(nrow=1e3, ncol=1e3, crs=proj4string(SPDF))
r[] <- 1:length(r)

## crop and mask
r2 <- crop(r, extent(SPDF))
r3 <- mask(r2, SPDF)

## Check that it worked
plot(r3)
plot(SPDF, add=TRUE, lwd=2)

    L'écrêtage raster à l'aide de shapefile dans R, mais en gardant la géométrie de la shapefile

    A parfaitement fonctionné! Je vous remercie.

    OriginalL'auteur Josh O'Brien