r.to.vect extracts vectors (aka, "arcs") from a raster map. These arcs may represent linear features (like roads or streams), or may represent area edge features (like political boundaries, or soil mapping units).
r.thin and r.to.vect may create excessive nodes at every junction, and may create small spurs or "dangling lines" during the thinning and vectorization process. These excessive nodes and spurs may be removed using v.clean.
A true vector tracing of the area edges might appear blocky, since the vectors outline the edges of raster data that are stored in rectangular cells. To produce a better-looking vector map, r.to.vect smoothes the corners of the vector data as they are being extracted. At each change in direction (i.e., each corner), the two midpoints of the corner cell (half the cell's height and width) are taken, and the line segment connecting them is used to outline this corner in the resultant vector map. (The cell's cornermost node is ignored.) Because vectors are smoothed by this program, the resulting vector map will not be "true" to the raster map from which it was created. The user should check the resolution of the geographic region (and the original data) to estimate the possible error introduced by smoothing.
r.to.vect extracts only area edges from the named raster input file. If the raster map contains other data (i.e., line edges, or point data) the output may be wrong.
Conversion of raster points to vector points:
Random sampling of points:
g.region raster=elevation -p # random sampling of points (note that r.random also writes vector points) r.random elevation raster_output=elevrand1000 n=1000 r.to.vect input=elevrand1000 output=elevrand1000 type=point # univariate statistics of sample points v.univar elevrand1000 column=value type=point # compare to univariate statistics on original full raster map r.univar elevation
Conversion of raster lines to vector lines:
Vectorization of streams in watershed basins map:
g.region raster=elevation -p r.watershed elev=elevation stream=elev.streams thresh=50000 r.to.vect -s input=elev.streams output=elev_streams type=line # drop "label" column which is superfluous in this example v.db.dropcolumn map=elev_streams column=label v.db.renamecolumn map=elev_streams column=value,basin_id # report length per basin ID v.report map=elev_streams option=length units=meters sort=asc
Conversion of raster polygons to vector polygons:
Vectorization of simplified landuse class map:
g.region raster=landclass96 -p # we smooth corners of area features r.to.vect -s input=landclass96 output=my_landclass96 type=area v.colors my_landclass96 color=random
Modified program for smoothed lines:
David Satnik,
Central Washington University
Updated 2001 by Andrea Aime, Modena, Italy
Update
Original r.to.sites, r.line and r.poly merged and updated to 5.7 by Radim Blazek
Available at: r.to.vect source code (history)
Latest change: Tuesday Apr 23 10:45:15 2024 in commit: f8115df1219e784a7136e7609f4c9bb16d928e2f
Main index | Raster index | Topics index | Keywords index | Graphical index | Full index
© 2003-2024 GRASS Development Team, GRASS GIS 8.4.1dev Reference Manual