如何从 HDF5 文件制作 RasterBrick? R
How to Make a RasterBrick from HDF5 Files? R
如何从多个 hdf5
文件中制作 Rasterbrick
in R
?通常,数据以 hdf5
格式提供,必须将其转换为更友好的格式以便于处理。
目前我知道 rhdf5 package
但我不确定如何获得 RasterBrick
。
source("http://bioconductor.org/biocLite.R")
biocLite("rhdf5")
library("rhdf5")
library("raster")
上访问多个 hdf5 文件
你可以用两个文件来说明。
谢谢!
AT.
一个选项是使用 gdalUtils
将 hdf5
文件转换为 GTiff
。一旦你这样做了,你就可以在一个堆栈中阅读它们。这是一个示例代码:
# list all the `hdf5` files
files <- list.files(path=".", pattern=paste(".*.h5",sep=""), all.files=FALSE, full.names=TRUE)
#choose the band that you want using the sds[] option and write GTiff files.
for (i in (files)) {
sds <- get_subdatasets(i)
r2 <- gdal_translate(sds[1], dst_dataset =paste(i,".tif",sep=""))}
如何从多个 hdf5
文件中制作 Rasterbrick
in R
?通常,数据以 hdf5
格式提供,必须将其转换为更友好的格式以便于处理。
目前我知道 rhdf5 package
但我不确定如何获得 RasterBrick
。
source("http://bioconductor.org/biocLite.R")
biocLite("rhdf5")
library("rhdf5")
library("raster")
上访问多个 hdf5 文件
你可以用两个文件来说明。
谢谢!
AT.
一个选项是使用 gdalUtils
将 hdf5
文件转换为 GTiff
。一旦你这样做了,你就可以在一个堆栈中阅读它们。这是一个示例代码:
# list all the `hdf5` files
files <- list.files(path=".", pattern=paste(".*.h5",sep=""), all.files=FALSE, full.names=TRUE)
#choose the band that you want using the sds[] option and write GTiff files.
for (i in (files)) {
sds <- get_subdatasets(i)
r2 <- gdal_translate(sds[1], dst_dataset =paste(i,".tif",sep=""))}