Computing SUV PET Dicom image
조회 수: 11 (최근 30일)
이전 댓글 표시
Hi, is there any reference for computing SUV from PET Dicom image?(matlab/c++)
댓글 수: 1
Wenhui Zhang
2021년 1월 26일
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.
답변 (2개)
Caroline Jordan
2018년 1월 2일
OsiriX uses the following formula:
SUV (g/ml) = Pixel value (Bq/ml) * Weight (kg) / Dose(Bq) * 1000 (g/kg)
댓글 수: 1
antogeo
2018년 4월 8일
Hello, Do you maybe know if "Dose" is corrected for time of acquisition or reads the dicom field: [0018, 1074]"Radionuclide Total Dose"?
Cata_Ray
2021년 8월 5일
편집: Cata_Ray
2021년 8월 16일
Using the formula:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/713447/image.png)
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
Cata_Ray
2021년 8월 16일
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
2021년 8월 20일
I included the dose calibration factor in my calculations becuase of the paper Shin et al., 2017. https://pubmed.ncbi.nlm.nih.gov/28551298/
참고 항목
카테고리
Help Center 및 File Exchange에서 DICOM Format에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!