View source: R/result.extract.interpolate.R
result.extract.interpolate | R Documentation |
Interpolates result values with given spatial information.
result.extract.interpolate(mask.grid, values, alt, x, y)
mask.grid |
The grid with spatial information the values are ordered by. |
values |
The values which should be interpolated. |
alt |
The related altitude for the gridcells of ‘mask.grid’. |
x |
The related Rechtswert (Gauss-Krueger-coordinates) for the gridcells of ‘mask.grid’. |
y |
The related Hochwert (Gauss-Krueger-coordinates) for the gridcells of ‘mask.grid’. |
Interpolates result values with given spatial information by external drift kriging.
A vector with the interpolated values.
Daniel Doktor, Maximilian Lange
Krige, D.G., 1951. A statistical approach to some basic mine valuation problems on the witwatersrand. Journal of the Chemical, Metallurgical and Mining Society of South Africa 52, 119-139. Pebesma, E.J., 2004. Multivariable geostatistics in S: the gstat package. Computers & Geosciences 30, 683-691.
result.extract.main
,data.main
## load preprocessed data data(dataFinal) ## load spatial information data(relatedGrid) ## set or load params params <- c(0, 0.058326, 0.109494, 0.039178, -10.34, -0.89, 18.11,-10.03, 28.61, 44.49) ## apply model result <- pim.solve(params, dataFinal, model.no=11, silent=FALSE, out2File=FALSE) ## resolve outlier information outliers <- result$outlier.bb + result$outlier.lc outliers.na <- which(is.na(outliers)==TRUE) outliers[outliers.na] <- rep(0, length(outliers.na)) mask.grid <- relatedGrid ## extract valid modelled values values.model <- result.extract.sub(mask.grid=mask.grid, result$doy.bb.pim, result$gk4.x, result$gk4.y, outliers=outliers, silent=FALSE, withOutliers=FALSE)$values ## interpolate result values with spatial informations of mask.grid values.model <- result.extract.interpolate(mask.grid=mask.grid, values=values.model, alt=mask.grid$alt, x=mask.grid$x, y=mask.grid$y)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.