dicomwrite creates a black and white image, why?

Hi there,

I have a problem in writing and reading a dicom CT image , I explained in the code below, I think I miss something, would you please let me know what I miss. Thanks a lot.

x=dicomread('t3.dcm'); % a grayscale image

metadata=dicominfo('t3.dcm');

... process x data

imshow(x,[]); % a correct grayscale image

dicomwrite (x,'t1.dcm',metadata);

x1=dicomread('t1.dcm');

imshow(x1,[]); % it is a black and white image, why?

Answer 1

[venkat vasu]

info = dicominfo('CT-MONO2-16-ankle.dcm'); I = dicomread(info); dicomwrite(I,'ankle.dcm',info);

this code will help you...

Answer 2

[Sean de Wolski]

What is the class of x after you read it in?

class(x) %?

And after this:

max(x(:))
min(x(:))

I would guess you are writing the file out and not accounting for the class of the data, causing it to scale wierdly or saturate at the edges.

[MYimage]

thank you for your answer, I corrected it by converting the data from double to uint16.