PERFORMANCE DIFFERENCE

FIELD

[0001] The present invention generally relates to a computer image processing field, particularly to a method for removing streak from detector cell with performance difference.

BACKGROUND

[0002] A computed tomography (CT) is a common endoscopy detection technology in the modern medicine. The CT technology obtains an image of a detected region of a human body by mainly rotating X-ray beams and a detector around the human body, and continuously carrying out section scanning, receiving attenuated X-ray information penetrating the human body by the detector and inputting said information into a computer, reconstructing an image through the computer based on the received attenuated X-ray information during each scanning.

[0003] The image reconstruction in the CT technology is an important factor that would affect the detection result. How to clearly and accurately reproduce an image of the detected region is a problem that has always been studied by the person skilled in the art.

[0004] From an imaging principle of a CT system, there are multiple factors that would affect a quality of a CT image.For example, in the existing CT system, due to reasons on the structure and manufacturing process of the detector itself, there would inevitably be some performance difference between detector cells on the detector , i.e., there would be some detector cells with performance difference. Projections sensed by those detector cells with performance difference cannot accurately reflect the image of the detected region, thus cannot be applied to the CT image reconstruction directly.

[0005] An existing processing method is: in the case that the detector cells with performance difference are determined to be present, usually by an interpolation method, a projection on the detector cell with performance difference is estimated from projections sensed on two detector cells adjacent to the detector cell with performance difference.

However, there would be some estimation error between a correct projection of the detector cell with performance difference and the estimated projection. Therefore, during the CT image reconstruction, when an image is reconstructed with projections using a filter backprojection method, the estimation error will be back-projected onto the reconstructed image, thereby forming a streak on the reconstructed image, thus affecting the quality of the CT reconstructed image.

[0006] Accordingly, it is necessary to provide an improved method for solving at least one of the technical problems as mentioned above.

SUMMARY

[0007] One aspect of the present invention is to provide a method for removing streaks from detector cells with performance difference, the method comprising the steps of:

al) estimating a projection on a detector cell with performance difference for each of one or more views;

a2) reconstructing one or more initial images with the estimated projection on the detector cell with performance difference for each view;

a3) estimating a projection error on the detector cell with performance difference for each view from the reconstructed one or more initial images having streaks; and

a4) reconstructing one or more output images with the estimated projection error on the detector cell with performance difference for each view to remove the streaks from the reconstructed one or more output images.

[0008] In one embodiment of the present invention, in the method as mentioned above, the step al) comprises estimating a projection on a detector cell with performance difference using an interpolation method.

[0009] In the method as mentioned above, the step a4) comprises:

subtracting the estimated projection error on the detector cell with performance difference from the estimated projection on the detector cell with performance difference to obtain a corrected projection on the detector cell with performance difference; and

reconstructing one or more output images with the corrected projection on the detector cell with performance difference.

[0010] In one embodiment of the present invention, in the method as mentioned above, the step a2) comprises reconstructing one or more initial images with the estimated projection on the detector cell with performance difference using a filter backprojection method, and the step a4) comprises reconstructing one or more output images with the estimated projection error on the detector cell with performance difference error using a filter backprojection method.

[0011] In one embodiment of the present invention, in the method as mentioned above, the step a3) comprises:

extracting an image stripe along a projection ray of the detector cell with performance difference from the reconstructed one or more initial images; and

estimating the projection error on the detector cell with performance difference using one or more gray values of one or more pixels in the extracted image stripe.

[0012] In one embodiment of the present invention, in the method as mentioned above, the step a3) further comprises:

forming a first straight line between a position of an X-ray source and a position of the detector cell with performance difference;

selecting second and third straight lines parallel to the first straight line, wherein a distance from the second straight line to the first straight line is the same as a distance from the third straight line to the first straight line; and

estimating the projection error on the detector cell with performance difference using one or more gray values of one or more corresponding pixels of the first, the second and the third straight lines.

[0013] In one embodiment of the present invention, in the method as mentioned above, the step a3) further comprises:

dividing the first straight line into M segments;

dividing respectively the second and the third straight lines into corresponding M segments as well;

computing the projection error on the detector cell with performance difference using the following equation:

