CC BY
59
DOI: 10.14529/ctcr180203
REVIEW OF IMAGE RECONSTRUCTION METHODS IN X-RAY COMPUTED TOMOGRAPHY WITH CONE-BEAM GEOMETRY
E.N. Simonov, [email protected],
M.V. Avramov, [email protected]
South Ural State University, Kyshtym, Russian Federation
The article reflects the main task of the X-ray computed tomography and its mathematical description. Direct and inverse Radon transform are given. The methods of image reconstruction in X-ray computed tomography are presented. Their brief classification are given. The author have reviewed classical problem of two-dimensional computed tomography and basics approaches to its solution. Emphasis is placed on back projection algorithm with filtering convolution. The derivation of the algorithm for parallel and fan-beam reconstruction are given. The analysis of the problem of three-dimensional reconstruction are presented. The author describes the additional conditions imposed on the projection data, the computational efficiency of the algorithms and the quality of the images. The basic trajectory of the X-ray source, providing the condition Tuy, are considered. The article gives an overview of existing methods of three-dimensional reconstruction with cone-beam geometry, their advantages, disadvantages, clinical applications. Their brief classification are given. Approximate algorithms of three-dimensional reconstruction are presented. The Feldkamp algorithm, the extended parallel backprojection, and the advanced single-slice rebinning are described. The author raises the question of developing approaches and methods for obtaining images with three-dimensional reconstruction for cone-beam spiral CT.
Keywords: computed tomography, three-dimensional image reconstruction.
Introduction
The main task of X-ray computed tomography (CT) is to reconstruct the internal structure of the objects of study by probing penetrating radiation. From a mathematical point of view, this problem can be considered as a problem of restoring some function in two-dimensional or three-dimensional space by known values of integrals along certain lines or planes. In CT, the desired function is a linear attenuation coefficient, and the projection data are known values of integrals.
The relationship between the projection data and the desired function of the linear attenuation coefficient is described by Radon transformation, which in the two-dimensional case is a linear integral of the function along the line. The representation of a function through its Radon transformation is called inversion formulas [1].
If the line L: x) = l is given in the normal form, the integral of the function f (x) along the line L is
+<»
P(6,l)= | f (lcos0-x'2sin9,lsin9 + x'2 cos9)dx'2, (1)
where
x'2 is the y-axis after the rotation by the angle 0.
The case formula will look like
I ^ 1
Iff 1 CP (9, l')
f ( x) =--2 ff---d9dl', (2)
W 4tc2 — 0 l — l' Cl'
-co 0
where l' is a curve passing through a specific point (xb x2) [1].
Thus, the main mathematical problem of X-ray computed tomography can be reduced to the integral equation of the form
T| = f |(x,y)dl = P(l,9). (3)
L(l ,9)
Equation (3), where T is the transformation operator, can be investigated by methods of integrated geometry of Radon. However, Radon methods allow us to obtain an exact solution of the basic KT equation provided that the initial data P(l, 0) are given precisely. But in practice, projective data
P(l, 9) contain different kinds of errors and are not exact. I mean ... it is not possible to accurately restore an image from P(l, 9) projections. In practice the regularization methods of equation (3) [2, 3] are applied.
All the variety of methods used for image restoration can be divided into two main groups: analytical and iterative methods [1, 2].
The implementation of any reconstruction algorithm depends on the X-ray generation scheme implemented in the tomograph. There are 3 basic schemes: parallel, fan, cone.
1. Reconstruction algorithm for parallel geometry of beams
The type of function ^(x, y), denoting the distribution of the desired physical quantity, a priori unknown. However, it is known that it is zero outside the reconstruction area Q, which is a circle of radius T with a center at the origin (Fig. 1).
Fig. 1. Parallel scan geometry
The position of the X-ray beam is determined by its distance l from the origin and angle 9. Projection defined by the formula (1) can be written in the form
+t
P(l,e)= J f (lcos0-1sin 0,lsin0 +1cos0)dt. (4)
—T
The next step in constructing the reconstruction algorithm is to find the inversion formula for (4). Moreover, as shown above, the expression (2) cannot be used directly.
One of the approaches of the derivation of the reconstruction approach is based on the use of the projection theorem. The first States that the operation on the projection at angle 9 is equivalent to a certain operation on the initial object ^(x, y). This theorem is called the generalized projection theorem. If the operation on ^(x, y) is reversible, such as the Fourier transform, then this implies a way of finding ^(x, y) given by P(l, 9). The second theorem, on the Central section, is formulated as follows. Let R(v, 9) is a one-dimensional Fourier transform of the function P(l, 9) according to the first argument
T
R(v,0)=J P (l,0)e—2nvldl, (5)
—T
where v = 1/2А/ is the spatial frequency; A! - the distance between the counts of the rays. Then it would be true equality
R (v, e) = ц (v cos 0, v sin 0), (6)
where ц is a two-dimensional Fourier transform of the function ц(х, y). In other words, the one-dimensional Fourier transform of projection data for variable / is equivalent to the two-dimensional Fourier transform of an object function expressed through polar coordinates.
It is obvious that the inverse two-dimensional Fourier transform from ц gives a reconstruction formula the function ц(х, y) on the set P(/, 0). However, in the practical implementation of this method of action there are significant difficulties associated with sampling errors and inaccuracy of projection data. The result is a very noisy image of the subject object's. To avoid this, it is necessary to take special measures [1, 2].
One of the ways solution to the problem is the use in (8) the Fourier transform of the function ц, and it will work on the window function W(v):
ц (v cos 0, v sin 0) = ц (v cos 0, v sin 0)W (v). (7)
Then, substituting (7) in (6), and taking the inverse Fourier transform, you can get
n T
ц = JJ P(/,0)g (/'- /)d/d0, (8)
0 -T
where
12 A!
g (/)= J |v|W (v) exp (i2nv!) dv, (9)
-1/2А/
/' = xcos0 + ysin0 - ray passing through a specific point P.
An internal integral in expression (8) is called a collapsed projection and is a convolution by / of projection P(/, 0) at the angle 0 and function g(!) given by expression (9). Function g(!) is called a collapsing function, and its appearance depends on the window function W(v) (Fourier function g (v) = |v|W (v)). The outer integral represents the inverse operation of projection. Formula (8), (9) form the basis of the method the reconstruction method of convolution and back-projection.
2. The Algorithm of reconstruction in fan geometry of the rays
To date, widespread scheme of scanning with a fan beam geometry, which is more preferable from the point of view of performance. Each beam is considered as a set of divergent rays, which defines the position of the radiation source - and у - the angle between the beam and the line connecting the source and the center of the circle (Fig. 2).
U is the position of the radiation source on a circle with radius D and center 0; в is the angle of inclination of the straight line 0U to the y-axis; 0 is the angle between the perpendicular 0R omitted from the origin of beam, and the x-axis; 0R equal to the length /; S - distance between points U and P; у is the angle between 0U and UP; r - length 0P; q(y, в) is a linear integral along the beam for fan geometry.
From formulas (7) to (9) for parallel geometry of the rays should
1 да T 2n
ц(rcosф,rsinф) = — J J J P(/,0)exp(i2nv[rcos(0-ф)-/])W(v)|v|d0d/dv.
2 -да-Т 0
Going from the variables (/, 0) variables (у, в), the operation of the reconstruction can be represented in two stages:
У m
q (у ',P)= J q (у,Р) cos у-g (S sin (у '-у)) dy; (10)
—ут
D 2?~
ц(x,y) = - J q (у',p)dp, (11)
20
where
00
g(Ssin(y'-y))= J \v\W(v)exp(i2%vSsin(y'-y))dv.
—00
If you accept y = y'— y,v' = vSsiny/y , you can write
g
(S sin у ) :
У
e2 • 2 ~ S sin у
J |v'exp(i2nv'у)dv'.
(12)
(13)
Fig. 2. Fan-shaped scanning geometry
The integral in (13) represents the kernel with the argument y obtained for parallel beams , and is the inverse Fourier transform of the spatial frequency module |v '|. The kernel for the fan-shaped geometry of the beams can be determined
■(У )■■
Y
-g
(Y).
o2 • 2
S sin y
As a result, formula (11) will take the form of a fan scan geometry
(14)
~ 1 1 r - -
H(*,y)- J dp J q (y, P)g (y '—y)dy, (15)
2 0 S —y„
where q (y,P) = q (y, P) D cos y is the modified projection;
y 2
g (y'-y) = g (y'-y)--is the modified kernel.
sin y
3. The reconstruction algorithm when the conical geometry of the beams
The above algorithms for parallel and fan geometry of rays formed the basis of the classical ap-
proach to solving the problem of computed tomography, in which measurements are made in the same plane with the help of linear detectors, and according to these data, the layer of the object is reconstructed.
ym
2
2
In practice, there are cases when the information in one slice of the studied object is not enough and it is necessary to understand the internal structure of the whole extended object. To obtain the necessary three-dimensional image, it is possible to carry out a series of measurements by the classical method and then combine the information obtained about the individual layers. This approach to solving a volume problem significantly increases the time of tomographic examination, which is not always acceptable in practice.
Another way to obtain a three-dimensional image is a three-dimensional reconstruction. In contrast to the classical two-dimensional CT method of measurement method of three-dimensional CT is produced in a cone of X-rays using a two-dimensional detector, which allows to obtain the projection of multiple layers in a single measurement. This approach provides some advantages over the classical research method: increases the photon flux, reduces the time of projections and makes it possible to increase their geometric.
In order to accurately restore the image along cone projections, it is necessary that each plane passing through the studied object at least once crosses the trajectory of the radiation source (Kirillov -Smith - Thuy condition [4, 5]). Examples of scanning schemes with the trajectory of the source that satisfy this condition: spiral, two circles, a circle and a straight line, two circles lying in mutually perpendicular planes [4].
Both algebraic and analytical methods were developed to solve the three-dimensional CT problem. These methods can be divided into two classes. The first one includes algorithms that divide a three-dimensional problem into a series of two-dimensional ones. From a set of cone projections by interpolation is obtained by a set of fan projections, on which a two-dimensional reconstruction of parallel layers of the object is performed. This approach is called the method of repackaging.
The second class includes algorithms that perform a complete three-dimensional reconstruction. Two-dimensional cone projections are filtered in some way, and then a reverse projection operation is performed. Such algorithms are often referred to as the FBP algorithms (filtered backprojection) or algorithms of the inverse projection with filtering [4, 5].
In [6] Feldkamp with co-authors proposed one of the first algorithms of three-dimensional reconstruction with cone geometry of rays. The article presents an algorithm for a source moving along a circle. This trajectory does not satisfy the data completeness condition because the planes parallel to the planes in which the source moves intersect the study object, but do not intersect the trajectory of the source. Therefore, accurate image recovery is not possible and the algorithm is approximate.
The Feldkamp algorithm belongs to the class of FBP algorithms. It performs a line-by-line convolution filtration, and the full geometry of the the inverse projection calculation.
Let a stationary object be specified for the study. The source-detector system revolves around the object, describing the circle. The plane containing the source is called media her plane. The intersection of the middle plane and the axis of rotation is taken as the origin. The distance from the source to the detector is D, from the source to the axis of rotation - d, from the axis of rotation to the detector - d' (Fig. 3).
With d ^ w, the algorithm passes into a layer-by-layer inversion formula for the parallel case. In addition, the algorithm gives the exact value of the density function in the Central section (z = 0) at any distance from the source.
For Fig. 4 /-rotation angle, Incoordinate along the detector that determines the recovery point.
cone scheme is taken into account only at the stage of
Fig. 3. The geometry of the scanning cone of rays
Fig. 4. The geometry at the mid-plane
Communication (l, 9) with (Y, Ф)
/ 1 /
l = Yd/ (d2 + Y2)2, Y = ldj (d2 -12)2, 9 = Ф + ^ + а.
(16)
The projection in cylindrical coordinates can be written as p (l, 0) = J rdr (f f ( r, $)8[r cos (0-^)-1 ] d f
Using the formula of inversion (3) in polar coordinates, when the point x = (r, and l = rcos(9 -and presenting the projection p (l, 9) through the Fourier transform (7), you can write the desired density function as
1
f (r,ф) = ^Re (f d0jvdv J p(l,0)eJvrM9-«-1 ]dl.
0
Performing the transition from the coordinates (l, 9) to (Y, O) by expressions (16) and from v to v' by
v' = v [d + r cos (<fr-O)]/(d 2 + Y 2 )2,
we can write down the final formula for the fan bundle in the median plane
f (r-ф)=f
d
2
[d + r cos (ф-Ф)]'
Pф [Y (r,ф)] dФ,
(17)
where
Y (r, ф) = dr sin (ф - Ф)/^ + r cos (ф - Ф)],
рф (Y И
d
ТРф (Y') g (Y - Y') dY'
g
-( d 2 + Y '2 ) 2
vy0
(Y) = Re J v exp ( jvY) dv.
(18)
(19)
Expression (17) is an inverse projection operation, (19) is a folding function, (18) is a rolled-up projection. Further, to determine the formula for three-dimensional reconstruction, it is necessary to determine the additional contribution of 5/ to the reconstructed density by the projection data for a small increment of the rotation angle 5®. From the projection data along the intersection of the plane of the detector, with the middle plane, we can calculate the contribution at points lying on the median plane. Projections that cross the plane of the detector along a line parallel to the median plane, but not lying in it,
themselves determine the plane. This plane is considered as an average plane-in a tilted coordinate system. It is necessary to correct the difference between the actual rotation of 5® along the vertical axis and an equivalent rotation of SO along the normal to the plane.
Further, the distance from the source to the detector in the inclined plane is substituted into the Radon's transformatuion. The density at the point r is taken equal to the sum of the additional contributions of all planes (one for each rotation angle) passing through r. Feldkamp's algorithm became one of the first algorithms of three-dimensional reconstruction, suitable for practical implementation. The algorithm gives an approximate solution for a circular trajectory, provides high quality of restoration, but does not take into account the cone angle (with the error of reconstruction grows with the angle), although it surpasses the classical approach and some iterative methods, it takes a lot of time for the computer account. In this connection, several modifications of the Feldkamp algorithm were proposed. The main Managing is given to the creation of high-speed algorithms and to the investigation of possible general the complexity of the method of collecting projection data [2, 5, 6].
So in [7] an algorithm was proposed, which was called the generalized parallel inverse projection (EPBP - extended parallel backprojection). From the classical Feldkamp algorithm, this the method is characterized by the use of repackaging and the transition to parallel data, which allows correctly take into account the redundancy of data, use an arbitrary pitch for spiral scanning, perform a reconstruction of the images of the heart.
An alternative to algorithms based on the Feldkamp method can be the single-layer repacking algorithm (ASSR - advanced single-slice rebinning). Term "One-layer repacking" implies the separation of the trajectory of the source into almost flat segments, each of which can be approximated by a flat curve. The oblique plane, holding such a curve, is called the reconstruction plane. Cone data transfom to these slanted planes that results in a set of fan projections. Behind the reconstructed flat tomograms are transformed into the required three-dimensional shape [2, 8].
Algorithms such as ASSR and EPBP give an approximate solution to the problem of three which is expressed in artifacts on the image. Development of an exact analytical algorithm ensuring optimal image quality is a complex task that has been the subject of intensive research in the world for several decades and has worthwhile time. The exact methods currently being developed are under investigation.
Conclusions
Algorithms of classical two-dimensional reconstruction of tomographic image for parallel and fan geometry of beams are considered. The analysis of the main approaches to the solution of the three-dimensional problem of KT is carried out, common approximate methods for reconstruction of three-dimensional images on a two-dimensional projection obtained when the conical geometry of the beams to planar and non-planar scan trajectories are considered.
The three-dimensional reconstruction of an image from cone beam projections is one of the most recent generalizations of CT. To obtain the exact solution it is necessary to use a flat trajectory that satisfies the completeness condition: every plane through the object under examination at least one time crosses the trajectory of the motion of the radiation source. Of all the promising trajectories providing accurate three-dimensional reconstruction, the most suitable for practical applications is a spiral covering the object both in radius and height. Spiral CT scan is a mode of volumetric tomography with non-planar geometry of radiation, which is performed by continuous longitudinal movement of the object of study.
The development of approaches and methods for obtaining tomographic images for three-dimensional reconstruction in the cone of rays for spiral CT is an urgent task of CT, which has become the subject of intensive research of the last few decades for the developers of reconstruction algorithms.
References
1. Simonov E.N. Fizika vizualizatsii izobrazheniy v rentgenovskoy komp'yuternoy tomogra/ii [Phisics Imaging in X-ray Computed Tomography]. Chelyabinsk, South Ural St. Univ. Publ., 2013. 479 p.
2. Kalender V. Komp'yuternaya tomogra/iya: osnovy, tekhnika, kachestvo izobrazheniy i oblasti kli-nicheskogo ispol'zovaniya. [Computed Tomography: Fundamentals, System Technology, Image Quality, Applications]. Moscow, Technosfera Publ., 2006. 344 p.
3. Tikhonov A.N., Arsenin V.Ya., Timonov A.A. Matematicheskie zadachi komp'yuternoy tomogra/ii [Mathematical Problems of Computed Tomography]. Moscow, Science Publ., 1987. 160 p.
4. Pikalov V.V., Lihachev A.V. [Comparison of Algorithms for Helical Tomography]. Numerical Methods and Programming, 2004, vol. 5, pp. 170-183. (in Russ.)
5. Likhachev A.V. [Comparison of Feldkamp Algorithm with Fourier Synthesis Algorithm for Three-Dimensional Imaging]. Avtometriya, 2006, no. 1, pp. 88-102. (in Russ.)
6. Feldkamp L.A., Davis L.C., Kress J.W. Practical Cone-Beam Algorithm. Journ. Opt. Soc. Amer. A., 1984, vol. 1, iss. 6, pp. 612-619. DOI: 10.1364/J0SAA.1.000612
7. Kachelrie M., Knaup M., Kalender W.A. Extended Parallel Backprojection for Standard Three-Dimensional and Phase-Correlated Four-Dimensional Axial and Spiral Cone-Beam CT with Arbitrary Pitch, Arbitrary Cone-Angle, and 100% Dose Usage. Medical Physics, June 2004, vol. 31, no. 6, pp. 1623-1641. DOI: 10.1118/1.1755569
8. Kachelrie M., Schaller S., Kalender W.A. Advanced Single-Slice Rebinning in Cone-Beam Spiral CT. Medical Physics, April 2000, vol. 27, no. 4, pp. 1033-1041.
Received 12 February 2018
УДК 519.6 DOI: 10.14529/ctcr180203
К ВОПРОСУ РАЗРАБОТКИ МЕТОДОВ РЕКОНСТРУКЦИИ ИЗОБРАЖЕНИЙ В РЕНТГЕНОВСКОЙ КОМПЬЮТЕРНОЙ ТОМОГРАФИИ С КОНУСНЫМ ПУЧКОМ ИЗЛУЧЕНИЯ
Е.Н. Симонов, М.В. Аврамов
Южно-Уральский государственный университет, филиал в г. Кыштыме, г. Кыштым, Россия
Отражена основная задача рентгеновской компьютерной томографии и ее математическое описание. Приведено прямое и обратное преобразование Радона. В статье представлены методы реконструкции изображений в рентгеновской компьютерной томографии. Дана их краткая классификация. Автором рассмотрена классическая задача двумерной компьютерной томографии, основные подходы к ее решению. Основное внимание уделено алгоритму обратного проецирования с фильтрацией сверткой. Приведен вывод алгоритма для параллельной и веерной схем сканирования. Проанализирована задача трехмерной реконструкции изображений. Описаны дополнительные условия, накладываемые на полноту проекционных данных, вычислительную эффективность алгоритмов, качество получаемых изображений. Рассмотрены основные траектории движения источника рентгеновского излучения, обеспечивающие выполнение условия Кириллова - Смита - Туя. Основное внимание уделено спиральной траектории. Приведен анализ существующих методов трехмерной реконструкции с конусной геометрией лучей, их достоинства, недостатки, области применения. Дана их краткая классификация. Приведены алгоритмы, дающие приближенное решение задачи трехмерной реконструкции. Описаны алгоритм Фельдкампа, алгоритм обобщенной обратной проекции, модернизированный алгоритм однослойной перепаковки. Затронут вопрос разработки подходов и методов получения изображений при трехмерной реконструкции в конусе лучей для спиральной компьютерной томографии.
Ключевые слова: компьютерная томография, трехмерная реконструкция.
Литература
1. Симонов, Е.Н. Физика визуализации изображений в рентгеновской компьютерной томографии /Е.Н. Симонов. - Челябинск: Издат. центр ЮУрГУ, 2013. - 479 с.
2. Календер, В. Компьютерная томография основы, техника, качество изображений и об-
ласти клинического использования / В. Календер; пер. с англ. А.В. Кирюшина, А.Е. Соловченко; под ред. В.Е. Синицына. - М.: Техносфера, 2006. - 344 с.
3. Тихонов, А.Н. Математические задачи компьютерной томографии / А.Н. Тихонов, В.Я. Ар-сенин, А.А. Тимонов. - М.: Наука, 1987. - 160 с.
4. Пикалов, В.В. Сравнение алгоритмов спиральной томографии /В.В. Пикалов, А.В. Лихачев // Вычислительные методы и программирование. - 2004. - Т. 5. - С. 170-183.
5. Лихачев, А.В. Сравнение алгоритма Фельдкампа с алгоритмом синтеза Фурье для трехмерной томографии /А.В. Лихачев //Автометрия. - 2006. - Т. 42, № 1 - С. 88-102.
6. Feldkamp, L.A. Practical cone-beam algorithm / L.A. Feldkamp, L.C. Davis, J.W. Kress // J. Opt. Soc. Am. A. - 1984. - Vol. 1, iss. 6. - P. 612-619. DOI: 10.1364/J0SAA.1.000612
7. Kachelrie, M. Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage /M. Kachelrie, M. Knaup, W.A. Kalender // Medical Physics. - 2004. - Vol. 31, no. 6. - P. 1623-1641. DOI: 10.1118/1.1755569
8. Kachelrie, M. Advanced single-slice rebinning in cone-beam spiral CT / M. Kachelrie, S. Schaller, W.A. Kalender //Medical Physics. - 2001. - Vol. 27, no. 4. - P. 1033-1041.
Симонов Евгений Николаевич, д-р техн. наук, профессор, заместитель директора по науке, Южно-Уральский государственный университет, филиал в г. Кыштыме, г. Кыштым; e.n.simonov@ yandex.ru.
Аврамов Максим Витальевич, аспирант кафедры экономики, управления и информационных технологий, Южно-Уральский государственный университет, филиал в г. Кыштыме, г. Кыш-тым; [email protected].
Поступила в редакцию 12 февраля 2018 г.
ОБРАЗЕЦ ЦИТИРОВАНИЯ
FOR CITATION
Simonov, E.N. Review of Image Reconstruction Methods in X-Ray Computed Tomography with Cone-Beam Geometry / E.N. Simonov, M.V. Avramov // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2018. - Т. 18, № 2. -С. 29-37. DOI: 10.14529/ctcr180203
Simonov E.N., Avramov M.V. Review of Image Reconstruction Methods in X-Ray Computed Tomography with Cone-Beam Geometry. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2018, vol. 18, no. 2, pp. 29-37. DOI: 10.14529/ctcr180203
