Как отрезать кирпич полигоном?

Question

У меня есть кирпичный файл переменных bioclim, кирпич был объединен с четырьмя кирпичами из 30sec, поэтому он немного велик. Я хотел бы получить кирпичный файл моей исследовательской области, разрезав его, используя полигон в качестве границы. Что мне делать? В противном случае, если это невозможно сделать с кирпичом, могу ли я сделать это с растром?

Заранее спасибо ~

Marco

Answer 1

ЗАКАНЧИВАТЬ extent(), если вы хотите, чтобы обрезать кирпич до меньшего прямоугольника. Также drawExtent(), если вы предпочитаете выбирать, нажав.

EDIT: Поскольку вы использовали термины «вырезать» и «маску», я не уверен, что правильно понял, но вот два способа, которые могут помочь. Вы могли бы использовать и то, и другое.

# an example with dimensions: 77, 101, 3 (nrow, ncol, nlayers) 
myGrid_Brick <- brick(system.file("external/rlogo.grd", package="raster")) 

# a simple polygon within those dimensions 
myTriangle_P <- Polygon(cbind(c(10, 80, 50, 10), c(10, 20, 65, 10))) 
myTriangle_Ps <- Polygons(list(myTriangle_P), "fubar") 
myTriangle_SP <- SpatialPolygons(list(myTriangle_Ps)) 
myTriangle_Ras <- rasterize(myTriangle_SP, myBrick) 

# this will crop a brick to minimal rectangle that circumscribes the polygon 
# extent(myCrop) is smaller than extent(myGrid) but no values are changed 
myCrop_Brick <- crop(myGrid_Brick, myTriangle_SP) 

# while this converts every coordinate that is NA in 
# the mask to become NA in the returned brick 
# while leaving the brick extent unchanged 
myMask_Brick <- mask(myGrid_Brick, myTriangle_Ras)