Проблема преобразования данных в формат hdf в формат tiff - PullRequest
Новая
       96

Проблема преобразования данных в формат hdf в формат tiff

0 голосов
/ 11 февраля 2020

У меня есть данные hdf5, и я хочу преобразовать их в растр, но я получаю сообщение об ошибке в последней строке кода. Ошибка как «Ошибка в (функция (классы, fdef, mtable)): невозможно найти унаследованный метод для функции« растр »для сигнатуры« массив »« »

library(raster)
library(gdalUtils)
library(sp)
library(sf)
library(gdalUtilities)
library(rgdal)      
library(tiff)
library(ncdf4)
library(hdf5r)
library(hdfqlr)
#library(BiocManager)
#BiocManager::install("rhdf5")
#source("http://bioconductor.org/biocLite.R")
#biocLite("rhdf5")
library(rhdf5)
library(gdalUtils)
library(raster)
library(sf)
library(rasterVis)
library(emmeans)
library()
library(lsmeans)
library(biomod2)
library(reshape)
library(reshape2)
setwd("D://2R")
options(max.print = 99999999)

# set path and filename
hdf_path <- "D:/2R/"
hdf_name <- "1C.F19.SSMIS.XCAL2016-V.20150906-S011918-E030110.007356.V05A.hdf5"
hdf_file <- paste(hdf_path, hdf_name, sep="")
h5ls(hdf_file)
# read the attributes of Temperature
h5readAttributes(hdf_file, name = "S1/Tc")
# get lons and lats
lon <- h5read(hdf_file, "/S1/Longitude")
nlon <- length(lon)
nlon
head(lon); tail(lon)
lat <- h5read(hdf_file, "/S1/Latitude")
nlat <- length(lat)
nlat
head(lat); tail(lat)
# read the temperature data
h1 <- h5read(hdf_file, "S1/Tc")
class(h1); str(h1)
t<- raster(h1)
...