Planar tomosynthesis reconstruction in a parallel-beam framework via virtual object reconstruction

Size: px
Start display at page:

Download "Planar tomosynthesis reconstruction in a parallel-beam framework via virtual object reconstruction"

Transcription

1 Planar tomosynthesis reconstruction in a parallel-beam framework via virtual object reconstruction Brian E. Nett a,shuaileng a and Guang-Hong Chen a,b a Department of Medical Physics, University of Wisconsin-Madison, WI 53704; b Department of Radiology, University of Wisconsin-Madison, WI ABSTRACT A framework for image reconstruction from planar tomosynthesis trajectories (i.e. all source positions reside in a single plane) is presented. The parallel beam geometry is a convenient starting point in deriving reconstruction algorithms, both analytic and iterative, as the relation between frequency space and image space is well known. We present a method for utilizing parallel beam reconstruction algorithms in an internally consistent manner. The concept of a virtual image object is utilized. This virtual image object has the property that cone-beam projections through the real object are directly related to parallel-beam projections of the virtual object. The virtual object may then be reconstructed using any algorithm derived for parallel beam projections. Finally, an affine transform may be applied to the virtual image object in order to generate the reconstruction result. In the implementation described here the backprojection operation is performed such that the real image object is reconstructed without introducing a rebinning in image space. Image reconstruction comparisons are given for a standard filtered backprojection (FBP) type algorithm where parallel projections were assumed in the algorithm derivation. Numerical simulations were performed for a C-arm type geometry and parallel beam FBP reconstructions using the virtual object are compared with the standard backprojection algorithm. Finally, a comparison was made between the new parallel beam reconstruction and the standard approximation where the cone-beams are assumed to be approximately parallel beams and a cone-beam backprojection is employed. A reduction in streaking artifacts was observed using the new algorithm compared with the standard approximation. 1. INTRODUCTION In a tomosynthesis acquisition typically several cone-beam projections of the object are taken from different source and detector positions. In the literature on tomosynthesis parallel beam geometry is often assumed when deriving reconstruction algorithms. 1 4 These algorithms which have been derived in the parallel beam limit are then applied directly to the cone-beam case, without a strong theoretical foundation. There are significant advantages of assuming parallel beam geometry as the relationship between the projection data and the Fourier transform of the image object is well known. Thus, there is a motivation to utilize a parallel-beam geometry. Since it is not feasible to acquire parallel beam data from the tomosynthesis trajectories in practice. In this paper, we utilize the concept of a virtual object to convert cone-beam projections of a real object into parallel beam projections of a virtual object. Although methods have been proposed to populate frequency space 5 from divergent beam projections via summing many shifted local Fourier transforms, the projection from a single view no longer corresponds to a single localized slice in frequency space. In the framework presentedherewedemonstrate, for planar tomosynthetic trajectories, a local transform that preserves the one-to-one correspondence between a measured projection and a slice in frequency space, as this correspondence is crucial in algorithm development. This framework for tomosynthesis reconstruction may be utilized for any planar tomosynthesis acquisition (i.e. the source moves in a single plane during the acquisition). When the detector is oriented parallel to the source plane the conversion algorithm may be applied directly to the acquired projection data. Another possible acquisition for planar tomosynthesis is acquiring cone-beam projections from a circular tomosynthesis acquisition using an interventional C-arm system. 6 This type of acquisition may be of interest for providing interventional fluoroscopy and/or interventional perfusion measurements. In order to utilize the virtual object concept to perform the reconstruction rebinning of the projection data is utilized. After the data has been rebinned to a All scientific correspondence should be addressed to G.-H. Chen via gchen7@wisc.edu V. 7 (p.1 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

2 virtual detector that is parallel to the source trajectory the reconstruction may be carried out using the virtual object concept. Simulation results included here demonstrate: the concept of the virtual object, the dependence of the virtual object scaling on acquisition geometry, and the correspondence between cone-beam projections of the real object and parallel beam projections virtual object. Additionally, one of the reconstruction algorithms derived in the parallel beam framework is implemented here using the virtual object concept. The filtered backprojection based algorithm proposed by Pelc 1 that compensates for the density of frequency space sampling of the tomosynthesis trajectory was implemented here. The filter function in this algorithm was derived in a parallel beam framework. A comparison between the implementation using parallel projections through the virtual object and cone-beam projections through the real object is presented. A reduction in streaking artifacts was observed using the new algorithm compared with the standard approximation. 2. METHODS In this work we propose an internally consistent reconstruction algorithm for planar tomosynthesis (i.e. all source positions reside in a single plane that is parallel to the reconstructed image planes). This will be accomplished via reconstruction of a virtual image object which has the property that for our proposed acquisition geometry the cone-beam projections of the real image object are directly related to parallel beam projections of the virtual image object. In this framework the virtual image object need not be directly reconstructed as the conversion from real object to virtual object as well as the scaling operation may be taken account in the backprojection step. Note that the virtual object must be invariant as the x-ray source moves to different positions. In planar tomosynthesis the x-ray source is constrained to a given plane (e.g. linear tomosynthesis, circular tomosynthesis, etc.). Thus, the mapping relating the real object to the virtual object may be dependent on the distance between the source and detector plane, but must be independent of the source s position within the source plane. The affine transform which we have utilized in this work to map between the real and virtual object satisfies this criteria. Additionally, we note that in some acquisition geometries the projection data is acquired on a detector that is not parallel to the source-plane. In this case the projection data may be re-binned to a corresponding virtual detector that is parallel to the source plane. The virtual object framework may then be applied directly using this virtual detector. We stress that for planar tomosynthesis we will use a source position independent (within a given plane) mapping between the real object space and the virtual object space Three Dimensional Virtual Object The mapping from the real object to the virtual object is presented below. This mapping is derived so that the measurements on the detector of cone-beam projections may be directly related to parallel beam projections of a virtual object. The concept of treating divergent X-ray projections as parallel projections was previously proposed by Edhom and Danielsson. 7 Here we follow a similar formalism but introduce an additional parameter in our system parameterization and the virtual object definition utilized here is slightly different than that previously proposed. The direction chosen here for the parallel rays is parallel to the iso-ray which passes from the source point to the origin of the object space. The geometry utilized in this derivation is given in Figure 1. Where a given point source of divergent x-rays is denoted by S(x s,y s,r) and resides in a given plane a height R above the x y axis. In the 3D case the source is free to move along any planar trajectory. The measurements of the divergent rays are then made on a detector which is separated from the planar source by a distance D and lies parallel to the source plane. For a sample point in the object, P (x p,y p,z p ), the divergent ray that originates at the source point intersects the detector plane at point A(x A,y A, (R D)). The first step in generating a virtual object is to establish a mapping with one-to-one correspondence between the real image object and the virtual image object. The iso-ray is defined here as the ray which originates from the source point and passes through the origin of the real image space. The attenuation along ray SA provides the divergent beam measurement detected at a detector element centered at point A. Thismeasurementatpoint A may also be associated with the attenuation of the virtual object along a parallel ray. The direction of choice here is rays which are parallel to the iso-ray. In this case the ray VAis parallel to the ray S0 and thus we would like to relate the attenuation along the ray VAin the virtual object to the attenuation along SA. Therefore, we V. 7 (p.2 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

3 need a one-to-one mapping which relates each point along the ray SA to a point along the ray VA. For simplicity the mapping we choose here is to draw a line from the origin though the image point and the point at which this line intersects the virtual ray is the virtual point associated with this real point (i.e. P (x, y, z) V (x v,y v,z v )). Thus, in order to derive the relationship between the real and virtual object we must simply relate (x, y, z) to (x v, y v,z v ). This mapping is derived using the similar triangles of Figure 1. z z v = PP 0 VV 0 = OP OV = z v = 1+ PV z OP VA SO = VV D SS 0 VA SO = PV OP z v z 1 = z v + D R R z v = = z v + D R R OP OP + PV = 1 1+ PV OP D R z z (6) z v = Q z, where Q = D R z. (7) Based on our mapping choice it is clear that the scaling factor is identical for all three axes. Therefore, the virtual object coordinates are related to the real object coordinates as given below. x v = Q x, y v = Q y, z v = Q z, where Q = D R z. (8) The second step in establishing the virtual object concept is to account for the difference in projection value for the rays traversing the real object and the rays traversing the virtual object. Namely, we examine the path lengths in the real and virtual object spaces: Proj = dlf(x, y, z), and Proj v = dl f v (x v,y v,z v ), (9) PA where dl = VA dz, and dl = dz v cos θ 0 cos θ = DRdz (R z) 2 cos θ, where each geometric parameter is demonstrated in Figure 1, and the parameters of Proj and Proj v have been suppressed. After substituting the relationships for dl and dl back into the integral equation one obtains the following relationships dz Q 2 Rdz Proj = f(x, y, z), and Proj v = cos θ 0 D cos θ f v(x v,y v,z v ). (10) By examining the above equations one may make an informed choice of the virtual object, f v (x v,y v,z v )= 1 f(x, y, z). (11) Q2 (1) (2) (3) (4) (5) Given this choice of the virtual object we may substitute this definition back into our expression for a given virtual projection V. 7 (p.3 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

4 z S( x, y, R) S S V( xv, yv, zv) R 0 S0 ( xs, ys,0) y P( x, y, z) V0 ( xv, yv,0) D O P 0 ( x, y,0) x V ( x, y, R D) D v v (D-R) A( x, y, R D) A A (a) Figure 1. The geometrical configuration assumed here in deriving the virtual object mapping. Proj v = Rdz f(x, y, z) (12) D cos θ By comparing the expression for Proj in Eq. 10 to the expression for Proj v in Eq. 12 we obtain the relationship between the projections of the real object and the projections of the virtual object. Proj v = R D cos θ 0 Proj (13) cos θ We may then incorporate the definitions of the cos θ and cos θ 0 into this relationship, as provided below cos θ = R D and cos θ 0 = x 2 s + ys 2 + R 2 (xa x s ) 2 +(y A y s ) 2 + D, (14) 2 Proj v = x 2 s + y 2 s + R2 (xa x s ) 2 +(y A y s ) 2 + D 2 Proj. (15) This multiplicative factor that relates the parallel beam projections of the virtual object and the cone-beam projections of the real object will be referred to here as the virtual projection scaling Geometric Detector Description of C-arm Type System First, we will introduce the geometry of consideration here as shown in Figure 2. The angle t is measured with respect to the x axis and positive angles are measured from x in the direction of y. The angle θ is measured V. 7 (p.4 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

5 t z x R y D û SSD û vˆ ŵ vˆ ŵ Figure 2. Demonstration of the three separate coordinate systems which will be employed here. The geometrical parameters of interest in this geometry are: R the distance from the source plane to the origin, D the distance from the source plane to the virtual detector plane, û, ˆv, ŵ the unprimed local unit vectors of the real detector, û, ˆv, ŵ the primed local unit vectors of the virtual detector, ˆx, ŷ, ẑ are the lab unit vectors, SSD distance from the source to the detector surface, SOD distance from the source to the origin, t the view angle parameter, and θ the half tomo angle for this geometry. with respect to the z axis. In the cone-beam geometry the iso-ray originates from the focal point and passes through the origin and intersects the center of the detector. The source to surface distance of the detector will be referred to as the SSD, and likewise the source to origin distance is referred to as the SOD. These distances are measured along the iso-ray. When utilizing a standard interventional C-arm system each projection is taken in the unprimed coordinate system (Figure 2). The virtual object concept derived above was done under the condition that the detector plane is parallel to the source plane. Therefore, in order to apply the virtual projection scaling we introduce the primed coordinate system which is simply a translation of the world coordinate system (Figure 2). The relation between the unprimed and the primed coordinates is given below. We introduce several coordinate systems in order that we may easily derive associated relations to connect the primed and unprimed systems. The global three dimensional Cartesian coordinate system enables one to express a given position, r, in terms of its three components (x, y, z) r = xî + yĵ + zˆk, (16) where î, ĵ, ˆk are the unit vectors along the x, y, z axes respectively. Equally, well we may express r in terms of a coordinate system which rotates along with a detector that is perpendicular to the iso-ray (where the iso-ray originates at the focal point and passes through the origin). Note that in general these two coordinate systems do not share a common origin. Thus, we may represent r in this rotating coordinate system as, r = uû + vˆv + wŵ + O, (17) r = uû + vˆv + wŵ (SSD SOD)ŵ, (18) r = uû + vˆv +(w SSD + SOD)ŵ, (19) V. 7 (p.5 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

6 where if one stands upright at the source position directed toward the detector û points from right to left on the detector, ˆv points from bottom to top in the detector, and w points toward the detector from the focal point. In addition to this detector configuration we are also interested in the configuration in which the detector plane remains parallel to the x y axis so that the magnification is unaltered from view to view. Conveniently, the coordinate system located at the center of this detector is simply a translation of the global coordinate system: r = u û + v ˆv + w ŵ + O, (20) r = u î + v ĵ + w ˆk + O, (21) where the origin of this primed system is shared with the origin of the unprimed system. Now that we have defined our different coordinate systems we may use them to derive the quantities of interest. The unit vectors for the unprimed system are related to those of the global and primed systems as follows. sin t cos t 0 î û ˆv = cos t cos θ sin t cos θ sin θ ĵ (22) ŵ cos t sin θ sin t sin θ cos θ ˆk This relationship will be used directly to determine the parallel rebinning relationship from the primed system to the unprimed detector system. The parallel beam geometry here assumes rays which are parallel to w. If we have measurements of ray sums in the primed system (u,v,t,θ) and we would like to relate these to measurements in the unprimed system (u, v, t, θ) the following relation may be utilized: ( ( )( ) u sin t cos t u = v) cos t cos θ sin t cos θ v (23) Conversely, for parallel rays if we would like to determine a point in the primed detector from the unprimed detector we may use the following relationship: ( ) ( )( ) u sin t cos t v = cos θ u cos t sin t (24) cos θ v Note, that as θ goes to π/2 two of the matrix values in Eq. 24 go to infinity. However, when θ is π/2 the geometry corresponds to a cone-beam computed tomography acquisition, and thus in this limit we would use a cone-beam CT reconstruction formula rather than this planar tomosynthesis framework. 8, 9 Due to space constraints the cone-beam rebinning relationship between the two planes is not presented here but may be derived straightforwardly using this system geometry Reconstruction Steps In the case of planar tomosynthesis the acquisition will be of cone-beam projection data. Here we discuss applications where the X-ray source will remain in one given plane throughout the acquisition. Under these geometrical conditions we may design a unique virtual object which is an affine transformation of our original object. This conversion to a virtual object may be utilized with any of the reconstruction algorithms assuming parallel beams. 1 4 With the virtual object definition and the virtual projection scaling known we will outline the reconstruction steps below. The most basic implementation of the planar tomosynthesis reconstruction in the parallel beam framework is given in Figure 3 (a). This flow chart is provided for one to become familiar with this basic framework. This example is given for an acquisition where a flat panel detector lies parallel to the source trajectory. The virtual projection scaling may be applied directly on the acquired data in this case without the need for a re-binning operation (Eq. 14). One may then apply any filter derived assuming a parallel beam geometry. The data may then be backprojected via a parallel-beam backprojection in order to reconstruct the virtual object. After the virtual object is reconstructed one may re-grid from the virtual object space to the real object space using Eq. 8 and Eq V. 7 (p.6 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

7 Analytic Projections onto the unprimed Detector Cone-beam forward project onto a flat panel detector lying parallel to the source trajectory Cone-beam Re-binning from unprimed to primed detector Apply the virtual projection scaling Apply the virtual projection scaling Apply a filter if desired Parallel-beam backprojection to form a virtual object Re-grid from the virtual object to the real object Parallel-beam Re-binning from primed to unprimed detector Apply Parallel Beam Based Filter Parallel-beam backprojection from unprimed detector For each recon point Transform from real to virtual coordinates Find the backprojection location of the virtual point onto the detector Bi-linear interpolate among the filtered values Multiply the result by Q 2 end (a) (b) Figure 3. (a) Basic implementation steps. (b) Implementation steps for C-arm circular tomosynthesis. After describing the fundamental steps necessary for reconstruction in this framework we will present the required steps to perform reconstruction from the C-arm circular tomosynthesis geometry defined above. In this case the cone-beam projections are made directly onto the unprimed frame (real detector). A cone-beam rebinning maps the acquired cone-beam projections measured in the unprimed frame to the effective cone-beam projections which would have been measured in the primed frame (virtual detector). The virtual projection scaling may then be applied directly on the data in the virtual detector using Eq. 14. A parallel rebinning using Eq. 23 maps the scaled data to the real detector. On the real detector a parallel beam based filter is then applied. 1 3 Here we have implemented the kernel introduced by Pelc 1 : G θ (k u,k v )=( π sin θ ) k u W (k) (25) m θ where k u and k v are the conjugate variables of u and v respectively, k = k 2 u + k2 v, m θ is the number of projections and W (k) is a general window function added to control the tradeoff between resolution and the noise properties of the reconstruction. Finally, a parallel beam backprojection operation is utilized here. In this implementation we do not require an image space rebinning operation. The voxel driven backprojection algorithm utilized here directly reconstructs a point in the real object by: (1) finding the corresponding point in the virtual object (Eq. 8), (2) performing a bilinear interpolation along parallel rays, and (3) multiplying by Q 2 according to Eq. 11. Note this backprojection formalism via a virtual object has also been utilized in cone-beam CT reconstruction Parallel Beam Approximation The conventional approach which is employed in many tomosynthesis reconstruction algorithms is to assume parallel beams during the algorithm development and then approximately apply the derived algorithm to conebeam data. 2 4 The virtual object concept discussed here enables one to circumvent this approximation as parallel beams through the virtual object are directly generated in this framework. Thus, we will provide a comparison here using a filtered backprojection algorithm implemented both with the virtual object concept utilizing a parallel beam filter, and an approximation in which the parallel-beam filter is followed by a cone-beam backprojection V. 7 (p.7 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

8 2.5. Reconstruction Parameters Numerical simulations have been conducted assuming cone-beam measurements are taken in the geometry described in Figure 2. Analytic projection data was generated through two phantoms: a small uniform sphere placed at the origin and a high contrast Shepp-Logan phantom. The radius of the small sphere was 0.05, and the parameters for the ellipsoids which make up the high contrast Shepp-Logan phantom are given in the Appendix. The scanning geometry was parameterized as: detector Size = 4 4, detector Sampling = , R=5, D=8, θ =15, and view Sampling = 40/(2π). The reconstruction of the small sphere was conducted with an in-plane image matrix of and a total of 256 slices were reconstructed. The high contrast Shepp-Logan phantom was reconstructed with an in-plane image matrix of and a total of 128 slices were reconstructed. No roll-off window was applied in the filtering step. 3. RESULTS 3.1. Virtual Object Concept For demonstration purposes the transformation of a simple geometrical object is included here to graphically visualize this operation. A simple cube of uniform linear attenuation is the sample object, and Figure 4 (a) demonstrates one given slice through this cube. The same slice is shown in Figure 4 (b) after the object has been transformed to a virtual object. Following this operation we apply the inverse transformation to arrive at the result given in Figure 4(c). Thus, as demonstrated in (d) after applying the forward and inverse transform the sample object is well recovered, with the exception of a slight overshoot at the boundary caused by slight interpolation error. Therefore, we have confidence that images reconstructed in the virtual object space may be transformed to the real object space in a straightforward manner. It is important to note that this example of image space rebinning is given to provide insight into the virtual object concept and in the actual implementation no interpolation is performed in image space (Figure 3 (b)) (a) (b) (c) (d) Figure 4. The central x z slice of a centered uniform cube where each side has dimension 2 (a), the central x v z v slice of the virtual object (b), by applying the inverse transform to the image in (b) the central x z slice is recovered (c), a plot of the central vertical line through (a) and (c) is shown in (d) where the asterisk corresponds to the true value from (a) and the line corresponds to the recovered value from (c). After demonstrating the forward and inverse transform to the virtual object space we assess the effect of the magnification parameter on the virtual object. The phantom is demonstrated in Figure 5 (a). The phantom shown here is a cube with uniform attenuation value and a length of 0.5 on each side. The virtual object for three different magnification values will be investigated. By fixing D and scaling R the effective magnification of the system is varied here. It is clear that the ratio of R/D effects both the size and shape of the virtual object (Figure 5 (b-d)). As the ratio of R/D increases the shape of the virtual object more closely resembles the shape of the original object since the incident cone-beams are closer to the parallel beam limit. In the third set of numerical experiments we have demonstrated here the equivalence between the cone-beam projection of the original image object and the scaled parallel beam projection of the virtual image object Figure 6(a-c). The respective images of forward projections through the uniform cube using the standard Siddon ray tracing algorithm 11 are displayed in Figure 6(a) and Figure 6(b). A sample plot (Figure 6 (c)) demonstrates V. 7 (p.8 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

9 (a) (b) (c) (d) Figure 5. The x z cross section of the phantom (a) and three separate virtual objects for different locations of the source and detector plane (b-d): (b) R=1.1 and D=3, (c) R=1.5 and D=3, and (d) R=2 and D=3. the agreement between the projection values calculated for cone-beam projections of the original image object when compared with parallel beam projections of the virtual image object. Since we now have equivalent parallel projection measurements through the virtual object we can directly sample the Fourier space of the virtual object by utilizing the 3D projection slice theorem Ray Sum (a) (b) Projection Index (c) Figure 6. (a) The cone-beam projection of the original object where the focal point is (0,0,2.2). (b) Weighted parallel beam projection through the virtual object where the rays are parallel to the z axis. (c) A comparison plot of the projection onto the detector row v = 0 where the line corresponds to the cone-beam projection from (a) and the red circles to the parallel beam projection of the virtual object in (b) Reconstruction Results Images have been reconstructed here using the virtual object concept in a filtered backprojection framework proposed by Pelc. 1 Since this algorithm utilizes a filter kernel derived in the parallel beam framework it is a good candidate to demonstrate the virtual object framework. However, we note that any algorithm derived with the parallel beam assumption may be utilized for planar tomosynthesis reconstruction. The two phantoms reconstructed were a uniform sphere and the high contrast Shepp-Logan phantom. A relatively small half tomographic angle of 15 was employed in this study and thus significant signal leakage between planes is expected. The filtered backprojection algorithm utilizing the virtual object concept is first compared with the standard backprojection reconstruction. The comparison for the case of a uniform sphere is given in Figure 7. Next the central slice of the high contrast Shepp-Logan phantom was reconstructed (Figure 8). As noted by Pelc 1 we expect that this FBP result will accentuate the high frequency content of the image. The density of sampling in frequency space has been normalized in the sampled regions. However, the circular tomosynthesis sampling in frequency space leaves large regions of frequency space unsampled. Thus, this tomosynthetic sampling effectively V. 7 (p.9 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

10 populates more of the high frequency components in frequency space. Therefore, the high spatial frequency components have been accentuate in Figure 7 and Figure 8. (a) (b) Figure 7. Reconstruction results for a uniform sphere placed at the origin: (a) the axial slices and (b) the coronal slices. The respective images: phantom (left), backprojection (center), and parallel beam filtered backprojection with the virtual object (right). Figure 8. Reconstruction results for the central slice of the high contrast Shepp-Logan phantom. The respective images: phantom (left), backprojection (center), and parallel beam filtered backprojection with the virtual object (right). The second set of reconstruction comparisons provided here demonstrate the difference between the parallel beam filtered backprojection algorithm utilizing the virtual object concept and the approximation to the algorithm of using the parallel beam filter and cone-beam backprojection. The approximation has the advantage that no rebinning is required on the detector plane, while the parallel beam reconstruction with virtual object has the advantage that the filter and backprojection operations are completely consistent. The reconstruction comparison for the small sphere are given in Figure 9. Additionally, the high contrast Shepp-Logan phantom has been reconstructed with both of these approaches (Figure 10). One may note an increase in streaking artifacts V. 7 (p.10 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

11 from the high contrast objects in Figure 9 (b) and Figure 10 when the parallel-beam filter is combined with the cone-beam backprojection operation. Based on the simulation results the approximation of treating the conebeam projections approximately as parallel beam projections in filter derivation is well justified. However, the streak artifact pattern is reduced using our new fully consistent algorithm. Quantitative comparisons between these two approaches are not provided here as the spatial resolution has not yet been matched. (a) (b) Figure 9. Reconstruction results for a uniform sphere placed at the origin: (a) the axial slices and (b) the coronal slices. The respective images: parallel beam filtered backprojection with the virtual object (left) and parallel beam filter with cone-beam backprojection. Figure 10. Reconstruction results for the central slice of the high contrast Shepp-Logan phantom. The respective images: parallel beam filtered backprojection with the virtual object (left) and parallel beam filter with cone-beam backprojection. 4. CONCLUSIONS A planar tomosynthesis framework has been proposed using a virtual object concept in the reconstruction algorithm. This framework may be utilized with any existing tomosynthesis algorithm derived assuming parallel beams. Additionally, new avenues for new algorithm development have been paved such as frequency space based iterative algorithms. 12 Reconstruction simulations for an FBP type algorithm demonstrate reduction in streaking artifacts from high contrast objects, while confirming the previous assertion that the cone-beams could be treated approximately as parallel beams. Future research includes quantitative resolution and noise measurements. Acknowledgments We would like to acknowledge our funding support: NIH R01 EB , #T32CA , and the Herman I. Shapiro Graduate Fellowship. Thanks to Professor Norbert Pelc for providing our group with a copy of his PhD thesis. Additionally, thanks to Charles Mistretta for helpful discussions, and Orhan Unal for computer network support V. 7 (p.11 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

12 Appendix The definition for the high contrast Shepp-Logan phantom used in the simulations presented here are given below. (x,y,z) β a b c ρ (0,0,.25) (0,0,.25) (-.22,0,0) (.22,0,0) (0,0.35,0) (0,.1,0) (-.08,-.65,0) (.06,-.650,0) (.06,-.105,.875) (0,.1,.875) Table 1. The definition of the high contrast Shepp-Logan phantom where the conventions for the ellipsoids are taken from Kak and Slaney 13. REFERENCES 1. N. Pelc, A Generalized filtered backprojection algorithm for three dimensional reconstruction. PhD thesis, Harvard University, Boston, G. Lauritsch and W. H. Harer, A theoretical framework for filtered backprojection in tomosynthesis, Proc. SPIE 3338, pp , G. M. Stevens, R. Fahrig, and N. J. Pelc, Filtered backprojection for modifying the impulse response of circular tomosynthesis, Med. Phys. 28, pp , H. Matsuo, A. Iwata, I. Horiba, and N. Suzumura, Three-dimensional image reconstruction by digital tomo-synthesis using inverse filtering, IEEE Tran. Med. Imaging 12, pp , G.-H. Chen, S. Leng, and C. Mistretta, A novel extension of the parallel-beam projection-slice theorem to the divergent fan-beam and cone-beam projections, Med. Phys. 32, pp , B. Nett, J. Zambelli, C. Riddel, B. Belanger, and G.-H. Chen, Circular tomosynthesis implemented with a clinical interventional flat-panel based C-Arm, SPIE Medical Imaging , P. R. Edholm and P.-E. Danielsson, Divergent X-ray projections may under certain conditions be treated as parallel projections, Computer Methods and Programs in Biomedicine 57, pp , L. A. Feldkamp, L. C. Davis, and J. W. Kress, Practical cone-beam algorithm, J. Opt. Soc. Am. A 1, pp , B. E. Nett, T.-L. Zhuang, and G.-H. Chen, A cone-beam FBP reconstruction algorithm for short-scan and super-short scan source trajectories., Fully 3D Reconstruction in Radiology and Nuclear Medicine, J. Zambelli, B. Nett, S. Leng, C. Riddell, B. Belanger, and G.-H. Chen, Novel C-arm based cone-beam CT using a source trajectory of two concentric arcs, SPIE Medical Imaging , R. L. Siddon, Fast calculation of the exact radiological path for a three-dimensional CT array, Med. Phys. 12, pp , J. Velikina, S. Leng, and G.-H. Chen, Limited-view-angle tomographic image reconstruction via total variation minimization, SPIE Medical Imaging , A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging, IEEE, Bellingham, WA, pp. 102, V. 7 (p.12 of 12) / Color: No / Format: Letter / Date: 3/3/2007 8:44:49 AM

2D Fan Beam Reconstruction 3D Cone Beam Reconstruction

2D Fan Beam Reconstruction 3D Cone Beam Reconstruction 2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel

More information

Introduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University

Introduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University Introduction to Medical Imaging Cone-Beam CT Klaus Mueller Computer Science Department Stony Brook University Introduction Available cone-beam reconstruction methods: exact approximate algebraic Our discussion:

More information

Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset

Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,

More information

Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D.

Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced

More information

Feldkamp-type image reconstruction from equiangular data

Feldkamp-type image reconstruction from equiangular data Journal of X-Ray Science and Technology 9 (2001) 113 120 113 IOS Press Feldkamp-type image reconstruction from equiangular data Ben Wang a, Hong Liu b, Shiying Zhao c and Ge Wang d a Department of Elec.

More information

2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner

2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner 2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner Moscow-Bavarian Joint Advanced Student School 2006 March 19 2006 to March 29 2006 Overview 2D Fan Beam Reconstruction Shortscan Reconstruction

More information

TESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT

TESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT SCA2016-080 1/7 TESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT Alexander Katsevich 1, 2 and Michael Frenkel 1 1 itomography Corp., 2 University of Central Florida (UCF) This paper

More information

Novel C-arm based cone-beam CT using a source trajectory of two concentric arcs

Novel C-arm based cone-beam CT using a source trajectory of two concentric arcs Novel C-arm based cone-beam CT using a source trajectory of two concentric arcs Joseph Zambelli a, Brian E. Nett a,shuaileng a, Cyril Riddell c, Barry Belanger d,guang-hong Chen a,b a Department of Medical

More information

Improvement of Efficiency and Flexibility in Multi-slice Helical CT

Improvement of Efficiency and Flexibility in Multi-slice Helical CT J. Shanghai Jiaotong Univ. (Sci.), 2008, 13(4): 408 412 DOI: 10.1007/s12204-008-0408-x Improvement of Efficiency and Flexibility in Multi-slice Helical CT SUN Wen-wu 1 ( ), CHEN Si-ping 2 ( ), ZHUANG Tian-ge

More information

Background. Outline. Radiographic Tomosynthesis: Image Quality and Artifacts Reduction 1 / GE /

Background. Outline. Radiographic Tomosynthesis: Image Quality and Artifacts Reduction 1 / GE / Radiographic Tomosynthesis: Image Quality and Artifacts Reduction Baojun Li, Ph.D Department of Radiology Boston University Medical Center 2012 AAPM Annual Meeting Background Linear Trajectory Tomosynthesis

More information

Scaling Calibration in the ATRACT Algorithm

Scaling Calibration in the ATRACT Algorithm Scaling Calibration in the ATRACT Algorithm Yan Xia 1, Andreas Maier 1, Frank Dennerlein 2, Hannes G. Hofmann 1, Joachim Hornegger 1,3 1 Pattern Recognition Lab (LME), Friedrich-Alexander-University Erlangen-Nuremberg,

More information

DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE

DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan

More information

Central Slice Theorem

Central Slice Theorem Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial

More information

An approximate cone beam reconstruction algorithm for gantry-tilted CT

An approximate cone beam reconstruction algorithm for gantry-tilted CT An approximate cone beam reconstruction algorithm for gantry-tilted CT Ming Yan a, Cishen Zhang ab, Hongzhu Liang a a School of Electrical & Electronic Engineering, Nanyang Technological University, Singapore;

More information

Non-Stationary CT Image Noise Spectrum Analysis

Non-Stationary CT Image Noise Spectrum Analysis Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de

More information

Cover Page. The handle holds various files of this Leiden University dissertation

Cover Page. The handle   holds various files of this Leiden University dissertation Cover Page The handle http://hdl.handle.net/1887/8289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 2017-0-13

More information

Advanced Image Reconstruction Methods for Photoacoustic Tomography

Advanced Image Reconstruction Methods for Photoacoustic Tomography Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic

More information

THE FAN-BEAM scan for rapid data acquisition has

THE FAN-BEAM scan for rapid data acquisition has 190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 26, NO. 2, FEBRUARY 2007 Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full Partial Scans Jiangsheng You*, Member, IEEE, Gengsheng L. Zeng, Senior

More information

Implementation of a backprojection algorithm on CELL

Implementation of a backprojection algorithm on CELL Implementation of a backprojection algorithm on CELL Mario Koerner March 17, 2006 1 Introduction X-ray imaging is one of the most important imaging technologies in medical applications. It allows to look

More information

Approximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach

Approximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach Fundamenta Informaticae 135 (2014) 1 19 1 DOI 10.3233/FI-2014-1109 IOS Press Approximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach Linda Plantagie Centrum Wiskunde

More information

Image Reconstruction 3 Fully 3D Reconstruction

Image Reconstruction 3 Fully 3D Reconstruction Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction

More information

Translational Computed Tomography: A New Data Acquisition Scheme

Translational Computed Tomography: A New Data Acquisition Scheme 2nd International Symposium on NDT in Aerospace 2010 - We.1.A.3 Translational Computed Tomography: A New Data Acquisition Scheme Theobald FUCHS 1, Tobias SCHÖN 2, Randolf HANKE 3 1 Fraunhofer Development

More information

FDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW)

FDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW) Biomedical Imaging Volume 26, Article ID 83983, Pages 8 DOI.55/IJBI/26/83983 FDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW) Dong Yang and Ruola Ning

More information

An FDK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated

An FDK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated Home Search Collections Journals About Contact us My IOPscience An FK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated projections acquired using a circular trajectory This content

More information

Comparison of Probing Error in Dimensional Measurement by Means of 3D Computed Tomography with Circular and Helical Sampling

Comparison of Probing Error in Dimensional Measurement by Means of 3D Computed Tomography with Circular and Helical Sampling nd International Symposium on NDT in Aerospace - We..A. Comparison of Probing Error in Dimensional Measurement by Means of D Computed Tomography with Circular and Helical Sampling Jochen HILLER, Stefan

More information

Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction

Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Klaus Erhard a, Michael Grass a, Sebastian Hitziger b, Armin Iske b and Tim Nielsen a a Philips Research Europe Hamburg,

More information

An Iterative Approach to the Beam Hardening Correction in Cone Beam CT (Proceedings)

An Iterative Approach to the Beam Hardening Correction in Cone Beam CT (Proceedings) Marquette University e-publications@marquette Biomedical Engineering Faculty Research and Publications Engineering, College of 1-1-1999 An Iterative Approach to the Beam Hardening Correction in Cone Beam

More information

Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam

Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam Home Search Collections Journals About Contact us My IOPscience Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam projections This content has been downloaded from IOPscience.

More information

Medical Image Reconstruction Term II 2012 Topic 6: Tomography

Medical Image Reconstruction Term II 2012 Topic 6: Tomography Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross

More information

NIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.

NIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13. NIH Public Access Author Manuscript Published in final edited form as: Med Phys. 2008 February ; 35(2): 660 663. Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic

More information

Computed Tomography. Principles, Design, Artifacts, and Recent Advances. Jiang Hsieh THIRD EDITION. SPIE PRESS Bellingham, Washington USA

Computed Tomography. Principles, Design, Artifacts, and Recent Advances. Jiang Hsieh THIRD EDITION. SPIE PRESS Bellingham, Washington USA Computed Tomography Principles, Design, Artifacts, and Recent Advances THIRD EDITION Jiang Hsieh SPIE PRESS Bellingham, Washington USA Table of Contents Preface Nomenclature and Abbreviations xi xv 1 Introduction

More information

The n-pi-method for Helical Cone-Beam CT

The n-pi-method for Helical Cone-Beam CT 848 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 9, SEPTEMBER 2000 The n-pi-method for Helical Cone-Beam CT R. Proksa*, Th. Köhler, M. Grass, and J. Timmer Abstract A new class of acquisition schemes

More information

Material for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material

Material for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material Material for Chapter : Integral Equations in Visual Computing Material Basic Principles of Tomography c 00 Bernhard Burgeth 0 Source: Images Figure : Radon Transform: ttenuation http://en.wikimedia.org/wiki/image:radon_transform.png

More information

Acknowledgments and financial disclosure

Acknowledgments and financial disclosure AAPM 2012 Annual Meeting Digital breast tomosynthesis: basic understanding of physics principles James T. Dobbins III, Ph.D., FAAPM Director, Medical Physics Graduate Program Ravin Advanced Imaging Laboratories

More information

CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION

CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology

More information

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation

More information

A Novel Two-step Method for CT Reconstruction

A Novel Two-step Method for CT Reconstruction A Novel Two-step Method for CT Reconstruction Michael Felsberg Computer Vision Laboratory, Dept. EE, Linköping University, Sweden mfe@isy.liu.se Abstract. In this paper we address the parallel beam 2D

More information

Investigating Oblique Reconstructions with Super-Resolution in Digital Breast Tomosynthesis

Investigating Oblique Reconstructions with Super-Resolution in Digital Breast Tomosynthesis Investigating Oblique Reconstructions with Super-Resolution in Digital Breast Tomosynthesis Raymond J. Acciavatti, Stewart B. Mein, and Andrew D.A. Maidment University of Pennsylvania, Department of Radiology,

More information

A practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT

A practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT Journal of X-Ray Science and Technology 13 (2005) 9 21 9 IOS Press A practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT Hui Li a, Zhaotian Zhang

More information

Accelerated C-arm Reconstruction by Out-of-Projection Prediction

Accelerated C-arm Reconstruction by Out-of-Projection Prediction Accelerated C-arm Reconstruction by Out-of-Projection Prediction Hannes G. Hofmann, Benjamin Keck, Joachim Hornegger Pattern Recognition Lab, University Erlangen-Nuremberg hannes.hofmann@informatik.uni-erlangen.de

More information

Graphics and Interaction Transformation geometry and homogeneous coordinates

Graphics and Interaction Transformation geometry and homogeneous coordinates 433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation

More information

RECENTLY, biomedical imaging applications of

RECENTLY, biomedical imaging applications of 1190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 9, SEPTEMBER 2005 A General Exact Reconstruction for Cone-Beam CT via Backprojection-Filtration Yangbo Ye*, Shiying Zhao, Hengyong Yu, Ge Wang, Fellow,

More information

Radon Transform and Filtered Backprojection

Radon Transform and Filtered Backprojection Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department

More information

Suitability of a new alignment correction method for industrial CT

Suitability of a new alignment correction method for industrial CT Suitability of a new alignment correction method for industrial CT Matthias Elter 1, Nicole Maass 1, Peter Koch 2 1 Siemens AG, Healthcare Sector, Erlangen, Germany, e-mail: matthias.elter@siemens.com,

More information

X-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes

X-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.

More information

A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT

A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,

More information

Iterative CT Reconstruction Using Curvelet-Based Regularization

Iterative CT Reconstruction Using Curvelet-Based Regularization Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in

More information

MULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM

MULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM MULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM M. Simon, C. Sauerwein, I. Tiseanu, S. Burdairon Hans Wälischmiller GmbH Klingleweg 8, D-88709 Meersburg, Germany e-mail: ms@hwm.com ABSTRACT A new flexible

More information

AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT

AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT 1102 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 40, NO. 4, AUGUST 1993 AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT Grant. T. Gullberg and Gengsheng L. Zeng, Department of Radiology,

More information

Adaptive algebraic reconstruction technique

Adaptive algebraic reconstruction technique Adaptive algebraic reconstruction technique Wenkai Lua) Department of Automation, Key State Lab of Intelligent Technology and System, Tsinghua University, Beijing 10084, People s Republic of China Fang-Fang

More information

Basics of treatment planning II

Basics of treatment planning II Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Dose calculation algorithms! Correction based! Model based 1 Dose calculation algorithms!

More information

A Backprojection-Filtration Algorithm for Nonstandard. Spiral Cone-beam CT with an N-PI Window

A Backprojection-Filtration Algorithm for Nonstandard. Spiral Cone-beam CT with an N-PI Window A Backprojection-Filtration Algorithm for Nonstandard Spiral Cone-beam CT with an N-PI Window Hengyong Yu, Yangbo Ye,, Shiying Zhao, Ge Wang, CT/Micro-CT Laboratory, Department of Radiology, Department

More information

Aberrations in Holography

Aberrations in Holography Aberrations in Holography D Padiyar, J Padiyar 1070 Commerce St suite A, San Marcos, CA 92078 dinesh@triple-take.com joy@triple-take.com Abstract. The Seidel aberrations are described as they apply to

More information

Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013

Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013 2455-5 Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage 29 April - 3 May, 2013 Lecture NoteBasic principles of X-ray Computed Tomography Diego Dreossi Elettra, Trieste

More information

Image Reconstruction from Projection

Image Reconstruction from Projection Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital

More information

Answers to practice questions for Midterm 1

Answers to practice questions for Midterm 1 Answers to practice questions for Midterm Paul Hacking /5/9 (a The RREF (reduced row echelon form of the augmented matrix is So the system of linear equations has exactly one solution given by x =, y =,

More information

USING cone-beam geometry with pinhole collimation,

USING cone-beam geometry with pinhole collimation, IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 56, NO. 3, JUNE 2009 687 A Backprojection-Based Parameter Estimation Technique for Skew-Slit Collimation Jacob A. Piatt, Student Member, IEEE, and Gengsheng L.

More information

Modern CT system generations Measurement of attenuation

Modern CT system generations Measurement of attenuation CT reconstruction repetition & hints Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 4, 16 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering

More information

Reconstruction in CT and hints to the assignments

Reconstruction in CT and hints to the assignments Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 24, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging

More information

GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY

GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Engineering

More information

Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT

Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig

More information

GE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018

GE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018 GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB

More information

MATHEMATICS Grade 7 Advanced Standard: Number, Number Sense and Operations

MATHEMATICS Grade 7 Advanced Standard: Number, Number Sense and Operations Standard: Number, Number Sense and Operations Number and Number Systems A. Use scientific notation to express large numbers and numbers less than one. 1. Use scientific notation to express large numbers

More information

Cover Page. The handle holds various files of this Leiden University dissertation

Cover Page. The handle   holds various files of this Leiden University dissertation Cover Page The handle http://hdl.handle.net/887/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 207-04-3

More information

A prototype table-top inverse-geometry volumetric CT system

A prototype table-top inverse-geometry volumetric CT system A prototype table-top inverse-geometry volumetric CT system Taly Gilat Schmidt a Department of Radiology, Stanford University, Stanford, California 94305 Josh Star-Lack NexRay, Inc., Los Gatos, California

More information

Estimation of Missing Fan-Beam Projections using Frequency Consistency Conditions

Estimation of Missing Fan-Beam Projections using Frequency Consistency Conditions Estimation of Missing Fan-Beam Projections using Frequency Consistency Conditions Marcel Pohlmann, Martin Berger, Andreas Maier, Joachim Hornegger and Rebecca Fahrig Abstract Reducing radiation dose is

More information

Index. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting

Index. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms

More information

CHAPTER 4 RAY COMPUTATION. 4.1 Normal Computation

CHAPTER 4 RAY COMPUTATION. 4.1 Normal Computation CHAPTER 4 RAY COMPUTATION Ray computation is the second stage of the ray tracing procedure and is composed of two steps. First, the normal to the current wavefront is computed. Then the intersection of

More information

GPU-based cone-beam CT reconstruction with extended FOV

GPU-based cone-beam CT reconstruction with extended FOV GPU-based cone-beam CT reconstruction with extended FOV Tamás Huszár 1, Gábor Jakab 12, and Attila Rácz 1 1 Mediso Medical Equipment Developing and Service Ltd. Budapest, 1022 Hungary http://www.mediso.hu,

More information

Research Article Reconstruction Method for Nonconcyclic Dual-Source Circular Cone-Beam CT with a Large Field of View

Research Article Reconstruction Method for Nonconcyclic Dual-Source Circular Cone-Beam CT with a Large Field of View Mathematical Problems in Engineering Volume 2013, Article ID 526247, 7 pages http://dx.doi.org/10.1155/2013/526247 Research Article Reconstruction Method for Nonconcyclic Dual-Source Circular Cone-Beam

More information

Revisit of the Ramp Filter

Revisit of the Ramp Filter IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 62, NO. 1, FEBRUARY 2015 131 Revisit of the Ramp Filter Gengsheng L. Zeng, Fellow, IEEE Abstract An important part of the filtered backprojection (FBP) algorithm

More information

Biophysical Techniques (BPHS 4090/PHYS 5800)

Biophysical Techniques (BPHS 4090/PHYS 5800) Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University

More information

Artifact Mitigation in High Energy CT via Monte Carlo Simulation

Artifact Mitigation in High Energy CT via Monte Carlo Simulation PIERS ONLINE, VOL. 7, NO. 8, 11 791 Artifact Mitigation in High Energy CT via Monte Carlo Simulation Xuemin Jin and Robert Y. Levine Spectral Sciences, Inc., USA Abstract The high energy (< 15 MeV) incident

More information

CIRCULAR scanning trajectory has been widely used in

CIRCULAR scanning trajectory has been widely used in IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 7, JULY 2006 869 Region of Interest Reconstruction From Truncated Data in Circular Cone-Beam CT Lifeng Yu, Yu Zou, Emil Y. Sidky, Charles A. Pelizzari,

More information

8/2/2016. Measures the degradation/distortion of the acquired image (relative to an ideal image) using a quantitative figure-of-merit

8/2/2016. Measures the degradation/distortion of the acquired image (relative to an ideal image) using a quantitative figure-of-merit Ke Li Assistant Professor Department of Medical Physics and Department of Radiology School of Medicine and Public Health, University of Wisconsin-Madison This work is partially supported by an NIH Grant

More information

Enhancement Image Quality of CT Using Single Slice Spiral Technique

Enhancement Image Quality of CT Using Single Slice Spiral Technique Enhancement Image Quality of CT Using Single Slice Spiral Technique Doaa. N. Al Sheack 1 and Dr.Mohammed H. Ali Al Hayani 2 1 2 Electronic and Communications Engineering Department College of Engineering,

More information

3D X-ray Laminography with CMOS Image Sensor Using a Projection Method for Reconstruction of Arbitrary Cross-sectional Images

3D X-ray Laminography with CMOS Image Sensor Using a Projection Method for Reconstruction of Arbitrary Cross-sectional Images Ke Engineering Materials Vols. 270-273 (2004) pp. 192-197 online at http://www.scientific.net (2004) Trans Tech Publications, Switzerland Online available since 2004/08/15 Citation & Copright (to be inserted

More information

AP Physics: Curved Mirrors and Lenses

AP Physics: Curved Mirrors and Lenses The Ray Model of Light Light often travels in straight lines. We represent light using rays, which are straight lines emanating from an object. This is an idealization, but is very useful for geometric

More information

Reflection & Mirrors

Reflection & Mirrors Reflection & Mirrors Geometric Optics Using a Ray Approximation Light travels in a straight-line path in a homogeneous medium until it encounters a boundary between two different media A ray of light is

More information

Beam Attenuation Grid Based Scatter Correction Algorithm for. Cone Beam Volume CT

Beam Attenuation Grid Based Scatter Correction Algorithm for. Cone Beam Volume CT 11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic Beam Attenuation Grid Based Scatter Correction Algorithm for More Info at Open Access Database

More information

An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography

An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography 4th International Symposium on NDT in Aerospace 2012 - Tu.3.A.3 An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography Jonathan HESS *, Patrick KUEHNLEIN *, Steven OECKL

More information

Cover Page. The handle holds various files of this Leiden University dissertation

Cover Page. The handle   holds various files of this Leiden University dissertation Cover Page The handle http://hdl.handle.net/188/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 201-04-13

More information

INTRODUCTION TO MEDICAL IMAGING- 3D LOCALIZATION LAB MANUAL 1. Modifications for P551 Fall 2013 Medical Physics Laboratory

INTRODUCTION TO MEDICAL IMAGING- 3D LOCALIZATION LAB MANUAL 1. Modifications for P551 Fall 2013 Medical Physics Laboratory INTRODUCTION TO MEDICAL IMAGING- 3D LOCALIZATION LAB MANUAL 1 Modifications for P551 Fall 2013 Medical Physics Laboratory Introduction Following the introductory lab 0, this lab exercise the student through

More information

Reconstruction in CT and relation to other imaging modalities

Reconstruction in CT and relation to other imaging modalities Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound

More information

SPECIAL TECHNIQUES-II

SPECIAL TECHNIQUES-II SPECIAL TECHNIQUES-II Lecture 19: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Method of Images for a spherical conductor Example :A dipole near aconducting sphere The

More information

Theoretically-exact CT-reconstruction from experimental data

Theoretically-exact CT-reconstruction from experimental data Theoretically-exact CT-reconstruction from experimental data T Varslot, A Kingston, G Myers, A Sheppard Dept. Applied Mathematics Research School of Physics and Engineering Australian National University

More information

Maths Year 11 Mock Revision list

Maths Year 11 Mock Revision list Maths Year 11 Mock Revision list F = Foundation Tier = Foundation and igher Tier = igher Tier Number Tier Topic know and use the word integer and the equality and inequality symbols use fractions, decimals

More information

Coordinate Transformations in Advanced Calculus

Coordinate Transformations in Advanced Calculus Coordinate Transformations in Advanced Calculus by Sacha Nandlall T.A. for MATH 264, McGill University Email: sacha.nandlall@mail.mcgill.ca Website: http://www.resanova.com/teaching/calculus/ Fall 2006,

More information

A method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals

A method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals More Info at Open Access Database www.ndt.net/?id=15210 A method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals Subodh P S 1,a, Reghunathan

More information

Grangeat-type helical half-scan computerized tomography algorithm for reconstruction of a short object

Grangeat-type helical half-scan computerized tomography algorithm for reconstruction of a short object Grangeat-type helical half-scan computerized tomography algorithm for reconstruction of a short object Seung Wook Lee a) CT/Micro-CT Laboratory, Department of Radiology, University of Iowa, Iowa City,

More information

Using Algebraic Geometry to Study the Motions of a Robotic Arm

Using Algebraic Geometry to Study the Motions of a Robotic Arm Using Algebraic Geometry to Study the Motions of a Robotic Arm Addison T. Grant January 28, 206 Abstract In this study we summarize selected sections of David Cox, John Little, and Donal O Shea s Ideals,

More information

Reconstruction in CT and relation to other imaging modalities

Reconstruction in CT and relation to other imaging modalities Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound

More information

Content. Coordinate systems Orthographic projection. (Engineering Drawings)

Content. Coordinate systems Orthographic projection. (Engineering Drawings) Projection Views Content Coordinate systems Orthographic projection (Engineering Drawings) Graphical Coordinator Systems A coordinate system is needed to input, store and display model geometry and graphics.

More information

Digital Image Processing

Digital Image Processing Digital Image Processing Image Restoration and Reconstruction (Image Reconstruction from Projections) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering

More information

Medical Image Processing: Image Reconstruction and 3D Renderings

Medical Image Processing: Image Reconstruction and 3D Renderings Medical Image Processing: Image Reconstruction and 3D Renderings 김보형 서울대학교컴퓨터공학부 Computer Graphics and Image Processing Lab. 2011. 3. 23 1 Computer Graphics & Image Processing Computer Graphics : Create,

More information

Quality control phantoms and protocol for a tomography system

Quality control phantoms and protocol for a tomography system Quality control phantoms and protocol for a tomography system Lucía Franco 1 1 CT AIMEN, C/Relva 27A O Porriño Pontevedra, Spain, lfranco@aimen.es Abstract Tomography systems for non-destructive testing

More information

APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.

APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R. APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.China Abstract: When Industrial Computerized Tomography (CT)

More information

FRED Slit Diffraction Application Note

FRED Slit Diffraction Application Note FRED Slit Diffraction Application Note The classic problem of diffraction through a slit finds one of its chief applications in spectrometers. The wave nature of these phenomena can be modeled quite accurately

More information

Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction

Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction Hung Nien and Jeffrey A. Fessler Department of Electrical Engineering and Computer Science University of Michigan, Ann

More information

A sufficient condition for spiral cone beam long object imaging via backprojection

A sufficient condition for spiral cone beam long object imaging via backprojection A sfficient condition for spiral cone beam long object imaging via backprojection K. C. Tam Siemens Corporate Research, Inc., Princeton, NJ, USA Abstract The response of a point object in cone beam spiral

More information