joint ictp-twas workshop on portable x-ray analytical...
TRANSCRIPT
![Page 1: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/1.jpg)
2455-5
Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage
Diego Dreossi
29 April - 3 May, 2013
Elettra, Trieste Italy
Lecture NoteBasic principles of X-ray Computed Tomography
![Page 2: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/2.jpg)
Diego DreossiSincrotrone Trieste, S.S. 14 km 163.5
AREA Science Park 34012 Basovizza (TS), Italy
Basic principles ofX-ray Computed Tomography
Workshop on Portable X-ray Analytical Instruments for Cultural HeritageICTP, Trieste, 29 april - 3 may, 2013
![Page 3: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/3.jpg)
Thanks to the big CT community for all the informationand material available on the web!
Acknowledgements
![Page 4: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/4.jpg)
Tomography
The word tomography is derived from the Greek ("part" or "section") and ("to write")
Projected image
Sections
Object
![Page 5: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/5.jpg)
4 7
5 3
11 11
8 8
11
8
4 7
5 3
75
7
8+715
8+513
11+718
11+718
8+715
8+513
11+718
11+718
![Page 6: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/6.jpg)
35
74
35
74
9 109 10
15+1025
13+922
18+1028
18+927
15+1025
13+922
18+1028
18+927
4 7
5 3
3
12
4
25+328
22+1234
28+1240
27+431
25+328
22+1234
28+1240
27+431
![Page 7: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/7.jpg)
31 40
34 28
12 2115 9
4 75 3
-19
/3
![Page 8: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/8.jpg)
![Page 9: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/9.jpg)
![Page 10: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/10.jpg)
Basic interactions
![Page 11: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/11.jpg)
x
No
xio eNN
xio eNN )( 321
Ni
x
Attenuation
NoNi
![Page 12: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/12.jpg)
o
i
kk N
Nx lno
i
NNdxx ln)(
x
iok
k
eNN
x
NoNi
Line integral
![Page 13: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/13.jpg)
The projection function is also determined by the angle of view .
x
y
t
p (t)
t1
A
B
ds
![Page 14: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/14.jpg)
x
y
A
B
y1
x1
t1
The equation of line AB is x cosy sin t
![Page 15: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/15.jpg)
The projection function can be written as
dydxtyxyxtp
dsyxtplinet
)sincos(),()(
),()(),(
p(t) is known as the Radon Transform of the function x,y).
A projection is formed by combining a set of line integrals. The simplest projection is a collection of parallel ray integrals as given by p(t) for a constant . This is known as parallel projection.
detector
Laminar beam
object
![Page 16: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/16.jpg)
Image reconstruction algorithms are derived to constructx,y) from p(t).
Classification of Algorithms
Backprojection
Fourier Domain Approach
Filtered Backprojection
Iterative Methods
Algebraic Reconstruction Technique (ART)
Iterative Least Squares
Simultaneous IterativeReconstruction Technique (SIRT)
![Page 17: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/17.jpg)
![Page 18: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/18.jpg)
![Page 19: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/19.jpg)
The above given process can be expressed mathematically. The reconstructed back-projection image b (x,y) at a particular view is :
dttyxtpyxb )sincos()(),(
Adding up the images at all angles (0-)
ddttyxtp
dyxbyxfb
0
0
)sincos()(
),(),(
![Page 20: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/20.jpg)
dudvevuFvuFFyxf
dxdyeyxfyxfFvuF
vyuxj
vyuxj
)(21
)(2
),(),(),(
),(),(),(
f(x,y) F(u,v)
ImageSpace
FourierSpace
![Page 21: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/21.jpg)
![Page 22: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/22.jpg)
•Thus the Fourier Transform of a projection at angle forms a line in the 2-D Fourier plane at this angle.•After filling the entire F(,) plane with the transforms of the projections at all angles, the object can be reconstructed using 2-D Inverse Fourier Transform
![Page 23: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/23.jpg)
Convolution Back-projection Algorithm
The back-projected function can rewritten in space domain as follows:
)(
)()( 111
11
tc
FtptpFF
Thus, instead of filtering in the frequency domain, p (t) can be convolved with a function c (t) and then back-projected.
![Page 24: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/24.jpg)
-0.2
-0.1
0
0.1
0.2
0.3
-10 -5 0 5 10
RamLakShepp-Logangen-Hamming
ampi
ezza
(u.a
.)
pixel
0
0.1
0.2
0.3
0.4
0.5
0 0.1 0.2 0.3 0.4 0.5
rampa
Shepp-Logan
gen-Hamming
ampi
ezza
(u.a
.)
frequenza (1/step)
• Ramp filter (RamLak): enhancement of high frequencies noise
• Gen-Hamming, Shepp-Logan: enhancement ofintermediate frequencies
• Convolution theorem convolution in the direct space asan alternative to multiplication in the Fourier space
![Page 25: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/25.jpg)
![Page 26: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/26.jpg)
SinogramImage
ProjectionProjection
t
t
![Page 27: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/27.jpg)
Sinogram
Back
Pro
ject
ion
Filtration
Image
![Page 28: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/28.jpg)
![Page 29: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/29.jpg)
Reconstruction in frequency domain
![Page 30: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/30.jpg)
ART
x y
z w
12
8
791113
![Page 31: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/31.jpg)
First generation
• EMI Mark I (Hounsfield), “pencil beam” or parallel-beam scanner• 180° - 240° rotation angle, angular step ~1 °• Scan time 5 min, reconstruction time 20 min• Resolution: 80 x 80 pixels (ea. 3 x 3 mm2),• Slice thickness 13 mm
![Page 32: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/32.jpg)
Second generation
• Hybrid system: Fan beam + linear array (~30 elements)• traslation and rotation• total scan time ~30 s• more complex agorithms (“fan” geometry)
![Page 33: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/33.jpg)
Third generation• the fan beam is covering all the sample• 500-700 elements (ionizating chambers or scintillators)• No traslations• total scan time ~ seconds
• reconstruction time ~ seconds
![Page 34: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/34.jpg)
Fourth generation
• Stationary ring of detectors (600 – 4800 scintillators)• Rotating X-ray source • total scan time ~ seconds• reconstruction time ~ seconds • Slice thickness 1mm
3D volumes constructed as a series of 2D slices
Fifth: electron beam scanner
![Page 35: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/35.jpg)
Helical (sixth generation)
![Page 36: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/36.jpg)
Conventional CT Helical CT
scanning N 360° scans at positions z1 to zn
One scan of n-360°from positions z1 to zn
Pre-processing corrections corrections
intermediate Z-interpolation
reconstruction Convolution and backprojection
Convolution and backprojection
result N images at positions z1 to zn
Images at arbitrary positions from z1 to zn
![Page 37: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/37.jpg)
Acq. time
Data 360°
Data helical
Matrix
Power
Slice thick.
Spatial res.
Contrast res
![Page 38: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/38.jpg)
Status3D medical CT:
• Helical trajectory• Similar to 3rd generation CT but with multiple rows of detectors (4, 8, 16, 32, 64, now even up to 640 rows)• FDK-like approximate reconstruction
3D lab-based CT:• 3D cone-beam micro-CT using circular trajectory• 5122, 10242, 20482 (flat-panel – CCD – II detectors-…)• FDK approximate reconstruction
![Page 39: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/39.jpg)
SourceEnergy Current Focal spot sizeFocal spot stability
Magnificationx-ray fluxPenumbrafield of view
DetectorScreen: efficiency, spatial resolutionPixel size: spatial resolution,signal, FovDynamic range, daq time
Design concepts
![Page 40: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/40.jpg)
Source
Focal spot sizeEmission angleMaterialTechnology
≈5° - 20°
![Page 41: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/41.jpg)
Source
FSM
MUFS1
Heat dissipationStabilityPhase effects
unsharpness
![Page 42: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/42.jpg)
Source
1.50 mm AlMean energy 46.2 keV1st HVL 2.38 mm Al, 2nd HVL 4.20 mm Al, HVL1/HVL2 = 0.567
0.25 mm CuMean energy 57.3 keV1st HVL 6.46 mm Al, 2nd HVL 8.00 mm Al, HVL1/HVL2 = 0.807
CharacteristicRadiation lines
Internal filtering
100 kVp, = 15°, 0.2 mm Be window500 mm air
![Page 43: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/43.jpg)
Source
![Page 44: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/44.jpg)
1.E-01
1.E+00
1.E+01
1.E+02
1.E+03
1.E+04
1.E+05
0 50 100E (keV)
c
m
Bone
Calcium Carbonate
paraffin wax
1 2
Lower energy Higher energy
1.0
0.2
x
Rel.
inte
nsity
0.1 – 0.2 transmission > 5 n
x
iok
k
eNN
![Page 45: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/45.jpg)
Detector
Xray
Detectorscreen lightguide
Direct (aSe, Si)PowdersStructured
FovHigh energy
![Page 46: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/46.jpg)
Detector CCD Area, sCMOS
Flat Panel (large area CMOS)
Flat Panel (aSi/aSe)
Photodiode array
Single photon counting
GOS (+FOP)CsI (+FOP)
Direct conversion
CsICeramic
MediPix, Pixirad, ….
![Page 47: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/47.jpg)
Detector
Dynamic range 72 db ~ 4000:1 -> 12 bit
![Page 48: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/48.jpg)
Detector
LAG: not suitable for CT!
![Page 49: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/49.jpg)
Reconstruction
Volume 150 x 150 x 150 mm3
Isotropic voxel size 50 mdataset = 30003 x 4 = 108 GB
DATA I/O / PROCESSING / STORAGE
Isotropic voxel size 60 mdataset = 25003 x 4 = 62.5 GB
![Page 50: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/50.jpg)
Reconstruction
Cone beamFanHelicalHelical+conepCTHBCTLocal Area CT
Artifacts reduction
Pre-processing Post-processing
![Page 51: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/51.jpg)
Geometry
Intensity~1/d2
FDK filtered back-projection is reasonable up to ~10°FDK filtered back-projection is reasonable up to ~10°
With microfocus sources phase contrast imaging possible
![Page 52: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/52.jpg)
Unsharpness
PSM
UPS2
22PSFSTOT UUU
3 333 )6.1()1(1 PSFSMM
UTOT
ASTM E 2698
Pixel Size FSM
MUFS1
Focal spot
ExamplesFS = 8 m, PS = 25 m, M =5
UTOT= 9.2 m PSeq = 5 m
FS = 16 m, PS = 50 m, M =1.25
UTOT= 64 m PSeq = 40 m
![Page 53: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/53.jpg)
Number of projections
K ≈ /2 N Reconstruction of an ellipse
Spatial and angular samplingon Fourier space
K Projections
N S
ampl
es
64 128 256 512
512
25
6
128
6
4
![Page 54: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/54.jpg)
NatINpPSpn
1
pn pixel noise (standard deviation evaluated on 2D slices)
k = constant dependend on rec algorithm (smooth kernel reduces noise, but also spatial resolution)
PS = pixel size (mm2) Np = number of projectionsI = source current (mA)t = integration time of the detector (s)Na = image averaging number
![Page 55: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/55.jpg)
Artifacts
Sample 1- motion2- metallic materials3- dimension > scan field
Physics1- Beam Hardening
Cupping & Streaks2- Undersampling3- Photon starvation4- Partial Volume
Scanner1- rings2- geometry3- stability
Technique1- helical in transverse plane2- helical in multislice3- multiplanar and 3D reconstrucion4- cone beam effect
![Page 56: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/56.jpg)
Artifacts simulation
![Page 57: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/57.jpg)
Artifacts simulation
![Page 58: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/58.jpg)
PIXEL
I
single pixel traversed by individual rays:
•will underestimate the attenuation
Motion artifactsrod moved during acquisition
Artifacts simulation
![Page 59: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/59.jpg)
More realistic cases
Beam hardening: cuppingBeam hardening: streaks
Photon starvationUndersampling
![Page 60: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/60.jpg)
More realistic cases
Rings
Metal artifacts
![Page 61: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/61.jpg)
More realistic cases
Dimensions > FOV
![Page 62: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/62.jpg)
Imaging RequirementsSomething to know:
1. The size of the object you want to scan
2. The material the object is made of
3. The level of fine detail or density differences you want to see
4. How fast you want to see the results
5. Where you want to place the instrument
6. Who will use it
![Page 63: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/63.jpg)
New high-speed Computed Tomography system for 3D mass production process control
![Page 64: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/64.jpg)
![Page 65: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/65.jpg)
![Page 66: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/66.jpg)
![Page 67: Joint ICTP-TWAS Workshop on Portable X-ray Analytical ...indico.ictp.it/event/a12177/session/12/contribution/7/material/0/0.pdf · X-ray Computed Tomography Workshop on Portable X-ray](https://reader033.vdocument.in/reader033/viewer/2022053101/606089ceb9169b17cb2dca3d/html5/thumbnails/67.jpg)