Image Reconstruction from Fan-Beam Projections without Back-Projection Weight in a 2-D Dynamic CT : Compensation of Time-Dependent Rotational , Uniform Scaling and Translational Deformations

In a dynamic CT, the acquired projections are corrupted due to strong dynamic nature of the object, for example: lungs, heart etc. In this paper, we present fan-beam reconstruction algorithm without position-dependent backprojection weight which compensates for the time-dependent translational, uniform scaling and rotational deformations occurring in the object of interest during the data acquisition process. We shall also compare the computational cost of the proposed reconstruction algorithm with the existing one which has position-dependent weight. To accomplish the objective listed above, we first formulate admissibility conditions on deformations that is required to exactly reconstruct the object from acquired sequential deformed projections and then derive the reconstruction algorithm to compensate the above listed deformations satisfying the admissibility conditions. For this, 2-D time-dependent deformation model is incorporated in the fan-beam FBP reconstruction algorithm with no backprojection weight, assuming the motion parameters being known. Finally the proposed reconstruction algorithm is evaluated with the motion corrupted projection data simulated on the computer.


Introduction
The object of interest being imaged may be dynamic in nature, hence it is very important to develop algorithms and techniques to compensate the motion artefacts in tomographic imaging.There are several techniques and algorithms developed to tackle the challenges posed by the dynamic nature of the object such as heart, lungs, etc.One straightforward approach would be to minimize or limit the data acquisition time using ultra fast scaffolds.But in practice, dedicated softwares must be implemented for good results.Different techniques have been developed for compensation of motion artefacts, namely, 1) Gating method by Kachelriess et al. and Flohr et al. [1,2].2) Integrating motion models in reconstruction process using algebraic and analytical formulations [3][4][5][6].
3) Reconstruction method for general type of motion with moderate amplitude [7,8].4) Method to correct the sinogram through simple motion models as global translations, magnifications or rotations [9].5) Reconstruction formula for compensation of translation and magnification [8].
Roux et al. [10] have proposed a new frame-work, wherein the real object, at a particular time, is mapped to a reference object at reference time through some transformation, thus giving rise to a virtual acquisition geometry.The key point to observe is that, the value of the Radon transform of the reference object, in the virtual geometry, at reference time is linked to the Radon transform of the transformed object, at a particular time.Thus giving rise to the image reconstruction formula from de-formed projections; involving parameters of the virtual acquisition geometry but the aforementioned reconstruction formula is convoyed with a position-dependent backprojection weight.Since large amount of data needs be processed in dynamic studies, the implementation of the FBP formula for dynamic CT proposed by Roux et al. [10] will consume more time than it's counterpart (static case).The solution to eliminate the position-dependent back-projection weight in the fan-beam FBP formula, proposed by Noo.et al. [11] inspired by [12] in a static 2-D CT was given in [13].In this paper, we reduce the computational cost of the fan-beam FBP reconstruction formula for a 2-D dynamic CT proposed by Roux.et al., by eliminating the position-dependent back-projection weight, by choosing the weighting function similar to the one proposed in [13] for compensating time-dependent rotational, uniform scaling and translational deformations.By uniform scaling, we mean that the scaling transformation scales the object with the same magnitude along both the orthogonal reference co-ordinates.Thus for the above mentioned deformations, our formula will consume less time than the one proposed by Roux et al., due to the elimination of computationally expensive backprojection weight.
We follow the sequential approach listed below to address the issues raised in the paper.1) Admissibility condition to reconstruct the object from the projections acquired during it's deformation (affine).2) Reconstruction formula for a dynamic CT, compensating the time-dependent affine deformation.3) Incorporating the techniques in the above formula to eliminate back-projection weight.4) Simulating the idea presented in the paper.

