在R中绘制VIIRS栅格数据(.h5)

时间:2018-02-09 14:56:37

标签: r raster hdf

我正在研究如何在R中处理.h5(hdf)VIIRS DNB栅格数据。我想绘制它或将其导出为地理位置。

到目前为止,我可以使用

读取文件
hdf<- h5file("mypath/GDNBO-SVDNB_npp_d20171101_t0110370_e0116174_b31151_c20180208224859630066_nobc_ops.h5", mode = "r")

但我无法弄清楚如何进入“光辉”乐队,更不用说绘制它了。我已经阅读了“h5”软件包文档而没有任何进展..

作为另一种选择,我尝试使用

将h5文件导出到geotiff
sds <- 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?

中处理这些类型的文件

0 个答案:

没有答案