How to multiply all the pixel number of each slice with their own rescale slope(dicominfo) in each slice

Dear all,

I have the file as attached. the format is dicom (.dcm)

How to multiply all the pixel number of each slice with their own rescale slope(in dicominfo) in each slice?

**that means each slice have their own rescale slope value.

so that the new pixel value will come out when I view this images.

[Simon Chan]

I am afraid there is only one rescale slope value for all the images belongs to the same DICOM file.

[mohd akmal masud]

Dear @Simon Chan,

Yah, I want to check the rescale slope for each slice, is it same or different.

I used this command

% Let say want to know the dicominfo for slice number 49, 50
S49 = dicominfo('jaszakspect.dcm', 'Frames',49);
S50 = dicominfo('jaszakspect.dcm', 'Frames',50);

But I got error

Error using dicominfo>parseInputs

Unrecognized parameter name "frames".

Error in dicominfo (line 13)

args = parseInputs(varargin{:});

[Rik]

What leads you to believe the rescale slope and intercept might be different for each slice?

[mohd akmal masud]

I'm, not sure that sir. That I guess

If just only one rescale slope for all slices, then how to multiply all the pixel number of each slice with rescale slope(in dicominfo)?

Answer 1

[Rik]

In general you should remember to include the intercept as well, although apparently that is 0 for this dataset.

unzip jaszakspect.zip
IM = dicomread('jaszakspect.dcm');
size(IM)

ans = 1×4

64 64 1 64

info = dicominfo('jaszakspect.dcm');
[min(IM(:)) max(IM(:))]

ans = 1×2

0 32767

IM = double(IM)*info.RescaleSlope + info.RescaleIntercept;
[min(IM(:)) max(IM(:))]

ans = 1×2

0 0.0029