Simulation of CT-images: from X-ray spectrum to CT-numbers Simulation of CT-images: from X-ray spectrum to CT-numbers S. Qamhiyeh 1) , A.Wysocka-Rabin 2 ) , M. Ellerbrock 1) , C. Penssel 3) , M. Kachelriess 3) and O. Jaekel 1) 1) German Cancer Research Center (DKFZ), Medical Physics, Heidelberg, Germany 2) Andrzej Soltan Institute for Nuclear Studies, Department of Accelerator Physics, Swierk, Poland 3) Friedrich-Alexander-University, Institute of Medical Physics, Erlangen-Nuremberg, Germany MCNEG / 2007 NPL, Teddington 28.03.2007 [email protected][email protected]
19
Embed
Simulation of CT-images: from X-ray spectrum to CT …resource.npl.co.uk/docs/science_technology/ionising radiation/clubs... · Simulation of CT-images: from X-ray spectrum to CT-numbers
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Simulation of CT-images:from X-ray spectrum to CT-numbers
Simulation of CT-images:from X-ray spectrum to CT-numbers
S. Qamhiyeh1), A.Wysocka-Rabin2), M. Ellerbrock1), C. Penssel3), M. Kachelriess3) and O. Jaekel1)
1) German Cancer Research Center (DKFZ), Medical Physics, Heidelberg, Germany2)Andrzej Soltan Institute for Nuclear Studies, Department of Accelerator Physics, Swierk, Poland3) Friedrich-Alexander-University, Institute of Medical Physics, Erlangen-Nuremberg, Germany
Uncertainties of CT-numbers can be up to 300 HU for bone-like materials. This can lead to uncertainties in the range of ions (position of Bragg-peak)
Changes of the CT-numbers of bone-like materials affect the calibration curve and cause up to 6 mm change in the range of Carbon ions. Up to 80% of the dose is allocated to the critical organs.
What influences CT-numbers ?What influences CT-numbers ?Design parameters of the CT-scanner (filters, collimators, generation)
Voltage of the X-ray tube
Size of Field of view (FOV)
Reconstruction algorithms
Parameters of the scanned object (material of phantom, size, geometry, homogenity, position in the FOV, movement)
Artefacts, especially metal artefacts
The dependance of CT-numbers on the above parameters is only known emprically.
Aim ...Aim ...
Develop a method to calculate CT-numbers under different conditions.
Analyze the effects of measurement conditions on calculated CT-numbers.
Estimate range uncertainties due to different measurement conditions of CT-numbers.
Cylindrical phantom from PMMA (16 cm diameter) and Gammex-
RMI substitutes (Middletown, WI).
Siemens Emotion CT equipped with a DURA352 X-ray tube.
Photon transport simulation throughthe filters,collimators and phantom
in CT scanner using EGSnrc/BEAMnrc
Reconstruction
Measurements of CT spectrum
using SpectroX
X-ray-tube
simulation
Responsefunction
of CT detectorsMeasurements of
CT-numbers
Projections and images
HU (Normalization ?)
Single projection Matrix – image
1
10
100
1000
1 65 129 193 257 321 385 449 513 577 641
Detector channels (672)
Flue
nce
per d
etec
tor c
hann
el
Fluence after phantom with Cortical bone insert
Initial fluence
Input CT-spectra calculated for the Emotion CT-scanner with DURA352 X-ray tube for 80, 110 and 130 kVInput CT-spectra calculated for the Emotion CT-scanner with DURA352 X-ray tube for 80, 110 and 130 kV
0.E+00
2.E-04
4.E-04
6.E-04
8.E-04
1.E-03
0 20 40 60 80 100 120 140Energy (KeV)
Cou
nts
(arb
itrar
y un
its)
80kV
110kV
130kV
0.E+00
1.E-02
2.E-02
3.E-02
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
energy (MeV)
fluen
ce/M
eV/in
ciden
t par
ticle
(cm
- 2 M
eV- 1
)
130 kV110 kV80 kV
0.00E+00
1.00E-03
2.00E-03
3.00E-03
4.00E-03
5.00E-03
6.00E-03
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
energy (MeV)
fluen
ce/M
eV/in
cide
nt p
artic
le (c
m- 2
MeV
- 1)
130 kV110 kV80 kV
Spectra in the center of FOV
Spectra at the edge of FOV
Input spectra
The effect of filters (Aluminium and Teflon) on the fluence and mean energy of the spectrum were investigated using 110 kV spectrum
The effect of filters (Aluminium and Teflon) on the fluence and mean energy of the spectrum were investigated using 110 kV spectrum
Effect of voltage, phantom size and phantom material on Carbon rangeEffect of voltage, phantom size and phantom material on Carbon range
0.00.20.40.60.81.01.21.41.61.82.0
-1000 -500 0 500 1000 1500 2000CT-numbers (HU)
Rel
ativ
e ra
nge
(WED
/T)
110 kV, PMMA, 16 cm
120 kV, PMMA, 16 cm
100 kV, PMMA, 16 cm
110 kV, PMMA, 30 cm
110 kV, Air, 16 cm
ConclusionsConclusions
It is POSSIBLE to calculate CT-numbers from MC simulation
The presented method was used to calculate the effects of different measurement parameters on CT-numbers and Carbon range
It is necessary to improve the effeciency of the MC simulation to use more than one projection for reconstruction (time needed foreach projection is 1-3 hrs of CPU time) to decrease the uncertainties of calculted CT-numbers
The simulation can not be used to calculate CT-numbers of irregular or off-center objects
Scheme of CT-scanner and phantom Scheme of CT-scanner and phantom
1) CT1-part was used to simulate X-ray source, filters and collimators. 1010
photons with a given energy spectrum. 1 week of CPU time.
2) CT2-part was used to simulate phantom and iserts. 5x108 photons. Input is PHSP from CT1. 1-3 hrs of CPU time.
Thank you for your attentionThank you for your attention
Questions ???Questions ???
Clinical Ion Therapy facility in Heidelberg (HIT)University Heidelberg / GSI / DKFZ1000 patients / yearp, He, C, O