DICOM, pixel value

how to proceed to find specific pixel value range (suppose 1000 to 1500) in a 512x512 16 uint DICOM image and color this range with yellow or any color ?

[KALYAN ACHARJYA]

Have you tried with dicom read?

[Rik]

Note that in the case of CT you need to account for the RescaleSlope and RescaleIntercept DICOM attributes to scale the data to HU. If you want to do that, you will need to use either the double data type, or int16 to account for negative values.

Answer 1

[Walter Roberson]

img = dicomread(FileName);
mask = 1000 <= img & img <= 1500;
R = img; R(mask) = 0;                   %yellow is 0, max, max
G = img; G(mask) = intmax(class(img));
B = img; B(mask) = intmax(class(img));
RGB = cat(3, R, G, B);
image(RGB)

[Deepa S]

thankyou, it works.

Answer 2

[Selva Karna]

clc
clear all;
close all;
your_dcm=dcmread('dcm/root');
ms=your_dcm>=a&your_dcm>=b;
rl=your_dcm;
fhol=find(ms==1);
rl(fhol)=0;

RGB_IMAGE(:,:,1)=rl;
RGB(:,:,2)=your_dcm;
RGB(:,:,3)=your_dcm;