Convert pixel value to HU

Question

0 votes

How can I change a voxel value of a Micro CT image series to HU units ?

i need to change the voxel value and the header in dicom to use the new dicom in the treatment planning system

the linear equation i will be used to change the voxel is

y = 0.7966x-568.5

1 Comment
Show -1 older comments Hide -1 older comments

Ahmad Abbas on 7 Aug 2020

Open in MATLAB Online

% this is the code i used to create histogram, which help me to create the equation (tmp11) which then i convert the pixel value to HU ( i am not sure if it is correct or not ) 
% i am trying to save the result with the metadata and changes, please can someone help
clear
myFolder = 'C:\Users\Admin\Desktop\MATLAB\Crop cochlea';
filePattern = fullfile(myFolder, '*dcm');
theFiles = dir(filePattern);
crop=zeros(209,218,394);
for k = 1 : length(theFiles)
    baseFileName = theFiles(k).name;
    fullFileName = fullfile(myFolder, baseFileName);
    fprintf(1, 'Now reading %s\n',fullFileName);
    x = dicomread(fullFileName);
    crop(:,:,k)=x;
end
%histogram(crop)
tmp11 = min(max (0.7966*crop -568.5, -1000),4096); 
hu = uint16(tmp11);