WO2019083014A1 - 光子計数型のx線検出データを処理する方法及び装置、並びに、x線装置 - Google Patents
光子計数型のx線検出データを処理する方法及び装置、並びに、x線装置Info
- Publication number
- WO2019083014A1 WO2019083014A1 PCT/JP2018/039884 JP2018039884W WO2019083014A1 WO 2019083014 A1 WO2019083014 A1 WO 2019083014A1 JP 2018039884 W JP2018039884 W JP 2018039884W WO 2019083014 A1 WO2019083014 A1 WO 2019083014A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- ray
- energy
- effective atomic
- atomic number
- bin
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims description 76
- 238000000034 method Methods 0.000 title claims description 52
- 238000001514 detection method Methods 0.000 title description 20
- 238000012937 correction Methods 0.000 claims abstract description 202
- 239000000126 substance Substances 0.000 claims description 64
- 238000003672 processing method Methods 0.000 claims description 33
- 238000001228 spectrum Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 16
- 238000010606 normalization Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 14
- 238000002083 X-ray spectrum Methods 0.000 claims description 12
- 238000007689 inspection Methods 0.000 claims description 10
- 238000012889 quartic function Methods 0.000 claims description 6
- 238000002360 preparation method Methods 0.000 claims description 5
- 230000001747 exhibiting effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 12
- 230000006870 function Effects 0.000 description 61
- 238000010521 absorption reaction Methods 0.000 description 24
- 239000013598 vector Substances 0.000 description 21
- 238000003384 imaging method Methods 0.000 description 18
- 230000005540 biological transmission Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 6
- 230000004907 flux Effects 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 238000009607 mammography Methods 0.000 description 4
- 239000013558 reference substance Substances 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000013480 data collection Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 238000005316 response function Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 239000004925 Acrylic resin Substances 0.000 description 2
- 229920000178 Acrylic resin Polymers 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 239000013256 coordination polymer Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 210000005075 mammary gland Anatomy 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229910004613 CdTe Inorganic materials 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 101100208721 Mus musculus Usp5 gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- 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
-
- 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/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- 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/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/083—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
- G01N23/087—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays using polyenergetic X-rays
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/18—Investigating the presence of flaws defects or foreign matter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1647—Processing of scintigraphic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/20—Sources of radiation
- G01N2223/206—Sources of radiation sources operating at different energy levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/423—Imaging multispectral imaging-multiple energy imaging
Definitions
- the present invention relates to a method and apparatus for processing X-ray detection data transmitted through an object, and an X-ray apparatus, and in particular, X-ray detection of photon counting type by irradiating the object with X-rays having a continuous spectrum. And apparatus for processing the collected data and photon counting X-ray apparatus.
- Patent Document 1 Japanese Patent Laid-Open No. 2010-091483: "Foreign material detection method and apparatus according to the name of the invention"
- This patent is based on a so-called dual energy method (or subtraction method).
- This inspection method utilizes the fact that, when X-rays of two types of energy (ie, two types of X-rays having different wavelengths) pass through a substance, there is a difference in the X-ray transmission information.
- the following processing is basic. First, two types of X-ray images of low energy and high energy are simultaneously created, and their images are mutually subtracted. Further, the component image of the foreign object is extracted from the difference image, and the component image is subjected to threshold processing to detect the foreign object.
- a beam hardening phenomenon (a ray hardening phenomenon) occurs.
- the low energy component is more easily absorbed than the high energy component, and the component ratio after material transmission shifts to the high energy side, and the effective energy is on the high energy side. It refers to the phenomenon of shifting.
- Patent Document 1 The beam hardening phenomenon which is inevitable in this X-ray apparatus is described, for example, in Patent Document 1 and Patent Document 2.
- Patent Document 1 an X-ray generator for irradiating a subject with pulsed X-rays having a continuous spectrum (multicolor) and an X-ray transmitted through the subject as photons are considered.
- An X-ray device is described which comprises a photon counting detector which outputs an electrical signal according to the energy of the particles.
- processing is performed to reconstruct an image in which the occurrence of an artifact due to a beam hardening phenomenon is suppressed.
- Patent Document 2 discloses a process for correcting the influence of the beam hardening phenomenon on the detection signal of the detector in an X-ray apparatus provided with an X-ray source and an X-ray detector. It is illustrated. As this processing, the thickness of the simulated subject is changed to create a plurality of projection data, the projection data for each transmission distance is plotted on a graph, and the plotted projection data and the theoretical value are associated to perform beam hardening correction. I have created a function. Therefore, for example, in the radiation detection apparatus described in Patent Document 1, in order to reduce the beam hardening phenomenon, a correction function for beam hardening correction described in Patent Document 2 is created, and the object to be examined is corrected with this correction function. It can be considered to correct the measured values of transmitted X-rays of
- the beam hardening correction technique using a correction function described in Patent Document 2 has a very narrow range of effective atomic numbers of an object to which one correction function can be applied.
- the effective atomic number Zeff referred to here is the average atomic number of the atomic numbers of the plural types of elements (substances) present in the path through which the X-ray beam passes.
- the application range that is, a certain accuracy
- the degree of the beam hardening phenomenon also changes depending on the thickness t of the material through which the X-ray passes and the magnitude of the X-ray energy. For this reason, with the correction function experimentally or simply estimated based on some substance, although convenient beam hardening correction can be performed, accurate information can not be obtained.
- the present inventors use photon counting X-ray detectors and collect photon counts of X-rays for each of a plurality of BINs (energy ranges) set in the spectrum of continuous X-rays, From this collected data, develop a technology to identify (specify, estimate, and evaluate) changes in the types and properties of substances present in the path without depending on the thickness t of the object (length of the X-ray path) There is.
- this technology generically substance identification. In order to identify this substance, it is necessary to consider what range of atomic number of substance is to be subjected to X-ray detection in the field (medical, foreign substance detection, etc.) (for example, Patent Document 3) reference).
- correction data to replace with is calculated.
- the correction data is calculated for each X-ray energy BIN.
- the actual measurement value that is, the X-ray attenuation amount ⁇ t before correction is corrected for each energy BIN and each pixel as needed.
- This correction is performed for each energy BIN as a process for each of a plurality of pixels forming a focused imaging region, so that the disadvantage that the beam hardening phenomenon changes depending on the amount of X-ray energy can be improved.
- the range of the substance that can be used to perform appropriate beam hardening correction is the effective atomic number In terms of Z eff , in practice, Z eff is about ⁇ 2.
- the present invention can be applied to X-ray photon measurement in which a substance having an element with a wider effective atomic number Z eff is to be measured, and the measurement value can be reduced with less calculation load and beam hardening with higher accuracy.
- Apparatus and method capable of presenting various quantitative image information and / or more accurate substance identification information based on the amount of X-ray attenuation that can be corrected and the beam hardening correction is performed, and the apparatus and method thereof To provide an X-ray apparatus to which
- a beam-like X-ray having energy of continuous X-ray spectrum is irradiated to an object, and the X-ray transmitted through the object is detected
- a method is provided for counting the number of photons of the X-ray for each pixel area consisting of at least one pixel at each of three or more preset X-ray energy BINs, and processing data indicating the count value.
- This method comprises the operation step of calculating, for each X-ray energy BIN and for each pixel area, count data represented by the ratio of the count value in both the absence state and the existence state of the object.
- beam hardening correction for correcting the beam hardening phenomenon received when the X-ray passes through the object is applied to the count data for each pixel area and for each X-ray energy BIN.
- a correction step for obtaining an X-ray attenuation amount ( ⁇ Low t, ⁇ Middle t, ⁇ High t) ( ⁇ is a linear attenuation coefficient and t is a thickness of the object along the X-ray projection direction);
- the low-energy side normalized X-ray attenuation amount ( ⁇ Low / ( ⁇ Low t, ⁇ Middle t) is normalized to normalize the X-ray attenuation amounts ( ⁇ Low t, ⁇ Middle t) ⁇ Low 2 + ⁇ High energy side normalized X-ray attenuation amount (( ⁇ High / (( High )) by normalizing the X-ray attenuation amounts ( ⁇ Middle t, ⁇ High t) of Middle 2 ) 1/2 )
- the pre-processing step sets a desired range (Zmin to Zmax) of effective atomic numbers of elements constituting the composition of the object, and a lower limit thereof from the effective atomic number (Zm) within the desired range
- Zmin to Zmax a desired range of effective atomic numbers of elements constituting the composition of the object
- Zm effective atomic number
- ⁇ t ⁇ is the linear attenuation coefficient and t is the thickness in the X-ray path direction of the substance
- Setting a straight line on the two-dimensional coordinates as a target function when it is assumed that a substance consisting of an element having a number (for example, Zm 7) is irradiated with monochromatic X-rays, and the horizontal axis on the two-dimensional coordinates
- a beam of X-rays having energy of continuous X-ray spectrum is irradiated to an object, and the X-rays transmitted through the object And detecting the number of photons of the X-ray for each pixel area consisting of at least one pixel in each of three or more X-ray energy BINs set in advance, and providing a method of processing data indicating the count value .
- correction information for correcting an X-ray count value that the X-ray has received a beam hardening phenomenon when the X-ray passes through the object the mass of a plurality of substances having known atomic numbers Based on the characteristics of thickness (mass thickness: tt) and the amount of X-ray attenuation ( ⁇ t: ⁇ is the linear attenuation coefficient and t is the thickness of material in the X-ray path direction) at the effective energy of each X-ray energy BIN
- the X-ray attenuation amount is finally obtained by performing the correction for each pixel area using a pre-processing step of preparing in advance for each X-ray energy BIN and the correction information prepared by the pre-processing step. And an attenuation amount processing step of processing the determined X-ray attenuation amount.
- a beam-like X-ray having energy of continuous X-ray spectrum is irradiated to an object, and the X-ray transmitted through the object is detected and preset 3
- An X-ray apparatus is provided which counts the number of photons of the X-ray for each pixel area consisting of at least one pixel at each of one or more X-ray energy BIN and processes the count value.
- the X-ray apparatus comprises correction information for correcting an X-ray count value at which the X-ray has received a beam hardening phenomenon when the X-ray passes through the object, a plurality of substances having known atomic numbers.
- the attenuation amount processing means comprises: an X-ray image generation means for generating an X-ray image by the photon measurement based on the X-ray attenuation amount corrected and finally determined by the correction means; And X-ray image presentation means for presenting
- the X-ray apparatus is an X-ray medical diagnostic apparatus or an X-ray nondestructive inspection apparatus having a configuration for detecting the X-rays by photon counting.
- Zmin to Zmax a desired range
- Zm a desired effective atomic number
- a straight line on the two-dimensional coordinates when assuming that a substance consisting of an element is irradiated with monochromatic X-rays is set as a target function.
- a beam hardening correction function for correcting hardening is stored in advance in the storage unit as the correction information.
- the beam hardening correction function can be calculated according to the above-described procedure, as long as the information on the generalized target function and the residual of the designated effective atomic number in the predetermined effective atomic number range is held. Therefore, even if the range of the effective atomic number set in advance is wider, it is not necessary to obtain the amount of operation in proportion to the area in calculating the beam hardening correction function. That is, beam hardening correction can be performed with less calculation load for an object having a wider range of elements of effective atomic number Z eff .
- the X-ray attenuation ( ⁇ Low t, ⁇ Middle t) of the two energy BINs on the low energy side among the three or more X-ray energies BIN is normalized Low energy side normalized X-ray attenuation ( ⁇ Low / ( ⁇ Low 2 + ⁇ Middle 2 ) 1/2 ) and high energy side two energy BIN X-ray attenuation ( ⁇ Middle t, ⁇ High t) Is normalized to obtain high-energy-side normalized X-ray attenuation ((.mu. High / (.mu. Middle 2 + .mu. High 2 ) 1/2 ) for each pixel area (area consisting of one or more pixels).
- the low energy side normalized X-ray attenuation ( ⁇ Low / ( ⁇ Low 2 + ⁇ Middle 2) ) 1/2) and the high energy side normalized X-ray attenuation amount ( ⁇ High / ( ⁇ Middle 2 + ⁇ High 2) low energy corresponding to each of the 1/2) Is estimated and a side effective atomic number (Z Low) and the high energy side effective atomic number (Z High) for each pixel region.
- the low-energy side effective atomic number (Z Low) and the high energy side effective atomic number (Z High is compared with each other to determine the degree of coincidence of the effective atomic number (Z Low , Z High ), which is a value that can be regarded as coincident or coincident by this coincidence determination.
- the effective atomic number is determined.
- the reference information prepared in advance it may be regarded as a value close to or true from the effective atomic number (Z Low , Z High ) estimated from the low energy side and the high energy side. Since effective atomic numbers having possible values can be introduced, it is possible to estimate the effective atomic number image of the substance with higher accuracy. This makes it possible to more reliably identify the type and nature of the substance present on the X-ray flux transmission path.
- the attenuation characteristics of elements of various atomic numbers are generalized with theoretical atomic number Z as a parameter, and other attenuation characteristics for the designated attenuation characteristics are residual Based on pre-processing of holding as information. Furthermore, using the correction information obtained in the pre-processing, the actual acquired data is subjected to beam hardening correction for each pixel area and for each energy BIN to obtain the correction value. The effective energy numbers on the low energy side and the high energy side are estimated from the beam-hardening-corrected correction values, and the estimated values are arithmetically driven to a true value or a value that can be regarded as such.
- correction information can be obtained and stored relatively easily. Once this correction information has been obtained, it can be used clinically as long as it falls within the scope of its application. That is, since the correction information is not out of the application range and is versatile, it is not necessary to frequently acquire correction information in accordance with clinical conditions as in the prior art. Also, the correction information may actually have the fitting coefficient of the residual function that fills the residual described above in accordance with the effective atomic number Zeff .
- the memory capacity can be reduced as compared with a correction scheme in which a large number of beam hardening correction functions must be prepared for each characteristic.
- the fact that the effective atomic number Z eff can be estimated with high accuracy means that appropriate beam hardening correction can also be searched. Therefore, it is possible to determine the pixel value of the image with higher accuracy. This increases the quantitativity of the image and reduces the non-uniformity between pixels caused by the inherent variation of the pixels.
- FIG. 1 is a block diagram illustrating a schematic configuration of an embodiment of an X-ray apparatus equipped with an image processing apparatus (implementing an image processing method) according to the present invention
- Fig. 2 is a graph illustrating a continuous spectrum of polychromatic X-rays and three energy BINs set in this spectrum
- FIG. 3 is a flowchart illustrating pre-processing performed as part of image processing in an X-ray apparatus according to an embodiment
- FIG. 4 is a flow chart illustrating data collection processing performed in the X-ray apparatus and post processing that forms another part of the image processing
- FIG. 5 is a graph showing a continuous spectrum of polychromatic X-rays used to estimate the effective atomic number Z eff simulated in consideration of the response function of the detector.
- FIG. 7 illustrates the theoretical attenuation characteristics at the time of monochromatic X-ray irradiation (the same figure (A): solid line: linear characteristics) and the theoretical attenuation characteristics at the time of polychromatic X-ray irradiation, taking low energy BIN as an example.
- A): dotted line: generalization in the horizontal axis direction, target function (FIG. (B): solid line: linear characteristic of X Y) and beam hardening correction curve (FIG. (B): dotted line : A graph explaining the process of creating a curve characteristic)
- FIG. 8 shows the residuals along the vertical axis of the beam hardening correction curve for a specific atomic number Z and the other beam hardening correction curves among the generalized beam hardening correction curves for a plurality of atomic number substances.
- a graph that illustrates FIG. 9 is a graph illustrating the fitting coefficients of a fitting function for correcting the residual in the direction of the vertical axis at low energy BIN
- FIG. 10 is a graph illustrating fitting coefficients of a fitting function for correcting the residual in the vertical axis direction at medium energy BIN
- FIG. 11 is a graph illustrating a fitting coefficient of a fitting function for correcting the residual in the direction of the vertical axis in high energy BIN;
- FIG. 12 is a table illustrating fitting coefficients stored as part of beam hardening correction information
- FIG. 14 is a schematic view of an X-ray image on which an ROI for setting a range for performing beam hardening correction is superimposed;
- FIG. 15 is a diagram showing the ratio of incident X-ray count to output X-ray count (ratio of absorption to transmission) for each energy BIN;
- FIG. 17 is a schematic view of the beam hardening correction information stored in the memory and the estimated value of the effective atomic number Z eff for each pixel and for each energy BIN.
- FIG. 18 is a diagram for explaining a method of estimating the effective atomic numbers Z High and Z Low on the low energy side and the high energy side from the X-ray attenuation normalized for each pixel, and a comparison thereof.
- 19 is a view schematically explaining an example of an effective atomic number image based on an X-ray attenuation amount for each pixel after beam hardening correction and an example of displaying the image.
- FIG. 20 is a schematic flowchart showing an overview of image display processing executed in the second embodiment of the present invention in conjunction with the flow of processing of the first embodiment; FIG.
- FIG. 21 is a diagram schematically illustrating storage contents for image display formed in an image memory
- FIG. 22 is a view showing an example of an image unique to photon counting displayed in the second embodiment
- FIG. 23 is a view showing another example of the photon counting unique image displayed in the second embodiment
- FIG. 24 is a diagram showing another example of the photon counting unique image displayed in the second embodiment.
- FIG. 25 is a diagram illustrating the relationship between the low energy side effective atomic number (Z Low ), the high energy side effective atomic number (Z High ), and the effective atomic number Zm preset when performing beam hardening correction; It is.
- a method and apparatus for processing counting data in photon counting X-ray detection (hereinafter referred to as data processing method and data processing apparatus) according to the present invention will be described with reference to FIGS. 1 to 19.
- the data processing method and the data processing apparatus are, for example, a method and an apparatus mounted or mounted on an X-ray mammography apparatus, a medical X-ray apparatus such as a dental X-ray imaging apparatus, an X-ray apparatus for foreign substance inspection.
- FIG. 10 An example of the structure which comprises the principal part of the X-ray apparatus 10 is shown in FIG.
- the hardware configuration of the data processing method and the X-ray apparatus on which the data processing apparatus according to the first embodiment is mounted or mounted may be a known one, and only the main part will be described.
- the X-ray apparatus 10 includes an X-ray generator 21 which generates X-rays of a continuous spectrum, collimates the X-rays into a beam, and irradiates the object space OS.
- the X-ray generator 21 is disposed at the output side of the X-ray tube 22 driven by receiving a high voltage supply, and the X-ray tube 22 emits X-rays generated by the X-ray tube 22 in the form of a beam.
- a collimator 23 for collimating the light.
- the focal point diameter of the tube focus F of the X-ray tube 22 is, for example, 0.5 mm ⁇ , and this tube focus F can be regarded as an almost point-like X-ray source.
- the X-rays emitted from the tube focal point F are made of a flux of photons (photons) having various energies (X-ray energy), and have a continuous energy spectrum according to the tube voltage.
- the X-ray apparatus 10 further includes a detector 24 for detecting beam-shaped X-rays transmitted through the imaging target OB located in the object space OS.
- the detector 24 has a detection layer 25 having a semiconductor (such as CdTe or CZT) that converts X-rays directly into an electric signal directly under the incident window, and the detection layer 25 has a size of 200 ⁇ m ⁇ 200 ⁇ m, for example.
- a pixel group in which the pixels having the pixel are two-dimensionally arrayed is formed.
- the detector 24 further has, for example, a layered data acquisition circuit 26 built in an ASIC, which processes the detection signal of each pixel for each pixel, for example, on the opposite side of the tube focus F of the detection layer 25.
- the data acquisition circuit 26 is configured as a photon counting circuit capable of counting the number of X-ray photons (photons) incident on the pixel group of the detection layer 25 for each pixel.
- this circuit can divide the X-ray spectrum into a plurality of X-ray energy ranges (referred to as BIN) by setting the threshold for discriminating the X-ray energy, and can count photons for each pixel at each energy BIN It has become.
- the count data created by processing the electric pulse signal in response to the incidence of the X-ray photon is frame data (a set of count data of each pixel Output as The frame rate, for example, varies from 300 fps to 6,600 fps. Except for the superposition phenomenon of photons incident on one pixel, for example, one electric pulse is excited when one photon is incident, so that the count data of each pixel reflects the number of the electric pulses.
- this detector 24 is classified as a photon counting type detector from the viewpoint of the detection processing method. That is, the detector 12 considers the number of photons to be X-rays, assuming that X-rays having continuous energy spectra (polychromatic X-rays) are a set of photons having various energies. And each pixel (note that one or more pixels may be counted).
- this energy BIN as described later, for example, as shown in FIG. 2, three energy BINs: Bin Low to Bin High are set.
- the number of energy BIN: Bin may be four or five as long as it is three or more.
- the region below the lower threshold TH1 of the energy [keV] and the region above the upper threshold TH4 are regions which can not be measured or are not used. Therefore, the region between the threshold values TH1 to TH4 is divided into one (in this case, the threshold values are TH1 and TH4 only) or a plurality of energy BINs. For example, threshold values TH2 and TH3 are set as shown in FIG. 2, and three energy BINs are formed.
- the imaging target OB located in the object space OS is scanned by beam-shaped X-rays.
- a belt conveyor is arranged to pass through the object space OS .
- the imaging target OB is X-ray scanned by placing the imaging target OB on the belt conveyor.
- there is a dental panoramic X-ray imaging apparatus In this case, a patient jaw as an imaging object OB is positioned in the object space OS between the X-ray generator 21 and the detector 24. .
- the jaws are X-ray scanned by rotating the pair of the X-ray generator 21 and the detector 24 while facing each other.
- the imaging target OB may be scanned with relative movement between the pair of the X-ray generator 21 and the detector 24 and the imaging target OB.
- the measurement data of the digital quantity output from the detector 24 is processed by utilizing the superiority of the energy discrimination by the processing apparatus mounted on the X-ray apparatus 10 or the processing apparatus placed outside the X-ray apparatus 10.
- This processing includes image reconstruction by tomosynthesis, creation of an absorption vector length image (two-dimensional image) based on the reconstructed image, and creation of a three-dimensional scatter chart based on the reconstructed image. These processes are proposed by the international publication number WO2016 / 171186 A1 etc.
- the measurement data of the digital quantity output from the detector 24 is also subjected to processing unique to the present invention.
- beam hardening which makes correction result estimated by assumed beam hardening correction curve retrospectively asymptotically approach to true value, and copes with a smaller amount of calculation for wider effective atomic number Z eff of substance. Correction is included.
- the beam hardening correction curve is a curve used to match the attenuation characteristics of the material undergoing beam hardening with the target function (linear attenuation characteristics with respect to the mass thickness when a monochromatic X-ray is irradiated). It is used in the sense, and the measurement value is corrected according to the difference between the beam hardening correction curve and the target function.
- the beam hardening phenomenon (linear hardening phenomenon) is generated when the irradiated polychromatic X-ray passes through the substance, and the low energy component is higher than the high energy component. It is apt to be absorbed and scattered inside, and the component ratio after material permeation is a phenomenon in which the high energy side increases and the effective energy shifts to the high energy side.
- This beam hardening phenomenon can be summarized as what occurs physically as a result of the fact that the interaction between the molecule (atom) of the object and the X-ray photons is different due to the difference in the energy of the X-ray photons.
- the effective atomic number Z eff is an average value taking into consideration the amount of interaction with the X-ray of the atomic number possessed by each of plural kinds of elements (substances) present in the path through which the beam-like X-ray passes. It is an atomic number.
- the X-ray apparatus 10 includes a data processor 30.
- the data processing apparatus 30 is constituted by a computer CP as an example.
- the computer CP itself may be a computer having a known computing function, and includes an interface (I / O) 31 connected to the detector 24 via the communication line LN.
- the interface 31 includes a buffer memory 32, a read-only memory (ROM) 33, a random access memory (RAM) 34, a processor 35 having a central processing unit (CPU) 35A, and an image memory 36 via an internal bus B.
- the input unit 37 and the display unit 38 are communicably connected to each other via the bus B.
- the ROM 33 stores various programs for computer readable measurement value correction, substance identification and the like in advance. For this reason, the ROM 33 includes a storage area (which functions as a non-transitory computer recording medium) 33A in which those programs are stored in advance. Further, the ROM 33 also includes first and second storage areas 33B and 33C for storing beam hardening correction data (also referred to as calibration data) for beam hardening correction of measurement values, which will be described later.
- beam hardening correction data also referred to as calibration data
- the processor 35 (i.e., the CPU 35A) reads a necessary program from its storage area 33A of the ROM 33 into its own work area and executes it.
- the processor 35 is a CPU for image processing.
- the buffer memory 32 is used to temporarily store frame data sent from the detector 24.
- the RAM 34 is used to temporarily store data necessary for the operation when the processor 35 operates.
- the image memory 36 is used to store various image data and information processed by the processor 35.
- the input unit 37 and the display unit 38 function as a man-machine interface with the user, and the input unit 37 receives input information from the user.
- the display 38 can display an image or the like under the control of the data processor 35.
- the data processing device 30 may be provided as a device or inspection system integrated with the X-ray device 10. Further, when the data processing device 30 is communicably connected to the X-ray device 10 via the communication line LN as in the present embodiment, the data processing device 30 may always be connected online, or may be necessary. It may be possible to communicate only at certain times. Furthermore, the data processing device 30 may be provided in a stand-alone manner. Of course, the data processing device 30 may be configured by a hardware circuit that performs pipeline processing and the like. ⁇ Data processing apparatus and data processing method>
- the beam hardening correction according to the present invention is a wider range (Zmin to Zmax) of the effective atomic number Z eff exhibited by a plurality of substances (materials) present in the X-ray irradiation path, as compared with the conventional beam hardening correction. Another advantage is that it is possible to cope with correction data created in a single preparation, and beam hardening correction can be performed with higher accuracy.
- step S1 theoretical values of the energy spectrum of continuous X-rays (polychromatic X-rays) to be irradiated to the X-ray tube 22 are prepared (FIG. 3, step S1).
- This theoretical spectrum is shown in FIG. 5 (A).
- an energy spectrum is created (folded) in consideration of a response function based on the semiconductor material of the detector 12, its thickness, pixel size and the like by Monte Carlo simulation (step S2).
- step S2 The energy spectrum which folded the theoretical spectrum by the response function is illustrated to FIG. 5 (B).
- Three energy BINs energy range are set in this energy spectrum (step S3).
- lower energy side BIN Bin Low (15-23 keV (20.2 keV), middle energy BIN: Bin Middle (23-32 keV (27.6 keV), and higher energy side BIN: Bin High (32-50 keV)
- the numerical values of 20.2 keV, 27.6 keV, and 38.1 keV in parentheses are the effective energy in each BIN, and the mass attenuation coefficient ⁇ / ⁇ ( ⁇ : linear attenuation coefficient corresponding to this effective energy) , ⁇ : mass density) corresponds to the slope of the target function of each BIN.
- step S5 the effective energy of each of the three energy BIN: Bin Low to Bin High is calculated (step S5). Furthermore, the above-mentioned mass attenuation coefficient ⁇ / ⁇ corresponding to the effective energy of each energy BIN: Bin Low ( ⁇ Bin High ) and provisionally applied to the reference effective atomic number Zm is the target function It is calculated as a slope (step S6).
- step S7 an operation of plotting a beam hardening correction curve on two-dimensional coordinates is performed for each of the energy BIN: Bin Low to Bin High (step S7).
- This plot is made using data that is already physically known (theoretical operation value). Specifically, in this plot, as shown in FIGS. 6A to 6C, the abscissa represents the mass thickness (.rho.t), and the ordinate represents the X-ray attenuation amount (.mu.t to the effective energy at each energy BIN). ) Is made.
- the linear theoretical attenuation characteristics (corresponding to the target function) obtained when the
- the value of the horizontal axis ( ⁇ t) is multiplied by the inclination of the above-described target function: mass attenuation coefficient ⁇ / ⁇ to generalize the value of the horizontal axis (step S8).
- the mass attenuation coefficient ⁇ / ⁇ is a known value determined by the atomic number Z and the X-ray energy.
- FIG. 7 illustrates low energy BIN: Bin Low , but according to this generalized operation, the graph of [Before] in FIG. 7 (A) and the graph of [After] in (B). Are generated in operation.
- the atomic number Z (effective atomic number Z eff ) of the substance is different, the degree of occurrence of the beam hardening phenomenon also differs. Therefore, the atomic number Z is a variable for each curve (straight line).
- the processor 35 performs residual calculation and fitting processing (step S9).
- FIGS. 8A, 8B, and 8C show residual (ratio) curves of BIN: Bin Low , Bin Middle, and Bin High , respectively.
- the residual is greater for the low energy BIN than for the higher energy side BIN. This is because photons with lower X-ray energy receive more beam hardening.
- a quartic function is illustrated as a fitting function, it may not necessarily be a quartic function, and a higher order function can be used to fit a residual with higher accuracy. What degree of the function to use may be determined in balance with the amount of operation.
- the processor 35 expresses the residual curve as a two-variable function of “mass thickness ( ⁇ t) ⁇ ⁇ ” and atomic number Z (step S10).
- the coefficients (a0, a1, a2, a3, a4) of the quartic function relating to the fitting described above are expressed as a function of atomic number Z, and the beam hardening correction curve for any atomic number Z in atomic number Zmin to Zmax within the desired range Make it possible to estimate
- the two-variable function is illustrated in FIGS. 9 to 11 for each energy BIN.
- a fitting function f ( ⁇ t) that comprehensively represents the two-variable function (that is, the residual) described above is created as follows (step S11).
- f ( ⁇ t) a 0 + a 1 ⁇ ( ⁇ t) + a 2 ⁇ ( ⁇ t) 2 + a 3 ⁇ ( ⁇ t) 3 + a 4 ⁇ ( ⁇ t) 4 (1)
- An example of the fitting coefficients M j and a j is shown for each energy BIN in FIGS. 12 (A), (B) and (C).
- the function f (.rho.t) is multiplied (step S12).
- the beam hardening correction function as final beam hardening correction information is within the predetermined atomic number range for each energy BIN. It is obtained by
- the processor 35 stores the beam hardening correction function obtained as described above as information for beam hardening correction in, for example, the first storage area 33B (or the second storage area 33C) of the ROM 33 (step S13).
- the corresponding beam hardening correction function is read out by the processor 35 into the work area. Beam hardening correction is performed using the read correction function as in the conventional case. An example of this execution will be described later.
- the fitting coefficients M j and a j may be stored in advance in, for example, the first storage area 33B (or the second storage area 33C) of the ROM 33. In that case, if necessary, they are read out, a beam hardening correction function according to the required atomic number Z is calculated, and a beam hardening correction is performed using it.
- Part 2 Acquisition processing and subsequent processing (including beam hardening correction and creation of effective atomic number image)>
- the processor 35 can perform X-ray transmission data collection and subsequent processing as shown in FIG. 4 interactively with the user at any time.
- the processor 35 controls the high voltage generator (not shown) to drive the X-ray tube 22 while moving the pair of the X-ray generator 21 and the detector 24 relative to the object OB and driving the X-ray tube 22 24 is driven to scan the object OB with a beam of X-rays.
- the high voltage generator (not shown) to drive the X-ray tube 22 while moving the pair of the X-ray generator 21 and the detector 24 relative to the object OB and driving the X-ray tube 22 24 is driven to scan the object OB with a beam of X-rays.
- the pair of the X-ray generator 21 and the detector 24 is fixed, and the object OB is moved to pass through the object space OS.
- the pair of the X-ray generator 21 and the detector 24 are rotated opposite to each other around the patient's jaw as the object OB.
- a tomosynthesis method is applied to the acquired frame data to create an optimum focus image when the object OB is viewed along the X-ray irradiation path (step S22) .
- This optimal focus image may be an image along a cross section of a certain height (depth) of the object OB, or a collection of height (depth) pixels exhibiting the optimum focus in the X-ray irradiation pass for each pixel It may be an image. Of course, it may be a known vector-length image proposed by the present inventor. Alternatively, the transmission image may be simply a scanogram.
- the collected data used to create these images may be used by averaging any one or more of the three energy BIN: Bin Low , Bin Middle , and Bin High selectively or entirely.
- the processor 35 causes the display 38 to display the image created in step S22 (step S23), and interactively or automatically sets an ROI (region of interest) on the image (step S24).
- An example of this is schematically shown in FIG. According to the figure, the ROI designated by the user is shown on the image IM OB of the object OB displayed on the display 38.
- the range of this ROI is the object of beam hardening correction and effective atomic number image described later. For this reason, this range specification is not necessarily required, and the entire image area of the displayed object OB may be set as a target range by default. Also, the number of ROIs may be more than one.
- the processor 35 calculates ⁇ t (attenuation amount) from the actual number of emitted photons for each energy BIN based on the frame data collected by placing the object OB in the object space OS (step S25). This is calculated in advance using, as calibration data, frame data collected only in the air, without placing the object in the object space OS.
- X-ray energy BIN: Bin Low to Bin High shown in FIG. 2 is schematically shown on the horizontal axis of FIG. 15, and the vertical axis is counted at each energy BIN: Bin Low to Bin High.
- the counts of x-ray photons are shown as frequency.
- ⁇ Low , ⁇ Middle and ⁇ High are virtual mean line attenuation coefficients at each energy BIN: Bin Low to Bin High (that is, linear attenuation coefficients with respect to effective energy at each energy BIN), and t is X in the object It is the length (thickness) in the transmission direction of the wire bundle.
- the virtual average linear attenuation coefficients ⁇ Low ( ⁇ Middle , ⁇ High ) of each energy BIN: Bin Low ( ⁇ Bin High ) do not depend on the thickness t.
- the processor 35 designates, as an initial value, the effective atomic number Z eff of one or more elements which are presumed to constitute or be included in the object OB when performing the beam hardening correction.
- one or more types of objects OB that are currently targeted for data collection for example, one or more types of breasts that may be present in a designated portion (eg, a three-dimensional partial region) of the ROI.
- the processor 35 specifies the position of one or more pixels constituting the ROI (step S27), and specifies the initial pixel position P (step S28).
- the count value (count) in the first pixel to be processed is subjected to beam hardening correction for each energy BIN (step S31).
- the count value (count) in the first pixel is subjected to beam hardening correction.
- the corrected attenuation amount is stored in the image memory 36.
- This beam hardening correction is shown in FIG.
- the dimensions of the horizontal axis and the vertical axis are the same as in FIG. 13, and both are dimensions of the X-ray attenuation amount ⁇ t.
- the measured value (dimension of attenuation) ⁇ t n of the first pixel at low energy side BIN: Bin Low is taken on the vertical axis
- the measured value ⁇ t n is affected by beam hardening.
- This calculation may be performed by estimating the intersections of the curves as shown in FIG. 16, or may be estimated from the ratio of the measured value ⁇ t n to the difference ⁇ t n .
- the beam hardening correction counts ⁇ Low t, ⁇ Middle at the first pixel in the ROI at each of the three energies BIN Low , Bin Middle and Bin High. t, ⁇ High t is calculated. This calculated value is temporarily stored, for example, in the image memory 36.
- the storage data in the memory 36 is schematically shown in FIG.
- the processor 35 performs normalization using the corrected counts ⁇ Low t, ⁇ Middle t, and ⁇ High t (step S32). This normalization process is performed for each pixel using the following equation.
- Bin Low Bin Low
- Bin Middle the factor of the length of the X-ray path (the thickness t of the object) passing through the object OB is eliminated, and the attenuation ⁇ High independent of the thickness t -nor is required.
- normalization is performed using two BINs on the low energy side: Bin Low and Bin Middle attenuation amounts ⁇ Low t and ⁇ Middle t.
- the normalized attenuation amounts ⁇ High-nor and ⁇ Low-nor are also stored in the image memory 36 for each pixel forming the ROI (see FIG. 17B).
- the processor 35 determines whether or not the beam hardening correction of the pixel value (that is, the attenuation amount) has been completed for all the pixels in the ROI (step S33). When the determination is NO, it is recognized that the pixel to be subjected to beam hardening correction still remains, the pixel flag is updated (not shown), and the processes of steps S31 and S32 are repeated. By this repetition, the above-described beam hardening correction and normalization are performed on all pixels in the ROI.
- step S33 If the processor 35 determines YES in step S33, it recognizes that the beam hardening correction and normalization have been completed for all the pixels in the ROI, and then proceeds to estimate the effective atomic number Zeff .
- the processor 35 reads out to the work area a reference curve which defines the relation of “amount of attenuation: atomic number” prepared in advance by theoretical calculation in the ROM (step S34).
- This reference curve shows the relationship of the amount of attenuation theoretically obtained from X-ray transmission data collected by irradiating a material having a known atomic number with X-rays.
- the processor 35 estimates the atomic numbers Z High and Z Low using the reference curves related to the two linear attenuation amounts ⁇ High-nor and ⁇ Low-nor (Step S35). Specifically, of the two linear attenuation amounts ⁇ High-nor and ⁇ Low-nor , estimation is made from the high energy side by applying the linear attenuation amount ⁇ High-nor to the reference curve on the high energy side (upper in FIG. 18) Calculate the atomic number Z High . Similarly, the linear attenuation amount ⁇ Low -nor is applied to the reference curve on the low energy side (FIG. 18 lower) to obtain the atomic number Z Low estimated from the low energy side. The obtained atomic numbers Z High and Z Low are stored in the image memory 36 (see FIG. 17C).
- step S35 is performed for each of the pixels forming the ROI.
- step S40 the latest amount of attenuation that has been beam-hardened for each pixel with respect to each energy BIN is stored while being updated.
- step S37 the processor 35 again determines whether or not the work of determining the effective atomic number Z eff has been completed for all the pixels (or one or more pixel regions obtained by bundling a plurality of pixels) in the ROI. (Step S43,). If the determination is NO, the process of step S34 and subsequent steps for the next pixel is repeated after the instruction of step S44. As a result, for all pixels, the effective atomic number Z eff that can be regarded as the average value of the atomic numbers of one or more elements present in the path of the X-ray flux incident on the pixels is determined. Are stored while being updated (see FIG. 17D).
- the degree of beam hardening that a substance having that atomic number receives depends on the magnitude of the energy of the X-ray photon. For this reason, specifying the atomic number of the substance from the linear attenuation coefficient of the specific energy BIN may increase the error, and therefore, the beam hardening correction may not be accurately performed. Therefore, as described above, the normalized linear attenuation is calculated on the high energy side and the low energy side, and the difference between the two attenuations is used to add the atomic number to a true value or a value that can be regarded as a true value. The high effective atomic number Z eff is estimated.
- this background image may be an image other than the focused image by the tomosynthesis method, or the effective atomic number image IM Zeff may be displayed alone.
- the above procedure It can calculate beam hardening correction function. Therefore, even if the range of the effective atomic number set in advance is wider, it is not necessary to obtain the amount of operation in proportion to the area in calculating the beam hardening correction function. That is, beam hardening correction can be performed with less calculation load for an object having a wider range of elements of effective atomic number Zeff .
- the beam hardening correction for each pixel described above actually causes not only the beam hardening phenomenon in a narrow sense, but also factors such as X-ray attenuation such as heel effect, charge sharing, etc. From a measured value including an error due to a circuit factor, correction in a broad sense including the error etc. is made. As a result, the measurement value is corrected as if it had been calibrated from the beginning, and its reliability is enhanced. This means that when performing image reconstruction and object analysis based on the measured values, the processing becomes more stable and more reliable. In the case of identifying the type or property of a substance based on measured values, the identification accuracy becomes high. Second Embodiment
- This X-ray apparatus relates to another application of the beam-hardened corrected count value in the first embodiment described above. Therefore, the same or equivalent components as those described in the first embodiment are denoted by the same reference numerals, and the description thereof will be omitted or simplified.
- step S40 causes the image memory 36 to perform beam hardening correction for each energy BIN in each of all target pixels (or all pixel areas)
- the stored count value is stored.
- the processor 35 reads the count values for each of the beam-hardened energy BINs from the image memory 36 into the work area, and determines them as the final count values for X-ray image creation (see FIG. 20, step S61).
- the stored contents of the determined count value can be schematically shown, for example, as shown in FIG.
- the finally determined counts are ⁇ Low t, ⁇ Middle t, ⁇ High t for low energy BIN: Bin Low , medium energy BIN: Bin Middle , and high energy BIN: Bin High , respectively. It shall appear.
- the processor 35 determines what X-ray image is to be displayed (presented) interactively with the user (step S62).
- this X-ray image in addition to a focused image (including a panoramic image) optimized for each pixel by tomosynthesis based on the count value, that is, the X-ray attenuation amount ⁇ t, quantitativeity of pixel information based on photon measurement
- a three-dimensional scatter plot, an absorption vector length image, and an average absorption value image are provided as unique images in pursuit of. That is, the user selectively designates, by default or interactively, a focused image, a three-dimensional scatter plot, an absorption vector length image, and an average absorption value image.
- the X-ray attenuation amount ⁇ t may be simply included, for example, in a gray scale encoded X-ray transmission image in the option.
- linear attenuation value ⁇ t has three degrees of freedom. Because of this, 3D linear attenuation vector ( ⁇ Low t, ⁇ Middle t, ⁇ High t) Set its length, ie, linear attenuation value vector length (( ⁇ Low t) 2 + ( ⁇ Middle t) 2 + ( ⁇ High t) 2 ) 1/2
- a normalized three-dimensional source weak value vector (hereinafter referred to as a linear attenuation vector) having a denominator as a denominator is a component of the length (thickness t) of the X-ray flux passing through the substance disappearing from the linear attenuation vector itself ( ⁇ Low , ⁇ Middle , ⁇ High ) / ( ⁇ Low 2 + ⁇ Middle 2 + ⁇ High 2 ) 1/2 (5) Calculated as When three axes orthogonal to one another define ⁇ Low t, ⁇ Middle t, and
- this three-dimensional line attenuation vector is calculated for each pixel and mapped to the above three-dimensional coordinates, those end points are distributed around a predetermined point on the surface of the sphere in a certain range around it, including statistical errors
- These scattered points are distributed as a set of The inventors refer to a three-dimensional mapping diagram depicting such scatter points as a three-dimensional scatter diagram.
- An example of this three-dimensional scatter plot is shown in FIG. In the figure, reference symbol Vr indicates a three-dimensional line attenuation vector, and reference symbol DP indicates a scatter point.
- This spherical surface that is, the distribution of the end points of the line attenuation vectors in the three-dimensional scattergram, is unique to the type of substance that constitutes the object itself. That is, if the type of substance is different, the distribution position is theoretically different, so that the type of substance can be identified. ⁇ About the absorption vector length image>
- the linear attenuation value vector length at each pixel is t ( ⁇ Low 2 + ⁇ Middle 2 + ⁇ High 2 ) 1/2 (6)
- the present inventors call this scalar quantity the absorption vector length (or pseudo absorption value).
- a two-dimensional image can be created with this absorption vector length as a pixel value, and the present inventors call this two-dimensional image an absorption vector length image (or a pseudo absorption image).
- An example of this absorption vector length image is schematically shown in FIG. ⁇ About the average absorption value image>
- the virtual mean line attenuation coefficients at three energies BIN: Bin Low , Bin Middle , and Bin High ie, the linear attenuation coefficients for the effective energy of each energy BIN
- pixel value t (a 1 ⁇ Low + a 2 ⁇ Middle + a 3 ⁇ High ) / 3 ...
- the reason why 3 is allocated to the denominator is to calculate an average value over three energy BINs: Bin Low , Bin Middle , and Bin High , that is, all energy BINs.
- coefficients a 1 , a 2 , and a 3 may be values fixed in advance by default, or may be variable while the user or the like interprets an image.
- Condition of coefficient: a 1 + a 2 + a 3 3 is the case of taking a weighted average, and the condition of this coefficient may be removed when handling as a pixel value of real number multiple of weighted average value.
- the present inventors defined an image having the pixel value calculated in this manner in each pixel as an average absorption value image.
- This average absorption value image is schematically shown in FIG.
- Each pixel has a pixel value calculated by the equation (7) or (8) described above.
- the pixel value may be a value calculated by bundling a predetermined number of pixels around a certain pixel.
- the average absorption value image which concerns on this invention is not necessarily limited to what is produced when three X-ray energy BIN is cut out from a continuous X-ray spectrum.
- the number of the plurality of X-ray energies BIN may be two or four or more of a low energy range and a high energy range obtained by dividing the continuous X-ray spectrum according to the magnitude of the X-ray energy.
- the processor 35 converts the image data of the X-ray image designated from among the various types of X-ray images as mentioned above into a predetermined or desired display format (step S63), and displays it on the display 38 (Display) (step S64).
- step S65 the processor 35 interactively determines whether or not to display another aspect (step S65), and repeats the process of step S63 as necessary.
- a unique X-ray image that is, a three-dimensional image, is obtained by utilizing the superiority of photon counting and performing accurate beam hardening correction on the count value.
- the scatter diagram, the absorption vector length image, and the average absorption value image can be selectively displayed as appropriate.
- the beam hardening correction also has an effect of absorbing the fluctuation due to the heel effect on the X-ray tube side. Therefore, it is possible to provide very useful image information in identifying (specifying, estimating) with high accuracy and stably the type and property of the substance (element) forming the object using these images. This is very effective clinically.
- the estimation means and the coincidence determination means are functionally configured.
- pre-steps necessary for processing in the X-ray apparatus, first and second correction steps, attenuation amount processing steps, and X-ray image presentation steps, and pre-preparation means are functionally configured.
- the energy BIN of the X-ray is not necessarily limited to three, and the threshold of energy discrimination may be increased to set four or more energy BINs.
- the threshold of energy discrimination may be increased to set four or more energy BINs.
- two low energy BINs are used to calculate the normalized effective atomic number Z Low of the low energy side as described above, and two BINs of the high energy side are calculated.
- the normalized effective atomic number Z High of the high energy side may be calculated in the same manner as described above.
- the normalized X-ray attenuation amount on the low energy side and the high energy side is ⁇ ⁇ Middle
- a variant is also possible in which calculation is performed as / ( ⁇ Middle 2 + ⁇ High 2 ) 1/2 and ⁇ High / ( ⁇ Low 2 + ⁇ Middle 2 ) 1/2 .
- the unit for calculating the effective atomic number Z eff is not necessarily limited to the configuration in which physical pixels defined in the detection layer of the detector are the unit, and pixel signals of a plurality of such pixels are bundled. It is also possible to virtually use the pixel area PA as an operation unit as shown in FIG. Of course, it may be a detector with one pixel.
- the present invention is not limited to the above-described configuration, and can be implemented in combination with various conventionally known modes without departing from the scope of the present invention.
- the method of determining the effective atomic number Z eff is not limited to the method described above.
- the effective atomic number Z eff is determined by comparing one of the atomic numbers Z High and Z Low estimated in step S 35 with the effective atomic number Z m designated in advance when performing the beam hardening correction. It is also good.
- the atomic number Z eff can be determined. Since only one of the low energy side effective atomic number (Z Low ) or the high energy side effective atomic number (Z High ) needs to be performed, it is advantageous for improving the processing speed.
- the effective atomic number Z eff is calculated . You may make it the structure decided. It is advantageous to improve the accuracy of substance identification.
- the number of effective atomic numbers estimated in the estimation step may be three or more, and the processing speed decreases as the number of estimated effective atomic numbers increases, but the accuracy of substance identification can be improved.
- the number of atomic numbers can be increased.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Toxicology (AREA)
- Mathematical Physics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
Abstract
Description
<本発明の目的>
[第1の実施形態]
<データ処理装置及びデータ処理方法>
<その1:事前処理>
a0=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a1=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a2=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a3=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a4=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
…(2)
で表される。つまり、係数aj(j=0~4)は原子番号Zの関数であり、係数Mj(j=0~4)は、係数ajに依存する量である。このフィッティング係数Mj,ajの一例を図12(A),(B),(C)にエネルギーBIN毎に示す。
上述のように事前処理が済むと、プロセッサ35はユーザとの間でいつでもインターラクティブに図4に示すように、X線透過データの収集及びその後の処理を行うことができる。
CoLow=CILow(e(-μLowt)
CoMiddle=CIMiddle(e(-μMiddlet)
CoHigh=CIHigh(e(-μHight)
…(3)
で表される。
[第2の実施形態]
<3次元散布図について>
(μLowt、μMiddlet、μHight)
を設定し、その長さ、即ち、線減弱値ベクトル長
((μLowt)2+(μMiddlet)2+(μHight)2)1/2
を分母とする正規化された3次元線源弱値ベクトル(以下、線減弱ベクトルと呼ぶ)は、線減弱ベクトルそのものから物質中のX線束の通過パスの長さ(厚みt)の成分が消えて
(μLow,μMiddle,μHigh)/(μLow 2+μMiddle 2+μHigh 2)1/2 …(5)
として演算される。互いに直交する3つの軸がμLowt、μMiddlet、μHightをそれぞれ表す3次元座標を設定すると、この3次元線減弱ベクトルの始点はその3次元座標の原点に位置し、終点が半径1の例えば球体表面に位置する。この3次元線減弱ベクトルを各画素について演算し、上記3次元座標にマッピングすると、それらの終点はかかる球体表面の所定の一点を中心に、その周辺の一定範囲に分布する、統計誤差を含む点の集合として散布点される。このような散布点を描いた3次元マッピング図を、本発明者等は3次元散布図と呼んでいる。この3次元散布図の例を図22に示す。同図において、参照符号Vrが3次元線減弱ベクトルを示し、参照符号DPが散布点を示す。
<吸収ベクトル長画像について>
t(μLow 2+μMiddle 2+μHigh 2)1/2 …(6)
で演算できるので、本発明者等は、このスカラー量を吸収ベクトル長(absorption vector length)(又は擬似的吸収値)と呼んでいる。この吸収ベクトル長を画素値として2次元画像を作成でき、本発明者等は、この2次元画像を吸収ベクトル長画像(又は擬似的な吸収画像)と呼んでいる。この吸収ベクトル長画像の例を図23に模式的に示す。
<平均吸収値画像について>
画素値=t(μLow+μMiddle+μHigh)/3 ‥‥(7)
又は
画素値=t(a1μLow+a2μMiddle+a3μHigh)/3 ‥‥(8)
ここで、
a1,a2,a3: 0以上の正の実数からなる重み付け係数であり、
a1+a2+a3=3とする
の式に基づいて演算できる。即ち、厚さtに依存したスカラー量としての画素値を得ることができる。ここで、分母に3を充てているのは、3つのエネルギーBIN:BinLow、BinMiddle、BinHigh、即ち全エネルギーBINに渡る平均値を演算するためである。
<変形例>
21 X線発生器
22 X線管
24 検出器
25 検出層
26 データ収集回路
30 データ処理装置(コンピュータ)
33 ROM
33A、33B,33C 記憶領域
35 プロセッサ(CPUを搭載)
37 入力器
38 表示器
P 画素
PA 画素領域
OB 対象物
Claims (30)
- 連続X線スペクトルのエネルギーを有するビーム状のX線を対象物に照射して、当該対象物を透過した前記X線を検出して、予め設定した2つ以上のX線エネルギーBINそれぞれにおいて少なくとも1画素からなる画像領域毎に当該X線の光子数を計数して、その計数値を示すデータを処理する方法において、
前記対象物が存在しない状態と存在する状態の双方における前記計数値の比で示される計数データを前記X線エネルギーBIN毎に且つ前記画素領域毎に演算する演算ステップと、
前記計数データに、予め指定される実効原子番号に応じた補正情報に基づき、前記X線が前記対象物を透過するときに受けるビームハードニング現象を補正するビームハードニング補正を、前記画素領域毎に且つ前記X線エネルギーBIN毎に施してX線減弱量(μt)(μは線減弱係数、tは前記対象物の前記X線の投影方向に沿った厚み)を求める補正ステップと、
2つ以上の前記X線エネルギーBINのうち選択される2つの前記X線エネルギーBINのX線減弱量を規格化して少なくとも1つの規格化X線減弱量を前記画素領域毎に求める規格化ステップと、
前記規格化X線減弱量と元素の実効原子番号との理論的な対応関係を示す参照情報から、少なくとも1つの実効原子番号を前記画素領域毎に推定する推定ステップと、
前記推定ステップで推定される少なくとも1つの実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち少なくとも2つの実効原子番号を比較して互いに一致しているか又は一致していると見做すことができる状態か否かを判定する一致判定ステップと、
を有することを特徴とする処理方法。 - 前記X線エネルギーBINが予め3つ以上設定されていて、
前記規格化ステップは、前記X線エネルギーBINのうち選択される2つの前記X線エネルギーBINの組み合わせが異なる状態で2つ以上の前記規格化X線減弱量を求める請求項1に記載の処理方法。 - 前記推定ステップは、前記規格化ステップで求められる2つ以上の前記規格化X線減弱量から2つ以上の実効原子番号を求めて、
前記一致判定ステップは、前記推定ステップで求められる2つ以上の実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち少なくとも2つの実効原子番号を比較する請求項2に記載の処理方法。 - 前記一致判定ステップは、前記推定ステップで推定される少なくとも1つの実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち選択される二つの値の差を演算する演算ステップと、この差が所定の閾値以下であるか否かを判定する差判定ステップとを含む請求項1~3のいずれかに記載の処理方法。
- 前記差判定ステップにより前記差が前記閾値以下であると判定されたときに、前記推定ステップで推定される少なくとも1つの実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち選択される二つが一致又は略一致していると見做して、その一致又は略一致している値が真の値又は真の値に近い実効原子番号であると前記X線エネルギーBIN毎に且つ前記画素領域毎に提示する提示ステップを有する請求項4に記載の処理方法。
- 前記画素領域は1つの画素から成る領域である請求項1~5のいずれかに記載の処理方法。
- 前記X線エネルギーBINは、エネルギースペクトル上で互いに隣接する低エネルギーBIN、中位エネルギーBIN、及び高エネルギーBINから成る3つのX線エネルギーBINであり、
前記低エネルギー側の2つのエネルギーBINは前記低エネルギーBIN及び前記中位エネルギーBINから成り、
前記高エネルギー側の2つのエネルギーBINは前記中位エネルギーBIN及び前記高エネルギーBINから成る請求項1~6の何れかに記載の処理方法。 - 前記X線エネルギーBINは、エネルギースペクトル上で互いに隣接する又はディスクリートな低エネルギーから高エネルギーまでの4つ以上のX線エネルギーBINであり、
前記低エネルギー側の2つのエネルギーBINは、前記低エネルギー側の互いに異なる2つのエネルギーBINから成り、
前記高エネルギー側の2つのエネルギーBINは、前記高エネルギー側の互いに異なる2つのエネルギーBINから成る、
ことを特徴とする請求項1~6の何れかに記載の処理方法。 - 連続X線スペクトルのエネルギーを有するビーム状のX線を対象物に照射し、当該対象物を透過した前記X線を検出し、予め設定した2つ以上のX線エネルギーBINそれぞれにおいて少なくとも1画素から成る画素領域毎に当該X線の光子数を計数し、その計数値を示すデータを処理する装置において、
前記対象物が存在しない状態と存在する状態の双方における前記計数値の比で示される計数データを、前記X線エネルギーBIN毎に且つ前記画素領域毎に演算する演算手段と、
前記計数データに、予め指定される実効原子番号に応じた補正情報に基づき前記X線が前記対象物を透過するときに受けるビームハードニング現象を補正するビームハードニング補正を、前記画素領域毎に且つ前記X線エネルギーBIN毎に施してX線減弱量(μt)(μは線減弱係数、tは前記対象物の前記X線の投影方向に沿った厚み)を求める補正手段と、
2つ以上の前記X線エネルギーBINのうち選択される2つの前記X線エネルギーBINのX線減弱量を規格化して少なくとも1つの規格化X線減弱量をを前記画素領域毎に求める規格化手段と、
前記X線減弱量と元素の実効原子番号との理論的な対応関係を示す参照情報から、少なくとも1つの実効原子番号を前記画素領域毎に推定する推定手段と、
前記推定ステップで推定される少なくとも1つの実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち少なくとも2つの実効原子番号を比較して互いに一致しているか又は一致していると見做すことができる状態か否かを判定する一致判定手段と、
を備えたことを特徴とする処理装置。 - 連続X線スペクトルのエネルギーを有するビーム状のX線を対象物に照射し、当該対象物を透過した前記X線を検出し、予め設定した2つ以上のX線エネルギーBINそれぞれにおいて少なくとも1画素から成る画素領域毎に当該X線の光子数を計数し、その計数値を示すデータを処理する方法において、
前記X線が前記対象物を透過するときに当該X線がビームハードニング現象を受けたX線計数値を補正するための補正情報を、原子番号が既知の複数の物質の質量厚さ(mass thickness:ρt)と前記各X線エネルギーBINの実効エネルギーにおけるX線減弱量(μt:μは線減弱係数、tは物質のX線パス方向の厚さ)との特性に基づき、前記X線エネルギーBIN毎に事前に準備する事前処理ステップと、
前記事前処理ステップにより準備されている前記補正情報を用いて前記補正を前記画素
領域毎に行って最終的に前記X線減弱量を確定させるとともに、この確定された前記X線減弱量を処理する減弱量処理ステップと、
を有することを特徴とする処理方法。 - 前記減弱量処理ステップは、前記補正ステップにより補正され且つ最終的に確定された前記X線減弱量に基づき光子計測によるX線画像を作成するX線画像作成ステップと、このX線画像を提示するX線画像提示ステップとを有する、ことを特徴とする請求項10に記載の処理方法。
- 前記事前処理ステップは、
前記対象物の組成を成す元素の実効原子番号の所望範囲(Zmin~Zmax)を設定するステップと、
前記所望範囲実効原子番号(Zm)の実効原子番号から、その下限値及び上限値値を含む複数の実効原子番号を離散的に選択し、当該複数の実効原子番号を持つ元素に対して、質量厚さ(mass thickness:ρt)を横軸に採り且つ当該各X線エネルギーBINの実効エネルギーにおける線減弱量(μt:μは線減弱係数、tは物質のX線パス方向の厚さ)を縦軸に採った2次元座標上で当該各実効原子番号のグラフを理論的に推定するステップと、
前記所望範囲(Zmin~Zmax)の実効原子番号から所望の実効原子番号(Zm=7)を指定するステップと、
前記指定された実効原子番号の元素から成る物質に単色X線が照射されたと仮定したときの前記2次元座標上の直線を目標関数として設定するステップと、
前記2次元座標上で、前記横軸方向に前記目標関数の傾き(μ/ρ)を乗算して、前記複数の実効原子番号それぞれの複数の曲線を当該実効原子番号の変数であるとして一般化するステップと、
前記一般化された複数の曲線の中から、前記指定された実効原子番号の元素の曲線を指定し、この指定曲線と他の曲線との残差に基づく、前記ビームハードニングを補正するためのビームハードニング補正関数を前記補正情報として、記憶部に事前に記憶するステップと、
を有する、ことを特徴とする請求項10に記載の処理方法。 - 前記補正ステップは、
前記補正されたビームハードニング補正関数のデータを前記記憶部から読み出し当該計数値の前記補正を行うことを特徴とする請求項12に記載の処理方法。 - 前記減弱量処理ステップは、
前記対象物が存在しない状態と存在する状態の双方における前記計数値の比で示される計数データを、前記X線エネルギーBIN毎に且つ前記画素領域毎に演算する演算ステップと、
前記記憶部から所望の実効原子番号に対応した前記ビームハードニング補正関数を前記補正情報として読出し、この補正情報に基づき前記計数データに前記ビームハードニング補正を、前記画素領域毎に且つ前記X線エネルギーBIN毎に施してX線減弱量(μt)(μは線減弱係数、tは前記対象物の前記X線の投影方向に沿った厚み)を求める第1の補正ステップと、
2つ以上の前記X線エネルギーBINのうち選択される2つの前記X線エネルギーBINのX線減弱量を規格化して少なくとも1つの規格化X線減弱量を前記画素領域毎に求める規格化ステップと、
前記規格化されたX線減弱量と元素の実効原子番号との理論的な対応関係を示す参照情報から、少なくとも1つの実効原子番号を前記画素領域毎に推定する推定ステップと、
前記推定ステップで推定される少なくとも1つの実効原子番号と前記第1の補正ステップで予め指定される実効原子番号とのうち少なくとも2つの実効原子番号を相互に比較して当該実効原子番号の一致度が許容範囲内か否かを判定する一致判定ステップと、
前記一致の判定により前記一致度が許容範囲に入っているときには当該一致度を呈する当該実効原子番号を前記画素領域における真の実効原子番号(Zeff)であると推定する推定ステップと、
を有する請求項12又は13に記載の処理方法。 - 前記減弱量処理ステップは、
前記推定ステップで推定される少なくとも1つの実効原子番号と前記第1の補正ステップを行う際に予め指定される実効原子番号とのうち選択される二つが互いに一致した又は一致したと見做せるときの前記補正情報を取得する補正情報取得ステップと、
前記第1の補正ステップにおいて前記補正情報で補正された前記X線減弱量(μt)を前記画素領域毎の最終的なX線減弱量として確定する減弱量確定ステップと、
前記減弱量確定ステップにより確定された前記画素領域の最終的なX線減弱量に基づき前記対象物の光子計測に拠るX線画像を提示する画像提示ステップと、
を有する請求項14に記載の処理方法。 - 前記減弱量処理ステップは、
前記推定ステップにより推定される前記実効原子番号から実効原子番号画像を生成するステップを有する、請求項14に記載の処理方法。 - 前記減弱量処理ステップは、
前記一致の判定により前記一致度が許容範囲に入っていないと判断されときに、前記記憶部から別の実効原子番号に対応した別のビームハードニング補正関数を前記補正情報として読出し、この補正情報に基づき前記計数データに前記ビームハードニング補正を、前記画素領域毎に且つ前記X線エネルギーBIN毎に施してX線減弱量を求める第2の補正ステップと、
この第2の補正ステップを行った後、前記規格化ステップ、前記推定ステップ、及び、前記一致判定ステップを1回以上、繰り返し指令する繰り返し指令ステップと、を有する、
ことを特徴とする請求項13~16の何れか一項に記載の処理方法。 - 前記第2の補正ステップは、前記別のビームハードニング補正関数を読み出すときに、その前に読み出した実効原子番号とその隣のディスクリートな実効原子番号と間の選択値に応じて前記残差を比例配分した補正情報を演算する、ことを特徴とする請求項15に記載の処理方法。
- 前記残差に基づく前記ビームハードニング補正関数は4次関数で表され、その4次関数は、係数a0、a1、a2、a3、a4としたときに、
f(ρt)=a0+a1×(ρt)+a2×(ρt)2+a3×(ρt)3+a4×(ρt)4
ここで、Zを実効原子番号とし、係数Mj(j=0~4)としたとき、係数aj(j=0~4)は、
a0=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a1=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a2=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a3=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
a4=M0+M1×Z+M2×Z2+M3×Z3+M4×Z4
で表されることを特徴とする請求項12~18の何れか一項に記載の処理方法。 - 前記残差に基づく前記ビームハードニング補正関数は、前記実効原子番号Zで高次関数で表される係数aを持つ、前記質量厚さ(mass thickness:(t)の高次関数で表されることを特徴とする請求項12~18の何れか一項に記載の処理方法。
- 前記減弱量処理ステップは、
前記計数値に基づく画像をモニタに表示させる表示ステップと、
前記モニタ上で関心領域(ROI)を表示させて当該関心領域を構成する画素それぞれを前記画素領域として設定する画素設定ステップと、
を有することを特徴とする請求項11~20の何れか一項に記載の処理方法。 - 前記減弱量処理ステップは、
前記計数値に基づく画像をモニタに表示させる表示ステップと、
前記モニタ上で関心領域(ROI)を表示させて当該関心領域を構成する画素のうち、互いに隣接する複数の画素を束ね、その束ねた複数の画素から成る各領域を前記画素領域として設定する画素設定ステップと、
を有することを特徴とする請求項11~20の何れか一項に記載の処理方法。 - 前記一致判定ステップは、前記推定ステップで求められる1つ以上の実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち選択される二つの値の差を演算する演算ステップと、この差が所定の閾値以下であるか否かを判定する差判定ステップとを含む請求項14~22のいずれかに記載の処理方法。
- 前記差判定ステップにより前記差が前記閾値以下であると判定されたときに、前記推定ステップで求められる1つ以上の実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち選択される二つの値が一致又は略一致していると見做して、その一致又は略一致している値が真の値又は真の値に近い実効原子番号であると前記X線エネルギーBIN毎に且つ前記画素領域毎に提示する提示ステップを有する、ことを特徴とする請求項23に記載の処理方法。
- 2つ以上の前記X線エネルギーBINは、エネルギースペクトル上で互いに隣接する低エネルギーBIN、中位エネルギーBIN、及び高エネルギーBINから成る3つのX線エネルギーBINであり、
前記低エネルギー側の2つのエネルギーBINは前記低エネルギーBIN及び前記中位エネルギーBINから成り、
前記高エネルギー側の2つのエネルギーBINは前記中位エネルギーBIN及び前記高エネルギーBINから成る、
ことを特徴とする請求項11~24何れか一項に記載の処理方法。 - 2つ以上の前記X線エネルギーBINは、エネルギースペクトル上で互いに隣接する又はディスクリートな低エネルギーから高エネルギーまでの4つ以上のX線エネルギーBINであり、
前記低エネルギー側の2つのエネルギーBINは、前記低エネルギー側の互いに異なる2つのエネルギーBINから成り、
前記高エネルギー側の2つのエネルギーBINは、前記高エネルギー側の互いに異なる2つのエネルギーBINから成る、
ことを特徴とする請求項11~24の何れか一項に記載の処理方法。 - 連続X線スペクトルのエネルギーを有するビーム状のX線を対象物に照射し、当該対象物を透過した前記X線を検出し、予め設定した2つ以上のX線エネルギーBINそれぞれにおいて少なくとも1画素から成る画素領域毎に当該X線の光子数を計数し、その計数値を処理するX線装置において、
前記X線が前記対象物を透過するときに当該X線がビームハードニング現象を受けたX線計数値を補正するための補正情報を、原子番号が既知の複数の物質の質量厚さ(mass thickness:ρt)と前記各X線エネルギーBINの実効エネルギーにおけるX線減弱量(μt:μは線減弱係数、tは物質のX線パス方向の厚さ)との特性に基づき、前記X線エネルギーBIN毎に事前に準備する事前準備手段と、
前記事前準備手段により準備されている前記補正情報を用いて前記補正を前記画素領域毎に行って最終的に前記X線減弱量を確定させるとともに、この確定された前記X線減弱量を処理する減弱量処理手段と、
を備えたことを特徴とするX線装置。 - 前記減弱量処理手段は、前記補正手段により補正され且つ最終的に確定された前記X線減弱量に基づき光子計測によるX線画像を作成するX線画像作成手段と、このX線画像を提示するX線画像提示手段とを有する、ことを特徴とする請求項27に記載のX線装置。
- 前記X線装置は、前記X線を光子計数法で検出する構成を有するX線医療診断機器又はX線非破壊検査装置であることを特徴とする請求項27又は28に記載のX線装置。
- 連続X線スペクトルのエネルギーを有するビーム状のX線を対象物に照射し、当該対象物を透過した前記X線を検出し、予め設定した2つ以上のX線エネルギーBINそれぞれにおいて少なくとも1画素から成る画素領域毎に当該X線の光子数を計数し、その計数値を処理する方法に係るプログラムをコンピュータにより読出し可能に記憶する記録媒体であって、そのプログラムを実行することにより当該コンピュータは、
前記対象物が存在しない状態と存在する状態の双方における前記計数値の比で示される計数データを、前記X線エネルギーBIN毎に且つ前記画素領域毎に演算する演算ステップと、
前記計数データに、予め指定される実効原子番号に応じた補正情報に基づき、前記X線が前記対象物を透過するときに受けるビームハードニング現象を補正するビームハードニング補正を、前記画素領域毎に且つ前記X線エネルギーBIN毎に施してX線減弱量(μt)(μは線減弱係数、tは前記対象物の前記X線の投影方向に沿った厚み)を求める補正ステップと、
2つ以上の前記X線エネルギーBINのうち選択される2つの前記X線エネルギーBINのX線減弱量を規格化して少なくとも1つの規格化X線減弱量を前記画素領域毎に求める規格化ステップと、
前記規格化されたX線減弱量と元素の実効原子番号との理論的な対応関係を示す参照情報から、少なくとも1つの実効原子番号を前記画素領域毎に推定する推定ステップと、 前記推定ステップで推定される少なくとも1つの実効原子番号と前記補正ステップで予め指定される実効原子番号とのうち少なくとも2つの実効原子番号を比較して互いに一致しているか又は一致していると見做せる状態か否かを判定する一致判定ステップと、
を実行可能であることを特徴とするコンピュータ読出し可能な非遷移的実態的記録媒体。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP18870980.2A EP3702812A4 (en) | 2017-10-26 | 2018-10-26 | METHOD AND DEVICE FOR PROCESSING PHOTONCOUNTING X-RAY DETECTION DATA AND X-RAY DEVICE |
CN201880039187.2A CN110770607B (zh) | 2017-10-26 | 2018-10-26 | 处理光子计数型x射线检测数据的方法及装置、以及x射线装置 |
US16/627,451 US11099141B2 (en) | 2017-10-26 | 2018-10-26 | Method and apparatus for processing photon counting-type X-ray detection data and X-ray apparatus |
KR1020197036226A KR102330146B1 (ko) | 2017-10-26 | 2018-10-26 | 광자 계수형의 x선 검출 데이터를 처리하는 방법 및 장치, 및 x선 장치 |
JP2019550321A JP7197174B2 (ja) | 2017-10-26 | 2018-10-26 | 光子計数型のx線検出データを処理する方法及び装置、並びに、x線装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017207458 | 2017-10-26 | ||
JP2017-207458 | 2017-10-26 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2019083014A1 true WO2019083014A1 (ja) | 2019-05-02 |
Family
ID=66247932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2018/039884 WO2019083014A1 (ja) | 2017-10-26 | 2018-10-26 | 光子計数型のx線検出データを処理する方法及び装置、並びに、x線装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US11099141B2 (ja) |
EP (1) | EP3702812A4 (ja) |
JP (1) | JP7197174B2 (ja) |
KR (1) | KR102330146B1 (ja) |
CN (1) | CN110770607B (ja) |
WO (1) | WO2019083014A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021210442A1 (ja) * | 2020-04-16 | 2021-10-21 | 浜松ホトニクス株式会社 | 放射線検査方法、放射線検査装置、放射線検査システム、及び放射線検査プログラム |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3508887A1 (en) * | 2018-01-09 | 2019-07-10 | Koninklijke Philips N.V. | Charge sharing calibration method and system |
US11977037B2 (en) | 2018-10-22 | 2024-05-07 | Rapiscan Holdings, Inc. | Insert for screening tray |
JP7271209B2 (ja) * | 2019-02-06 | 2023-05-11 | キヤノン株式会社 | 放射線撮像装置、放射線撮像システムおよび放射線撮像装置の制御方法 |
DE102019111567A1 (de) * | 2019-05-03 | 2020-11-05 | Wipotec Gmbh | Verfahren und Vorrichtung zur Röntgeninspektion von Produkten, insbesondere von Lebensmitteln |
JP7440271B2 (ja) * | 2020-01-10 | 2024-02-28 | 富士フイルムヘルスケア株式会社 | 放射線撮像装置および光子計数型検出器の較正方法 |
US12019035B2 (en) * | 2021-07-16 | 2024-06-25 | Rapiscan Holdings, Inc. | Material detection in x-ray security screening |
CN114063138B (zh) * | 2021-11-16 | 2023-07-25 | 武汉联影生命科学仪器有限公司 | 扫描成像系统有效能量的测定方法、设备和扫描成像系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030072409A1 (en) * | 2001-10-12 | 2003-04-17 | Kaufhold John Patrick | Methods and apparatus for estimating a material composition of an imaged object |
JP2006101926A (ja) | 2004-09-30 | 2006-04-20 | M & C:Kk | 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法 |
JP2010091483A (ja) | 2008-10-09 | 2010-04-22 | Anritsu Sanki System Co Ltd | 異物検出方法および装置 |
WO2010061810A1 (ja) | 2008-11-27 | 2010-06-03 | 株式会社 日立メディコ | 放射線撮像装置 |
WO2015111728A1 (ja) | 2014-01-23 | 2015-07-30 | 株式会社ジョブ | X線検査装置及びx線検査方法 |
WO2016171186A1 (ja) | 2015-04-20 | 2016-10-27 | 株式会社ジョブ | X線検査用のデータ処理装置及びデータ処理方法、並びに、その装置を搭載したx線検査装置 |
WO2017069286A1 (ja) | 2015-10-23 | 2017-04-27 | 株式会社ジョブ | X線装置、データ処理装置及びデータ処理方法 |
JP2017127638A (ja) * | 2016-01-18 | 2017-07-27 | 東芝メディカルシステムズ株式会社 | X線ct装置、情報処理装置、および情報処理方法 |
JP2017528735A (ja) * | 2014-08-16 | 2017-09-28 | エフ・イ−・アイ・カンパニー | 材料特性評価のための断層撮影再構成 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009245583A (ja) | 2003-06-18 | 2009-10-22 | Konica Minolta Opto Inc | 光ピックアップ装置用の集光素子 |
CN100365664C (zh) * | 2006-03-02 | 2008-01-30 | 西安交通大学 | 基于hl一致性条件的ct投影数据射束硬化效应校正方法 |
US7298812B2 (en) * | 2006-03-31 | 2007-11-20 | General Electric Company | Image-based material decomposition |
WO2011036436A1 (en) * | 2009-09-22 | 2011-03-31 | Isis Innovation Limited | X-ray imaging |
CN103900931B (zh) * | 2012-12-26 | 2017-04-26 | 首都师范大学 | 一种多能谱ct成像方法及成像系统 |
-
2018
- 2018-10-26 KR KR1020197036226A patent/KR102330146B1/ko active IP Right Grant
- 2018-10-26 US US16/627,451 patent/US11099141B2/en active Active
- 2018-10-26 WO PCT/JP2018/039884 patent/WO2019083014A1/ja unknown
- 2018-10-26 JP JP2019550321A patent/JP7197174B2/ja active Active
- 2018-10-26 EP EP18870980.2A patent/EP3702812A4/en active Pending
- 2018-10-26 CN CN201880039187.2A patent/CN110770607B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030072409A1 (en) * | 2001-10-12 | 2003-04-17 | Kaufhold John Patrick | Methods and apparatus for estimating a material composition of an imaged object |
JP2006101926A (ja) | 2004-09-30 | 2006-04-20 | M & C:Kk | 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法 |
JP2010091483A (ja) | 2008-10-09 | 2010-04-22 | Anritsu Sanki System Co Ltd | 異物検出方法および装置 |
WO2010061810A1 (ja) | 2008-11-27 | 2010-06-03 | 株式会社 日立メディコ | 放射線撮像装置 |
WO2015111728A1 (ja) | 2014-01-23 | 2015-07-30 | 株式会社ジョブ | X線検査装置及びx線検査方法 |
JP2017528735A (ja) * | 2014-08-16 | 2017-09-28 | エフ・イ−・アイ・カンパニー | 材料特性評価のための断層撮影再構成 |
WO2016171186A1 (ja) | 2015-04-20 | 2016-10-27 | 株式会社ジョブ | X線検査用のデータ処理装置及びデータ処理方法、並びに、その装置を搭載したx線検査装置 |
WO2017069286A1 (ja) | 2015-10-23 | 2017-04-27 | 株式会社ジョブ | X線装置、データ処理装置及びデータ処理方法 |
JP2017127638A (ja) * | 2016-01-18 | 2017-07-27 | 東芝メディカルシステムズ株式会社 | X線ct装置、情報処理装置、および情報処理方法 |
Non-Patent Citations (2)
Title |
---|
KIMOTO, NATSUMI ET AL.: "Precise material identification method based on a photon counting technique with correction of the beam hardening effect in X-ray spectra", APPLIED RADIATION AND ISOTOPES, vol. 124, 9 March 2017 (2017-03-09), pages 16 - 26, XP029979981, DOI: doi:10.1016/j.apradiso.2017.01.049 * |
See also references of EP3702812A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021210442A1 (ja) * | 2020-04-16 | 2021-10-21 | 浜松ホトニクス株式会社 | 放射線検査方法、放射線検査装置、放射線検査システム、及び放射線検査プログラム |
Also Published As
Publication number | Publication date |
---|---|
JPWO2019083014A1 (ja) | 2020-12-24 |
CN110770607A (zh) | 2020-02-07 |
KR20200004392A (ko) | 2020-01-13 |
JP7197174B2 (ja) | 2022-12-27 |
EP3702812A1 (en) | 2020-09-02 |
US20200249179A1 (en) | 2020-08-06 |
CN110770607B (zh) | 2024-02-13 |
EP3702812A4 (en) | 2021-07-21 |
US11099141B2 (en) | 2021-08-24 |
KR102330146B1 (ko) | 2021-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102330146B1 (ko) | 광자 계수형의 x선 검출 데이터를 처리하는 방법 및 장치, 및 x선 장치 | |
JP7217020B2 (ja) | X線装置、x線検査方法、及びデータ処理装置 | |
JP6590381B2 (ja) | X線装置、データ処理装置及びデータ処理方法 | |
US7643866B2 (en) | Method for producing a computed tomography display of tissue structures by applying a contrast medium | |
US8150131B2 (en) | Computer tomography method for determining an object region with the aid of projection images | |
EP2243021B1 (en) | System and method for quantitative imaging of chemical composition to decompose multiple materials | |
US20190180482A1 (en) | Image reconstruction device, x-ray ct device, and image reconstruction method | |
US20140133719A1 (en) | System and method for multi-material correction of image data | |
US20170231589A1 (en) | X-ray ct apparatus, upsampling method of projection data, and image reconstruction method | |
EP3977935A1 (en) | Data processing device and data processing method for processing x-ray detection data and x-ray inspection apparatus provided with device and method | |
US20080226019A1 (en) | Multiple Scatter Correction | |
EP3367086B1 (en) | X-ray device, data processing device, and data processing method | |
JP7348642B2 (ja) | データ処理装置、そのデータ処理装置を搭載したx線装置、及びデータ処理方法 | |
JP5920770B2 (ja) | X線コンピュータ断層撮影方法及び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: 18870980 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 20197036226 Country of ref document: KR Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 2019550321 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2018870980 Country of ref document: EP Effective date: 20200526 |