wherein represents the projection error on the detector cell with performance difference, represents a gray value of the segment on the first straight line, represents a gray value of the segment on the second straight line, represents a gray value of the segment on the third straight line, M represents the number of the segments divided, represents a constant mapping a streak level to the

projection error on the detector cell with performance difference, and a represents the number of non-zero elements in

[0014] In one embodiment of the present invention, in the method as mentioned above, the step a3) further comprises:

determining whether the segment on individual of the first, the second and the third straight lines is in a pixel range of the reconstructed one or more initial images;

when it is determined that the segment on individual of the first, the second and the third straight lines is not in the pixel range of the reconstructed one or more initial images, the gray value of the segment on individual of the first, the second and the third straight lines is zero; and

when it is determined that the segment on individual of the first, the second and the third straight lines is in the pixel range of the reconstructed one or more initial images, the gray value of the segment on individual of the first, the second and the third straight lines is interpolated from the reconstructed one or more initial images.

[0015] In one embodiment of the present invention, in the method as mentioned above, the distance from the second straight line to the first straight line and the distance from the third straight line to the first straight line are determined by a size of a detector cell of a detector, a distance from the X-ray source to the detector, a distance from the X-ray source to a rotating center of the X-ray source and a reconstruction kernel.

[0016] The method for removing streaks from detector cells with performance difference in accordance with the above embodiments of the present invention can reduce an estimation error by estimating a projection error on the detector cell with performance difference from the reconstructed one or more initial images having streaks, subtracting the estimated projection error on the detector cell with performance difference from the estimated projection on the detector cell with performance difference to further obtain a corrected projection on the detector cell with performance difference, and can remove the streaks from the reconstructed one or more output images by reconstructing one or more output images with the corrected projection on the detector cell with performance difference, improving quality of CT reconstructed images.

[0017] Another aspect of the present invention is to provide a method for removing streaks from detector cells with performance difference, the method comprising the steps of: bl) estimating a projection on a detector cell with performance difference for each of one or more views;

b2) reconstructing one or more initial images with the estimated projection on the detector cell with performance difference for each view;

b3) computing a possible streak direction for each of one or more pixels in the reconstructed one or more initial images;

b4) identifying a streak feature for each pixel along the possible streak direction;

b5) when identifying a pixel is a streak pixel, smoothing the streak pixel in the streak direction so as to smooth a difference image; and

b6) updating one or more output images with the smoothed difference image.

[0018] In one embodiment of the present invention, in the method as mentioned above, the step b3) comprises:

forming a first straight line by an X-ray source and each of one or more pixels in the reconstructed one or more initial images; and

defining the possible streak direction for each of one or more pixels in the reconstructed one or more initial images by the first straight line.

[0019] In one embodiment of the present invention, in the method as mentioned above, the step b4) comprises:

selecting several pixel points near the pixel along the possible streak direction in the reconstructed one or more initial images;

selecting a second straight line passing each of the several pixel points and perpendicular to the possible streak direction;

for each of the several pixel points, finding two pixel points adjacent to the pixel point in the second straight line;

comparing a gray value of the pixel point to gray values of its adjacent two pixel points; and

identifying a streak feature for the pixel using the compared result.

[0020] In one embodiment of the present invention, in the method as mentioned above, the step b4) further comprises:

if the gray value of the pixel point is always smaller than the gray values of its adjacent two pixel points for each of the several pixel points, regarding the pixel as a streak pixel; or

if the gray value of the pixel point is always bigger than the gray values of its adjacent two pixel points for each of the several pixel points, regarding the pixel as a streak pixel.

[0021] In one embodiment of the present invention, in the method as mentioned above, the step b5) comprises:

smoothing the streak pixel in the streak direction using the following equation:

wherein represent weighting factors, each of the reconstructed one or more initial images is an image, represents a gray value of the streak pixel respectively represent gray values of the two pixels adjacent to the streak pixel in a second straight line passing the streak pixel and perpendicular to the streak direction, and represents a correct gray value of the streak pixel

[0022] In one embodiment of the present invention, in the method as mentioned above, the step b6) comprises:

obtaining an output image using the following equation:

wherein represents an output image, represents a reconstructed initial image, represents a corrected image obtained from smooth represents a smooth operator and represents a parameter for mapping a difference level to after the smooth operator.

