Как использовать начало координат для проецирования растра из файла h5 в R?

Я хотел бы проецировать data.frame из файла h5 в проектируемый растр. Файл h5, который я использую, является сеткой, эквивалентной снежно-водной, для Сьерра-Невады, Калифорния, а начало координат находится на одном из подуровней файла. Я не могу понять, как создать спроектированную растровую сетку, используя источник. Ниже приведен мой код, который загружает данные и извлекает структуру файла h5. Я использую R с библиотекой rhdf5

Как мне взять координированную информацию и использовать ее для проецирования моего datSWE, чтобы он представлял собой проецируемый растр?

Вот несколько дополнительных веб-ссылок, которые я использовал для создания этого первого раздела кода.

https://rpubs.com/fenclmar/98163

https://github.com/NEONScience/NEON-Data-Skills/blob/master/tutorials-in-development/R/R-RS-data-institute/hyperspectral/subset-h5-file-R.Rmd

# Call the R HDF5 Library, install from BioConductor repo
library(rhdf5)

#download the data
theUrl<-"ftp://ftp.snow.ucsb.edu/pub/org/snow/products/reconstruction/sierra/reconstruction_sierra_500m_CY2001.h5"
theFile<-"reconstruction_sierra_500m_CY2001.h5"
try(curl::curl_download(theUrl, theFile, quiet=TRUE))

# View structure of file
f <- "reconstruction_sierra_500m_CY2001.h5"
h5ls(f)

#look at how many "levels" of nesting 
h5ls(f,recursive=5)
fiu_struct <- h5ls(f,recursive=5)

# have a look at the structure.
fiu_struct
fiu_struct[3,1]
## Let's view the metadata for the coordinate info
coordInfo  <- h5readAttributes(f,fiu_struct[3,1])
coordInfo

##Extract the snow water equivalent data
datSWE <- h5read(f, "Grid/swe")
str(datSWE)
day1<-datSWE[,,1] #extract first day of the time series grid

0 ответов

Другие вопросы по тегам