back to first page ..
back to Stream order

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

• 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.