Re: large scale universal kriging for ppt data

[Edzer J. Pebesma]

Hi Femke,

What you need for each time step is a column file with
x, y coordinates, precipitation, and elevation. You can
figure the command file from the manual. If you fix the
variogram, it'll be easy; if you want to refit the variogram
each time step, more work will be involved.
Looping can be done with a shell script, or, what
I would do, is use the gstat package within R (or the
library within S-PLUS); no asciigrids needed there.
Hope this helps,
--
Edzer
femke wrote:

Hello gstat users,
 
I have a question regarding some large scale geostatistics.  Briefly, 
I have 1000 ascii grids of hourly precipitation data (49 points, the 
rest of it nodata), which I want to interpolate.  Given the comment in 
the message archives on driftless kriging, I figure I want to do 
Universal Kriging with a DEM (in asciigrid format) as a mask.  I want 
to iterate over these grids, fit the best model given the weighted 
least squares statistic, interpolate, and spit out a surface (asciigrid).
 
Reading the user's manual it isn't entirely clear to me how gstat 
would accomplish this.  Could someone give me some insight into the 
best way to do this?   
 
Cheers,
 
femke