.TH r.li.mps 1grass "" "GRASS 6.4.4" "Grass User's Manual"
.SH NAME
\fI\fBr.li.mps\fR\fR  - Calculates mean patch size index on a raster map, using a 4 neighbour algorithm
.SH KEYWORDS
raster, landscape structure analysis, patch index
.SH SYNOPSIS
\fBr.li.mps\fR
.br
\fBr.li.mps help\fR
.br
\fBr.li.mps map\fR=\fIname\fR \fBconf\fR=\fIname\fR \fBoutput\fR=\fIname\fR  [\-\-\fBoverwrite\fR]  [\-\-\fBverbose\fR]  [\-\-\fBquiet\fR] 
.SS Flags:
.IP "\fB\-\-overwrite\fR" 4m
.br
Allow output files to overwrite existing files
.IP "\fB\-\-verbose\fR" 4m
.br
Verbose module output
.IP "\fB\-\-quiet\fR" 4m
.br
Quiet module output
.PP
.SS Parameters:
.IP "\fBmap\fR=\fIname\fR" 4m
.br
Name of input raster map
.IP "\fBconf\fR=\fIname\fR" 4m
.br
Configuration file
.IP "\fBoutput\fR=\fIname\fR" 4m
.br
Name for output raster map
.PP
.SH DESCRIPTION
\fIr.li.mps\fR (mean patch size) calculates the mean size of the 
patches in the sampling area as:
.br
\fC
.DS
.br
MPS = A / Npatch
.br
.DE
\fR
with:
.RS
.IP
\fBA\fR:sampling area size without null cells
.IP
\fBNpatch\fR: number of patches
.RE
The unit of the mean size is hectare.
.PP
This index is calculated using a 4 neighbour algorithm, diagonal cells 
are ignored when tracing a patch.
.SH NOTES
Do not use absolute path names for the \fBoutput\fR map/file.
If the "moving window" method was selected in \fBr.li.setup\fR, then the
output will be a raster map, otherwise an ASCII file will be generated in
the $HOME/.r.li/output/ folder.
.PP
If the input raster contains NULL value cells, \fIr.li.mps\fR
returns 0 for these cells. 
.br
If you want to change these 0 values to NULL, run subsequently on the resulting map:
\fC
.DS
.br
r.null setnull=0 map=my_map
.br
.DE
\fR
after index calculation.
.SH EXAMPLES
To calculate mean path size index on map \fImy_map\fR, using
\fImy_conf\fR configuration file (previously defined with
\fIr.li.setup\fR) and saving results in \fImy_out\fR, run:
\fC
.DS
.br
r.li.mps map=my_map conf=my_conf output=my_out
.br
.DE
\fR
.PP
Forest map (Spearfish sample dataset) example:
\fC
.DS
.br
g.region rast=landcover.30m \-p
.br
r.mapcalc "forests = if(landcover.30m >= 41 && landcover.30m <= 43,1,null())"
.br
r.li.mps map=forests conf=movwindow7 out=forests_mps_mov7
.br
r.univar forests_mps_mov7
.br
.DE
\fR
.SH SEE ALSO
\fI
r.li - package overview
.br
r.li.setup
\fR
.SH 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)
.SH AUTHORS
Serena Pallecchi student of Computer Science University of Pisa (Italy).
.br
Commission from Faunalia Pontedera (PI), Italy (www.faunalia.it)
.PP
\fILast changed: $Date: 2014-02-19 01:26:04 +0100 (Wed, 19 Feb 2014) $\fR
.PP
Full index
.PP
© 2003-2014 GRASS Development Team