1r.li.patchdensity(1) Grass User's Manual r.li.patchdensity(1)
2
3
4
6 r.li.patchdensity - Calculates patch density index on a raster map,
7 using a 4 neighbour algorithm
8
10 raster, landscape structure analysis, patch index
11
13 r.li.patchdensity
14 r.li.patchdensity help
15 r.li.patchdensity map=name conf=string output=name [--overwrite]
16 [--verbose] [--quiet]
17
18 Flags:
19 --overwrite
20 Allow output files to overwrite existing files
21
22 --verbose
23 Verbose module output
24
25 --quiet
26 Quiet module output
27
28 Parameters:
29 map=name
30 Name of input raster map
31
32 conf=string
33 Configuration file
34
35 output=name
36 Name for output raster map
37
39 r.li.patchdensity is a patch index for r.li suite. It calculates the
40 function:
41 f(sample_area) = (Patch_Number/Area) * 1000000
42 that is 1000000 by number of patch for area unit. This index is cal‐
43 culated using a 4 neighbour algorithm.
44
46 To calculate patch density index on map my_map, using my_conf configu‐
47 ration file (previously defined with r.li.setup) and saving results in
48 the raster map my_out, run:
49 r.li.patchdensity map=my_map conf=my_conf out=my_out
50
51
52 Example for Spearfish forest areas:
53 g.region rast=landcover.30m -p
54 # extract forested areas:
55 r.category landcover.30m
56 r.mapcalc "forests=if(landcover.30m >= 41 && landcover.30m
57 <=43,1,null())"
58 # patch density (7x7 moving window defined in r.li.setup):
59 r.li.patchdensity forests conf=movwindow7 out=forests_p_dens7
60 d.rast.leg forests_p_dens7
61 r.to.vect forests out=forests feature=area
62 d.vect forests type=boundary
63
64
66 A map of NULL values is considered to have zero patches.
67 If you want to have null values instead run
68
69 r.null setnull=0 map=my_out
70
71 after index calculation.
72 If raster area is 0, r.li.patchdensity returns -1. This is only possi‐
73 ble
74 if the raster is masked.
75
77 r.li package overview
78 r.li.daemon
79 r.li.setup
80
82 Claudio Porta and Lucio Davide Spano, students of Computer Science
83 University of Pisa (Italy).
84 Commission from Faunalia Pontedera (PI) (www.faunalia.it)
85
87 Please send bugs reports to
88 spano@cli.di.unipi.it,
89 porta@cli.di.unipi.it
90
91
92 Full index
93
94 © 2003-2008 GRASS Development Team
95
96
97
98GRASS 6.3.0 r.li.patchdensity(1)