I would be happy to look into this, but the file you point to requires authentification.
On 8/22/19 5:52 AM, Frederico Faleiro wrote: > Dear list, > > I am using some netCDF files from the CMIP5 climate models in raster > package, but I am having some issues with one model. > The netCDF file from the GFDL-ESM2G model (e.g. > http://aims3.llnl.gov/thredds/fileServer/css03_data/cmip5/output1/NOAA-GFDL/GFDL-ESM2G/rcp45/mon/atmos/Amon/r1i1p1/v20120412/pr/pr_Amon_GFDL-ESM2G_rcp45_r1i1p1_204101-204512.nc > ) > have the error message as in bellow example. > > #Example > library(raster) > s1 <- stack("pr_Amon_GFDL-ESM2G_rcp45_r1i1p1_204101-204512.nc") > Error in .rasterObjectFromCDF(x, type = objecttype, band = band, ...) : > cells are not equally spaced; you should extract values as points > > # I check some solutions that recomend force read the file with the > argument "stopIfNotEqualSpaced = F" as bellow. > sf <- stack("pr_Amon_GFDL-ESM2G_rcp45_r1i1p1_204101-204512.nc", > *stopIfNotEqualSpaced > = F*) > bf <- brick("pr_Amon_GFDL-ESM2G_rcp45_r1i1p1_204101-204512.nc", > *stopIfNotEqualSpaced > = F*) > > I performed some tests and only "works" with brick. However I did not find > any solution to check where is the problem and fix it in the file. > > How can I check if it is really an issue and fix it? > > sessionInfo() > R version 3.5.1 (2018-07-02) > Platform: x86_64-w64-mingw32/x64 (64-bit) > Running under: Windows 10 x64 (build 18362) > > Matrix products: default > > locale: > [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 > [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C > [5] LC_TIME=Portuguese_Brazil.1252 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] raster_2.8-19 sp_1.3-1 > > loaded via a namespace (and not attached): > [1] compiler_3.5.1 rgdal_1.4-3 Rcpp_1.0.1 codetools_0.2-15 > ncdf4_1.16.1 > [6] grid_3.5.1 lattice_0.20-35 > > Best regards, > > Frederico > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > -- Edzer Pebesma Institute for Geoinformatics Heisenbergstrasse 2, 48151 Muenster, Germany Phone: +49 251 8333081
pEpkey.asc
Description: application/pgp-keys
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo