我正在研究如何在R中处理.h5(hdf)VIIRS DNB栅格数据。我想绘制它或将其导出为地理位置。
到目前为止,我可以使用
读取文件hdf<- h5file("mypath/GDNBO-SVDNB_npp_d20171101_t0110370_e0116174_b31151_c20180208224859630066_nobc_ops.h5", mode = "r")
但我无法弄清楚如何进入“光辉”乐队,更不用说绘制它了。我已经阅读了“h5”软件包文档而没有任何进展..
作为另一种选择,我尝试使用
将h5文件导出到geotiffsds <- get_subdatasets("subset_1_d20171101_t0110370.h5")
gdal_translate(sds[17], dst_dataset = "radiance.tif", overwrite=T) #subset 17 is the radiance band i need
然而,在转换过程中,当在例如地理位置打开地理位置时,文件已经失去了地理位置。 QGIS,它没有任何位置或投影。
任何人都知道如何在R?
中处理这些类型的文件