[GRASS-dev] On i.histo.match (Re: On (Landsat) imagery naming patterns)

Nikos Alexandris nik at nikosalexandris.net
Thu Feb 21 01:18:59 PST 2013


MM:
> > >> For rescaling to 0-1000, you could try something like
> > >> r.mapcalc "bandX_int = round(1000 * bandX)"

NA:
> > > Nice!  But I wonder why round-ing is required?

MM:
> > rounding is not required, just a bit more accurate than truncating
> > with int(1000 * bandX).

Just FYI, results look nice!  I even convert back to 0-1.0 via 

r.mapcalc --o "${HistoMatchedMap} = ${HistoMatchedMap} / 1000.0"

Thnx, Nikos

ps- I wonder if *10000 is *better* for higher precision?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 230 bytes
Desc: This is a digitally signed message part.
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20130221/a8864e8f/attachment.pgp>


More information about the grass-dev mailing list