[SOLVED] Polygon Block Kriging

Here you can suggest new features to be added to RGeostats package.

[SOLVED] Polygon Block Kriging

Postby Maxime Beauchamp » Wed Nov 06, 2013 2:32 pm

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
Maxime Beauchamp
 
Posts: 29
Joined: Thu Jan 17, 2013 2:48 pm

Re: Polygon Block Kriging

Postby Didier Renard » Sun Nov 10, 2013 9:05 pm

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.
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm


Return to Global suggestions

Who is online

Users browsing this forum: No registered users and 7 guests

cron