netProp-deprecated: Propagate values along a river network. This function does...

In rtop: Interpolation of Data with Variable Spatial Support

Description

Pass values along a river network when the river network has more segments than the prediction polygons.

Usage

netProp(network, from = "FROMJCT", to = "TOJCT", pred = "pred", 
        iprint = 1)

Arguments

network	object of class SpatialLinesDataFrame describing the river network
from	name of the column giving the endpoint ID of each line segment
to	name of the column giving the start ID of each line segment
pred	name of column with predictions
iprint	if iprint >= 1 the function will give some information about the convergence of the value propagation. Use iprint = 0 to suppress this output.

Value

The function will propagate the predictions upwards along the river network. The result is a SpatialLinesDataFrame with predictions for all line segments, which can easier be plotted.

Note

This function works when the topology of the river network is similar to the example here, that the from-column is always the upstream part of a river segment, and that all segments are actually connected.

Author(s)

Jon Olav Skoien

Examples

## Not run: 
library(rgdal)
rpath = system.file("extdata",package="rtop")
observations = readOGR(rpath,"observations")
predictionLocations = readOGR(rpath,"predictionLocations")
observations$obs = observations$QSUMMER/observations$AREASQKM

# Setting some parameters 
params = list(geoDist = TRUE, rresol = 25, cloud = FALSE, model = "Sph")
# Build an object
rtopObj = createRtopObject(observations,predictionLocations, 
              formulaString = obs~1, params = params)
# Fit a variogram (function also creates it)
rtopObj = rtopFitVariogram(rtopObj)
# Check the variogram fit
rtopObj = checkVario(rtopObj, cloud = TRUE, identify = TRUE)
# Predicting at prediction locations
rtopObj = rtopKrige(rtopObj)
# Cross-validation
rtopObj = rtopKrige(rtopObj,cv=TRUE)
cor(rtopObj$predictions$observed, rtopObj$predictions$var1.pred)

rnet = readOGR(".", "rnet")
pred = rtopObj$predictions
rnet$pred = 
    pred$var1.pred[match(rnet$TOJCT, pred$JCTID)]

# will only plot for a few discontinous river segments
spplot(rnet, "pred", col.regions = bpy.colors())
rnet = netProp(rnet)
# will show a prediction for all segments
spplot(rnet, "pred", col.regions = bpy.colors())

 
## End(Not run)

rtop documentation built on May 2, 2019, 6:48 p.m.