[GRASS-user] Filtering high "outliers" in Landsat reflectance imagery?

Glynn Clements glynn at gclements.plus.com
Thu Aug 1 19:21:39 PDT 2013


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)
> 
> from python (input and output names
> changed for simplicity) but instead of replacing the values greater
> than 1 with 1, it replaces them with NaNs.  Is there something wrong
> with my mapcalc expession?

No.

Is there anything unusual about the actual map names? Names containing
any of the characters

	^ # ( ) [ ] + - % > < ! & | ? : ; ~

need to quoted.

-- 
Glynn Clements <glynn at gclements.plus.com>


More information about the grass-user mailing list