readGDAL() and HDF5 files
On Mon, 26 Jan 2009 10:33:54 +0100 (CET),
Roger Bivand <Roger.Bivand at nhh.no> wrote:
On Sun, 25 Jan 2009, Sebastian P. Luque wrote:
On Sun, 25 Jan 2009 21:35:12 +0100 (CET), Roger Bivand <Roger.Bivand at nhh.no> wrote:
On Sun, 25 Jan 2009, Roger Bivand wrote:
On Sat, 24 Jan 2009, Sebastian P. Luque wrote:
[...]
Hi,
Using the lower left geographical coordinate:
ll <- c(project(matrix(c(-80.7299, 33.9755), nrow=1), "+init=epsg:3411")) grd <- GridTopology(ll, c(10000,10000), c(760, 1120))
and so on, I get:
which isn't bad. Maybe we need (ll+5000) instead, but it's hard to tell. The assumption in project() is no datum shift but +ellps=WGS84, so that might need attention too.
Should we draw these problems to the data providers' attention (it is the Norwegian Met. Office, and there are R people there too)? Perhaps they have a FAQ?
Yes, in fact I sent them the relevant info over the weekend, to which they replied today and will investigate. And sure enough, they pointed me to an R package they built that works with a custom HDF5 library for their files: http://saf.met.no/p/tools.html The package hasn't been updated for 2 years, so they will see what needs to change, based on our exchanges here. Cheers,
Seb