Limiting Intensity Values of an Image

Hi,

I have about 200 CT scans that I need to limit their intensity values between -2048 and 2048. I tried histeq and imadjust but they did not work. I get the best result when I use imshow(image, [-2048,2048]). But I need to save the result data of this imshow.

Original image-> http://img850.imageshack.us/img850/268/39979786.jpg

Imshow result -> http://img820.imageshack.us/img820/3909/76610719.jpg

What can I do?

Best

Answer 1

[Walter Roberson]

See the FEX contribution "freezeColors"

Answer 2

[Image Analyst]

Of course histeq() is not what you want. But I don't know why imadjust() didn't work - you must have called it incorrectly. Alternatively you could do

    scaledImage = uint8(255 * (single(yourImage) + 2048)/4096);
    imwrite(scaledImage, fullFileName);

That will scale your image between -2048 and +2048 and look like what imshow() does.