Dicom file not displaying properly.

I am trying to read and display dicom(.dcm) file in matlab and save it later. However, the file is not displaying properly. Please find my code below:

info=dicominfo("IMG-0012-00218.dcm");
I=dicomread(info);

figure, imshow(I,[])
h = drawfreehand; %draw something 
M = ~h.createMask();
I(M) = 0;
imshow(I,[]);
imwrite(I,'Image1.png','mode','loseless')

I am unable to attach the .dcm file here as it is not allowing me to do it. I am getting the image below which is not correct. I want the image to be displayed in the way defined in 'IMG-0012-00218.png'.

Answer 1

[Harsha]

Hi Warid Islam

To get the output as shown in IMG-0012-00218.png you have to change the window level in the raw data.

Refer to this MATLAB answer by Simon Chan to resolve your issue.

https://www.mathworks.com/matlabcentral/answers/1435937-dicomwrite-to-save-dicom-image-with-window-level-setting

[Warid Islam]

Hi @Harsha,

Thank you for your response. I have a query. Your code solves my above question. I want to do further operation on that image. I want to manually segment the image and save the resultant image in the .png format. I tried the code below. However, I didn't get the desired result.The saved image looks like below.

I=dicomread('IMG-0012-00219.dcm');
metadata = dicominfo('IMG-0012-00219.dcm');
WC = metadata.WindowCenter;                                         % original [40; 40]
WW = metadata.WindowWidth; 
figure, imshow(I,[WC(1)-WW(1), WC(1)+WW(1)]); 

h = drawfreehand; %draw something 
M = ~h.createMask();
I(M) = 0;
imshow(I,[WC(1)-WW(1), WC(1)+WW(1)]);
imwrite(im2double(I),'Image3.png','mode','loseless')