On Thu, 21 Aug 2003, Edzer J. Pebesma wrote: I tried to use the area option mentioned in appendix A.3. But it turned out to be extremely slow. Probably I made a mistake.
Below is my cmd file. It took gstat 2 hours to initialise the maps whereas it does take usually a few seconds. 'mask' and 'area map' are both raster maps on the same grid cell support. For the area map I chose integers classifying each polygon. Each polygon (area) has a different size. There are 31 in total. --------------------------------------------------------- # # gstat command file, Linux version 2.4.2 (26 April 2003) # Fri Jul 25 00:23:02 2003 # data(Pox0t30): 'proplskPox0t30.gstat', x=1, y=2, v=23, mv=-9999, min=10, max=120; variogram(Pox0t30): 95.6177 Nug(0) + 89.5835 Exp(13000.3); mask: 'mask'; predictions(Pox0t30): 'map/Pox0t30_OK_block_LEI.pred'; variances(Pox0t30): 'map/PoxOT30_OK_block_LEI.var'; area: 'LEI'; # blocksize; #do I have to specify block kriging here? ---------------------------------------------------------- Ulrich > Ulrich, is block kriging for non-rectangular blocks what you want? If yes, > appendix A.3 should give you just enough information. Verify your results > with a rectangular block kriging, approximately the size & location of > the non-rectangular area. > -- > Edzer > > Ulrich Leopold wrote: > > >Dear list, > > > >I would like to aggregate point data to non-regular supports, e.g. > >polygons. Is it possible to specify a polygon file or a rasterized > >polygon file via 'mask' to obtain kriging estimates per polygon area? I > >could not figure it quite out from the manual what to do. > > > >Ulrich > > > > > > -- __________________________________________________ Ulrich Leopold MSc. Department of Physical Geography Institute for Biodiversity and Ecosystem Dynamics Faculty of Science University of Amsterdam Nieuwe Achtergracht 166 NL-1018WV Amsterdam Phone: +31 20 525 7456 (7451 Secretary) Fax: +31 20 525 7431 Mobile: +31 64 220 3028 Email: [EMAIL PROTECTED] URL: www.science.uva.nl/ibed/research/Research_Fields/fg/stafffg/index.html
