[R-sig-Geo] INTAMAP--->GeoTIFF

[Jon Olav Skoien]

Hi Tarek,
> output = interpolate(observations = e,
> predictionLocations = grd,
> outputWhat = list(mean=T, variance=T, excprob = 7),
> maximumTime = 30,
> methodName = "automatic",
> optList = list() )
> plot(output)
> #how generate a GeoTIFF ?
You can generate the GeoTIFF with a call to writeGDAL with the results 
of the interpolation, which will be in output$predictions:
 > writeGDAL(output$predictions, fname = "result.tif", drivername = "GTiff")
You will then get a GeoTIFF with two bands - the interpolation and the 
kriging variance. If you are just interested in the interpolation 
result, you can instead call:
 > writeGDAL(output$predictions[,"var1.pred"], fname = "result.tif", 
drivername = "GTiff")


Cheers,
Jon