On Fri, Apr 3, 2009 at 9:07 AM, Roger Bivand <[email protected]> wrote: > On Fri, 3 Apr 2009, Jarek Jasiewicz wrote: > >> Hi >> >> I during batch processing of couple of maps I found strange behaviour of >> readRAST6 when I import cell map with value 999. >> Instead of 999 Na is imported (?? _ >> >> It is for me completely strange, so I test it with folowing: > > OK. This is a consequence of trying to guess the NODATA value to pass to > r.out.gdal,
The C version of r.out.gdal has a nodata parameter. > and I don't see it in 0.5-19 (my local slightly changed from > CRAN latest 0.5-18). For certain range values retrieved from r.info, NODATA > is assigned the value 999 (R/bin_link.R, about line 113). Please try 0.5-18, > I suspect that the condition was changed. Alternatively, set NODATA to a > suitable value outside the range of the data. Note that I don't have a 6.5 > GRASS, and if r.info output has changed format, the parsing will be wrong, It didn't change the format to my knowledge. BUT: To get the range, use r.info -r test_map min=999 max=999 The parsing of plain r.info will fail in non-English locales (or enforce C locale). ... > I guess that there is a smarter way of getting the raster min/max too, isn't > there? Yes, please always use the script output flags which are not subject to locale and way easier to parse. If missing in a module, we are happy to add it. It also exists in GRASS 6.3: http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_3/raster/r.info/main.c#L74 Best Markus _______________________________________________ grass-stats mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-stats
