Instrumentation#
X-ray production#
Energy distribution#
The principle of an x-ray tube is the following. Electrons in a heated filament, with a current of a few mA typically, are extracted using a very high voltage \(U\) (several tens of kV typically). Their kinetic energy is close to 0 in the filament (they can be considered at rest). Therefore their kinetic energy reach \(U\) keV just before impinging onto the anode target. From electron-matter interactions in the target follow X-ray mostly produced by Bremsstrahlung but also to a lesser extent Fluorescence. Refer to section Electron-matter interactions.
A monoenergetic electron beam impacting a target thin metal produces a uniform energy spectrum of Bremsstrahlung x-rays. For a thick anode, the radiant-energy spectrum of emitted photons is a function linearly decreasing.
Thin-target model of the radiant-energy spectrum (excerpt from Attix DOI):
Consider a beam of electrons of kinetic energy \(T_o\) entering the page perpendicularly, and each passing the nucleus at some distance (impact parameter) \(b\) (see section Electron-matter interactions). The differential interaction cross-section when \(b = b_1\) is proportional to the area \(dA_1 = 2\pi b_1 db\). For \(b_2 = 2b_1\), \(dA_2 = 2 dA_1\). Thus twice as many photons \(N_2\) come from interactions in \(dA_2\), as the \(N_1\) from \(dA_1\). If the magnitude of the interaction (i.e., the x-ray quantum energy \(hv\) produced) is assumed to be proportional to \(1/b\), then \(hv_1 = 2 hv_2\), Therefore \(N_1hv_1 = N_2hv_2\), and the x-ray radiant-energy spectrum should be flat, as it is observed to be.
The high voltage accelerating electrons conditions the maximum energy value of the emitted photons, while the current only changes the photon rate without modify the maximum value.
The linear coefficient of the materials being more important at low energies, the spectrum exiting the anode target no longer has its low-energy components (spectrum B of Fig. 90).
X-ray spectra emitted at 60 and 120 kV high voltage. We sees the characteristic lines of the anode fluorescence in tungsten following ionizations by the electron beam.
Beam hardening#
When x-rays pass through matter, the low-energy components are more attenuated than the high-energy ones because the linear attenuation coefficient globally decreases when the photon energy is increased. Figure Fig. 92 shows how a 60kV x-ray spectrum is affected by a 150-micron-thick cupper foil.
The signal measured \(I\left(p\right)\) by the flat panel detector is proportional to the deposited energy in the active area of the pixel \(p=\left(x,y\right)\). Assuming there are \(M\) different materials crossed by the x-rays from the source to pixel \(p\), the image \(I\left(p\right)\) can be modeled as follows:
where \(E\) is the photon energy, \(D\left(E,p\right)\) is the detector response, \(N_0\left(E,p\right)\) the number of photons emitted by the source, \(\mu_m\left(E\right)\) the linear attenuation coefficient of material \(m\) and \(X_m\left(p\right)\) the corresponding material thickness. If we have an homogeneous material \(M\) this boils down to
We can see that since the primary spectrum is not monochromatic the measurement is the sum of several attenuation laws. We could however compute the effective monochromatic attenuation that would give the same measurement:
with
where $E_{\mathrm{eff}} corresponds to the monochromatic energy that would give the same attenuation than the one measured.
Illustration of beam hardening by adding successively layers of 1/2 attenuation (\ emph {HVL}) in aluminum. We notice that we need a \ emph {HVL} more and more important to stop half the photons, and that the average energy of the photons in the spectrum also increases.
When the material has no attenuation discontinuity in the energy range of the spectrum, such as aluminum previously, the transmitted spectrum is deformed but remains continued. On the other hand, the presence of discontinuity photoelectric, as here for tin, goes strongly modify the shape of the transmitted spectrum.
Standard characteristics#
The radiative efficiency of the generator is very bad and the vast majority of electrical power is converted into heat at the anode. Summary of typical characteristics of X-rays.
HT: \([20-800]\) kV
Power: \(1\) W – \(5\) kW
Focus: \(1\) micron – \(5\) mm
Absorbed power \(= \text{kV} \times \text{mA}\)
Radiative power \(\approx 10^{-6} \times Z \times \text{kV}^2 \times \text{mA}\) (NB : Bremsstrahlung slope \(\propto Z\))}
Yield \(\approx 10^{-6} \times Z \times \text{kV}\) E.g. @ 60 kV (W anode): 99.5% heat, and 0.5% radiative
mA \textit{vs} e-
Anode temperature (1 calorie = the amount of energy required to warm one gram of water from 19.5 to 20.5°C at standard atmospheric pressure)
The standard components of the generator are: the high-voltage transformer(s), the x-ray tube, the high-voltage cable and the controller (for the setting of high-voltage and current).
Linear accelerator#
Beyond 450 kV, it becomes difficult to accelerate electrons by a single high voltage. So it is necessary to accelerate them in successive stages for the bring to the desired energy. This is the principle of the linear accelerator which can deliver x-rays of several MeV. These energies are needed as soon as the thickness to be traversed exceeds a few centimeters of metal.
Radioactive sources#
Detectors#
Films and Image Intensifiers#
Analog films were used initially. Flexible, they can tackle the pipes, but require a chemical process of development. More versions modern non-sensitive to ambient light are used for quality assurance of beams in radiotherapy.
The photostimulable plates have the same function as the silver films but are reusable and do not require no chemical process.
Image intensifiers have long allowed to perform dynamic real-time imagery. They present very good resistance to radiation.
Flat panels#
Flat screens are now the benchmark in digital X-ray imaging. The reading electronics being in the field irradiation, special precautions should be taken to protect it. Two technologies exist: either indirect by the conversion into optical photons in the so-called scintillator material (thanks to a particular doping which allows to have a visible fluorescence following the effect photoelectric), or direct using the electrons expelled as a result of interactions in the detector (but this requires electrodes at each pixel). The vast majority of detectors are in integration, without the possibility of photon counting (because this requires a spectrometric chain at each pixel and a ultra fast electronics).
Standard characteristics#
Characteristics:
Mode: integrating
Readout: 1-100 FPS
Pixel sizes: 1-1000 microns
Dynamic range: 8-20 bits
\(< 20\) MPixels
Mode: counting
Countrate: 2-200 Mcps / mm\(^2\)
Pixel sizes: 50-1000 microns
No. of threshold: 1-100 (typ. 4)
\(< 1\) MPixels
Geometry: flat panel or linear
The typical energy response of some flat screen in the market is given in Fig. 109. At low energy, interactions are very probable and the detector has a response close to 100%, ie total absorption. When the energy is too great, X-rays pass through without interacting in the detector, which drop efficiency. Note the discontinuity of the CsI scintillator screen absorption (\(K_{\alpha}\) edge is \(33\) keV for I and \(36\) keV for Cs) and of the lead-glass shielding (\(K_{\alpha}\) edge is \(88\) keV for Pb).
Optical coupling#
The protection of the reading electronics against radiation requires optical coupling solutions remote, with a tandem optic (as in Fig. 110) or a bundle of fibers (as in Fig. 111), but this lowers the detection efficiency compared to a flat screen for which the solid angle for the photodiodes is close to \(2\pi\) sr.
Fig. 112 shows the simplified architecture at the level of a pixel to get the x-ray energy spectrum and not a single integral scalar value.
Photon counting#
The temporal response of a single photon is relatively long (several tens of ns), which can quickly pose stacking problems if the photon flux is too much important (see Fig. 113): this phenomenom is called pile-up. Charge sharing of a single photon deposit between neighboring pixels may also occur, which results in the detection of several photons of lower energy.
A photon does not necessarily deposit all of its energy: this is only the case if there is interaction by effect photoelectric. During a Compton interaction only one fraction of its energy is deposited: how to go up then to the incident energy of the photon if we want to have the energy distribution? X-ray “color” imaging is not easy to master. Fig. 114 shows the detected radiation spectrum of \(^{22}_{11}\mathsf{Na} \Rightarrow ~^{22}_{10}\mathsf{Ne} + (180.8{\%}) \beta^+ + (99.9{\%}) \gamma_{1.27\mathsf{MeV}}\).
The detected spectrum would be the incident emitted one if all photons had deposited all their energy by photoelectric effect. Instead, many photons did interact by Compton scattering and only part of their energy was deposied in the detector: the maximum energy deposition by the Compton electron is when \(\theta=\pi\), which is called the Compton edge (see Fig. 114).
Flat Field#
X-ray detectors have the same characteristics and limits than standard imaging systems: afterglow, the dynamics of digitization, spatial heterogeneities and temporal… Normalization by an objectless image (flat field) corrects a part of the artifacts.
Scatter sensitivity \(\Rightarrow\) slot-scanning vs flat panel
Veiling glare due to secondary energy carriers (optical photons or e\(^-\)) \(\Rightarrow\) drop at low frequencies
Dynamic range
Spatial artifacts \(\Rightarrow\) defective pixel map
Temporal artifacts \(\Rightarrow\) lag (residual signal), flicking pixels
Temporal stability \(\Rightarrow\) varying performances over time
The image formation process may therefore be approximated as follows for a mono-material and mono-energy:
where \(p\) is the pixel, \(I_{\mathsf{Source}}\) the x-rays emitted by the source, \(R_{\mathsf{Detector}}\) the detector response, \(\mu\) the effective attenuation coefficient, \(X\) the material thickness, and \(I_{\mathsf{Dark}}\) the dark current. The flattening makes it possible to cope with those spatial variations of x-ray emission and detector response.
Flattening
Magnification & Unsharpness#
The x-ray image is a projection of the attenuation and the size of the projection image is larger than the size of the object: there is a magnification \(M \geq 1\) which can be measured as the ration of the source-to-detector distance over the source-to-object distance (see image Fig. 116):
The blurring in image X has two causes: the detector itself (the inherent unsharpnness) and the non-punctuality of the source of radiation (geometric unsharpness). If the focal task of the source is oriented, the geometric unsharpness may depend on the position on the detector.
Total unsharpness \(u^2_{\mathsf{total}}\)
The geometric unsharpness is directly related to the source size \(f\) up to the magnification \(M\)as Fig. 117 shows. Note that an angled anode generated a space-variant geometric unsharpness.
where, given the notation of figure Fig. 117, the magnification \(M\) is given by
The inherent unsharpness, also named point-spread function (PSF), depends on the detector technology (direct or inidrect), the scintillator (if indirect) thickness, the coating and mounting materials, the incident photon energy…
Image primitives#
SNR, CNR and SPR#
The contrast between two areas of the RX image is directly linked to the law of attenuation. A change in current or the generator voltage will directly impact it (modification of \(I_0\) and \(\mu\). Alone, contrast is not enough not to characterize the detectability, it is necessary to take a measurement related to noise in the radiograph. The presence of secondary radiation from interactions (scattered, fluorescence…) will also degrade the dynamics.
Image primitives
Absolute contrast
Relative contrast
Counting process: Poisson noise (mean = variance) for the nr of photons
SPR: Scatter-to-Primary Ratio
Nota Bene: remind that the generator settings (mA, s & kV) influence those primitives
MTF and DQE#
The contrast in the image also depends on the content spatial frequency of the object. The transfer function of modulation (MTF) makes it possible to characterize the evolution of contrast with the spatial frequency.
Fig. 123 shows a typical bar-pattern test object to assess the spatial resolution. It ia usually made of lead and PMMA. The higher the spatial frequency the lower the intensity contrast.
The drop of intensity contrast with respect to the spatial frequency is observed in Fig. 124.
The MTF value at frequency \(f\) is simply computed from the contrast ratio with respect to the contrast when the frequency tends to \(0\).
The ratio of the MTF to the noise power spectrum (NPS) gives the quantum detection efficiency (DQE), which ultimately represents a contrast on noise as a function of the spatial frequency.
we have
with \(N_0\) in [#phots~mm\(^{-2}\)] for a given exposure. Then
Nota Bene: the patient (scattering, magnification…) is now being taken into account to compute effective figures of merit.
Moiré or aliasing can also occur in X-ray imaging. It is directly related to the spatial sampling of the detector, ie the size of the pixel.