Translating Multi-Otsu Threshold from java

[Alexandre Fioravante de Siqueira]

Hi everyone,

I was thinking about using multi-Otsu thresholding for some research. I am used to use the ImageJ available plugin, http://imagej.net/Multi_Otsu_Threshold.

Could you give me a hand on writing such a function which could do this for python? Maybe putting it into scikit-image. Or else, could you point me a better alternative?

Thank you very much!

Alex

[Josh Warner]

We do have Otsu, and I think it would be within the scope of scikit-image. 

I think I'd find use for this, personally, and depending on the implementation you may be able to reuse code from simple binary Otsu internally.

Seems reasonable - PR would be welcome!

Josh

[Alexandre Fioravante de Siqueira]

Hi Josh,

OK, I'm on it. I'll try to implement from the paper.

Alex

[Alexandre Fioravante de Siqueira]

Hi all,

the preliminary code is ready, but it could be improved a lot. For instance, there's a function named maxsigma() which has a lot of nested for loops.

Here is the code:

https://gist.github.com/alexandrejaguar/fd767149b3a667f8be125734893f020d
Could you give me a hand to improve that?

Thank you very much!

Kind regards,

Alex

[Alexandre Fioravante de Siqueira]

Done (for now): https://github.com/scikit-image/scikit-image/pull/2076.

[Alexandre Fioravante de Siqueira]

Hi all,

I am having a bad time trying to correct the first version of this code.

The values from my function and Tosa's does not match. In my tests, my maxSigma is higher than his; I believe I am having a problem with the values returned from the histogram.

Could you give me a hand on that? This is the code so far:

https://github.com/alexandrejaguar/software/blob/master/threshold_multiotsu.py

Thanks for your time!

Kind regards,

Alex