[0023] The method for removing streaks from detector cells with performance difference in accordance with the above embodiments of the present invention can remove the streaks from the reconstructed one or more initial images directly by identifying pixels in the streak direction from the reconstructed one or more initial images having streaks and smoothing the streak pixels, thus improving quality of CT reconstructed images.

[0024] A further aspect of the present invention is to provide a method for removing streaks from detector cells with performance difference, the method comprising the steps of: cl) estimating a projection on a detector cell with performance difference for each of one or more views when an X-ray source is located in each of one or more positions;

c2) estimating a projection error on a conjugate detector cell for a conjugate view when the X-ray source is located in a conjugate position relative to the each position, wherein the each position and its conjugate position are located in a connecting line between the X-ray source and the detector cell with performance difference when the X-ray source is located in the each position; and

c3) reconstructing one or more output images using the estimated projection on the detector cell with performance difference and the estimated projection error on the conjugate detector cell.

[0025] In one embodiment of the present invention, in the method as mentioned above, the step cl) comprises:

computing projections on detector cells adjacent to the detector cell with performance difference; and

estimating a projection on the detector cell with performance difference with the computed projections on the detector cells adjacent to the detector cell with performance difference using an interpolation method.

[0026] In one embodiment of the present invention, in the method as mentioned above, the step c2) comprises:

c21) computing a conjugate projection on the conjugate detector cell for a conjugate view when the X-ray source is located in a conjugate position relative to the each position;

c22) estimating a conjugate projection on the conjugate detector cell for the conjugate view when the X-ray source is located in the conjugate position relative to the each position; and

c23) subtracting the computed conjugate projection from the estimated conjugate projection to estimate the projection error, on the conjugate detector cell.

[0027] In one embodiment of the present invention, in the method as mentioned above, the step c22) further comprises:

computing projections on detector cells adjacent to the conjugate detector cell; and estimating a conjugate projection on the conjugate detector cell with the computed projections on the detector cells adjacent to the conjugate detector cell using an interpolation method.

[0028] In one embodiment of the present invention, in the method as mentioned above, the step c3) comprises:

subtracting the estimated projection error on the conjugate detector cell from the estimated projection on the detector cell with performance difference to obtain an improved projection; and

reconstructing one or more output images with the improved projection.

[0029] The method for removing streaks from detector cells with performance difference in accordance with the above embodiments of the present invention can remove the streaks from detector cells with performance difference from the reconstructed one or more output images by utilizing a conjugate projection ray of the X-ray source and utilizing a conjugate projection ray sensed by a conjugate detector cell to the detector cell with performance difference, estimating an estimation error that is possible to occur on the detector cell with performance difference with a projection error that is possible to occur on the conjugate detector cell, reducing the estimation error caused by low sampling of the detector cell with performance difference, and reconstructing one or more output images with the improved projection, thus improving quality of CT reconstructed images.

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] These and other features, aspects and advantages of the present invention can be better understood when the following detailed description is read with reference to the accompanying drawings, in which the same reference signs represent the same components in the whole drawings, in which:

Fig. 1 is a schematic diagram of one type of CT system when detecting a patient;

Fig. 2 is a flow chart of a method for removing streaks from detector cells with performance difference in accordance with a first embodiment of the present invention;

Fig. 3 is a schematic diagram used by the method as shown in Fig. 2 when estimating a projection error on a detector cell with performance difference;

Fig. 4 is a flow chart of a method for removing streaks from detector cells with performance difference in accordance with a second embodiment of the present invention;

Fig. 5 is a schematic diagram used by the method as shown in Fig. 4 when estimating a possible streak direction; Fig. 6 is a flow chart of a method for removing streaks from detector cells with performance difference in accordance with a third embodiment of the present invention;

Fig. 7a shows a simplified illustration when an X-ray source is located in one position; and

Fig. 7b shows a simplified illustration when the X-ray source is located in a conjugate position of the position as shown in Fig. 7a.

DETAILED DESCRIPTION

[0031] In order to help the person skilled in the art to exactly understand the subject matters claimed by the present invention, embodiments of the present invention will be described in detail hereinafter with reference to the accompanying drawings. In the following detailed description for those embodiments, some known functions or structures will not be described in detail in the Description, to avoid disclosure of the present invention to be affected by unnecessary details.

