how can I obtain hounsfield units(HU) from a .dcm image?

Question

sara 2014-9-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/153570-how-can-i-obtain-hounsfield-units-hu-from-a-dcm-image

评论： Walter Roberson 2022-5-24

hi

if I have a CT image how can I obtain hounsfield units from this... I know that we can read .dcm Image with dicomread(img)..

and I should say when I use dicomread(img) and then I want to imshow this I just see a black screen...

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

drummer 2020-5-7

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/153570-how-can-i-obtain-hounsfield-units-hu-from-a-dcm-image#answer_431061

在 MATLAB Online 中打开

After loading your image:

yourImage = dicomread('yourImage.dcm');

You can go straight to

info = dicominfo('yourImage.dcm');

Now, you have to find your DICOM attribute that corresponds to the Hounsfield units.

This attribute leads you the the linear correlation between the voxel value in the image and the Hounsfield units.

The DICOM attribute I´m talking about is (0028, 1053), or the 'Rescale Slope' attibute.

H_unit(x,y) = voxel_value(x,y)*RescaleSlope + Intercept.

It is exactly like y = a.x + b.

You can do like this:

rSlope = info.RescaleSlope;
for j = 1 : size(yourImage, 1) % This loop multiply each voxel value by the rescale slope
    for i = 1 : size(yourImage, 2)
        hounsfieldImage(i,j) = yourImage(i,j)*rSlope;
    end
end
figure
imshow(hounsfieldImage, 'DisplayRange', []);

Cheers

7 个评论
显示 5更早的评论隐藏 5更早的评论

yildiz 2022-5-24

thank you for your kindly reply. ı can not see resclae slope and intercept from matlab but in microdicom viewer rescale slope is 1 and rescale intercept is -8192. I wrote this value manually in my code but ı think this value is wrong. Can you help me ?

Walter Roberson 2022-5-24

has the image already had water suppressed?

请先登录，再进行评论。

Answer 2

Marco Castro 2017-10-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/153570-how-can-i-obtain-hounsfield-units-hu-from-a-dcm-image#answer_285754

编辑：Marco Castro 2017-10-14

在 MATLAB Online 中打开

Hi, well, first you need to do the reading of the image:

Image1=dicomread('Image.dcm');

After you convert the dicom image in uint16:

 Image1 = uint16(Image1); % This is very important, if you don't declare this, can make you troubles

After you make a vector of the size of the image:

[m,n] = size(Image1);

After that you need two matrices of zeros of the size of the image(after this you will see the use of both)

B1 = zeros(m,n);
S1 = zeros(m,n);

Then you make a vector with the values of the Hounsfield units, example

 Bone = [200,300,3000]; % The information of this i take it from wikipedia https://en.wikipedia.org/wiki/Hounsfield_scale

And then you do two for loops, these will be repeated according to the size of the image, of rows and columns

for i=1:m
  for j=1:n
      P=Image1(i,j); %This contain the information of the image
      if(P>Bone(1,1) && P<Bone(1,2)) % Make the condition that allows compare the information in the variable "P"
          B1(i,j) = 1; % Make a bit image and change the values of zeros for ones
          S1(i,j) = P; % This transform the value of P that contain the image information
      end
  end
end
B1=uint8(B1); % Transforms the double-to-single bit image
S1=uint8(S1*255);  % This is the same that the line above but multiplies the result for 255
And then for show it:
figure(1)
subplot(2,1,1);
imshow(S1);
title('Image treated');
subplot(2,1,2);
imshow(Image1*255); % This multiplies the values of the image just to observe them
title('Original Image');

I tried and it work out, but is in the bones from a CT in the skull, I hope it helps you.

2 个评论
显示无隐藏无

Walter Roberson 2017-10-14

This does not look anything like I would expect? Generally calculation of HU involves looking at the dicom headers and doing a linear transformation. The process is described at https://www.mathworks.com/matlabcentral/answers/114123-how-to-calculation-hounsfield-unit#answer_127872

boloroo bolor 2018-4-3