[R] help on rgdal package

[Hana Sevcikova]

I started to use the rgdal package in order to work with DEM files. Since
these files are based on non-rectangular quads, there are lots of
NA-values in the corresponding datasets. Does anybody know how to remove
the NA-values from an object of class GDALDataset in order to get a
rectangular dataset?

Also, the function saveDataset doesn't seem to work for these kind of data
(or at least for me). I get the following error:

> x <- new("GDALDataset", "xxx.dem")
> getDriverLongName(getDriver(x))
[1] "USGS Optional ASCII DEM"

> saveDataset(x,"xxx.data")
Error in .local(.Object, ...) :
        GDAL Error 6: GDALDriver::Create() ... no create method implemented for this format.

Or should the dataset be raster data? Then I get:

> rd<-getRasterData(x)
> saveDataset(rd,"test.data")
Error in .assertClass(dataset, "GDALReadOnlyDataset") :
        Object is not a member of class GDALReadOnlyDataset

Does anybody see what am I doing wrong? Thanks.

Hana Sevcikova