[GRASS-user] r.mapcalc appears to have changed behavior
Jim Mochel
jim.mochel.hivol at gmail.com
Tue Sep 27 06:45:58 EDT 2011
Thank you. !
I am not sure how I did that but I am sure it seemed like a good idea at the
time :-)
Jim
On Tue, Sep 27, 2011 at 1:02 AM, Michael Barton <michael.barton at asu.edu>wrote:
> It looks like for the 2nd calculation, your computational region is set to
> extents of 0W - 1E, 0S - 1N, with a resolution of 1.
>
> So r.mapcalc is only operating on a single 1x1 cell that is outside the
> area of your raster map.
>
> Michael
> ____________________
> C. Michael Barton
> Director, Center for Social Dynamics & Complexity
> Professor of Anthropology, School of Human Evolution & Social Change
> Arizona State University
>
> voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC)
> fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC)
> www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu
>
>
>
>
>
>
>
>
>
>
>
> On Sep 27, 2011, at 11:50 PM, <grass-user-request at lists.osgeo.org> <
> grass-user-request at lists.osgeo.org> wrote:
>
> > Date: Mon, 26 Sep 2011 23:19:11 -0400
> > From: Jim Mochel <jim.mochel.hivol at gmail.com>
> > Subject: [GRASS-user] r.mapcalc appears to have changed behavior
> > To: grass-user at lists.osgeo.org
> > Message-ID:
> > <CAE=Y9-Q7tFZ=cJ4TC6nceYA5B=
> R9E7RRwZbVfF_fYbMNTEpC6A at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > I had made a number of scripts that worked fine in the 6.4.0 build 6 but
> > appear to have problems with the same code in later revisions:
> >
> > The method used to take about 6-8 secs and produce about what you would
> > expect on the screen. Now the same calculation takes 0 seconds the file
> is
> > radically changed. See the before and after below.
> >
> > Any help would be greatly appreciated.
> > Thanks
> > Jim Mochel
> >
> > r.mapcalc na_dem_sub_20 = na_dem_input - 20
> >
> >
> >
> +----------------------------------------------------------------------------+
> > | Layer: na_dem Date: Sat Sep 24 17:37:59 2011
> > |
> > | Mapset: PERMANENT Login of Creator: jmochel
> > |
> > | Location: wgs84
> > |
> > | DataBase: /home/jmochel/grassdata
> > |
> > | Title: ( na_dem )
> > |
> > | Timestamp: none
> > |
> >
> |----------------------------------------------------------------------------|
> > |
> > |
> > | Type of Map: raster Number of Categories: 65535
> > |
> > | Data Type: CELL
> > |
> > | Rows: 8384
> > |
> > | Columns: 9102
> > |
> > | Total Cells: 76311168
> > |
> > | Projection: x,y
> > |
> > | N: 4384500 S: -3999500 Res: 1000
> > |
> > | E: 4639500 W: -4462500 Res: 1000
> > |
> > | Range of data: min = 0 max = 65535
> > |
> > |
> > |
> > | Data Description:
> > |
> > | generated by r.in.gdal
> > |
> > |
> > |
> > | Comments:
> > |
> > | r.in.gdal input="/home/jmochel/gis/raw-data/hydo1k/na_dem.bil"
> outpu\
> > |
> > | t="na_dem"
> > |
> > |
> > |
> >
> +----------------------------------------------------------------------------+
> >
> >
> > r.info map=na_dem_sub_20
> >
> >
> +----------------------------------------------------------------------------+
> > | Layer: na_dem_sub_20 Date: Mon Sep 26 23:06:03 2011
> > |
> > | Mapset: PERMANENT Login of Creator: jmochel
> > |
> > | Location: wgs84
> > |
> > | DataBase: /home/jmochel/grassdata
> > |
> > | Title: ( na_dem_sub_20 )
> > |
> > | Timestamp: none
> > |
> >
> |----------------------------------------------------------------------------|
> > |
> > |
> > | Type of Map: raster Number of Categories: 610
> > |
> > | Data Type: CELL
> > |
> > | Rows: 1
> > |
> > | Columns: 1
> > |
> > | Total Cells: 1
> > |
> > | Projection: x,y
> > |
> > | N: 1 S: 0 Res: 1
> > |
> > | E: 1 W: 0 Res: 1
> > |
> > | Range of data: min = 610 max = 610
> > |
> > |
> > |
> > | Data Description:
> > |
> > | generated by r.mapcalc
> > |
> > |
> > |
> > | Comments:
> > |
> > | na_dem_input at PERMANENT - 20
> > |
> > |
> > |
> >
> +----------------------------------------------------------------------------+
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20110927/5c3fad27/attachment-0001.html
More information about the grass-user
mailing list