search for: savedataset

...gdal 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(.Ob...

R-geographers... I'm trying to solve a problem to implement a line-by-line tiled processing using RGDAL (read 1 line of an image, process the one line, write one line of the image to a binary file). Everything except for the final step I'm able to do using a combination of RGDAL and r-base commands. Below is the basic structure, the input file "elev" can be any image file