ポリゴンを使ってレンガをカットする方法は？

Question

私はbioclim変数のレンガファイルを持っています、レンガは4つの30秒タイルレンガからマージされていましたので、少し大きいです。研究領域のレンガファイルをポリゴンを境界として切断して取得したいと考えています。私は何をすべきか？それ以外の場合は、レンガではできない場合は、ラスターでもできますか？あなたは小さな長方形にレンガをトリミングしたい場合は、事前に

感謝〜

マルコ

Answer 1

はextent()をチェックしてください。また、drawExtent()をクリックして選択する場合は、

EDIT：「カット」と「マスク」という用語を使用しているため、正しく理解しているかどうかはわかりませんが、ここでは役立つ2つの方法があります。両方を使用することもできます。

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