2013-03-26 68 views
2

我有4個單波段的光柵(相同的分辨率,相同的程度...),我想結合起來有一個四波段的光柵。任何人都可以告訴我如何用R來做到這一點? 謝謝。結合一個帶光柵獲得多波段光柵

+2

你看着'光柵:: stack'? – johannes 2013-03-26 16:02:27

+0

是與R的矩陣,sp的SpatialGridDataFrame,柵格的柵格,還是什麼? – mdsumner 2013-03-26 20:00:19

回答

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 
+0

堆棧的問題是我獲得了4層的柵格堆棧,而我想要的是帶有4個帶的單個柵格層。具體來說,我有一個只有紅色樂隊,另一個藍色樂隊和第三個綠色樂隊。我想將它們組合成一個RGB光柵圖層。 – 2013-03-27 09:32:43

+0

堆棧中的圖層可以像光柵一樣工作,就像任何多波段光柵一樣,您可以在arcgis中打開tif堆棧,併爲RGB通道分配波段進行測試。在R中,如果要顯示某個堆棧的某些波段作爲RGB組合,則可以在光柵包中使用plotRGB函數:plotRGB(stack_raster,r = 1,g = 2,b = 3,stretch ='hist' );數字代表堆棧層索引。 – 2013-03-27 17:37:14

+0

謝謝,它的作品! – 2013-03-28 09:12:48