intamap (version 1.3-10)

blockPredict: Spatial block prediction

Description

blockPredict is a generic method for prediction of spatially aggregated variables within the intamap-package package.

Usage

blockPredict(object, ...)

Arguments

object

a list object of the type described in intamap-package

...

other arguments that will be passed to the requested interpolation method. See the individual interpolation methods for more information. The following arguments from object can be overrun through ...: [object Ob

Value

• a list object similar to object, but extended with predictions at a the set of locations defined object.

code

gstat

Details

The function blockPredict is a wrapper around the spatialPredict.block function within the intamap-package package, to simplify the calls for block predictions. Block predictions are spatial predictions assumed to be valid for a certain area. The blocks can either be given by passing SpatialPolygons as the predicitonLocations or by passing the block-argument through the parameters of the object or through the ...-argument. There are esentially two ways to solve the problems of block predictions.

• analytical

{block predictions can be found directly by block kriging} numerical{block predictions can be found through numerical simulations over a set of points within the block, the requested output is found by averaging over these simulations }

References

http://www.intamap.org/

See Also

gstat,autoKrige

Examples

# This example skips some steps that might be necessary for more complicated
# tasks, such as estimateParameters and pre- and postProcessing of the data
data(meuse)
coordinates(meuse) = ~x+y
meuse$value = log(meuse$zinc)
data(meuse.grid)
gridded(meuse.grid) = ~x+y
proj4string(meuse) = CRS("+init=epsg:28992")
proj4string(meuse.grid) = CRS("+init=epsg:28992")

# set up intamap object:
obj = createIntamapObject(
	observations = meuse,
	predictionLocations = meuse.grid,
	targetCRS = "+init=epsg:3035",
	class = "automap"
)

# do interpolation step:
obj = conformProjections(obj)
obj = estimateParameters(obj)  # blockPredict
obj = blockPredict(obj,block=c(100,100))  # blockPredict