我有一个光栅堆栈,stk由R中的三个光栅图像组成.这是一个简单的例子
# set up a raster stack with three layers
> library(raster)
> r <- raster(nrows=10,ncols=10)
> r[] <- rnorm(100)
> stk <- stack(r,r,r)
# layer names are set by default
> names(stk)
[1] "layer.1" "layer.2" "layer.3"
我为栅格图层指定了名称:
# set layer names to "one", "two" and "three"
> names(stk) <- c('one','two','three')
> names(stk)
[1] "one" "two" "three"
当我使用以下命令将RasterStack写入GeoTiff(多层)时:
writeRaster(stk,"myStack.tif", format="GTiff")
根据文件名重命名图层(见> names(stk)下文).
当我读入光栅堆栈时:
> stk <- stack("myStack.tif")
# the layer names have been set automatically based …