Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Cone-beam mammo-computed tomography from data along two tilting arcs Kai Zeng,a兲 Hengyong Yu,b兲 Laurie L. Fajardo,c兲 and Ge Wangd兲 CT/Micro-CT Laboratory, Department of Radiology, University of Iowa, Iowa City, Iowa 52242 共Received 23 December 2005; revised 19 July 2006; accepted for publication 20 July 2006; published 13 September 2006兲 Over the past several years there has been an increasing interest in cone-beam computed tomography 共CT兲 for breast imaging. In this article, we propose a new scheme for theoretically exact cone-beam mammo-CT and develop a corresponding Katsevich-type reconstruction algorithm. In our scheme, cone-beam scans are performed along two tilting arcs to collect a sufficient amount of information for exact reconstruction. In our algorithm, cone-beam data are filtered in a shiftinvariant fashion and then weighted backprojected into the three-dimensional space for the final reconstruction. Our approach has several desirable features, including tolerance of axial data truncation, efficiency in sequential/parallel implementation, and accuracy for quantitative analysis. We also demonstrate the system performance and clinical utility of the proposed technique in numerical simulations. © 2006 American Association of Physicists in Medicine. 关DOI: 10.1118/1.2336510兴 Key words: cone-beam CT, mammography, exact reconstruction, Katsevich algorithm I. INTRODUCTION Breast cancer is ranked as the second leading cause of cancer death in women in the United States. It has been recognized that mass screening and early treatment are extremely important to reduce the mortality of breast cancer. Due to its specificity and sensitivity, x-ray mammography has been the method of choice for screening and diagnosis.1,2 However, x-ray mammography is far from being perfect because up to 17% of breast cancers are not identified with mammography, and normal breasts are associated with 70%–90% of mammograms suspicious of cancers.3 A major limitation of x-ray mammography is its projective nature, while the real anatomy and pathology is really in three dimensions 共3D兲. To address this problem, x-ray tomosynthesis and cone-beam computed tomography 共CT兲 are two compelling solutions. Tomosynthesis is a three-dimensional 共3D兲 imaging technique to reconstruct a series of images from a limited number of projections.4 Since its introduction in 1972, the area of tomosynthesis has been significantly advanced largely due to the development of the area detectors.5 A primary application of tomosynthesis is for breast imaging.6–8 The tomosynthetic algorithms are either analytic or iterative. The analytic algorithms are straightforward and efficient, such as self-masking,9 selective plane removal,10 and matrix inversion tomosynthesis.11 While the iterative algorithms are robust against noisy data and flexible to integrate prior knowledge, as it is done using algebraic reconstruction techniques,12,13 expectation-maximization,14 etc. none of these algorithms can avoid the inherent drawback of tomosynthesis due to the data incompleteness. Technically speaking, a breast volume should be imaged very well by cone-beam CT. Since more information of the object is acquired, the image quality of CT is much better than tomosynthesis, in terms of contrast resolution, geometrical distortion, etc. The concept of breast CT was proposed two decades ago,15 but little progress had been made initially because of compromised image quality and involved radia3621 Med. Phys. 33 „10…, October 2006 tion exposure. Again, thanks to the advancement in the digital detector technology, a number of groups investigated the feasibility and prototypes of cone-beam mammo-CT.16,17 Nevertheless, the algorithms for breast CT are still based on the traditional Feldkamp-type algorithms,18 and reconstruct images approximately with various artifacts. The fundamental classic results on exact cone-beam CT reconstruction were achieved by Grangeat,19 Smith,20 and Tuy.21 The recent breakthroughs on exact cone-beam CT algorithms were reviewed by Zhao et al.22 Up to now, there are a number of accurate and efficient cone-beam CT algorithms for various scanning trajectories, such as a helix,23–25 an arc-plus-line,26,27 a circle-plus-arc,28,29 and a saddle curve.30,31 Also, there are several algorithms which allow exact image reconstruction in the case of general trajectories.32–36 To improve image quality with cone-beam mammo-CT, we are motivated to design a cone-beam scanning mode that allows theoretically exact image reconstruction. In this article, we propose a novel scheme for cone-beam mammo-CT, which is theoretically exact, and develop a corresponding Katsevich-type reconstruction algorithm. In our scheme, cone-beam scans are performed along two tilting arcs to collect a sufficient amount of information for exact reconstruction. We derive our corresponding algorithm in the framework established by Katsevich,29,33 which may handle axial truncation of cone beam data. In what follows, we describe our system setup and derive the algorithm in Sec. II, describe numerical simulation results in Sec. III, and discuss relevant issues and conclude the article in Sec. IV. II. METHODS AND MATERIALS A. Cone-beam mammo-CT system In the proposed cone-beam mammo-CT system 共Fig. 1兲, a patient lays down on a table with one breast hanging through a hole. The x-ray tube and a flat-panel camera are fixed to a 0094-2405/2006/33„10…/3621/13/$23.00 © 2006 Am. Assoc. Phys. Med. 3621 3622 Zeng et al.: Cone-beam mammo-CT 3622 B. Cone-beam reconstruction method 1. Review of general cone-beam reconstruction formula First, we need to review Katsevich’s general cone beam reconstruction framework briefly.33 Let the scanning locus ⌳ be a finite union of smooth curves defined in R3: I ª 艛 关al,bl兴 → R3, l s 苸 I → y共s兲 苸 R3,兩ẏ共s兲兩 ⫽ 0 on I, 共1兲 where −⬁ ⬍ al ⬍ bl ⬍ ⬁ and ẏ共s兲 ª dy / ds. Assume f is smooth, compactly supported and identically equals zeros in a neighborhood of the locus ⌳, the cone beam transform of f along ⌳ is defined as FIG. 1. Proposed 3D mammography cone-beam CT system. 共a兲 System configuration and 共b兲 the cross section along the dashed line in 共a兲. rigid frame such as a C-arm to produce cone-beam projections. Then, a beast can be scanned twice along two tilting arcs respectively 关Figs. 1共b兲 and 2共a兲兴 by rotating the C-arm in different scanning planes. The other breast can be scanned in a similar fashion. Compared with the existing breast CT systems that only perform approximate reconstruction, our system for the first time aims at theoretically exact conebeam mammo-CT with fewer artifacts and better accuracy. Unlike other C-arm based algorithms,27,29 our scanning trajectory is more symmetric, which helps improve image quality in general. D f 共y, 兲 ª 冕 ⬁ f共y + t兲dt, y 苸 ⌳,  苸 S2 , 共2兲 0 where S2 is the unit sphere in R3. For given x 苸 R3 and 苸 R3 − 0, we also define 共s,x兲 ª x − y共s兲 , 兩x − y共s兲兩 x 苸 R3 \ ⌳, s 苸 I, ⌸共x, 兲 ª 兵z 苸 R3:共z − x兲 · = 0其. 共3兲 With y共s j兲 denoting points of intersection of ⌸共x , 兲 and ⌳ and ⬜ denoting the great circle 兵␣ 苸 S2 : ␣ ·  = 0其 which consists of unit vectors perpendicular to , we introduce the sets Crit共s,x兲 ª 兵␣ 苸 ⬜共s,x兲:⌸共x, ␣兲 is tangent to ⌳ or contains an end point of ⌳其 Ireg共x兲 ª 兵s 苸 I:Crit共s,x兲 傺 ⬜共s,x兲其 Crit共x兲 ª 艛 Crit共s,x兲. s苸I FIG. 2. Scanning geometry for exact cone-beam reconstruction. 共a兲 3D view of the tilting arcs and the breast and 共b兲 the exact reconstruction region U. Medical Physics, Vol. 33, No. 10, October 2006 共4兲 To reconstruct f at any fixed x 苸 R3, Katsevich suggested the following main assumptions about the locus ⌳.33 Property C1. (Tuy’s condition Ref. 21). Any plane through x intersects ⌳ at least at one point. Property C2. There exists a constant N1 such that the number of directions in Crit共s , x兲 does not exceed N1 for any s 苸 Ireg共x兲. Property C3. There exists a constant N2 such that the number of points in ⌸共x , ␣兲 艚 ⌳ does not exceed N2 for any ␣ 苸 S2 − Crit共x兲. If we only use a subset ⌳共x兲 傺 ⌳ for image reconstruction at x, ⌳ and I are replaced by ⌳共x兲 and I共x兲 for the earlier mentioned formulas, respectively. Let us consider a weighting function n共s , x , ␣兲, s 苸 Ireg共x兲, ␣ 苸 共⬜共s , x兲 − Crit共s , x兲兲. If n共s , x , ␣兲 is the piecewise constant and 兺 j:y共s j兲苸⌳艚⌸共x,␣兲n共s j , x , ␣兲 = 1 for almost all ␣ 苸 S2, we arrive at Katsevich’s general inversion formula33 as follows: 3623 Zeng et al.: Cone-beam mammo-CT f共x兲 = − ⫻ 1 42 冏冕 冕 3623 c 共s,x兲 m 兺 I共x兲 m 兩x − y共s兲兩 2 0 D f 共y共q兲,cos ␥共s,x兲 q + sin ␥␣⬜共s,x, m兲兲 冏 d␥ ds, sin ␥ q=s 共5兲 ␣⬜共s,x, 兲 ª 共s,x兲 ⫻ ␣共s,x, 兲, ␣共s,x, 兲 苸 ⬜共s,x兲, 共6兲 where is a polar angle in the plane perpendicular to 共s , x兲 and m 苸 关0 , 兲 are the points where 共s , x , 兲 are discontinuous, and cm共s , x兲 are values of the jumps cm共s,x兲 ª lim 共共s,x, m + 兲 − 共s,x, m − 兲兲, →0+ 共7兲 共s,x, 兲 ª sgn共␣ · ẏ共s兲兲n共s,x, ␣兲, ␣ = ␣共s,x, 兲 苸 ⬜共s,x兲. 共8兲 2. Existence of PI segments For our two-tilting-arcs scanning locus 关Fig. 2共a兲兴, we describe our geometry first, and then validate every property required by Katsevich’s general reconstruction framework. Let the three components of x 苸 R3 are x1 , x2 , x3, respectively. f共x兲 苸 C⬁0 is a smooth function inside the cylinder x21 + x22 ⬍ r2. Our scanning orbit ⌳ = Ca 艛 Cb defined on the real interval I = I1 艛 I2 consists two tilting arcs 关see Fig. 2共a兲兴 is parameterized by Ca ª 兵y共s兲:y 1 = R cos共s + smx兲cos t,y 2 = R sin共s + smx兲,y 3 = R cos共s + smx兲sin t,s 苸 I1其, Cb ª 兵y共s兲:y 1 = R cos共s − smx兲cos t,y 2 = R sin共s − smx兲,y 3 = − R cos共s − smx兲sin t,s 苸 I2其 共9兲 where I1 = 关− / 2 − 2smx , / 2兴, I2 = 共 / 2 , 3 / 2 + 2smx兴, R ⬎ r is the radius of the arcs, t is the tilting angle satisfying R cos t ⬎ r, and smx = arcsin共r / R cos t兲. smx is the fan angle. The cone beam projection data of f is defined as in 共Eq. 共2兲兲. As defined in Appendix A, the reconstruction region Ū can be decomposed into two kinds of reconstruction zones Ū1 and Ū2 关Figs. 2共b兲 and 14兴. There are exactly two PI segments within Ū2 and there is only one PI segment within Ū1. Here PI means “,” and a PI segment of x is the line segment containing x and its two end points on the scanning orbit ⌳. Our problem is to reconstruct f inside a region U which is the intersection of Ū and the cylindrical support 关Fig. 2共b兲兴. In practical clinical applications, the region in the chest should be avoided and only the breast part is concerned. Based on the definition of Ū, for any fixed x 苸 Ū, there exist at least one and at most two PI segments. For each PI segment, one end point must be on arc A while the other Medical Physics, Vol. 33, No. 10, October 2006 FIG. 3. Illustration of Pl lines for 共a兲 x 苸 Ū2 and 共b兲 x 苸 Ū1. must be on arc B if we ignore those points on the scanning arcs’ planes of zero measure. We denote the corresponding angular variable as sa = sa共x兲 for arc A and sb = sb共x兲 for arc B. The existence and uniqueness of these PI segments can be divided in the two cases: 共1兲 R sin共sa + smx兲 ⬎ x2 ⬎ R sin共sb − smx兲 and 共2兲 R sin共sa + smx兲 ⬍ x2 ⬍ R sin共sb − smx兲. Moreover, for any x 苸 Ū2 关Fig. 3共a兲兴, there exist exactly two PI segments with our symmetric imaging geometry. While, for any x 苸 Ū1 关Figs. 3共b兲 and 14共d兲兴, there exists exactly one PI segment, either in case 1 or case 2. With our scanning geometry, for a pair of given sa共x兲 and sb共x兲 for x 苸 Ū2 关Fig. 3共a兲兴, the PI segment may correspond to two PI interand vals Ī1共x兲 = 关− / 2 − smx , sa共x兲兴 艛 关sb共x兲 , 3 / 2 + smx兴 mx mx Ī2共x兲 = 关sa共x兲 , / 2 − s 兴 艛 共 / 2 + s , sb共x兲兴 yielding two PI arcs ⌳1 = ⌳1共x兲 = ⌳共Ī1共x兲兲 and ⌳2 = ⌳2共x兲 = ⌳共Ī2共x兲兲 关Fig. 3共a兲兴. 3. Reconstruction method To derive a theoretically exact image reconstruction algorithm in Katsevich’s framework,33 we need to check the properties of the locus ⌳ and the weighting function in our TABLE I. Weighting functions n共s , x , ␣兲 used for our Katsevich-type reconstruction. Case Value of n共s , x , ␣兲 1 IP, s1 苸 I1 1 IP, s1 苸 I2 3 IPs, s1 , s2 苸 I1 ; s3 苸 I2 1 1 n共s1 , x , ␣兲 = n共s2 , x , ␣兲 = 1, n共s3 , x , ␣兲 = −1 n共s1 , x , ␣兲 = n共s3 , x , ␣兲 = 1, n共s2 , x , ␣兲 = −1 3 IPs, s1 苸 I1 ; s2 , s3 苸 I2 , s2 ⬍ s3; x is above the plane 关in Ū2, Fig. 2共b兲兴 3 IPs, s1 苸 I1 ; s2 , s3 苸 I2 , s2 ⬍ s3; x is below the plane 关in Ū1, Fig. 2共b兲兴 n共s1 , x , ␣兲 = n共s2 , x , ␣兲 = 1, n共s3 , x , ␣兲 = −1 3624 Zeng et al.: Cone-beam mammo-CT 3624 particular geometry. Let us only consider the PI-arcs ⌳1 to reconstruct x 苸 Ū2 in the following, while ⌳2 can be similarly utilized. As there exists at least one PI segment for any x 苸 Ū, property C1, Tuy’s condition, is satisfied. Also, for y共s兲 苸 ⌳1共x兲 there are a limited number of planes ⌸共x , ␣兲, ␣ 苸 ⬜共s , x兲, which satisfies ⌸共x , ␣兲 contains an endpoint of ⌳1 or is tangent to ⌳1. Thus property C2 is satisfied. Moreover, apart from a set of zero measure, any plane ⌸共x , ␣兲, ␣ 苸 S2, through x intersects ⌳1共x兲 at an odd number of intersection points 共IP兲. And the number of IP is either one or three which validates the property C3. Overall, ⌳1 satisfies the three properties C1 – C3 required in Ref. 33. Following the work by Katsevich,29 we define a similar weighting function n共s , x , ␣兲 in Table I, where s j, j = 1 , 2 , 3 are the angular parameters satisfying y共s j兲 苸 共⌳1 艚 ⌸共x , ␣兲兲. Thus, we have 兺 j:y共s j兲苸⌳1艚⌸共x,␣兲n共s j , x , ␣兲 = 1 for almost all ␣ 苸 S2, and n共s , x , ␣兲 is the piece-wise constant. Therefore, our scanning arcs and weighting function are qualified to be fitted into Katsevich’s general reconstruction framework. Then, we need to find the discontinuities in of 共2.8兲 for each s 苸 Ī1共x兲, and determine the filtering directions. Similar to Katsevich’s treatment,29 we can analyze the discontinuities by projecting the source trajectory ⌳ = Ca 艛 Cb onto the detector plane DP共s0兲 共Fig. 4兲, which is defined as the plane through the origin O and perpendicular to the line through y共s0兲 and O. In this article, all the variables are denoted with hat signify the objects on a projection plane. The detector coordinates are defined along the directional vector u and v. For the sources on the arc A, that is, s0 苸 关−共 / 2 + 2smx兲 , / 2兴, u = 共sin共s0 + smx兲cos t,− cos共s0 + smx兲,sin共s0 + smx兲sin t兲, v = 共− sin t,0,cos t兲. 共10兲 And for the source on the arc B, that is, s0 苸 关 / 2 , 共3 / 2 + 2smx兲兴 FIG. 4. Projections of a Pl line on a detector plane for the configuration of Fig. 3共a兲. 共a兲 The projected trajectory with the source on Ca, 共b兲 the projected trajectory with the source on Cb, and 共c兲 an impossible trajectory projection. u = 共sin共s0 + smx兲cos t,− cos共s0 − smx兲,− sin共s0 − smx兲sin t兲, v = 共sin t,0,cos t兲. With the source on Ca, we can describe the projection of Cb on the detector plane as u共s兲 = R sin共s0 + smx兲cos共s − smx兲cos 2t − cos共s0 + smx兲sin共s − smx兲 , cos共s − smx兲cos共2t兲cos共s0 + smx兲 + sin共s − smx兲sin共s0 + smx兲 − 1 v共s兲 = R cos共s − smx兲sin 2t , cos共s − s 兲cos共2t兲cos共s0 + smx兲 + sin共s − smx兲sin共s0 + smx兲 − 1 mx where s0 苸 关− / 2 − 2smx , / 2兴 indicates the source position on the arc A, s 苸 共 / 2 , 3 / 2 + 2smx兲 and u共s兲 and v共s兲 are the detector coordinates. With the source on Cb, by geometric symmetry we can immediately obtain the same formula except for a changed range of s. When the source is on Ca 关Fig. 4共a兲兴, the discontinuities may occur at polar angle 1, 2, and 3 which corresponding to the plane P1, P2, and P3. It is easy to compute that the magnitudes of the jumps at these three Medical Physics, Vol. 33, No. 10, October 2006 共11兲 共12兲 positions are 2, 0, and 0, respectively. Thus, we only have a jump along direction parameterized by 1, which is the same filtering direction used in the Feldkamp algorithm.18 When the source is on Cb 关Fig. 4共b兲兴, the projection of a PI segment can only have one IP with the projection of ⌳1共x兲, giving the tangential direction 1 with a jump of magnitude 2. The last case 关Fig. 4共c兲兴 is simply impossible, as shown in Appendix B. Thus, the filtering directions should be along horizontal 3625 Zeng et al.: Cone-beam mammo-CT 3625 TABLE II. Linear attenuation coefficients at 38 keV. Material Attenuation coefficient 共cm−1兲 Skin Breast Calcifications Mass Fibrous 0.2963 0.2465 3.2822 0.2667 0.2667 1 22 f¯I2共x兲 = − ⫻ 冕 2 0 冕 兺 ¯I 共x兲 m=1,2 2 1 兩x − y共s兲兩 冏 D f 共y共q兲,cos ␥共s,x兲 q + sin ␥␣⬜共s,x, m兲兲 冏 d␥ ds. sin ␥ q=s 共14兲 Combining these two formulas, we have f共x兲 = 共f¯I1共x兲 + f¯I2共x兲兲/2. 共15兲 Consequently, our reconstruction algorithm can be divided into four steps: FIG. 5. Uncompressed breast phantom which contains masses, calcifications, and fibroses of different sizes and densities. lines 关Fig. 4共a兲兴 or tangential lines 关Fig. 4共b兲兴, depending on m determined based on our above analysis. This is exact the basic same result obtained by Katsevich.29 Therefore, the reconstruction formula can be expressed as f¯I1共x兲 = − ⫻ 1 22 冕 2 0 冕 1 兺 ¯I 共x兲 m=1,2 兩x − y共s兲兩 1 冏 D f 共y共q兲,cos ␥共s,x兲 q + sin ␥␣⬜共s,x, m兲兲 冏 共1兲 differentiate the projection data collected along PI-arc ⌳1共x兲 to obtain / qD f ; 共2兲 filter every projection along the direction defined by m. Projection data on arc A are horizontally filtered, and projection data on arc B are tangentially filtered; 共3兲 backproject the filtered data to reconstruct an image; and 共4兲 repeat the above three steps for the projection data collected along the other PI-arc ⌳2共x兲, and average the two images for the final reconstruction. For x 苸 Ū1 关Fig. 3共b兲兴, we can derive the formula and steps in a similar fashion as that for x 苸 Ū2 关Fig. 3共a兲兴. Since our reconstruction formula is derived in Katsevich’s framework,29,33, it can be regarded as a variant of the conebeam reconstruction algorithm in the case of circle-and-arc.29 C. Uncompressed breast phantom d␥ ds. sin ␥ q=s 共13兲 On the other hand, we can also reconstruct x from ⌳2共x兲. Similarly, we have A 3D mathematical mammography phantom 共Fig. 5兲 was designed in reference to several commercially available mammography phantoms, such as the American College of Radiology 共ACR兲 accreditation phantom, Mammography Imaging Screening Trial phantom and Uniform phantom.37,38 TABLE III. Dimensions and linear attenuation coefficients of the fibroses. Fibrous No. 1 2 3 4 5 6 7 8 9 10 Diameter 共mm兲 Height 共mm兲 AC 共cm−1兲 1.56 10 0.267 1.12 10 0.267 0.89 10 0.267 0.75 10 0.267 0.54 10 0.267 0.4 10 0.267 1.12 10 0.267 1.12 10 0.253 1.12 10 0.249 1.12 10 0.248 Medical Physics, Vol. 33, No. 10, October 2006 3626 Zeng et al.: Cone-beam mammo-CT 3626 TABLE IV. Dimensions and linear attenuation coefficients of masses. Mass No. 1 2 3 4 5 6 7 8 9 10 Diameter 共mm兲 AC 共cm−1兲 8.00 0.267 6.00 0.267 3.00 0.267 2.00 0.267 1.00 0.267 0.5 0.267 3.00 0.267 3.00 0.253 3.00 0.249 3.00 0.248 The breast was modeled as half an ellipsoid. While the existing phantoms contain important structures simulating mass, fibrous and calcification, they mimic a compressed breast for x-ray mammography, being suboptimal for our 3D CT simulation. To address this limitation, our phantom targeted an uncompressed breast with representative anatomical and pathological features, including structures of different sizes and contrasts. A detailed description of the phantom is as follows. The three ellipsoidal semiaxes were set to 50, 50, and 100 mm. The skin thickness was set to 2.5 mm. While the fibroses were modeled as cylinders, the calcifications and mass as balls 共Fig. 5兲. The phantom was positioned in the nonnegative space, attached to the chest wall defined on z = 0 mm. Fibroses were placed on the planes of z = 22.5 mm and z = 62.5 mm. Masses were centered on the planes of z = 35 mm and z = 75 mm. Calcifications were scattered on the plane of z = 42.5 mm. The linear attenuation coefficients 共AC兲 of the breast structures were specified for 38 keV x-rays16,39 as listed in Table II. The dimensions of the features in the phantom were specified in reference to those in the existing mammography phantoms, especially the ACR phantom. Tables III–V enlist the sizes and attenuation properties of the fibroses in the planes z = 22.5 and 62.5 mm, the parameters of the masses in the planes z = 35.0 and 75.0 mm, and the characteristics of the calcifications at z = 42.5 mm, respectively. III. NUMERICAL SIMULATION To corroborate the correctness of our proposed algorithm and demonstrate its clinical utility, we implemented it in C⫹⫹ for numerical tests. These results were also compared with a circular cone-beam CT scan using the FDK algorithm. The system setup and geometric parameters were summarized in Tables VI and VII, respectively. In our simulation, the reconstruction region U was delimited by the chest wall 关Fig. 6共a兲兴. The FDK algorithm can only approximately reconstruct the volume above the scanning plane 关Fig. 6共b兲兴, while our proposed method can exactly reconstruct a larger volume. In all the cases, the reconstruction results obtained using our algorithm were excellent except for the streak artifacts primarily caused by the high-density calcifications 共ten times higher than the background兲 共Figs. 6–8兲, which were also observed in the FDK results. The other structures, such as fibrous and mass, were reconstructed very well. Although the detector element size was comparable or even larger than the fine structures in our phantom, these features were clearly visible in the reconstructed images. As far as the calcifications are concerned, despite that they were substantially smaller than one pixel, they were revealed in the reconstruction. The low-contrast structures were well revealed 共Figs. 7 and 8兲 using our algorithm. Our algorithm generally produced more accurate results than the FDK algorithm, because there were significant density drops in the FDK results, especially far away from the central plane 共Figs. 7–10兲. Also, due to the approximate nature of the FDK algorithm, significant streak artifacts were observed 共Figs. 7 and 8兲 in the FDK results. Besides, our algorithm preserved the shape of the breast phantom, but there were substantial shape distortions in the FDK results especially near the top of the phantom 共Figs. 7 and 8兲. A number of numerical tests with noisy projection data were performed to simulate the practical imaging process. As the image quality is closely related to radiation dose, which is proportion to the total number of involved x-ray photons given other conditions being equal, the same total number of photons was used in each test for image quality comparison. In our study, the same radiation dose was uniformly distributed to every detector cell in each projection view for simplicity. Then, the Poisson noise was added to the projection data. The image quality difference between our algorithm and the FDK algorithm was similar to that in the noise-free case, except that in the noisy case our algorithm showed a better noise tolerance than the FDK algorithm.40 Specifically, the images obtained using our algorithm were smoother than the FDK results, and had no geometric distortion that was shown in the FDK counterparts 共Figs. 7–11兲. In Figs. 10 and 11, the image noise associated with the proposed algorithm appeared to be larger than that with the FDK algorithm. The reason was that the density drop inherent in the FDK reconstruction, made many pixel values go outside the fixed display window. As a result, a large portion of noise became invisible in the displayed image. The aforementioned image artifacts in our simulation were around the high contrast structures. The profiles across their boundaries can be modeled as piece-wise constant func- TABLE V. Dimensions and linear attenuation coefficients of calcifications. Calcification No. 1 2 3 4 5 6 Diameter 共mm兲 AC 共cm−1兲 0.7 3.282 0.54 3.282 0.4 3.282 0.32 3.282 0.24 3.282 0.16 3.282 Medical Physics, Vol. 33, No. 10, October 2006 3627 Zeng et al.: Cone-beam mammo-CT TABLE VI. Simulated imaging parameters for tilting arc based mammo-CT scan. Number of projections Detector array size Number of detector elements Detector element size Number of photons per view 共for noisy projections兲 Arc radius Arc tilting angle Arc angular range Reconstruction grid Voxel size 500⫻ 2 132⫻ 155 mm2 256⫻ 300 0.57⫻ 0.57 mm2 1e7 ⫻ 256⫻ 300 330 mm 20° 226° 256⫻ 256⫻ 320 0.43⫻ 0.43⫻ 0.43 mm3 3627 TABLE VII. Simulated system setup parameters of a circular scan for FDK algorithm. Number of projections Two-dimensional detector array size Detector elements Detector element size Radius of circle Total photon number per view 共for noisy projections兲 Reconstruction grid Voxel size 1000 132⫻ 133 mm2 256⫻ 256 0.57⫻ 0.57 mm2 330 mm 1e7 ⫻ 256⫻ 300 256⫻ 256⫻ 320 0.43 mm3 flawless, precisely as we expected, but the results using the FDK algorithm contained significant density dropping artifacts 共Figs. 12兲. tions, which do not satisfy the assumption f 苸 C⬁0 共R3兲 that was made for the formulation of Katsevich’s general exact cone-beam reconstruction scheme.33 Therefore, it should not be surprising that such high-contrast nondifferentiable features had caused intensity fluctuations in the reconstructed images, as we previously discussed.41 To confirm the correctness of the implementation of our proposed reconstruction formula, we also applied it to reconstruct the traditional 3D Shepp-Logan phantom. The reconstruction results looked IV. DISCUSSIONS AND CONCLUSION Our reconstruction algorithm was derived from Katsevich’s general scheme and Katsevich’s circle-and-arc algorithm, which can survive the longitude data truncation. As a result, no x-rays need to be sent into the chest wall. Also, this algorithm can be efficiently implemented in parallel similar to what was recently reported of the parallel implementation of the Katsevich algorithm.42 Additionally, our scanning or- FIG. 6. Representative phantom slices. 共a兲 The ideal phantom slice at x2 = 0.0 mm, 共b兲 x1 = 0.0 mm, 共c兲 x3 = 35 mm, and 共d兲 x3 = 75.0 mm. The dotted lines shows the exact reconstruction region. The display window is 关0.244, 0.255兴. Medical Physics, Vol. 33, No. 10, October 2006 3628 Zeng et al.: Cone-beam mammo-CT FIG. 7. Reconstructed images at x2 = 0.0 mm using the proposed algorithm and the FDK algorithm. The images in the left column are from the proposed algorithm, and that in the right column are from the FDK algorithm. The images in the top row are from noise free projection data, and that in the bottom row from noisy data with the same number of photons. The display window is 关0.244, 0.255兴. bit is symmetric, and may be numerically more stable than other less symmetric scanning trajectories, such as circleplus-line and circle-plus-arc.27,29 Although our algorithm requires a longer scanning trajectory than a circular scan, it does not mean that our algorithm require more radiation dose. In our imaging protocol, we can distribute the same dose as that with the circular scan to each view and reconstruct 3D images from these projections. As our algorithm can utilize all the information from two arcs, the final image shall not be inherently noisier than the circular scan reconstruction. Our numeric experiments have indicated that our results were actually less noisy than that from the circular scan using the FDK algorithm. To compare the noise levels in the images reconstructed using our algorithm and the FDK algorithm respectively, the standard deviations of pixel values within a homogeneous region-of-interest were computed 共Fig. 11兲. Generally speaking, the variance of a stochastic process cannot be computed from a single realization of that process. To compare the noise properties of the reconstructed images, the ensemble averages should be used, instead of the averages based on a single image. These averages are equivalent only in the case of an ergodic stochastic process. However, it is not common in either the quality-control practice or the reconstruction literature to do the scan and reconstruction many times for evaluation of the noise level.40,43–45 Hence, in this project the noise image has been approximately considered as being ergodic. In conclusion, we have proposed a novel cone-beam mammo-CT scheme based on two tilting arcs and formulated Medical Physics, Vol. 33, No. 10, October 2006 3628 FIG. 8. Reconstructed images at x1 = 0.0 mm using the proposed algorithm and the FDK algorithm. The images in the left column are from the proposed algorithm, and that in the right column are from the FDK algorithm. The images in the top row are from noise free projection data, and that in the bottom row from noisy data with the same number of photons. The display window is 关0.244, 0.255兴. a Katsevich-type reconstruction algorithm. The numerical simulation results are consistent with the ideal 3D image volume within the numerical error. Compared with other existing mammo-CT algorithms, our work promises better diagnostic performance for breast imaging, and may have a FIG. 9. Reconstructed images at x3 = 35.0 mm using the proposed algorithm and the FDK algorithm. The images in the left column are from the proposed algorithm, and that in the right column are from the FDK algorithm. The images in the top row are from noise free projection data, and that in the bottom row from noisy data with the same number of photons. The display window is 关0.244, 0.255兴. 3629 Zeng et al.: Cone-beam mammo-CT 3629 x1 = 共R cos s̄a cos t兲 + 共1 − 兲共R cos s̄b cos t兲, 共A1a兲 x2 = 共R sin s̄a兲 + 共1 − 兲共R sin s̄b兲, 共A1b兲 x3 = 共R cos s̄a sin t兲 + 共1 − 兲共− R cos s̄b sin t兲. 共A1c兲 We can rewrite 共A1a兲–共A1c兲 as FIG. 10. Reconstructed images at x3 = 75.0 mm using the proposed algorithm and the FDK algorithm. The images in the left column are from the proposed algorithm, and that in the right column are from the FDK algorithm. The images in the top row are from noise free projection data, and that in the bottom row from noisy data with the same number of photons. The display window is 关0.244, 0.255兴. commercial potential. Our scheme also can be used to other similar tomographic imaging applications, such as singlephoton emission computed tomography. ACKNOWLEDGMENTS This work is partially supported by the NIH/NIBIB Grant Nos. EB002667 and EB004287. APPENDIX A: DEFINITION OF REGION Ū, Ū1, and Ū2 1. Monotony of x3 with respect to sa and sb First, we prove that x3 is monotonic for the first case R sin共sa + smx兲 ⬎ x2 ⬎ R sin共sb − smx兲 with respect to sb as what Katsevich did.29 Without loss of generality, let us parameterize x by 0 ⬍ ⬍ 1 with s̄b = sb − smx and s̄a = sa + smx. Hence, we have x1 = cos s̄a + 共1 − 兲cos s̄b , R cos t 共A2a兲 x2 = sin s̄a + 共1 − 兲sin s̄b , R 共A2b兲 x3 = cos s̄a − 共1 − 兲cos s̄b . R sin t 共A2c兲 Then, let us fix x1 and x2, and let , s̄a be functions of s̄b, i.e., 共s̄b兲 and s̄a共s̄b兲. Differentiating 共A2a兲 and 共A2b兲, we have 共cos s̄a − cos s̄b兲 共sin s̄a − sin s̄b兲 d ds̄b d ds̄b − sin s̄a + cos s̄a ds̄a ds̄b ds̄a ds̄b = 共1 − 兲sin s̄b , 共A3a兲 = − 共1 − 兲cos s̄b . 共A3b兲 Solving 共A3兲 for d / ds̄b and ds̄a / ds̄b, we obtain d ds̄b ds̄a ds̄b = = 共1 − 兲sin共s̄b − s̄a兲 1 − cos共s̄a − s̄b兲 共A4a兲 , 共1 − 兲 . 共A4b兲 Also, differentiating 共A2c兲 with respect to s̄b, we have 1 dx3 d = 共cos s̄a + cos s̄b兲 R sin t ds̄b ds̄b − sin s̄a ds̄a ds̄b + 共1 − 兲sin s̄b . 共A5兲 Noting FIG. 11. Reconstructed images at x3 = 46.8 mm from the proposed algorithm and the FDK algorithm with the same number of photons. The noise level is measured as the standard deviation from the homogeneous region indicated by the dotted rectangle. 共a兲 The noise from the proposed algorithm is 2.03e-0.3, and 共b兲 the noise from the FDK algorithm is 2.92e-03. The display window is 关0.244, 0.255兴. Medical Physics, Vol. 33, No. 10, October 2006 3630 Zeng et al.: Cone-beam mammo-CT 3630 FIG. 12. Reconstructed Shepp-Logan phantom images at x2 = −0.25 and profiles along the dash lines. 共a兲 is reconstructed from the proposed algorithm and 共b兲 from the FDK algorithm. The images are displayed with the window 关1.0, 1.05兴. 2. Upper and lower limits of x3 given „x1 , x2… to admits a PI line sin共s̄b − s̄a兲共cos s̄a + cos s̄b兲 − 共sin s̄b − sin s̄a兲cos共s̄b − s̄a兲 = 共sin s̄b cos s̄a − cos s̄b sin s̄a兲共cos s̄a + cos s̄b兲 − 共sin s̄b − sin s̄a兲共cos s̄a cos s̄b + sin s̄a sin s̄b兲 = cos2 s̄a sin s̄b − sin s̄a cos2 s̄b − sin s̄a sin2 s̄b 共A6兲 + sin2 s̄a sin s̄b = sin s̄b − sin s̄a and substituting 共A4兲 and 共A5兲, we have dx3 ds̄b = 共1 − 兲R sin t 冋 共cos s̄a + cos s̄b兲sin共s̄b − s̄a兲 册 1 − cos共s̄b − s̄a兲 + 共sin s̄b − sin s̄a兲 = 2共1 − 兲R sin t 共sin s̄b − sin s̄a兲 1 − cos共s̄b − s̄a兲 . 共A7兲 Because 0 ⬍ ⬍ 1 and R sins̄a ⬎ x2 ⬎ R sins̄b, dx3 / ds̄b does not change its sign. Hence, x3 is monotonic decreasing with s̄b 共or sb兲. As s̄a 共or sa兲 rotates together with s̄b 共or sb兲 in the same direction, x3 is also monotonic decreasing with s̄a 共or sa兲. Similarly, for the second case R sin共sa + smx兲 ⬍ x2 ⬍ R sin共sb − smx兲, we can prove that x3 is monotonous increasing with s̄b 共or sb兲. As s̄a 共or sa兲 rotates together with s̄b 共or sb兲 in the same direction, x3 is also monotonic increasing with s̄a 共or sa兲. It is easy to see that if we project all the x that admit at least one PI line onto the plane x3 = 0, the results must be within an ellipse 兵共x1 , x2兲 兩 x21 / cos2 t + x22 = R2其. Moreover, given a 共x1 , x2兲 within this ellipse, we can readily get the limit of x3 with respect to the assumption of cases 1 and 2 based on the aforementioned monotony property. If we do so for every point 共x1 , x2兲 within the ellipse, we can obtain an image showing the appearance of the region. Based on assumption for the first case, we have x2共y共sa兲兲 艌 x2 艌 x2共y共sb兲兲 where x2共y共sb兲兲 represents the second component of y共sb兲. For any fixed 共x1 , x2兲 within the ellipse and all the points admitting a PI line 共Fig. 13兲, we have the monotony of x3 with respect to sa or sb. Therefore, we can determine the upper and lower limits of x3 for any given 共x1 , x2兲 as follows: Let sa−共x2兲 , sa+共x2兲 be two particular sa± on arc A 关Fig. 13共a兲兴, at which x2共yArcA共sa±兲兲 = x2, x1共yArcA共sa+兲兲 ⬎ 0 and x1共yArcA共sa−兲兲 ⬍ 0. Note that for 兩x2兩 ⬍ 兩R cos共smx兲兩 there does not exist sa−. Let sb−共x2兲, sb+共x2兲 be two particular sb± on arc B 关Fig. 13共b兲兴, at which x2共yArcB共sb±兲兲 = x2, x1共yArcB共sb+兲兲 ⬎ 0 and x1共yArcB共sb−兲兲 ⬍ 0. Note that for some 兩x2兩 ⬍ 兩R cos共smx兲兩 there dose not exist sb+. According to the monotony, i.e., x3 is monotonous decreasing with sa, the upper limit of x3 can be reached at the FIG. 13. Projection of the scanning arcs, the Pl line 共case 1兲 and X onto the traverse plane. 共a兲 The projection of ArcA, and 共b兲 the projection of ArcB. The dotted curve shows the range of sa and sb. Medical Physics, Vol. 33, No. 10, October 2006 3631 Zeng et al.: Cone-beam mammo-CT 3631 minimum sa, i.e., sa = sa+. As far as the lower limit of x3 case1共x1 , x2兲 is concerned, it is more complicated to calculate the maximum sa for a particular 共x1 , x2兲. We must check all the possible boundaries from arc A 关the dotted line in Fig. 14共a兲兴, such as sa− and sa max. Besides, as y共sa兲, x, and y共sb兲 are on the same line, the maximum sa is also limited by the boundaries on arc B 关the dotted line in Fig. 14共b兲兴, such as sb+ and sb max. Hence, the lower limit of x3 is the maximum one of those corresponding x3s at those four critical positions. For example, in Fig. 14, the maximum sa is achieved when the PI line intersects arc B on y共sb max兲. Finally, the strategy to determine the limits can be described as follows: Up共x3 case 1共x1,x2兲兲 = x3共sa+兲 = x3共sb−兲 冦 x3共sa−兲, if sa− exists 冉 冊 2 Low共x3 case 1共x1,x2兲兲 = max x3共sb+兲, if sb+ exists 3 x3 sb max = + 2sm 2 x3 sa max = 冉 冊 冧 where Up共x3 case1共x1 , x2兲兲 is the upper limit of x3, Low共x3 case1共x1 , x2兲兲 is the lower limit of x3. And for the second case’s assumption, we can analysis it similarly. 3. Definition of Ū, Ū1, and Ū2 Now, we can define the region Ū as 共A9a兲, which is the disjunction of the regions determined in cases 1 and 2. Based on the earlier results, it is easy to see that for any x 苸 Ū, it admits at least one PI line and at most two PI lines. The upper and lower surfaces of the region Ū are shown in Fig. 14. Furthermore, the two PI-segments region Ū2 is defined as in 共A9b兲, which admits two PI segments and has the same 共A8兲 upper surface as region Ū. Ū2’s lower surface is shown in Fig. 14共c兲. Ū1 is defined as 共A9c兲, which admits only one PI line. Figure 14共d兲 shows the relative position between Ū2 and Ū1, FIG. 14. Geometry of the exact reconstruction region Ū. 共a兲 The upper surface of the region Ū with respect to two tilting arcs, 共b兲 the lower surface of the region Ū, 共c兲 the lower surface of the region that allows 2 Pl lines, and 共d兲 the 2 Pl lines region vs the 1 Pl line region 共the 1 Pl line region in green, and the 2 Pl line region in red兲. Medical Physics, Vol. 33, No. 10, October 2006 3632 Zeng et al.: Cone-beam mammo-CT 3632 冦冨 冦冨 冧 冧 x21 + x22 ⬍ R2, and given a 共x1,x2兲, Ū = x cos2 t , x3 苸 关min共Lowcase 1,Lowcase 2兲,max共Upcase 1,Upcase 2兲兴 共A9a兲 x21 2 2 2 + x2 ⬍ R , and given a 共x1,x2兲, cos t , Ū2 = x x3 苸 关max共Lowcase 1,Lowcase 2兲,min共Upcase 1,Upcase 2兲兴 共A9b兲 Ū1 = Ū − Ū2 . 共A9c兲 APPENDIX B: PROOF OF THE IMPOSSIBLE CASE †FIG. 4„c…‡ Again, let us denote s̄b = sb − smx, s̄ = s − smx, and s̄t = st + smx. If the case exists 关Fig. 4共c兲兴, as pointed out by Katsevich29 for the circle-and-arc trajectory, there must be a plane tangent to Ca at s̄t. Thus, we must have 冨 − sin s̄t cos t cos s̄t − sin s̄t sin t 共cos s̄ − cos s̄t兲cos t sin s̄ − sin st − 共cos s̄t + cos s̄兲sin t 共cos s̄b − cos s̄t兲cos t sin s̄b − sin st − 共cos s̄t + cos s̄b兲sin t where 3 / 2 ⬎ s̄ ⬎ s̄b ⬎ , sin s̄b ⬍ sin s̄t. Simplifying 共B1兲, we have 2 cos t sin t关sin s̄t sin共s̄ − s̄b兲 + 共cos s̄ − cos s̄b兲兴 = 0. 共B2兲 Since 0 ⬍ t ⬍ / 2, cos t sin t ⫽ 0, we have sin s̄t sin共s̄ − s̄b兲 + 共cos s̄ − cos s̄b兲 = 0. 共B3兲 On the other hand, we obtain that cos s̄ − cos s̄b + sin s̄t sin共s̄ − s̄b兲 =cos s̄ − cos s̄b − sin s̄b sin共s̄b − s̄兲 + 共sin s̄b − sin s̄t兲sin共s̄b − s̄兲 =cos s̄ − cos s̄b − sin2 s̄b cos s̄ + sin s̄b cos s̄b sin s̄ + 共sin s̄b − sin s̄t兲sin共s̄b − s̄兲 =− cos s̄b共1 − sin s̄b sin s̄兲 + cos2 s̄b cos s̄ + 共sin s̄b − sin s̄t兲sin共s̄b − s̄兲 =− cos s̄b共1 − cos共s̄b − s̄兲兲 + 共sin s̄b − sin s̄t兲sin共s̄b − s̄兲 ⬎ 0, 共B4兲 where we have used the fact that for 3 / 2 ⬎ s̄ ⬎ s̄b ⬎ , sin s̄b ⬍ sin s̄t. Therefore, we arrive at a contradiction. a兲 Electronic mail: [email protected] Electronic mail: [email protected] c兲 Electronic mail: [email protected] d兲 Corresponding author. Electronic mail: [email protected] 1 L. T. Niklason, “Current and future developments in digital mammography,” Eur. J. Cancer 38, S14–S15 共2002兲. 2 J. M. Lewin, C. J. D’Orsi, and R. E. Hendrick, “Digital mammography,” Radiol. Clin. North Am. 42, 871–84 共2004兲. 3 A. I. Mushlin, R. W. Kouides, and D. E. Shapiro, “Estimating the accub兲 Medical Physics, Vol. 33, No. 10, October 2006 冨 = 0, 共B1兲 racy of screening mammography: A meta-analysis,” Am. J. Prev. Med. 14, 143–153 共1998兲. 4 D. G. Grant, “Tomosynthesis: A three-dimensional radiographic imaging technique,” IEEE Trans. Biomed. Eng. 19, 20–28 共1972兲. 5 J. T. Dobbins III and D. J. Godfrey, “Digital x-ray tomosynthesis: Current state of the art and clinical potential,” Phys. Med. Biol. 48, R65–R106 共2003兲. 6 E. D. Pisano and C. A. Parham, “Digital mammography, sestamibi breast scintigraphy, and positron emission tomography breast imaging,” Radiol. Clin. North Am. 38, 861–986 共2000兲 x. 7 L. T. Niklason et al., “Digital tomosynthesis in breast imaging,” Radiology 205, 399–406 共1997兲. 8 G. M. Stevens et al., “Circular tomosynthesis: Potential in imaging of breast and upper cervical spin—Preliminary phantom and in vitro study,” Radiology 228, 569–575 共2003兲. 9 D. P. Chakraborty et al., “Self-masking subtraction tomosynthesis,” Radiology 150, 225–229 共1984兲. 10 D. N. Ghosh Roy et al., “Selective plane removal in limited angle tomographic imaging,” Med. Phys. 12, 65–70 共1985兲. 11 D. J. Godfrey, A. Rader, and J. T. Dobbins III, “Practical strategies for the clinical implementation of matrix inversion tomosynthesis 共MITS兲,” Medical Imaging 2003: Physics of Medical Imaging, Feb 16–18 共The International Society for Optical Engineering, San Diego, CA, 2003兲. 12 R. Gordon, R. Bender, and G. T. Herman, “Algebraic reconstruction techniques 共ART兲 for three-dimensional electron microscopy and x-ray photography,” J. Theor. Biol. 29, 471–481 共1970兲. 13 P. Bleuet et al., “An adapted fan volume sampling scheme for 3-D algebraic reconstruction in linear tomosynthesis,” IEEE Trans. Nucl. Sci. 49, 2366–2372 共2002兲. 14 T. Wu et al., “A comparison of reconstruction algorithms for breast tomosynthesis,” Med. Phys. 31, 2636–2647 共2004兲. 15 C. H. Chang et al., “Computed tomographic evaluation of the breast,” AJR, Am. J. Roentgenol. 131, 459–464 共1978兲. 16 B. Chen and R. Ning, “Cone-beam volume CT breast imaging: Feasibility study,” Med. Phys. 29, 755–770 共2002兲. 17 J. Zhong, R. Ning, and D. Conover, “Image denoising based on multiscale singularity detection for cone beam CT breast imaging,” IEEE Trans. Med. Imaging 23, 696–703 共2004兲. 18 L. A. Feldkamp, L. C. Davis, and J. W. Kress, “Practical cone-beam algorithm,” J. Opt. Soc. Am. A 1, 612–619 共1984兲. 19 P. Grangeat, “Mathematical frame work of cone beam 3D reconstruction via the first derivative of the radon transform,” Lect. Notes Math. 1497, 66–97 共1991兲. 3633 Zeng et al.: Cone-beam mammo-CT 20 B. D. Smith, “Image reconstruction from cone-beam projections: Necessary and sufficient conditions and reconstruction methods,” IEEE Trans. Med. Imaging MI-4, 14–25 共1985兲. 21 H. K. Tuy, “An inverse formula for cone-beam reconstruction,” SIAM J. Appl. Math. 43, 546–552 共1983兲. 22 S. Zhao, H. Yu, and G. Wang, “A unified framework for exact cone-beam reconstruction formulas,” Med. Phys. 32, 1712–1721 共2005兲. 23 A. Katsevich, “An improved exact filtered backprojection algorithm for spiral computed tomography,” Adv. Appl. Math. 32, 681–697 共2004兲. 24 F. Noo, J. Pack, and D. Heuscher, “Exact helical reconstruction using native cone-beam geometries,” Phys. Med. Biol. 48, 3787–3818 共2003兲. 25 H. Yu and G. Wang, “Studies on implementation of the katsevich algorithm for spiral cone-beam CT,” J. X-Ray Sci. Technol. 12, 97–116 共2004兲. 26 H. Hu, “Exact regional reconstruction of longitudinally unbounded objects using the circle-and-line cone-beam tomographic system,” Medical Imaging 1997: Physics of Medical Imaging 共SPIE, Newport Beach, CA, 1997兲. 27 A. Katsevich, “Image reconstruction for the circle and line trajectory,” Phys. Med. Biol. 49, 5059–5072 共2004兲. 28 X. Tang and R. Ning, “A cone beam filtered backprojection 共CB-FBP兲 reconstruction algorithm for a circle-plus-two-arc orbit,” Med. Phys. 28, 1042–1055 共2001兲. 29 A. Katsevich, “Image reconstruction for the circle-and-arc trajectory,” Phys. Med. Biol. 50, 2249–2265 共2005兲. 30 J. D. Pack, F. Noo, and H. Kudo, “Investigation of saddle trajectories for cardiac CT imaging in cone-beam geometry,” Phys. Med. Biol. 49, 2317–2336 共2004兲. 31 H. Yu et al., “Exact BPF and FBP algorithms for nonstandard saddle curves,” Med. Phys. 32, 3305 共2005兲. 32 Y. Yangbo et al., “A general exact reconstruction for cone-beam CT via backprojection-filtration,” IEEE Trans. Med. Imaging 24, 1190 共2005兲. 33 A. Katsevich, “A general scheme for constructing inversion algorithms Medical Physics, Vol. 33, No. 10, October 2006 3633 for cone beam CT,” Int. J. Math. Math. Sci. 21, 1305–1321 共2003兲. J. D. Pack, F. Noo, and R. Clackdoyle, “Cone-beam reconstruction using the backprojection of locally filtered projections,” IEEE Trans. Med. Imaging 24, 70–85 共2005兲. 35 Y. Ye and G. Wang, “Filtered backprojection formula for exact image reconstruction from cone-beam data along a general scanning curve,” Med. Phys. 32, 42–48 共2005兲. 36 Y. Zou and X. Pan, “An extended data function and its generalized backprojection for image reconstruction in helical cone-beam CT,” Phys. Med. Biol. 49, N383–N387 共2004兲. 37 J. Law, “A new phantom for mammography,” Br. J. Radiol. 64, 116–120 共1991兲. 38 ICRU, “Phantoms and computational models in therapy, diagnosis and protection,” ICRU Report No. 48, 1992. 39 ICRU, “Tissue substitutes in radiation dosimetry and measurement,” ICRU Report No. 44, 1989. 40 G. Wang and M. W. Vannier, “Low-contrast resolution in volumetric x-ray CT—analytical comparison between conventional and spiral CT,” Med. Phys. 24, 373–376 共1997兲. 41 H. Yu, S. Zhao, and G. Wang, “A differentiable Shepp-Logan phantom and its applications in exact cone-beam CT,” Phys. Med. Biol. 50, 5583– 5595 共2005兲. 42 J. Deng et al., “A parallel implementation of the Katsevich algorithm for 3-D CT image reconstruction,” J. Supercomput. 共in press兲. 43 L. Yu and X. Pan, “Half-scan fan-beam computed tomography with improved noise and resolution properties,” Med. Phys. 30, 2629–2637 共2003兲. 44 T. Fuchs et al., “Spiral interpolation algorithms for multislice spiral CT. II. Measurement and evaluation of slice sensitivity profiles and noise at a clinical multislice system,” IEEE Trans. Med. Imaging 19, 835 共2000兲. 45 A. Polacin, W. A. Kalender, and G. Marchal, “Evaluation of section sensitivity profiles and image noise in spiral CT,” Radiology 185, 29–35 共1992兲. 34