← Back to branch summary

~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to raster/r.fill.dir/description.html

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
raster/r.fill.dir/description.html
1
 
<h2>DESCRIPTION</h2>
2
 
 
3
 
 
4
 
<em>r.fill.dir</em> filters and generates a depressionless
5
 
elevation map and a flow direction map from a given raster elevation map.
6
 
 
7
 
<h2>NOTES</h2>
8
 
 
9
 
The <b>type</b> parameter is the type of format at which the user wishes to create
10
 
the flow direction map. The <i>agnps</i> format gives category values from
11
 
1-8, with 1 facing north and increasing values in the clockwise direction.
12
 
The <i>answers</i> format gives category values from 0-360 degrees, with 0
13
 
(360) facing east and values increasing in the counter clockwise direction
14
 
at 45 degree increments. The <i>grass</i> format gives the same category
15
 
values as the <em><a href="r.slope.aspect.html">r.slope.aspect</a></em>
16
 
program.<p>
17
 
 
18
 
The method adopted to filter the elevation map and rectify it is
19
 
based on the paper titled "Software Tools to Extract Structure from Digital
20
 
Elevation Data for Geographic Information System Analysis" by S.K. Jenson
21
 
and J.O. Domingue (1988).
22
 
<p>
23
 
 
24
 
The procedure takes an elevation layer as input and initially fills all the
25
 
depressions with one pass across the layer. Next, the flow direction
26
 
algorithm tries to find a unique direction for each cell. If the watershed
27
 
program detects areas with pothholes, it delineates this area from the rest
28
 
of the area and once again the depressions are filled using the neighborhood
29
 
technique used by the flow direction routine. The final output will be a
30
 
depressionless elevation layer and a unique flow direction layer.
31
 
<p>
32
 
 
33
 
This (D8) flow algorithm performs as follows: At each raster cell the code
34
 
determines the slope to each of the 8 surrounding cells and assigns the flow
35
 
direction to the highest slope out of the cell.  If there is more than one
36
 
equal, non-zero slope then the code picks one direction based on preferences
37
 
that are hard-coded into the program.  If the highest slope is flat and in
38
 
more than one direction then the code first tries to select an alternative
39
 
based on flow directions in the adjacent cells. <em>r.fill.dir</em> iteratates that process,
40
 
effectively propagating flow directions from areas where the directions are
41
 
known into the area where the flow direction can't otherwise be resolved.
42
 
 
43
 
<p>
44
 
The flow direction map can be encoded in either ANSWERS (Beasley et.al,
45
 
1982) or AGNPS (Young et.al, 1985) form, so that it can be readily used as
46
 
input to these hydrologic models. The resulting depressionless elevation
47
 
layer can further be manipulated for deriving slopes and other attributes
48
 
required by the hydrologic models.
49
 
<p>
50
 
 
51
 
In case of local problems, those unfilled areas can be stored optionally.
52
 
Each unfilled area in this maps is numbered. The <b>-f</b> flag
53
 
instructs the program to fill single-cell pits but otherwise to just find
54
 
the undrained areas and exit. With the <b>-f</b> flag set the program
55
 
writes an elevation map with just single-cell pits filled, a direction map
56
 
with unresolved problems and a map of the undrained areas that were found
57
 
but not filled. This option was included because filling DEMs was often not
58
 
the best way to solve a drainage problem. These options let the user get a
59
 
partially-fixed elevation map, identify the remaining problems and fix the
60
 
problems appropriately.
61
 
<p>
62
 
 
63
 
<em>r.fill.dir</em> is sensitive to the current window setting. Thus 
64
 
the program can be used to generate a flow direction map for any 
65
 
sub-area within the full map layer. Also, <em>r.fill.dir</em> is
66
 
sensitive to any <em>mask</em> in effect.
67
 
<p>
68
 
 
69
 
In some cases it may be necessary to run r.fill.dir repeatedly (using output
70
 
from one run as input to the next run) before all of problem areas are
71
 
filled.
72
 
 
73
 
<h2>EXAMPLE</h2>
74
 
 
75
 
<div class="code"><pre>
76
 
r.fill.dir input=ansi.elev elevation=ansi.fill.elev direction=ansi.asp
77
 
</pre></div>
78
 
 
79
 
<p>
80
 
will create a depressionless (sinkless) elevation map ansi.fill.elev and a flow
81
 
direction map ansi.asp for the type "grass".
82
 
 
83
 
<h2>SEE ALSO</h2>
84
 
 
85
 
<em>
86
 
<a href="r.fillnulls.html">r.fillnulls</a>,
87
 
<a href="r.slope.aspect.html">r.slope.aspect</a></em>
88
 
 
89
 
<p>
90
 
Beasley, D.B. and L.F. Huggins. 1982. ANSWERS (areal nonpoint source watershed environmental response simulation): User's manual. U.S. EPA-905/9-82-001, Chicago, IL, 54 p.
91
 
<p>
92
 
Jenson, S.K., and J.O. Domingue. 1988. Extracting topographic structure from
93
 
digital elevation model data for geographic information system analysis. Photogram. Engr. and Remote Sens. 54: 1593-1600.
94
 
<p>
95
 
Young, R.A., C.A. Onstad, D.D. Bosch and W.P. Anderson. 1985. Agricultural nonpoint surface pollution models (AGNPS) I and II model documentation. St. Paul: Minn. Pollution control Agency and Washington D.C., USDA-Agricultural Research
96
 
Service.
97
 
<p>
98
 
 
99
 
<h2>AUTHOR</h2>
100
 
Fortran version:
101
 
Raghavan Srinivasan, Agricultural Engineering Department, Purdue
102
 
University<br>
103
 
Rewrite to C with enhancements:
104
 
Roger S. Miller
105
 
 
106
 
<p><i>Last changed: $Date: 2008-05-15 20:59:22 +0200 (Thu, 15 May 2008) $</i>