[0032] Unless defined otherwise, the technical or scientific terms used in the Claims and the Description should have meanings as commonly understood by one of ordinary skilled in the art to which the present disclosure belongs. The terms "first", "second" and the like in the Description and the Claims do not mean any sequential order, quantity or importance, but are only used for distinguishing different components. The terms "a", "an" and the like do not denote a limitation of quantity, but denote the existence of at least one. The terms "comprises", "comprising", "includes", "including" and the like mean that the element or object in front of the "comprises", "comprising", "includes" and "including" covers the elements or objects and their equivalents illustrated following the "comprises", "comprising", "includes" and "including", but do not exclude other elements or objects. The term "coupled" or "connected" or the like is not limited to being connected physically or mechanically, but may comprise electric connection, no matter directly or indirectly.

[0033] Fig. 1 illustrates a schematic diagram of one type of CT system when detecting a patient. As shown in Fig. 1, the CT system 10 generally comprises a rotatable gantry 12 and a supporting table 15 disposed in a hollow imaging area 14 in the rotatable gantry 12 for supporting a patient 30. The rotatable gantry 12 includes an X-ray source S and a detector

18 disposed opposite to the X-ray source S, in which the detector 18 includes a plurality of independent detector cells D arranged in an array. When the rotatable gantry 12 is located in a certain scanning position, the X-ray source S emits a sector-shaped beam of X-rays 162 towards a direction of the detector 18, and the plurality of detector cells D sense the X-rays attenuated through the patient 30 respectively, thus a set of projections are sensed by the detector cells D to obtain a corresponding view. As the rotatable gantry 12 rotates, the X-ray source S and the detector 18 rotate along with it around a rotating center O. The CT system 10 carries out a plurality of scans. During each of the plurality of scans, each corresponding view may be sensed by all the detector cells D. In the case that the detector cells D are normal, one or more images may be reconstructed directly using each corresponding view. However, when there is some performance difference between the detector cells D on the detector 18, i.e., in the case that the detector 18 has a detector cell with performance difference (in which ° is a column number of the detector cell wlth performance difference in the array of the detector 18, is a row number of the detector cell with performance difference in the array of the detector 18), the projection for each view sensed by the detector cell with performance difference cannot correctly reflect soft tissue features of the patient 30, therefore, the projection for each view sensed by the detector cell with performance difference cannot be used. It is needed to estimate a projection on the detector cell with performance difference for each view as accurately as possible by the following method of the present invention.

[0034] Fig. 2 illustrates a flow chart of a method for removing streaks from detector cells with performance difference in accordance with a first embodiment of the present invention.Now referring to Fig. 2, the method for removing streaks from detector cells with performance difference in accordance with a first embodiment of the present invention comprises the following steps:

[0035] In Step sll, a projection on a detector cell with performance difference for each of one or more views is estimated, in which is a column number of the detector cell with performance difference, is a row number of the detector cell with performance difference, is a view number of the projection. In the present embodiment, the projection _on the detector cell

D(n, row) with performance difference may be estimated using a known interpolation method.The known interpolation method may be any suitable interpolation method.For example, the projection on the detector cell with performance difference is estimated with the projections measured on two detector cells adjacent to the detector cell with performance difference by a known interpolation method.

[0036] In Step sl2, one or more initial images I are reconstructed with the projection on the detector cell with performance difference for each view estimated through Step sl l . In the present embodiment, one or more initial images I may be reconstructed with the estimated projection on the detector cell with performance difference using a filter backprojection method. Since there is an estimation error between the estimated projection of the detector cell

with performance difference in Step s ll and a correct projection of the detector cell with performance difference, there exist streaks in the one or more initial images I reconstructed in Step si 2, and the estimation error is just the

reason that causes the existence of streaks on the reconstructed one or more initial images I. Therefore, if the estimation error can be reduced to a very small value, even approaching zero, then the streaks on the reconstructed one or more initial images I can be removed. Accordingly, the following steps are just to attempt to reduce the estimation error so as to remove the streaks on the reconstructed one or more initial images I.

[0037] In Step sl3, a projection error on the detector cell with performance difference for each view is estimated from the reconstructed one or more initial images I having streaks.

[0038] In one embodiment of the present invention, Step s l3 comprises:

