Computing SUV PET Dicom image
Mostra commenti meno recenti
Hi, is there any reference for computing SUV from PET Dicom image?(matlab/c++)
1 Commento
Wenhui Zhang
il 26 Gen 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.
Risposte (2)
Caroline Jordan
il 2 Gen 2018
0 voti
OsiriX uses the following formula:
SUV (g/ml) = Pixel value (Bq/ml) * Weight (kg) / Dose(Bq) * 1000 (g/kg)
1 Commento
antogeo
il 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"?
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 Commenti
Sophia Lerebours
il 16 Ago 2021
Is there a particular reason that you did not include the dose calibration factor in your calculations?
Cata_Ray
il 16 Ago 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
il 20 Ago 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/
il 8 Ott 2015
il 20 Ago 2021
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 