Page 1 of 1

[SOLVED] Polygon Block Kriging

PostPosted: Wed Nov 06, 2013 2:32 pm
by Maxime Beauchamp
Hi,

I am back on the forum. I'd like to know if it is possible to perform block kriging on polygons. For instance, I would like to perfom a kriging of pollution levels on French regions. I have shapefiles that give me the polygons associated to each region and I'd like to know their predicted value. Is it possible with RGeoS anf if not, is it planned to implement this kind of block kriging in future developments ?

Kind regards,

Maxime beauchamp
INERIS
DRC/MOCA

Re: Polygon Block Kriging

PostPosted: Sun Nov 10, 2013 9:05 pm
by Didier Renard
Hi Maxime

Welcome back to the forum.
There is a procedure (already implemented) which enables calculations of the mean value of a target variable over a polygon.
I must first define what is a POLYGON in RGeoS. A POLYGON is a set of POLYSET(s). The number of polysets in a polygon is simple positive.
Then for each polyset, the mean value of the variable is calculated using a global Kriging (rather than a block kriging). In particular, it also gives access to the variance of the estimation error. As far as the data are concerned, it offers the possibility to select all the samples included inside the polyset, or to use all available data instead.
Finally this global estimation is carried out using a discretization of the polygon for the calculation of the covariance terms. Thsi discretization can be performed regularly (discertization points located on the nodes of a regular grid, with a mesh fixed by the user) or a stratified grid (where samples are generated randomly within cells of a regular grid defined by the user).

Hope this will serve.