Step s l31, extracting an image stripe along a projection ray of the detector cell D(n, row) with performance difference from the reconstructed one or more initial images I; and

Step si 32, estimating the projection error on the detector cell with performance difference using one or more gray values of one or more pixels in the extracted image stripe.

[0039] Referring to Fig. 3, a first straight line SD is formed between the position of the X-ray source S and the position of the detector cell with performance difference. Two straight lines parallel to the first straight line SD are selected at two sides respectively adjacent to the first straight line SD, which are referred to as a second straight line LI and a third straight line L2 respectively, and a distance wl from the second straight line LI to the first straight line SD is equal to a distance w2 from the third straight line L2 to the first straight line SD. Specifically, the distance wl from the second straight line LI to the first straight line SD and the distance w2 from the third straight line L2 to the first straight line SD are determined by a size of a detector cell D of the detector 18, a distance from the X-ray source S to the detector 18, a distance from the X-ray source S to the rotating center O of the X-ray source S and a reconstruction convolution kernel.

[0040] Specifically, Step si 32 further comprises estimating the projection error on the detector cell with performance difference using one or more gray values of one or more corresponding pixels of the first, the second and the third straight lines SD, LI and L2.

[0041] Continuing to refer to Fig. 3, first, the first straight line SD is divided into M segments, in which SDk in Fig. 3 represents the segment on the first straight line SD. Then, the second straight line LI and the third straight line L2 are respectively divided into corresponding M segments as well, in which in Fig. 3 represents the segment on the second straight line LI, and L2kin Fig. 3 represents the segment on the third straight line L2.

[0042] The projection error on the detector cell with performance difference may be calculated using the following equation:

Wherein represents a gray value of the segment On the first straight line

, represents a gray value of the segment on the second straight line LI, represents a gray value of the segment on the third straight line L2. First, it is needed to determine whether the segments on individual of the first, the second and the third straight lines SD, LI and L2 are in a pixel range of the reconstructed one or more initial images I, wherein the individual positions may be determined based on a distance from the X-ray source S to the detector cell with performance difference, a fan angle Ύ (an angle between the connecting line SD from the X-ray source S to the detector cell with performance difference and a connecting line SO from the X-ray source S to its rotating center O) of the X-ray source S and a projection angle β (an angle between the connecting line SD from the X-ray source S to the detector cell with performance difference and a horizontal x-axis) of the X-ray source S. If it is determined that the segments, on individual of the first, the second and the third straight lines SD, LI and L2 are within the pixel range of the reconstructed one or more initial images I, the gray values segments,

On individual of the first, the second and the third straight lines SD, LI and L2 may be interpolated from the reconstructed one or more initial images I based on the individual positions; if it is determined that the segments, On individual of the first, the second and the third straight lines SD, LI and L2 are not within the pixel range of the reconstructed one or more initial images I, the gray values and of the segments and L2k On individual of the first, the second and the

third straight lines SD, LI and L2 are equal to zero.

[0043] In addition, M represents the number of the segments divided, a represents the number of non-zero elements in in Equation (2), represents a constant mapping a streak level to the projection error on the detector cell with performance difference.

[0044] In Step s l4, one or more output images I0 are reconstructed with the estimated projection error on the detector cell with performance difference for each

view, to remove the streaks from the reconstructed one or more output images I0. In the present embodiment, one or more output images I0 are reconstructed with the estimated projection error on the detector cell with performance difference using a filter backprojection method.

[0045] Furthermore, Step s l4 comprises:

Step s l41, subtracting the estimated projection error on the detector cell with performance difference from the estimated projection on the detector cell with performance difference to obtain a corrected projection on the detector cell with performance difference, i.e., as shown in the following equation:

Step sl42, reconstructing one or more output images I0 with the corrected projection on the detector cell with performance difference obtained through Step sl41.

[0046] The method for removing streaks from detector cells with performance difference of the first embodiment of the present invention can reduce an estimation error by estimating a projection error on the detector cell with performance difference from the reconstructed one or more initial images I having streaks, subtracting the estimated projection error on the detector cell with performance difference from the estimated projection on the detector cell with performance difference to further obtain a corrected projection ) on the detector cell with performance difference, and can remove the streaks from the reconstructed one or more output images I0 by reconstructing one or more output images I0 with the corrected projection P(ii, row, view) on the detector cell with performance difference, improving quality of CT reconstructed images.

