Watershed.Tributary-methods: S4 Method for Function 'Watershed.Tributary'

In Watersheds: Spatial Watershed Aggregation and Spatial Drainage Network Analysis

Description

S4 Method for function Watershed.Tributary. Definition of the order of tributary zhyd watersheds of the current zhyd watershed.

Value

The method returns a list of 4 objects:

c2c3	A list of length 2 with objects SpatialPolygonsDataFrame of length 1 ordered that represents the greater watershed and second tributary of the current zhyd watershed object.
c2	An object SpatialPolygonsDataFrame of length 1 ordered that represents the greater watershed tributary of the current zhyd watershed object.
c3	An object SpatialPolygonsDataFrame of length 1 ordered that represents the second watershed tributary of the current zhyd watershed object.
node_trib	An object SpatialPointsDataFrame of length 2 that represents the station points of the tributary watershed objects.

Methods

signature(x = "SpatialPointsDataFrame", xo = "SpatialPointsDataFrame", y = "SpatialLinesDataFrame", z = "SpatialPointsDataFrame", zhyd = "SpatialPolygonsDataFrame", c1 = "SpatialPolygonsDataFrame")

See Also

See Also the class Watershed and the methods Watershed.Order and Watershed.Order2.

Examples

library(Watersheds)
data(WatershedsData)

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

station1 = SpatialPoints(coords=cbind(4328448.74, 3118576.86), 
	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]]
c1_inlet = a[[2]]
c1_outlet = a[[3]]
sb1 = a[[7]]
riverIO = a[[8]]
nodeIO = a[[9]]			
c1_river = a[[10]]
c1_node = a[[11]]	

a = Watershed.Tributary(x=c1_inlet,xo= c1_outlet,y=riverIO,z=nodeIO,zhyd=zhyd1, c1=c1)
c2c3 = a[[1]]
c2 = a[[2]]
c3 = a[[3]]
node_trib = a[[4]]
	
bbox1 = slot(c2c3, "bbox")
bbox = matrix(0,2,2)
bbox[,1] = bbox1[,1]*.998
bbox[,2] = bbox1[,2]*1.002

plot(c1, xlim=bbox[1,], ylim=bbox[2,],col="gray50")			
plot(c2, col="gray75", add=TRUE)
plot(c3, col="gray85", add=TRUE)
plot(slot(watershed,"station"),pch=24, bg="blue",add= TRUE)
plot.PolyLineAttribute(c1, "order", 450, 0.8)
plot.PolyLineAttribute(c2, "order", 450, 0.8)
plot.PolyLineAttribute(c3, "order", 450, 0.8)				
plot(c1_river, col="blue", add=TRUE)
plot(node_trib,pch=21,bg="red",cex=.8,add=TRUE)
plot.PointAttribute(node_trib,"ELEV",600,0.7)
title(main="Current zhyd watershed (1)",
	sub="First order tributary nodes (1.1, 1.2)")	

Watersheds documentation built on May 2, 2019, 12:10 p.m.