[R-sig-Geo] readGDAL() and HDF5 files

[Sebastian P. Luque]

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:

> http://spatial.nhh.no/misc/hudson_st.png

> 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