[0047] Fig. 4 illustrates a flow chart of a method for removing streaks from detector cells with performance difference in accordance with a second embodiment of the present invention. Now referring to Fig. 4, the method for removing streaks from detector cells with performance difference in accordance with a second embodiment of the present invention comprises the following steps:

[0048] In Step s21, a projection on a detector cell with performance difference for each of one or more views is estimated, in which is a column number of the detector cell with performance difference, is a row number of the detector cell with performance difference, is a view number of the projection. In the present embodiment, the projection on the detector cell D(n, row) with performance difference may be estimated using a known interpolation method. The known interpolation method may be any suitable interpolation method. For example, the projection on the detector cell with performance difference is estimated with the projections sensed on two detector cells adjacent to the detector cell with performance difference by a known interpolation method.

[0049] In Step s22, one or more initial images I are reconstructed with the projection on the detector cell with performance difference for each view estimated through Step sll, each of the reconstructed one or more initial images I being an image. In the present embodiment, one or more initial images I may be reconstructed with the estimated projection on the detector cell with performance difference using a filter backprojection method. Since there is an estimation error

between the estimated projection of the detector cell with performance difference in Step s21 and a correct projection of the detector cell with performance difference, there exist streaks in the reconstructed one or more initial images I in Step s22. The streaks present on the reconstructed one or more initial images I are generally distributed along a straight line, while there is very little straight line structure in images of a human tissue. Therefore, if a straight line structure can be identified from the reconstructed one or more initial images I and pixels on said straight line can be smoothed, then streaks can be removed from the reconstructed one or more initial images I. Accordingly, the following steps are to remove the streaks on the reconstructed one or more initial images I in light of such ideas of invention.

[0050] In Step s23, a possible streak direction for each of one or more pixels in the reconstructed one or more initial images I is calculated.

[0051] Referring to Fig. 5, the position of the X-ray source S and each pixel P(i, j) of one or more pixels in the reconstructed one or more initial images I form a first straight line SP, wherein i=[0, N] , j=[0, N].

[0052] The first straight line SP defines the possible streak direction for each pixel P(i, j) of one or more pixels in the reconstructed one or more initial images I, i.e., the possible streak direction is a direction along the first straight line SP.

[0053] In Step s24, a streak feature for each pixel P(i, j) along the possible streak direction SP is identified.

[0054] Continuing to refer to Fig. 5, first, several pixel points adjacent

to the pixel P(i, j) are selected along the possible streak direction SP in the reconstructed one or more initial images I.

[0055] A straight line is selected to pass each pixel point of the several pixel points and perpendicular to the possible streak direction SP, which is referred to as a second straight line LI.

[0056] For each pixel point of the several pixel points, find two pixel points

PQ _k 2 adjacent to the pixel point in the second straight line LI, compare a gray value of the pixel point PQ _k to gray values of its adjacent two pixel points and identify a

streak feature for the pixel P(i, j) using the compared results.

[0057] In Step s25, identify whether a pixel P(i, j) is a streak pixel? In the case of identifying that the pixel P(i, j) is a streak pixel, the process proceeds to Step s25; in the case of identifying that the pixel P(i, j) is not a streak pixel, the process returns back to Step s24.

[0058] In one embodiment of the present invention, if there is always the gray value of the pixel point PQ _k that is smaller than the gray values of its adjacent two pixel points

PQ _k 2 for each pixel point PQ _k, k=[l, R] of the several pixel points, the pixel P(i, j) is regarded as a streak pixel. Or, alternatively, in another embodiment of the present invention, if there is always the gray value of the pixel point PQ _k that is bigger than the gray values of its adjacent two pixel points for each pixel point PQ _k, k=[l, R] of the several pixel points, the

pixel P(i, j) is regarded as a streak pixel.

[0059] In Step s26, when a pixel P(i, j) is identified as a streak pixel, the streak pixel P(i, j) in the streak direction is smoothed so as to smooth a difference image.

[0060] Referring to Fig. 5, similarly, a second straight line LI is selected in a direction passing the streak pixel P(i, j) and perpendicular to the streak direction SP. Then, two pixels PI, P2 adjacent to the streak pixel P(i, j) are found in the second straight line LI.

