How to use The Minimum Error Thresholding method for CT images using Matlab?

I want to implement the method of "Minimum Error Threshold - By J. Kittler and J. Illingworth" for CT images to finding the best threshold for segmentation in Matlab. So, I found the following code that is written using sundus. However, in CT images typical Hounsfield Unit values are between -1000 : >1000.

So, How Can I set the range of the threshold for CT images using that function and is there any code of "Kittler" for CT images?

Answer 1

[Matt J]

You could pre-map your input image to have the intensity range expected by your code using imadjust,

    J = imadjust(I,[low_in high_in],[low_out high_out])