[R] How to open grib file in R?
b.rowlingson at lancaster.ac.uk
Mon Jan 14 17:52:50 CET 2013
On Mon, Jan 14, 2013 at 4:36 PM, Jonsson <amen.alyaari at bordeaux.inra.fr> wrote:
> but I got this error:
> C:\Users\aalyaari\Desktop\Vol_025_H14_2010060700.grib has GDAL driver GRIB
> and has 721 rows and 1440 columns Error in validityMethod(as(object,
> superClass)) :
> Geographical CRS given to non-conformant data: -90.125 90.125
> Any idea why I am getting this error?
At a guess, and this is just a guess, I'd say the system tried to
give a Geographical CRS to non-conformant data...
So, what's a Geographical CRS? One with latitude and longitude.
What's "non-conformant data"? Data that can't conform to latitude and
longitude. Why is -90.125 not conformable to latitude-longitude?
Because it's further south than the south pole.
Why have you got data that goes further south than the south pole and
further north than the north pole? I don't know. I suspect you have
points at the poles and the software thinks these are grid cells that
are 0.250 lat-long degrees wide. Ask your source.
For what its worth, the raster package will ignore the error, but
Warning: Inside GRIB2Inventory, Message # 5
ERROR: Ran out of file reading SECT0
There were 76 trailing bytes in the file.
class : RasterLayer
band : 1
dimensions : 721, 1440, 1038240 (nrow, ncol, ncell)
resolution : 0.25, 0.25 (x, y)
extent : -0.125, 359.875, -90.125, 90.125 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +a=6367470 +b=6367470 +no_defs
data source : /nobackup/rowlings/Downloads/Vol_025_H14_2010060700.grib
names : Vol_025_H14_2010060700
I suspect there's an off-by-one error in the code that created this.
Remember there are only three possible mistakes in computing - off by
one errors and divide by zero errors.
More information about the R-help