3
3
#========================
4
4
# North Carolina tests
7
G7RLI=$HOME/.grass7/r.li
9
# switch GRASS GIS version here:
12
echo "SAMPLINGFRAME 0|0|1|1
13
SAMPLEAREA 0.0|0.0|1.0|1.0" > $G7RLI/landcover_whole_whole # GRASS7
15
echo "SAMPLINGFRAME 0|0|1|1
16
SAMPLEAREA 0.0|0.0|1.0|1.0" > $G6RLI/history/landcover_whole_whole # GRASS6
6
G_RLI=$HOME/.grass7/r.li
9
echo "SAMPLINGFRAME 0|0|1|1
10
SAMPLEAREA 0.0|0.0|1.0|1.0" > $G_RLI/landcover_whole_whole # GRASS7
19
12
# North Carolina location:
20
13
export GRASS_OVERWRITE=1
21
g.region rast=landclass96 -p
14
g.region raster=landclass96 -p
22
15
RASTER_MAP=landclass96
23
#r.to.vect in=basin_50K out=basin_50K feature=area
16
#r.to.vect in=basin_50K output=basin_50K feature=area
25
18
echo "-------------------------"
26
r.li.dominance $RASTER_MAP conf=landcover_whole_whole out=forests_dominance_whole
19
r.li.dominance $RASTER_MAP conf=landcover_whole_whole output=forests_dominance_whole
27
20
cat $G_RLI/output/forests_dominance_whole
29
22
echo "-------------------------"
30
r.li.edgedensity $RASTER_MAP conf=landcover_whole_whole out=forests_edgedens_whole
23
r.li.edgedensity $RASTER_MAP conf=landcover_whole_whole output=forests_edgedens_whole
31
24
cat $G_RLI/output/forests_edgedens_whole
33
26
echo "-------------------------"
34
r.li.mpa $RASTER_MAP conf=landcover_whole_whole out=forests_mpa_whole
27
r.li.mpa $RASTER_MAP conf=landcover_whole_whole output=forests_mpa_whole
35
28
cat $G_RLI/output/forests_mpa_whole
37
30
echo "-------------------------"
38
r.li.mps $RASTER_MAP conf=landcover_whole_whole out=forests_mps_whole
31
r.li.mps $RASTER_MAP conf=landcover_whole_whole output=forests_mps_whole
39
32
cat $G_RLI/output/forests_mps_whole
41
34
echo "-------------------------"
42
r.li.padcv $RASTER_MAP conf=landcover_whole_whole out=forests_padcv_whole
35
r.li.padcv $RASTER_MAP conf=landcover_whole_whole output=forests_padcv_whole
43
36
cat $G_RLI/output/forests_padcv_whole
45
38
echo "-------------------------"
46
r.li.padrange $RASTER_MAP conf=landcover_whole_whole out=forests_padrange_whole
39
r.li.padrange $RASTER_MAP conf=landcover_whole_whole output=forests_padrange_whole
47
40
cat $G_RLI/output/forests_padrange_whole
49
42
echo "-------------------------"
50
r.li.padsd $RASTER_MAP conf=landcover_whole_whole out=forests_padsd_whole
43
r.li.padsd $RASTER_MAP conf=landcover_whole_whole output=forests_padsd_whole
51
44
cat $G_RLI/output/forests_padsd_whole
53
46
echo "-------------------------"
54
r.li.patchdensity $RASTER_MAP conf=landcover_whole_whole out=forests_p_dens_whole
47
r.li.patchdensity $RASTER_MAP conf=landcover_whole_whole output=forests_p_dens_whole
55
48
cat $G_RLI/output/forests_p_dens_whole
57
50
echo "-------------------------"
58
r.li.patchnum $RASTER_MAP conf=landcover_whole_whole out=forests_patchnum_whole
51
r.li.patchnum $RASTER_MAP conf=landcover_whole_whole output=forests_patchnum_whole
59
52
cat $G_RLI/output/forests_patchnum_whole
61
54
echo "-------------------------"
62
r.li.pielou $RASTER_MAP conf=landcover_whole_whole out=forests_pielou_whole
55
r.li.pielou $RASTER_MAP conf=landcover_whole_whole output=forests_pielou_whole
63
56
cat $G_RLI/output/forests_pielou_whole
65
58
echo "-------------------------"
66
r.li.renyi $RASTER_MAP conf=landcover_whole_whole out=forests_renyi_whole7_a06 alpha=0.6
59
r.li.renyi $RASTER_MAP conf=landcover_whole_whole output=forests_renyi_whole7_a06 alpha=0.6
67
60
cat $G_RLI/output/forests_renyi_whole7_a06
69
62
echo "-------------------------"
70
r.li.richness $RASTER_MAP conf=landcover_whole_whole out=forests_richness_whole
63
r.li.richness $RASTER_MAP conf=landcover_whole_whole output=forests_richness_whole
71
64
cat $G_RLI/output/forests_richness_whole
73
66
echo "-------------------------"
74
r.li.shannon $RASTER_MAP conf=landcover_whole_whole out=forests_shannon_whole
67
r.li.shannon $RASTER_MAP conf=landcover_whole_whole output=forests_shannon_whole
75
68
cat $G_RLI/output/forests_shannon_whole
77
70
echo "====== End of r.li tests ========================"