I've been trying to write a script creating binary maps:
a=raw_input ("enter a threshold")
b=float(a)
grass.run_command('r.mapcalculator', formula="**'base_map'>=b**", outfile='map'+str(a), overwrite=True)
When I'm trying to execute it there's an error message :
Invalid map
Parse error
Obviously, when the relevant part of the code is changed (e.g. formula="'base_map'>=**0.7**"
) everything is working fine. However, my goal was to write a script creating a map automatically when the "a" (the threshold) is provided.
I would be grateful for any hints.
Cheers