Eric Goddard wrote:
> Hi, I need to do the same thing to some WV-2 imagery. I'm using the command
>
> grass.mapcalc("$output = if($input_rast>1.0, 1.0, $input_rast)",
> output=wv2_out, input_rast=wv2_in)
The expression seems fine -- in general. I also see in my *bash* history
entries like
--%<---
r.mapcalc "TopoCorr.B.Trimmed.DOS1.ToCR.1 = if( TopoCorr.B.Trimmed.DOS1.ToCR.1
> 1.0, 1.0, TopoCorr.B.Trimmed.DOS1.ToCR.1 )"
or
for MAP in `g.mlist rast pat=TopoCorr.B.Trimmed.DOS1.ToCR.?`; do r.mapcalc
"${MAP} = if( ${MAP} > 1.0, 1.0, ${MAP} )"; done
--->%--
Back to yours:
# in python, you instruct
grass.mapcalc(
# and put code inside
"$output = if($input_rast>1.0, 1.0, $input_rast
# then close it
)"
Maybe it should be (not tested :-p):
grass.mapcalc( "$output = if( $input_rast > 1.0, 1.0, $input_rast)" )
?
Nikos
[rest deleted]
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
