[GRASS-user] r.mapcalc-problem

GRASS GRASS user
1

hi,

i test a little bit the examples in

http://skagit.meas.ncsu.edu/~helena/grasswork/grassbookdat07/ncexternal/grasstestrast.txt

on grass64 from osgeo4w on a Win Vista 32 box.

I copy simply the command from the file above in the command line of the wx-gui

first example:

# find developed areas close to lakes
r.buffer lakes out=lakes_buff dist=60,120,240,500 --o
d.rast lakes_buff
r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
d.rast developed_lake
d.vect streets_wake co=grey
d.rast -o lakes

following error message for the r.mapcalc-command:

(Thu Mar 12 22:32:47 2009)
r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
syntax error, unexpected $end, expecting '='
Parse error
(Thu Mar 12 22:32:48 2009) Command finished (0 sec)

second example:

r.mapcalc "streets_travtime=1./streets_speed"
r.info streets_travtime

following error message for the r.mapcalc-command:

(Thu Mar 12 22:34:39 2009)
r.mapcalc "streets_travtime=1./streets_speed"
syntax error, unexpected $end, expecting '='
Parse error
(Thu Mar 12 22:34:39 2009) Command finished (0 sec)

best regards
helli
_______________________________________________________________________
DSL zum Nulltarif + 20 Euro Extraprämie bei Online-Bestellung über die
DSL Freundschaftswerbung! http://dsl.web.de/?ac=OM.AD.AD008K15279B7069a

2

On Thu, Mar 12, 2009 at 10:38 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

hi,

i test a little bit the examples in

http://skagit.meas.ncsu.edu/~helena/grasswork/grassbookdat07/ncexternal/grasstestrast.txt

on grass64 from osgeo4w on a Win Vista 32 box.

I copy simply the command from the file above in the command line of the wx-gui

first example:

# find developed areas close to lakes
r.buffer lakes out=lakes_buff dist=60,120,240,500 --o
d.rast lakes_buff
r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
d.rast developed_lake
d.vect streets_wake co=grey
d.rast -o lakes

following error message for the r.mapcalc-command:

(Thu Mar 12 22:32:47 2009)
r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
syntax error, unexpected $end, expecting '='
Parse error
(Thu Mar 12 22:32:48 2009) Command finished (0 sec)

I get the same on Linux. Apparently something in the GUI
is eating some characters of the formula.

I have opened a ticket for this:
http://trac.osgeo.org/grass/ticket/528

Markus

second example:

r.mapcalc "streets_travtime=1./streets_speed"
r.info streets_travtime

following error message for the r.mapcalc-command:

(Thu Mar 12 22:34:39 2009)
r.mapcalc "streets_travtime=1./streets_speed"
syntax error, unexpected $end, expecting '='
Parse error
(Thu Mar 12 22:34:39 2009) Command finished (0 sec)

best regards
helli
_______________________________________________________________________
DSL zum Nulltarif + 20 Euro Extraprämie bei Online-Bestellung über die
DSL Freundschaftswerbung! http://dsl.web.de/?ac=OM.AD.AD008K15279B7069a

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

3

Hi,

2009/3/13 Markus Neteler <neteler@osgeo.org>:

r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
syntax error, unexpected $end, expecting '='
Parse error

try without quotation marks, e.g.

r.mapcalc developed_lake = if(landuse96_28m==1 || landuse96_28m==2,
lakes_buff, null())

Martin

--
Martin Landa <landa.martin gmail.com> * http://gama.fsv.cvut.cz/~landa

4

On Fri, Mar 13, 2009 at 9:00 AM, Martin Landa <landa.martin@gmail.com> wrote:

Hi,

2009/3/13 Markus Neteler <neteler@osgeo.org>:

r.mapcalc "developed_lake=if(landuse96_28m==1 || landuse96_28m==2, lakes_buff, null())"
syntax error, unexpected $end, expecting '='
Parse error

try without quotation marks, e.g.

r.mapcalc developed_lake = if(landuse96_28m==1 || landuse96_28m==2,
lakes_buff, null())

That works. Do you see a chance to also accept it quoted
as required for the true cmd line?

BTW: a subsequent
d.rast developed_lake
adds an entry to the manager but in the form
d.rast flags=developed_lake

which should be
d.rast map=developed_lake

Markus