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

[Sebastian P. Luque]

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,

>> If you have the cell centre offset for the lower left cell, then:

>> is.na(ice$band1) <- ice$band1 < 0 vice <- ice$band1 mice <-
>> matrix(vice, ncol=760, byrow=TRUE) mice1 <- t(mice[1120:1,]) grd <-
>> GridTopology(c(-4941217, -4941217), c(10000,10000), c(760, 1120))

> With:

> grd <- GridTopology(c(-((760/2)*10000) + 5000, -((1120/2)*10000) +
> 5000), c(10000,10000), c(760, 1120))

> that is assuming that the input grid is centred at the pole, and that
> the cell sizes are 10km by 10km, we get much closer, but the whole
> grid is shifted South. Once we establish the metadata by reverse
> engineering, setting up your workflow ought to be feasible.

Hi Roger,

Yes, that is correct, and I should've posted this to begin with, from
OSI-SAF's documents:

-------------- next part --------------
A non-text attachment was scrubbed...
Name: OSI-SAF_nh_grid.png
Type: image/png
Size: 26705 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20090125/cf6cb95b/attachment.png>
-------------- next part --------------

and:

-------------- next part --------------
A non-text attachment was scrubbed...
Name: OSI-SAF_nh_hdf5.png
Type: image/png
Size: 44635 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20090125/cf6cb95b/attachment-0001.png>
-------------- next part --------------

but I'm unable to access the Bx and By components in the HDF5 files.


Cheers,

-- 
Seb