On Thu, Aug 1, 2013 at 9:21 PM, Glynn Clements <gl...@gclements.plus.com> wrote: > > 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.
No special characters. An example image name is x2NOV04172233_M3DS_R6C2_052823926030_01_P001_TOAR.1. The mapcalc evaluates successfully (and correctly!) on my machine at home. My work machine is using the following build: version=7.0.svn date=2013 revision=56943 build_date=2013-08-01 ./configure --prefix=/usr/local/grass70 --with-postgres=yes --with-postgres-includes=/usr/include/postgresql --enable-largefile --with-sqlite --with-freetype=yes --with-freetype-includes=/usr/include/freetype2 --with-proj-share=/usr/share/proj --with-python --with-cxx --with-wxwidgets --with-tcltk-includes=/usr/include/tcl8.5 --with-odbc=yes --enable-64bit --with-geos=yes --with-openmp=yes --with-liblas=yes --with-cairo=yes --with-opencl=yes --with-pthread I'll have to bost my home machine's version later today. I compiled the latest svn yesterday to see if that would resolve the issue, but I still get null values. Not sure which version I was using before. Any recommendations on how I may be able to track the error down? I'm using Ubuntu 13.04 64-bit. I wonder if my issue could be linked to http://osgeo-org.1560.x6.nabble.com/seemingly-random-errors-in-r-mapcalc-td5070327.html ? Unfortunately my mapcalc operation doesn't actually raise an error so I don't think Soren's recommendation would help in my case. _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user