The gstat manual (appendix A.3) mentions the possibility to do block kriging
for irregularly shaped blocks, using the area command/option.
It mentions that the 'file' should contain the points that describe the
shape. I'm not sure what it expects. Is it :
1). A "mask" map that has a value for all cells within the irregular block
and No Data for all other cells
2). The vertices of the irregular shape specified by x,y coordinates in a
GeoEas type data file.
In the latter case, should the number of discretisation points within the
block be specified by using the set nblockdiscr=x option?
Marc Hoogerwerf
ALTERRA
Green World Research
POBox 47, 6700 AA Wageningen
Tel:+31 (317) 474649
Fax: +31 (317) 419000
Web site: http://www.alterra.wageningen-ur.nl
Visiting address: Droevendaalsesteeg 3, Gebouwnummer 101
Email: [EMAIL PROTECTED]
