Identifies all cells belonging to a (sub)catchment and returns them in a vector
[C, BoundingBox [, INLETS] ] = GetCatchmentCells ( handle_fdr , outlet [, POT_INLETS] )
integer: handle on the flow direction grid defining the stream network
integer: ID of the outlet pixel, as returned by XY2rec
(optional) vector of potential inlet IDs. By default, the function returns all cells upstream the outlet pixel. If you want to exclude some subcatchments then provide a vector of inlet IDs defining those subcatchments. The recursive accumulation will not go further upstream whenever one of these inlets is met.
output: vector of (sub)catchement cells.
output, matrix [4 x 1]: bounding box of the (sub)catchment in the form [colmin,colmax,rowmin,rowmax]
output: vector of IDs of the inlet sites which have actually been encountered in the extraction
Identifies all cells belonging to a (sub)catchment and return them in a vector