Computing SUV PET Dicom image

35 views (last 30 days)
Alex marco
Alex marco on 8 Oct 2015
Commented: Sophia Lerebours on 20 Aug 2021
Hi, is there any reference for computing SUV from PET Dicom image?(matlab/c++)
  1 Comment
Wenhui Zhang
Wenhui Zhang on 26 Jan 2021
Hi, do you know the answer now? I am quite new to this area and if you know now, could you please tell me the answer.

Sign in to comment.

Answers (2)

Caroline Jordan
Caroline Jordan on 2 Jan 2018
OsiriX uses the following formula:
SUV (g/ml) = Pixel value (Bq/ml) * Weight (kg) / Dose(Bq) * 1000 (g/kg)
  1 Comment
antogeo
antogeo on 8 Apr 2018
Hello, Do you maybe know if "Dose" is corrected for time of acquisition or reads the dicom field: [0018, 1074]"Radionuclide Total Dose"?

Sign in to comment.


Cata_Ray
Cata_Ray on 5 Aug 2021
Edited: Cata_Ray on 16 Aug 2021
Using the formula:
The Dose should be corrected in relation to the injection and acquisition time. That information can be found in the DICOM Header and access via "dicominfo".
In this example I call the DICOM Header structure "header"
PET = %load PET imgage;
% Calculate SUV factor from DICOM header
delta_time = (str2double(header.AcquisitionTime) - str2double(header.RadiopharmaceuticalStartTime)) / 100; % [min]
half_life = header.RadionuclideHalfLife/ 60; % [min]
corrected_dose = header.RadionuclideTotalDose * exp(- delta_time * log(2) / half_life); % [Bq]
SUV_factor = (header.RescaleSlope * str2double(header.patient_weight) * 1000)/ (corrected_dose); % [g/Bq] = [] * [Kg]* 1000[g/kg] / [Bq]
% Create SUV image
PET_SUV = double(PET) * double(SUV_factor) ; %[g/ml] = [Bq/ml] * [g/Bq]
  3 Comments
Cata_Ray
Cata_Ray on 16 Aug 2021
It is not needed for this method of SUV calculation. You can see the formula at the begining of the answer.
Is there a praticular reason for you to want to use it?
Sophia Lerebours
Sophia Lerebours on 20 Aug 2021
I included the dose calibration factor in my calculations becuase of the paper Shin et al., 2017. https://pubmed.ncbi.nlm.nih.gov/28551298/

Sign in to comment.

Categories

Find more on DICOM Format in Help Center and File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!