[0061] In one embodiment of the present invention, in Step s26, the steak pixel P(i, j) in the streak direction is smoothed using the following equation:

wherein represent weighting factors, represents a gray value of

the streak pixel respectively represent gray values of the two pixels PI,

P2 adjacent to the streak pixel , and represents a correct gray value of the streak pixel

[0062] The positions of each pixel point PQ _k and its adjacent two pixel points and two pixels PI, P2 adjacent to the streak pixel may be determined based on a distance from the X-ray source S to the streak pixel a fan angle Ύ (an angle between the connecting line SP from the X-ray source S to the streak pixel and the connecting line SO from the X-ray source S to its rotating center O) of the X-ray source S and a projection angle β (an angle between the connecting line SP from the X-ray source S to the streak pixel and the positive direction of the horizontal x-axis) of the X-ray source S and the like. The individual gray values may be determined based on the individual positions in the reconstructed one or more initial images I.

[0063] A difference image is obtained by the following equation:

wherein represents a reconstructed initial image, Ϊ represents a corrected image obtained from

[0064] After obtaining the difference image , the difference image is smoothed to remove a high-frequency component from the difference image

[0065] In Step s27, one or more output images are updated with the smoothed difference image, such that the output images will become more natural, without flicker.

[0066] In one embodiment of the present invention, the output image is obtained using the following equation:

wherein represents a smooth operator and represents a parameter for mapping a difference level to the difference image after the smooth operator.

[0067] The method for removing streaks from detector cells with performance difference of the second embodiment of the present invention can remove the streaks from the reconstructed one or more initial images I directly by identifying pixels in the streak direction from the reconstructed one or more initial images I having streaks and smoothing the streak pixels, thus improving quality of CT reconstructed images.

[0068] The projections on the detector cells with performance difference may be estimated using the projections on the adjacent detector cells by a known interpolation method. However, tiny structures in the human body may cause a larger estimation error, which cannot be reduced by an advanced interpolation method. Therefore, the present invention provides a method for removing streaks from detector cells with performance difference of the third embodiment. The reason for causing said estimation error is a low sampling rate due to the detector cell with performance difference. In the present method, the estimation error may be reduced using redundant information from a sine diagram. Since the conjugate projection ray of the X-ray source S is most likely to pass similar tiny structures, and since a taper angle of the beam 162 of X-rays cannot use a conjugate projection of the conjugate projection ray directly, the present method utilizes a high sampling rate of a conjugate detector cell at a conjugate position to generate correct date, thereby obtaining the estimation error caused by the low sampling rate of the detector cell D(n, row) with performance difference by an interpolation method. the method for removing streaks from detector cells with performance difference of the third embodiment of the present invention is described in detail in the following.

[0069] Fig. 6 illustrates a flow chart of the method for removing streaks from detector cells with performance difference in accordance with the third embodiment of the present invention. Now referring to Fig. 6, the method for removing streaks from detector cells with performance difference in accordance with the third embodiment of the present invention comprises the following steps:

[0070] In Step s31, a projection on a detector cell with performance difference for each of one or more views when the X-ray source S is located in each of one or more positions is estimated.

[0071] When the rotatable gantry 12 is located in a certain scanning position, the X-ray source S is located in a position as shown in Fig. 7a, in which projections on detector cells adjacent to the detector cell with performance difference may be calculated from the sensed information on the detector cells adjacent to the detector cell D(n, row) with performance difference. In one embodiment of the present invention, by the sensed information on four detector cells adjacent to the detector cell with performance difference, projections

and on the four detector cells adjacent to the detector cell with performance difference in the position may be calculated respectively. According to different scanning methods, the row numbers of the adjacent four detector cells may be the same as or different from the row number of the detector cell with performance difference.

[0072] Then, using any suitable known interpolation method, the projection on the detector cell with performance difference may be estimated with the calculated projections

and on the detector cells adjacent to the detector cell with performance difference, i.e., as shown in the following equation:

[0073] However, the present invention is not limited thereto. In another embodiment of the present invention, the projection on the detector cell with performance difference may be estimated with the projections on two or more detector cells adjacent to the detector cell with performance difference.

