r - ポリゴンを使用してレンガをカットするには?

Question

私は bioclim 変数のブリック ファイルを持っています。ブリックは 4 つの 30 秒タイル ブリックからマージされたので、少し大きくなっています。ポリゴンを境界にして研究領域のブリックファイルを切り取りたいです。私は何をすべきか？そうでなくて、ブリックでは無理ならラスターでもできるのでしょうか？

よろしくお願いします〜

マルコ

Answer 1

extent()レンガを小さな長方形に切り取る場合は、チェックしてください。drawExtent()クリックして選択する場合も同様です。

編集：「カット」と「マスク」という用語を使用したため、正しく理解できたかどうかはわかりませんが、役立つ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)