[GRASSLIST:2189] Re: How to Display USGS DEM Files?

Glynn Clements glynn.clements at virgin.net
Fri Jul 27 18:45:48 EDT 2001

Tom Russo wrote:

> As I recall, when I used r.in.gdal to do exactly what you're doing, I followed
> the r.in.gdal run with r.out.ascii, then used a sed script to replace the
> things like -32768 with 0.

r.mapcalc *might* be more efficient, e.g.:

	r.mapcalc 'result=if(source == -32768, null(), source)'

> Now, there may well be a way to specify how the nulls are handled by r.in.gdal

There isn't. With r.in.bin, you can specify a NULL value, but then you
have to specify the resolution and boundaries manually (r.in.bin
doesn't read .hdr files or similar).

Glynn Clements <glynn.clements at virgin.net>

More information about the grass-user mailing list