how can I obtain hounsfield units(HU) from a .dcm image?
80 次查看(过去 30 天)
显示 更早的评论
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 个评论
采纳的回答
drummer
2020-5-7
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 个评论
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 ?
更多回答(1 个)
Marco Castro
2017-10-14
编辑:Marco Castro
2017-10-14
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
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 3-D Volumetric Image Processing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!