R: Watershed inlet and outlet nodes: case 2

Watershed.IOR2

R Documentation

Watershed inlet and outlet nodes: case 2

Description

The function determines the inlet and outlet nodes for zhyd watershed objects. This case 2 is for those watersheds that its river inlet and outlet object is length 2 (length(riverIO)=2).

Usage

Watershed.IOR2(x, dist, node)

Arguments

`x`	An object `"SpatialPointsDataFrame"` or `"SpatialPoints"` as are described in package `sp` over the function will search the inlet and outlet nodes of the watershed.
`dist`	A vector with the distances of each point in `x` to the current zhyd boundary.
`node`	An object `"SpatialPointsDataFrame"` as are described in package `sp` over the function will search the inlet and outlet nodes of the watershed. It must be the entire node search object.

Value

A list of length 2:

`inlet`	A `"SpatialPointsDataFrame"` that represents the inlet node of the current `zhyd`.
`outlet`	A `"SpatialPointsDataFrame"` that represents the outlet node of the current `zhyd`.

Note

If there are not inlet or outlet node of the current zhyd is returned 0.

Author(s)

J.A. Torres

Examples

library(Watersheds)
data(WatershedsData)

station1 = WatershedsData$station
subbasin1 = WatershedsData$subbasin
zhyd1 = WatershedsData$zhyd
river1 = WatershedsData$river
node1 = WatershedsData$node	

station1 = SpatialPoints(coords=cbind(4330341.36,3284797.06), 
	proj4string=slot(subbasin1,"proj4string"))
watershed = new("Watershed",station=station1,subbasin=subbasin1,
	zhyd=zhyd1,river=river1,c1=subbasin1,node=node1)

a = Watershed.Order(watershed)
c1 = a[[1]]
nodeIO = a[[9]]			
c1_river = a[[10]]
c1_node = a[[11]]	
		
# determining inlet and outlet watershed nodes
	# determining distances of nodeIO to c1
	boundary = gBoundary(c1)
	dist = gDistance(nodeIO, boundary, byid =TRUE)
	a = Watershed.IOR2(x=nodeIO, dist=dist, node=c1_node)
	str(a)		
	c1_inlet = a$inlet; c1_inlet
	c1_outlet = a$outlet; c1_outlet
	
plot(c1,col="gray60")			
plot(station1,pch=24, bg="blue",add= TRUE)
plot(c1_river, col="blue", add=TRUE)
plot(c1_outlet,pch=21, bg="red",add= TRUE)
plot.PointAttribute(c1_outlet,"ELEV",700,0.8)
title(main="Watershed outlet, case II")