femke
Tue, 10 Feb 2004 11:42:52 -0800
|
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
|