Re: [GRASS-user] Problems with grass for OSX
Dear Jose, when did your problem occur? After installation/updating or haphazardly? I had a similar problem with the Python GUI but after I updated all frameworks it worked fine. In case you use the GRASS version from Kyngchaos can you update the GDAL complete framework, free type and cairo and report whether the error still occurs? Regards Ralf Am 10.12.2011 um 05:20 schrieb grass-user-requ...@lists.osgeo.org: > Message: 3 > Date: Fri, 9 Dec 2011 17:04:12 -0200 > From: Jos? Carlos Guerrero Ant?nez > Subject: [GRASS-user] Problems with grass for OSX > To: grass-user@lists.osgeo.org > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > Hi, > > I can not boot properly grass and this is the message that appears in Terminal > window. > > OSX 10.6.8 > > Thanks in advance, > > > torre-b-024:~ joseca_IMac$ > '/Applications/GRASS-6.4.app/Contents/MacOS/grass.sh'; exit > Rebuilding Addon HTML manual pages index... > Rebuilding Addon menu... > Python 2.6.1 found. > > WELCOME TO GRASS Version 6.4.1 2011 > > 1) Have at your side all available GRASS tutorials > > 2) When working on your location, the following materials > are extremely useful: > - A topo map of your area > - Current catalog of available computer maps > > 3) Check the GRASS webpages for feedback mailinglists and more: > http://www.grass-gis.org > http://grass.osgeo.org > > Hit RETURN to continue > > Starting GRASS ... > Traceback (most recent call last): > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 641, in OnSelectLocation >self.listOfLocations[self.lblocations.GetSelection()])) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 590, in UpdateMapsets >self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 629, in GetListOfMapsets > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 651, in EncodeString > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 459, in getdefaultlocale >return _parse_localename(localename) > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 391, in _parse_localename >raise ValueError, 'unknown locale: %s' % localename > ValueError: unknown locale: UTF-8 > Traceback (most recent call last): > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 858, in >GRASSStartUp = StartUp(0) > File > "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", > line 7981, in __init__ > File > "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", > line 7555, in _BootstrapApp > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 843, in OnInit >StartUp.OnBrowse(None) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 703, in OnBrowse >self.OnSetDatabase(event) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 692, in OnSetDatabase >self.OnSelectLocation(None) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 641, in OnSelectLocation >self.listOfLocations[self.lblocations.GetSelection()])) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 590, in UpdateMapsets >self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 629, in GetListOfMapsets > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 651, in EncodeString > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 459, in getdefaultlocale >return _parse_localename(localename) > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 391, in _parse_localename >raise ValueError, 'unknown locale: %s' % localename > ValueError: unknown locale: UTF-8 > Error in GUI startup. If necessary, please > report this error to the GRASS developers. > Switching to text mode now. > Hit RETURN to continue... > > > > -- > > Dr. José Carlos Guerrero > Laboratorio de Desarrollo Sustentable y Gestión Ambiental del Territorio > (LDSGAT) > Instituto de Ecología y Ciencias Ambientales (IECA) > Facultad de Ciencias > Universidad de la República > Montevideo-URUGUAY >
Re: [GRASS-user] How to find the second largest value in a series of rasters.
what about computing the first largest value, making your categorical map, and using it as a mask ? if you mask all cell containing the first largest value, then the second largest value will become the first largest value of the masked raster serie ... my 2 ct Sylvain 2011/12/10 Marcello Gorini > > > Marcello: > > >> > My problem is that I also need to find the *second largest* value and >> the >> > corresponding raster number which contains the second largest value. >> > >> > I am doing that by iterating over all classes through a shell script, >> but it >> > obviously takes much more time than using a simple r.series command. >> > >> > Can anyone share any idea on how to accomplish that using r.series or >> > something similar? >> >> > Glynn: > > >> There isn't an efficient way to do this using existing tools. If you >> need the efficiency, I suggest adding a "second largest" aggregate to >> lib/stats (based upon c_max.c and c_maxx.c), and updating r.series to >> use it. >> >> > Thanks for pointing me the way. I think it is a little beyond my skills, > but I might give it a try. > > Cheers, > Marcello. > > ___ > 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
Re: [GRASS-user] Problems with grass for OSX
This looks like you do not have the correct version of gdal installed. __ C. Michael Barton Director, Center for Social Dynamics & Complexity Professor of Anthropology, School of Human Evolution & Social Change Arizona State University Tempe, AZ 85287-2402 USA voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC) fax: 480-965-7671(SHESC), 480-727-0709 (CSDC) www:http://csdc.asu.edu, http://shesc.asu.edu http://www.public.asu.edu/~cmbarton On Dec 9, 2011, at 9:20 PM, grass-user-requ...@lists.osgeo.org wrote: > Date: Fri, 9 Dec 2011 17:04:12 -0200 > From: Jos? Carlos Guerrero Ant?nez > Subject: [GRASS-user] Problems with grass for OSX > To: grass-user@lists.osgeo.org > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > Hi, > > I can not boot properly grass and this is the message that appears in Terminal > window. > > OSX 10.6.8 > > Thanks in advance, > > > torre-b-024:~ joseca_IMac$ > '/Applications/GRASS-6.4.app/Contents/MacOS/grass.sh'; exit > Rebuilding Addon HTML manual pages index... > Rebuilding Addon menu... > Python 2.6.1 found. > > WELCOME TO GRASS Version 6.4.1 2011 > > 1) Have at your side all available GRASS tutorials > > 2) When working on your location, the following materials > are extremely useful: > - A topo map of your area > - Current catalog of available computer maps > > 3) Check the GRASS webpages for feedback mailinglists and more: > http://www.grass-gis.org > http://grass.osgeo.org > > Hit RETURN to continue > > Starting GRASS ... > Traceback (most recent call last): > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 641, in OnSelectLocation >self.listOfLocations[self.lblocations.GetSelection()])) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 590, in UpdateMapsets >self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 629, in GetListOfMapsets > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 651, in EncodeString > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 459, in getdefaultlocale >return _parse_localename(localename) > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 391, in _parse_localename >raise ValueError, 'unknown locale: %s' % localename > ValueError: unknown locale: UTF-8 > Traceback (most recent call last): > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 858, in >GRASSStartUp = StartUp(0) > File > "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", > line 7981, in __init__ > File > "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", > line 7555, in _BootstrapApp > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 843, in OnInit >StartUp.OnBrowse(None) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 703, in OnBrowse >self.OnSetDatabase(event) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 692, in OnSetDatabase >self.OnSelectLocation(None) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 641, in OnSelectLocation >self.listOfLocations[self.lblocations.GetSelection()])) > File > "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line > 590, in UpdateMapsets >self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 629, in GetListOfMapsets > File > "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", > line 651, in EncodeString > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 459, in getdefaultlocale >return _parse_localename(localename) > File > "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", > line 391, in _parse_localename >raise ValueError, 'unknown locale: %s' % localename > ValueError: unknown locale: UTF-8 > Error in GUI startup. If necessary, please > report this error to the GRASS developers. > Switching to text mode now. > Hit RETURN to continue... > > > > -- > > Dr. Jos? Carlos Guerrero > Laboratorio de Desarrollo Sustentable y Gesti?n Ambiental del Territorio > (LDSGAT) > Instituto de Ecolog?a y Ciencias Ambientales (IECA) > Facultad de Ciencias > Univ
Re: [GRASS-user] How to find the second largest value in a series of rasters.
Marcello: > > My problem is that I also need to find the *second largest* value and the > > corresponding raster number which contains the second largest value. > > > > I am doing that by iterating over all classes through a shell script, > but it > > obviously takes much more time than using a simple r.series command. > > > > Can anyone share any idea on how to accomplish that using r.series or > > something similar? > > Glynn: > There isn't an efficient way to do this using existing tools. If you > need the efficiency, I suggest adding a "second largest" aggregate to > lib/stats (based upon c_max.c and c_maxx.c), and updating r.series to > use it. > > Thanks for pointing me the way. I think it is a little beyond my skills, but I might give it a try. Cheers, Marcello. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to initialize raster maps.
Glynn: > Note that using the same map as both input and output results in > undefined behaviour. Even if it happens to work, there's no guarantee > that it will continue to work in future versions. > > If you need to implement an iterative algorithm, you should generate a > new map, then replace the original with g.remove+g.rename afterwards. > > Now, you scared me, because my ongoing work is pretty much based on this kind of code. I use r.mapcalc to create maps that serve as counters such as the command i++ in C. I initialize them as zero maps and then based on some logic I keep incrementing them, pretty much like the following pseudo-code to implement a sum: sum = 0 for i = 1 to 10 sum = sum + i end But I do more that that. For instance, if certain conditions are met in other maps at a given iteration and the initialized map is zero, then update the initialized map with the value of "i". If not, don't change it. In this way, my map is updated every iteration until no more cells meet the condition or some threshold is met and processing stops. I even use d.rast to display the map as it is being updated. How should I proceed to avoid the "undefined behaviour" and still get what I want? (I hope I was able to make myself understood :) ) Thanks for your comments. Glynn: > There is no point in referring to your original map in the r.mapcalc > expression; the result will always be created with the bounds and > resolution of the current region. AFAIK, the only modules which don't > behave this way are r.in.*, which create the output to match the size > of the input file. > > Yes, Moritz's suggestion showed me that. I perfectly understand it now, thanks. Glynn: I'm sure that there are other ways to create a map consisting of > zeros, but there's no reason to suspect that any of them will be > faster than r.mapcalc. > > There really isn't. I just suspected it was taking too long and found out in my little experiment that was only my fault, putting a logic expression where there was no need of it. Setting the region and r.mapcalc zero_map = 0 does the job. Thank you. Best, Marcello. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to find the second largest value in a series of rasters.
Marcello Gorini wrote: > I have to find the maximum value for each cell in a series of rasters and > generate both a maximum value map and a categorical map with the raster > number which contains the maximum value. > > R.series is, therefore, perfect to accomplish that using > method=maximum,max_raster. > > My problem is that I also need to find the *second largest* value and the > corresponding raster number which contains the second largest value. > > I am doing that by iterating over all classes through a shell script, but it > obviously takes much more time than using a simple r.series command. > > Can anyone share any idea on how to accomplish that using r.series or > something similar? There isn't an efficient way to do this using existing tools. If you need the efficiency, I suggest adding a "second largest" aggregate to lib/stats (based upon c_max.c and c_maxx.c), and updating r.series to use it. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to initialize raster maps.
Marcello Gorini wrote: > >> Inside a loop in a shell script, I keep incrementing the values of a > >> raster map using r.mapcalc with some logical expressions. Note that using the same map as both input and output results in undefined behaviour. Even if it happens to work, there's no guarantee that it will continue to work in future versions. If you need to implement an iterative algorithm, you should generate a new map, then replace the original with g.remove+g.rename afterwards. > >> Ideally, I would like to initialize it as a raster with the same > >> dimensions as my base map, but containing only zeros. > >> > >> In order to do that, I use the following expression: > >> > >> r.mapcalc "initialized_map=if(base_map>**99,1,0)" There is no point in referring to your original map in the r.mapcalc expression; the result will always be created with the bounds and resolution of the current region. AFAIK, the only modules which don't behave this way are r.in.*, which create the output to match the size of the input file. > But there is no way to avoid r.mapcalc, is there? I'm sure that there are other ways to create a map consisting of zeros, but there's no reason to suspect that any of them will be faster than r.mapcalc. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problems with grass for OSX
Hi, I can not boot properly grass and this is the message that appears in Terminal window. OSX 10.6.8 Thanks in advance, torre-b-024:~ joseca_IMac$ '/Applications/GRASS-6.4.app/Contents/MacOS/grass.sh'; exit Rebuilding Addon HTML manual pages index... Rebuilding Addon menu... Python 2.6.1 found. WELCOME TO GRASS Version 6.4.1 2011 1) Have at your side all available GRASS tutorials 2) When working on your location, the following materials are extremely useful: - A topo map of your area - Current catalog of available computer maps 3) Check the GRASS webpages for feedback mailinglists and more: http://www.grass-gis.org http://grass.osgeo.org Hit RETURN to continue Starting GRASS ... Traceback (most recent call last): File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 641, in OnSelectLocation self.listOfLocations[self.lblocations.GetSelection()])) File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 590, in UpdateMapsets self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) File "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", line 629, in GetListOfMapsets File "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", line 651, in EncodeString File "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", line 459, in getdefaultlocale return _parse_localename(localename) File "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", line 391, in _parse_localename raise ValueError, 'unknown locale: %s' % localename ValueError: unknown locale: UTF-8 Traceback (most recent call last): File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 858, in GRASSStartUp = StartUp(0) File "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", line 7981, in __init__ File "/Users/Shared/unix/wxpython-snow/lib/python2.6/site-packages/wx-2.8-mac-unicode/wx/_core.py", line 7555, in _BootstrapApp File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 843, in OnInit StartUp.OnBrowse(None) File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 703, in OnBrowse self.OnSetDatabase(event) File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 692, in OnSetDatabase self.OnSelectLocation(None) File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 641, in OnSelectLocation self.listOfLocations[self.lblocations.GetSelection()])) File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 590, in UpdateMapsets self.listOfMapsets = utils.GetListOfMapsets(self.gisdbase, location) File "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", line 629, in GetListOfMapsets File "/Users/Shared/src/GRASS/grass-6.4.1-snow/dist.i386-apple-darwin10.7.0/etc/wxpython/gui_modules/utils.py", line 651, in EncodeString File "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", line 459, in getdefaultlocale return _parse_localename(localename) File "/System/Library/Frameworks/Python.framework/Versions/2.6/lib/python2.6/locale.py", line 391, in _parse_localename raise ValueError, 'unknown locale: %s' % localename ValueError: unknown locale: UTF-8 Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... -- Dr. José Carlos Guerrero Laboratorio de Desarrollo Sustentable y Gestión Ambiental del Territorio (LDSGAT) Instituto de Ecología y Ciencias Ambientales (IECA) Facultad de Ciencias Universidad de la República Montevideo-URUGUAY ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.rast.stats (r.category) fail?
Hi Annalisa, On Fri, Dec 9, 2011 at 6:22 PM, Annalisa Minelli wrote: > Hi all, > I'm experiencing some problems with v.rast.stats and (I suppose) v.category. > I have a lines vector file and I want to obtain some statistics on each > feature. > I try to use the module but it soon exits with the message: > > v.rast.stats vector=ord1_pl_3d_cat raster=dem colprefix=stats percentile=90 > --verbose -c > ERROR: No categories found in raster map > > evenif I have some feature (with categories and "cat" column too on > associated dbf) in my file. could you please tell us - which operating system - which GRASS version > Then I tried r.category by itself and I discovered that sometimes it works > but sometimes not. > Align region on input raster map doesn't affect the result. > > Any hint? If possible, please send a North Carolina data set example for easy replication. cheers Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.rast.stats (r.category) fail?
Hi all, I'm experiencing some problems with v.rast.stats and (I suppose) v.category. I have a lines vector file and I want to obtain some statistics on each feature. I try to use the module but it soon exits with the message: v.rast.stats vector=ord1_pl_3d_cat raster=dem colprefix=stats percentile=90 --verbose -c ERROR: No categories found in raster map evenif I have some feature (with categories and "cat" column too on associated dbf) in my file. Then I tried r.category by itself and I discovered that sometimes it works but sometimes not. Align region on input raster map doesn't affect the result. Any hint? Best regards, Annalisa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to initialize raster maps.
Marcello: > > But there is no way to avoid r.mapcalc, is there? > > Hamish: > I doubt you'd get very much faster or more efficient even by > writing your own C module. (it would probably only take you a > minute to create a r.zero module out of the doc/raster/r.example > code; just change "return x;" to "return 0;" and compile. > > Well, out of curiosity, I made the changes to r.example and compiled it. I then made a small script that initializes a map with the dimensions of one of my maps ( 677 x 1114) inside a loop using the three approaches mentioned in this thread: 1 - r.mapcalc "initialized_map=if(base_map>99,1,0)" # my very smart method :) 2 - r.mapcalc initialized_map=0 # Moritz method 3 - r.example input=initialized_map output=zero_map # Hamish idea Since I initialize rasters 374 times in my original script, I made the loop in this way. Results: 1 - 49 seconds 2 - 19 seconds 3 - 58 seconds Bottom line: - Thanks Moritz for the tip, it will save me at least 30 seconds from now on. - Thanks Hamish, because I never knew r.example existed and it was fun to test the different approaches. But, no, r.example should stay just as that an example :) Thanks to all. Cheers, Marcello. P.S.: Hamish: > tip: to avoid all your scripts breaking when you want to use > grass7, always add a space around the '=' in r.mapcalc. > > r.mapcalc "initialized_map = 0" > > I will keep it in mind from now on. Tks. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to find the second largest value in a series of rasters.
Marcello: > I have to find the maximum value for each cell in a series > > of rasters and generate both a maximum value map and a > > categorical map with the raster number which contains the > > maximum value. > > Hamish: > see also: > > http://grass.osgeo.org/wiki/Time_series#Common_legends_for_many_raster_maps > > it sounds like you have done something similar already, but > maybe you could adapt it somehow, or it is faster. > > I think you'll have to use g.copy + r.null (or better r.reclass) > to hide the highest values. > > Thanks for the ideas. I still haven't quite figured out how to do it without a loop, but I will keep trying. Anyway, I found very useful this kind of expression: for MAP in `g.mlist rast pattern=$MAP_PATTERN` ; do eval `r.info -r $MAP` I didn't know it was possible. Thanks. Cheers, Marcello. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] errors in grass7 svn installation 49631
Hi, 2011/12/9 Pankaj Kr Sharma : > 33, in > from vdigit.wxdisplay import DisplayDriver, GetLastError > File "/usr/local/grass-7.0.svn/etc/gui/wxpython/vdigit/wxdisplay.py", line > 70, in > errtype = CFUNCTYPE(UNCHECKED(c_int), String, c_int) > NameError: name 'CFUNCTYPE' is not defined > GRASS 7.0.svn (oy3):~ > hopefully fixed in r49632. Martin -- Martin Landa * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] errors in grass7 svn installation 49631
Today, I updated grass7 and make resulted in a long list of modules having error during compilation. Even then , I ran make install as message during installation was that error may be fixed during run-time. GRASS7 installed and ran also, but with following message during startup. Is is OK or should I revert back to earlier version. rdcgis1@rdcgis1:/usr/svn/gis/grass7$ svn info Path: . Working Copy Root Path: /usr/svn/gis/grass7 URL: https://svn.osgeo.org/grass/grass/trunk Repository Root: https://svn.osgeo.org/grass Repository UUID: 15284696-431f-4ddb-bdfa-cd5b030d7da7 Revision: 49631 Node Kind: directory Schedule: normal Last Changed Author: martinl Last Changed Rev: 49629 Last Changed Date: 2011-12-09 02:04:16 +0530 (Fri, 09 Dec 2011) Message during startup: __ ___ _____ / / __ \/ | / ___/ ___/ / / _/ ___/ / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \ / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ / \/_/ |_/_/ |_/// \/___/// Welcome to GRASS 7.0.svn (2011) GRASS homepage: http://grass.osgeo.org This version running through:Bash Shell (/bin/bash) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui wxpython When ready to quit enter:exit Launching 'wxpython' GUI in the background, please wait... To run a command as administrator (user "root"), use "sudo ". See "man sudo_root" for details. GRASS 7.0.svn (oy3):~ > Traceback (most recent call last): File "/usr/local/grass-7.0.svn/etc/gui/wxpython/wxgui.py", line 34, in from lmgr.frame import GMFrame File "/usr/local/grass-7.0.svn/etc/gui/wxpython/lmgr/frame.py", line 44, in from lmgr.layertreeimport LayerTree File "/usr/local/grass-7.0.svn/etc/gui/wxpython/lmgr/layertree.py", line 35, in from mapdisp.frame import MapFrame File "/usr/local/grass-7.0.svn/etc/gui/wxpython/mapdisp/frame.py", line 38, in from vdigit.toolbarsimport VDigitToolbar File "/usr/local/grass-7.0.svn/etc/gui/wxpython/vdigit/toolbars.py", line 24, in from vdigit.mainimport VDigit File "/usr/local/grass-7.0.svn/etc/gui/wxpython/vdigit/main.py", line 18, in from vdigit.wxdigit import IVDigit, GV_LINES File "/usr/local/grass-7.0.svn/etc/gui/wxpython/vdigit/wxdigit.py", line 33, in from vdigit.wxdisplay import DisplayDriver, GetLastError File "/usr/local/grass-7.0.svn/etc/gui/wxpython/vdigit/wxdisplay.py", line 70, in errtype = CFUNCTYPE(UNCHECKED(c_int), String, c_int) NameError: name 'CFUNCTYPE' is not defined GRASS 7.0.svn (oy3):~ > Please help. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
