<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">
<TITLE>Exotic patch based statistics, a solution...</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>
<P><FONT SIZE=2>Hi list,<BR>
<BR>
just want to share some code that could be useful to some of you.<BR>
<BR>
I was looking for a way to do other patch based statistics than those provided in r.statistics. Using masks and a loop in GRASS was terribly slow (certainly for very patchy areas (in my case some 100000). So this is the fastest I could make it without programming a module myself. I use R, but not the spgrass6 link, from within a bash script.<BR>
<BR>
The script uses a map with all numbered patches as can be provided by r.le.patch or r.le.stats and a cover map which contains the values u want to evaluate with some exotic statistic or function.<BR>
<BR>
In my setup I formatted the output as a rules file so the results could be translated back into an actual map using the patch number map as common factor. Other approaches are also possible, depending on the need of the user.<BR>
<BR>
Hope some people can benefit from it.<BR>
<BR>
Cheers,<BR>
Koen<BR>
<BR>
-- THE BASH SCRIPT ----><BR>
<BR>
#!/bin/bash<BR>
#<BR>
###########################################################################################################################<BR>
#<BR>
# Grass bash script to evaluate the content of patches based upon a patch and cover map...<BR>
#<BR>
# !!! Only works within the GRASS terminal !!!<BR>
#<BR>
#<BR>
# Calculate the a map with numbered patches with r.le.patch or r.le.trace. Use this as input for the patches.number.map.<BR>
# The cover.map map is the map you want to evaluate based upon the delineated patches with another function than those<BR>
# provided by r.statistics.<BR>
#<BR>
# use: sh myscript.sh patch.number.map cover.map (name of the script depends on how you save it ofcourse)<BR>
#<BR>
###########################################################################################################################<BR>
<BR>
r.out.xyz input=$1 output=- fs=' ' | awk '{print $3}' > patches<BR>
r.out.xyz input=$2 output=- fs=' ' | awk '{print $3}' > cover<BR>
<BR>
#<BR>
# Make a two column matrix with first the patch number than the associated cover value<BR>
#<BR>
<BR>
paste patches cover > Rinput<BR>
<BR>
<BR>
###########################################################################################################################<BR>
#<BR>
# R code starts here!<BR>
#<BR>
###########################################################################################################################<BR>
R --slave --vanilla --quiet --no-save <<EEE<BR>
<BR>
<BR>
rdata <- read.table("Rinput")<BR>
<BR>
for (i in 1:nrpatches){<BR>
<BR>
rdatasub <- subset(rdata[,2],rdata[,1]==i)<BR>
<BR>
#<BR>
# complicated patch based statistics go here...<BR>
#<BR>
<BR>
<BR>
}<BR>
EEE<BR>
<BR>
#<BR>
# Close R session within bash GRASS script...<BR>
#<BR>
<BR>
###########################################################################################################################<BR>
<BR>
<BR>
#<BR>
# The output of the R script is best printed according to the rules specifications so you can use the results<BR>
# directly as a reclassify rule within GRASS.<BR>
#<BR>
<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>