QGIS User Guide, Release 2.6
Outputs
D8 Contributing Area Grid [raster] A gridof contributing area values calculated as the cells own con-
tribution plus the contribution fromupslope neighbors that drain in to it according to the D8 flowmodel.
Consoleusage
processing.runalg(’taudem:d8contributingarea’-p, -o, -wg, -nc, -ad8)
See also
D8 Flow Directions
Description
Creates2 grids. The firstcontains theflowdirectionfromeachgridcelltoone ofitsadjacentordiagonalneighbors,
calculated using the direction of steepest descent. The second contain the slope, as evaluated in the direction of
steepest descent, and is reported as drop/distance, i.e. tan of the angle. Flow direction is reported as NODATA for
any grid cell adjacent to the edge of the DEM domain, or adjacent to a NODATA value in the DEM. In flat areas,
flow directions are assigned away from higher ground and towards lower ground using the method of Garbrecht
andMartz (1997). The D8 flow direction algorithm may be applied toa DEM that has not had its pits filled,but it
will then result in NODATA values for flow direction and slope at the lowest point of each pit.
D8 Flow Direction Coding:
• 1 — East
• 2 — Northeast
• 3 — North
• 4 — Northwest
• 5 — West
• 6 — Southwest
• 7 — South
• 8 — Southeast
The flow direction routingacross flat areas is performedaccordingto the methoddescribedbyGarbrecht,J.andL.
W.Martz,(1997),“The AssignmentofDrainageDirectionOverFlatSurfaces in RasterDigital Elevation Models”,
Journal of Hydrology,193: 204-213.
Parameters
Pit Filled Elevation Grid [raster] A grid of elevation values. This is usually the output of the “Pit
Remove” tool, in which case it is elevations with pits removed. Pits are low elevation areas in digital
elevation models (DEMs) that are completely surrounded by higher terrain. They are generally taken to
be artifacts of the digitation process that interfere with the processing of flow across DEMs. So they are
removedbyraisingtheirelevationtothe pointwherethey just drainoffthe domain. This step is not essential
ifyouhave reason to believe that the pits inyour DEMare real. If a fewpits actuallyexistandsoshould not
be removed, while at the same time others are believed to be artifacts that need to be removed, the actual
pits should have NODATA elevation values inserted at their lowest point. NODATA values serve to define
edges of the domain in the flow field, and elevations are only raised to where flow is off an edge, so an
internal NODATA value will stopa pit frombeing removed, if necessary.
18.8. TauDEM algorithm provider
595
Pdf files thumbnail preview - Draw thumbnail images for PDF in C#.net, ASP.NET, MVC, Ajax, WinForms, WPF
Support Thumbnail Generation with Various Options for Quick PDF Navigation
pdf thumbnail; .pdf printing in thumbnail size
Pdf files thumbnail preview - VB.NET PDF Thumbnail Create SDK: Draw thumbnail images for PDF in vb.net, ASP.NET, MVC, Ajax, WinForms, WPF
Support Thumbnail Generation with Various Options for Quick PDF Navigation
pdf thumbnails in; pdf thumbnails
QGIS User Guide, Release 2.6
Outputs
D8 Flow Direction Grid [raster] A grid of D8 flow directions which are defined, for each cell, as the
direction of the one of its eight adjacent ordiagonal neighbors with the steepest downward slope.
D8 Slope Grid [raster] A grid giving slope in the D8 flow direction. This is measured as drop/distance.
Consoleusage
processing.runalg(’taudem:d8flowdirections’-fel, -p, -sd8)
See also
D-Infinity Contributing Area
Description
Calculates a grid of specific catchment area which is the contributing area per unit contour length using the
multiple flow direction D-infinity approach. D-infinity flow direction is defined as steepest downward slope on
planar triangular facets on a block centered grid. The contribution at each grid cell is taken as the grid cell length
(or when the optional weight grid input is used, from the weight grid). The contributing area of each grid cell is
then taken as its own contribution plus the contribution from upslope neighbors that have some fraction draining
to it according to the D-infinity flow model. The flow fromeach cell either all drains to one neighbor, ifthe angle
falls along acardinal(0,㔋/2,㔋,3㔋/2)orordinal (㔋/4, 3㔋/4,5㔋/4,7㔋/4)direction,or is onan angle fallingbetween
the direct angle to two adjacent neighbors. In the latter case the flow is proportioned between these two neighbor
cells according to how close the flow direction angle is to the direct angle to those cells. The contour length used
here is the grid cell size. The resulting units of the specific catchment area are length units the same as those of
the grid cell size.
When the optional weight grid is not used, the result is reported in terms of specific catchment area, the upslope
area per unit contour length, taken here as the number of cells times grid cell length (cell area divided by cell
length). This assumes that grid cell length is the effective contour length, in the definition of specific catchment
area and does not distinguish any difference in contour length dependent upon the flow direction. When the
optional weight grid is used,the result is reported directly as a summation of weights, without any scaling.
If the optional outlet point shapefile is used, only the outlet cells and the cells upslope (by the D-infinity flow
model) ofthemare in the domain to be evaluated.
By default, the tool checks foredge contamination. This is defined as the possibility that a contributing area value
may be underestimated due to grid cells outside of the domain not being counted. This occurs when drainage
is inwards from the boundaries or areas with “no data” values for elevation. The algorithm recognizes this and
reports “no data” for the contributing area. It is common to see streaks of “no data” values extending inwards
from boundaries alongflow paths that enterthe domain at a boundary. This is the desired effect and indicates that
contributing area forthese grid cells is unknown due to it being dependent on terrain outside ofthe domainof data
available. Edge contamination checking may be turned off in cases where you know it is not an issue or want to
ignore these problems, iffor example, the DEM has beenclipped along a watershed outline.
Parameters
D-Infinity Flow Direction Grid [raster] A grid of flow directions based on the D-infinity flow
method using the steepest slope of a triangular facet. Flow direction is determined as the direction of
the steepest downward slope on the 8 triangular facets of a 3x3 block centered grid. Flow direction is en-
coded as an angle in radians, counter-clockwise fromeast as a continuous (floating point) quantitybetween
0and 2㔋. The resulting flow in a grid is then usually interpreted as being proportioned between the two
neighboring cells that define the triangular facet with the steepest downward slope.
596
Chapter 18. Processing providers and algorithms
How to C#: Preview Document Content Using XDoc.Word
With the SDK, you can preview the document content according to the preview thumbnail by the ways as following. C# DLLs for Word File Preview. Add references:
no pdf thumbnails in; disable pdf thumbnails
How to C#: Preview Document Content Using XDoc.PowerPoint
With the SDK, you can preview the document content according to the preview thumbnail by the ways as following. C# DLLs: Preview PowerPoint Document.
view pdf thumbnails; print pdf thumbnails
QGIS User Guide, Release 2.6
Outlets Shapefile [vector: point] Optional.
Apoint shapefile defining the outlets of interest. If this input file is used, only the cells upslope of these
outlet cells are considered to be within the domain being evaluated.
Weight Grid [raster] Optional.
Agrid giving contribution to flow for each cell. These contributions (also sometimes referred to as weights
orloadings)areusedin the contributingarea accumulation. If this input fileis not used,theresult is reported
in terms of specific catchment area (the upslope area per unit contour length) taken as the number of cells
times grid cell length (cell area divided by cell length).
Check for edge contamination [boolean] A flag that indicates whether the tool should check for edge
contamination. Edge contamination is defined as the possibility that a contributing area value may be
underestimated due to the fact that grid cells outside of the domain have not been evaluated. This occurs
when drainage is inwards from the boundaries or areas with NODATA values for elevation. The algorithm
recognizes this and reports NODATAfor the impated cells. It is common to see streaks of NODATA values
extendinginwards fromboundaries along flowpaths that enter the domain ata boundary. This is the desired
effect and indicates that contributingarea forthesegridcells is unknowndue toitbeing dependentonterrain
outside of the domain ofavailable data. Edge contaminationcheckingmay be turned off incases where you
knowthis is not anissue,or wantto ignore these problems,iffor example,the DEM has been clipped along
awatershedoutline.
Default: True
Outputs
D-Infinity Specific Catchment Area Grid [raster] Agridofspecific catchmentarea whichisthe
contributing area per unit contour length using the multiple flow direction D-infinity approach. The con-
tributing area of each grid cell is then taken as its own contribution plus the contribution from upslope
neighbors that have some fraction draining to it according to the D-infinity flow model.
Consoleusage
processing.runalg(’taudem:dinfinitycontributingarea’-ang, -o, -wg, -nc, -sca)
See also
D-Infinity Flow Directions
Description
Assigns a flow direction based on the D-infinity flow method using the steepest slope of a triangular facet (Tar-
boton, 1997, “A New Method for the Determination of Flow Directions and Contributing Areas in Grid Digital
Elevation Models”, Water Resources Research, 33(2): 309-319). Flow direction is defined as steepest down-
ward slope on planar triangular facets on a block centered grid. Flow direction is encoded as an angle in radians
counter-clockwise from east as a continuous (floating point)quantity between 0 and 2㔋. The flow direction angle
is determined as the direction of the steepest downward slope on the eight triangular facets formed in a 3 x 3 grid
cell window centered on the grid cell of interest. The resulting flow in a grid is then usually interpreted as being
proportioned betweenthe two neighboring cells that define the triangular facet with the steepest downward slope.
Ablock-centered representation is used with each elevation value taken to represent the elevation of the center
of the corresponding grid cell. Eight planar triangular facets are formed between each grid cell and its eight
neighbors. Each of these has a downslope vector whichwhen drawn outwards from the center may be at an angle
that lies within or outside the 45 degree (㔋/4 radian) angle range ofthefacet at the center point. Ifthe slope vector
18.8. TauDEM algorithm provider
597
How to C#: Preview Document Content Using XDoc.excel
document in memory. With the SDK, you can preview the document content according to the preview thumbnail by the ways as following.
pdf file thumbnail preview; how to create a thumbnail of a pdf document
How to C#: Overview of Using XDoc.PowerPoint
Tell C# users how to: create a new PowerPoint file and load PowerPoint; merge, append, and split PowerPoint files; insert, delete, move, rotate Create Thumbnail.
pdf thumbnail generator; how to show pdf thumbnails in
QGIS User Guide, Release 2.6
angle is within the facet angle, it represents the steepest flow direction on that facet. If the slope vector angle is
outside a facet, the steepest flow direction associated with that facet is taken along the steepest edge. The slope
and flow direction associated with the grid cell is taken as the magnitude and direction of the steepest downslope
vectorfromall eight facets. Slope is measured as drop/distance, i.e. tan of the slope angle.
In the case where no slope vectors are positive (downslope),theflowdirectionis set usingthe method of Garbrecht
andMartz(1997)forthedetermination offlowacross flat areas. This makes flatareas drainawayfromhighground
andtowards low ground. The flow path grid toenforce drainage along existing streams is an optionalinput,and if
used, takes precedence over elevations forthe setting of flow directions.
The D-infinity flow direction algorithm may be applied to a DEM that has not had its pits filled, but it will then
result in“no data” values for the D-infinity flow direction and slope associated with the lowest point ofthe pit.
Parameters
Pit Filled Elevation Grid [raster] A grid of elevation values. This is usually the output of the “Pit
Remove” tool, in which case it is elevations with pits removed.
Outputs
D-Infinity Flow Directions Grid [raster] A grid of flow directions based on the D-infinity flow
method using the steepest slope of a triangular facet. Flow direction is determined as the direction of
the steepest downward slope on the 8 triangular facets of a 3x3 block centered grid. Flow direction is en-
coded as an angle in radians, counter-clockwise fromeast as a continuous (floating point) quantitybetween
0and 2㔋. The resulting flow in a grid is then usually interpreted as being proportioned between the two
neighboring cells that define the triangular facet with the steepest downward slope.
D-Infinity Slope Grid [raster] A grid of slope evaluated using the D-infinity method described in Tar-
boton, D. G., (1997), “A New Method for the Determination of Flow Directions and Contributing Areas in
Grid Digital Elevation Models”, Water Resources Research, 33(2): 309-319. This is the steepest outwards
slope on one of eight triangular facets centered at each grid cell, measured as drop/distance, i.e. tan of the
slope angle.
Consoleusage
processing.runalg(’taudem:dinfinityflowdirections’-fel, -ang, -slp)
See also
Grid Network
Description
Creates 3 grids that contain for each grid cell: 1) the longest path, 2) the total path, and 3) the Strahler order
number. These values are derived from the network defined by the D8 flow model.
The longest upslope length is the length of the flow path from the furthest cell that drains to each cell. The total
upslope path length is the length of the entire grid network upslope of each grid cell. Lengths are measured
between cell centers taking into account cell size and whetherthe direction is adjacent or diagonal.
Strahler order is defined as follows: A network of flow paths is defined by the D8 Flow Direction grid. Source
flow paths have a Strahler order number of one. When two flow paths of different order join the order of the
downstream flow path is the order of the highest incoming flow path. When two flow paths of equal order join
the downstreamflow path order is increased by 1. When more than two flowpaths join the downstreamflow path
order is calculated as the maximum of the highest incoming flow path order or the second highest incoming flow
path order +1. This generalizes the common definition to cases where more than two flow paths join at a point.
598
Chapter 18. Processing providers and algorithms
How to C#: Overview of Using XDoc.Word
Tell C# users how to: create a new Word file and load Word from pdf; merge, append, and split Word files; insert, delete, move, rotate Create Thumbnail.
pdf files thumbnails; cannot view pdf thumbnails in
How to C#: Set Image Thumbnail in C#.NET
Convert Jpeg to PDF; Merge PDF Files; Split PDF Document; Remove Password from VB.NET How-to, VB.NET PDF, VB.NET Word, VB How to C#: Set Image Thumbnail in C#.NET
create thumbnail jpeg from pdf; pdf thumbnail viewer
QGIS User Guide, Release 2.6
Where the optional mask grid and threshold value are input, the function is evaluated only considering grid cells
that lie in the domain with mask grid value greater than or equal to the threshold value. Source (first order) grid
cells are taken as those that do not have any other grid cells from inside the domain draining in to them, and
only when two of these flow paths join is order propagated according to the ordering rules. Lengths are also only
evaluated counting paths within the domain greater than or equal to the threshold.
If the optional outlet point shapefile is used, only the outlet cells and the cells upslope (by the D8 flow model) of
them are in the domain tobe evaluated.
Parameters
D8 Flow Direction Grid [raster] A grid of D8 flow directions which are defined, for each cell, as the
direction of the one of its eight adjacent ordiagonal neighbors with the steepest downward slope. This grid
can be obtained as the output ofthe “D8 FlowDirections”tool.
Outlets Shapefile [vector: point] Optional.
Apoint shapefile defining the outlets of interest. If this input file is used, only the cells upslope of these
outlet cells are considered to be within the domain being evaluated.
Mask Grid [raster] Optional.
Agrid that is used to determine the domain do be analyzed. If the mask grid value >= mask threshold (see
below), then the cell will be included in the domain. While this tool does not have an edge contamination
flag, if edge contamination analysis is needed, then a mask grid from a function like “D8 Contributing
Area” that does support edge contamination can be used to achieve the same result.
Mask Threshold [number] This input parameteris used in thecalculationmaskgridvalue>=maskthreshold
to determine if the grid cell is in the domain to be analyzed.
Default: 100
Outputs
Longest Upslope Length Grid [raster] A grid that gives the length ofthelongest upslope D8flow path
terminating at each grid cell. Lengths are measured between cell centers taking into account cell size and
whetherthe direction is adjacent or diagonal.
Total Upslope Length Grid [raster] The total upslope path length is the length of the entire D8 flow
grid network upslope of each grid cell. Lengths are measured between cell centers taking into account cell
size and whether the direction is adjacent or diagonal.
Strahler Network Order Grid [raster] A grid giving the Strahler order number for each cell. A net-
workofflowpaths is definedbytheD8FlowDirectiongrid. Source flow paths have aStrahlerordernumber
of one. When two flow paths of different order join the order of the downstream flow path is the order of
the highest incoming flow path. When two flow paths of equal order join the downstream flow path order
is increased by 1. When more than two flow paths join the downstream flow path order is calculated as the
maximumofthe highest incomingflow path order or the secondhighest incoming flowpathorder+ 1. This
generalizes the common definition to cases where more thantwo flow paths join at a point.
Consoleusage
processing.runalg(’taudem:gridnetwork’, d8_flow_dir_grid, outlets_shape, mask_grid, threshold, longest_len_grid, total_len_grid, strahler_grid)
18.8. TauDEM algorithm provider
599
How to C#: Generate Thumbnail for Excel
Document. Conversion. Convert Excel to PDF. Convert Excel File: Merge Excel Files. File: Split Excel Document. Insert Image. Thumbnail Create. Thumbnail Create. |
create thumbnail jpg from pdf; pdf files thumbnail preview
How to C#: Generate Thumbnail for Word
Images. Convert Word to ODT. Convert PDF to Word. a Word File. File: Merge Word Files. File: Split Text Search. Insert Image. Thumbnail Create. Thumbnail Create.
show pdf thumbnails; thumbnail view in for pdf files
QGIS User Guide, Release 2.6
See also
Pit Remove
Description
Identifies all pits in the DEM and raises theirelevationto the level of the lowest pourpoint around theiredge. Pits
are low elevation areas indigitalelevationmodels (DEMs)thatare completely surroundedby higher terrain. They
are generally taken to be artifacts that interfere with the routing of flow across DEMs, so are removed by raising
their elevation to the point where they drain off the edge of the domain. The pour point is the lowest point on the
boundary of the “watershed” draining to the pit. This step is not essential if you have reason to believe that the
pits in yourDEM are real. If a few pits actuallyexistand so shouldnot be removed, while at the same time others
are believed to be artifacts thatneedto be removed, the actual pits shouldhave NODATA elevation values inserted
at theirlowestpoint. NODATA values serve to define edges in the domain, and elevations are only raised towhere
flow is offan edge, so an internal NODATA value will stop a pit from being removed,if necessary.
Parameters
Elevation Grid [raster] A digital elevation model (DEM) grid to serve as the base input for the terrain
analysis and stream delineation.
Outputs
Pit Removed Elevation Grid [raster] Agridofelevationvalues withpits removedsothatflowis routed
off of the domain.
Consoleusage
processing.runalg(’taudem:pitremove’-z, -fel)
See also
.
18.8.2 Specialized Grid Analysis
D8 Distance To Streams
Description
Computes the horizontal distance to stream foreach grid cell, moving downslope according tothe D8 flowmodel,
until a stream grid cell is encountered.
Parameters
D8 Flow Direction Grid [raster] This input is a grid of flow directions that are encoded using the D8
method where all flow from a cells goes to a single neighboring cell in the direction of steepest descent.
This grid can be obtained as the output of the “D8 Flow Directions” tool.
600
Chapter 18. Processing providers and algorithms
How to C#: Overview of Using XDoc.Excel
Tell C# users how to: create a new Excel file and load Excel; merge, append, and split Excel files; insert, delete, move, rotate, copy and Create Thumbnail.
pdf thumbnails in; create thumbnail jpeg from pdf
How to C#: Generate Thumbnail for PowerPoint
Conversion. Convert PowerPoint to PDF. Convert PowerPoint to File. File: Merge PowerPoint Files. File: Split PowerPoint. Text Search. Insert Image. Thumbnail Create
can't see pdf thumbnails; pdf thumbnail preview
QGIS User Guide, Release 2.6
Stream Raster Grid [raster] A grid indicating streams. Such a grid can be created by several of the tools
in the “Stream NetworkAnalysis” toolset. However, the tools in the “Stream NetworkAnalysis” toolset
only create grids with a value of 0 for no stream, or 1 for stream cells. This tool can also accept grids
with values greater than 1, which can be used in conjunction with the Threshold parameter to determine
the location of streams. This allows Contributing Area grids to be used to define streams as well as the
normal StreamRaster grids. This grid expects integer (long integer) values and any non-integer values will
be truncated toan integerbefore being evaluated.
Threshold [number] This value acts as threshold on the Stream Raster Grid to determine the location
of streams. Cells with a Stream Raster Grid value greater than or equal to the Threshold value
are interpreted as streams.
Default: 50
Outputs
Output Distance to Streams [raster] A grid giving the horizontal distance along the flow path as de-
fined by the D8 Flow Directions Grid to the streams in the Stream Raster Grid.
Consoleusage
processing.runalg(’taudem:d8distancetostreams’-p, -src, -thresh, -dist)
See also
D-Infinity Avalanche Runout
Description
Identifies an avalanche’s affected area and the flow path length to each cell in that affacted area. All cells downs-
lope from each source area cell, up to the point where the slope from the source to the affected area is less than
athreshold angle called the Alpha Angle can be in the affected area. This tool uses the D-infinity multiple flow
directionmethod for determining flowdirection. This will likelycause verysmall amounts offlow to be dispersed
to some downslope cells that might overstate the affected area, so a threshold proportion can be set to avoid this
excess dispersion. The flow path length is the distance from the cell in question to the source cell that has the
highest angle.
All points downslope from the source area are potentially in the affected area, but not beyond a point where the
slope from the source to the affected area is less than a threshold angle called the Alpha Angle.
Slope is to be measured using the straight line distance from source point to evaluation point.
It makes morephysicalsensetomefortheangle tobe measuredalongtheflowpath. Nevertheless itisequallyeasy
tocode straight line angles as angles along the flow path,so an option that allows switching will be provided. The
most practical way toevaluate avalanche runout is to keep track of the source pointwith the greatest angle toeach
point. Then the recursive upslope flow algebra approach will look at a grid cell and all its upslope neighbors that
flow to it. Information from the upslope neighbors will be used to calculate the angle to the grid cell in question
andretain it in the runout zone ifthe angle exceeds the alpha angle. This procedure makes the assumptionthatthe
maximumangle at a grid cell will be from the set of cells that have maximum angles to the inflowing neighbors.
This willalways be true ofangle is calculatedalong a flow path,but Ican conceive ofcases where flow paths bend
back on themselves where this would not be the case for straight line angles.
TheD-infinity multiple flow direction fieldassigns flow fromeachgridcellto multiple downslope neighbors using
proportions (Pik) that vary between 0 and 1 and sum to 1 for all flows out of a grid cell. It may be desirable to
specify a threshold T that this proportion has to exceed before a grid cell is counted as flowing to a downslope
18.8. TauDEM algorithm provider
601
QGIS User Guide, Release 2.6
grid cell, e.g. Pik > T (=0.2 say) to avoid dispersion to grid cells that get very little flow. T will be specified as
auser input. If all upslope grid cells are to be used T may be input as 0.
Avalanche source sites are to be input as a short integer grid (name suffix
*
ass, e.g. demass) comprised of
positive values where avalanches may be triggered and 0 values elsewhere.
The following grids are output:
• rz — A runout zone indicator with value 0 to indicate that this grid cell is not in the runout zone and value
>0 to indicate that this grid cell is in the runout zone. Since there may be information in the angle to the
associated source site,this variable will be assigned the angle to the source site (in degrees)
• dm —Along flow distance fromthe source site that has the highest angle to the point in question
Parameters
D-Infinity Flow Direction Grid [raster] A grid giving flow direction by the D-infinity method.
Flow direction is measured in radians, counter clockwise from east. This can be created by the tool “D-
Infinity Flow Directions”.
Pit Filled Elevation Grid [raster] This input is a grid of elevation values. As a general rule, it is
recommended that you use a grid of elevation values that have had the pits removed for this input. Pits are
generally taken to be artifacts that interfere with the analysis offlow across them. This grid can be obtained
as the output of the “Pit Remove” tool, in which case it contains elevation values where the pits have been
filled tothe point where they just drain.
Avalanche Source Site Grid [raster] This is a grid of source areas for snow avalanches that are com-
monly identified manually using a mix of experience and visual interpretation of maps. Avalanche source
sites are to be input as a short integer grid (name suffix
*
ass, e.g. demass) comprised of positive values
where avalanches may be triggered and 0 values elsewhere.
Proportion Threshold [number] This value is a threshold proportion that is used to limit the disperson
of flow caused by using the D-infinity multiple flow direction method for determining flow direction. The
D-infinity multiple flow direction method often causes very small amounts of flow to be dispersed to some
downslope cells that might overstate the affected area, so a threshold proportion can be set to avoid this
excess dispersion.
Default: 0.2
Alpha Angle Threshold [number] This value is the threshold angle, called the Alpha Angle, that is used
to determine which of the cells downslope from the source cells are in the affected area. Only the cells
downslope from each source area cell, up to the point where the slope from the source to the affected area
is less than a threshold angle are in the affected area.
Default: 18
Measure distance along flow path [boolean] This option selects the method used to measure the
distance used to calculate the slope angle. If option is True then measure it along the flow path, where
the False option causes the slope to be measure along the straight line distance from the source cell to the
evaluation cell.
Default: True
Outputs
Runout Zone Grid [raster] This grid Identifies the avalanche’s runout zone (affected area) using a runout
zone indicator with value 0 to indicate that this grid cell is not in the runout zone and value > 0 to indicate
that this grid cell is in the runout zone. Since there may be information in the angle to the associated source
site,this variable will be assigned the angle to the source site (in degrees).
Path Distance Grid [raster] This is a grid of the flow distance from the source site that has the highest
angle to each cell.
602
Chapter 18. Processing providers and algorithms
QGIS User Guide, Release 2.6
Consoleusage
processing.runalg(’taudem:dinfinityavalancherunout’-ang, -fel, -ass, -thresh, -alpha, -direct, -rz, -dfs)
See also
D-Infinity Concentration Limited Accumulation
Description
This function applies to the situation where an unlimited supply ofa substance is loaded into flow at a concentra-
tion orsolubility threshold Csol over a region indicated by an indicator grid (dg). It a grid ofthe concentration of
asubstance at each location in the domain, where the supply of substance from a supply area is loaded into the
flow at a concentration or solubility threshold. The flow is first calculated as a D-infinity weighted contributing
area of an input Effective Runoff Weight Grid (notionally excess precipitation). The concentation of substance
overthe supply area (indicatorgrid) is at the concentrationthreshold. As the substance moves downslope withthe
D-infinity flow field,it is subject to firstorder decay in moving fromcell to cell as well as dilution due to changes
in flow. The decay multiplier grid gives the fractional (first order) reduction in quantity in moving fromgrid cell
xto the next downslope cell. If the outlets shapefile is used, the tool only evaluates the part of the domain that
contributes flow to the locations given by the shapefile. This is useful for a tracking a contaminant or compound
from an area with unlimited supply of that compound that is loaded into a flow at a concentration or solubility
threshold overa zone and flow fromthe zone may be subject to decay orattenuation.
The indicator grid (dg) is used to delineate the area of the substance supply using the (0, 1) indicator function
i(x). A[] denotes the weighted accumulation operator evaluated using the D-Infinity Contributing Area func-
tion. The Effective Runoff Weight Grid gives the supply to the flow (e.g. the excess rainfall if this is overland
flow) denoted as w(x). The specific discharge is then given by:
Q(x)=A[w(x)]
This weighted accumulation Q(x) is output as the Overland Flow Specific Discharge Grid. Over the substance
supply area concentration is at the threshold (the threshold is a saturation orsolubility limit). If i(x) = 1, then
C(x) = Csol, and L(x) = Csol Q(x),
where L(x) denotes the load being carried by the flow. At remaining locations, the load is determined by load
accumulation and the concentration by dilution:
Here d(x) = d(i, j) is a decay multiplier giving the fractional (first order) reduction in mass in moving
from grid cell x to the next downslope cell. If travel (or residence) times t(x) associated with flow between
cells are available d(x) may be evaluated as exp(-k t(x)) where k is a first order decay parameter. The
Concentration grid output is C(x). If the outlets shapefile is used, the tool only evaluates the part of the domain
that contributes flow to the locations given by the shapefile.
Useful for a tracking a contaminant released or partitioned to flow at a fixed threshold concentration.
Parameters
D-Infinity Flow Direction Grid [raster] A grid giving flow direction by the D-infinity method.
Flowdirectionis measured inradians,counter clockwise from east. This grid canbe created bythefunction
“D-Infinity FlowDirections”.
Disturbance Indicator Grid [raster] A grid that indicates the source zone of the area of substance
supply and must be 1 inside the zone and 0 orNODATAover the rest ofthe domain.
18.8. TauDEM algorithm provider
603
QGIS User Guide, Release 2.6
Decay Multiplier Grid [raster] Agridgivingthefactorby which flow leavingeachgridcell is multiplied
before accumulation on downslope grid cells. This maybe usedto simulate the movement ofan attenuating
or decaying substance. Iftravel (orresidence)times t(x) associated with flow between cells are available
d(x) may be evaluated as exp(-k t(x)) where k is a first order decayparameter.
Effective Runoff Weight Grid [raster] Agrid givingtheinput quantity(notionallyeffective runoff or
excess precipitation) to be used in the D-infinity weighted contributing area evaluation of Overland Flow
Specific Discharge.
Outlets shapefile [vector: point] Optional.
This optional input is a point shapefile defining outlets of interest. If this file is used, the tool will only
evaluate the area upslope of these outlets.
Concentration Threshold [number] The concentration or solubility threshold. Over the substance sup-
ply area,concentration is at this threshold.
Default: 1.0
Check for edge contamination [boolean] This option determines whether the tool should check for
edge contamination. Edge contamination is defined as the possibility that a value may be underestimated
due to grid cells outside of the domain not being considered when determining contributing area.
Default: True
Outputs
Concentration Grid [raster] A grid giving the resulting concentration of the compound of interest in the
flow.
Consoleusage
processing.runalg(’taudem:dinfinityconcentrationlimitedaccumulation’-ang, -dg, -dm, -q, -o, -csol, -nc, -ctpt)
See also
D-Infinity Decaying Accumulation
Description
The D-Infinity Decaying Accumulation tool creates a grid of the accumulated quantity at each location in the
domain where the quantity accumulates with the D-infinity flow field,but is subject to first order decay in moving
from cell to cell. By default, the quantity contribution of each grid cell is the cell length to give a per unit width
accumulation, but can optionally be expressed with a weight grid. The decay multiplier grid gives the fractional
(first order) reduction in quantity in accumulating from grid cell x to the next downslope cell.
Adecayed accumulationoperatorDA[.] takes as inputa mass loading fieldm(x) expressedat each grid location
as m(i, j) that is assumed to move with the flow field but is subject to first order decay in moving from cell to
cell. The output is the accumulated mass at each location DA(x). The accumulation of m at each grid cell can be
numerically evaluated.
Here d(x) = d(i ,j) is adecaymultipliergiving the fractional(firstorder)reductioninmass in moving from
grid cell x to the next downslope cell. If travel (or residence) times t(x) associated with flow between cells are
available d(x) may be evaluated as exp(-k t(x)) where k is a first order decay parameter. The weight grid
is used to represent the mass loading m(x). If not specified this is taken as 1. If the outlets shapefile is used the
function is only evaluated on that part ofthe domain that contributes flow to the locations given by the shapefile.
604
Chapter 18. Processing providers and algorithms
Documents you may be interested
Documents you may be interested