R'de şekil dosyası kullanarak kırpma raster'i kırpma, ancak şekil dosyası

Question

'un geometrisini koruyorum Düzensiz bir şekil dosyasına (Amazon biome) dayalı olarak bir raster'i kırpmak (veya kesmek) için {raster} kullanıyorum ancak çıktı her zaman dikdörtgen bir boyuta sahip. Ancak, çıktıya şekil dosyasının tam olarak aynı geometrisinde ihtiyacım var. Herhangi bir ipucu? Şerefe.

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

Answer 1

Seçeneklerden biri 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)