I need to visualize some area using values provided in HDF-EOS version 5 file and create map using those values and other GIS layers. Unfortunately, my QGIS do not recognize the HDF5 file.
Do I really need to compile GDAL manually as described at: SciGeo?
Maybe there's an easier way to do that? I thought about conversion to the GeoTIFF. However HEG tool crashes during conversion.
I opened my product in HDFView, but I don't see any option to generate GeoTIFF in this application.
Last, but not least, I've tried R! I read appropriate data set, but I don't know what to do next. Is there any possibility to save this data in different format? Or maybe I am able to (how? I can't find this in h5 package documentation)
My current code is:
pacman::p_load(h5)
# download HDF-5 file
destfile <- 'SMAP.h5'
if(!file.exists(destfile))
{
download.file('ftp://n5eil01u.ecs.nsidc.org/SAN/SMAP/SPL2SMAP.003/2015.07.07/SMAP_L2_SM_AP_02293_D_20150707T084111_R13080_001.h5',destfile)
}
# read file
file <- h5file(destfile ,'r')
dataset <- file["Soil_Moisture_Retrieval_Data"]["soil_moisture"]
# @Spacedman: Does your R have a GDAL that supports these files?
subset(rgdal::gdalDrivers(),grepl("HDF",long_name))
# it outputs:
# name long_name create copy
# 51 HDF4Image HDF4 Dataset TRUE FALSE
# 53 HDF5Image HDF5 Dataset FALSE FALSE
No comments:
Post a Comment