组合一个波段栅格以获得多波段光栅

时间:2013-03-26 15:55:57

标签: r raster

我有4个单波段的栅格(相同的分辨率,相同的范围...),我想要结合使用四波段的光栅。任何人都可以告诉我如何用R做到这一点? 谢谢。

1 个答案:

答案 0 :(得分:2)

要创建具有相同范围和分辨率的栅格的栅格堆栈,请使用栅格包中的堆栈功能:

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