Я хочу получить подмножество данных со спутника CALIPSO.
Я скачал плитку из центра обработки данных Langley и нашел эту ссылку о том, как читать в формате HDF, в R: http://geoscripting -wur.github.io / IntroToRaster / hdf.html и эту запись в stackoverflow Чтение hdf файлов в R и преобразование их в растры geoTIFF
library(gdalUtils)
library(raster)
# Provides detailed data on hdf4 files
tile <- "somepath/CAL_LID_L1-ValStage1-V3-40.2019-04-01T09-04-47ZD.hdf"
test <- gdalinfo(tile)
Часть продукции:
...
[84] " SUBDATASET_21_NAME=HDF4_SDS:UNKNOWN:\"somepath/CAL_LID_L1-ValStage1-V3-40.2019-04-01T09-04-47ZD.hdf\":20"
[85] " SUBDATASET_21_DESC=[63255x583] Total_Attenuated_Backscatter_532 (32-bit floating-point)"
...
sds <- get_subdatasets(tile)
gdalinfo дал мне более 130 записей, тогда как для sds я получаю только 67. Тем не менее, gdalinfo обычно дают странный и «реальный» номер, имя которого совпадает с 67 записями в sds.
Теперь ссылка выше говорит, что я должен преобразовать его в TIF:
name <- sds[21]
outfilename <- "somepath/CAL_LID_L1-ValStage1-V3-40.2019-04-01T09-04-47ZD.tif"
gdal_translate(sds[1], dst_dataset = outfilename)
# Load the Geotiff created into R
r <- raster(outfilename)
#or more straight away:
r <- raster(sds[20])
Я тоже пытался
plot(r)
однако, вывод не имеет никакого смысла, и размеры не соответствуют. Я также пробовал некоторые другие подмножества
Прежде всего мне нужно уметь правильно читать данные.
Я хотел бы установить тайлы таким образом, чтобы у меня были только измерения с определенного периода времени.
Возможно, кто-то работал с CALIOP Data раньше или сталкивался с подобной проблемой и может помочь?