v.random(1)

NAME
v.random - Generates random 2D/3D vector points.

KEYWORDS
vector, sampling, statistics, random, point pattern, stratified random sampling, level1

SYNOPSIS
v.random
v.random --help
v.random [-zab] output=name npoints=integer [restrict=name] [layer=string] [cats=range] [where=sql_query] [zmin=float] [zmax=float] [seed=integer] [column=name] [column_type=string] [--overwrite] [--help] [--verbose] [--quiet] [--ui]

  Flags:

-z
    Create 3D output

-a
    Generate n points for each individual area (requires restrict parameter)

-b
    Do not build topology
    Advantageous when handling a large number of points

--overwrite
    Allow output files to overwrite existing files

--help
    Print usage summary

--verbose
    Verbose module output

--quiet
    Quiet module output

--ui
    Force launching GUI dialog

  Parameters:

output=name [required]
    Name for output vector map

npoints=integer [required]
    Number of points to be created

restrict=name
    Name of input vector map
    Restrict points to areas in input vector

layer=string
    Layer number or name (’-1’ for all layers)
    A single vector map can be connected to multiple database tables. This number determines which table to use. When used with direct OGR access this is the layer name.
    Default: -1

cats=range
    Category values
    Example: 1,3,7-9,13

where=sql_query
    WHERE conditions of SQL statement without ’where’ keyword
    Example: income < 1000 and population >= 10000

zmin=float
    Minimum z height (needs -z flag or column name)
    Default: 0.0

zmax=float
    Maximum z height (needs -z flag or column name)
    Default: 0.0

seed=integer
    The seed to initialize the random generator. If not set the process ID is used

column=name
    Name of column for z values
    Writes z values to column

column_type=string
    Type of column for z values
    Options: integer, double precision
    Default: double precision

DESCRIPTION
v.random randomly generates vector points within the current region using the selected random number generator.

v.random can generate also 3D vector points or write random value to the attribute table. Point height range or attribute value range is controlled by specifying zmin and zmax options. Both z values are included in range (zmin <= z <= zmax). Generated random attribute value type can be controlled by column_type. Use integer column type for integers and double precision for floating point numbers. Integer values are calculated by rounding random floating point number.

To produce repeatable results a random seed can be set using the option seed.

  Restriction to vector areas

If an restrict vector map with areas is specified, the location of random points is restricted to the selected areas. By default, the requested number of points are distributed across all areas.

If the -a flag is given, the requested number of points is generated for each individual area. For example, if 20 points should be generated and the input map has 100 individual areas, 2000 points will be generated in total.

Attributes attached to restrict vector map are also transferred if the layer parameter is defined > 0, see example below.

NOTES
Importantly, attributes will only be transferred if layer > 0 (e.g., layer=1).

EXAMPLES
All examples are based on the North Carolina sample dataset.

  Generating random points in 2D
Generate 20 random points with binary attributes (only 0 or 1):

v.random output=binary_random npoints=20 zmin=0 zmax=1 column=’binary’
v.db.select binary_random
cat|binary
1|0.63495
2|0.233421
3|0.489302
4|0.748264
5|0.505556
6|0.32975
[...]
v.univar -d binary_random
Calculating geometric distances between 20 primitives...
[...]
minimum: 148.515
maximum: 16572.8
[...]

  Generating random points in 2D with binary attributes
Generate 20 random points with binary attributes (only 0 or 1):

v.random output=binary_random npoints=20 zmin=0 zmax=1 column=’binary’ column_type=integer
v.db.select binary_random
cat|binary
1|0
2|0
3|0
4|0
5|1
6|0
[...]

  Generating random points in 3D
Generate 20 random 3D points using a specific random seed:

v.random seed=52 output=height_random npoints=40 zmin=110 zmax=170 -z
v.univar -d height_random
Calculating geometric distances between 40 primitives...
[...]
minimum: 334.889
maximum: 18351.9
range: 18017
sum: 5.38425e+06
mean: 7266.2
mean of absolute values: 7266.2
population standard deviation: 3563.95
[...]
skewness: 0.34703

