Information
-
Patent Grant
-
6373919
-
Patent Number
6,373,919
-
Date Filed
Monday, June 12, 200024 years ago
-
Date Issued
Tuesday, April 16, 200222 years ago
-
Inventors
-
Original Assignees
-
Examiners
Agents
-
CPC
-
US Classifications
Field of Search
US
- 378 4
- 378 62
- 378 98
- 378 901
-
International Classifications
-
Abstract
In order to provide a projection data correction method and apparatus that prevents artifact generation in the direction which provides a low count of penetrating radiation, and a radiation tomographic imaging method employing such a projection data correction method and a radiation tomographic imaging apparatus comprising the projection data correction apparatus, projection data Projij is individually corrected by a weighted addition (604, 606, 608) with an average value (602) containing the proximate data. The weight SGij of the projection data is decreased as the count of the penetrating radiation from which the projection data is derived becomes smaller.
Description
BACKGROUND OF THE INVENTION
The present invention relates to a projection data correction method and apparatus and a radiation tomographic imaging method and apparatus, and more particularly to a method and apparatus for individually correcting data constituting a projection of an object to be imaged by penetrating radiation, and a radiation tomographic imaging method employing such a projection data correction method and a radiation tomographic imaging apparatus comprising the projection data correction apparatus.
In a radiation tomographic imaging apparatus employing an X-ray, i.e., an X-ray CT (computerized tomography) apparatus, an X-ray beam having a width containing a range to be imaged and having a thickness in a direction orthogonal to the width is emitted from an X-ray emitting apparatus. The thickness of the X-ray beam can be changed by regulating the openness of an X-ray passing aperture of a collimator, and the slice thickness for imaging is thus adjusted.
An X-ray detecting apparatus comprises a multi-channel X-ray detector to detect an X-ray, and the detector has a multiplicity (e.g., of the order of 1,000) of X-ray detector elements arranged in an array in the direction of the X-ray width.
The X-ray emitting/detecting apparatus is rotated (or scans) around an object to be imaged to generate X-ray projection data of the object in a plurality of view directions around the object, and a tomographic image is produced (or reconstructed) by a computer based on the projection data.
The SNR (signal-to-noise ratio) of the projection data is varied with the “count” of the penetrating X-ray, and a smaller count gives a lower SNR. If the object to be imaged provides significantly different counts, i.e., significantly different SNR's of the projection data, depending on the view direction, the reconstructed image suffers from false linear images, or so-called streaking artifacts, due to the effect of a noise pattern in the direction which provides a poor SNR.
Thus, when the thorax or abdomen is imaged with both the arms put on the sides of the body, for example, the streaking artifacts inevitably appear in the direction which connects the arms, leading to a problem that an image having a good quality cannot be obtained.
SUMMARY OF THE INVENTION
The present invention is directed to solving the aforementioned problem, and the object of the invention is to provide a projection data correction method and apparatus that prevents artifact generation in the direction which provides a low count of penetrating radiation, and a radiation tomographic imaging method employing such a projection data correction method and a radiation tomographic imaging apparatus comprising the projection data correction apparatus.
(1) In accordance with a first aspect of the invention to solve the aforementioned problem, there is provided a projection data correction method comprising the steps of: individually calculating an average value of data and its proximate data, the data constituting a projection of an object to be imaged by penetrating radiation; calculating a first product by multiplying the data by a first weighting factor that corresponds to a count of the penetrating radiation from which the data is derived; calculating a second product by multiplying the average value by a second weighting factor that is a complement of the first weighting factor with respect to one; and adding the first and second products.
(2) In accordance with a second aspect of the invention to solve the aforementioned problem, there is provided a projection data correction apparatus comprising: average value calculating means for individually calculating an average value of data and its proximate data, said data constituting a projection of an object to be imaged by penetrating radiation; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; and adder means for adding said first and second products.
(3) In accordance with a third aspect of the invention to solve the aforementioned problem, there is provided a radiation tomographic imaging apparatus comprising: radiation emitting/detecting means for emitting radiation toward an object to be imaged in a plurality of view directions and detecting penetrating radiation, and image producing means for producing an image using a projection based on signals detected by said radiation emitting/detecting means, wherein said image producing means comprises: average value calculating means for individually calculating an average value of data and its proximate data, said data constituting said projection; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; adder means for adding said first and second products; and image reconstruction means for reconstructing an image using a projection comprised of data resulted from said addition.
(4) In accordance with a fourth aspect of the invention to solve the aforementioned problem, there is provided the radiation tomographic imaging apparatus as described in (3), wherein an X-ray is employed as the radiation.
In any one of (1)-(4), it is preferred that the first weighting factor be given by the following equation:
1−e
(−C/α)
,
insofar as artifacts are effectively prevented, wherein
C: a count of the penetrating radiation from which the data is derived, and
α: a constant determined by the area of the projection. (EFFECT)
According to the present invention, the projection data is individually corrected by a weighted addition with an average value containing the proximate data. The weight of the data is made to correspond to the count of penetrating radiation from which the data is derived. Thus, data with a lower count has a smaller weight of that data and a greater weight of the average value, which can compensate for reduction in SNR.
Thus, the present invention can provide a projection data correction method and apparatus that prevents artifact generation in the direction which provides a low count of penetrating radiation, and a radiation tomographic imaging apparatus comprising such a projection data correction apparatus.
Further objects and advantages of the present invention will be apparent from the following description of the preferred embodiments of the invention as illustrated in the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1
is a block diagram of an apparatus in accordance with an embodiment of the present invention.
FIG. 2
is a schematic diagram showing the detector array in the apparatus shown in FIG.
1
.
FIGS. 3 and 4
are schematic diagrams showing the X-ray emitting/detecting apparatus in the apparatus shown in FIG.
1
.
FIG. 5
is a block diagram of the central processing apparatus with respect to projection data correction function.
FIG. 6
illustrates a concept of a projection of an object to be imaged.
FIG. 7
is a graphical representation illustrating a relationship between a count of a penetrating X-ray and a weighting factor.
FIG. 8
is a block diagram of the central processing apparatus with respect to the projection data correction function.
FIGS. 9 and 10
are flow charts of the operation of the apparatus shown in FIG.
1
.
DETAILED DESCRIPTION OF THE INVENTION
Several embodiments of the present invention will now be described in more detail with reference to the accompanying drawings.
FIG. 1
shows a block diagram of an X-ray CT apparatus, which is an embodiment of the present invention. The configuration of the apparatus represents an embodiment of the apparatus in accordance with the present invention, and the operation of the apparatus represents an embodiment of the method in accordance with the present invention.
As shown in
FIG. 1
, the apparatus comprises a scan gantry
2
, an imaging table
4
and an operator console
6
. The scan gantry
2
is an embodiment of the radiation emitting/detecting means of the present invention. The scan gantry
2
has an X-ray tube
20
. An X-ray (not shown) emitted from the X-ray tube
20
is formed into, for example, a fan-shaped X-ray beam, i.e., a fan beam, by a collimator
22
, and cast upon a detector array
24
. The detector array
24
has a plurality of X-ray detector elements arranged in an array in the direction of the fan-shaped X-ray beam spread. The configuration of the detector array
24
will be described later.
The X-ray tube
20
, collimator
22
and detector array
24
together constitute an X-ray emitting/detecting apparatus, which will be described later. The detector array
24
is connected with a data acquisition section
26
for acquiring data detected by the individual X-ray detector elements in the detector array
24
.
The emission of the X-ray from the X-ray tube
20
is controlled by an X-ray controller
28
. The connection relationship between the X-ray tube
20
and the X-ray controller
28
is omitted in the drawing. The collimator
22
is controlled by a collimator controller
30
. The connection relationship between the collimator
22
and the collimator controller
30
is omitted in the drawing.
The above-described components from the X-ray tube
20
through the collimator controller
30
are supported on a rotating section
34
of the scan gantry
2
. The rotation of the rotating section
34
is controlled by a rotation controller
36
. The connection relationship between the rotating section
34
and the rotation controller
36
is omitted in the drawing.
The imaging table
4
is configured to carry an object to be imaged (not shown in
FIG. 1
) into and out of an X-ray irradiation space in the scan gantry
2
. The relationship between the object to be imaged and the X-ray irradiation space will be described later.
The operator console
6
has a central processing apparatus
60
, which is comprised of, for example, a computer. The central processing apparatus
60
is connected with a control interface
62
, which is in turn connected with the scan gantry
2
and the imaging table
4
. The central processing apparatus
60
controls the scan gantry
2
and the imaging table
4
via the control interface
62
.
The data acquisition section
26
, X-ray controller
28
, collimator controller
30
and rotation controller
36
in the scan gantry
2
are controlled via the control interface
62
. The individual connections between these sections and the control interface
62
are omitted in the drawing.
The central processing apparatus
60
is also connected with a data acquisition buffer
64
, which is in turn connected with the data acquisition section
26
in the scan gantry
2
. Data acquired at the data acquisition section
26
is input to the data acquisition buffer
64
, and the buffer
64
temporarily stores the input data.
The central processing apparatus
60
performs data correction for projections of a plurality of views acquired via the data acquisition buffer
64
, and performs image reconstruction using the corrected projection data. The central processing apparatus
60
is an embodiment of the image producing means, and is also an embodiment of the image reconstruction means of the present invention. The image reconstruction is performed using a filtered back projection technique, for example.
The central processing apparatus
60
is also connected with a storage device
66
for storing several data, reconstructed images, programs and so forth. The central processing apparatus
60
is moreover connected with a display device
68
that displays the reconstructed image and other information output from the central processing apparatus
60
, and an operating device
70
that is operated by a human operator to supply several instructions and information to the central processing apparatus
60
.
FIG. 2
schematically illustrates a configuration of the detector array
24
. The detector array
24
is a multi-channel X-ray detector having a multiplicity of X-ray detector elements
24
(
i
) arranged in a row. The multiplicity of the X-ray detector elements
24
(
i
) together form an X-ray impinging surface curved into a cylindrical concave surface. Reference symbol ‘i’ designates a channel index and ‘i’=1-1,000, for example.
Each X-ray detector element
24
(
i
) is formed of a combination of a scintillator and a photodiode, for example. It will be easily recognized that the X-ray detector element
24
(
i
) is not limited thereto but may be a semiconductor X-ray detector element using cadmium telluride (CdTe) or the like, or an ionization chamber X-ray detector element using xenon (Xe) gas, for example.
FIG. 3
illustrates a relationship among the X-ray tube
20
, collimator
22
and detector array
24
in the X-ray emitting/detecting apparatus. FIG.
3
(
a
) is a view from the front of the scan gantry
2
and (b) is a view from the side of the scan gantry
2
. As shown, the X-ray emitted from the X-ray tube
20
is formed into a fan-shaped X-ray beam
400
by the collimator
22
, and cast upon the detector array
24
.
In FIG.
3
(
a
), the spread of the fan-shaped X-ray beam
400
, i.e., the width of the X-ray beam
400
is illustrated. The direction of the width of the X-ray beam
400
is identical to the direction of the channel row in the detector array
24
. In FIG.
3
(
b
), the thickness of the X-ray beam
400
is illustrated.
As exemplarily shown in
FIG. 4
, an object to be imaged
8
is rested on the imaging table
4
, and is carried into the X-ray irradiation space with the object's body axis intersecting the fan surface of the X-ray beam
400
. The scan gantry
2
has a cylindrical structure containing the X-ray emitting/detecting apparatus therein.
The X-ray irradiation space is formed in the internal space of the cylindrical structure of the scan gantry
2
. An image of the object to be imaged
8
sliced by the X-ray beam
400
is projected on the detector array
24
, and the detector array
24
detects the X-ray after passing through the object
8
. The slice thickness ‘th’ of the X-ray beam
400
impinging upon the object to be imaged
8
is regulated by the openness of an aperture of the collimator
22
.
FIG. 5
shows a block diagram of the central processing apparatus
60
with respect to the projection data correction. The respective function of each block shown in
FIG. 5
is implemented by, for example, a computer program. The central processing apparatus
60
shown in
FIG. 5
is an embodiment of the projection data correction apparatus of the present invention. The configuration of the apparatus represents an embodiment of the apparatus in accordance with the present invention, and the operation of the apparatus represents an embodiment of the method in accordance with the present invention.
As shown, the central processing apparatus
60
has a moving averaging unit
602
, which is an embodiment of the average value calculating means of the present invention. The moving average unit
602
individually calculates a moving average value in the direction of the channel ‘i’ for the projection data Proj
ij
(wherein ‘i’=channel index and ‘j’=view index) constituting a projection of the object to be imaged
8
as exemplarily shown in FIG.
6
. Data for use in the moving averaging are of proximate seven channels, for example.
The projection data Proj
ij
is multiplied by a factor SG
ij
at a multiplier
604
, which is an embodiment of the first product calculating means of the present invention. The factor SG
ij
is given by the following equation:
SG
ij
=1−e
(−C
ij
/α)
, (1)
wherein C
ij
is a count of the penetrating X-ray from which the projection data Proj
ij
is derived, and a is a function of the area of the projection ‘s’, which is given by the following equation:
α=f(s). (2)
A relationship between the count C
ij
and the factor SG
ij
by Eq. (1) is graphically shown in
FIG. 7. A
certain range having a small count C
ij
gives a factor SG
ij
that is approximately proportional to C
ij
. Moreover, α is a function of the area of the projection ‘s’, and is varied with the area of the projection. Thus, a factor can be obtained which is adapted for the cross-sectional area of the object to be imaged.
The moving average value of the projection data Proj
ij
is multiplied by a factor 1−SG
ij
at a multiplier
606
, which is an embodiment of the second product calculating means of the present invention. Then, an adder
608
calculates the sum of output data from the multipliers
604
and
606
. The adder
608
is an embodiment of the adder means of the present invention.
As a whole, the block diagram shown in
FIG. 5
performs a calculation as follows:
Proj′
ij
=SG
ij
·Proj
ij
+(1−SG
ij
)SM7{circle around (x)}Proj
ij
, (3)
wherein
SM7{circle around (x)}Proj
ij
represents moving averaging of data of seven channels proximate to the projection data Proj
ij
.
Equation (3) can be rewritten as follows:
Proj′
ij
=SM7{circle around (x)}Proj
ij
+SG
ij
(Proj
ij
−SM7{circle around (x)}Proj
ij
. (4)
A block diagram of an apparatus for performing the calculation given by Eq. (4) is shown in
FIG. 8
, in which similar parts to those shown in
FIG. 5
are denoted by similar symbols and description thereof will be omitted hereinafter. A subtracter
612
calculates the difference between the projection data Proj
ij
and its moving average value. The difference is multiplied by the factor SG
ij
at a multiplier
614
. An adder
616
calculates the sum of output data from the multiplier
614
and the moving averaging unit
602
.
Both of Equations (3) and (4) indicate that a weighted addition is individually performed on the projection data and a moving average value containing the proximate data. The weight of the projection data in the result of the addition is decreased as the count of the penetrating X-ray from which the projection data is derived becomes smaller.
On the contrary, the weight of the moving average value is complementarily increased. Therefore, the weight of the moving average value is increased as the count becomes smaller. Thus, although projection data based on a small count has a poor SNR, the reduction in SNR can be compensated for by increasing the weight of the moving average value for that projection data as above.
Now the operation of the present apparatus will be described.
FIG. 9
shows a flow chart of the operation of the present apparatus. As shown, a human operator inputs a scan plan via the operating device
70
at Step
912
. The scan plan includes the condition under which the X-ray is emitted, slice thickness, slice position and so forth. Thereafter, the present apparatus operates under control of the operation by the human operator and the central processing apparatus
60
following the input scan plan.
At Step
914
, scan positioning is performed. Specifically, the operator manipulates the operating device
70
to move the imaging table
4
and align the center of a site to be imaged of the object
8
with the center of rotation (i.e., isocenter) of the X-ray emitting/detecting apparatus.
After the scan positioning, a scan is performed at Step
916
. Specifically, the X-ray emitting/detecting apparatus is rotated around the object to be imaged
8
and collects projections for, for example, 1,000 views per rotation, into the data acquisition buffer
64
. After or along with the scan, data correction and image reconstruction are performed at Steps
918
and
920
, respectively.
A detailed flow chart of the data correction at Step
918
is shown in FIG.
10
. As shown, an offset correction is performed on the scan data at Step
932
for correcting an offset appended to the count of the penetrating X-ray.
Next, a reference correction is made at Step
934
for correcting the X-ray intensity by dividing the count of the penetrating X-ray by a count of the X-ray measured at a reference channel.
Then, a factor SG
ij
is calculated at Step
936
. The factor SG
ij
is calculated by Eq. (1) using the offset-corrected and reference-corrected count C
ij
, and is saved in the storage device
66
.
Next, a negative logarithm of the count C
ij
is calculated at Step
938
, a sensitivity correction is performed at Step
940
, and a beam-hardening correction is performed at Step
942
to generate projection data Proj
ij
.
The projection data correction by aforementioned Eq. (3) or (4) is performed on the projection data at Step
944
. The factor SG
ij
for the correction calculation is read from the storage device
66
for use.
The image reconstruction at Step
920
is performed by, for example, a filtered back projection technique using projections for a plurality of views comprised of such corrected data.
Since the projection data correction at Step
944
compensates for SNR reduction in a portion having a low count, a directional noise pattern will not occur when the thorax or abdomen is imaged with both the arms put on the sides of the body, thereby preventing false images such as streaking artifacts from being contained in the reconstructed image. The reconstructed tomographic image is displayed on the display device
68
at Step
922
. The display image is a tomographic image having a good quality not containing streaking artifacts or the like.
While the preceding description is made on a case in which an X-ray is employed as radiation, the radiation is not limited to the X-ray but may be other types of radiation such as a y-ray. However, the X-ray is preferred at the present time insofar as practical means for generation, detection, control and so forth of the X-ray are most well-improved.
Many widely different embodiments of the invention may be configured without departing from the spirit and the scope of the present invention. It should be understood that the present invention is not limited to the specific embodiments described in the specification, except as defined in the appended claims.
Claims
- 1. A projection data correction method comprising the steps of:individually calculating an average value of data and its proximate data, said data constituting a projection of an object to be imaged by penetrating radiation; calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; and adding said first and second products.
- 2. A projection data correction method comprising the steps of:individually calculating an average value of data and its proximate data, said data constituting a projection of an object to be imaged by penetrating radiation; calculating a first product by multiplying said data by a first weighting factor that is approximately proportional to a count of the penetrating radiation from which said data is derived in a range below a predetermined value of the count; calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; and adding said first and second products.
- 3. A projection data correction method as defined in claim 1 or 2, wherein said first weighting factor is given by the following equation:1−e(−C/α), whereinC: a count of the penetrating radiation from which said data is derived, and α: a constant determined by the area of said projection.
- 4. A projection data correction apparatus comprising:average value calculating means for individually calculating an average value of data and its proximate data, said data constituting a projection of an object to be imaged by penetrating radiation; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; and adder means for adding said first and second products.
- 5. A projection data correction apparatus comprising:average value calculating means for individually calculating an average value of data and its proximate data, said data constituting a projection of an object to be imaged by penetrating radiation; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that is approximately proportional to a count of the penetrating radiation from which said data is derived in a range below a predetermined value of the count; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; and adder means for adding said first and second products.
- 6. A projection data correction apparatus as defined in claim 4 or 5, wherein said first weighting factor is given by the following equation:1−e(−C/α), whereinC: a count of the penetrating radiation from which said data is derived, and α: a constant determined by the area of said projection.
- 7. A radiation tomographic imaging method for emitting radiation toward an object to be imaged in a plurality of view directions and detecting penetrating radiation, and producing an image using a projection based on signals detected by said step of emitting/detecting radiation, wherein said step of producing an image is performed by:individually calculating an average value of data and its proximate data, said data constituting said projection; calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; adding said first and second products; and reconstructing an image using a projection comprised of data resulted from said addition.
- 8. A radiation tomographic imaging method as defined in claim 7, wherein an X-ray is employed as said radiation.
- 9. A radiation tomographic imaging method for emitting radiation toward an object to be imaged in a plurality of view directions and detecting penetrating radiation, and producing an image using a projection based on signals detected by said step of emitting/detecting radiation, wherein said step of producing an image is performed by:individually calculating an average value of data and its proximate data, said data constituting said projection; calculating a first product by multiplying said data by a first weighting factor that is approximately proportional to a count of penetrating radiation from which said data is derived in a range below a predetermined value of the count; calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; adding said first and second products; and reconstructing an image using a projection comprised of data resulted from said addition.
- 10. A radiation tomographic imaging method as defined in claim 9, wherein an X-ray is employed as said radiation.
- 11. A radiation tomographic imaging method as defined in claim 7, 8, 9 or 10, wherein said first weighting factor is given by the following equation:1−e(−C/α), whereinC: a count of the penetrating radiation from which said data is derived, and α: a constant determined by the area of said projection.
- 12. A radiation tomographic imaging apparatus comprising:radiation emitting/detecting means for emitting radiation toward an object to be imaged in a plurality of view directions and detecting penetrating radiation, and image producing means for producing an image using a projection based on signals detected by said radiation emitting/detecting means, wherein said image producing means comprises: average value calculating means for individually calculating an average value of data and its proximate data, said data constituting said projection; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that corresponds to a count of the penetrating radiation from which said data is derived; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; adder means for adding said first and second products; and image reconstruction means for reconstructing an image using a projection comprised of data resulted from said addition.
- 13. A radiation tomographic imaging apparatus as defined in claim 12, wherein an X-ray is employed as said radiation.
- 14. A radiation tomographic imaging apparatus comprising:radiation emitting/detecting means for emitting radiation toward an object to be imaged in a plurality of view directions and detecting penetrating radiation, and image producing means for producing an image using a projection based on signals detected by said radiation emitting/detecting means, wherein said image producing means comprises: average value calculating means for individually calculating an average value of data and its proximate data, said data constituting said projection; first product calculating means for calculating a first product by multiplying said data by a first weighting factor that is approximately proportional to a count of the penetrating radiation from which said data is derived in a range below a predetermined value of the count; second product calculating means for calculating a second product by multiplying said average value by a second weighting factor that is a complement of said first weighting factor with respect to one; adder means for adding said first and second products; and image reconstruction means for reconstructing an image using a projection comprised of data resulted from said addition.
- 15. A radiation tomographic imaging apparatus as defined in claim 14, wherein an X-ray is employed as said radiation.
- 16. A radiation tomographic imaging apparatus as defined in claim 12, 13, 14 or 15, wherein said first weighting factor is given by the following equation:1−e(−C/α), whereinC: a count of the penetrating radiation from which said data is derived, and α: a constant determined by the area of said projection.
Priority Claims (1)
Number |
Date |
Country |
Kind |
11-226295 |
Aug 1999 |
JP |
|
US Referenced Citations (7)
Foreign Referenced Citations (1)
Number |
Date |
Country |
9833361 |
Jul 1998 |
WO |