X-ray computed tomography operates by using an X-ray generator that rotates around the object; X-ray detectors are positioned on the opposite side of the circle from the X-ray source.
A visual representation of the raw data obtained is called a sinogram, yet it is not sufficient for interpretation. Once the scan data has been acquired, the data must be processed using a form of tomographic reconstruction, which produces a series of cross-sectional images. In terms of mathematics, the raw data acquired by the scanner consists of multiple "projections" of the object being scanned. These projections are effectively the Radon transformation of the structure of the object. Reconstruction essentially involves solving the inverse Radon transformation.
Structure
In conventional CT machines, an X-ray tube and detector are physically rotated behind a circular shroud (see the image above right). An alternative, short lived design, known as electron beam tomography (EBT), used electromagnetic deflection of an electron beam within a very large conical X-ray tube and a stationary array of detectors to achieve very high temporal resolution, for imaging of rapidly moving structures, for example the coronary arteries. Systems with a very large number of detector rows, such that the z-axis coverage is comparable to the xy-axis coverage are often termed cone beam CT, due to the shape of the X-ray beam (strictly, the beam is pyramidal in shape, rather than conical). Cone-beam CT is commonly found in medical fluoroscopy equipment; by rotating the fluoroscope around the patient, a geometry similar to CT can be obtained, and by treating the 2D X-ray detector in a manner similar to a CT detector with a massive number of rows, it is possible to reconstruct a 3D volume from a single rotation using suitable software.
Contrast media
Contrast mediums used for X-ray CT, as well as for plain film X-ray, are called radiocontrasts. Radiocontrasts for X-ray CT are, in general, iodine-based.[2] This is useful to highlight structures such as blood vessels that otherwise would be difficult to delineate from their surroundings. Using contrast material can also help to obtain functional information about tissues. Often, images are taken both with and without radiocontrast.
Schematic configuration and motion
In this section, the schematic configuration and motion of the parallel beam irradiation optical system configured to obtain the p(s,θ) of above-mentioned (eq. 5) will be explained. In this section, how to obtain the p(s,θ) of (eq.5) by utilizing parallel beam irradiation optical system will also be explained. Configuration and motions of parallel beam irradiation optical system, referring Fig. 3.
Statements
Numbers (1)–(7) shown in Fig. 3 (see the numbers within the parentheses) respectively indicate: (1) = an object; (2) = the parallel beam light source; (3) = the screen; (4) = transmission beam; (5) = the datum circle (a datum feature); (6) = the origin (a datum feature); and (7) = a fluoroscopic image (a one-dimensional image; p (s, θ)).
Two datum coordinate systems xy and ts are imagined in order to explain the positional relations and movements of features (0)–(7) in the figure. The xy and ts coordinate systems share the origin (6) and they are positioned on the same plane. That is, the xy plane and the ts plane are the same plane. Henceforth, this virtual plane will be called “the datum plane”. In addition, a virtual circle centered at the abovementioned origin (6) is set on the datum plane (it will be called “the datum circle” henceforth). This datum circle (5) will be represents the orbit of the parallel beam irradiation optical system. Naturally, the origin (6), the datum circle (5), and the datum coordinate systems are virtual features which are imagined for mathematical purposes.
The μ(x,y) is absorption coefficient of the object (3) at each (x,y), p(s,θ) (7) is the collection of fluoroscopic images.
Motion of parallel beam irradiation optical system
The parallel beam irradiation optical system is the key component of a CT scanner. It consists of a parallel beam X-ray source (2) and the screen (3). They are positioned so that they face each other in parallel with the origin (6) in between, both being in contact with the datum circle (6).
These two features ((2) and (3)) can rotate counterclockwise [Notes 1] around the origin (6) together with the ts coordinate system while maintaining the relative positional relations between themselves and with the ts coordinate system (so, these two features ((2) and (3)) are always opposed each other). The ts plane is positioned so that the direction from a collimated X-ray source (2) to the screen (3) matches the positive direction of the t-axis while the s-axis parallels these two features. Henceforth, the angle between the x- and the s-axes will be indicated as θ. That is, parallel beam irradiation optical system where the angle between the object and the transmission beam equals θ. This datum circle (6) will be represents the orbit of the parallel beam irradiation optical system.
On the other hand, the object (1) will be scanned by CT scanner is fixed to xy coordination system. Hence, object (1) will not be moved while the parallel beam irradiation optical system are rotated around the object (1). The object (1) must be smaller than datum circle.
Increment/Table speed
The distance the table moves for every 360° rotation of the X-ray generator is called the increment or table feed for axial scan modes. For helical scan modes, it is called table speed.[3] Setting an increment that is smaller than the slice thickness results in overlap between the slices. A beneficial effect of this is a smoother transition between images when scrolling through the stack.[4]
Obtaining transmission image ‘s’
During the above-mentioned motion (that is pivoting around the object(1)) of parallel beam irradiation optical system, the collimated X-ray source (2) emits transmission beam (4) which are effectively “parallel rays” in a geometrical optical sense. The traveling direction of each ray of the transmission beam (4) is parallel to the t-axis. The transmission beam (4), emitted by the X-ray source (2), penetrates the object and reaches the screen (3) after attenuation due to absorption by the object.
Optical transmission can be presumed to occur ideally. That is, transmission beam penetrates without diffraction, diffusion, or reflection although it is absorbed by the object and its attenuation is assumed to occur in accordance with the Beer-Lambert law.
Consequently, a fluoroscopic image (7) is recorded on the screen as a one-dimensional image (one image is recorded for every θ corresponding to all s values). When the angle between the object and transmission beam is θ and if the intensity of transmission beam (4) having reached each "s" point on the screen is expressed as p(s, θ), it expresses a fluoroscopic image (7) corresponding to each θ.
Tomographic reconstruction
The technique of filtered back projection is one of the most established algorithmic techniques for this problem. It is conceptually simple, tunable and deterministic. It is also computationally undemanding, with modern scanners requiring only a few milliseconds per image.However, this is not the only technique available: the original EMI scanner solved the tomographic reconstruction problem by linear algebra, but this approach was limited by its high computational complexity, especially given the computer technology available at the time. More recently, manufacturers have developed iterative physical model-based maximum likelihood expectation maximization techniques. These techniques are advantageous because they use an internal model of the scanner's physical properties and of the physical laws of X-ray interactions. Earlier methods, such as filtered back projection, assume a perfect scanner and highly simplified physics, which leads to a number of artifacts, high noise and impaired image resolution. Iterative techniques provide images with improved resolution, reduced noise and fewer artifacts, as well as the ability to greatly reduce the radiation dose in certain circumstances.[5] The disadvantage is a very high computational requirement, but advances in computer technology and high-performance computing techniques, such as use of highly parallel GPU algorithms or use of specialized hardware such as FPGAs or ASICs, now allow practical use.
Basic principle
In this section, the basic principle of tomography in the case that especially uses tomography utilizing the parallel beam irradiation optical system will be explained.
Tomography is a technology that uses a tomographic optical system to obtain virtual 'slices' (a tomographic image) of specific cross section of a scanned object, allowing the user to see inside the object without cutting. There are several types of tomographic optical system including the parallel beam irradiation optical system. Parallel beam irradiation optical system may be the easiest and most practical example of a tomographic optical system therefore, in this article, explanation of "How to obtain the Tomographic image" will be based on "the parallel beam irradiation optical system". The resolution in tomography is typically described by the Crowther criterion.
Fig. 3 is intended to illustrate the mathematical model and to illustrate the principle of tomography. In Fig.3, absorption coefficient at a cross-sectional coordinate (x, y) of the subject is modeled as μ(x, y). Consideration based on the above assumptions may clarify the following items. Therefore, in this section, the explanation is advanced according to the order as follows:
- (1)Results of measurement, i.e. a series of images obtained by transmitted light are expressed (modeled) as a function p (s,θ) obtained by performing radon transform to μ(x, y), and
- (2)μ(x, y) is restored by performing inverse radon transform to measurement results.
(1) The Results of measurement of p(s,θ) in a parallel beam irradiation optical system
Consider the mathematical model where the absorption coefficient of the object at each point is represented by the function and suppose that the transmission beam penetrates without diffraction, diffusion, or reflection. Also assume the beam is absorbed by the object and its attenuation occurs in accordance with the Beer-Lambert law.What we want to know then is the values of the function . What we can measure will be the values of the function .
When the attenuation conforms to the Beer-Lambert law, the relation between and is given by equation (1) and the absorbance along the light beam path is given by equation (2). Here is the intensity of the light beam before transmission, while is the beam intensity after transmission.
Here, the direction from the light source toward the screen is defined as the direction and that perpendicular to the direction and parallel with the screen is defined as the direction. (Both the and coordinate systems are chosen such that they are reflections of each other without mirror-reflective transformation.)[clarification needed]
By using a parallel beam irradiation optical system, one can experimentally obtain the series of fluoroscopic images (these are one-dimensional images of a specific cross section of a scanned object) for each angle between the object and the transmitted light beam. In Fig.3, the plane rotates counter clockwise.[Notes 1] around the point of origin in the plane in such a way “to keep mutual positional relationship between the light source (2) and screen (7) passing through the trajectory (5).” Rotation angle of this case is same as above-mentioned θ.
The beam having an angle is the collection of lines , represented by equation (3) below.
The function is defined by equation (4). That is equal to the line integral of μ(x,y) along of (eq. 3) as the same manner of (eq.2). This means that, of following (eq. 5) is a resultant of Radon transformation of μ(x,y).