Random points with different X, Y, and Z coordinates

  Generating random points in selected areas
Generate 3 random points only in selected areas ("RALEIGH" related ZIP code areas):

v.random restrict=zipcodes_wake output=zipcodes_local_random_n3 npoints=3 where="ZIPNAME = ’RALEIGH’" -a
# visualization
d.mon wx0
d.vect zipcodes_wake
d.vect zipcodes_wake fcolor=yellow where="ZIPNAME = ’RALEIGH’"
d.vect zipcodes_local_random_n3 color=red icon=basic/circle

Fixed number of random points generated in selected areas

  Generating random adjacent areas
To generate random adjacent areas, first the centroids are generated as points, then a triangulation is run (North Carolina sample dataset:

g.region vector=nc_state
v.random output=randpoints6k npoints=6000
v.voronoi input=randpoints6k output=randareas6k
v.info -t randareas6k
v.category randareas6k option=print
# plot vector areas
d.mon wx0
d.vect randareas6k -c

Random adjacent areas from random points (here: used as centroids)

To eventually obtain isolated areas, selected areas can be extracted with v.extract.

These vector areas can also be rasterized:

# rasterize areas
# note: rastermaps must result in at least 6k pixel in this example
g.region vector=nc_state res=500 -p -a
v.to.rast randareas6k out=randareas6k use=cat
r.colors randareas6k color=random
d.rast randareas6k

  Random sampling from raster map
Generate 20 random samples from a raster map:

g.region -p raster=elevation
v.random output=random_samples npoints=20
v.db.addtable map=random_samples columns=’cat INTEGER, sample DOUBLE PRECISION’
v.what.rast map=random_samples raster=elevation column=sample
v.db.select random_samples
cat|sample
1|103.9935
2|129.1266
3|96.01388
[...]

  Random sampling from vector map
Generate 20 random points and sample attribute data from geology (vector) map:

g.region -p vector=geology
v.random output=random_samples npoints=20
v.db.addtable map=random_samples columns=’cat integer, geology varchar(100)’
v.what.vect map=random_samples column=geology query_map=geology query_layer=1 query_column=GEO_NAME
v.db.select random_samples
cat|geology
1|PzZm
2|
3|Zatm
[...]

  Stratified random sampling: Random sampling from vector map by attribute
Generate 20 random points restricted to forested areas:

g.region -p raster=landclass96
r.to.vect -v input=landclass96 output=landclass96 type=area
v.random restrict=landclass96 output=random_samples npoints=20 where="label = ’forest’" layer=1
v.db.select map=random_samples
cat|landclass96_cat|landclass96_label
1|5|forest
2|5|forest
3|5|forest
...

Random points only sampled in forested areas (stratified random sampling)

  Stratified random sampling: Random sampling from vector map with spatial constraints
Generating n points for each individual area: in this example two random points in each water body:

g.region -p raster=landclass96
r.to.vect -v input=landclass96 output=landclass96 type=area
v.random restrict=landclass96 output=random_samples npoints=2 where="label = ’water’" layer=1 -a

Two random points sampled in each individual water body (stratified random sampling)

SEE ALSO
g.region, r.random, v.db.addtable, v.perturb, v.sample, v.univar, v.what.rast, v.what.vect

SQL support in GRASS GIS

AUTHOR
James Darrell McCauley <darrell@mccauley-usa.com>,
when he was at: Agricultural Engineering Purdue University

SOURCE CODE

Available at: v.random source code (history)

Latest change: Mon Oct 26 11:54:36 2020 in commit: 3b6d257bdfc18a58dd42c5ab06c69ada99c56a24

Main index | Vector index | Topics index | Keywords index | Graphical index | Full index

© 2003-2022 GRASS Development Team, GRASS GIS 7.8.8dev Reference Manual