r.li.pielou
Calculates Pielou's diversity index on a raster map
r.li.pielou input=name config=name output=name [--overwrite] [--verbose] [--quiet] [--qq] [--ui]
Example:
r.li.pielou input=name config=name output=name
grass.script.run_command("r.li.pielou", input, config, output, overwrite=False, verbose=False, quiet=False, superquiet=False)
Example:
gs.run_command("r.li.pielou", input="name", config="name", output="name")
Parameters
input=name [required]
Name of input raster map
config=name [required]
Configuration file
output=name [required]
Name for output raster map
--overwrite
Allow output files to overwrite existing files
--help
Print usage summary
--verbose
Verbose module output
--quiet
Quiet module output
--qq
Very quiet module output
--ui
Force launching GUI dialog
input : str, required
Name of input raster map
Used as: input, raster, name
config : str, required
Configuration file
Used as: input, file, name
output : str, required
Name for output raster map
Used as: output, raster, name
overwrite: bool, optional
Allow output files to overwrite existing files
Default: False
verbose: bool, optional
Verbose module output
Default: False
quiet: bool, optional
Quiet module output
Default: False
superquiet: bool, optional
Very quiet module output
Default: False
DESCRIPTION
r.li.pielou calculates the "Pielou's diversity index" as:
with:
- J: Pielou eveness index
- H': Shannon index
- S: number of classes
NOTES
Do not use absolute path names for the config and output
file/map parameters. If the "moving window" method was selected in
g.gui.rlisetup, then the output will be a raster map, otherwise an
ASCII file will be generated in the folder
C:\Users\userxy\AppData\Roaming\GRASS8\r.li\output\
(MS-Windows) or
$HOME/.grass8/r.li/output/
(GNU/Linux).
If the input raster contains NULL value cells, r.li.pielou returns NULL for these cells.
EXAMPLES
To calculate Pielou's diversity index on map my_map, using my_conf configuration file (previously defined with g.gui.rlisetup) and saving results in my_out, run:
r.li.pielou input=my_map conf=my_conf output=my_out
Forest map (Spearfish sample dataset) example:
g.region raster=landcover.30m -p
r.mapcalc "forests = if(landcover.30m >= 41 && landcover.30m <= 43,1,null())"
r.li.pielou input=forests conf=movwindow7 out=forests_pielou_mov7
r.univar forests_pielou_mov7
Forest map (North Carolina sample dataset) example:
g.region raster=landclass96 -p
r.mapcalc "forests = if(landclass96 == 5, 1, null() )"
r.li.pielou input=forests conf=movwindow7 out=forests_pielou_mov7
# verify
r.univar forests_pielou_mov7
r.to.vect input=forests output=forests type=area
d.mon wx0
d.rast forests_pielou_mov7
d.vect forests type=boundary
SEE ALSO
r.li (package overview), g.gui.rlisetup
REFERENCES
McGarigal, K., and B. J. Marks. 1995. FRAGSTATS: spatial pattern analysis program for quantifying landscape structure. USDA For. Serv. Gen. Tech. Rep. PNW-351. (PDF)
AUTHORS
Luca Delucchi and Duccio Rocchini, Fondazione E. Mach (Italy), based on the r.li.shannon code developed by Serena Pallecchi, student of Computer Science University of Pisa (Italy)
SOURCE CODE
Available at: r.li.pielou source code
(history)
Latest change: Friday Feb 07 19:16:09 2025 in commit a82a39f