Re: [GRASS-user] r.mapcalc and producing a max value map

2008-07-14 Thread Michael Barton



C. Michael Barton, Professor of Anthropology
Director of Graduate Studies
School of Human Evolution & Social Change
Center for Social Dynamics & Complexity
Arizona State University

Phone: 480-965-6262
Fax: 480-965-7671
www: 



On Jul 14, 2008, at 6:32 AM, <[EMAIL PROTECTED]> wrote:


Date: Mon, 14 Jul 2008 15:22:03 +0200
From: Christian Schwartze <[EMAIL PROTECTED]>
Subject: [GRASS-user] (no subject)
To: grass-user@lists.osgeo.org
Message-ID: <[EMAIL PROTECTED]>
Content-Type: text/plain; charset=ISO-8859-1

Dear GRASS users,

could somebody help me to find some solution for the following task?
I have a raster map representing many distinct areas (MAP_1) and a  
second int
raster with specific values (MAP_2). The result map should contain  
the same
areas of MAP_1 but with the max MAP_2 value of each area as IDs. I  
tried some
approaches with r.mapcalc but I don't know how to code the region  
aspect of

MAP_1 into r.mapcalc.

MAP_1   MAP_2   ---> MAP_MAX
-   -   ---

32  5266723974  95
111222  1429784155  ---> 777555
111322  4152763923  777955
12  4125718443  75



Christian,

I think what you need is the following:

r.statistics base=map_1 cover=map_2 method=max output=map_max

Michael
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.mapcalc and producing a max value map

2008-07-14 Thread G. Allegri
Take a loot to r.statistics:
http://grass.osgeo.org/grass64/manuals/html64_user/r.statistics.html

2008/7/14 Christian Schwartze <[EMAIL PROTECTED]>:
> Dear GRASS users,
>
> could somebody help me to find some solution for the following task?
> I have a raster map representing many distinct areas (MAP_1) and a second int
> raster with specific values (MAP_2). The result map should contain the same
> areas of MAP_1 but with the max MAP_2 value of each area as IDs. I tried some
> approaches with r.mapcalc but I don't know how to code the region aspect of
> MAP_1 into r.mapcalc.
>
>MAP_1   MAP_2   --->MAP_MAX
>-   -   ---
>
>32  5266723974  95
>111222  1429784155  --->777555
>111322  4152763923  777955
>12  4125718443  75
>
> Thanks a lot for any help!
> Regards,
> Christian.
>
>
>
> 
> This mail was sent through http://webmail.uni-jena.de
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.mapcalc and producing a max value map

2008-07-14 Thread Christian Schwartze
Dear GRASS users,

could somebody help me to find some solution for the following task?
I have a raster map representing many distinct areas (MAP_1) and a second int
raster with specific values (MAP_2). The result map should contain the same
areas of MAP_1 but with the max MAP_2 value of each area as IDs. I tried some
approaches with r.mapcalc but I don't know how to code the region aspect of
MAP_1 into r.mapcalc.

MAP_1   MAP_2   --->MAP_MAX
-   -   ---

32  5266723974  95
111222  1429784155  --->777555
111322  4152763923  777955
12  4125718443  75

Thanks a lot for any help!
Regards,
Christian.




This mail was sent through http://webmail.uni-jena.de
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user