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