如何从多个Rasterbrick
文件中R
hdf5
进行hdf5
?通常,数据以rhdf5 package
格式提供,并且必须将其转换为更友好的格式以便于处理。
目前我知道RasterBrick
,但如何获得source("http://bioconductor.org/biocLite.R")
biocLite("rhdf5")
library("rhdf5")
library("raster")
是我不确定的。
import static org.junit.Assert.assertTrue;
您可以在此链接http://mirador.gsfc.nasa.gov/cgi-bin/mirador/cart.pl?C1=GPM_3IMERGHH&CGISESSID=fb3b45e091f081aba8823f3e3f85a7d9&LBT_THRESHOLD=4000000上访问多个hdf5文件。
您可以使用两个文件进行说明。
谢谢!
AT
答案 0 :(得分:0)
一种方法是使用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=""))}