如何从HDF5文件制作RasterBrick? [R

时间:2016-02-10 18:49:37

标签: r maps time-series raster hdf5

如何从多个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

1 个答案:

答案 0 :(得分:0)

一种方法是使用gdalUtilshdf5文件转换为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=""))}