[GRASS-SVN] r34710 - in grass-addons: display/d.frame.split raster/r.colors.tools

svn_grass at osgeo.org svn_grass at osgeo.org
Thu Dec 4 05:10:30 EST 2008


Author: hamish
Date: 2008-12-04 05:10:30 -0500 (Thu, 04 Dec 2008)
New Revision: 34710

Added:
   grass-addons/display/d.frame.split/d.frame.split.README
Removed:
   grass-addons/raster/r.colors.tools/r.colors.stddev
Log:
mo' move to main cleanup

Added: grass-addons/display/d.frame.split/d.frame.split.README
===================================================================
--- grass-addons/display/d.frame.split/d.frame.split.README	                        (rev 0)
+++ grass-addons/display/d.frame.split/d.frame.split.README	2008-12-04 10:10:30 UTC (rev 34710)
@@ -0,0 +1 @@
+module moved into main GRASS source tree as d.split.frame 12/2008

Deleted: grass-addons/raster/r.colors.tools/r.colors.stddev
===================================================================
--- grass-addons/raster/r.colors.tools/r.colors.stddev	2008-12-04 10:06:59 UTC (rev 34709)
+++ grass-addons/raster/r.colors.tools/r.colors.stddev	2008-12-04 10:10:30 UTC (rev 34710)
@@ -1,145 +0,0 @@
-#!/bin/sh
-############################################################################
-#
-# MODULE:       r.colors.stddev
-# AUTHOR:       M. Hamish Bowman, Dept. Marine Science, Otago Univeristy,
-#                 New Zealand
-# PURPOSE:      
-#
-# COPYRIGHT:    (c) 2007 Hamish Bowman, and the GRASS Development Team
-#               This program is free software under the GNU General Public
-#               License (>=v2). Read the file COPYING that comes with GRASS
-#               for details.
-#
-#############################################################################
-
-#%Module
-#% description: Set color rules based on stddev from a map's mean value.
-#% keywords: raster
-#%End
-#% option
-#% key: input
-#% type: string
-#% gisprompt: old,cell,raster
-#% key_desc: name
-#% description: Name of input raster map 
-#% required: yes
-#%end
-#%flag
-#% key: b
-#% description: Color using standard deviation bands
-#%end
-#%flag
-#% key: z
-#% description: Force center at zero
-#%end
-
-
-if  [ -z "$GISBASE" ] ; then
-    echo "You must be in GRASS GIS to run this program." 1>&2
-    exit 1
-fi
-
-if [ "$1" != "@ARGS_PARSED@" ] ; then
-    exec g.parser "$0" "$@"
-fi
-
-
-eval `r.univar -g "$GIS_OPT_INPUT"`
-# $? is result of the eval not r.univar (???)
-#if [ $? -ne 0 ] ; then
-#   echo "ERROR: Problem running r.univar" 1>&2
-#   exit 1
-#fi
-
-
-if [ $GIS_FLAG_Z -eq 0 ] ; then
-
-   MEAN_MINUS_2STDEV=`echo "$mean $stddev" | awk '{print $1 - 2*$2}'`
-   MEAN_PLUS_2STDEV=`echo "$mean $stddev" | awk '{print $1 + 2*$2}'`
-
-   if [ $GIS_FLAG_B -eq 0 ] ; then
-     # smooth free floating blue/white/red
-     r.colors "$GIS_OPT_INPUT" color=rules << EOF
-       0% blue
-       $MEAN_MINUS_2STDEV blue
-       $mean white
-       $MEAN_PLUS_2STDEV red
-       100% red
-EOF
-   else
-     # banded free floating  black/red/yellow/green/yellow/red/black
-     MEAN_MINUS_1STDEV=`echo "$mean $stddev" | awk '{print $1 - $2}'`
-     MEAN_MINUS_3STDEV=`echo "$mean $stddev" | awk '{print $1 - 3*$2}'`
-     MEAN_PLUS_1STDEV=`echo "$mean $stddev" | awk '{print $1 + $2}'`
-     MEAN_PLUS_3STDEV=`echo "$mean $stddev" | awk '{print $1 + 3*$2}'`
-
-     # reclass with labels only works for category (integer) based maps
-     #r.reclass input="$GIS_OPT_INPUT" output="${GIS_OPT_INPUT}.stdevs" << EOF
-
-     # >3 S.D. outliers colored black so they show up in d.histogram w/ white background
-     r.colors "$GIS_OPT_INPUT" color=rules << EOF
-       0% black
-       $MEAN_MINUS_3STDEV black
-       $MEAN_MINUS_3STDEV red
-       $MEAN_MINUS_2STDEV red
-       $MEAN_MINUS_2STDEV yellow
-       $MEAN_MINUS_1STDEV yellow
-       $MEAN_MINUS_1STDEV green
-       $MEAN_PLUS_1STDEV green
-       $MEAN_PLUS_1STDEV yellow
-       $MEAN_PLUS_2STDEV yellow
-       $MEAN_PLUS_2STDEV red
-       $MEAN_PLUS_3STDEV red
-       $MEAN_PLUS_3STDEV black
-       100% black
-EOF
-   fi
-
-
-else 
-   # data centered on 0  (e.g. map of deviations)
-   r.mapcalc "r_col_stdev_abs_$$ = abs($GIS_OPT_INPUT)"
-   eval `r.info -r "r_col_stdev_abs_$$"`
-
-   # current r.univar truncates percentage to the base integer
-   STDDEV2=`r.univar -eg "r_col_stdev_abs_$$" perc=95.4500 | grep ^percentile | cut -f2 -d=`
-
-   if [ $GIS_FLAG_B -eq 0 ] ; then
-     # zero centered smooth blue/white/red
-     r.colors "$GIS_OPT_INPUT" color=rules << EOF
-       -$max blue
-       -$STDDEV2 blue
-       0 white
-       $STDDEV2 red
-       $max red
-EOF
-   else
-     # zero centered banded  black/red/yellow/green/yellow/red/black
-
-     # current r.univar truncates percentage to the base integer
-     STDDEV1=`r.univar -eg "r_col_stdev_abs_$$" perc=68.2689 | grep ^percentile | cut -f2 -d=`
-     STDDEV3=`r.univar -eg "r_col_stdev_abs_$$" perc=99.7300 | grep ^percentile | cut -f2 -d=`
-
-     # >3 S.D. outliers colored black so they show up in d.histogram w/ white background
-     r.colors "$GIS_OPT_INPUT" color=rules << EOF
-       -$max black
-       -$STDDEV3 black
-       -$STDDEV3 red
-       -$STDDEV2 red
-       -$STDDEV2 yellow
-       -$STDDEV1 yellow
-       -$STDDEV1 green
-       $STDDEV1 green
-       $STDDEV1 yellow
-       $STDDEV2 yellow
-       $STDDEV2 red
-       $STDDEV3 red
-       $STDDEV3 black
-       $max black
-EOF
-   fi
-
-   g.remove rast="r_col_stdev_abs_$$" --quiet
-fi
-



More information about the grass-commit mailing list