[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