Suite

Lire les sous-variables netcdf et les convertir en grille raster dans R ?


je veux ouvrir lefsvariable qui est sousRécupérationcatégorie et convertissez-la en grille raster.

fichier de bibliothèque (ncdf4) <- "http://co2web.jpl.nasa.gov/thredds/fileServer/OCO-2/B7101Ar_r02/2014/09/06/LtCO2/oco2_LtCO2_140906_B7101Ar_150909171225s.nc4" var <- nc_open(file," Récupération/fs")

Le serveur à https://co2web.jpl.nasa.gov/thredds/catalog/OCO-2/B7101Ar_r02/2014/09/06/LtCO2/catalog.html?dataset=nasa.jpl.oco2/B7101Ar_r02/2014/09 /06/LtCO2/oco2_LtCO2_140906_B7101Ar_150909171225s.nc4 ne semble pas les servir via OpenDAP de sorte que ncdf4 puisse les lire à distance. Essayez de le copier localement et de faire l'extraction avec :

fn <- 'oco2_LtCO2_140906_B7101Ar_150909171225s.nc4' nc <- nc_open(fn) vn <- 'Retrieval/fs' ncatt_get(nc,vn) x <- ncvar_get(nc,vn) > system(paste(')md5',fn MD5 (/tmp/oco2_LtCO2_140906_B7101Ar_150909171225s.nc4) = aba9389fe3bb0cf9dfe321c24fc211b4 > résumé(x) Min. 1er Qu. Médiane Moyenne 3e Qu. Max. NA 0,000 0,074 0,165 0,270 0,317 11,600 17782

Un extrait de la sortie 'nc' :

… 65 variables (hors variables de dimension) : float latitude[sounding_id] (stockage contigu) unités : degrés_north long_name : latitude missing_value : -999999 commentaire : latitude centrale de la mesure float longitude[sounding_id] (stockage contigu) unités : degrés_east long_name : longitude missing_value : -999999 commentaire : longitude centrale de la mesure double time[sounding_id] (stockage contigu) unités : secondes depuis 1970-01-01 00:00:00 long_name : heure missing_value : -999999…

S'ils le servaient en tant qu'OpenDAP, cela pourrait fonctionner à distance comme :

uri <- 'http://geoport.whoi.edu/thredds/dodsC/examples/bora_feb.nc' nc2 <- nc_open(uri) nc2