1r.carve(1) Grass User's Manual r.carve(1)
2
6 r.carve - Generates stream channels.
7 Takes vector stream data, transforms it to raster and subtracts depth
8 from the output DEM.
9
11 raster, hydrology
12
14 r.carve
15 r.carve --help
16 r.carve [-n] raster=name vector=name output=name [points=name]
17 [width=float] [depth=float] [--overwrite] [--help] [--verbose]
18 [--quiet] [--ui]
19 Flags:
21 -n
22 No flat areas allowed in flow direction
23 --overwrite
25 Allow output files to overwrite existing files
26 --help
28 Print usage summary
29 --verbose
31 Verbose module output
32 --quiet
34 Quiet module output
35 --ui
37 Force launching GUI dialog
38 Parameters:
40 raster=name [required]
41 Name of input raster elevation map
42 vector=name [required]
44 Name of input vector map containing stream(s)
45 Or data source for direct OGR access
46 output=name [required]
48 Name for output raster map
49 points=name
51 Name for output vector map for adjusted stream points
52 width=float
54 Stream width (in meters)
55 Default is raster cell width
56 depth=float
58 Additional stream depth (in meters)
59
61 r.carve accepts vector stream data as input, transforms them to raster,
62 and subtracts a default-depth + additional-depth from a DEM. If the
63 given width is more than 1 cell, it will carve the stream with the
64 given width. With the -n flag it should eliminate all flat cells within
65 the stream, so when and if the water gets into the stream it will flow.
66 The points option generates x,y,z for points which define the stream
67 with the z-value of the bottom of the carved-in stream. These points
68 can then be combined with contours to interpolate a new DEM with better
69 representation of valleys.
70
72 r.carve does not create a depressionless DEM because many depressions
73 are in flat areas and not in the streams.
74
76 North Carolina sample dataset:
77 # set computational region
78 g.region raster=elev_lid792_1m -p
79 # visualize original data
80 d.mon wx0
81 d.rast elev_lid792_1m
82 d.vect streams
83 # carve
84 r.carve rast=elev_lid792_1m vect=streams out=carved_dem width=3 depth=0.5
85 # visualize resulting carved DEM map
86 d.rast carved_dem
87 # visualize
88 r.relief input=elev_lid792_1m output=elev_lid792_1m_shaded
89 r.relief input=carved_dem output=carved_dem_shaded
90 d.rast elev_lid792_1m_shaded
91 d.erase
92 d.rast carved_dem_shaded
93 # flow accumulation
94 r.watershed elevation=elev_lid792_1m accumulation=elev_lid792_1m_accum
95 r.watershed elevation=carved_dem accumulation=carved_dem_accum
96 d.rast elev_lid792_1m_accum
97 d.erase
98 d.rast carved_dem_accum
99 Fig: Original 1m LiDAR based DEM with vector streams map on Fig: Original 1m LiDAR based DEM shown as shaded terrain
101 top
102 Fig: Carved 1m LiDAR based DEM Fig: Carved 1m LiDAR based DEM shown as shaded terrain
104 Fig: Flow accumulation in original 1m LiDAR based DEM Fig: Flow accumulation in carved 1m LiDAR based DEM
106
109 The module does not operate yet in latitude-longitude locations. It
110 has not been thoroughly tested, so not all options may work properly -
111 but this was the intention.
112
114 Terrain modeling and Soil Erosion Simulations for Fort Hood and Fort
115 Polk test areas, by Helena Mitasova, Lubos Mitas, William M. Brown,
116 Douglas M. Johnston, GMSL (Report for CERL 1999)
117
119 r.flow, r.fill.dir, r.watershed
120
122 Bill Brown (GMSL)
123 GRASS 6 update: Brad Douglas
124
126 Available at: r.carve source code (history)
127 Main index | Raster index | Topics index | Keywords index | Graphical
129 index | Full index
130 © 2003-2019 GRASS Development Team, GRASS GIS 7.8.2 Reference Manual
132GRASS 7.8.2 r.carve(1)