У меня есть растровый стек 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 on the filename
# they should be "one", "two" and "three"
> names(stk)
[1] "myStack.1" "myStack.2" "myStack.3"
Знаете ли вы, как сохранить имена слоев при записи RasterStacks в R? Я попытался записать стек в форматы GeoTIFF и NetCDF.
Спасибо, Кевин