[R-sig-Geo] Strange spatial reference system netCDF [stars and sf functions not working now]

maurizio marchi m@ur|z|o@m@rch| @end|ng |rom |bbr@cnr@|t
Tue Feb 18 15:08:05 CET 2020


Dear all,
Thanks for your time. Unfortunately the solution the community (i.e. 
Edzer Pebesma, https://github.com/r-spatial/stars/issues/175) is not 
working anymore. Does anyone knows why the code is now broken? Here the 
output of the old R code...

/f<-"tasmax_rcp85_land-rcm_eur_12km_01_mon_198012-208011.nc"//
//library(stars)//
//  Loading required package: abind//
//  Loading required package: sf//
//  Linking to GEOS 3.8.0, GDAL 3.0.2, PROJ 6.2.1//
//r<-read_ncdf(f, ncsub = cbind(start = c(1, 1, 1, 1), count = c(418, 
406, 3, 1)), eps=1e-3)//
//  no 'var' specified, using tasmax//
//  other available variables://
//  rotated_latitude_longitude, ensemble_member, time, time_bnds, 
grid_latitude, grid_latitude_bnds, grid_longitude, grid_longitude_bnds, 
ensemble_member_id, month_number, year, yyyymm//
//  Error in UseMethod("GPFN") : //
//  no applicable method for 'GPFN' applied to an object of class 
"rotated_latitude_longitude"//
//rx<-read_stars(f, proxy = TRUE) # only for the crs!//
//  Warning messages://
//  1: In CPL_read_gdal(as.character(x), as.character(options), 
as.character(driver),  ://
//    GDAL Message 1: The dataset has several variables that could be 
identified as vector fields, but not all share the same primary 
dimension. Consequently they will be ignored.//
//  2: In CPL_read_gdal(as.character(x), as.character(options), 
as.character(driver),  ://
//    GDAL Message 1: dimension #1 (time) is not a Time dimension.//
//  3: In CPL_read_gdal(as.character(x), as.character(options), 
as.character(driver),  ://
//    GDAL Message 1: dimension #0 (ensemble_member) is not a Time 
dimension.//
//st_crs(r)<-st_crs(rx)//
//  Warning messages://
//  1: In CPL_crs_from_proj4string(x) ://
//    GDAL Error 1: PROJ: proj_as_wkt: DerivedGeographicCRS can only be 
exported to WKT2//
//  2: In CPL_crs_from_proj4string(x) ://
//    GDAL Error 1: PROJ: proj_as_wkt: DerivedGeographicCRS can only be 
exported to WKT2//
//r0 = stars:::st_transform_proj.stars(r, 4326)//
//  Error: invalid crs: , reason: generic error of unknown origin//
//  proj_create: unrecognized format / unknown name//
//r0//
//Error: object 'r0' not found//
/
All my bests
-- 
*Maurizio Marchi,
PhD Forest Science - Ecological Mathematics*
Researcher
CNR - Institute of Biosciences and BioResources (IBBR), Florence 
division (Italy)
SkypeID: maurizioxyz
http://ibbr.cnr.it/ibbr/info/people/maurizio-marchi
#####------#####
EUFGIS National Focal Point for Italy (www.eufgis.org)
Scopus Author ID: 57188626512
ResearcherID: T-3813-2019
http://b4est.eu/ project
/...CPUs don't need to sleep.../

	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list