It extracts the channel network from the drainage directions in five possible ways:
- by using a threshold value on the contributing areas, then only
the pixels with contributing area greater than the threshold are the
channel heads (mode 0)
- by using a threshold value of the estimated stress tangential to the bottom (mode 1)
- by using a threshold value on the estimated stress tangential to the bottom only in convergent sites (mode 2)
Individuated the beginning of the channel incision, the points upriver
are considered as canalized.
Inputs
- the file containing the matrix of the drainage directions (obtained with DrainDir or Markoutlets)
- the file containing the matrix of contributing areas (calculated
with DrainDir, Tca or Multitca)
- if we choose mode 1, the file containing the matrix of the slopes (calculated with Gradient or Slope)
- if we choose mode 2, the file containing the matrix with the classes aggregated by GC
Output
- the file containing the matrix which individuates the channel
network
- the file containing the total contributing areas calculated with the new drainage directions
- if we choose to create the vector layer of the extracted network, the vector map with the channel network ("filePath" is the path
plus the file name of the vector layer)