On Mon, 12 Nov 2007 23:29:59 +0100 (CET),
Roger Bivand <[EMAIL PROTECTED]> 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

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to