| merge {GSIF} | R Documentation |
Merges objects of class "SpatialPredictions" or "RasterBrickSimulations" and produces average predictions where the two objects overlap spatially. If the predictions are available at different resolutions, then it downscales all other grids to the smallest grid cell size using bicubic splines (for predictions) i.e. nearest neighbor algorithm (for simulations). Weigths can be passed via the RMSE.l argument, otherwise they will be estimated from validation slot (if objects are of the class "SpatialPredictions").
## S4 method for signature 'SpatialPredictions,SpatialPredictions' merge(x, y, ..., RMSE.l = NULL, silent = TRUE)
x |
object of class |
y |
object of class |
... |
additional objects of class |
RMSE.l |
numeric; list of mean prediction errors for each object (these are used as weights during the averaging) |
silent |
logical; specifies whether to print out the progress and used RMSE's |
Returns an object of type "SpatialPixelsDataFrame" or "RasterBrickSimulations" that contains only the merged values.
Merging of multiple spatial predictions using weighted averaging is a heuristic approach to mapping. This method assumes that the predictions are completely independent (independent covariates, independent models), but this not might be the case and hence the merged predictions will be sub-optimal. Merging multiple predictions is however attractive for situations where the predictions do not have the same extent, so that spatial predictions with larger coverage can be used to fill in the gaps in locally produced predictions.
Tomislav Hengl and Gerard B.M. Heuvelink
Heuvelink, G.B.M., Bierkens, M.F.P. (1992) Combining soil maps with interpolations from point observations to predict quantitative soil properties. Geoderma 55(1-2): 1-15.