back to first page [..][[BR]]
back to ["Stream order"] [[BR]]

== Shreve Algortithm ==
#met plus d'une journée à tourner pour l'ensemble de la CCM (France).

{{{
data<-ccm@data
data<-data[order(data$wso_id,data$pfafstette),]
head(data)

Ndata<-dim(data)[1]
shreve<-rep(0,Ndata)
data$shreve<-shreve

for(j in as.numeric(levels(as.factor(data$wso_id)))){
  for(i in 1:length(data$wso_id[data$wso_id==j])){
    if(data$pfafstette[data$wso_id==j][i]<max(data$pfafstette[data$wso_id==j])){
      data$shreve[data$wso_id==j][i]<-sum(data$strahler[data$wso_id==j][i:length(data$wso_id[data$wso_id==j])]==1)
    }
    if(data$strahler[data$wso_id==j][i]==1){
      data$shreve[data$wso_id==j][i]<-1
    }
  }
}
data
}}}

'''Shreve order posant problème'''[[BR]]
- Les pfafstette maximum pour un wso_id donné devraient correspondre à un ordre de strahler de 1, or ce n'est pas le cas pour ces différents wso_id.
- problème avec les bassins versants ayant plusieurs sous bassins versants.