Учитывая растровые временные ряды (спутниковые изображения в пределах одного и того же пространственного экстента, каждое из которых сделано в разную дату), как я могу добавить дату в качестве имени для каждого слоя / полосы в (Geo) Tiff?Возможно, он недостаточно указан?
Используя R в качестве примера:
library(rgdal)
library(raster)
## example data (taken from https://rspatial.org/raster/rs/1-introduction.html#data)
dir.create('/tmp/data', showWarnings = FALSE)
if (!file.exists('/tmp/data/rs/samples.rds')) {
download.file('https://biogeo.ucdavis.edu/data/rspatial/rsdata.zip', dest = '/tmp/data/rsdata.zip')
unzip('/tmp/data/rsdata.zip', exdir='/tmp/data')
}
## just two raster layers
layers <- c("/tmp/data/rs/LC08_044034_20170614_B1.tif", "/tmp/data/rs/LC08_044034_20170614_B2.tif")
lc_stack <- raster::stack(layers)
## set arbitrary timestamps
lc_stack_z <- setZ(lc_stack,c("2019-09-20", "2019-09-21"))
timeseries_rs <- writeRaster(lc_stack_z, "/tmp/data/rs/timeseries-raster.tiff", "GTiff")
Вывод gdalinfo
выглядит следующим образом (нет информации об имени группы):
[...]
Band 1 Block=1497x1 Type=Float32, ColorInterp=Gray
Min=0.096 Max=0.735
Minimum=0.096, Maximum=0.735, Mean=nan, StdDev=nan
NoData Value=-3.39999999999999996e+38
Metadata:
STATISTICS_MAXIMUM=0.73462820053101
STATISTICS_MEAN=nan
STATISTICS_MINIMUM=0.096417911350727
STATISTICS_STDDEV=nan
Band 2 Block=1497x1 Type=Float32, ColorInterp=Undefined
Min=0.075 Max=0.718
Minimum=0.075, Maximum=0.718, Mean=nan, StdDev=nan
NoData Value=-3.39999999999999996e+38
Metadata:
STATISTICS_MAXIMUM=0.71775615215302
STATISTICS_MEAN=nan
STATISTICS_MINIMUM=0.074839904904366
STATISTICS_STDDEV=nan
Существует Библиотека тегов Tiff , но там я не могу найти какой-либо (семантически подходящий) тег, который бы соответствовал моим потребностям.
Есть ли полезная практика, как добавлять имена слоев / групп в(Geo) Tiff - оптимально, который может быть интерпретирован / прочитан другим программным обеспечением, таким как QGis.