Theory
Consider an object having strong dynamic nature and thus it deforms as a function of time during the data acquisition process.Let the object be denoted as t f at time and 0 t f at time 0 , where 0 is the reference time, where the object is assumed to be stationery.t t t  be the transformation that maps the real object at particular time to the reference time 0 t as shown in Figure 1.Therefore t transforms the object from "real state" to "reference state".
where and Thus, each real X-ray path on the evolving object is associated with a virtual X-ray path on the reference object, giving rise to a virtual acquisition geometry.Since affine transformations preserve intersection, each fan beam projection is transformed into another fan beam projection.
In a fan beam geometry, the parameter  uniquely determines one set of projection and also gives an understanding of the source location.Since  is varied with respect to time for acquiring all the projections, therefore can be replaced with   , g    be the dynamic fan-beam projection data collected with a circular scan trajectory given by where is the radius of the source trajectory.[14] to get a relation between the Hilbert Transform of the acquired fan-beam data and the Hilbert transform of the reference object projections.The result is given below: where OH is the Hilbert Filtered Parallel beam Projections of a ray perpendicular to n and the distance of s from the origin along n, thus Hilbert transform is in the direction of n whereas p H g  is the Hilbert filtered fanbeam projections along the direction T A  and n   0  a is the virtual source trajectory given by: . On modification of the condition stated by Noo et al. [11] (Noo's formula will be briefed in the next section) for static case, the following admissibility condition is deduced for the dynamic CT: Region of interest (ROI) in a dynamic object can be successfully reconstructed if all lines passing through the neighborhood of ROI intersects the virtual trajectory   0  a [10].

Fan-Beam Reconstruction Formula Based on Hilbert Filtering
A new class of algorithm based on Hilbert FBP (HFBP) has been derived [11] from the relation given in [14] for equi-angular static fan-beam projections, which is given by where   is the angle characterizing the ray that diverges from  0  a and contains , the projection data is characterized by  and are defined through parametrization with respect to n  and  cos sin and sin cos ,  e  (7)

 
, w   is the redundancy weight and where 0 is the point, to be reconstructed, on the reference object.
x Two unit orthogonal vectors u and given in the above equation are described as follows: weight in ramp FBP algorithm [15].
Noise non uniformity is appeared in the reconstructed image since ramp FBP algorithm is implemented using shift invariant filter [16].Noise uniformity is achieved with shift variant filter [17].Recently, two fan-beam algorithms [13,18] with no backprojection weight have been derived to address the position or distance dependent non-uniformity of resolution and noise variance by eliminating the backprojection weight.Extension of these two algorithms for cone-beam tomography has been given in [19].

Fan-Beam Reconstruction Formula with No Backprojection Weight
As discussed earlier, fan-beam FBP reconstruction formula ( 4) is convoyed with a position dependent weight, thus increasing the computational cost of the overall algorithm.The solution to eliminating this weight has been given in [13] by choosing a weighting function proportional to the back-projection weight x a (9) and such that x (11) This is the fan-beam FBP reconstruction algorithm with no back-projection weight for equi-angular curved detector geometry.

Incorporation of Motion Model in Static
Fan-Beam FBP Reconstruction Formula In this section, we shall brief the readers about the motion compensation algorithm given in [10] to tackle motion artefacts arising due to affine deformations during the data acquisition process.This formula was derived by incorporating motion model or deformation model into the fan-beam FBP reconstruction formula given in [11].
Since the static formula is convoyed with a position dependent back-projection weight, therefore the affine deformation compensation formula for dynamic CT is also computationally expensive due to the back-brojection weight inhereted from formula given in [11].We present a technique to eliminate this position dependent backprojection weight for translational, uniform scaling and rotational deformation by choosing the weighting func-tion similar to the one chosen in [13] for the static case.

Reconstruction Formula with Backprojection Weight: Compensation of Affine Deformation
Let   , g    be the deformed projection data.A fan- beam ray is parametrized by the parameters  and  , where  is the angle between a ray of interest and the line connecting source position to the origin.With reference to the dynamic object, the normal of a ray, passing through a point to be reconstructed, is denoted as x   ,    , where  is the angle of the ray connecting source and point to be reconstructed with in the fan: v .For each real x-ray, diverging from the vertex point, there exist a virtual x-ray path given by,

 
, s x a is the unit vector along the real x-ray path and x a is the unit vector along the virtual x-ray path.Consider as the normal to the virtual ray, joining the source and the point to be reconstructed on the reference object.
With Reference [10], we directly write the following relation, And the fan-beam reconstruction formula from deformed projection data is given by x a x (13) where , , where, the set and where H h is the impulse response of Hilbert Filter.For static fan-beam case, the virtual trajectory and the real trajectory is the same, similarly the reference object is the same as dynamic object.Thus , where is the Identity matrix.

Reconstruction Formula with No Backprojection Weight: Compensation for Translational, Uniform Scaling and Rotational Deformations
In this section, we use the weighting function (9), and eliminate the position dependent back-projection weight in the image reconstruction formula given by Equation (13).The weighting Equations ( 9), (10), the above Equation ( 19) is valid as it normalizes to one, the contribution from any line through within the FOV.The back projection weight in the virtual acquisition geometry is given by: (20) where Substituting the result obtained in Equation ( 21) to the back-projection weight we get, The key point to observe from Equation ( 22) is that the back-projection weight is independent of translational 0 such that, Open Access OJMI A. V. NARASIMHADHAN ET AL. 140 and rotational deformation, it only depends on scaling and shearing.The formula (13) compensates for time dependent affine deformation but is convoyed with a position dependent weight as discussed earlier.Therefore, if the object being imaged is subjected to translational or rotational deformation or both, then the term     a x in the weighting function (19) does not change (i.e.

 
  0 0 ), but scaling and shearing affects the back-projection weight as suggested by Equation (21).We propose to eliminate the position dependent back-projection weight for rotational, uniform scaling and translational deformations or their combinations.
Concretely, for a translational deformation, the matrix is chosen as and the  b can be chosen any vector since it does not effect the modification process of weighting term (22).In case of only rotational deformation, the transformation is done with a matrix given by therefore 1 So, back-projection weight does not change in case of rotational deformation.Now we shall discuss about uniform scaling wherein the back-projection weight changes.For uniform scaling deformation, therefore, Thus, by employing the weighting function (19) and the Equation (26) in Equation ( 13) we deduce the following fan-beam FBP reconstruction formula: where set is defined as: where  1, 2, 3, k   and c is the initial value of  .Equation ( 27) gives us the exact image reconstruction formula for compensating rotational and translational deformation during acquisition of equi-angular fan beam projections without any position dependent back-projection weight.

Salient Features of Proposed Algorithm
Since the technique involves derivative of the Hilbert filtered projection data, we have an advantage in terms of noise performance [20,21].The implementation of modified formula exhibits less variance compared to the fanbeam algorithm with no backprojection weight since more implementation blocks in modified formula.It can be mentioned here that the modified formula gives better computational efficiency than the unmodified for dynamic case.However, we shall give execution time form both algorithm in the simulation section.

Simulation Results
We simulate projection of low contrast Shep Logan phantom by considering the static as well as rotational and scaled deformed object.The scan parameters which we have used for simulation have been given Table 1.We have chosen A  or t A as diagonal matrix with following parameters is chosen as constant term with respect to time. Figure 2 shows the sinogram of original Shep Logan and deformed Shep Logan phantom.Firstly, the Shep Logan phantom is reconstructed using the fan-beam algorithm with no backprojection weight for a static object, which inturn yields an exact reconstruction of the Shep Logan Phantom (Figure 3(b)).Secondly, we reconstruct the Shep Logan Phantom from rotationally deformed data using fan-beam algorithm with no backprojection weight and our presented algorithm (Figures 3(c) and (d)).The observations completely abide with the theoretical and mathematical formulations and the Shep Logan phantom is exactly reconstructed.It is clearly observed from the Figure 3 that our proposed reconstruction algorithm removes the rotational motion artefacts.

Computational Complexity
The only difference between the Equations ( 27) and ( 13)   time than the reconstruction using (13).However, we have implemented proposed formula for dynamic projection using "C" complier.From the Table 2, it is clearly observed that computational efficiency is achieved significantly using modified fan-beam algorithm with no backprojection weight than existing exact fan-beam reconstruction for dynamic projections data.

Conclusion
In this paper, we ha struction algorithm compensating time-dependent rotational deformation.Since large amount of data is needed for dynamic CT, elimination of position-dependent back-projection weight from dynamic fan-beam reconstruction algorithm given by Roux et al., facilitates more advantages in terms of computational efficiency than the same in static CT.The presented algorithm inherits the better noise and resolution characteristics same as fan-beam reconstruction algorithm with no backprojection weight for static CT.It is very clear from the above simulation results that the proposed algorithm compensates for time-dependent rotational, translational and scaling deformation and exactly reconstructs the object, under consideration without any motion artefacts.

2  2 . 1 .
Admissibility Conditions on the Deformation for Fan-Beam Case Roux et al. have proposed a derivation of admissibility conditions based on generalization of Hamaker et al.
which is tangent to the source trajectory at  divergent beam reconstruction algorithms have spatially more uniform resolution and noise performance, and reduced divergent beam artefacts as a result of


Thus the value of radon transform is same for the rays characterized by x  .Therefore, by choosing the weight- ing function(9), the position dependent back-projection weight is eliminated, yielding the below given formula: