INTAMAP--->GeoTIFF

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