如何设置此医学 DICOM 图像的显示范围

How to set the display range for this medical DICOM images

这是我用来显示 DICOM 图像的代码,当我给出完整的显示范围时,它显示像左图那样的毛刺图像,而当我增加较低的显示范围时,图像看起来更清晰。

[img, map] = dicomread('D:\Work 2017\Mercy\trial\trial4\Export0000\MG0001.dcm');
info = dicominfo('D:\Work 2017\Mercy\trial\trial4\Export0000\MG0001.dcm' );
mini = min(img(:));
maxi = max(img(:));

figure,
subplot(131), imshow(img, [mini maxi]); title('Full display range')
subplot(132), imshow(img, [maxi*0.7 maxi]); title('70% and above display range')
subplot(133), imshow(img, [maxi*0.8 maxi]); title('80% and above display range')

我希望始终看到类似于右侧图像的图像,而不给出我在上面代码中使用的显示范围

通常 DICOM 将具有 WindowCenter and WindowWidth 标签,用于指定推荐的 window/level 设置。您可以通过以下方式将这些转换为颜色限制

% Get the DICOM header which contains the WindowCenter and WindowWidth tags
dcm = dicominfo(filename);

% Compute the lower and upper ranges for display
lims = [dcm.WindowCenter - (dcm.WindowWidth / 2), ...
        dcm.WindowCenter + (dcm.WindowWidth / 2)];

% Load in the actual image data
img = dicomread(dcm);

% Display with the limits computed above
imshow(img, lims);

或更简短

lims = dcm.WindowCenter + [-0.5 0.5] * dcm.WindowWidth;

如果这些值不可接受,那么最好提供 user-adjustable window/level(例如 imtool 中的对比工具),因为不太可能有任何方法可靠地获得 "acceptable" 对比度,因为它是主观的。