WEIGHTED ANALYTIC FILTERED BACK PROJECTION RECONSTRUCTION METHOD AND SYSTEM FOR ASYMMETRIC CONE ANGLE ARTIFACTS

Abstract

Disclosed in the present invention are a weighted analytic filtered back projection reconstruction method and system for asymmetric cone angle artifacts. The method comprises the following steps: dividing a reconstruction area into a plurality of weight regions on the basis of relative positions of a ray source ring and a detector ring; acquiring the projection data volume of voxel points in each weight area irradiated by X-rays; according to the projection data volume of the voxel points in each weight area irradiated by the X-rays, assigning a different initial weight to each weight area; performing smooth transition on the initial weight of each weight area by means of a transition weight to form a final weight assigned to each weight area; and according to different final weights of the weight regions, performing final weighted analytic reconstruction on projection data p (, , ) to acquire a back projection image.

Claims

1. A weighted analytic filtered back projection reconstruction method for asymmetric cone angle artifacts, comprising the following steps: dividing a reconstruction region into a plurality of weight regions on the basis of relative positions of a ray source ring and a detector ring, wherein the ray source ring and the detector ring are mutually staggered to form asymmetric cone angle artifacts; acquiring the projection data volume of voxel points in each weight region irradiated by X-rays; according to the projection data volume of the voxel points in each weight region irradiated by the X-rays, assigning a different initial weight to each weight region; performing smooth transition on the initial weight of each weight region by means of a transition weight, to form a final weight assigned to each weight region; and according to the different final weights of the weight regions, performing final weighted analytic reconstruction on projection data p(, , ) collected under large cone beam opening angle geometry, to acquire a back projection image.

