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

2009-01-26 Thread Sebastian P. Luque
On Mon, 26 Jan 2009 10:33:54 +0100 (CET), Roger Bivand wrote: > On Sun, 25 Jan 2009, Sebastian P. Luque wrote: >> On Sun, 25 Jan 2009 21:35:12 +0100 (CET), >> Roger Bivand wrote: >> On Sun, 25 Jan 2009, Roger Bivand wrote: >>> On Sat, 24 Jan 2009, Sebastian P. Luque wrote: [...] > Hi, > Usin

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

2009-01-26 Thread Roger Bivand
On Sun, 25 Jan 2009, Sebastian P. Luque wrote: On Sun, 25 Jan 2009 21:35:12 +0100 (CET), Roger Bivand 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)

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

2009-01-25 Thread Roger Bivand
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 <- GridTop

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

2009-01-25 Thread Roger Bivand
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(1,1)

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

2009-01-24 Thread Sebastian P. Luque
On Sat, 24 Jan 2009 12:56:13 +0100 (CET), Roger Bivand wrote: > On Fri, 23 Jan 2009, Sebastian P. Luque wrote: >> On Sat, 24 Jan 2009 08:11:15 +1100, >> Michael Sumner wrote: >> Hi Sebastian, When you say the geographic coordinates are not a >>> regular grid - is it that the actual grid is Merc

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

2009-01-24 Thread Roger Bivand
On Fri, 23 Jan 2009, Sebastian P. Luque wrote: On Sat, 24 Jan 2009 08:11:15 +1100, Michael Sumner wrote: Hi Sebastian, When you say the geographic coordinates are not a regular grid - is it that the actual grid is Mercator but the NetCDF file stores an X and Y vector separately for each uniqu

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

2009-01-23 Thread Sebastian P. Luque
On Sat, 24 Jan 2009 08:11:15 +1100, Michael Sumner wrote: > Hi Sebastian, When you say the geographic coordinates are not a > regular grid - is it that the actual grid is Mercator but the NetCDF > file stores an X and Y vector separately for each unique longitude and > latitude? I've seen this ma

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

2009-01-23 Thread Michael Sumner
Hi Sebastian, When you say the geographic coordinates are not a regular grid - is it that the actual grid is Mercator but the NetCDF file stores an X and Y vector separately for each unique longitude and latitude? I've seen this many times, but never with enough metadata to determine that with

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

2009-01-23 Thread Sebastian P. Luque
On Fri, 23 Jan 2009 11:01:21 +0100 (CET), Roger Bivand wrote: > On Fri, 16 Jan 2009, Sebastian P. Luque wrote: >> Hi, >> I'm trying to use HDF5 files with this structure: > Sorry for the delay. Having now set up a reprodicable route (OSGeo4W > GDAL and drivers, rgdal built against OSGeo4W as in

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

2009-01-23 Thread Roger Bivand
On Fri, 16 Jan 2009, Sebastian P. Luque wrote: Hi, I'm trying to use HDF5 files with this structure: Sorry for the delay. Having now set up a reprodicable route (OSGeo4W GDAL and drivers, rgdal built against OSGeo4W as in a posting a couple of days ago), I can see the problem, but have no r

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

2009-01-16 Thread Sebastian P. Luque
Hi, I'm trying to use HDF5 files with this structure: ------ $ gdalinfo conc_200901011200.hdf Driver: HDF5/Hierarchical Data Format Release 5 Files: conc_200901011200.hdf Size is 512, 512 Coordinate System is `' Subdatasets: