Dicom file not displaying properly.

5 次查看（过去 30 天）
Warid Islam
Warid Islam 2022-11-15
评论： Warid Islam 2022-11-18
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'.

请先登录，再回答此问题。

回答（1 个）

Harsha
Harsha 2022-11-18
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
  1 个评论
Warid Islam
Warid Islam 2022-11-18
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')

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Convert Image Type 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by