2. The weighted analytic filtered back projection reconstruction method according to claim 1, wherein the dividing a reconstruction region into a plurality of weight regions on the basis of relative positions of a ray source ring and a detector ring specifically comprises: constructing a space coordinate system by using a center of the detector ring as a circle center O, using an axis of the detector ring as a Z axis, using a horizontal radial direction of the detector ring as an X axis, and using a vertical radial direction of the detector ring as a Y axis, wherein in the space coordinate system, cone beam opening angles of the X-rays received by the detector ring fall within a range [.sub.2, .sub.1], and a rotation angle corresponding to coordinates v=(x, y, z) of a voxel point satisfies x=r sin and y=r cos ; and dividing the reconstruction region into a total of eight weight regions A to H by using the cone beam opening angles of the X-rays on two sides of the Z axis, wherein a range (v, ) of angles at which the eight weight regions are irradiated by the X-rays is as follows: $\left(v,\right)=\{\begin{array}{cc}\left(\mathrm{empty}\mathrm{set}\right)& \mathrm{if}r<c_{1}z-R_S\left(\mathrm{Region}A\right)\\ \left[-\frac{}{2}-{}_1,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}r{c}_{1}z-R_{S}0\left(\mathrm{Region}B\right)\\ \left[-\frac{}{2}-{}_1,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}{c}_{2}z-R_{S}r{R}_S-c_{1}z>0\left(\mathrm{Region}C\right)\\ \left[-,+\right]& \mathrm{else}\mathrm{if}r<R_S-c_{1}z\mathrm{and}r{c}_{2}z-R_S\left(\mathrm{Region}D\right)\\ \left[-\frac{}{2}-{}_1,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}r{R}_S-c_{1}z>0\mathrm{and}r>c_{2}z-R_S>0\left(\mathrm{Region}E\right)\\ \left[-,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\right]& \mathrm{else}\mathrm{if}{R}_S-c_{1}z>r>c_{2}z-R_S>0\left(\mathrm{Region}F\right)\\ \left[-,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\right]& \mathrm{else}\mathrm{if}r{c}_{2}z-R_{S}0\left(\mathrm{Region}G\right)\\ \left(\mathrm{empty}\mathrm{set}\right)& \mathrm{else}\left(\mathrm{Region}H\right)\end{array}$ where c.sub.1=cot .sub.1, c.sub.2=cot .sub.2, and the definitions of .sub.1 and .sub.2 are as follows: ${}_1=\{\begin{array}{cc}\arcsin \frac{R_S-c_{1}z}{r}-\arctan \frac{\sqrt{r^2-{\left(R_S-c_{1}z\right)}^2}}{c_{1}z}& \mathrm{if}r>.Math.R_S-c_{1}z.Math.\\ \mathrm{sgn}\left(R_S-c_{1}z\right)\frac{}{2}& \mathrm{else}\end{array}$ ${}_2=\{\begin{array}{cc}\arcsin \frac{R_S-c_{2}z}{r}-\arctan \frac{\sqrt{r^2-{\left(R_S-c_{2}z\right)}^2}}{c_{2}z}& \mathrm{if}r>.Math.R_S-c_{2}z.Math.\\ \mathrm{sgn}\left(R_S-c_{2}z\right)\frac{}{2}& \mathrm{else}\end{array}$ where R.sub.S represents a vertical distance from a ray source to a center of an FOV; z represents a Z-axis coordinate of a voxel point v; r represents the length of a vector (x, y) in an XY plane; and represents an angle formed between the vector (x, y) in the XY plane and the Y axis.

3. The weighted analytic filtered back projection reconstruction method according to claim 2, wherein the projection data volume (v, ) of the voxel points in each weight region irradiated by the X-rays is calculated by the following formula: $\left(v,\right)=\{\begin{array}{cc}0& \mathrm{if}r<c_{1}z-R_S\left(\mathrm{Region}A\right)\\ +2{}_1& \mathrm{else}\mathrm{if}r{c}_{1}z-R_{S}0\left(\mathrm{Region}B\right)\\ +2{}_1& \mathrm{else}\mathrm{if}{c}_{1}z-R_{S}r{R}_S-c_{1}z>0\left(\mathrm{Region}C\right)\\ 2& \mathrm{else}\mathrm{if}r<R_S-c_{1}z\mathrm{and}r{c}_{2}z-R_S\left(\mathrm{Region}D\right)\\ 2{}_1-2{}_2& \mathrm{else}\mathrm{if}r{R}_S-c_{1}z>0\mathrm{and}r>c_{2}z-R_S>0\left(\mathrm{Region}E\right)\\ -2{}_2& \mathrm{else}\mathrm{if}{R}_S-c_{1}z>r>c_{2}z-R_S>0\left(\mathrm{Region}F\right)\\ -2{}_2& \mathrm{else}\mathrm{if}r{R}_S-c_{2}z0\left(\mathrm{Region}G\right)\\ 0& \mathrm{else}\left(\mathrm{Region}H\right)\end{array}$

4. The weighted analytic filtered back projection reconstruction method according to claim 3, wherein the assigning a different initial weight to each weight region specifically comprises: according to the projection data volume (v, ) of the voxel points in each weight region irradiated by the X-rays, dividing the eight weight regions A to H into: a non-irradiated region, a partially irradiated region, a fully irradiated region, and a non-fixed irradiated region; assigning an initial weight 0 to the non-irradiated region; assigning an initial weight 0<W.sub.PS(, )<1 to the partially irradiated region and the non-fixed irradiated region; and assigning an initial weight W.sub.FS(, )=1 to the fully irradiated region.

5. The weighted analytic filtered back projection reconstruction method according to claim 4, wherein the weight region A and the weight region H are divided as non-irradiated regions; the weight region B and the weight region G are divided as regions irradiated at an angle less than 180; the weight region C and the weight region F are divided as regions irradiated at an angle equal to or greater than 180 and less than 360; the weight region D is divided as a 360 fully irradiated region; and the weight region E is divided as a non-fixed irradiated region; if the region irradiated at an angle equal to or greater than 180 and less than 360 is the weight region C, i.e., c.sub.2zR.sub.S>R.sub.Sc.sub.1z is satisfied, then $W_{PS}\left(v,\right)=\{\begin{array}{cc}0& \mathrm{if}<-/2-{}_1\\ 1+s\left(\frac{-+/2}{{}_1}\right)& \mathrm{else}\mathrm{if}<-/2+{}_1\\ 2& \mathrm{else}\mathrm{if}<+/2-{}_1\\ 1-s\left(\frac{-+/2}{{}_1}\right)& \mathrm{else}\mathrm{if}<+/2+{}_1\\ 0& \mathrm{else}\end{array}$ if the region irradiated at an angle equal to or greater than 180 and less than 360 is the weight region F, i.e., R.sub.Sc.sub.1zc.sub.2zR.sub.S is satisfied, then $W_{PS}\left(v,\right)=\{\begin{array}{cc}2& \mathrm{if}<-/2-{}_2\\ 1+s\left(\frac{-+/2}{{}_2}\right)& \mathrm{else}\mathrm{if}<-/2-{}_2\\ 0& \mathrm{else}\mathrm{if}<+/2+{}_2\\ 1-s\left(\frac{--/2}{{}_2}\right)& \mathrm{else}\mathrm{if}<+/2-{}_2\\ 2& \mathrm{else}\end{array}$ where a function s(t)=sin(t/2); for the region irradiated at an angle less than 180 and/or the non-fixed irradiated region, $W_{PS}\left(v,\right)=\{\begin{array}{cc}\min \left(2,\frac{2}{\left(v,\right)}\right)& \mathrm{if}\left(v,\right)\\ 0& \mathrm{else}\end{array}$ for the 360 fully irradiated region, the same weight W.sub.FS(, )=1 is assigned to all data.

6. The weighted analytic filtered back projection reconstruction method according to claim 5, wherein the performing smooth transition on the initial weight of each weight region by means of a transition weight, to form a final weight assigned to each weight region specifically comprises: introducing the transition weight W.sub.T(, ) to perform smooth transition on the initial weight of each weight region, to form the final weight W.sub.C(, )=(1W.sub.T(, ))W.sub.FS(, )=W.sub.T(, )W.sub.PS(, ), wherein for a non-360 fully irradiated region, W.sub.T(, )=1; and for the 360 fully irradiated region, $W_T\left(v,\right)=\frac{1}{2}\{\begin{array}{cc}0& \mathrm{if}r<r_0-r\\ 1+s\left(\frac{r-r_0}{r}\right)& \mathrm{else}\mathrm{if}r{r}_0+r\\ 2& \mathrm{else}\end{array}$ where r=R.sub.M.sup.2/2R.sub.S, r.sub.0=min(R.sub.Sc.sub.1z, c.sub.2zR.sub.S)r, and R.sub.M is the radius of the FOV.

7. The weighted analytic filtered back projection reconstruction method according to claim 6, wherein final weighted analytic reconstruction is performed by the following formula: $f\left(v\right)=\frac{1}{2}{W}_C\left(v,\right)h\left(\right)*p\left(,,\right)d$ where p(, , ) is the parallel beam geometric projection data obtained after weighted rearrangement of the p(, , ), is a projection image horizontal index variable after the rearrangement, and h() is a filter kernel.

8. A back projection reconstruction system, comprising a processor and a memory, the processor reading a computer program in the memory, and performing the following operations: dividing a reconstruction region into a plurality of weight regions on the basis of relative positions of a ray source ring and a detector ring, wherein the ray source ring and the detector ring are mutually staggered to form asymmetric cone angle artifacts; according to the projection data volume of voxel points in each weight region irradiated by the X-rays, assigning a different initial weight to each weight region; performing smooth transition on the initial weight of each weight region by means of a transition weight, to form a final weight assigned to each weight region; and according to the different final weights of the weight regions, performing final weighted analytic reconstruction on projection data p(, , ) collected under large cone beam opening angle geometry, to acquire a back projection image.

9. The back projection reconstruction system according to claim 8, wherein the ray source ring is formed by a plurality of X-ray sources surrounding in a circle, the detector ring is formed by a plurality of X-ray detectors surrounding in a circle, and the detector ring is arranged on an inner side of the ray source ring.

10. A computer-readable storage medium comprising a program instruction, the program instruction, when executed by a processor, implementing the weighted analytic filtered back projection reconstruction method according to claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0047] FIG. 1 is a schematic diagram of a geometric configuration to which an xFDK weighting method is applicable in the related art;

[0048] FIG. 2 is a schematic diagram of a geometric configuration to which a bixFDK weighting method is applicable according to an embodiment of the present disclosure;

[0049] FIG. 3 is a flowchart of a weighted analytic filtered back projection reconstruction method for asymmetric cone angle artifacts provided in an embodiment of the present disclosure;

[0050] FIG. 4 is a slice view of a ray source ring and a detector ring in a space coordinate system according to an embodiment of the present disclosure;

[0051] FIG. 5 is a sagittal view of a ray source ring and a detector ring in a space coordinate system according to an embodiment of the present disclosure;

[0052] FIG. 6 is a schematic diagram of regions obtained by dividing a reconstruction region according to an embodiment of the present disclosure;

[0053] FIG. 7 is a schematic diagram of a projection image reconstructed by a generic FDK algorithm;

[0054] FIG. 8 is a schematic diagram of a projection image reconstructed by a generic xFDK weighting algorithm;

[0055] FIG. 9 is a schematic diagram of a projection image reconstructed by a bixFDK weighting algorithm according to an embodiment of the present disclosure; and

[0056] FIG. 10 is schematic structural diagram of a back projection reconstruction system for asymmetric cone angle artifacts provided in an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0057] The technical contents of the present disclosure will be specifically described in detail below with reference to the accompanying drawings and specific embodiments.

[0058] As shown in FIG. 2, an embodiment of the present disclosure first provides a weighted analytic filtered back projection reconstruction method for asymmetric cone angle artifacts, which compensate for a cone angle shadow caused by incomplete data based on an analytic derivation back projection weighting factor calculation method, to separately process a near-end cone angle and a far-end cone angle of the asymmetric cone angle artifacts, so that both ends can obtain correct weights, and thereby the cone angle is correctly compensated for.

[0059] Compared with an existing cone angle weighted reconstruction algorithm, the method has advantages of accurate compensation, higher calculation speed, and introduction of no other artifacts.

[0060] As shown in FIG. 3, the weighted analytic filtered back projection reconstruction method for asymmetric cone angle artifacts according to the embodiment of the present disclosure specifically includes steps S1 to S4:

[0061] S1: A reconstruction region is divided into a plurality of weight regions on the basis of relative positions of a ray source ring and a detector ring.

[0062] Specifically, in this embodiment, the ray source ring and the detector ring are mutually staggered to form asymmetric cone angle artifacts. Moreover, as shown in FIG. 4 and FIG. 5, a space coordinate system is constructed by using a center of the detector ring as a circle center O, using an axis of the detector ring as a Z axis, using a horizontal radial direction of the detector ring as an X axis, and using a vertical radial direction of the detector ring as a Y axis.

[0063] In the space coordinate system, cone beam opening angles of the X-rays received by the detector ring fall within a range [.sub.2, .sub.1], and a rotation angle corresponding to coordinates v=(x, y, z) of a voxel point satisfies x=r sin and y=r cos .

[0064] As shown in FIG. 6, in this embodiment, the reconstruction region is divided into a total of eight weight regions A to H by using the cone beam opening angles of the X-rays on two sides of the Z axis.

[0065] A range (v, ) of angles at which the eight weight regions are irradiated by the X-rays is as follows:

[00009]$\left(v,\right)=\{\begin{array}{cc}\left(\mathrm{empty}\mathrm{set}\right)& \mathrm{if}r<c_{1}z-R_s\left(\mathrm{Region}A\right)\\ \left[-\frac{}{2}-{}_1,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}r<c_{1}z-R_{s}0\left(\mathrm{Region}B\right)\\ \left[-\frac{}{2}-{}_1,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}{c}_{2}z-R_{s}r{R}_s-c_{1}z>0\left(\mathrm{Region}C\right)\\ \left[-,+\right]& \mathrm{else}\mathrm{if}r<R_s-c_{1}z\mathrm{and}r{c}_{2}z-R_s\left(\mathrm{Region}D\right)\\ \left[-\frac{}{2}-{}_1,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\frac{}{2}+{}_1\right]& \mathrm{else}\mathrm{if}r{R}_s-c_{1}z>0\mathrm{and}r>c_{2}z-R_s>0\left(\mathrm{Region}E\right)\\ \left[-,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\right]& \mathrm{else}\mathrm{if}{R}_s-c_{1}z>r>c_{2}z-R_s>0\left(\mathrm{Region}F\right)\\ \left[-,-\frac{}{2}-{}_2\right].Math.\left[+\frac{}{2}+{}_2,+\right]& \mathrm{else}\mathrm{if}r{c}_{2}z-R_{s}0\left(\mathrm{Region}G\right)\\ \left(\mathrm{empty}\mathrm{set}\right)& \mathrm{else}\left(\mathrm{Region}H\right)\end{array}$ [0066] where c.sub.1=cot .sub.1, c.sub.2=cot .sub.2, and the definitions of .sub.1 and .sub.2 are as follows:

[00010]${}_1=\{\begin{array}{cc}\arcsin \frac{R_S-c_{1}z}{r}-\arctan \frac{\sqrt{r^2-{\left(R_S-c_{1}z\right)}^2}}{c_{1}z}& \mathrm{if}r>.Math.R_S-c_{1}z.Math.\\ \mathrm{sgn}\left(R_S-c_{1}z\right)\frac{}{2}& \mathrm{else}\end{array}$${}_2=\{\begin{array}{cc}\arcsin \frac{R_S-c_{2}z}{r}-\arctan \frac{\sqrt{r^2-{\left(R_S-c_{2}z\right)}^2}}{c_{2}z}& \mathrm{if}r>.Math.R_S-c_{2}z.Math.\\ \mathrm{sgn}\left(R_S-c_{2}z\right)\frac{}{2}& \mathrm{else}\end{array}$ [0067] where R.sub.S represents a vertical distance from a ray source to a center of an FOV; z represents a Z-axis coordinate of a voxel point v; r represents the length of a vector (x, y) in an XY plane; and represents an angle formed between the vector (x, y) in the XY plane and the Y axis.

[0068] S2: The projection data volume of voxel points in each weight region irradiated by the X-rays is acquired.

[0069] Specifically, in this embodiment, the projection data volume (v, ) of the voxel points in each weight region irradiated by the X-rays is calculated by the following formula:

[00011]$\left(v,\right)=\{\begin{array}{cc}0& \mathrm{if}r<c_{1}z-R_S\left(\mathrm{Region}A\right)\\ +2{}_1& \mathrm{else}\mathrm{if}r{c}_{1}z-R_{S}0\left(\mathrm{Region}B\right)\\ +2{}_1& \mathrm{else}\mathrm{if}{c}_{1}z-R_{S}r{R}_S-c_{1}z>0\left(\mathrm{Region}C\right)\\ 2& \mathrm{else}\mathrm{if}r<R_S-c_{1}z\mathrm{and}r{c}_{2}z-R_S\left(\mathrm{Region}D\right)\\ 2{}_1-2{}_2& \mathrm{else}\mathrm{if}r{R}_S-c_{1}z>0\mathrm{and}r>c_{2}z-R_S>0\left(\mathrm{Region}E\right)\\ -2{}_2& \mathrm{else}\mathrm{if}{R}_S-c_{1}z>r>c_{2}z-R_S>0\left(\mathrm{Region}F\right)\\ -2{}_2& \mathrm{else}\mathrm{if}r{R}_S-c_{2}z0\left(\mathrm{Region}G\right)\\ 0& \mathrm{else}\left(\mathrm{Region}H\right)\end{array}$

[0070] S3: According to the projection data volume of the voxel points in each weight region irradiated by the X-rays, a different initial weight is assigned to each weight region.

[0071] In this embodiment, according to the projection data volume of the voxel points in each weight region irradiated by the X-rays, the eight weight regions A to H are divided into four types of regions: a non-irradiated region, a partially irradiated region, a fully irradiated region, and a non-fixed irradiated region. An initial weight 0 is assigned to the non-irradiated region. An initial weight 0<W.sub.PS(, )<1 is assigned to the partially irradiated region and the non-fixed irradiated region. An initial weight W.sub.FS(, )=1 is assigned to the fully irradiated region.

[0072] Specifically, the weight region A and the weight region H are divided as non-irradiated regions; the weight region B and the weight region G are divided as regions irradiated at an angle less than 180; the weight region C and the weight region F are divided as regions irradiated at an angle equal to or greater than 180 and less than 360; the weight region D is divided as a 360 fully irradiated region; and the weight region E is divided as a non-fixed irradiated region. The regions irradiated at an angle less than 180 and the regions irradiated at an angle equal to or greater than 180 and less than 360 are all partially irradiated regions, but the specific initial weight W.sub.PS(, ) assigned thereto has certain differences as follows:

[0073] If the region irradiated at an angle equal to or greater than 180 and less than 360 is the weight region C, i.e., c.sub.2zR.sub.S>R.sub.Sc.sub.1z is satisfied, then

[00012]$W_{PS}\left(v,\right)=\{\begin{array}{cc}0& \mathrm{if}<-/2-{}_1\\ 1+s\left(\frac{-+/2}{{}_1}\right)& \mathrm{else}\mathrm{if}<-/2+{}_1\\ 2& \mathrm{else}\mathrm{if}<+/2-{}_1\\ 1-s\left(\frac{-+/2}{{}_1}\right)& \mathrm{else}\mathrm{if}<+/2+{}_1\\ 0& \mathrm{else}\end{array}$

[0074] If the region irradiated at an angle equal to or greater than 180 and less than 3600 is the weight region F, i.e., R.sub.Sc.sub.1zc.sub.2zR.sub.S is satisfied, then

[00013]$W_{PS}\left(v,\right)=\{\begin{array}{cc}2& \mathrm{if}<-/2-{}_2\\ 1+s\left(\frac{-+/2}{{}_2}\right)& \mathrm{else}\mathrm{if}<-/2-{}_2\\ 0& \mathrm{else}\mathrm{if}<+/2+{}_2\\ 1-s\left(\frac{--/2}{{}_2}\right)& \mathrm{else}\mathrm{if}<+/2-{}_2\\ 2& \mathrm{else}\end{array}$ [0075] where a function s(t)=sin(t/2).

[0076] For the region irradiated at an angle less than 1800 and/or the non-fixed irradiated region, i.e., the weight region B, the weight region G, and/or the weight region E,

[00014]$W_{PS}\left(v,\right)=\{\begin{array}{cc}\min \left(2,\frac{2}{\left(v,\right)}\right)& \mathrm{if}\left(v,\right)\\ 0& \mathrm{else}\end{array}$

[0077] For the 360 fully irradiated region, i.e., the weight region D, the same weight W.sub.FS(, )=1 is assigned to all data.

[0078] S4: Smooth transition is performed on the initial weight of each weight region by means of a transition weight, to form a final weight assigned to each weight region.

[0079] Specifically, in this embodiment, the transition weight W.sub.T(, ) is introduced to perform smooth transition on the initial weight of each weight region, to form the final weight W.sub.C(, )=(1W.sub.T(, ))W.sub.FS(, )=W.sub.T(, )W.sub.PS(, ) assigned to each weight region.

[0080] For a non-360 fully irradiated region (i.e., the other seven weight regions than the weight region D), W.sub.T(, )=1.

[0081] For the 360 fully irradiated region (i.e., the weight region D),

[00015]$W_T\left(v,\right)=\frac{1}{2}\{\begin{array}{cc}0& \mathrm{if}r<r_0-r\\ 1+s\left(\frac{r-r_0}{r}\right)& \mathrm{else}\mathrm{if}r{r}_0+r\\ 2& \mathrm{else}\end{array}$ [0082] where r=R.sub.M.sup.2/2R.sub.S, r.sub.0=min(R.sub.Sc.sub.1z, c.sub.2zR.sub.S)r, and R.sub.M is the radius of the FOV (as shown in FIG. 4).

[0083] S5: According to the different final weights of the weight regions, final weighted analytic reconstruction is performed on the projection data p(, , ) collected under large cone beam opening angle geometry, to acquire a back projection image.

[0084] Specifically, in this embodiment, final weighted analytic reconstruction is performed by the following formula:

[00016]$f\left(v\right)=\frac{1}{2}{W}_C\left(v,\right)h\left(\right)*p\left(,,\right)d$ [0085] where p(, , ) is the parallel beam geometric projection data obtained after weighted rearrangement of the p(, , ), is a projection image horizontal index variable after the rearrangement, and h() is a filter kernel.

[0086] As shown in FIG. 7 to FIG. 9, FIG. 7 is a projection image reconstructed by a generic FDK algorithm, FIG. 8 is a projection image reconstructed by a generic xFDK weighting algorithm, and FIG. 9 is a projection image reconstructed by a bixFDK weighting algorithm in this embodiment. From the figures, there are obvious cone angle shadows at the upper and lower edges of the reconstruction result obtained by the generic FDK algorithm. By the xFDK algorithm, only a shadow at the upper edge can be reduced, and there is still a shadow in a region below a dashed line at a lower edge. By the bixFDK algorithm in this embodiment, shadows at the upper and lower edges are significantly reduced.

[0087] Therefore, by the steps S1 to S5, an asymmetric cone angle artifact removal method for analytic reconstruction in static CT is designed and implemented. The method can effectively estimate and compensate for asymmetric cone angle artifacts caused by staggering of ray sources and detectors. After analytic reconstruction is weighted according to the method, a reconstructable range of a single axial scan in static CT can be greatly improved, so that the uniformity and accuracy of a CT value at an axial scan analytic FBP reconstruction edge layer can be significantly improved. The utilization rate of a ray dose in static CT is indirectly improved by expanding the reconstructable range, and the ray dose required for static CT imaging is effectively reduced. In addition, it can be understood that in this embodiment, a weighting formula only applicable to a symmetric cone angle configuration system is extended to be also applicable to an asymmetric cone angle configuration system, so that the method can be applied to geometric configuration in static CT.

[0088] Based on the weighted analytic filtered back projection reconstruction method, the present disclosure further provides a back projection reconstruction system. As shown in FIG. 10, the back projection reconstruction system includes one or more processors 21 and a memory 22. The memory 22 is coupled to the processors 21 and is configured to store one or more programs. When the one or more programs are executed by the one or more processors 21, the one or more processors 21 implement the weighted analytic filtered back projection reconstruction method as described in the above embodiment.

[0089] The processors 21 are configured to control overall operations of the back projection reconstruction system, to perform all or some steps of the weighted analytic filtered back projection reconstruction method. The processors 21 may be central processing units (CPU), graphic processing units (GPU), field programmable gate arrays (FPGA), application-specific integrated circuits (ASIC), digital signal processing (DSP) chips, or the like. The memory 22 is configured to store various types of data to support operations on the back projection reconstruction system. The data may include, for example, any application program or an instruction of method for operations on the back projection reconstruction system, and data related to the application program. The memory 22 may be implemented by any type of volatile or non-volatile storage device or a combination thereof, for example, a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM), an erasable programmable read-only memory (EPROM), a programmable read-only memory (PROM), a read-only memory (ROM), a magnetic memory, a flash memory, or the like.

[0090] In an exemplary embodiment, the back projection reconstruction system may be specifically implemented by a computer chip or an entity, or implemented by a product having a certain function, to perform the weighted analytic filtered back projection reconstruction method, and achieve a technical effect consistent with the method. A typical embodiment is a computer. Specifically, the computer may be, for example, a personal computer, a laptop, an in-car human-machine interaction device, a cellular phone, a camera phone, a smartphone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or any combination of the devices.

[0091] In another exemplary embodiment, the present disclosure further provides a computer-readable storage medium including a program instruction. The program instruction, when executed by a processor, performs the steps of the weighted analytic filtered back projection reconstruction method in any of the above embodiments. For example, the computer-readable storage medium may be the aforementioned memory including the program instruction, and the program instruction may be executed by a processor of a back projection reconstruction system to perform the weighted analytic filtered back projection reconstruction method, and achieve a technical effect consistent with the method.

[0092] Based on the above, the weighted analytic filtered back projection reconstruction method and system for asymmetric cone angle artifacts provided in the embodiments of the present disclosure have the following beneficial effects:

[0093] 1. The method and system can effectively estimate and compensate for asymmetric cone angle artifacts caused by staggering of ray sources and detectors. After analytic reconstruction is weighted according to the method, a reconstructable range of a single axial scan in static CT can be greatly improved, and the ray dose required for static CT imaging is effectively reduced.

[0094] 2. A weighting formula only applicable to a symmetric cone angle configuration system is extended to be also applicable to an asymmetric cone angle configuration system, so that the method can be applied to geometric configuration in static CT.

[0095] The weighted analytic filtered back projection reconstruction method and system for asymmetric cone angle artifacts provided by the present disclosure are described in detail above. Any obvious change made by a person of ordinary skill in the art to the present disclosure without departing from the essence of the present disclosure shall constitute an infringement of the patent right of the present disclosure, and the person shall bear corresponding legal responsibility.