SpatialGridDataFrame to netcdf grid
On Mon, 12 Nov 2007 23:29:59 +0100 (CET),
Roger Bivand <Roger.Bivand at nhh.no> wrote:
[...]
Then the fallback is to write a GTiff, and use gdal_translate to convert it (one of the gdal utilities that build with the GDAL shared object). If the same problem occurs, you may need to install netCDF and GDAL from source. If that doesn't resolve it, then you'll need to reshape the data to the required format and use functions in the ncdf package to create and write to the file.
I found out that the Debian packages with the gdal_translate binary is
gdal-bin, and trying to writeGDAL() to GTiff and then gdal_translate as
suggested earlier gives the same error message about cdfid not being
valid.
Just so I understand how the SpatialGridDataFrame object is structured:
R> gridparameters(locs.sub.ts)
cellcentre.offset cellsize cells.dim
1 -69.12 0.06109 113
2 63.64 0.03424 122
R> bbox(locs.sub.ts)
min max
coords.x1 -69.15 -62.25
coords.x2 63.62 67.80
How are the data (e.g. z values) aligned with respect to these values? I
would like to know what as.data.frame(locs.sub.ts) returns: the
coordinates and data for the midpoint of the cell, or the data for some
other reference point in the cell. Thanks again for any help.
Seb