[GRASS-user] Setting Language
Greetings I'm doing some work in GRASS with a couple of friends from different countries (/ different languages) and our GRASS has some words translated in in our own languages but this can be quite difficult when, for example, I'm running something on a laptop of an italian :). So I want to define languages to English. I know that I must add a set language in etc/init.bat right? but which parameter shall i define? set LANGUAGE= set LANG= set LC_MESSAGES= set LC_ALL= And how do I call English Language? ENG? Thanks Katrin PS- Nice work on some translations (e.g. German and Spanish) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] snap point to line/break line at given points
Hello, I have a line-vector (river) and some points on that line. First I am straightening that line with v.generalize (method: boyle). Now, there is the problem that the points are not exactly on the line anymore. So I want to snap (move) the points to the new generalized line again. Which tool is able to do that? After that I want to split the line at these points and want to give the attribute value to the newly created segments. These segment values should represent the number/value of the upstream point (which was used to split). How can this be done in GRASS? Do you have any ideas to solve these two problems/challenges? thank you /johannes -- GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit gratis Handy-Flat! http://portal.gmx.net/de/go/dsl ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Extrac values with r.mapcalc
Hi Gaspar, r.mapcalc "new.map = if(flat_map==0.0234,flat_map,null())" -> if you want to preserve original values or r.mapcalc "new.map = if(flat_map==0.0234,1,null())" -> if you prefer to assign a value of 1 to those values equal to 0.0234 in the original raster map It should work. Lluís El dc 13 de 04 de 2011 a les 18:21 -0500, en/na Gaspar Reyes Póndigo va escriure: > Thanks Markus. > > I did test with r.mapcalc novo_raster = if"(mi flat_map==0.0234)", but > the map is not appropiate. I know the area of value (0.0234), but is > not is equal in the novo_raster. > > > El 13/04/2011, a las 17:27, Markus Neteler escribió: > > > > 2011/4/13 Gaspar Reyes Póndigo : > > > > > Hello, Dear grass_user. My first forum participation in the > > > grass_user list > > > email. > > > > > > Welcome! > > > > > > > How to extract values of map raster with r.mapcalc? I need make > > > new map > > > raster based on specific value of other raster. > > > > > > Please take a look at the if() conditions: > > > > http://grass.osgeo.org/grass64/manuals/html64_user/r.mapcalc.html > > > > Therein is also a section EXAMPLES. > > > > Markus > > > > ___ > 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] Error with GRASS 6.4.1 on CentOS 5
>>> if if ignoredDelimiter is None: >>> ignoredDelimiter = '' >> >> Thanks a lot, I will try to patch it. > > probably you will find more such statements in the code (see eg. > r45962). In that case, could you send a patch which I could apply in > SVN? I applied your fixes from yesterday (rev 45959 to 45964 from the 6.4 release branch) and the wxPython UI now starts properly (an updated version is now available in the ELGIS testing repository). I did various searches for other ternary operators but could not find any other. So maybe the issue will happen again when executing some actions, but at least the UI is now usable (I quickly opened a map set). Thanks again for your quick help! Mathieu PS: for reference, the description of this feature available in Python 2.5 but not in Python 2.4: http://www.python.org/dev/peps/pep-0308/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [El] Re: [GRASS-user] Error with GRASS 6.4.1 on CentOS 5
On Thu, Apr 14, 2011 at 1:02 PM, Mathieu Baudier wrote: ... > I applied your fixes from yesterday (rev 45959 to 45964 from the 6.4 > release branch) and the wxPython UI now starts properly (an updated > version is now available in the ELGIS testing repository). Thanks for your efforts, Mathieu, and the quick packaging! Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Setting Language
On Thu, Apr 14, 2011 at 9:27 AM, katrin eggert wrote: > Greetings > I'm doing some work in GRASS with a couple of friends from different > countries (/ different languages) and our GRASS has some words translated in > in our own languages but this can be quite difficult when, for example, I'm > running something on a laptop of an italian :). So I want to define > languages to English. I know that I must add a set language in etc/init.bat > right? > but which parameter shall i define? This is for English: LANG=en_US LANGUAGE=en_US LC_MESSAGES=en_US And this for Italian: LANG=it_IT LANGUAGE=it_IT LC_ALL=it_IT etc. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Raster addition
Hello. I am trying to add some raster maps. I thought I could do this in r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where the maps overlap and ignoring areas where only on of the maps have values. I would like to do something similar to "v.overlay operator=or". Anyone who knows how to do this? Martin Ytre-Eide___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster addition
On 14/04/11 13:54, Martin Album Ytre-Eide wrote: Hello. I am trying to add some raster maps. I thought I could do this in r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where the maps overlap and ignoring areas where only on of the maps have values. This is due to the way r.mapcalc handles NULL values. See the man page, section NULL support for an explanation. I would like to do something similar to "v.overlay operator=or". Anyone who knows how to do this? It depends on what you want the new map to contain as values. If you want to add the value of old_1 to old_2 in places where they overlap, one way would be to transform your null values to 0 using r.null null= and then applying your formula. If you just want a (any) value where one of the maps is non-null, then you can use if. Something like this (untested): new_map = if(!isnull(old_1) || !isnull(old_2), 1, null()) which should give you 1 where any of the two contains a value and a NULL value elsewhere. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
SV: [GRASS-user] Raster addition
I would like the NULL values that overlap with none NULL values to be treated as 0, and still have NULL values where both maps have NULL values. The combined map would be a large map with NULL values and added values. Maybe the best way is to use r.null and set NULL values to 0. Add the maps and then convert the 0 values to NULL values? Or could one set the value in the if statement to be a the sum of maps somehow? I guess my biggest problem is summing with NULL values - maybe there is a way to treat NULL values as 0? I'll have a look at the man page Thanks, Martin -Opprinnelig melding- Fra: Moritz Lennert [mailto:mlenn...@club.worldonline.be] Sendt: 14. april 2011 14:15 Til: Martin Album Ytre-Eide Kopi: grass-user@lists.osgeo.org Emne: Re: [GRASS-user] Raster addition On 14/04/11 13:54, Martin Album Ytre-Eide wrote: > Hello. > > I am trying to add some raster maps. I thought I could do this in > r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where > the maps overlap and ignoring areas where only on of the maps have > values. This is due to the way r.mapcalc handles NULL values. See the man page, section NULL support for an explanation. > I would like to do something similar to "v.overlay operator=or". > > Anyone who knows how to do this? It depends on what you want the new map to contain as values. If you want to add the value of old_1 to old_2 in places where they overlap, one way would be to transform your null values to 0 using r.null null= and then applying your formula. If you just want a (any) value where one of the maps is non-null, then you can use if. Something like this (untested): new_map = if(!isnull(old_1) || !isnull(old_2), 1, null()) which should give you 1 where any of the two contains a value and a NULL value elsewhere. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] assign heigght value to pixel(x,y)
Hi all, how to assign height value(z) to a particular pixel having x and y coordinates in an orthphoto using grass. Thanks in advance, nani ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: SV: [GRASS-user] Raster addition
On 14/04/11 14:31, Martin Album Ytre-Eide wrote: I would like the NULL values that overlap with none NULL values to be treated as 0, and still have NULL values where both maps have NULL values. The combined map would be a large map with NULL values and added values. Maybe the best way is to use r.null and set NULL values to 0. Add the maps and then convert the 0 values to NULL values? I think this is probably the easiest. Or could one set the value in the if statement to be a the sum of maps somehow? I guess my biggest problem is summing with NULL values - maybe there is a way to treat NULL values as 0? I'll have a look at the man page If you chose to do it all directly in r.mapcalc, then I guess the only option is nested if statements new_map = if(!isnull(old_1) && !isnull(old_2), old_1+old_2, if(!isnull(old_1) && isnull(old_2), old_1, etc) Moritz Thanks, Martin -Opprinnelig melding- Fra: Moritz Lennert [mailto:mlenn...@club.worldonline.be] Sendt: 14. april 2011 14:15 Til: Martin Album Ytre-Eide Kopi: grass-user@lists.osgeo.org Emne: Re: [GRASS-user] Raster addition On 14/04/11 13:54, Martin Album Ytre-Eide wrote: Hello. I am trying to add some raster maps. I thought I could do this in r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where the maps overlap and ignoring areas where only on of the maps have values. This is due to the way r.mapcalc handles NULL values. See the man page, section NULL support for an explanation. I would like to do something similar to "v.overlay operator=or". Anyone who knows how to do this? It depends on what you want the new map to contain as values. If you want to add the value of old_1 to old_2 in places where they overlap, one way would be to transform your null values to 0 using r.null null= and then applying your formula. If you just want a (any) value where one of the maps is non-null, then you can use if. Something like this (untested): new_map = if(!isnull(old_1) || !isnull(old_2), 1, null()) which should give you 1 where any of the two contains a value and a NULL value elsewhere. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: SV: [GRASS-user] Raster addition
Original-Nachricht > Datum: Thu, 14 Apr 2011 14:31:32 +0200 > Von: Martin Album Ytre-Eide > An: \'Moritz Lennert\' > CC: "grass-user@lists.osgeo.org" > Betreff: SV: [GRASS-user] Raster addition > I would like the NULL values that overlap with none NULL values to be > treated as 0, and still have NULL values where both maps have NULL values. The > combined map would be a large map with NULL values and added values. that is possible with the r.mapcalc and a if condition like: if(map A &&& map B, 0, null()) something like that (untested) might work. the null () provides NULL values where the condition is false. The triple &&& is the logic and (treats NULL values different from double &&). see: http://grass.fbk.eu/gdp/html_grass64/r.mapcalc.html for more information. /j > > Maybe the best way is to use r.null and set NULL values to 0. Add the maps > and then convert the 0 values to NULL values? > Or could one set the value in the if statement to be a the sum of maps > somehow? I guess my biggest problem is summing with NULL values - maybe there > is a way to treat NULL values as 0? I'll have a look at the man page > > Thanks, > > Martin > > -Opprinnelig melding- > Fra: Moritz Lennert [mailto:mlenn...@club.worldonline.be] > Sendt: 14. april 2011 14:15 > Til: Martin Album Ytre-Eide > Kopi: grass-user@lists.osgeo.org > Emne: Re: [GRASS-user] Raster addition > > On 14/04/11 13:54, Martin Album Ytre-Eide wrote: > > Hello. > > > > I am trying to add some raster maps. I thought I could do this in > > r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where > > the maps overlap and ignoring areas where only on of the maps have > > values. > > This is due to the way r.mapcalc handles NULL values. See the man page, > section NULL support for an explanation. > > > I would like to do something similar to "v.overlay operator=or". > > > > Anyone who knows how to do this? > > It depends on what you want the new map to contain as values. If you want > to add the value of old_1 to old_2 in places where they overlap, one way > would be to transform your null values to 0 using r.null null= and then > applying your formula. > > If you just want a (any) value where one of the maps is non-null, then you > can use if. Something like this (untested): > > new_map = if(!isnull(old_1) || !isnull(old_2), 1, null()) > > which should give you 1 where any of the two contains a value and a NULL > value elsewhere. > > Moritz > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user -- GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit gratis Handy-Flat! http://portal.gmx.net/de/go/dsl ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: SV: [GRASS-user] Raster addition
Martin Album Ytre-Eide wrote: > > I would like the NULL values that overlap with none NULL values to be > treated as 0, and still have NULL values where both maps have NULL values. > The combined map would be a large map with NULL values and added values. > hy, i guess r.patch can do that for you ... http://grass.osgeo.org/grass70/manuals/html70_user/r.patch.html r.patch ... hope that i´m not wrong cheers, tim -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Raster-addition-tp6272457p6272716.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Adding labels to Classes
Greetings I have a Land Cover map and I want to substitute Land Cover codes by numbers. Besides r.reclass (1 = 1 Agriculture ; 2 = 2 Forest 3 = 3 Water) what other method does exists in GRASS? Thanks ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
SV: SV: [GRASS-user] Raster addition -Solved
Don't know who of you led me in the right direction, but eventually the line: "C=A + if(isnull(B),0,B)" in the manual solved my issue. New_map = if(isnull(old_1),0,old_1) + if(isnull(old_2,0,old_2)) + etc and a following new_map=if(new_map>0,new_map,null()) I had four maps in my case, and that is why "new_map = if(!isnull(old_1) && !isnull(old_2), old_1+old_2, if(!isnull(old_1) && isnull(old_2), old_1, etc)" would become to long (4!=24 statements?) Thank you both Johannes and Moritz. Martin -Opprinnelig melding- Fra: Johannes Radinger [mailto:jradin...@gmx.at] Sendt: 14. april 2011 15:13 Til: Martin Album Ytre-Eide; mlenn...@club.worldonline.be Kopi: grass-user@lists.osgeo.org Emne: Re: SV: [GRASS-user] Raster addition Original-Nachricht > Datum: Thu, 14 Apr 2011 14:31:32 +0200 > Von: Martin Album Ytre-Eide > An: \'Moritz Lennert\' > CC: "grass-user@lists.osgeo.org" > Betreff: SV: [GRASS-user] Raster addition > I would like the NULL values that overlap with none NULL values to be > treated as 0, and still have NULL values where both maps have NULL > values. The combined map would be a large map with NULL values and added > values. that is possible with the r.mapcalc and a if condition like: if(map A &&& map B, 0, null()) something like that (untested) might work. the null () provides NULL values where the condition is false. The triple &&& is the logic and (treats NULL values different from double &&). see: http://grass.fbk.eu/gdp/html_grass64/r.mapcalc.html for more information. /j > > Maybe the best way is to use r.null and set NULL values to 0. Add the > maps and then convert the 0 values to NULL values? > Or could one set the value in the if statement to be a the sum of maps > somehow? I guess my biggest problem is summing with NULL values - > maybe there is a way to treat NULL values as 0? I'll have a look at > the man page > > Thanks, > > Martin > > -Opprinnelig melding- > Fra: Moritz Lennert [mailto:mlenn...@club.worldonline.be] > Sendt: 14. april 2011 14:15 > Til: Martin Album Ytre-Eide > Kopi: grass-user@lists.osgeo.org > Emne: Re: [GRASS-user] Raster addition > > On 14/04/11 13:54, Martin Album Ytre-Eide wrote: > > Hello. > > > > I am trying to add some raster maps. I thought I could do this in > > r.mapcalc, but "new_map= old_1 + old_2" seems to produce a map where > > the maps overlap and ignoring areas where only on of the maps have > > values. > > This is due to the way r.mapcalc handles NULL values. See the man > page, section NULL support for an explanation. > > > I would like to do something similar to "v.overlay operator=or". > > > > Anyone who knows how to do this? > > It depends on what you want the new map to contain as values. If you > want to add the value of old_1 to old_2 in places where they overlap, > one way would be to transform your null values to 0 using r.null null= > and then applying your formula. > > If you just want a (any) value where one of the maps is non-null, then > you can use if. Something like this (untested): > > new_map = if(!isnull(old_1) || !isnull(old_2), 1, null()) > > which should give you 1 where any of the two contains a value and a > NULL value elsewhere. > > Moritz > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user -- GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit gratis Handy-Flat! http://portal.gmx.net/de/go/dsl ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Extrac values with r.mapcalc
Hi! Gaspar, you can do it one way as follows. 1) create a vector of desired area that you would like to extract values for. 2) use v.to rast for converting it into a raster use value = 1 for that raster 3)open mapcalc and put input raster as A and new created raster as B 4) enter formula if(B == 1), A, 0), also name the new raster 5) the new raster shall have the values from the input raster I hope I got it right for you. rgs, Ambrish Dhaka On 4/14/2011 3:16 AM, Gaspar Reyes Póndigo wrote: Hello, Dear grass_user. /My first/ forum /participation/ in the grass_user list email. How to extract values of map raster with r.mapcalc? I need make new map raster based on specific value of other raster. ___ 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] Setting Language
Thanks Markus Just to confirm: I have to place that in etc/Init.bat right? It doesn't matter where? Thanks 2011/4/14 Markus Neteler > On Thu, Apr 14, 2011 at 9:27 AM, katrin eggert > wrote: > > Greetings > > I'm doing some work in GRASS with a couple of friends from different > > countries (/ different languages) and our GRASS has some words translated > in > > in our own languages but this can be quite difficult when, for example, > I'm > > running something on a laptop of an italian :). So I want to define > > languages to English. I know that I must add a set language in > etc/init.bat > > right? > > but which parameter shall i define? > > This is for English: > > LANG=en_US > LANGUAGE=en_US > LC_MESSAGES=en_US > > And this for Italian: > LANG=it_IT > LANGUAGE=it_IT > LC_ALL=it_IT > > etc. > > Markus > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while running v.db.addcol
Greetings I0'm using winGRASS. When I run v.db.addcol in my WinXP I get this error: Default driver / database set to: driver: dbf database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ 0 [main] sh 5240 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 5436 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump c:/grass6/scripts/v.db.addcol: [: : integer expression expected 0 [main] sh 2264 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3416 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump ERROR: Column not found But when I run in Win7 or other Win Machine (XP also) i do not get this error. Any explanation for this weird incident? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error while running v.db.addcol
On Thu, Apr 14, 2011 at 6:13 PM, Jenny Turner wrote: > Greetings > I0'm using winGRASS. Please always indicate the version... > When I run v.db.addcol in my WinXP I get this error: > Default driver / database set to: > driver: dbf > database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ > 0 [main] sh 5240 open_stackdumpfile: Dumping stack > trace to sh.exe.stackdump > 0 [main] sh 5436 open_stackdumpfile: Dumping stack > trace to sh.exe.stackdump > c:/grass6/scripts/v.db.addcol: [: : integer expression > expected -> problem > 0 [main] sh 2264 open_stackdumpfile: Dumping stack > trace to sh.exe.stackdump > 0 [main] sh 3416 open_stackdumpfile: Dumping stack > trace to sh.exe.stackdump > ERROR: Column not found > But when I run in Win7 or other Win Machine (XP also) i do not get this > error. Any explanation for this weird incident? Some fixes have been applied recently http://trac.osgeo.org/grass/log/grass/branches/releasebranch_6_4/scripts/v.db.addcol which are available now in GRASS 6.4.1. If you can, please update to this bugfix release. But maybe others have an idea, too. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Setting Language
On Thu, Apr 14, 2011 at 6:12 PM, katrin eggert wrote: > Thanks Markus > Just to confirm: I have to place that in etc/Init.bat right? It doesn't > matter where? > Thanks Well it may matter... Maybe put it in the beginning above the line: set GIS_LOCK=1 Please let us know how it goes. I hope that we can add a GUI switch at the Prague Community Sprint. Markus > 2011/4/14 Markus Neteler >> >> On Thu, Apr 14, 2011 at 9:27 AM, katrin eggert >> wrote: >> > Greetings >> > I'm doing some work in GRASS with a couple of friends from different >> > countries (/ different languages) and our GRASS has some words >> > translated in >> > in our own languages but this can be quite difficult when, for example, >> > I'm >> > running something on a laptop of an italian :). So I want to define >> > languages to English. I know that I must add a set language in >> > etc/init.bat >> > right? >> > but which parameter shall i define? >> >> This is for English: >> >> LANG=en_US >> LANGUAGE=en_US >> LC_MESSAGES=en_US >> >> And this for Italian: >> LANG=it_IT >> LANGUAGE=it_IT >> LC_ALL=it_IT >> >> etc. >> >> Markus > > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user