WO2011036961A1 - X線ct装置 - Google Patents
X線ct装置 Download PDFInfo
- Publication number
- WO2011036961A1 WO2011036961A1 PCT/JP2010/063714 JP2010063714W WO2011036961A1 WO 2011036961 A1 WO2011036961 A1 WO 2011036961A1 JP 2010063714 W JP2010063714 W JP 2010063714W WO 2011036961 A1 WO2011036961 A1 WO 2011036961A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- ray
- detector
- subject
- incidence rate
- focal position
- Prior art date
Links
- 238000012937 correction Methods 0.000 claims abstract description 62
- 230000035945 sensitivity Effects 0.000 claims abstract description 60
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 20
- 238000007781 pre-processing Methods 0.000 claims description 16
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 description 45
- 238000000034 method Methods 0.000 description 24
- 238000010586 diagram Methods 0.000 description 17
- 238000013480 data collection Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 239000002184 metal Substances 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4021—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4411—Constructional features of apparatus for radiation diagnosis the apparatus being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/587—Alignment of source unit to detector unit
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
Definitions
- the present invention relates to an X-ray CT apparatus that corrects a sensitivity change caused by focal movement.
- X-ray CT (Computed Tomography) is a method for obtaining a tomographic image of a subject by calculating an X-ray absorption coefficient at each point on the tomographic plane by computer reconstruction from multi-directional X-ray projection data of the subject. Characteristic X-rays and continuous X-rays generated when an electron beam is incident on the target metal are used as the X-rays applied to the subject.
- the device target metal, etc.
- the focal point which is the X-ray generation position, moves in the body axis direction.
- the focal point moves, the X-ray irradiation range also moves.
- a problem associated with the movement of the focal position is an artifact that occurs in the image of the edge slice when the detector is displaced in the body axis direction separately from the focal movement.
- the focal position is fixed, the change in the X-ray incidence rate due to the detector position shift can be corrected by the sensitivity correction data.
- the X-ray incidence rate is the ratio of the X-ray dose actually incident on the detector to the X-ray dose irradiated to the region where the detector is viewed from the focal point in the absence of the subject. For example, if a collimator is present in the region where the detector is viewed from the focal point, the X-ray incidence rate is less than 100%.
- an X-ray tube is generally used that generates X-rays by irradiating an anode (target metal) with electrons accelerated at a high voltage.
- the X-ray generation efficiency is low with respect to the energy used for accelerating electrons, and most of the energy is changed to heat. Therefore, the X-ray generation position (focal position) is caused by the thermal expansion of the heated anode (target metal). ) May fluctuate.
- the X-ray tube is provided in a rotating gantry, there is a possibility that the position of the X-ray tube is shifted and the focal position is changed by rotating the gantry.
- Another problem associated with the movement of the focal position is an artifact that occurs in all slices of an image when a two-dimensional grid that removes scattered rays from the subject is used.
- a plurality of X-ray detectors smaller than the collection pitch in the slice direction are arranged between the grids, and the insensitive area corresponds to the focal position obtained by measurement. There is a way to change.
- JP-A-6-269443 Japanese Patent Laid-Open No. 10-234724 JP 2000-93418 A
- the present invention has been made in view of the above points, and it is an object of the present invention to provide an X-ray CT apparatus that estimates sensitivity correction data in a state where the focal position is moved with high accuracy and simple processing.
- the present invention provides one or more detections having an X-ray tube for irradiating X-rays, a collimator for limiting the X-ray irradiation area, and a plurality of detectors for detecting the X-rays.
- Sensitivity correction data is created based on a module, a grid installed in the detector, which removes scattered rays from the subject when the subject is irradiated with X-rays, a storage device, and the output of the detector
- a pre-processing device for correcting the output of the detector using the sensitivity correction data, and reconstructing the image of the subject based on the corrected output.
- An X-ray CT apparatus including an image display device configured to display a configured image, wherein the pre-processing device changes a focal position as a position of the X-ray tube in a state where the subject is not present, and a focal position
- An X-ray CT apparatus comprising means for obtaining sensitivity correction data for each focal position using an X-ray incidence rate at the time of subject imaging.
- the present invention it is possible to estimate the sensitivity correction data in a state where the focal position is moved with high accuracy and simple processing.
- FIG. 1 is a schematic configuration diagram of an X-ray CT apparatus according to a first embodiment of the present invention. It is a figure which shows the module comprised from the detector and 1-dimensional grid which concern on 1st Embodiment. It is a figure which shows an example of the circular arc arrangement
- FIG. 1 is a schematic configuration diagram of an X-ray CT apparatus 100.
- the X-ray CT apparatus 100 includes a detector 1, an X-ray tube 2, a collimator 3, a one-dimensional grid 4, a data acquisition device 6, a preprocessing device 7, a storage device 8, an image reconstruction processing device 9, and an image display device 10. It is the composition which includes.
- the X-ray tube 2 emits X-rays.
- the X-ray tube 2 generates X-rays by irradiating the anode (target metal) with electrons accelerated by a high voltage.
- the X-ray tube 2 is also referred to as a focal point, and the position is defined as a focal position.
- the collimator 3 limits the irradiated X-rays to a predetermined irradiation range (for example, a beam shape). X-rays whose irradiation range is narrowed by the collimator 3 pass through the subject 5 and a bed (not shown).
- a predetermined irradiation range for example, a beam shape
- the one-dimensional grid 4 removes scattered rays in the x direction generated from the subject 5.
- the detector 1 detects X-rays that have passed through the subject 5.
- the data collection device 6 converts the output from the detector 1 from an analog signal to a digital signal.
- the preprocessing device 7 estimates sensitivity correction data corresponding to the focal position.
- the storage device 8 stores the relationship between the X-ray incidence rate in the reference detector that detects X-rays that do not pass through the subject 5 and the X-ray incidence rate in other detectors.
- the image reconstruction processing device 9 performs image reconstruction calculation after performing various correction processes such as sensitivity correction and CT value adjustment on the collected data.
- the image display device 10 displays a reconstructed image.
- the subject 5 is imaged as follows. X-rays generated from the X-ray tube 2 are detected by the detector 1 after passing through the subject 5. From the incident X-ray dose at this time and the X-ray dose detected by the detector 1, the attenuation rate D of the X-ray subject 5 is obtained.
- the attenuation rate D of X-rays passing through the subject 5 is expressed as follows by an exponential function.
- ⁇ is an attenuation coefficient of the subject 5 and the exponential line integration is performed on a straight line connecting the position of the X-ray tube 2 (X-ray generation position) and the detector 1. Therefore, by taking the logarithm of the attenuation rate D, the profile of the integral value of the attenuation coefficient ⁇ in a certain direction can be obtained. Further, by measuring the X-ray while rotating the X-ray tube 2 and the detector 1 around the subject 5, a profile of the integral value of the attenuation coefficient ⁇ in all directions is obtained. By reconstructing an image from these profiles, the distribution of the attenuation coefficient ⁇ in the subject 5 can be acquired as an image.
- FIG. 2 is a diagram showing a module composed of the detector 1 and the one-dimensional grid 4. As shown in FIG. 2, one detector module 11 has N z detectors 1 arranged in the z direction and N x detectors 1 in the x direction.
- FIG. 3 is a diagram showing an example of module arrangement in the X-ray CT apparatus 100. As shown in FIG. 3, the detector module 11 of the X-ray CT apparatus 100 is assumed in a row N m pieces in a circular arc shape in the x direction.
- the X-ray dose measured by the detector 1 is expressed by the following equation.
- I is the intensity of the irradiated X-ray
- S is the intrinsic sensitivity of the detector 1
- R is the X-ray incidence rate.
- m is the module number of the detector module 11 (1 ⁇ m ⁇ N m
- z is the detector number of the detector 1 in the z direction among the detector modules 11 whose module number is m (1 ⁇ z ⁇ N z).
- X is the detector number (1 ⁇ x ⁇ N x ) of the detector 1 in the x direction among the detector modules 11 whose module number is m.
- F represents the focal position in the z direction.
- FIG. 4 is a diagram showing an X-ray irradiation range in N s slice acquisition.
- the focal point 12 in FIG. 4 is the X-ray tube 2 in FIG.
- image reconstruction is performed using the output from the detector 1 in which the detector numbers in the z direction are (N z / 2 ⁇ N s / 2 + 1) to (N z / 2 + N s / 2). . If the position of the focal point 12 moves in the z direction, the collimator 3 covers the irradiation area. Therefore, in the detector 1 of the detector number (N z / 2 + N s / 2) in the z direction that becomes an end slice.
- the X-ray incidence rate R becomes smaller than 100%.
- FIG. 5 is a diagram for explaining that the X-ray incidence rate changes as the irradiation range moves.
- the z-direction shift of the detector 1 in the end slice does not occur. Since the detector 1 of the module number m changes in the z-ray incidence rate because the detector 1 in the partial slice is displaced in the z direction, an artifact is generated in the reconstructed image as a result. Further, as shown in FIG.
- the X-ray incidence rate R is a function that does not depend on the position in the x direction in the detector module 11, and the expression (3) is as follows.
- the reference detector is a detector for directly detecting irradiated X-rays. Therefore, the detector 1 at the end in the x direction that is less likely to be covered by the subject 5 is used as a reference detector.
- x 1 to 4
- the opening width of the collimator 3 is set to full open, and the intrinsic sensitivity correction data X air of the detector 1 is collected without the subject 5.
- the intrinsic sensitivity correction data X air is expressed by the following equation.
- the focal position f of the output value X is omitted.
- the output value X f of the detector 1 is expressed by the following equation.
- the focus 12 itself which is the position of the X-ray tube 2
- the focus 12 is moved by moving the collimator 3 in the z direction to change the irradiation area. It is possible to acquire the same collected data.
- the X-ray incidence rate R can be obtained from the following equation from the equations (4) and (5).
- the X-ray incidence rate R may be obtained by averaging in the x direction in the detector module 11.
- f a polynomial is approximated as follows. Note that f represents the focal position in the z direction.
- Coefficients a n in the formula (7) is a coefficient indicating the relationship between the X-ray incidence rate of X-ray incidence rate of the detector module 11 is a reference signal reference detector is present with each detector module 11. That is, if looking for the coefficients a n advance the least squares method or the like, even during imaging an object 5, the X-ray incidence rate R of each detector module 11 from the X-ray incidence rate R in the reference detector Formula ( 7) can be estimated.
- the X-ray incidence rate R in the reference detector can be obtained by dividing the output value X of the detector 1 by the intrinsic sensitivity S and the X-ray intensity I.
- FIG. 6 is a diagram showing an irradiation range in the 2-slice acquisition mode.
- the aperture width of the collimator 3 is set to the 2-slice acquisition mode
- the X-ray incidence rate R of the reference detector at a certain focal position f is obtained by the equation (6)
- the integral value R sum in the z direction is obtained.
- This value R sum is a constant that does not depend on the focal position f, and is stored in the storage device 8.
- the value I sum integrated in the z direction after correcting the intrinsic sensitivity of the output of the reference detector at this time is (I patient / I) ⁇ R sum. .
- the X-ray incidence rate R in the reference detector is obtained by dividing the output of the reference detector after the intrinsic sensitivity correction by the value obtained by dividing I sum by R sum (I patient / I).
- Sensitivity correction data X ′ air for each is obtained as follows.
- the image reconstruction processing device 9 performs various corrections such as sensitivity correction and CT value adjustment using the sensitivity correction data X ′ air created by the preprocessing device 7 on the collected data detected by the detector 1 at each focal position. After that, the image is reconstructed, and the image display device 10 displays the image.
- the effect of distortion due to the position shift or rotation of the collimator 3 is also included in the change in the X-ray incidence rate.
- Correction data can be estimated.
- the measurement of X-ray incidence rate of the formula (6), i.e., calculation of the coefficients a n may be carried out at time of installation of the apparatus, measured every scan speed and tilt angle parameter. Further, when the detector module 11 and the X-ray tube 2 are exchanged, the X-ray incidence rate of the equation (6) is also measured.
- the collimator aperture width is set to full open and measurement is performed only once. Can be greatly reduced.
- the present invention can be combined with a method of moving the collimator 3 in accordance with the focal position as disclosed in Japanese Patent Laid-Open No. 4-227238.
- the value of the X-ray incidence rate for each detector module 11 is approximated by the value of the reference detector in the same slice as shown in Equation (7), but may be approximated as follows: Good.
- Equation (9) it is also possible to approximate with the sum of a plurality of reference detectors in the z direction.
- the X-ray incidence rate of the reference detector at the time of imaging the subject 5 is also obtained by the sum of the plurality of detectors.
- Statistical accuracy can be improved by using the X-ray incidence rate added in the z direction.
- the above-described embodiment can be applied even in the maximum slice acquisition mode, and the shadow of the collimator 3 can be made with respect to the reference detector by shifting the detector module 11.
- the X-ray incidence rate of the reference detector and the X-ray incidence rate of other detector modules 11 can be related.
- FIG. 8 is a schematic configuration diagram of an X-ray CT apparatus 200 according to the second embodiment.
- the X-ray CT apparatus 200 includes a detector 1, an X-ray tube 2, a collimator 3, a one-dimensional grid 4, a data acquisition device 6, a preprocessing device 7, a storage device 8, an image reconstruction processing device 9, an image display device 10, The focus position detector 13 and the focus position measuring device 14 are included.
- the detector 1, the X-ray tube 2, the collimator 3, the one-dimensional grid 4, the data acquisition device 6, the image reconstruction processing device 9, and the image display device 10 have the same configuration as that of the X-ray CT apparatus 100 shown in FIG. It has a function.
- the focal position detector 13 detects the position of the X-ray tube 2 (focal point 12), and the focal position measuring device 14 obtains focal position coordinates.
- the storage device 8 stores the relationship between the focal position measured by the focal position measurement device 14 and the X-ray incidence rate in the detector 1.
- the detector 1 of the X-ray CT apparatus 200 is mounted in units of modules. Similarly to the case of the X-ray CT apparatus 100, the detector 1 is arranged in one detector module 11 as shown in FIG. The detector module 11 is arranged in an arc shape as shown in FIG. In X-ray CT apparatus 200, a case where deviation for each detector module 11 for N s slice collection mode has occurred.
- the X-ray incidence rate for each detector module 11 is expressed by the X-ray incidence rate in the reference detector as shown in Expression (7), but in the second embodiment, detection is performed.
- the X-ray incidence rate R for each instrument module 11 is represented by a focal position f. That is, the focal position measured by the focal position detector 13 and the focal position measuring device 14 becomes the reference signal.
- the opening width of the collimator 3 is set to be fully open, and the intrinsic sensitivity of the detector 1 is collected without the subject 5 (formula (4)). Then, set the collimator 3 in the opening width of the N s slice collection mode, collected in the absence of the subject 5 while changing the focal position f (Equation (5)). At this time, the focal position f is simultaneously measured by the focal position detector 13 and the focal position measuring device 14. Thus, equation (6) is obtained from equations (4) and (5).
- the X-ray incidence rate R in each detector module 11 is approximated by a polynomial as follows according to the focal position f.
- the coefficient b n in the expression (10) is a coefficient indicating the relationship between the focal position f as a reference signal and the X-ray incidence rate R for each detector module 11. That is, if the coefficient b n is obtained in advance by the least square method or the like and stored in the storage device 8, the X-ray incidence rate of each detector module 11 from each focal position f even during imaging of the subject 5. R can be estimated from equation (10).
- the coefficient b n is obtained for each detector module 11 in the end slice.
- Sensitivity correction data corresponding to the focal position is created by multiplying the calculated X-ray incidence rate R by equation (4).
- the collected data is reconstructed after various corrections such as sensitivity correction and CT value adjustment by sensitivity correction data corresponding to the focal position f in the image reconstruction processing device 9, and an image is displayed on the image display device 10.
- the accuracy of the resulting sensitivity correction data X ′ air is also increased.
- the polynomial in equation (10) is actually a second order approximation, and the number of coefficients b n that must be held in the storage device 8 can be reduced.
- the coefficient b n may be obtained when the apparatus is installed.
- sensitivity correction data for each view corresponding to the focal position is obtained.
- the influence of the deviation of the irradiation area due to this control delay can be reduced by performing sensitivity correction using sensitivity correction data for each view after data collection is completed.
- FIG. 9 is a schematic configuration diagram of an X-ray CT apparatus 300 according to the third embodiment.
- the X-ray CT apparatus 300 includes a detector 1, an X-ray tube 2, a collimator 3, a data acquisition device 6, a preprocessing device 7, a storage device 8, an image reconstruction processing device 9, an image display device 10, a two-dimensional grid 15,
- the configuration includes an X-ray intensity detector 16 and an X-ray intensity measurement device 17.
- the detector 1, the X-ray tube 2, the collimator 3, the data collection device 6, the storage device 8, the image reconstruction processing device 9, and the image display device 10 have the same configuration and function as those of the X-ray CT apparatus 100 shown in FIG. Have
- the two-dimensional grid 15 removes scattered rays in the x and z directions generated from the subject 5.
- the X-ray intensity detector 16 detects the irradiated X-ray intensity, and the X-ray intensity measuring device 17 measures the X-ray intensity.
- FIG. 10 is a diagram showing a detector module 11 including a detector and a two-dimensional grid. As shown in FIG. 10, one detector detector module 11 has N z detectors 1 arranged in the z direction and N x detectors 1 in the x direction. A two-dimensional grid 15 for removing scattered rays in the x and z directions generated from the subject 5 is mounted on the front surface of the detector 1. Further, the detector module 11 in the X-ray CT apparatus 300, as shown in FIG. 3, and shall be N m pieces arranged in a circular arc shape in the x-direction.
- the X-ray CT apparatus 300 a case where deviation occurs in each detector module 11 in the imaging of the N s slice collection mode.
- a method of obtaining will be described.
- the aperture width of the collimator 3 is set to fully open, and the intrinsic sensitivity of the detector 1 is collected without the subject 5.
- the output X air of the detector is a function depending on the focal position f as in the following equation.
- f air is the focal position when the intrinsic sensitivity S of the detector is collected, and the X-ray incidence rate at that position is R (z, m, f air ). Then, set the collimator 3 in the opening width of the N s slice acquisition mode, when collected in the absence of the subject 5 while changing the focal position, the output X f of the detector 1 at the focal position f is as follows.
- R (z, m, f) is the X-ray incidence rate at the focal position f.
- R ′ is expressed as the following equation.
- Coefficients c n in the formula (14) is a coefficient indicating the relationship between the X-ray incidence rate of X-ray incidence rate of the detector module 11 is a reference signal reference detector is present with each detector module 11. That is, if seeking the coefficients c n advance the least squares method or the like, even during imaging an object 5, the 'X-ray incidence rate of each detector module 11 from R' X-ray incidence rate R in the reference detector It can be estimated from equation (14). Coefficients c n are determined for each slice in each detector module 11, stored in the storage device 8. In other slice acquisition modes, the coefficient c n is obtained for each slice in each detector module 11.
- the X-ray incidence rate R ′ in the reference detector is obtained by dividing the output value X of the detector by the intrinsic sensitivity S and the X-ray intensity I.
- the X-ray intensity I is measured by the X-ray intensity detector 16 and the X-ray intensity measuring device 17.
- Sensitivity correction data X ′ air for each focal position is obtained as follows.
- the image reconstruction processing device 9 performs various corrections such as sensitivity correction and CT value adjustment using the sensitivity correction data X ′ air created by the preprocessing device 7 on the collected data detected by the detector 1 at each focal position. After that, the image is reconstructed, and the image display device 10 displays the image.
- FIG. 11 is a diagram showing a detector module 11 composed of a two-dimensional grid from which the detector and a part of the grid are removed.
- FIG. 12 is a schematic configuration diagram of an X-ray CT apparatus 400 according to the fourth embodiment.
- the X-ray CT apparatus 400 includes a detector 1, an X-ray tube 2, a collimator 3, a data acquisition device 6, a preprocessing device 7, a storage device 8, an image reconstruction processing device 9, an image display device 10, a two-dimensional grid 15,
- the configuration includes an X-ray intensity / focus position detector 18 and an X-ray intensity / focus position measurement device 19.
- the detector 1, the X-ray tube 2, the collimator 3, the data acquisition device 6, the image reconstruction processing device 9, and the image display device 10 have the same configuration and function as those of the X-ray CT apparatus 100 shown in FIG.
- the two-dimensional grid 15 removes scattered rays in the x direction and the z direction that are generated when the subject 5 is irradiated with X-rays, and has a structure shown in FIG.
- the X-ray intensity / focus position detector 18 detects the position (focus position) and X-ray intensity of the X-ray tube 2, and the X-ray intensity / focus position measurement device 19 detects the position (focus position) of the X-ray tube 2.
- the X-ray intensity irradiated from the X-ray tube 2 is measured.
- the detector 1 of the X-ray CT apparatus 400 is mounted in units of modules. Similarly to the detector module 11 shown in FIG. 10, one detector module 11 of the X-ray CT apparatus 400 is provided with N z detectors 1 in the z direction and N x detectors 1 in the x direction. A two-dimensional grid 15 for removing scattered rays in the x and z directions generated from the subject 5 is mounted on the front surface of the detector 1. Further, the detector module 11 in the X-ray CT apparatus 400, as shown in FIG. 3, and shall be N m pieces arranged in a circular arc shape in the x-direction.
- the X-ray CT apparatus 400 similarly to the X-ray CT apparatus 300, since the shadow is caused by the two-dimensional grid 15 due to the focus movement, the X-ray incidence rate in all slices is obtained.
- the X-ray CT apparatus 400 consider the case where deviation for N s slice collection mode for each detector module 11 has occurred, determining the X-ray incidence rate for each slice from the measured focal position in each of the detector modules 11 A method will be described.
- the opening width of the collimator 3 is set to be fully open, and the intrinsic sensitivity of the detector 1 is collected without the subject 5 (formula (11)). Then, set the collimator 3 in the opening width of the N s slice collection mode, collected in the absence of the subject 5 while changing the focal position f (Equation (12)). At this time, the focal position f is simultaneously measured by the X-ray intensity / focus position detector 18 and the X-ray intensity / focus position measuring device 19. Thus, equation (13) is obtained from equations (11) and (12). Next, the X-ray incidence rate R ′ in each detector module 11 is approximated by a polynomial as follows according to the focal position f.
- Coefficients d n in formula (16) is a coefficient indicating the relationship between the X-ray incidence rate R for each and a reference signal the focus position f detector module 11. That is, to previously obtain the coefficients d n advance the least squares method or the like, be stored in the storage device 8, even during imaging an object 5, X-ray incidence rate of each detector module 11 from each focus position f R ′ can be estimated from equation (16).
- Coefficients d n is determined for each slice in each detector module 11, stored in the storage device 8. Also, in other slice acquisition modes, the coefficient d n is obtained for each slice in each detector module 11. By multiplying the value of the X-ray incidence rate R ′ obtained by the equation (16) by the equation (11), sensitivity correction data X ′ air for each focal position is obtained as follows.
- the X-ray intensity I is measured by the X-ray intensity / focus position detector 18 and the X-ray intensity / focus position measuring device 19.
- the collected data detected by the detector 1 at each focal position f is corrected as sensitivity correction data X ′ air by the preprocessing device 7 according to the above-described processing, and is output to the image reconstruction processing device 9.
- the image reconstruction processing device 9 performs various corrections such as CT value adjustment on the sensitivity correction data X ′ air , then reconstructs the image, and the image display device 10 displays the image.
- the X-ray incidence rate corresponding to the focal position is obtained by a continuous function such as Expression (16)
- the accuracy of the sensitivity correction data X ′ air obtained as a result is obtained. Also gets higher.
- the polynomial in Equation (16) is sufficient in practice, approximations about the secondary, only a number of coefficients d n must be held in the storage device 8 is also small. The coefficient d n is measured, such as during installation of the device may be put and stored.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
焦点の移動による入射X線量の変化を高精度かつ簡便に補正するX線CT装置を提供すること。 前処理装置(7)は、被写体(5)が無い状態において、X線管(2)の位置となる焦点位置を変え、焦点位置毎における各検出モジュールのX線入射率を取得し、その各検出モジュールのX線入射率を被写体(5)がない状態で求めたリファレンス検出器が存在するモジュールのX線入射率で近似し、近似多項式の係数を記憶装置(8)に記憶する。被写体(5)の撮像時において、リファレンス検出器が存在するモジュールのX線入射率と記憶した係数を用いて、各検出モジュールの被写体撮像時のX線入射率を算出し、算出された被写体撮像時のX線入射率を用いて、焦点位置毎の感度補正データを求める。
Description
本発明は、焦点移動により生じる感度変化を補正するX線CT装置に関する。
X線CT(Computed Tomography)は、被写体の多方向のX線投影データからコンピュータによる再構成により断層面の各点のX線吸収係数を算出し、被写体の断層像を得る手法である。被写体に照射するX線には、電子線をターゲット金属に入射したときに発生する特性X線と連続X線が用いられる。このとき、電子線のエネルギ付与により機器(ターゲット金属など)は熱伸を生じ、X線の発生位置である焦点は体軸方向へ移動してしまう。焦点の移動に伴ってX線の照射範囲も移動してしまう。
この焦点位置の移動に伴う問題として、焦点移動とは別に、検出器が体軸方向にずれを生じている場合に、端部スライスの画像において発生するアーチファクトが挙げられる。もしも焦点位置が固定であるならば、検出器位置ずれによるX線入射率の変化は感度補正データによって補正することが可能である。ここで、X線入射率とは、被写体が無い状態において、焦点から検出器を見込む領域に照射されるX線量に対する、実際に検出器に入射するX線量の割合である。例えば、焦点から検出器を見込む領域にコリメータが存在すれば、X線入射率は100%よりも小さくなる。
ここで、X線源には、通常、高電圧で加速された電子を陽極(ターゲット金属)に照射し、X線を発生させるX線管が用いられる。しかし、電子の加速に使われたエネルギに対してX線の発生効率は低く、ほとんどのエネルギは熱に変わるため、加熱された陽極(ターゲット金属)の熱膨張によりX線の発生位置(焦点位置)が変動する現象が生じるおそれがある。また、X線管は回動するガントリに備えられるため、ガントリが回動することによりX線管の位置がずれ、焦点位置が変化するおそれがある。そのため、感度補正データ収集時における焦点位置と患者(被検体、被写体)データ収集時における焦点位置が異なっていると、検出器のX線入射率が変化するため、感度補正は不完全となり画像にアーチファクトが生じる。コリメータの開口幅を広く設定することでアーチファクトを低減できるが、無効被ばくが生じる。
この対策として、特許文献1或いは特許文献2に記載されるように、焦点移動範囲をいくつかに分割し、それぞれの小範囲毎に予めX線検出器感度補正データを収集、保管しておき、被検体のデータ収集時の焦点位置に対応するX線検出器感度補正データを用いて,X線検出器感度を補正する方法がある。
ここで、X線源には、通常、高電圧で加速された電子を陽極(ターゲット金属)に照射し、X線を発生させるX線管が用いられる。しかし、電子の加速に使われたエネルギに対してX線の発生効率は低く、ほとんどのエネルギは熱に変わるため、加熱された陽極(ターゲット金属)の熱膨張によりX線の発生位置(焦点位置)が変動する現象が生じるおそれがある。また、X線管は回動するガントリに備えられるため、ガントリが回動することによりX線管の位置がずれ、焦点位置が変化するおそれがある。そのため、感度補正データ収集時における焦点位置と患者(被検体、被写体)データ収集時における焦点位置が異なっていると、検出器のX線入射率が変化するため、感度補正は不完全となり画像にアーチファクトが生じる。コリメータの開口幅を広く設定することでアーチファクトを低減できるが、無効被ばくが生じる。
この対策として、特許文献1或いは特許文献2に記載されるように、焦点移動範囲をいくつかに分割し、それぞれの小範囲毎に予めX線検出器感度補正データを収集、保管しておき、被検体のデータ収集時の焦点位置に対応するX線検出器感度補正データを用いて,X線検出器感度を補正する方法がある。
焦点位置の移動に伴うもう一つの問題として、被写体の散乱線を除去する2次元グリッドを用いた場合に、画像の全てのスライスにおいて発生するアーチファクトが挙げられる。
前記の問題の対策として、特許文献3に記載されるように、スライス方向の収集ピッチよりも小さい複数のX線検出器をグリッド間に配置し、不感領域を測定によって求めた焦点位置に対応して変化させる方法がある。
前記の問題の対策として、特許文献3に記載されるように、スライス方向の収集ピッチよりも小さい複数のX線検出器をグリッド間に配置し、不感領域を測定によって求めた焦点位置に対応して変化させる方法がある。
しかしながら、特許文献1及び特許文献2に記載された方法では、X線源である焦点移動が連続的であるのに対して感度補正データは離散的であるため、補正の精度が低下してしまうという問題がある。
また、焦点位置に対応した感度補正データが検出器毎に必要であり、X線CT装置が有する記憶装置内に記憶しなければならないデータ数は膨大となる。
また、特許文献3に記載された方法では、実際の収集ピッチよりも小さい複数の検出器を用いる必要があり、それに付随して処理しなければならない信号の数は増大し、全体の装置構成が複雑となる。
また、焦点位置に対応した感度補正データが検出器毎に必要であり、X線CT装置が有する記憶装置内に記憶しなければならないデータ数は膨大となる。
また、特許文献3に記載された方法では、実際の収集ピッチよりも小さい複数の検出器を用いる必要があり、それに付随して処理しなければならない信号の数は増大し、全体の装置構成が複雑となる。
本発明は、以上の点に鑑みてなされたものであり、焦点位置が移動した状態における感度補正データを高い精度で且つ簡便な処理により推定するX線CT装置を提供することにある。
前述した目的を達成するために本発明は、X線を照射するX線管と、前記X線の照射領域を制限するコリメータと、前記X線を検出する複数の検出器を有する1以上の検出モジュールと、前記検出器に設置され、前記X線が被写体に照射された場合に前記被写体からの散乱線を除去するグリッドと、記憶装置と、前記検出器の出力を元に感度補正データを作成する前処理装置と、前記感度補正データを用いて前記検出器の出力を補正して、前記補正した出力に基づき前記被写体の画像を再構成する画像再構成装置と、前記画像再構成装置で再構成した画像を表示する画像表示装置と、を備えたX線CT装置であって、前記前処理装置は、前記被写体が無い状態において、前記X線管の位置となる焦点位置を変え、焦点位置毎における各検出モジュールのX線入射率を取得する手段と、前記被写体が無い状態において、参照信号と前記取得した各検出モジュールのX線入射率との関係を求める手段と、前記求めた関係を前記記憶装置に記憶する手段と、前記被写体の撮像時において、被写体撮像時の参照信号と前記記憶された関係を用いて、各検出モジュールの被写体撮像時のX線入射率を算出する手段と、前記算出された被写体撮像時のX線入射率を用いて、焦点位置毎の感度補正データを求める手段と、を有することを特徴とするX線CT装置である。
本発明によれば、焦点位置が移動した状態における感度補正データを高い精度で且つ簡便な処理により推定することが可能である。
≪第1の実施形態≫
以下に、添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。なお、以下の説明および添付図面において、略同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略することとする。
以下に、添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。なお、以下の説明および添付図面において、略同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略することとする。
図1は、X線CT装置100の概略構成図である。
X線CT装置100は、検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10を含む構成である。
X線管2はX線を照射する。X線管2は、高電圧で加速された電子を陽極(ターゲット金属)に照射することにより、X線を発生させる。以下の説明ではX線管2を焦点とも呼び、その位置を焦点位置とする。コリメータ3は、照射されたX線を所定の照射範囲(例えば、ビーム形状)に制限する。コリメータ3によって照射範囲を絞られたX線は、被射体5や寝台(図示せず)を透過する。
X線CT装置100は、検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10を含む構成である。
X線管2はX線を照射する。X線管2は、高電圧で加速された電子を陽極(ターゲット金属)に照射することにより、X線を発生させる。以下の説明ではX線管2を焦点とも呼び、その位置を焦点位置とする。コリメータ3は、照射されたX線を所定の照射範囲(例えば、ビーム形状)に制限する。コリメータ3によって照射範囲を絞られたX線は、被射体5や寝台(図示せず)を透過する。
1次元グリッド4は、被写体5から発生するx方向の散乱線を除去する。検出器1は、被写体5を透過したX線を検出する。データ収集装置6は、検出器1からの出力をアナログ信号からディジタル信号に変換する。前処理装置7は、焦点位置に対応した感度補正データを推定する。記憶装置8は、被写体5を通過しないX線を検出するリファレンス検出器におけるX線入射率と他の検出器におけるX線入射率との関係を記憶する。画像再構成処理装置9は、収集データに対して感度補正やCT値調整などの様々な補正処理を実行後、画像再構成演算を実行する。画像表示装置10は、再構成画像を表示する。
X線CT装置100では次のように被写体5の撮像が行われる。
X線管2から発生するX線は被写体5を通過した後、検出器1によって検出される。このときの入射X線量と検出器1で検出されるX線量から、X線の被写体5による減衰率Dを求める。被写体5を通過するX線の減衰率Dは指数関数により以下のように表される。
X線管2から発生するX線は被写体5を通過した後、検出器1によって検出される。このときの入射X線量と検出器1で検出されるX線量から、X線の被写体5による減衰率Dを求める。被写体5を通過するX線の減衰率Dは指数関数により以下のように表される。
ここで、μは被写体5の減弱係数であり、指数関数の線積分はX線管2の位置(X線の発生位置)と検出器1を結ぶ直線上で行われる。したがって、この減衰率Dの対数をとることによって、ある方向における減弱係数μの積分値のプロファイルを求めることができる。さらに、X線管2と検出器1を被写体5の周りを回転させながらX線を測定することで、全方向における減弱係数μの積分値のプロファイルを求める。これらのプロファイルから画像再構成することによって、被写体5内の減弱係数μの分布を画像として取得することができる。
次に、前処理装置7で実行する、焦点位置に対応した感度補正データの推定方法について説明する。
本実施形態におけるX線CT装置100の検出器1は、モジュール単位で実装される。図2は、検出器1と1次元グリッド4から構成されるモジュールを示す図である。図2に示すように、1つの検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置されている。
図3は、X線CT装置100におけるモジュールの配置の一例を示す図である。図3に示すように、X線CT装置100の検出器モジュール11は、x方向に円弧状にNm個並んでいるものとする。
ここで、被写体5が無い状態においてX線を照射した場合、検出器1で測定されるX線量は次式で表される。
本実施形態におけるX線CT装置100の検出器1は、モジュール単位で実装される。図2は、検出器1と1次元グリッド4から構成されるモジュールを示す図である。図2に示すように、1つの検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置されている。
図3は、X線CT装置100におけるモジュールの配置の一例を示す図である。図3に示すように、X線CT装置100の検出器モジュール11は、x方向に円弧状にNm個並んでいるものとする。
ここで、被写体5が無い状態においてX線を照射した場合、検出器1で測定されるX線量は次式で表される。
この式(2)で、Iは照射するX線の強度、Sは検出器1の固有感度、RはX線入射率である。mは検出器モジュール11のモジュール番号(1≦m≦Nm)、zはモジュール番号がmである検出器モジュール11の内のz方向の検出器1の検出器番号(1≦z≦Nz)、xはモジュール番号がmである検出器モジュール11の内のx方向の検出器1の検出器番号(1≦x≦Nx)である。また、fはz方向の焦点位置を表す。
ここで、X線CT装置100による撮像をx方向に並んだ検出器1をNs列用いて画像を収集するNsスライス収集モードで撮像するとする。図4は、Nsスライス収集におけるX線照射範囲を示す図である。図4における焦点12は図1におけるX線管2である。Nsスライス収集モードでは、z方向の検出器番号が(Nz/2-Ns/2+1)~(Nz/2+Ns/2)における検出器1での出力を用いて画像再構成を行う。
焦点12の位置がz方向に移動してしまうと、コリメータ3が照射領域を覆ってしまうため、端部スライスとなるz方向の検出器番号(Nz/2+Ns/2)の検出器1におけるX線入射率Rは100%よりも小さくなってしまう。
焦点12の位置がz方向に移動してしまうと、コリメータ3が照射領域を覆ってしまうため、端部スライスとなるz方向の検出器番号(Nz/2+Ns/2)の検出器1におけるX線入射率Rは100%よりも小さくなってしまう。
さらに、検出器1がz方向に互いに位置ずれを生じている場合、感度補正データ収集時の焦点位置と被写体5の撮像時の焦点位置とが異なると、画像にアーチファクトが生じてしまう。
図5は、照射範囲が移動することでX線入射率が変化することを説明するための図である。例えば、図5(a)に示すモジュール番号(m+1)の検出器モジュール11では、端部スライスにおける検出器1のz方向のずれは生じていないが、モジュール番号mの検出器モジュール11では、端部スライスにおける検出器1のz方向のずれが生じるため、モジュール番号mの検出器モジュール11ではX線入射率が変わってしまい、結果、再構成画像にアーチファクトが生ずることになる。また、図5(b)に示すように、焦点12の位置がz方向に移動し(図4参照)、かつ、モジュール番号mの検出器モジュール11がz方向に位置ずれを生じている場合にもX線入射率が変わってしまい、再構成画像にアーチファクトが生ずることになる。
このように、実際の検出器1のずれはモジュール単位で起こることが多いため、本実施形態においても検出器モジュール11毎にずれが発生している場合を考える。
従ってX線入射率Rは検出器モジュール11内のx方向の位置に依存しない関数となり、式(3)は次式のようになる。
図5は、照射範囲が移動することでX線入射率が変化することを説明するための図である。例えば、図5(a)に示すモジュール番号(m+1)の検出器モジュール11では、端部スライスにおける検出器1のz方向のずれは生じていないが、モジュール番号mの検出器モジュール11では、端部スライスにおける検出器1のz方向のずれが生じるため、モジュール番号mの検出器モジュール11ではX線入射率が変わってしまい、結果、再構成画像にアーチファクトが生ずることになる。また、図5(b)に示すように、焦点12の位置がz方向に移動し(図4参照)、かつ、モジュール番号mの検出器モジュール11がz方向に位置ずれを生じている場合にもX線入射率が変わってしまい、再構成画像にアーチファクトが生ずることになる。
このように、実際の検出器1のずれはモジュール単位で起こることが多いため、本実施形態においても検出器モジュール11毎にずれが発生している場合を考える。
従ってX線入射率Rは検出器モジュール11内のx方向の位置に依存しない関数となり、式(3)は次式のようになる。
即ち、端部スライスにおけるX線入射率Rが検出器モジュール11毎に求まれば、アーチファクトを生じないデータを取得することができる。
次に、リファレンス検出器におけるX線入射率Rから他の検出器モジュール11におけるX線入射率Rを求める方法について説明する。リファレンス検出器は、照射されたX線を直接検出するための検出器である。そのため、被写体5に覆われる可能性が少ないx方向の端部の検出器1がリファレンス検出器として用いられる。
次に、リファレンス検出器におけるX線入射率Rから他の検出器モジュール11におけるX線入射率Rを求める方法について説明する。リファレンス検出器は、照射されたX線を直接検出するための検出器である。そのため、被写体5に覆われる可能性が少ないx方向の端部の検出器1がリファレンス検出器として用いられる。
本実施例では、モジュール番号m=1における検出器モジュール11のx方向の検出器番号1~4の検出器をリファレンス検出器に設定する。例えば、Nsスライス収集モードにおける端部スライス(z=Nz/2+Ns/2)におけるリファレンス検出器の出力として、m=1、z=Nz/2+Ns/2、x=1~4の4つの検出器1(図2に示す検出器モジュール11がモジュール番号m=1であるとすると、8スライス(Nz=8)の場合、枠A内の検出器)の平均値を用いる。
はじめに、コリメータ3の開口幅を全開に設定し、被写体5の無い状態で検出器1の固有感度補正データXairを収集する。このとき、X線照射率Rは100%となるため、固有感度補正データXairは次式のようになる。
ただし、検出器1の出力は焦点位置に依存しないとして、出力値Xの焦点位置fを省略した。
次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置を変えながら被写体5の無い状態で収集すると、検出器1の出力値Xfは次式のようになる。
次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置を変えながら被写体5の無い状態で収集すると、検出器1の出力値Xfは次式のようになる。
出力値Xfの収集において、X線管2の位置である焦点12そのものを移動させながら収集させてもよいが、コリメータ3をz方向に移動させて照射領域を変えることで、焦点12を移動させた場合と同じ収集データを取得することができる。
X線入射率Rは式(4)と式(5)とから次式のように求められる。
X線入射率Rは式(4)と式(5)とから次式のように求められる。
ここで、統計精度を向上させるために、検出器モジュール11内のx方向における平均によりX線入射率Rを求めてもよい。
次に、検出器モジュール11毎のX線入射率R(z,m,f)をリファレンス検出器が存在する検出器モジュール11(モジュール番号m=1)のX線入射率R(z,1,f)によって以下のように多項式で近似する。なお、fはz方向の焦点位置を表す。
次に、検出器モジュール11毎のX線入射率R(z,m,f)をリファレンス検出器が存在する検出器モジュール11(モジュール番号m=1)のX線入射率R(z,1,f)によって以下のように多項式で近似する。なお、fはz方向の焦点位置を表す。
式(7)における係数anは、参照信号であるリファレンス検出器が存在する検出器モジュール11のX線入射率と検出器モジュール11毎のX線入射率との関係を示す係数である。即ち、この係数anを予め最小自乗法等で求めておけば、被写体5を撮像中においても、リファレンス検出器におけるX線入射率Rから各検出器モジュール11のX線入射率Rを式(7)より推定することが可能である。係数anは、端部スライス(z=Nz/2+Ns/2)において検出器モジュール11毎に求め、記憶装置8に記憶される。また、他のスライス収集モードにおいても、係数anは、端部スライスにおいて検出器モジュール11毎に求める。
被写体5を撮像している場合、リファレンス検出器におけるX線入射率Rは、検出器1の出力値Xを固有感度SおよびX線強度Iで除することで求められる。図4に示すようにNsスライス収集モードであれば、X線強度Iは端部スライス{(z=Nz/2-Ns/2+1),(z=Nz/2+Ns/2)}の内側のスライス{(z=Nz/2-Ns/2+2)から(z=Nz/2+Ns/2-1)}では焦点12の位置ずれや検出器モジュール11の位置ずれの影響を受けない。このため、被写体5を撮像している場合のリファレンス検出器におけるX線強度Iは内側のスライス{(z=Nz/2-Ns/2+2)から(z=Nz/2+Ns/2-1)}の検出器1の出力から求められる。
ただし、2スライス収集モード(Ns=2)の場合は、2つのスライス{(z=Nz/2),(z=Nz/2+1)}が端部スライスとなり、内側のスライスが存在せず内側のスライスの検出器出力が得られないことから、次のようにX線強度Iを求める。
ただし、2スライス収集モード(Ns=2)の場合は、2つのスライス{(z=Nz/2),(z=Nz/2+1)}が端部スライスとなり、内側のスライスが存在せず内側のスライスの検出器出力が得られないことから、次のようにX線強度Iを求める。
図6は、2スライス収集モードにおける照射範囲を示す図である。
はじめに、コリメータ3の開口幅を2スライス収集モードに設定し、ある焦点位置fにおけるリファレンス検出器のX線入射率Rを式(6)により求め、それらのz方向への積分値Rsumを求める。この値Rsumは、焦点位置fによらない定数であり、記憶装置8に記憶される。被写体5の撮像時のX線強度をIpatientとすると、このときのリファレンス検出器の出力を固有感度補正した後、z方向へ積分した値Isumは(Ipatient/I)×Rsumとなる。ここで、Iは固有感度補正データ収集時のX線の強度である。したがって、IsumをRsumで割った値(Ipatient/I)で、固有感度補正後のリファレンス検出器の出力を除することで、リファレンス検出器におけるX線入射率Rを求める。
はじめに、コリメータ3の開口幅を2スライス収集モードに設定し、ある焦点位置fにおけるリファレンス検出器のX線入射率Rを式(6)により求め、それらのz方向への積分値Rsumを求める。この値Rsumは、焦点位置fによらない定数であり、記憶装置8に記憶される。被写体5の撮像時のX線強度をIpatientとすると、このときのリファレンス検出器の出力を固有感度補正した後、z方向へ積分した値Isumは(Ipatient/I)×Rsumとなる。ここで、Iは固有感度補正データ収集時のX線の強度である。したがって、IsumをRsumで割った値(Ipatient/I)で、固有感度補正後のリファレンス検出器の出力を除することで、リファレンス検出器におけるX線入射率Rを求める。
以上のようにして、リファレンス検出器におけるX線入射率Rを求めれば、各検出器モジュール11におけるX線入射率も求めることが可能であり、その値を式(4)に乗じることで焦点位置毎の感度補正データX'airが次式のように得られる。
画像再構成処理装置9は、各焦点位置において検出器1によって検出された収集データを、前処理装置7によって作成された感度補正データX'airによる感度補正やCT値調整などの各種補正を行った後に画像を再構成し、画像表示装置10はその画像を表示する。
本発明は、焦点位置に対応したX線入射率を式(7)のような連続的な関数によって求めているため、結果として得られる式(8)の感度補正データX'airの精度も高くなる。
また、式(7)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数anの数も少なくてすむ。特許文献1及び特許文献2に記載されている方法では、焦点位置毎に各検出器の補正データを記憶する必要があるが、本発明では係数anを記憶装置8に保持しておけばよい。
また、式(7)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数anの数も少なくてすむ。特許文献1及び特許文献2に記載されている方法では、焦点位置毎に各検出器の補正データを記憶する必要があるが、本発明では係数anを記憶装置8に保持しておけばよい。
また、前記のリファレンス検出器を参照信号とする方法では、コリメータ3の位置ずれや回転による歪みの効果もX線入射率の変化に含まれるため、焦点移動だけでなくこれらの効果も補正できる感度補正データの推定が可能である。
また、式(6)のX線入射率の測定、即ち、係数anの算出は装置の据付時に行えばよく、スキャン速度とチルト角度のパラメータ毎に測定する。また、検出器モジュール11やX線管2の交換を行った場合にも式(6)のX線入射率の測定を行う。
また、式(6)のX線入射率の測定、即ち、係数anの算出は装置の据付時に行えばよく、スキャン速度とチルト角度のパラメータ毎に測定する。また、検出器モジュール11やX線管2の交換を行った場合にも式(6)のX線入射率の測定を行う。
従来の感度補正データはスライス収集モード毎に開口幅を変えて測定する必要があったが、本発明ではコリメータ開口幅を全開に設定して1度だけ測定すればよく、日々の感度補正データ収集の頻度を大幅に低減できる。
本発明は、特開平4-227238号公報に示すような焦点位置に対応してコリメータ
3を移動させる方法と組み合わせることも可能である。
本発明は、特開平4-227238号公報に示すような焦点位置に対応してコリメータ
3を移動させる方法と組み合わせることも可能である。
次に、本実施の形態に関する変形例について説明する。
前記の実施形態では、検出器モジュール11毎のX線入射率の値は、式(7)に示すように、同一スライスにおけるリファレンス検出器の値により近似したが、以下のように近似してもよい。
前記の実施形態では、検出器モジュール11毎のX線入射率の値は、式(7)に示すように、同一スライスにおけるリファレンス検出器の値により近似したが、以下のように近似してもよい。
式(9)に示すように、z方向における複数のリファレンス検出器についての和で近似することも可能である。このとき、被写体5の撮像時におけるリファレンス検出器のX線入射率も、複数の検出器の和で求める。z方向に加算したX線入射率を用いることで、統計精度を向上することができる。
また、別の変形例として、被写体5の撮像時におけるリファレンス検出器のX線入射率をビュー方向、即ち、検出器1の回転方向への加算平均によって求める方法がある。焦点の移動が無視できる範囲のビュー数においてX線入射率の平均を求めることで、統計精度を向上することができる。
また、前記の実施形態では、リファレンス検出器が存在する検出器モジュール11のモジュール番号はm=1としてきたが、m=Nmにおける検出器モジュール11の検出器1をリファレンス検出器に用いてもよい。どちらの検出器モジュール11におけるリファレンス検出器を用いるかについては、検出器1の出力値により判断する。つまり、リファレンス検出器に入射したX線が被写体を通過したかを判断するために閾値を設定し、この閾値よりも出力が高い(減衰が小さい)場合のリファレンス検出器の出力を用いる。
モジュール番号m=1、m=Nmの両方の検出器モジュール11におけるリファレンス検出器が有効であるならば、それぞれのリファレンス検出器におけるX線入射率から、他の検出器モジュール11のX線入射率を求め、それらの加算平均を用いて感度補正データを計算してもよい。
モジュール番号m=1、m=Nmの両方の検出器モジュール11におけるリファレンス検出器が有効であるならば、それぞれのリファレンス検出器におけるX線入射率から、他の検出器モジュール11のX線入射率を求め、それらの加算平均を用いて感度補正データを計算してもよい。
また、前記の実施形態を利用して最大スライス収集モードにおける補正処理を行うことも可能である。
図7は、検出器モジュール11を装置に実装する際に、モジュール番号m=1とm=Nmの検出器モジュール11をスライス中心軸よりも意図的にずらして配置したモジュール配列を示す図である。図7に示すように配置することで、最大スライス収集モードにおいても前記実施形態は適用可能となり、検出器モジュール11をずらすことにより、リファレンス検出器に対してコリメータ3の影を作ることができるため、リファレンス検出器のX線入射率と他の検出器モジュール11のX線入射率を関係付けることができる。図7に示す例では、z方向の検出器番号z=1、z=NzにおけるX線入射率は、それぞれモジュール番号m=Nm、m=1におけるリファレンス検出器の出力を用いる。
図7は、検出器モジュール11を装置に実装する際に、モジュール番号m=1とm=Nmの検出器モジュール11をスライス中心軸よりも意図的にずらして配置したモジュール配列を示す図である。図7に示すように配置することで、最大スライス収集モードにおいても前記実施形態は適用可能となり、検出器モジュール11をずらすことにより、リファレンス検出器に対してコリメータ3の影を作ることができるため、リファレンス検出器のX線入射率と他の検出器モジュール11のX線入射率を関係付けることができる。図7に示す例では、z方向の検出器番号z=1、z=NzにおけるX線入射率は、それぞれモジュール番号m=Nm、m=1におけるリファレンス検出器の出力を用いる。
≪第2の実施形態≫
次に、第2の実施形態について説明する。
図8は、第2の実施形態に係るX線CT装置200の概略構成図である。
X線CT装置200は、検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、焦点位置検出器13、焦点位置計測装置14を含む構成である。
検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
次に、第2の実施形態について説明する。
図8は、第2の実施形態に係るX線CT装置200の概略構成図である。
X線CT装置200は、検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、焦点位置検出器13、焦点位置計測装置14を含む構成である。
検出器1、X線管2、コリメータ3、1次元グリッド4、データ収集装置6、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
焦点位置検出器13は、X線管2(焦点12)の位置を検知し、焦点位置測定装置14は焦点の位置座標を求める。
記憶装置8は、焦点位置測定装置14によって測定された焦点位置と検出器1におけるX線入射率との関係を記憶する。
記憶装置8は、焦点位置測定装置14によって測定された焦点位置と検出器1におけるX線入射率との関係を記憶する。
次に、前処理装置7で実行する、焦点位置に対応した感度補正データの推定方法について説明する。
X線CT装置200の検出器1はモジュール単位で実装されており、X線CT装置100の場合と同様に、1つの検出器モジュール11には図2に示すように検出器1が配置され、検出器モジュール11は図3に示すように円弧状に配置される。
X線CT装置200において、Nsスライス収集モードについて検出器モジュール11毎にずれが発生した場合を考える。
X線CT装置100の場合では、式(7)に示すように、検出器モジュール11毎のX線入射率をリファレンス検出器におけるX線入射率で表したが、第2の実施形態では、検出器モジュール11毎のX線入射率Rを焦点位置fによって表す。即ち、焦点位置検出器13、焦点位置測定装置14によって測定された焦点位置が参照信号となる。
X線CT装置200の検出器1はモジュール単位で実装されており、X線CT装置100の場合と同様に、1つの検出器モジュール11には図2に示すように検出器1が配置され、検出器モジュール11は図3に示すように円弧状に配置される。
X線CT装置200において、Nsスライス収集モードについて検出器モジュール11毎にずれが発生した場合を考える。
X線CT装置100の場合では、式(7)に示すように、検出器モジュール11毎のX線入射率をリファレンス検出器におけるX線入射率で表したが、第2の実施形態では、検出器モジュール11毎のX線入射率Rを焦点位置fによって表す。即ち、焦点位置検出器13、焦点位置測定装置14によって測定された焦点位置が参照信号となる。
まず、コリメータ3の開口幅を全開に設定し、被写体5の無い状態で検出器1の固有感度を収集する(式(4))。次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置fを変えながら被写体5の無い状態で収集する(式(5))。このとき、焦点位置fを焦点位置検出器13及び焦点位置測定装置14によって同時に測定する。こうして式(4)と式(5)から、式(6)が得られる。
次に、各検出器モジュール11におけるX線入射率Rを焦点位置fによって以下のように多項式で近似する。
次に、各検出器モジュール11におけるX線入射率Rを焦点位置fによって以下のように多項式で近似する。
式(10)における係数bnは、参照信号である焦点位置fと検出器モジュール11毎のX線入射率Rとの関係を示す係数である。即ち、この係数bnを予め最小自乗法等で求めておき、記憶装置8に記憶しておけば、被写体5を撮像中においても、各焦点位置fから各検出器モジュール11のX線入射率Rを式(10)より推定することが可能である。
他のスライス収集モードにおいても、係数bnは、端部スライスにおいて検出器モジュール11毎に求める。
求めたX線入射率Rを式(4)に乗じることで、焦点位置に対応した感度補正データを作成する。収集データは、画像再構成処理装置9において、焦点位置fに対応した感度補正データによる感度補正やCT値調整などの各種補正後に画像再構成され、画像表示装置10によって画像が表示される。
求めたX線入射率Rを式(4)に乗じることで、焦点位置に対応した感度補正データを作成する。収集データは、画像再構成処理装置9において、焦点位置fに対応した感度補正データによる感度補正やCT値調整などの各種補正後に画像再構成され、画像表示装置10によって画像が表示される。
第2の実施形態では、焦点位置に対応したX線入射率を式(10)のような連続的な関数によって求めているため、結果として得られる感度補正データX'airの精度も高くなる。
また、式(10)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数bnの数も少なくてすむ。
また、係数bnは装置の据付時などに求めればよい。
また、式(10)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数bnの数も少なくてすむ。
また、係数bnは装置の据付時などに求めればよい。
本第2の実施形態を、特開平4-227238号公報に記載されたような焦点位置に対応してコリメータ3を移動させる方法と組み合わせることも可能である。
式(10)に示す各検出器モジュール11のX線入射率Rを収集する際、コリメータ3の開口の中心は、スライス中心軸上に位置する。そのため、焦点移動に対応してコリメータ3を移動した場合には、測定した焦点位置fとコリメータ開口の中心位置から、コリメータ開口の中心がスライス中心軸上に存在したときに相当する焦点位置f’を算出し、f’を式(10)に代入して各検出器モジュール11におけるX線入射率Rを計算する。このX線入射率Rから、焦点位置に対応したビュー毎の感度補正データが求まる。
コリメータ3を移動させる方法では、焦点の移動を検出してからコリメータ3の移動の命令を発行し、移動が完了するまで遅延が生じる場合がある。この制御の遅延に起因した照射領域のずれの影響は、データ収集終了後、ビュー毎の感度補正データを用いて感度補正することで低減することができる。
式(10)に示す各検出器モジュール11のX線入射率Rを収集する際、コリメータ3の開口の中心は、スライス中心軸上に位置する。そのため、焦点移動に対応してコリメータ3を移動した場合には、測定した焦点位置fとコリメータ開口の中心位置から、コリメータ開口の中心がスライス中心軸上に存在したときに相当する焦点位置f’を算出し、f’を式(10)に代入して各検出器モジュール11におけるX線入射率Rを計算する。このX線入射率Rから、焦点位置に対応したビュー毎の感度補正データが求まる。
コリメータ3を移動させる方法では、焦点の移動を検出してからコリメータ3の移動の命令を発行し、移動が完了するまで遅延が生じる場合がある。この制御の遅延に起因した照射領域のずれの影響は、データ収集終了後、ビュー毎の感度補正データを用いて感度補正することで低減することができる。
≪第3の実施形態≫
次に、第3の実施形態について説明する。
図9は、第3の実施形態に係るX線CT装置300の概略構成図である。
X線CT装置300は、検出器1、X線管2、コリメータ3、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、2次元グリッド15、X線強度検出器16、X線強度測定装置17を含む構成である。
検出器1、X線管2、コリメータ3、データ収集装置6、記憶装置8、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
次に、第3の実施形態について説明する。
図9は、第3の実施形態に係るX線CT装置300の概略構成図である。
X線CT装置300は、検出器1、X線管2、コリメータ3、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、2次元グリッド15、X線強度検出器16、X線強度測定装置17を含む構成である。
検出器1、X線管2、コリメータ3、データ収集装置6、記憶装置8、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
2次元グリッド15は、被写体5から発生するx方向とz方向の散乱線を除去する。
X線強度検出器16は照射されるX線強度を検出し、X線強度測定装置17は、X線強度を測定する。
X線強度検出器16は照射されるX線強度を検出し、X線強度測定装置17は、X線強度を測定する。
次に、前処理装置7で実行する、焦点位置に対応した感度補正データの推定方法について説明する。
X線CT装置300の検出器1は、モジュール単位で実装される。図10は、検出器と2次元グリッドから構成される検出器モジュール11を示す図である。図10に示すように、1つの検出器検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置される。また、被写体5から発生するx方向とz方向の散乱線を除去するための2次元グリッド15が検出器1の前面に装着される。
また、X線CT装置300における検出器モジュール11は、図3に示すように、x方向に円弧状にNm個配置されるものとする。
X線CT装置300の検出器1は、モジュール単位で実装される。図10は、検出器と2次元グリッドから構成される検出器モジュール11を示す図である。図10に示すように、1つの検出器検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置される。また、被写体5から発生するx方向とz方向の散乱線を除去するための2次元グリッド15が検出器1の前面に装着される。
また、X線CT装置300における検出器モジュール11は、図3に示すように、x方向に円弧状にNm個配置されるものとする。
前記の第1の実施形態及び第2の実施形態では、コリメータ3の影が生じる端部スライスにおけるX線入射率のみを考慮していたが、第3の実施形態では、焦点移動によって2次元グリッド15による影も生じるため、全てのスライスにおけるX線入射率を求める。
以下、X線CT装置300において、Nsスライス収集モードの撮像において検出器モジュール11毎にずれが発生した場合を考える。
以下、X線CT装置300において、Nsスライス収集モードの撮像において検出器モジュール11毎にずれが発生した場合を考える。
第1の実施形態と同様に、モジュール番号m=1の検出器モジュール11にリファレンス検出器を設け、リファレンス検出器におけるX線入射率から他の検出器モジュール11におけるスライス毎のX線入射率を求める方法について説明する。
はじめに、コリメータ3の開口幅を全開に設定し、被写体5の無い状態で検出器1の固有感度を収集する。このとき、焦点位置が必ずしもスライス中心軸上にあるとは限らないため、検出器の出力Xairは次式のように焦点位置fに依存した関数となる。
はじめに、コリメータ3の開口幅を全開に設定し、被写体5の無い状態で検出器1の固有感度を収集する。このとき、焦点位置が必ずしもスライス中心軸上にあるとは限らないため、検出器の出力Xairは次式のように焦点位置fに依存した関数となる。
ここで、fairは検出器の固有感度Sを収集したときの焦点位置であり、その位置での
X線入射率をR(z,m,fair)とする。
次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置を変えながら被写体5の無い状態で収集すると、焦点位置fにおける検出器1の出力Xfは以下のようになる。
X線入射率をR(z,m,fair)とする。
次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置を変えながら被写体5の無い状態で収集すると、焦点位置fにおける検出器1の出力Xfは以下のようになる。
ここで、R(z,m,f)は焦点位置fにおけるX線入射率である。
式(11)と式(12)より、X線入射率R’は次式のように表される。
式(11)と式(12)より、X線入射率R’は次式のように表される。
次に、検出器モジュール11毎のX線入射率R'(z,m,f)をリファレンス検出器が存在するモジュール番号m=1の検出器モジュール11のX線入射率R'(z,1,f)によって以下の多項式で近似する。
式(14)における係数cnは、参照信号であるリファレンス検出器が存在する検出器モジュール11のX線入射率と検出器モジュール11毎のX線入射率との関係を示す係数である。即ち、この係数cnを予め最小自乗法等で求めておけば、被写体5を撮像中においても、リファレンス検出器におけるX線入射率R’から各検出器モジュール11のX線入射率R’を式(14)より推定することが可能である。係数cnは、それぞれの検出器モジュール11においてスライス毎に求め、記憶装置8に記憶される。また、他のスライス収集モードにおいても、係数cnは、それぞれの検出器モジュール11においてスライス毎に求められる。
被写体5を撮像している場合、リファレンス検出器におけるX線入射率R’は、検出器の出力値Xを固有感度SおよびX線強度Iで除することで求める。X線強度Iは、X線強度検出器16、X線強度測定装置17で測定される。
このように、リファレンス検出器におけるX線入射率R’が求まれば、各検出器モジュール11におけるX線入射率R’も求めることが可能であり、その値を式(11)に乗じることで焦点位置毎の感度補正データX'airが以下のように得られる。
このように、リファレンス検出器におけるX線入射率R’が求まれば、各検出器モジュール11におけるX線入射率R’も求めることが可能であり、その値を式(11)に乗じることで焦点位置毎の感度補正データX'airが以下のように得られる。
画像再構成処理装置9は、各焦点位置において検出器1によって検出された収集データを、前処理装置7によって作成された感度補正データX'airによる感度補正やCT値調整などの各種補正を行った後に画像を再構成し、画像表示装置10はその画像を表示する。
このように、第3の実施形態では、焦点位置に対応したX線入射率を式(14)のような連続的な関数によって求めているため、結果として得られる感度補正データX'airの精度も高くなる。
また、式(14)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数cnの数も少なくてすむ。
また、式(14)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数cnの数も少なくてすむ。
前記の方法では、被写体5を撮像している場合のX線強度IをX線強度検出器16とX線強度測定装置17によって求めたが、X線強度を他の方法で求めることも可能である。
図11は、検出器と一部のグリッドを取り除いた2次元グリッドから構成される検出器モジュール11を示す図である。
図11に示すように、モジュール番号m=1の検出器モジュール11におけるリファレンス検出器に対する散乱線除去用のグリッドを取り除くことで、X線強度Iをリファレンス検出器により測定することができる。この場合、図11に示すグリッドを一部取り除いたモジュール番号m=1のリファレンス検出器はX線強度Iを測定するために利用し、式(13)に示すX線入射率R’は、モジュール番号m=Nmにおけるリファレンス検出器から求めればよい。
図11は、検出器と一部のグリッドを取り除いた2次元グリッドから構成される検出器モジュール11を示す図である。
図11に示すように、モジュール番号m=1の検出器モジュール11におけるリファレンス検出器に対する散乱線除去用のグリッドを取り除くことで、X線強度Iをリファレンス検出器により測定することができる。この場合、図11に示すグリッドを一部取り除いたモジュール番号m=1のリファレンス検出器はX線強度Iを測定するために利用し、式(13)に示すX線入射率R’は、モジュール番号m=Nmにおけるリファレンス検出器から求めればよい。
≪第4の実施形態≫
次に、第4の実施形態について説明する。
図12は、第4の実施形態に係るX線CT装置400の概略構成図である。
X線CT装置400は、検出器1、X線管2、コリメータ3、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、2次元グリッド15、X線強度/焦点位置検出器18、X線強度/焦点位置測定装置19を含む構成である。
検出器1、X線管2、コリメータ3、データ収集装置6、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
次に、第4の実施形態について説明する。
図12は、第4の実施形態に係るX線CT装置400の概略構成図である。
X線CT装置400は、検出器1、X線管2、コリメータ3、データ収集装置6、前処理装置7、記憶装置8、画像再構成処理装置9、画像表示装置10、2次元グリッド15、X線強度/焦点位置検出器18、X線強度/焦点位置測定装置19を含む構成である。
検出器1、X線管2、コリメータ3、データ収集装置6、画像再構成処理装置9、画像表示装置10は、図1に示すX線CT装置100のそれと同様の構成、機能を有する。
2次元グリッド15は、X線が被写体5に照射される際に発生するx方向とz方向の散乱線を除去するものであり、図10に示す構造を有する。
X線強度/焦点位置検出器18は、X線管2の位置(焦点位置)とX線強度を検出し、X線強度/焦点位置測定装置19は、X線管2の位置(焦点位置)とX線管2から照射されるX線強度を測定する。
X線強度/焦点位置検出器18は、X線管2の位置(焦点位置)とX線強度を検出し、X線強度/焦点位置測定装置19は、X線管2の位置(焦点位置)とX線管2から照射されるX線強度を測定する。
次に、前処理装置7で実行する、焦点位置に対応した感度補正データの推定方法について説明する。
X線CT装置400の検出器1はモジュール単位で実装される。図10に示す検出器モジュール11と同様、X線CT装置400の1つの検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置される。また、被写体5から発生するx方向とz方向の散乱線を除去するための2次元グリッド15が検出器1の前面に装着される。
また、X線CT装置400における検出器モジュール11は、図3に示すように、x方向に円弧状にNm個配置されるものとする。
X線CT装置400の検出器1はモジュール単位で実装される。図10に示す検出器モジュール11と同様、X線CT装置400の1つの検出器モジュール11にはz方向にNz個、x方向にNx個の検出器1が配置される。また、被写体5から発生するx方向とz方向の散乱線を除去するための2次元グリッド15が検出器1の前面に装着される。
また、X線CT装置400における検出器モジュール11は、図3に示すように、x方向に円弧状にNm個配置されるものとする。
X線CT装置400では、X線CT装置300と同様に、焦点移動によって2次元グリッド15による影を生じるため、全てのスライスにおけるX線入射率を求める。
以下、X線CT装置400において、Nsスライス収集モードについて検出器モジュール11毎にずれが発生した場合を考え、測定した焦点位置からそれぞれの検出器モジュール11におけるスライス毎のX線入射率を求める方法について説明する。
以下、X線CT装置400において、Nsスライス収集モードについて検出器モジュール11毎にずれが発生した場合を考え、測定した焦点位置からそれぞれの検出器モジュール11におけるスライス毎のX線入射率を求める方法について説明する。
まず、コリメータ3の開口幅を全開に設定し、被写体5の無い状態で検出器1の固有感度を収集する(式(11))。次に、コリメータ3をNsスライス収集モードにおける開口幅に設定し、焦点位置fを変えながら被写体5の無い状態で収集する(式(12))。このとき、焦点位置fをX線強度/焦点位置検出器18、X線強度/焦点位置測定装置19によって同時に測定する。こうして式(11)と式(12)から、式(13)が得られる。
次に、各検出器モジュール11におけるX線入射率R’を焦点位置fによって以下のように多項式で近似する。
次に、各検出器モジュール11におけるX線入射率R’を焦点位置fによって以下のように多項式で近似する。
式(16)における係数dnは、参照信号である焦点位置fと検出器モジュール11毎のX線入射率Rとの関係を示す係数である。即ち、この係数dnを予め最小自乗法等で求めておき、記憶装置8に記憶しておけば、被写体5を撮像中においても、各焦点位置fから各検出器モジュール11のX線入射率R’を式(16)より推定することが可能である。
係数dnは、それぞれの検出器モジュール11においてスライス毎に求め、記憶装置8に記憶される。また、他のスライス収集モードにおいても、係数dnは、それぞれの検出器モジュール11においてスライス毎に求められる。
式(16)により求めたX線入射率R’の値を式(11)に乗じることで焦点位置毎の感度補正データX'airが以下のように得られる。
式(16)により求めたX線入射率R’の値を式(11)に乗じることで焦点位置毎の感度補正データX'airが以下のように得られる。
ここで、X線強度Iは、X線強度/焦点位置検出器18、X線強度/焦点位置測定装置19で測定される。
このように、各焦点位置fにおいて検出器1によって検出された収集データは、前処理装置7によって前記の処理に従って感度補正データX'airとして補正され、画像再構成処理装置9に出力される。画像再構成処理装置9は、感度補正データX'airに対してCT値調整などの各種補正を行った後、画像を再構成し、画像表示装置10は画像を表示する。
このように、各焦点位置fにおいて検出器1によって検出された収集データは、前処理装置7によって前記の処理に従って感度補正データX'airとして補正され、画像再構成処理装置9に出力される。画像再構成処理装置9は、感度補正データX'airに対してCT値調整などの各種補正を行った後、画像を再構成し、画像表示装置10は画像を表示する。
このように、第4の実施形態では、焦点位置に対応したX線入射率を式(16)のような連続的な関数によって求めているため、結果として得られる感度補正データX'airの精度も高くなる。
また、式(16)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数dnの数も少なくてすむ。
また、係数dnは装置の据付時などに測定し、記憶して置けばよい。
また、式(16)における多項式は実際には2次程度の近似で十分であり、記憶装置8に保持しなければならない係数dnの数も少なくてすむ。
また、係数dnは装置の据付時などに測定し、記憶して置けばよい。
1…検出器
2…X線管
3…コリメータ
4…1次元グリッド
5…被写体
6…データ収集装置
7…前処理装置
8…記憶装置
9…画像再構成処理装置
10…画像表示装置
11…検出器モジュール(検出モジュール)
12…焦点
13…焦点位置検出器
14…焦点位置測定装置
15…2次元グリッド
16…X線強度検出器
17…X線強度測定装置
18…X線強度/焦点位置検出器
19…X線強度/焦点位置測定装置
2…X線管
3…コリメータ
4…1次元グリッド
5…被写体
6…データ収集装置
7…前処理装置
8…記憶装置
9…画像再構成処理装置
10…画像表示装置
11…検出器モジュール(検出モジュール)
12…焦点
13…焦点位置検出器
14…焦点位置測定装置
15…2次元グリッド
16…X線強度検出器
17…X線強度測定装置
18…X線強度/焦点位置検出器
19…X線強度/焦点位置測定装置
Claims (6)
- X線を照射するX線管と、前記X線の照射領域を制限するコリメータと、前記X線を検出する複数の検出器を有する1以上の検出モジュールと、前記検出器に設置され、前記X線が被写体に照射された場合に前記被写体からの散乱線を除去するグリッドと、記憶装置と、前記検出器の出力を元に感度補正データを作成する前処理装置と、前記感度補正データを用いて前記検出器の出力を補正して、前記補正した出力に基づき前記被写体の画像を再構成する画像再構成処理装置と、前記画像再構成処理装置で再構成した画像を表示する画像表示装置と、を備えたX線CT装置であって、
前記前処理装置は、
前記被写体が無い状態において、前記X線管の位置となる焦点位置を変え、焦点位置毎における各検出モジュールのX線入射率を取得する手段と、
前記被写体が無い状態において、参照信号と前記取得した各検出モジュールのX線入射率との関係を求める手段と、
前記求めた関係を前記記憶装置に記憶する手段と、
前記被写体の撮像時において、被写体撮像時の参照信号と前記記憶された関係を用いて、各検出モジュールの被写体撮像時のX線入射率を算出する手段と、
前記算出された被写体撮像時のX線入射率を用いて、焦点位置毎の感度補正データを求める手段と、
を有することを特徴とするX線CT装置。 - 前記取得された焦点位置毎の各検出モジュールにおけるX線入射率は、前記参照信号の多項式で近似され、
前記関係は、前記多項式の係数である
ことを特徴とする請求の範囲第1項に記載のX線CT装置。 - 前記参照信号は、前記被写体を通過しないX線を検出する検出器の出力を用いる
ことを特徴とする請求の範囲第1項に記載のX線CT装置。 - 前記参照信号は、前記焦点位置を表す値を用いる
ことを特徴とする請求の範囲第1項に記載のX線CT装置。 - 前記グリッドとして、スライス方向に対して垂直な方向の散乱線を除去する1次元グリッドを用い、
前記被写体の複数スライス収集モードの撮像を行う場合、前記焦点位置毎の感度補正データは端部スライスにおいて求める
ことを特徴とする請求の範囲第1項に記載のX線CT装置。 - 前記グリッドとして、スライス方向及びスライス方向に対して垂直な方向の散乱線を除去する2次元グリッドを用い、
前記被写体の複数スライス収集モードの撮像を行う場合、前記焦点位置毎の感度補正データは全てのスライスにおいて求める
ことを特徴とする請求の範囲第1項に記載のX線CT装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/390,752 US8761483B2 (en) | 2009-09-28 | 2010-08-12 | Calibration system for focal spot shift in an X-ray CT device |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009-223335 | 2009-09-28 | ||
JP2009223335A JP5400546B2 (ja) | 2009-09-28 | 2009-09-28 | X線ct装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011036961A1 true WO2011036961A1 (ja) | 2011-03-31 |
Family
ID=43795721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2010/063714 WO2011036961A1 (ja) | 2009-09-28 | 2010-08-12 | X線ct装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8761483B2 (ja) |
JP (1) | JP5400546B2 (ja) |
WO (1) | WO2011036961A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110621985A (zh) * | 2017-07-03 | 2019-12-27 | 株式会社岛津制作所 | X线计算机断层装置 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5400546B2 (ja) * | 2009-09-28 | 2014-01-29 | 株式会社日立メディコ | X線ct装置 |
JP6316283B2 (ja) * | 2012-06-05 | 2018-04-25 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | X線ctイメージャの動き層分解較正 |
JP2014147418A (ja) * | 2013-01-31 | 2014-08-21 | Ge Medical Systems Global Technology Co Llc | 画像生成装置および放射線断層撮影装置並びにプログラム |
US20160199019A1 (en) * | 2015-01-13 | 2016-07-14 | Arineta Ltd. | Method and apparatus for focal spot position tracking |
CN107582089B (zh) * | 2017-09-29 | 2021-06-29 | 上海联影医疗科技股份有限公司 | 准直器、成像设备、焦点位置跟踪方法及校正方法 |
WO2019090541A1 (en) | 2017-11-08 | 2019-05-16 | Shenzhen United Imaging Healthcare Co., Ltd. | Systems and methods for correcting projection images in computed tomography image reconstruction |
CN109431528B (zh) * | 2017-11-08 | 2022-05-31 | 上海联影医疗科技股份有限公司 | 校正投影图像的***和方法 |
US10898159B2 (en) | 2019-01-11 | 2021-01-26 | General Electric Company | X-ray imaging system use and calibration |
CN111728632B (zh) * | 2020-07-31 | 2023-08-15 | 上海联影医疗科技股份有限公司 | 射线探测装置、射线探测方法和ct图像重建方法 |
CN113749676B (zh) * | 2020-10-20 | 2024-07-12 | 宽腾(北京)医疗器械有限公司 | 一种实现ct精确校直的方法 |
CN113017673B (zh) * | 2021-03-03 | 2022-08-26 | 明峰医疗***股份有限公司 | 一种模体的定位装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06269443A (ja) * | 1993-03-18 | 1994-09-27 | Toshiba Corp | X線ct装置 |
JP2000033085A (ja) * | 1998-07-17 | 2000-02-02 | Shimadzu Corp | X線ct装置 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5550886A (en) * | 1994-11-22 | 1996-08-27 | Analogic Corporation | X-Ray focal spot movement compensation system |
JPH105210A (ja) * | 1996-06-21 | 1998-01-13 | Toshiba Corp | X線ct装置及びそのミスアライメント補正方法 |
JP3527381B2 (ja) * | 1997-02-28 | 2004-05-17 | 株式会社島津製作所 | X線ct装置 |
JP4298011B2 (ja) * | 1998-08-06 | 2009-07-15 | 株式会社東芝 | X線ファンビーム位置移動検出システムおよびそれを適用したctスキャナ装置 |
JP2000093418A (ja) | 1998-09-22 | 2000-04-04 | Toshiba Corp | X線検出装置およびx線イメージング装置 |
US6094469A (en) * | 1998-10-21 | 2000-07-25 | Analogic Corporation | Computed tomography system with stable beam position |
US6979826B2 (en) * | 2002-07-29 | 2005-12-27 | Ge Medical Systems Global Technology Company Llc | Scintillator geometry for enhanced radiation detection and reduced error sensitivity |
US7616801B2 (en) * | 2002-11-27 | 2009-11-10 | Hologic, Inc. | Image handling and display in x-ray mammography and tomosynthesis |
US8063379B2 (en) * | 2006-06-21 | 2011-11-22 | Avraham Suhami | Radiation cameras |
JP2008036284A (ja) * | 2006-08-09 | 2008-02-21 | Toshiba Corp | 医用画像合成方法及びその装置 |
US7912270B2 (en) * | 2006-11-21 | 2011-03-22 | General Electric Company | Method and system for creating and using an impact atlas |
GB0708655D0 (en) * | 2007-05-04 | 2007-06-13 | Imp Innovations Ltd | Imaging of moving objects |
JP2008299982A (ja) * | 2007-06-01 | 2008-12-11 | Hitachi Ltd | 光ディスク、光ディスク装置および光ディスク装置の制御方法 |
JP5153348B2 (ja) * | 2008-01-09 | 2013-02-27 | 株式会社日立ハイテクノロジーズ | 荷電粒子ビーム軌道補正器及び荷電粒子ビーム装置 |
US9995507B2 (en) * | 2009-04-15 | 2018-06-12 | Richard Norman | Systems for cost-effective concentration and utilization of solar energy |
JP5400546B2 (ja) * | 2009-09-28 | 2014-01-29 | 株式会社日立メディコ | X線ct装置 |
US20120278055A1 (en) * | 2009-11-18 | 2012-11-01 | Koninklijke Philips Electronics N.V. | Motion correction in radiation therapy |
US9764109B2 (en) * | 2011-11-19 | 2017-09-19 | Yale University | Method of correlating brain activity |
KR20130069042A (ko) * | 2011-12-16 | 2013-06-26 | 삼성전자주식회사 | 자기 공명 영상 촬영 방법 및 그에 따른 자기 공명 영상 촬영 장치 |
-
2009
- 2009-09-28 JP JP2009223335A patent/JP5400546B2/ja active Active
-
2010
- 2010-08-12 US US13/390,752 patent/US8761483B2/en not_active Expired - Fee Related
- 2010-08-12 WO PCT/JP2010/063714 patent/WO2011036961A1/ja active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06269443A (ja) * | 1993-03-18 | 1994-09-27 | Toshiba Corp | X線ct装置 |
JP2000033085A (ja) * | 1998-07-17 | 2000-02-02 | Shimadzu Corp | X線ct装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110621985A (zh) * | 2017-07-03 | 2019-12-27 | 株式会社岛津制作所 | X线计算机断层装置 |
CN110621985B (zh) * | 2017-07-03 | 2022-03-11 | 株式会社岛津制作所 | X线计算机断层装置 |
Also Published As
Publication number | Publication date |
---|---|
US8761483B2 (en) | 2014-06-24 |
JP5400546B2 (ja) | 2014-01-29 |
US20120177272A1 (en) | 2012-07-12 |
JP2011067554A (ja) | 2011-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5400546B2 (ja) | X線ct装置 | |
EP2276408B1 (en) | X ray imaging system with scatter radiation correction and method of using same | |
EP2046203B1 (en) | X-ray detector gain calibration depending on the fraction of scattered radiation | |
US8699659B2 (en) | Systems and methods for focal spot motion correction | |
US7502440B2 (en) | Radiodiagnostic apparatus | |
CN106725568B (zh) | Ct扫描仪散焦强度测量方法 | |
CN106687045B (zh) | 数据处理装置、x射线ct装置以及参考校正方法 | |
KR20070011176A (ko) | X선 ct 장치 | |
CN102770077B (zh) | X射线摄像装置 | |
US20100158318A1 (en) | Focal spot size measurement with a movable edge located in a beam-shaping device | |
US6876718B2 (en) | Scatter correction methods and apparatus | |
JP2011167334A (ja) | 放射線撮影装置 | |
US7476025B2 (en) | Shadow mask for an X-ray detector, computed tomography unit having a shadow mask, and a method for adjusting a shadow mask | |
JP5526775B2 (ja) | 放射線撮像装置 | |
JP5220383B2 (ja) | 放射線撮像装置 | |
EP3835830B1 (en) | Systems and methods for estimating a focal spot motion and calculating a corresponding correction | |
JP4584550B2 (ja) | X線計測装置 | |
JP2006102299A (ja) | X線線量補正方法およびx線ct装置 | |
JPH11299768A (ja) | X線ct装置 | |
JP2018114211A (ja) | X線ct装置 | |
JP7207856B2 (ja) | Ct撮影装置 | |
JP2013009874A (ja) | X線ct装置 | |
JP5849838B2 (ja) | 放射線断層画像生成装置、放射線断層撮影装置および放射線断層画像生成プログラム | |
JP2017070590A (ja) | 画像処理装置及びその制御方法、コンピュータプログラム | |
JPH1189827A (ja) | X線コンピュータ断層撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 10818637 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13390752 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 10818637 Country of ref document: EP Kind code of ref document: A1 |