Description Value Methods See Also Examples
S4 Method for function Watershed.Order
. Definition of the properties of the current zhyd
watershed.
The method returns a list of 11 objects:
c1 |
An object |
c1_inlet |
An object |
c1_outlet |
An object |
c2 |
An object |
c3 |
An object |
node_trib |
An object |
sb1 |
An object |
riverIO |
An object |
nodeIO |
An object |
c1_river |
An object |
c1_node |
An object |
signature(x = "Watershed")
The function takes the object of class Watershed
and identifies the subbasin that contains the current station
(class SpatialPoints
) and subsets the zhyd
object to subbasin and extract the current zhy
object that contains station
. Posteriorly, identifies the inlet and outlet stretches and probable inlet and outlet nodes of the zhyd
. Then, runs the functions Watershed. ,IOR1, IOR2, IOR3
, or IOR4
for determining the actual inlet and outlet nodes. Finally, the method executes the S4 method Watershed.Tributary
for defining tributary nodes and tributary catchments of the current zhyd
watershed.
See Also the class Watershed
and the methods Watershed.Order2
and Watershed.Tributary
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 | 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]]
c2 = a[[4]]
c3 = a[[5]]
node_trib = a[[6]]
sb1 = a[[7]]
riverIO = a[[8]]
nodeIO = a[[9]]
c1_river = a[[10]]
c1_node = a[[11]]
bbox1 = slot(c1, "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(c1_node,pch=21,bg="blue",cex=.5,add=TRUE)
plot(nodeIO,pch=21,bg="blue",cex=.5,add=TRUE)
plot(c1_inlet, pch=21, bg="green",add= TRUE)
plot(c1_outlet,pch=21, bg="red",add= TRUE)
plot.PointAttribute(nodeIO,"ELEV",600,0.7)
title(main="Current zhyd watershed (1)",
sub="First order tributary watersheds (1.1, 1.2)")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.