TODO list
=========

# TODO: change the examples to North Carolina wherever possible:

# Spearfish tests
g.region rast=landcover.30m -p
r.mapcalc "forests = if(landcover.30m >= 41 && landcover.30m <= 43,1,null())"
d.mon x0
d.rast forests

# r.li.setup:
# created as described in EXAMPLES, moving window, in r.li.setup/description.html:
echo "SAMPLINGFRAME 0|0|1|1
SAMPLEAREA -1|-1|0.015021459227467811|0.011058451816745656
MOVINGWINDOW" > $HOME/.r.li/history/movwindow7

r.li.patchdensity forests conf=movwindow7 out=forests_p_dens7
d.rast.leg forests_p_dens7

r.to.vect forests out=forests feature=area
d.vect forests type=boundary

## CHECK THIS:

# MASK test
g.copy rast=fields,MASK
r.li.patchdensity forests conf=movwindow7 out=forests_p_dens7mask --o
d.rast.leg forests_p_dens7mask
# -> no negative values! but MASK is respected

# zero data test
r.mapcalc forests=0
r.li.patchdensity forests conf=movwindow7 out=forests_p_dens7 --o
d.rast.leg forests_p_dens7mask
r.info -r forests_p_dens7
min=19.841270
max=19.841270
# -> all cells are 19.841270 ?!


========================
North Carolina tests

# Tests with return an ASCII output file
sh ./r.li.testing_nc_asc.sh

---------
echo "SAMPLINGFRAME 0|0|1|1
SAMPLEAREA 0.0|0.0|1.0|1.0" > $HOME/.r.li/history/landsat_test

g.region rast=lsat7_2002_40 -p
r.li.shannon input=lsat7_2002_40 conf=landsat_test out=landsat_shannon
--> Result written to ASCII file <$HOME/.r.li/output/landsat_shannon>