Load DICOM image, binarize, and save a new DICOM image

3 次查看（过去 30 天）

显示更早的评论

georgina company 2024-5-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2113826-load-dicom-image-binarize-and-save-a-new-dicom-image

评论： georgina company 2024-5-2

采纳的回答： Karl

在 MATLAB Online 中打开

I load a DICOM Image and binarize it.

Then I save a new dicom image. However, the saved image is only black, and not black and white.

Here is the code:

for i = 1:length(d)
    na=d(i).name;
    TF = startsWith(na,'Tumour');
    if TF == 1
        info = dicominfo(na);
        greyImage = dicomread(info);
        BW=imbinarize(greyImage);
        fname1 = strcat('bin_',na)
        dicomwrite(BW, fname1,info);
    else
        if TF == 0
            fprintf('File do not match')
        end
    end
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

DGM 2024-5-1

编辑：DGM 2024-5-1

I can run this code on a carefully-selected DICOM file, and the binarized result can be read and viewed fine. Your written image is black. Without any information about your image file and how you're viewing it, we cannot know whether this is because:

something is going wrong with how you're viewing it
something is going wrong with how you're saving it
or something is going wrong with how you binarized it to begin with.

I'm going to guess that it's the latter; after all, imbinarize() is not magic. Unless some effort were taken, it might indeed be expected to give garbage results. Did you look at the range of your image data? A lot of DICOM images only use a tiny portion of their available dynamic range. If present, things like annotations, borders, or other decorations may present as features which lie far outside the actual dynamic range of the image data. For example, I have images where:

the available dynamic range is [-32768 32767]
the actual image data spans only [-59 101]
the borders and annotations span [-2048 0]

What do you think will happen when I blindly binarize that?

Again, this is just a guess.

Walter Roberson 2024-5-1

在 MATLAB Online 中打开

dicomwrite(BW, fname1,info);

That asks to write a dicom image using the attributes of the original image. But the attributes of the original image are potentially unsuitable for writing an image of type logical.

georgina company 2024-5-2

Der both,

I guess it is what Sir Roberson suggests. However, if I save the binarized image without the info when I construct the MedicalVolume it says that PatientPosition is needed.

Then, how can I construct the medical Volume?

请先登录，再进行评论。

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