4つの1バンドのラスター(同じ解像度、同じ範囲...)があり、それらを組み合わせて4つのバンドのラスターを作成します。誰でもRでこれを行う方法を教えてもらえますか?ありがとうございました。
3690 次
1 に答える
3
同じ範囲と解像度のラスターのラスタースタックを作成するには、ラスターパッケージのスタック関数を使用します。
s <- stack(raster1, raster2)
この場合、ラスターはラスターオブジェクトまたはラスターへのファイルパスにすることができます。
編集(具体例):
library(raster)
fn <- system.file("external/test.grd", package="raster")
s <- stack(raster(fn), raster(fn)*2) #here i had to create the raster object since I
#was multiplying one of the input rasters
s
class : RasterStack
dimensions : 115, 80, 9200, 2 (nrow, ncol, ncell, nlayers)
resolution : 40, 40 (x, y)
extent : 178400, 181600, 329400, 334000 (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:28992
+towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812 +proj=sterea
+lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000
+ellps=bessel +units=m +no_defs
names : test.1, test.2
min values : 128.434, 256.868
max values : 1805.78, 3611.56
于 2013-03-26T17:22:23.223 に答える