[0074] In Step s32, a projection error on a conjugate detector cell for a conjugate view when the X-ray source S is located in a conjugate position relative to each position is estimated, wherein each position and its conjugate position are located in a connecting line between the X-ray source S and the detector cell with performance difference when the X-ray source S is located in each position.

[0075] As the rotatable gantry 12 rotates, the X-ray source S and the detector 18 rotate around the rotating center O along with it. When the rotatable gantry 12 rotates to a conjugate position of the position of the X-ray source S as shown in Fig. 7a, i.e., the X-ray source S rotates to be located in the conjugate position as shown in Fig. 7b along with it, i.e., the conjugate position of the X-ray source S as shown in Fig. 7b corresponds to the position of the detector cell with performance difference when in the position of the X-ray source S as shown in Fig. 7a. When in the conjugate position of the X-ray source S as shown in Fig. 7b, the detector cell with performance difference rotates from the position as shown in Fig. 7a to the position as shown in Fig. 7b, and at this moment, the conjugate detector cell of the detector cell with performance difference just corresponds to the position of the X-ray source S as shown in Fig. 7a.

[0076] In one embodiment of the present invention, Step s32 comprises:

[0077] In Step s321, a conjugate projection on the conjugate detector cell for a conjugate view when the X-ray source S is located in a conjugate position relative to each position is calculated.

[0078] In the conjugate position as shown in Fig. 7b, the conjugate projection rays emitted from the X-ray source S will be sensed by the conjugate detector cell

of the detector cell with performance difference, and the conjugate projection on the conjugate detector cell may be calculated directly from the sensed information of the conjugate detector cell

[0079] In Step s322, a conjugate projection on the conjugate detector cell for a conjugate view when the X-ray source S is located in a conjugate position relative to each position is estimated.

[0080] When the X-ray source S is in the conjugate position as shown in Fig. 7b, the projections on the detector cells adjacent to the conjugate detector cell may be calculated by the sensed information of the detector cell in the conjugate position corresponding to the detector cells adjacent to the detector cell D(n, row) with performance difference. In one embodiment of the present invention, by the four detector cells adjacent to the conjugate detector cell D(n, row) , projections , and οn the four detector cells adjacent to the conjugate detector cell in the position may be computed respectively. According to different scanning methods, the row numbers , and of the

adjacent four detector cells may be the same as or different from the row number rowc of the conjugate detector cell

[0081] Then, using any suitable known interpolation method, the conjugate projection on the conjugate detector cell may be estimated with the computed projections

and ,

on the detector cells adjacent to the conjugate detector cell ) i.e., as shown in the following equation:

[0082] However, the present invention is not limited thereto. In another embodiment of the present invention the conj ·ugate projection on the conjugate detector cell may be estimated with the projections on two or more detector cells adjacent to the conjugate detector cell a known interpolation method.

[0083] In Step s323, the computed conjugate projection on the conjugate detector cell is subtracted from the estimated conjugate projection on the conjugate detector cell to estimate the projection error on the conjugate detector cell , i. e., as shown in the following equation:

[0084] In Step s33, one or more output images are reconstructed using the estimated projection on the detector cell with performance difference and the estimated projection error on the conjugate detector cell

[0085] Furthermore, Step s33 comprises:

Step s331, subtracting the estimated projection error on the conjugate detector cell from the estimated projection on the detector cell with performance difference to obtain an improved projection , i.e., as shown in the following equation: Step s332, reconstructing one or more output images with the improved projection.

[0086] The method for removing streaks from detector cells with performance difference of the third embodiment of the present invention can remove the streaks from the detector cells with performance difference from the reconstructed one or more output images Icby utilizing a conjugate projection ray of the X-ray source S and utilizing a conjugate projection ray sensed by a conjugate detector cell to the detector cell wlth performance difference, estimating an estimation error that is possible to occur on the detector cell with performance difference with a projection error that is possible to occur on the conjugate detector cell reducing the estimation error caused by low sampling of the detector cell with performance difference, and reconstructing one or more output images with the improved projection thus improving quality of CT reconstructed

images.

[0087] Although the present invention has been set forth in details in combination with specific embodiments, the person skilled in the art shall be understood that many modifications and variations may be made to the present invention. Therefore, it should be recognized that the intention of the claims is to cover all these modifications and variations within the real concept and range of the present invention.