[GRASS-dev] Landsat 5TM pre-processing - histogram matching - problem

Nikos Alexandris nik at nikosalexandris.net
Wed May 6 03:35:21 PDT 2015


joanna mardas:

> 
> > I was trying to convert to integer through *r.mapcalc* (the only other
> > way I've found) with the function *int(x)* and x was my map (I hope it
> > was ok to put map instead of x) but the result was a map with "one
> > shade" of grey, so maybe it was wrong to put a map instead of x...

Moritz:

> int() rounds, so if you maps varies between 0 and 1, rounding will make 
> you lose almost all information. You have to multiply your values by a 
> multiple of ten corresponding to the degree of precision you want before 
> converting it to int. I.e. for a value of 0.12345678, if you want to 
> keep the first five numbers after the decimal point, you have to do 
> something like this:
> 
> newMap = int(oldMap * 100000)


I wrote on "Mon, 4 May 2015 22:14:41 +0300" how I approached this in the
past.  will try to support, if I can, for the rest of questions as well.

Thanks, Nikos

[snip]


More information about the grass-dev mailing list