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

joanna mardas joanna.mardas at wp.pl
Sat May 9 17:28:02 PDT 2015


Hey!
> joanna, once again, the easy "other way" is posted in my first reply, I
> think.  You just need to multiply with 1000, perform the histogram
> matching, then divide by 1000.0 to get back to floats.

So I did that: int(oldmap*1000), then i.histo.match (band to band). I did also r.recode (oldlow:oldhigh:0:255), and then i.histo.match (band to band).  Histograms in the middle (results of int(oldmap*1000) and i.hist.match) look strange and the images have "holes"-places with no data, totally transparent. After r.recode  there are no places with no data or there are 1 or 2 very small (max 5 pixels). I don't know why it is like that. So I will either do the classification after r.recode and 
i.histo.match or without r.recode and i.histo.match.
BTW, does classification work with both float and integer?
Joanna
-------------- next part --------------
A non-text attachment was scrubbed...
Name: histograms.jpg
Type: image/jpeg
Size: 309471 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20150510/fdf91c88/attachment-0001.jpg>


More information about the grass-dev mailing list