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?

Attachment: signature.asc
Description: This is a digitally signed message part.

_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev

Reply via email to