WO2019117147A1 - X線装置および構造物の製造方法 - Google Patents

X線装置および構造物の製造方法 Download PDF

Info

Publication number
WO2019117147A1
WO2019117147A1 PCT/JP2018/045516 JP2018045516W WO2019117147A1 WO 2019117147 A1 WO2019117147 A1 WO 2019117147A1 JP 2018045516 W JP2018045516 W JP 2018045516W WO 2019117147 A1 WO2019117147 A1 WO 2019117147A1
Authority
WO
WIPO (PCT)
Prior art keywords
ray
rays
scattering
scattered
energy
Prior art date
Application number
PCT/JP2018/045516
Other languages
English (en)
French (fr)
Inventor
一明 鈴木
Original Assignee
株式会社ニコン
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2017239035A external-priority patent/JP2021028575A/ja
Priority claimed from JP2017245413A external-priority patent/JP2021028576A/ja
Application filed by 株式会社ニコン filed Critical 株式会社ニコン
Publication of WO2019117147A1 publication Critical patent/WO2019117147A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/02Investigating 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/04Investigating 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
    • G01N23/046Investigating 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 using tomography, e.g. computed tomography [CT]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/20Investigating 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 using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G01N23/201Investigating 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 using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials by measuring small-angle scattering

Definitions

  • the present invention relates to a method of manufacturing an x-ray device and structure.
  • a Compton camera which detects Compton scattering of incident radiation (X-ray or ⁇ -ray) and generates an image of a radiation source based on the result (for example, Patent Document 1).
  • X-ray or ⁇ -ray Compton scattering of incident radiation
  • Patent Document 1 Patent Document 1
  • the nondestructive inspection device which detects Compton scattering and generates an image of a to-be-measured object based on the result.
  • X-rays including the transmitted X-rays transmitted through the object to be measured or the scattered X-rays scattered inside the object to be measured are incident, and the incident X-rays are scattered
  • a first detector for detecting a scattering position a second detector for detecting an incident position of the incident X-ray upon which the X-ray scattered or transmitted by the first detector is incident; the scattering position;
  • the X-ray apparatus has a first energy which is lost by incidence and scattering of X-rays including transmitted X-rays transmitted through the object to be measured or scattered X-rays scattered inside the object to be measured A second detector for detecting the second energy of the X-rays scattered or transmitted by the first detector, and the transmission based on the first energy and the second energy. And a specific unit that specifies at least one of an X-ray and the scattered X-ray.
  • the X-ray apparatus is configured to transmit transmitted X-rays transmitted through the object to be measured or electrons generated by scattering after incidence of X-rays including scattered X-rays scattered inside the object to be measured
  • the transmitted X based on a detector for detecting information on the progression and information on the progression of the X-rays after the scattering, information on the progression of the electrons, and information on the progression of the X-rays
  • an identifying unit that identifies at least one of a ray and the scattered X-ray.
  • an X-ray apparatus includes: a detector for detecting X-rays including transmitted X-rays transmitted through a device under test or scattered X-rays scattered inside the device under test; An identifying unit configured to estimate an incident direction of the X-ray and to identify at least one of the transmitted X-ray and the scattered X-ray based on the estimated incident direction and the position of the X-ray source.
  • a method of manufacturing a structure creates design information related to the shape of the structure; creating the structure based on the design information; The X-ray apparatus according to any one of the fourth to the fourth aspects is used for measurement to acquire shape information, and the acquired shape information is compared with the design information.
  • FIG. 1 is a block diagram schematically showing a configuration of an X-ray apparatus according to a first embodiment. It is a figure which shows typically the path
  • FIG. 16 is a view schematically showing a positional relationship among an X-ray source, a detector, and an object to be measured, for explaining an example in the case where the maximum energy of X-rays is unknown in the third embodiment.
  • FIG. 13 is a view schematically showing a path of X-rays incident on the scattering portion and the absorption portion of the detector for explaining an example in the case where the maximum energy of the X-rays is known in the third embodiment.
  • FIG. 18 is a view schematically showing a positional relationship among an X-ray source, a detector, and an object to be measured, for explaining an example in the case where the maximum energy of X-rays is known in the third embodiment.
  • the X-ray apparatus destroys the internal information (for example, the internal structure) of the object to be measured by irradiating the object to be measured with X-rays and detecting X-rays transmitted through the object to be measured. Get without.
  • An X-ray apparatus intended for industrial parts such as mechanical parts and electronic parts is called an industrial X-ray inspection apparatus.
  • FIG. 1 is a diagram showing an example of the configuration of an X-ray apparatus 100 according to the present embodiment.
  • the X-ray apparatus 100 includes a housing 1, an X-ray source 2, a placement unit 3, a detector 4, and a control device 5.
  • the housing 1 is disposed such that the lower surface thereof is substantially parallel (horizontal) to a floor surface of a factory or the like. Inside the housing 1, an X-ray source 2, a placement unit 3, and a detector 4 are accommodated.
  • the housing 1 includes an X-ray shielding material to prevent X-rays from leaking to the outside of the housing 1.
  • lead is included as an X-ray shielding material.
  • the X-ray source 2 emits X-rays toward the Z-axis + direction along the optical axis Zr parallel to the Z-axis with the emission point P shown in FIG. 1 as the apex according to the control by the controller 5.
  • the emission point P coincides with the focal position of the electron beam propagating inside the X-ray source 2 described later. That is, the optical axis Zr is an axis connecting the emission point P which is the focal position of the electron beam of the X-ray source 2 and the center of the imaging region of the detector 4 described later.
  • the X-ray emitted from the X-ray source 2 may be any of a conically-spreading X-ray (so-called cone beam), a fan-like X-ray (so-called fan beam), and a linear X-ray (any pencil beam) May be.
  • cone beam conically-spreading X-ray
  • fan beam fan-like X-ray
  • linear X-ray any pencil beam
  • movement which moves a beam and the to-be-measured object S relatively.
  • the X-ray source 2 emits X-ray photons one by one when viewed on a short time scale.
  • the X-ray source 2 emits at least one of, for example, about 50 eV ultrasoft X-rays, about 0.1 to 2 keV soft X-rays, about 2 to 20 keV X-rays and about 20 to several MeV hard X-rays. . It is in the region of hard X-rays that the scattering is remarkable.
  • the mounting unit 3 includes a mounting table 31 on which the object to be measured S is mounted, and a manipulator unit 36 including a rotation driving unit 32, an X-axis moving unit 33, a Y-axis moving unit 34 and a Z-axis moving unit 35.
  • X-axis source 2 is provided on the Z-axis + side.
  • the mounting table 31 is rotatably provided by the rotation drive unit 32. As described later, when the rotation axis Yr by the rotation drive unit 32 moves in the X axis, Y axis, and Z axis directions, the mounting table 31 moves together.
  • the rotation drive unit 32 is configured by, for example, an electric motor, and rotates the mounting table 31 by the rotational force generated by the electric motor controlled and driven by the control device 5 described later.
  • the rotation axis Yr of the mounting table 31 is parallel to the Y axis and passes through the center of the mounting table 31.
  • the X-axis moving unit 33, the Y-axis moving unit 34, and the Z-axis moving unit 35 are controlled by the control device 5 to move the mounting table 31 in the X-axis direction, the Y-axis direction, and the Z-axis direction.
  • the Z-axis moving unit 35 is controlled by the control device 5 so that the distance from the X-ray source 2 to the object S is a distance according to the magnification of the object S in the image to be captured.
  • the table 31 is moved in the Z-axis direction.
  • the detector 4 shown in FIG. 1 is provided closer to the Z-axis + side than the mounting table 31. That is, the mounting table 31 is provided between the X-ray source 2 and the detector 4 in the Z-axis direction.
  • the detector 4 includes a scattering portion 41 and an absorption portion 42, and the transmitted X-rays transmitted through the object S (or the transmitted X-rays which go straight without being irradiated to the object S) and the scattered X-rays Detect at least one.
  • the scattering portion 41 is provided on the Z direction ⁇ side, that is, on the side closer to the object S than the absorbing portion 42.
  • the thickness of the scattering portion 41 in the Z direction is smaller than the thickness of the absorbing portion 42 in the Z direction.
  • the scattering unit 41 has an incident surface 411 parallel to the XY plane. X-rays emitted from the X-ray source 2 and transmitted through the object to be measured S mounted on the mounting table 31 or scattered in the object to be measured S enter the incident surface 411. In the following description, X-rays transmitted through the inside of the object S without being scattered and entering the detector 4 without being scattered are transmitted X-rays and scattered inside the object S and detected. X-rays incident on are called scattered X-rays.
  • the scattering unit 41 includes, for example, a scintillator unit including a known scintillation material, and a light receiving unit that amplifies and receives light emitted from the scintillator unit.
  • the X-rays incident on the incident surface 411 of the scattering unit 41 are Compton scattered at a certain probability in the scintillator unit.
  • fluorescence is emitted from the scintillator unit by the recoil electrons.
  • the fluorescence is amplified by the photomultiplier in the light receiving unit and output to the control device 5 as an electric signal.
  • the X-rays Compton-scattered in the scintillator portion change in the traveling direction, and then exit the scattering portion 41 and enter the absorbing portion 42 provided on the Z direction + side.
  • X-rays that have not been Compton scattered in the scintillator part exit the scattering part 41 without changing the traveling direction in the scattering part 41 and enter the absorbing part 42.
  • the absorbing portion 42 has an incident surface 421 parallel to the XY plane.
  • the X-rays arriving from the scattering section 41 are incident on the incident surface 421.
  • the absorbing unit 42 is configured of a scintillator unit including a known scintillation material, a light receiving unit that amplifies and receives light emitted from the scintillator unit, and the like.
  • the X-rays incident on the incident surface 421 of the scintillator unit are absorbed by the scintillator unit of the absorber unit 42 and emit fluorescence, and the light energy of the emitted fluorescence is amplified by the photomultiplier in the above-mentioned light receiving unit to become electrical energy It is converted and output to the control device 5 as an electric signal.
  • the scattering unit 41 and the absorption unit 42 of the detector 4 both have a structure in which a plurality of pixels of the scintillator unit and the light receiving unit are two-dimensionally arranged, and each of the plurality of pixels of the scintillator unit and the light receiving unit Are arranged to correspond to each other.
  • the control device 5 has a microprocessor and its peripheral circuits, etc., and reads and executes a control program stored in advance in a storage medium (for example, a flash memory etc.) (not shown). Control each part.
  • the control device 5 controls the operation of the X-ray source 2, the mounting table control unit 52 which controls the drive operation of the manipulator unit 36, and the object under test based on the electrical signal output from the detector 4. It has the image generation part 53 which produces
  • the identifying unit 54 determines whether the X-rays from the X-ray source 2 are scattered within the object S based on the electrical signals output from the scattering unit 41 and the absorbing unit 42.
  • the identifying unit 54 identifies, among X-rays from the object S, scattered X-rays scattered at the object S or identifies X-rays other than the scattered X-rays, based on the determination result. That is, the identifying unit 54 identifies whether the X-rays incident on the detector 4 are X-rays scattered by the object S.
  • the image generation unit 53 generates intensity distribution data from X-rays other than the scattered X-rays identified by the identification unit 54, that is, transmitted X-rays that are not scattered in the object S or X-rays that are partially absorbed. It selects as internal information of the to-be-measured object S at the time of producing
  • the image reconstruction unit 56 projects the X-ray irradiation direction relative to the object S relative to the plurality of X-ray intensity distribution data obtained by changing the X-ray irradiation direction relative to the known image reconstruction processing method. By using this, a reconstructed image of the object S is generated.
  • Cross-sectional image data and three-dimensional data, which are the internal structure (cross-sectional structure) of the object S, are generated by the image reconstruction process.
  • the cross-sectional image data includes structural data of the object S in a plane parallel to the XZ plane.
  • Image reconstruction processing includes back projection, filter-corrected back projection, successive approximation, and the like.
  • FIG. 2A shows a case where an electrical signal is output at a certain point in both the pixel 41 a which is one pixel of the scattering portion 41 and the pixel 42 a which is one pixel of the absorbing portion 42.
  • the scattering unit 41 detects the scattering position where the incident X-ray is scattered in the detector 4 according to the position of the pixel 41a
  • the absorbing unit 42 detects the incident position where the X-ray is incident according to the position of the pixel 42a.
  • the identifying unit 54 calculates the energy E1 and the energy E2, respectively, based on the electrical signal output from the pixel 41a and the electrical signal output from the pixel 42a.
  • the energy E1 corresponds to the energy lost due to the Compton scattering of X-rays in the pixel 41a
  • the energy E2 corresponds to the energy of the X-rays reaching the pixel 42a after being Compton scattered in the pixel 41a.
  • the identifying unit 54 can calculate the Compton scattering angle ⁇ 1 according to the following equation (1), assuming that the angle of change in the X-ray traveling direction due to Compton scattering in the pixel 41a, ie, the Compton scattering angle is ⁇ 1.
  • the path L2 of the X-ray emitted from the pixel 41a and incident on the pixel 42a Assuming that the path of the X-ray incident on the pixel 41a is a path L1, the path L1 and the path L2 intersect at a scattering point T1 to form a Compton scattering angle ⁇ 1. Accordingly, the path L1 corresponds to any one of the straight lines toward the scattering point T1 on the conical surface CS with an angle of ⁇ 1 between the central axis L2 and the generatrix with the scattering point T1 as the apex and the path L2 as the central axis. I understand that.
  • the identifying unit 54 determines that the X-ray emitted from the X-ray source 2 is incident on the pixel 41 a without changing the traveling direction inside the object S, and the traveling path is L1. it can. Even when the X-ray source 2 does not exist on the conical surface CS, when the distance between the X-ray source 2 and the conical surface CS is smaller than a predetermined length, the identifying unit 54 similarly to the above.
  • the specifying unit 54 specifies that the X-ray incident on the detector 4 is the transmitted X-ray transmitted without being scattered in the inside of the measured object S.
  • the predetermined length described above is preferably set as appropriate based on the size of the pixels 41a and 42a and the arrangement of the X-ray source 2, the object S, the detector 4 and the like in the X-ray apparatus 100. .
  • FIG. 2C shows the case where the X-ray source 2 does not exist on the conical surface CS and the distance between the X-ray source 2 and the conical surface CS is larger than a predetermined length.
  • the X-ray source 2 does not exist on L1, which is one of the straight lines on the conical surface CS toward the scattering point T1.
  • the identifying unit 54 can determine that the X-ray emitted from the X-ray source 2 is scattered inside the object S, changes the traveling direction, and enters the pixel 41a.
  • the X-ray emitted from the X-ray source 2 and traveling along the path L0 enters the object S and is scattered at the scattering point T3 inside the object to change the traveling direction to L1. It is identified as the scattered X-ray incident on the pixel 41a.
  • the identifying unit 54 determines that the X-ray is not scattered by the object S Determine and identify X-rays other than scattered X-rays.
  • the identifying unit 54 determines that the X-rays are scattered by the object S, and the scattered X-rays Identify what is.
  • the detection accuracy ⁇ 1 of the Compton scattering angle ⁇ 1 is determined as follows. That is, it is assumed that the energy ET of the X-ray incident on the scattering unit 41 is constant. When the energy ET is the sum of the energy E1 and the energy E2, the detection accuracy ⁇ 1 satisfies the relationship of the following equation (2) obtained by differentiating the equation (1).
  • the detection accuracy ⁇ 1 is calculated by the following equation (3).
  • the energy of the X-ray emitted from the X-ray source 2 and incident on the scattering section 41 is 511 keV
  • the energy resolution is 3 percent
  • the Compton scattering angle ⁇ 1 is 20 °
  • the detection accuracy ⁇ 1 is 5.3 °
  • the detection accuracy ⁇ 1 is 1.8 °.
  • the resolution of the scattering angle ⁇ 1 is 10 mrad.
  • the identifying unit 54 is used to calculate the energy E1 obtained by the scattering unit 41 of the detector 4, the position (XSi, YSi, ZSi) of the pixel that has output the electric signal used for calculating the energy E1, and the energy E1.
  • Each information of time Ti when the electrical signal is output from the scattering unit 41 is acquired.
  • the position (XSi, YSi, ZSi) of this pixel corresponds to the scattering point T1 in each of FIGS. 2 (a), 2 (b), and 2 (c).
  • the time Ti corresponds to the time when X-rays enter the scattering unit 41 and are Compton scattered.
  • the identifying unit 54 is used to calculate the energy E2 obtained by the absorbing unit 42 of the detector 4, the position (XAj, YAj, ZAj) of the pixel that has output the electric signal used for calculating the energy E2, and the energy E2.
  • Each information of time Tj when the electrical signal is output from the absorption unit 42 is acquired.
  • the position (XAj, YAj, ZAj) of this pixel corresponds to the absorption point T2 in each of FIGS. 2 (a), 2 (b), and 2 (c).
  • the time Tj corresponds to the time when the X-ray enters the absorber 42.
  • the identifying unit 54 receives the X-ray incident on the scattering unit 41 based on the time Ti at which the electrical signal used to calculate the energy E1 is output and the time Tj at which the electrical signal used to calculate the energy E2 is output. It is determined whether or not the X-ray incident on the absorption unit 42 is an X-ray emitted at a certain timing from the X-ray source 2. Specifically, the specifying unit 54 determines the energy E1 obtained by the scattering unit 41 and the energy obtained by the absorbing unit 42 when the time difference between the times Ti and Tj (that is, Tj-Ti) is within the predetermined time T1. It is determined that E2 is output for the same X-ray photon.
  • the predetermined time T1 is determined based on the time required for the X-ray to reach from the scattering unit 41 to the absorbing unit 42. For example, the predetermined time T1 is determined by calculation or various tests, simulations, etc., and is stored in advance in a memory (not shown).
  • Tij which is a time difference between time Ti and Tj, exceeds the predetermined time T1
  • the identifying unit 54 determines that the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42 have the same X value. It is determined that the line photon is not output.
  • the identifying unit 54 identifies the scattered X-ray of the energy E1 output from the scattering unit 41 and the energy E2 output from the absorbing unit 42. Do not use for processing.
  • the distance in the Z direction between the scattering portion 41 and the absorbing portion 42 is 20 mm
  • the time difference Tij between the time Ti and Tj is about 0.1 nsec.
  • the distance between the scattering unit 41 and the absorbing unit 42 may be determined as a preferable value based on each pixel size of the detector 4 and the detection accuracy ⁇ 1 represented by the above-described equation (3). it can.
  • the identifying unit 54 determines the Compton scattering based on the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42. It is determined whether the angle ⁇ 1 should be calculated. The case where the Compton scattering angle ⁇ 1 should not be calculated will be described next.
  • the readout circuit of the detector 4 outputs an electrical signal used for calculation of the energy E1 and the energy E2 from the scattering unit 41 and the absorbing unit 42 at predetermined time intervals. Referred to herein as the predetermined time and the energy output interval T O.
  • T o 100 nsec.
  • the above predetermined time T1 needs to be set to less than the energy output interval T O. That is, T1 ⁇ T O Set to meet
  • the electrical signal used for calculating the respective from the scattering portion 41 and the absorbent portion 42 energy E1 and the energy E2 may as output.
  • the X-ray source 2 emits X-ray photons one by one when viewed on a short time scale.
  • the relationship between the energy output interval T O for outputting an electric signal to be used from the X-ray source 2 and the average emission interval Te that emits X-ray photon from the detector 4 to calculate the energy will be described.
  • the X-ray photons emitted from the X-ray source 2 by one radiation enter the pixel 41a and the pixel 42a, respectively, and are used for calculating the energy from all the pixels of the scattering portion 41 and the absorbing portion 42 of the detector 4. The process of completing the signal output is called one event.
  • obtaining energy information for one frame from the detector 4 based on a single X-ray photon emitted from the X-ray source 2 is referred to as one event.
  • the Compton scattering angle ⁇ 1 in the scattering section 41 only one detection value of a single X-ray photon emitted from the X-ray source 2 is used for the processing described later during one event. This is because there is a possibility that the Compton scattering angle ⁇ 1 can not be correctly determined when X-ray photons of two or more radiations simultaneously enter the detector 4. This is called event separation.
  • the radiation interval Te is required to be greater than the energy output interval T O is the time required to obtain the energy information of one frame. That is, T O ⁇ Te It is necessary to meet
  • FIG. 3 shows the case where the radiation interval Te from the X-ray source 2 is smaller than the energy output interval TO. Note that FIG. 3 shows the case where two X-ray photons are emitted from the X-ray source 2 at substantially the same time in different directions, in order to simplify the description.
  • One X-ray photon (first X-ray photon) travels along the path L11 from the X-ray source 2, passes through the object S, and enters the pixel 41a of the scattering portion 41.
  • the first X-ray photons are Compton scattered at the scattering point T31 of the pixel 41a, the path changes to L12, and the light enters the pixel 42a of the absorbing unit 42.
  • the X-ray source 2 emits another X-ray photon (second X-ray photon) a short time after the first X-ray photon is emitted.
  • the second X-ray photon travels along a path L21 different from the path L11, passes through the object S, and enters the pixel 41b of the scattering portion 41 of the scattering portion 41.
  • the second X-ray photons are Compton scattered at the scattering point T32 of the pixel 41b, the path changes to L22, and the light enters the pixel 42a of the absorbing unit 42.
  • the identifying unit 54 determines whether or not the sum of the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42 exceeds a predetermined value ETL in a single event.
  • the value of ETL is set based on the energy of a single X-ray photon emitted from the X-ray source 2 in a normal state. For example, the energy resolution is added to the maximum value of the X-ray energy based on the energy of the X-ray source 2 determined by the driving acceleration voltage.
  • the specifying unit 54 determines that the energy E1 and the energy E2 of the plurality of X-ray photons are calculated as shown in FIG. . In this case, the specifying unit 54 does not use the energy E1 corresponding to the two pixels obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42 for specifying the scattered X-ray.
  • the sum of the energy E1 and the energy E2 is equal to or less than the predetermined value ETL, and each corresponds to a single pixel in each of the scattering portion 41 and the absorbing portion 42.
  • the specifying unit 54 uses the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42 to calculate the Compton scattering angle ⁇ 1 based on the above-mentioned equation (1). Then, as described with reference to FIG.
  • the identifying unit 54 does not use the energy E1 and the energy E2 for identifying the scattered X-ray.
  • the specifying unit 54 uses the calculated Compton scattering angle ⁇ 1 and the acquired positions (XSi, YSi, ZSi) and positions (XAj, YAj, ZAj) to form the conical surface CS shown in FIG. To determine whether the X-ray source 2 is present. Specifically, the identifying unit 54 determines the position (XSi, YSi, ZSi) and the position (XSi, YSi, ZSi), the position (XAj, YAj, ZAj), and the calculated Compton scattering angle ⁇ 1.
  • a position (XSi, YSi, ZSi) is a vertex, and a cone (straight cone) in which the angle between the generatrix and the central axis is ⁇ 1 is set.
  • the identifying unit 54 determines whether the X-ray source 2 is present on the conical surface CS of this cone or at a position where the distance from the conical surface CS is smaller than a predetermined distance.
  • the predetermined length is appropriately set based on the sizes of the pixels 41a and 42a and the arrangement of the X-ray source 2, the object S, the detector 4 and the like in the X-ray apparatus 100 as described above. Is preferred.
  • the identifying unit 54 determines that the X-ray source 2 is present on the conical surface CS or at a position where the distance from the conical surface CS is smaller than a predetermined distance, the X-rays emitted from the X-ray source 2 After being incident on the scattering portion 41 without being scattered inside the object S and Compton scattered by the scattering portion 41, it is judged that it has reached the absorbing portion 42, and it is specified that it is an X-ray other than the scattered X-ray Do.
  • the image generation unit 53 generates X-ray intensity distribution data based on the energy information (electric signal) obtained from the scattering unit 41 and the absorption unit 42, and generates X-ray projection image data.
  • the image generation unit 53 uses the internal information of the object S.
  • the identifying unit 54 determines that the X-ray source 2 does not exist on the conical surface CS or at a position where the distance from the conical surface CS is smaller than a predetermined distance, the X-ray emitted from the X-ray source 2
  • the light scatters inside the object S and enters the scattering part 41, and after Compton scattering by the scattering part 41, it is judged that it has reached the absorbing part 42, and the scattered X-ray is specified.
  • the image generation unit 53 does not generate X-ray intensity distribution data based on the energy information (electrical signal) obtained from the scattering unit 41 and the absorption unit 42. That is, the image generation unit 53 does not use it as internal information of the object S when generating the X-ray projection image data.
  • the identification unit 54 transmits the X-rays transmitted through the object S based on the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorption unit 42. At least one of the transmitted X-rays) which are straightly moved without being irradiated to the object S and the scattered X-rays scattered by the object S is specified.
  • the scattering portion 41 detects the scattering position at which the X-ray is scattered in the detector 4, and the absorbing portion 42 detects the incident position at which the X-ray is incident. Based on the position and the incident position, at least one of the transmitted X-ray and the scattered X-ray can be identified.
  • the X-ray projection image data generated by the image generation unit 53 is suppressed from blurring of the edge or the like of the object S due to the scattered X-rays that have caused scattering, etc. As a result, a high quality image of the object S can be generated.
  • a part of X-rays emitted from the X-ray source 2 and incident on the scattering unit 41 of the detector 4 may be incident on the absorbing unit 42 without being Compton scattered by the scattering unit 41.
  • the image generation unit 53 does not generate X-ray intensity distribution data based on the energy information (electrical signal) output from the absorption unit 42. That is, the image generation unit 53 does not use it as internal information of the object S when generating the X-ray projection image data.
  • the image reconstructing unit 56 is known for a plurality of X-ray projection image data generated for each predetermined measurement angle while changing the irradiation angle of X-rays irradiated from the X-ray source 2 to the object S. An image reconstruction process is performed to generate a reconstructed image. As a result, cross-sectional image data and three-dimensional data, which are the internal structure (cross-sectional structure) of the object S, are generated.
  • step S1 the X-ray control unit 51 of the control device 5 causes the X-ray source 2 to start irradiation of X-rays, and proceeds to step S2.
  • step S2 the specifying unit 54 of the control device 5 outputs the pixel at the position (XSi, YSi, ZSi) of the scattering unit 41 of the detector 4 and the pixel at the position (XAj, YAj, ZAj) of the absorbing unit 42.
  • the X-ray energies E1 and E2 are acquired, and the process proceeds to step S3.
  • step S3 the identifying unit 54 determines whether or not the time difference ⁇ Tij between the time Ti when the X-ray enters the scattering unit 41 and the time Tj when the X-ray enters the absorbing unit 42 is less than or equal to a predetermined value T1. . If the time difference ⁇ Tij between the time Tj and Ti is less than or equal to the predetermined value T1, that is, it is determined that the X-ray incident on the scattering portion 41 and the X-ray incident on the absorbing portion 42 are in the same event. If yes, step S3 is affirmed and the process proceeds to step S4.
  • step S3 is negative and the process proceeds to step S9 described later.
  • step S4 it is determined whether or not the X-ray incident on the absorbing unit 42 has passed through different positions of the object S.
  • the sum ET of the energy E1 from the scattering unit 41 and the energy E2 from the absorbing unit 42 exceeds the predetermined value ETL, it is determined that the X-ray has passed through different positions of the object S, and the determination in step S4 is negative. Then, the process proceeds to step S9 described later.
  • step S4 it is determined whether an electrical signal corresponding to the energy E1 or E2 is output from a single point (single pixel) of the detector 4.
  • step S5 When an electric signal corresponding to the energy E1 or E2 is output from a single place in each of the scattering unit 41 and the absorbing unit 42 of the detector 4, an affirmative determination is made in step S5, and the process proceeds to step S6. If electrical signals corresponding to the energy E1 and E2 are output from a plurality of locations (plural pixels) in at least one of the scattering portion 41 and the absorbing portion 42, the negative determination is made in step S5, and the process proceeds to step S9 described later.
  • step S6 using the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42, the specifying unit 54 calculates the Compton scattering angle ⁇ 1 based on the above-mentioned equation (1). Then, the process proceeds to step S7.
  • the identifying unit 54 determines whether or not X-rays have been scattered inside the object S. Specifically, from the calculated Compton scattering angle ⁇ 1, the position (XSi, YSi, ZSi) is taken as a vertex, with the straight line passing through the position (XSi, YSi, ZSi) and the position (XAj, YAj, ZAj) as the central axis.
  • the angle between the generatrix and the central axis sets a cone (straight cone) of ⁇ 1.
  • the identifying unit 54 determines whether the X-ray source 2 is present on the conical surface CS of this cone or at a position where the distance from the conical surface CS is smaller than a predetermined distance. If it is determined that the X-ray is not scattered inside the object S, an affirmative determination is made in step S7, and the process proceeds to step S8. On the other hand, when it is determined that the X-ray has been scattered inside the measured object S, a negative determination is made in step S7, and the process proceeds to step S9 described later.
  • step S8 the image generation unit 53 determines the position (XSi, YSi, ZSi) of the X-ray that has entered the scattering unit 41, the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42.
  • the information is stored as internal information to be an incident position, and the process proceeds to step S9. That is, in a process to be described later, the image generation unit 53 generates X-ray projection image data based on the position (XSi, YSi, ZSi) of the X-ray incident on the scattering unit 41.
  • step S9 the control device 5 determines whether or not the identification process has been performed for all the events of the detector 4.
  • step S9 If the identification process has been performed for all the events, an affirmative determination is made in step S9 and the process proceeds to step S10.
  • step S10 the image generation unit 53 generates X-ray projection image data of the object S using the generated internal information, and proceeds to step S12 described later.
  • step S9 the determination in step S9 is negative and the process proceeds to step S11.
  • the identifying unit 54 outputs the electric signal (pixel from the pixel at the position (XSi, YSi, ZSi) of the scattering unit 41 and the pixel at the position (XAj, YAj, ZAj) at the next event. That is, the electric signal corresponding to the energy E1 and E2 is acquired, and the process returns to step S3 and the same processing is performed thereafter.
  • step S12 the control device 5 determines whether generation of X-ray projection image data has been performed for all measurement angles. If generation of X-ray projection image data has been performed for all measurement angles, an affirmative determination is made in step S12, and the process proceeds to step S13.
  • step S13 the image reconstruction unit 56 generates three-dimensional data using a plurality of X-ray projection image data, and ends the process. The control device 5 can display this three-dimensional data on a monitor (not shown) or store it in a memory (not shown). If there is a measurement angle at which X-ray projection image data is not generated in step S12, the determination in step S12 is negative and the process proceeds to step S14.
  • step S14 the placement stand control unit 52 drives the manipulator unit 36 to rotationally drive the placement stand 31 to a predetermined measurement angle, and the process returns to step S2, and the same processing is performed thereafter.
  • the identifying unit 54 of the control device 5 determines the object S based on the energy loss E1 of the X-ray obtained by the scattering unit 41 and the energy E2 of the scattered X-ray obtained by the absorbing unit 42. Among the X-rays incident on the detector 4, the scattered X-rays scattered inside the object S are specified. Further, the identifying unit 54 of the control device 5 sets the scattering position of the X-ray in the detector 4 detected by the scattering unit 41 and the incident position to the absorbing unit 42 of the X-ray detected by the absorbing unit 42. Based on this, at least one of transmitted X-rays and scattered X-rays can be identified.
  • the identifying unit 54 determines the scattering angle ⁇ 1 of the X-rays scattered by the scattering unit 41, the position (scattering position) of the X-rays incident on the scattering unit 41, and the position of the scattered X-rays incident on the absorbing unit Scattered X-rays are identified based on the relationship satisfied with the incident position). As a result, it is possible to generate X-ray projection image data by excluding scattered X-rays among the X-rays incident on the detector 4.
  • the identifying unit 54 is based on the scattering angle ⁇ 1, the position of the X-ray incident on the scattering unit 41 (ie, the scattering position), and the position of the X-ray incident on the absorbing unit 42 (ie, the incident position).
  • the position where the presence of the X-ray source 2 is estimated (circumference C2) is calculated, and when the X-ray source 2 exists at the calculated position, the X-rays incident on the detector 4 are scattered X by the object S Identify as transmitted X-rays other than lines.
  • the image quality of the X-ray projection image data can be improved.
  • the specifying unit 54 has a predetermined time difference between the time Ti when the electric signal corresponding to the energy E1 is output from the scattering unit 41 and the time Tj when the electric signal corresponding to the energy E2 from the absorbing unit 42.
  • the scattered X-rays are identified using X-rays within the range T1.
  • the X-rays emitted from the X-ray source 2 at different timings and excluding the X-rays incident on the detector 4 can be identified to identify the scattered X-rays, so that the accuracy of identifying the scattered X-rays can be improved.
  • the identifying unit 54 identifies the scattered X-rays by irradiating the X-ray source 2, transmitting the different positions of the object S, and excluding the X-rays incident on the same pixel of the absorbing unit 42. Do. As a result, it is suppressed that the information of different positions inside the object S is included in the internal information as the information of the same position, so it is possible to suppress the deterioration of the image quality of the X-ray projection image data.
  • the specifying unit 54 determines that the sum of the energy E1 obtained by the scattering unit 41 and the energy E2 obtained by the absorbing unit 42 exceeds a value based on the energy E0 of the X-ray emitted from the X-ray source 2 And the different positions of the object S are transmitted to detect that the plurality of X-ray photons are incident on the same position of the absorber 42, and the scattered X-ray is specified excluding the detected X-ray photons Do. Thereby, X-rays transmitted through different positions of the object to be measured S and incident on the same pixel of the absorbing unit 42 can be excluded from the internal information of the object to be measured S.
  • the scattering portion 41 and the absorbing portion 42 of the detector 4 are described as an example of the scintillator portion containing the scintillation material. It is not limited to the example.
  • the scattering unit 41 and the absorbing unit 42 of the detector 4 may be semiconductor detectors that directly convert an incident X-ray into an electric signal without converting it into light energy and output the electric signal.
  • the scattering unit 41 and the absorbing unit 42 have CdTe, generate electrons in a number proportional to the energy of the incident X-ray, and output as an electric signal.
  • the thickness of the scattering portion 41 in the Z direction is smaller than the thickness of the absorbing portion 42 in the Z direction.
  • the image generation unit 53 generates X-ray projection image data based on the position (XSi, YSi, ZSi) of X-rays incident on the scattering unit 41
  • the explanation is given taking the example as an example.
  • the image generation unit 53 may generate X-ray projection image data based on the position (XAj, YAj, ZAj) of the X-ray incident on the absorption unit 42.
  • the image generation unit 53 assumes that the scattering unit 41 does not scatter X-rays based on the position of the X-ray source 2 and the position (XSi, YSi, ZSi) of X-rays incident on the scattering unit 41
  • the position (XAj ', YAj', ZAj ') on the absorbing portion 42, which is incident when it is input, is calculated.
  • the image generation unit 53 may generate X-ray projection image data from the energy E1 and the energy E2 based on the calculated positions (XAj ', YAj', ZAj ').
  • the structure manufacturing system of the present embodiment produces a molded article such as an electronic component provided with, for example, a door portion of an automobile, an engine portion, a gear portion, and a circuit board.
  • FIG. 5 is a block diagram showing an example of the configuration of a structure manufacturing system 400 according to the present embodiment.
  • the structure manufacturing system 400 includes the X-ray device 100 described in the first embodiment, a design device 410, a molding device 420, a control system 430, and a repair device 440.
  • the design device 410 is a device used by a user when creating design information on the shape of a structure, and performs design processing for creating and storing the design information.
  • the design information is information indicating the coordinates of each position of the structure.
  • the design information is output to the molding apparatus 420 and a control system 430 described later.
  • the forming apparatus 420 performs a forming process of forming and forming a structure using the design information generated by the design apparatus 410.
  • the forming apparatus 420 is also included in an aspect of the present invention for performing at least one of lamination processing, casting processing, forging processing, and cutting processing represented by 3D printer technology.
  • the X-ray apparatus 100 performs measurement processing for measuring the shape of the structure formed by the forming apparatus 420.
  • the X-ray apparatus 100 outputs, to the control system 430, information (hereinafter referred to as shape information) indicating the coordinates of the structure, which is the measurement result of measuring the structure.
  • shape information information (hereinafter referred to as shape information) indicating the coordinates of the structure, which is the measurement result of measuring the structure.
  • the control system 430 includes a coordinate storage unit 431 and an inspection unit 432.
  • the coordinate storage unit 431 stores the design information created by the design device 410 described above.
  • the inspection unit 432 determines whether the structure formed by the forming device 420 is formed according to the design information created by the design device 410. In other words, the inspection unit 432 determines whether or not the formed structure is non-defective. In this case, the inspection unit 432 performs an inspection process of reading the design information stored in the coordinate storage unit 431 and comparing the design information with the shape information input from the X-ray apparatus 100. The inspection unit 432 compares, for example, the coordinates indicated by the design information with the coordinates indicated by the corresponding shape information as the inspection process, and as a result of the inspection process, if the coordinates of the design information and the coordinates of the shape information match. It is determined that the non-defective product is molded according to the design information.
  • the inspection unit 432 determines whether the difference between the coordinates is within a predetermined range, and if within the predetermined range, restoration is possible. It is determined that the product is defective.
  • the inspection unit 432 outputs repair information indicating the defective portion and the amount of repair to the repair device 440.
  • the defective portion is a portion having coordinates of shape information that does not match the coordinates of design information, and the amount of repair is a difference between the coordinates of design information at the defective portion and the coordinates of shape information.
  • the repair device 440 performs repair processing to rework defective parts of the structure based on the input repair information. The repair device 440 performs again the same processing as the molding processing performed by the molding device 420 in the repair processing.
  • step S31 the design device 410 is used when the user designs a structure, and design information about the shape of the structure is created and stored by design processing, and the process proceeds to step S32.
  • the present invention is not limited to only the design information created by the design device 410, and in the case where design information is already input, one that acquires design information by inputting the design information is also included in an aspect of the present invention.
  • the forming apparatus 420 forms and structures a structure based on the design information by the forming process, and proceeds to step S33.
  • step S33 the X-ray apparatus 100 performs measurement processing, measures the shape of the structure, outputs shape information, and proceeds to step S34.
  • step S34 the inspection unit 432 performs an inspection process of comparing the design information generated by the design apparatus 410 with the shape information measured and output by the X-ray apparatus 100, and the process proceeds to step S35.
  • step S35 the inspection unit 432 determines whether or not the structure formed by the forming device 420 is non-defective based on the result of the inspection process. If the structure is non-defective, that is, if the difference between the coordinates of the design information and the coordinates of the shape information is within a predetermined range, an affirmative determination is made in step S35, and the process is ended.
  • step S35 is negative and the process proceeds to step S36.
  • step S36 the inspection unit 432 determines whether or not the defective portion of the structure can be repaired. If the defective portion is not recoverable, that is, if the difference between the coordinate of the design information and the coordinate of the shape information at the defective portion exceeds the predetermined range, the determination in step S36 is negative and the process is ended. If the defective portion can be repaired, that is, if the difference between the coordinate of the design information and the coordinate of the shape information in the defective portion is within the predetermined range, step S36 is positively determined, and the process proceeds to step S37. In this case, the inspection unit 432 outputs repair information to the repair device 440. In step S37, the repair device 440 performs a repair process on the structure based on the input repair information, and returns to step S33. As described above, the repair device 440 performs again the same processing as the molding processing performed by the molding device 420 in the repair processing.
  • the X-ray apparatus 100 of the structure manufacturing system 400 performs measurement processing for acquiring the shape information of the structure created by the forming device 420 based on the design processing of the design device 410, and the inspection unit of the control system 430
  • a step 432 performs an inspection process to compare the shape information acquired in the measurement process with the design information created in the design process. Therefore, the quality control of the structure can be performed because the defect of the structure and the internal information of the structure can be acquired by nondestructive inspection, and it can be determined whether the structure is a non-defective item created according to the design information.
  • the repair device 440 performs the repair process for performing the forming process on the structure again based on the comparison result of the inspection process. Therefore, when the defective portion of the structure can be repaired, the same process as the forming process can be applied to the structure again, which contributes to the production of a high quality structure close to the design information.
  • FIG. 7 is a diagram showing an example of the configuration of the X-ray apparatus 100A according to the second embodiment.
  • a coordinate system including an X axis, a Y axis, and a Z axis is set as in FIG. 1.
  • the X-ray apparatus 100A includes a housing 1, an X-ray source 2, a placement unit 3, a control device 5, and a detector 4A as in the first embodiment.
  • the detector 4A is provided on the Z axis + side with respect to the mounting table 31. That is, the mounting table 31 is provided between the X-ray source 2 and the detector 4 in the Z-axis direction.
  • the detector 4A includes a scattering unit 41A and an absorbing unit 42A.
  • the scattering portion 41A is provided on the Z direction ⁇ side, that is, the side closer to the object S than the absorbing portion 42A. X-rays emitted from the X-ray source 2 and transmitted through the object to be measured S placed on the mounting table 31 or scattered in the object to be measured S enter the scattering portion 41A.
  • X-rays transmitted through the inside of the object S without being scattered and entering the detector 4 without being scattered are transmitted X-rays, and are scattered inside the object S and the detector 4A X-rays incident on are called scattered X-rays.
  • the scattering portion 41A has, for example, the same configuration as that of the known Japanese Patent Application Laid-Open No. 2017-67601. That is, as shown in FIG. 8, the scattering section 41A has a chamber 411A, a pixel electrode 412A and a drift electrode 413A. Inside the chamber 411A, a mixed gas of a rare gas such as argon or xenon and a gas having a quenching action including alkane or carbon dioxide such as ethane or methane at normal temperature is enclosed.
  • the chamber 411A is not limited to one in which the mixed gas is sealed, and may be one in which a single gas is sealed.
  • the pixel electrode 412A is provided on the surface on the Z direction + side of the chamber 411A, and the drift electrode 413A is provided on the surface on the Z direction ⁇ side of the chamber 411.
  • the pixel electrode 412A has a plurality of pixels two-dimensionally arranged on the XY plane.
  • an electric field is formed between the pixel electrode 412A and the drift electrode 413A by the voltage applied to the pixel electrode 412A and the drift electrode 413A.
  • the scattering section 41 When X-rays from the object S are incident on the scattering section 41, the X-rays (photons) collide with electrons constituting the gas particles in the chamber 411A with a certain probability, and are Compton scattered.
  • the Compton-scattered X-ray photons exit the scattering portion 41A in a state in which the traveling direction changes, and enter the absorbing portion 42A provided on the Z direction + side.
  • X-rays not Compton scattered in the chamber 411A exit the scattering portion 41A without changing the traveling direction in the chamber 411A and enter the absorbing portion 42A.
  • the secondary electrons constituting the electron cloud are attracted to the Z direction + side by the electric field formed by the pixel electrode 412A and the drift electrode 413A, and enter the pixel electrode 412A.
  • the pixel electrode 412A has a plurality of pixels each including a cathode electrode extending in the X direction and an anode electrode disposed in a plurality of openings provided in the cathode electrode, and secondary electrons forming an electron cloud are incident thereon
  • the pixel outputs an electrical signal to the controller 5.
  • the scattering unit 41 can detect information on the progress of electrons generated by the scattering of the X-ray incident in the chamber 411A. That is, the shape of the electron cloud is obtained.
  • the information on the electron's travel includes the initial travel direction of the recoil electrons (the initial travel direction of the electrons collided with X-ray photons) and the travel trajectory of the recoil electrons represented by the distribution of many secondary electrons. At least one of (track) is included.
  • the absorbing unit 42A is configured of a scintillator unit including a known scintillation material, a light receiving unit that amplifies and receives light emitted from the scintillator unit, and the like.
  • the X-ray incident on the incident surface of the scintillator unit is absorbed by the scintillator unit of the absorbing unit 42A and emits fluorescence, and the light energy of the emitted fluorescence is amplified by the photomultiplier in the above-described light receiving unit to be electrical energy It is converted and output to the control device 5 as an electric signal.
  • the absorbing unit 42A of the detector 4A has a structure in which a plurality of pixels are two-dimensionally arranged in the scintillator unit and the light receiving unit, and the plurality of pixels in the scintillator unit and the light receiving unit correspond to each other Are arranged in By this, it is possible to collectively acquire the intensity data of X-rays emitted from the X-ray source 2 and transmitted, scattered or partially absorbed by the object S at a plurality of pixels, and the respective intensity data
  • the intensity distribution data can be obtained by integrating.
  • Each of the plurality of pixels of the absorbing portion 42A is arranged to correspond to each of the plurality of pixels of the pixel electrode 412A of the scattering portion 41A.
  • the scintillator portion of the absorbing portion 42A may not have a structure in which a plurality of pixels are two-dimensionally arranged.
  • the identifying unit 54 of the control device 5 estimates the incident direction of the X-ray to the detector 4A, and transmits X-ray (or X-ray) based on the estimated incident direction and the position of the X-ray source 2 At least one of transmitted X-rays) and scattered X-rays which go straight without being irradiated to the object S is specified. That is, the identifying unit 54 determines whether the X-rays from the X-ray source 2 are scattered inside the object S based on the electrical signals output from the scattering unit 41A and the absorbing unit 42A.
  • the identifying unit 54 identifies, among X-rays from the object S, scattered X-rays scattered at the object S or identifies X-rays other than the scattered X-rays, based on the determination result. That is, the identifying unit 54 identifies whether the X-ray incident on the detector 4A is the X-ray scattered by the object S.
  • FIG. 8 shows the case where the electrical signals are detected by the plurality of pixels 412a to 412c of the pixel electrode 412A of the scattering portion 41A, and the electrical signals are detected by one pixel 42a of the absorbing portion 42A. Note that it is assumed that the pixel electrode 412A detects an electrical signal in the order of the pixels 412c, 412b, and 412a.
  • the pixel electrode 412A secondary electrons constituting an electron cloud generated along a track of recoil electrons generated by Compton scattering of X-rays in the scattering portion 41A are drift electric fields in the chamber 411A.
  • the electrons are attracted to the Z direction + side, and when the electrons reach the pixel electrode 412A, the signal is detected.
  • an electron cloud which is information on the traveling of recoil electrons in the chamber 411A, according to the time when electrons are sequentially detected in the pixels 412c, 412b, and 412a and the two-dimensional position of the pixels 412c, 412b, and 412a.
  • the three-dimensional shape (track of electrons) of can be determined.
  • the procedure for obtaining the three-dimensional shape of the electron cloud will be specifically described.
  • the moving velocity Vd at which any one secondary electron contributing to the generation of the electron cloud moves in the Z direction along the drift electric field is determined by the strength of the drift electric field.
  • the moving velocity Vd is considered to be substantially constant in the chamber 411A.
  • the time for drifting from the drift electrode 413A to the pixel electrode 412A is Td. When the electrons reach the pixel electrode 412A, a signal is detected.
  • the time taken for the X-rays to reach the absorbing part 42A after Compton scattering and the time required for the generation of an electron cloud due to this are extremely short compared to Td, so the time of Compton scattering and the scattered X-rays There is no problem considering that the time when the light reaches the absorbing portion 42A and the time when the electron cloud is formed are substantially the same.
  • the thickness D of the chamber 411A of the scattering portion 41A in the Z direction is 300 mm
  • the electron causing the voltage change in the pixel (Xi, Yj) of the pixel electrode 412A starts drifting, that is, the position where the electron cloud is generated in the chamber 411A is the pixel electrode 412A as the Z direction reference It is considered that (Xi, Yj, Vd ⁇ t sig ).
  • the reference in the Z direction may be the position of the drift electrode 413A.
  • the position in the chamber 411A corresponding to the pixel whose voltage has changed is (Xi, Yj, D (Td ⁇ t sig ) / Td) It becomes.
  • the scattering section 41A functions as a first detector that detects the initial traveling direction F1 of recoil electrons generated by the Compton scattering of X-rays at the scattering point T1 in the chamber 411A. Accordingly, the scattering unit 41A detects the initial traveling direction F1 of recoil electrons as information on the progress of the electrons generated by the scattering of the X-rays in the chamber 411A.
  • the X-rays Compton scattered at the scattering point T1 change in the traveling direction and enter the pixel 42a of the absorbing portion 42A. That is, the pixel 42a of the absorbing portion 42A is the incident position of X-rays incident on the absorbing portion 42A, and the absorbing portion 42A receives the incident X-rays from the scattering portion 41A and detects the incident position of the incident X-rays. It functions as a second detector. In this way, it is possible to grasp the traveling direction of the Compton scattered X-ray and the shape of the electron cloud generated by recoil electrons due to Compton scattering.
  • the data regarding the correspondence of a range and energy shall be beforehand stored in memory (not shown) as energy data.
  • h is Planck's constant.
  • the X-rays incident on the chamber 411A, the Compton-scattered X-rays, and the recoil electrons are expressed by the following equations (4) to (6) using the momentum conservation law and the energy conservation law.
  • c is the speed of light in vacuum
  • is the Compton scattering angle
  • is the emission angle of the recoil electrons (azimuth angle in the initial traveling direction).
  • is a Lorentz factor and is expressed by the following equation (7).
  • the scattering angle ⁇ is represented by the following equation (8) according to the above equations (4) to (7).
  • the identifying unit 54 detects the scattering point T1 based on the positional relationship between the pixels 412a and 412b corresponding to a position near at least the projection position of the scattering point T1 in the Z direction among the pixels 412a to 412c of the pixel electrode 412A. It is possible to obtain the traveling direction (azimuth angle of the initial traveling direction) of the recoil electrons emitted from. Preferably, the azimuth may be obtained based on the output from a plurality of pixel groups linearly arranged corresponding to the position until the recoil electron performs the first Rutherford scattering. Can be improved.
  • the identifying unit 54 estimates the incident direction of the X-ray to the detector 4A, and based on the estimated incident direction and the position of the X-ray source 2, the transmitted X-ray (or the straight object without irradiating the object S At least one of transmitted X-rays and scattered X-rays. The details will be described below.
  • the path L2 is a path connecting the position of the scattering point T1 and the position of the pixel 42a by a straight line.
  • a plane (scattering surface) including the path L2 of the scattered X-rays and the initial traveling direction F1 of the recoil electrons is determined by the path L2 and the initial traveling direction F1 of the recoil electrons emitted from the scattering point T1.
  • the scattering direction which is the traveling direction of the X-ray from the scattering point T1 is determined based on the position of the scattering point T1 and the position of the pixel 42a, and the scattering surface based on the scattering direction and the initial traveling direction F1 of recoil electrons.
  • the scattering direction may be calculated based on the path L2. Also, although an example was shown where the scattering plane was determined using the initial traveling direction F1 of recoil electrons, the scattering plane was determined using the track P1 of recoil electrons and the path L2 of scattering X-rays or the scattering direction. May be Since the path L1 of the X-ray incident on the chamber 411A exists in the above-described scattering surface, the path L1 and the path L2 intersect at the scattering point T1 at the scattering point T1 at the scattering point T1.
  • FIG. 9A shows a case where electric signals are output from the pixels 412a to 412c of the scattering portion 41A and the pixel 42a of the absorbing portion 42A, as in FIG. In such a case, the X-ray source 2 is present at or near the extension of the path L1.
  • FIG. 9B shows the case where electric signals are output from the pixels 412a and 412f of the scattering portion 41A and the pixel 42a of the absorbing portion 42A.
  • the initial traveling direction of recoil electrons is F2
  • X-rays are scattered at a scattering point T3 inside the object S.
  • the X-ray source 2 is not present on or near the extension of the path L1.
  • the identifying unit 54 can calculate the azimuth angle ⁇ of the initial advance of the recoil electrons using the following equation (9) obtained from the above equations (4) to (7).
  • the azimuth (exit angle) ⁇ is calculated using the equation (9)
  • the path L1 has the scattering point T1 as a vertex and the initial traveling direction F1 of recoil electrons as a central axis, and the central axis and generatrix
  • the angle formed by is the cone of ⁇
  • the scattering point T1 is the apex
  • the path L2 is the central axis
  • the angle between the central axis and the generatrix corresponds to the tangent to the cone of ⁇ .
  • the path L1 of the X-ray incident on the chamber 411A can be determined, but both the scattering angle ⁇ and the azimuth angle ⁇ are calculated. Thus, it is possible to determine the course L1 with higher accuracy.
  • the identifying unit 54 determines that the X-ray is not scattered by the object S and is scattered. Identify X-rays other than X-rays. On the other hand, when the X-ray source 2 is not located on the path L1 and in the vicinity of the path L1, the identifying unit 54 determines that the X-rays are scattered by the object S and indicates that they are scattered X-rays. Identify.
  • the identifying unit 54 for identifying whether or not X-rays are scattered by the object S, and identifying scattered X-rays and X-rays other than the scattered X-rays (hereinafter referred to as identifying process Will be described.
  • the Compton scattering angle ⁇ calculated by the equation (8) satisfies the following equation (10).
  • the detection accuracy ⁇ is calculated by the following equation (12). Note that Shall be For example, assuming that the energy of X-rays emitted from the X-ray source 2 and incident on the scattering section 41A is 511 keV, the energy resolution is 3 percent, and the Compton scattering angle ⁇ is 20 °, the detection accuracy ⁇ is 5.3 ° Become. When the energy resolution is 1 percent, the detection accuracy ⁇ is 1.8 °.
  • the distance in the Z direction between the exit surface (pixel electrode 412 side) of the scattering portion 41A and the absorbing portion 42A is 20 mm, and the size of each pixel of the pixel electrode 412A of the scattering portion 41A and the absorbing portion 42A is 200 ⁇ m ⁇ 200 ⁇ m. Assuming that the resolution of the scattering angle ⁇ is 10 mrad.
  • the identifying unit 54 calculates the energy E1 calculated based on the output of the scattering unit 41A of the detector 4A, the position of the pixel of the pixel electrode 412A that outputs the electrical signal, and information of the time t sig when the electrical signal is output from the pixel To get The identifying unit 54 calculates a track having a three-dimensional shape of the electron cloud, based on the position of each pixel at which the electrical signal is output at the pixel electrode 412A and the time t sig .
  • the position (Xi, Yj, Vd ⁇ t sig0 ) determined by the above corresponds to the scattering point T1.
  • the scattering point T1 corresponds to a position at which the calculated distance Vd ⁇ t sig is apart from the position of the pixel 412a in the Z-axis direction.
  • the identifying unit 54 outputs the energy E2 obtained from the absorbing unit 42A of the detector 4A, the position (XAj, YAj, ZAj) of the pixel that has output the electrical signal corresponding to the energy E2, and the energy E2 from the absorbing unit 42A.
  • Each piece of information at time Tj is acquired.
  • the position (XAj, YAj, ZAj) of this pixel corresponds to the absorption point T2 in each of FIGS. 8, 9 (a), and 9 (b).
  • the time Tj corresponds to the time when the X-ray enters the absorber 42A.
  • the predetermined time T1 is determined by calculation or various tests, simulations, etc., and is stored in advance in a memory (not shown).
  • the time Tj + t sig is the time difference between the Tj t sig exceeds a predetermined time T1
  • the specific part 54, and the energy E2 obtained by the energy E1 calculated absorbing portion 42A the same X-ray photons It is determined that the output is not about.
  • the identifying unit 54 processes the energy E 1 calculated and the energy E 2 obtained by the absorbing unit 42 A to identify the scattered X-rays.
  • the distance in the Z direction between the scattering portion 41A and the absorbing portion 42A is 20 mm as described above, the time difference t sig between time Tj + t sig and Tj is about 0.1 nsec.
  • the distance between the scattering unit 41A and the absorbing unit 42A may be determined as a preferable value based on each pixel size of the detector 4A and the detection accuracy ⁇ represented by the above-described equation (12). it can.
  • the identifying unit 54 determines the Compton scattering angle ⁇ based on the calculated energy E1 and the energy E2 obtained by the absorbing unit 42A. Determine whether to calculate. The case where the Compton scattering angle ⁇ should not be calculated will be described next.
  • the readout circuit of the detector 4A outputs an electrical signal from the scattering unit 41A and outputs energy E2 from the absorbing unit 42A at predetermined time intervals. Referred to herein as the predetermined time and the energy output interval T O.
  • T o 100 nsec.
  • the above predetermined time T1 needs to be set to less than the energy output interval T O. That is, T1 ⁇ T O Set to meet
  • the energy output within the interval T O the electrical signal is output from the scattering portion 41A, it can be such that the electrical signal corresponding to the energy E2 from the absorbent portion 42A is output.
  • the X-ray source 2 emits X-ray photons one by one when viewed on a short time scale.
  • Energy information from the detector 4 and the average emission interval Te from the X-ray source 2 emits X-ray photons (i.e., electrical signals from the pixels 42a of the pixel electrode 412 and the absorption section 42A) the relationship between the output interval T O for outputting Will be explained.
  • the X-ray photons emitted from the X-ray source 2 by one radiation are respectively incident on the chamber 411A and the absorbing portion 42A, and energy information is calculated from all the pixels of the scattering portion 41A and the absorbing portion 42A of the detector 4A. The process of completing the output of the used electrical signal is called one event.
  • obtaining energy information for one frame from the detector 4 based on a single X-ray photon emitted from the X-ray source 2 is referred to as one event.
  • the radiation interval Te is required to be greater than the energy output interval T O is the time required to obtain the energy information of one frame. That is, T O ⁇ Te It is necessary to meet
  • Figure 10 is a radiation interval Te from the X-ray source 2 shows the case energy output interval T O less.
  • T O the case energy output interval
  • first X-ray photon travels along the path L11 from the X-ray source 2, transmits the object S, enters the chamber 411A of the scattering section 41A, and Compton scatters at the scattering point T31.
  • the path changes to L12, and the light enters the pixel 42a of the absorbing portion 42A.
  • a scattering point T31 exists in the Z-axis direction of the pixel 412a of the pixel electrode 412A.
  • another X-ray photon (second X-ray photon) is emitted from the X-ray source 2 immediately after (almost simultaneously) the emission of the first X-ray photon.
  • the second X-ray photon travels on a path L21 different from the path L11, passes through the object S, enters the chamber 411A, and Compton scatters at the scattering point T32 to change the path to L22.
  • To the pixel 42a of At the scattering point T32 recoil electrons whose initial traveling direction is F4 are generated, and an electrical signal is output from the pixel of the pixel electrode 412A to which the secondary electrons forming the electron cloud generated according to the track have arrived.
  • the scattering point T32 exists in the Z-axis direction of the pixel 412x of the pixel electrode 412A.
  • the identifying unit 54 determines whether the sum of the energy E1 calculated based on the output of the electrical signal from the scattering unit 41A and the energy E2 obtained by the absorbing unit 42A exceeds a predetermined value ETL in a single event. Determine The value of ETL is set based on the energy of a single X-ray photon emitted from the X-ray source 2 in a normal state. For example, the energy resolution is added to the maximum value of the X-ray energy based on the energy of the X-ray source 2 determined by the driving acceleration voltage.
  • the identifying unit 54 determines that the energy E1 and the energy E2 are obtained by the plurality of X-ray photons, as shown in FIG. Do. In this case, the identifying unit 54 does not use the energy E1 calculated based on the output of the electrical signal from the scattering unit 41A and the energy E2 obtained by the absorbing unit 42A for identifying the scattered X-ray.
  • the specifying unit 54 uses the energy E1 calculated based on the output of the electrical signal from the scattering unit 41A and the energy E2 obtained by the absorbing unit 42A to use the above-described equation (8).
  • the Compton scattering angle ⁇ is calculated, and as described with reference to FIGS.
  • the specifying unit 54 uses the calculated Compton scattering angle ⁇ , the initial traveling direction F1 of recoil electrons, the scattering point T1, and the position (XAj, YAj, ZAj) of the pixel 42a detected by the absorbing unit 42A. Then, it is determined whether or not the X-ray source 2 is present on or near the path L1 shown in FIG. Specifically, the identification unit 54 determines the recoil electrons from the initial traveling direction F1 of the recoil electrons, the position of the scattering point T1, the position (XAj, YAj, ZAj), and the calculated Compton scattering angle ⁇ .
  • a straight line passing through the position and position (XAj, YAj, ZAj) of the scattering point T1 on the plane including the initial traveling direction F1, the position of the scattering point T1, and the position (XAj, YAj, ZAj) is set as the path L1.
  • the identifying unit 54 determines whether the X-ray source 2 exists at a position on the path L1 or at a position at which the distance from the path L1 is smaller than a predetermined distance.
  • the predetermined distance is appropriately determined based on the size of each pixel of the pixel electrode 412A and the pixel 42a and the arrangement of the X-ray source 2, the object S, the detector 4A, etc. in the X-ray device 100 as described above. It is preferable to set.
  • the identifying unit 54 determines that the X-ray source 2 is present on the path L1 or at a position where the distance from the path L1 is smaller than a predetermined distance, the X-ray emitted from the X-ray source 2
  • the light is incident on the scattering portion 41A without being scattered inside S, and after Compton scattering by the scattering portion 41A, it is determined that it has reached the absorbing portion 42A, and it is specified that it is an X-ray other than the scattered X-ray.
  • the image generation unit 53 generates X-ray intensity distribution data based on the energy information (electrical signal) output from the absorption unit 42A, and generates an X-ray projection image data. Used as internal information of S.
  • the image generation unit 53 uses the internal information of the object S.
  • the identifying unit 54 determines that the X-ray source 2 does not exist on the path L1 or at a position where the distance from the path L1 is smaller than a predetermined distance
  • the X-ray emitted from the X-ray source 2 is After being scattered inside the object to be measured S and incident on the scattering portion 41A and Compton scattered by the scattering portion 41A, it is determined that it has reached the absorbing portion 42A, and the scattered X-ray is specified.
  • the image generation unit 53 does not generate X-ray intensity distribution data based on the energy information (electric signal) output from the absorption unit 42A. That is, the image generation unit 53 does not use it as internal information of the object S when generating the X-ray projection image data.
  • the X-ray projection image data generated by the image generation unit 53 is suppressed from blurring of the edge or the like of the object S due to the scattered X-rays that have caused scattering, etc. As a result, a high quality image of the object S can be generated.
  • a part of X-rays emitted from the X-ray source 2 and incident on the scattering portion 41A of the detector 4A may be incident on the absorbing portion 42A without being Compton scattered by the scattering portion 41A.
  • the image generation unit 53 does not generate X-ray intensity distribution data based on the energy information (electrical signal) output from the absorption unit 42A. That is, the image generation unit 53 does not use it as internal information of the object S when generating the X-ray projection image data.
  • the image reconstructing unit 56 is known for a plurality of X-ray projection image data generated for each predetermined measurement angle while changing the irradiation angle of X-rays irradiated from the X-ray source 2 to the object S. An image reconstruction process is performed to generate a reconstructed image. As a result, cross-sectional image data and three-dimensional data, which are the internal structure (cross-sectional structure) of the object S, are generated.
  • step S101 the X-ray control unit 51 of the control device 5 causes the X-ray source 2 to start irradiation of X-rays, and proceeds to step S102.
  • step S102 the specifying unit 54 of the control device 5 acquires the electrical signal output from the pixel electrode 412A of the scattering unit 41A of the detector 4A, detects the track of recoil electrons, and detects the maximum length of the track and the reverse.
  • step S103 the identifying unit 54 acquires (calculates) the energy E2 of the X-ray based on the electrical signal output from the pixel at the position (XAj, YAj, ZAj) of the absorbing unit 42A, and proceeds to step S104.
  • step S104 the identification unit 54 determines whether or not the time difference ⁇ Tij between the time Tj + t sig when the X-ray enters the scattering part 41A and the time Tj when the X-ray enters the absorbing part 42A is less than or equal to a predetermined value T1. Do. When the time difference ⁇ Tij between time Tj and Tj + t sig is less than or equal to a predetermined value T1, that is, it is determined that the X-ray incident on the scattering portion 41A and the X-ray incident on the absorbing portion 42A are in the same event. If YES in step S104, the process proceeds to step S105.
  • step S104 When the time difference ⁇ Tij between time Tj and Tj + t sig exceeds a predetermined value T1, that is, it is determined that the X-ray incident on the scattering portion 41A and the X-ray incident on the absorbing portion 42A are not in the same event If it is determined that the answer in step S104 is negative, the process proceeds to step S110 described later.
  • step S105 it is determined whether the X-rays incident on the absorbing portion 42A have passed through different positions of the object S.
  • the sum ET of the energy E1 calculated in step S102 and the energy E2 obtained by the absorbing unit 42A acquired in step S103 exceeds the predetermined value ETL, the X-rays are transmitted through different positions of the object S
  • the process proceeds to step S110 described later.
  • the sum ET of the energy E1 and the energy E2 is less than the predetermined value ETL, X-rays are normally emitted from the X-ray source 2, and energy information based on a single X-ray photon is emitted from the scattering portion 41A and the absorbing portion 42A.
  • step S105 a single track is detected by the detector 4A, and it is determined whether an electrical signal corresponding to the energy E2 is output from a single point (single pixel) of the detector 4A.
  • step S106 a single track is detected based on the output of the scattering unit 41A of the detector 4A and energy E2 is output from a single point of the absorbing unit 42A, step S106 is positively determined, and the step is performed. It progresses to S107.
  • step S106 is negatively determined. The process proceeds to step S110 described later.
  • step S107 the identifying unit 54 calculates the Compton scattering angle ⁇ based on the above-described equation (8) using the calculated energy E1 and the energy E2 output from the absorbing unit 42, and proceeds to step S108. move on.
  • step S ⁇ b> 108 the identifying unit 54 determines whether or not X-rays have been scattered inside the object S. Specifically, the identifying unit 54 determines the scattering point T1 and the position (XAj,) on the scattering plane including the initial traveling direction of recoil electrons, the position of the scattering point T1, and the position (XAj, YAj, ZAj).
  • a straight line passing the scattering point T1 is set as a path L1 with a Compton scattering angle ⁇ with respect to a straight line connecting YAj and ZAj).
  • the identifying unit 54 determines whether the X-ray source 2 exists at a position on the path L1 or at a position at which the distance from the path L1 is smaller than a predetermined distance. When it is determined that the X-ray is not scattered inside the object S, an affirmative determination is made in step S108, and the process proceeds to step S109. On the other hand, when it is determined that the X-ray has been scattered inside the measured object S, a negative determination is made in step S108, and the process proceeds to step S110 described later.
  • step S109 the identifying unit 54 calculates the position (expected incident position) on the absorbing unit 42A which is expected to be incident when X-rays incident on the chamber 411A do not cause the Compton scattering.
  • the specifying unit 54 calculates the expected incident position based on the Compton scattering angle ⁇ , the position of the scattering point T1, and the path L1.
  • the image generation unit 53 stores the energy E2 from the absorption unit 42A and the above-described expected incident position calculated by the identification unit 54 as internal information, and proceeds to step S110.
  • step S110 the control device 5 determines whether or not the identification process has been performed for all the events of the detector 4A.
  • step S110 If the identification process has been performed for all the events, an affirmative determination is made in step S110, and the process proceeds to step S111.
  • step S111 the image generation unit 53 generates X-ray projection image data of the object S using the generated internal information, and proceeds to step S113 described later.
  • step S110 If there is an event for which the identification process has not been performed in step S110, the determination in step S110 is negative and the process proceeds to step S112.
  • step S112 the specifying unit 54 acquires the electrical signal from the pixel electrode 412A of the scattering unit 41A and the electrical signal from the absorbing unit 42A in the next event, returns to step S104, and performs the same processing thereafter.
  • step S113 the control device 5 determines whether generation of X-ray projection image data has been performed for all measurement angles. If generation of X-ray projection image data has been performed for all measurement angles, an affirmative determination is made in step S113, and the process proceeds to step S114.
  • step S114 the image reconstruction unit 56 generates three-dimensional data using a plurality of X-ray projection image data, and ends the process. The control device 5 can display this three-dimensional data on a monitor (not shown) or store it in a memory (not shown). If there is a measurement angle at which X-ray projection image data is not generated in step S113, the determination in step S113 is negative, and the process proceeds to step S115.
  • step S115 the mounting table control unit 52 drives the manipulator unit 36 to rotationally drive the mounting table 31 to a predetermined measurement angle, and the process returns to step S102, and the same processing is performed thereafter.
  • the scattering portion 41A of the detector 4A transmits X-rays transmitted through the object S (or X-rays transmitted straight without being irradiated to the object S) or scattering scattered inside the object S It detects information on the progress of electrons generated by the incidence and scattering of X including X-rays.
  • the identifying unit 54 of the control device 5 transmits X-rays based on the information on the scattering position of X-rays incident on the detector 4A and scattered, the incident position of X-rays incident on the absorbing portion 42A, and the progress of electrons.
  • the scattering portion 41A of the detector 4A is a recoil electron generated by incidence X of transmitted X-rays transmitted through the object S or X including scattered X-rays scattered inside the object S Detect the direction of travel
  • the identification unit 54 of the controller 5 determines at least the transmitted X-rays and the scattered X-rays based on the traveling direction of the recoil electrons, the incident position of the X-rays detected by the absorber 42A, and the scattering angle ⁇ of the X-rays. Identify one. Thereby, it becomes possible to sort the X-rays scattered X-rays scattered by the object S and the X-rays not scattered.
  • When generating X-ray projection image data it becomes possible to exclude scattered X-rays that cause noise from the internal information of the object S, and the image quality of the X-ray projection image data can be improved. .
  • the identifying unit 54 calculates the energy E1 lost due to the X-rays being Compton scattered in the scattering unit 41A, and the traveling direction of the recoil electrons, the calculated energy E1 and the energy E2 obtained by the absorbing unit 42A. And at least one of the transmitted X-rays and the scattered X-rays is identified. Thereby, it becomes possible to sort the X-rays scattered X-rays scattered by the object S and the X-rays not scattered.
  • the scattering unit 41 detects the scattering point T1 where X-rays are Compton scattered in the chamber 411A, and the identifying unit 54 performs the Compton scattering in the chamber 411A and the scattering point T1, the initial traveling direction of recoil electrons, and At least one of the transmitted X-ray and the scattered X-ray is specified based on the relationship between the Compton scattering angle ⁇ of the X-ray and the incident position of the X-ray incident on the absorbing portion 42A. Thereby, it is possible to generate X-ray projection image data by excluding scattered X-rays among X-rays incident on the detector 4A.
  • the identifying unit 54 determines the presence of the X-ray source 2 based on the scattering point T1, the initial traveling direction of recoil electrons, the Compton scattering angle ⁇ , and the incident position of X-rays to the absorbing unit 42A.
  • the position to be estimated (pathway L1) is calculated, and when the X-ray source 2 exists at the calculated position, the X-ray is specified as a transmission X-ray.
  • the image quality of the X-ray projection image data can be improved.
  • the specifying unit 54 determines that the time difference between the time Tj + t sig at which the electrical signal corresponding to the energy E1 is output from the scattering unit 41A and the time Tj at which the electrical signal corresponding to the energy E2 is output from the absorbing unit 42A is predetermined.
  • the scattered X-rays are identified using X-rays within the range T1 of
  • the scattered X-rays can be identified by excluding the X-rays emitted from the X-ray source 2 at different timings and incident on the detector 4A, so that the accuracy of identifying the scattered X-rays can be improved.
  • the identifying unit 54 identifies the scattered X-rays by irradiating the X-ray source 2, transmitting different positions of the object S, and excluding the X-rays incident on the same pixel of the absorbing unit 42A. Do. As a result, it is suppressed that the information of different positions inside the object S is included in the internal information as the information of the same position, so it is possible to suppress the deterioration of the image quality of the X-ray projection image data.
  • the specifying unit 54 measures The different positions of the object S are transmitted to detect the plurality of X-ray photons incident on the same position of the absorber 42A, and the detected X-ray photons are excluded to specify the scattered X-ray. Thereby, X-rays transmitted through different positions of the object to be measured S and incident on the same pixel of the absorbing section 42A can be excluded from the internal information of the object to be measured S.
  • the X-ray apparatus 100A of the second embodiment described above can be modified as follows.
  • the scattering unit 41A may be configured by a known three-dimensional quantum image detector using an SOI.
  • FIG. 12A shows the configuration of the detector 4 in this case.
  • a three-dimensional quantum image detector using an SOI includes a Si sensor unit 500 having a thickness along the Z direction, and the Si sensor unit 500 detects a track of recoil electrons.
  • the Si sensor unit 500 has a silicon layer 501 and an integrated circuit 502 on the Z direction + side of the silicon layer 501.
  • the X-rays from the object S are transmitted through the silicon layer 501 or Compton scattered, and the integrated circuit 502 detects the three-dimensional shape of the track of recoil electrons by Compton scattering in the same manner as the embodiment.
  • the scattering unit 41A and the absorbing unit 42A may have a configuration in which a plurality of three-dimensional quantum image detectors are stacked along the Z direction.
  • FIGS. 12 (b) and 12 (c) show the configuration of the detector 4A in this case.
  • FIG. 12 (b) is a cross-sectional view in the YZ plane
  • FIG. 12 (c) is an external view in the XY plane.
  • the scattering unit 41A includes a plurality of three-dimensional quantum image detectors 511a to 511d (generally denoted by reference numeral 511) stacked along the Z direction, and a processing circuit 513 provided on the processing substrate 512. .
  • the absorbing unit 42A includes a plurality of three-dimensional quantum image detectors 511e to 511h (generally denoted by reference numeral 511) stacked along the Z direction, and a processing circuit 515 provided on the processing substrate 514. .
  • the scattering portion 41A and the absorption portion 42A each have four three-dimensional quantum image detectors 511a to 511d and 511e to 511h, respectively.
  • the number of three-dimensional quantum image detectors 511 included in the unit 41A and the absorbing unit 42A is not limited to the examples illustrated in FIGS. 12 (b) and 12 (c).
  • the scattering unit 41A and the absorbing unit 42A may have different numbers of three-dimensional quantum image detectors 511.
  • the three-dimensional quantum image detector 511 includes a substrate 520 and a semiconductor element 521.
  • the semiconductor element 521 is formed in a thin plate shape, and is provided on the substrate 520.
  • the semiconductor element 521 has a plurality of strip electrodes on the top surface (Z direction side).
  • the strip electrodes are connected by wires 523 to processing circuits 513 and 515 provided on the processing substrates 512 and 514, respectively.
  • the strip electrode outputs an electrical signal when X-rays or recoil electrons are incident.
  • the X-ray apparatus 100A described in the second embodiment and the modification described above can be applied to the structure manufacturing system 400 described with reference to FIG. 5 in the first embodiment.
  • the structure manufacturing system 400 has the X-ray apparatus 100A instead of the X-ray apparatus 100 shown in FIG. 5, and the X-ray apparatus 100A is created by the forming apparatus 420 based on the design process of the design apparatus 410.
  • the measurement processing for acquiring the shape information of the structure is performed, and the inspection unit 432 of the control system 430 performs inspection processing for comparing the shape information acquired in the measurement processing with the design information generated in the design processing.
  • the same function and effect as the function and effect obtained by the structure manufacturing system 400 of the first embodiment can be obtained.
  • the third embodiment has the same configuration as that of the X-ray apparatus 100A of the second embodiment shown in FIG. That is, the X-ray apparatus 100A of the third embodiment is the same as the first embodiment shown in the casing 1, the X-ray source 2, the placement unit 3 and the control device 5 shown in FIG. It has a detector 4A of the embodiment.
  • the specifying unit 54 of the control device 5 is the same as the second embodiment, based on the electrical signals output from the scattering unit 41A and the absorbing unit 42A, from the X-ray source 2 It is determined whether X-rays have been scattered within the object S.
  • the identifying unit 54 identifies, among X-rays from the object S, scattered X-rays scattered at the object S or identifies X-rays other than the scattered X-rays, based on the determination result. That is, the identifying unit 54 identifies whether the X-ray incident on the detector 4A is the X-ray scattered by the object S.
  • the identifying unit 54 When the identifying unit 54 identifies that the X-ray incident on the detector 4A is a scattered X-ray, the identifying unit 54 identifies a position at which the scattered X-ray is scattered in the object S.
  • a method of specifying the position (scattering position) at which the X-ray is scattered in the object to be measured S will be described in detail.
  • the identifying unit 54 receives the energy E1 that the X-ray has lost due to Compton scattering in the scattering unit 41A of the detector 4A and the X that has reached the absorbing unit 42A after Compton scattering.
  • the energy E2 of the line is added to calculate the energy E0 of the X-ray incident on the chamber 411A.
  • the identifying unit 54 calculates the scattering angle ⁇ and the azimuth angle ⁇ by using the above-described equation (8) or (9) to determine the path L1 of the X-ray incident on the chamber 411A (see FIG. 9). .
  • the specifying unit 54 determines a position at which the X-ray is scattered inside the object S (hereinafter, referred to as a scattering position).
  • a scattering position a position at which the X-ray is scattered inside the object S
  • the method of determining the scattering position when the X-ray emitted from the X-ray source 2 is monochromatic and when the X-ray emitted from the X-ray source 2 has an energy distribution and the maximum energy is unknown, X-ray The X-rays emitted from the source 2 have an energy distribution, and the maximum energy is known.
  • Es The energy of monochromatic X-rays emitted from the X-ray source 2 will be described as Es.
  • FIG. 13 shows a case where electrical signals are output from the pixels 412a and 412f of the scattering portion 41A and the pixel 42a of the absorbing portion 42A, as in FIG. 9B.
  • the path L1 of X-rays characterized based on the outputs of the scattering portion 41A and the absorbing portion 42A and the path L11 of X-rays from the X-ray source 2 intersect at the scattering angle ⁇
  • the point P1 to be determined is determined as the scattering position.
  • the identifying unit 54 determines the coordinates of the absorption point T2 (the position of the pixel that outputs the energy E2 in the absorbing unit 42A) and the coordinates of the scattering point T1 (the position of the pixel that outputs the electrical signal at the latest time of the pixel electrode 412A).
  • the coordinates of the scattering position P1 in the XYZ coordinate system are calculated based on the position determined by the position in the Z direction and the scattering angle ⁇ of the scattering section 41A, the position of the X-ray source 2, and the scattering angle ⁇ . .
  • the identifying unit 54 calculates the coordinates of the plurality of scattering positions P by performing the above-described process on the plurality of scattered X-rays identified as having been scattered. As a result, the frequency distribution of the scattering positions is obtained inside the object S. By this frequency distribution, it is possible to know the internal information of the object S, for example, the distribution state of the material having a large X-ray scattering rate inside the object S, and the like. Further, the specifying unit 54 can generate image data in which each of the calculated plurality of scattering positions P is superimposed on the image data corresponding to the object S. Thereby, it is possible to provide an image that visually represents the internal information of the object S.
  • the scattering angle ⁇ and the scattering position P1 obtained by the above equation (13) are closest to the detector 4A among multiple scattering.
  • the scattering angle and the scattering position when scattered on the side (Z direction + side) or the different scattering angle and the scattering position may be indicated.
  • the identification unit 54 performs the same process as described above each time the object S is rotated by a predetermined rotation amount, and indicates the distribution state of scattering inside the object S.
  • Image data may be generated. The data value of the scattering position different from the actual one is averaged by the multiple scattering by adding these plural image data to generate the composite image data, and the representation accuracy of the distribution state of the scattering inside the object S Deterioration of the
  • FIG. 14 schematically shows the positional relationship between the X-ray source 2, the detector 4A and the object S in the XZ plane, and FIG.
  • the measurement object S includes, as an example, a region S1 made of a substance having different scattering characteristics (scattering rate) inside.
  • the scattering characteristics of the object S are illustrated as known, but in actuality, both the size of the object S and the scattering characteristics are obtained by the following procedure.
  • FIG. 14A X-rays incident on incident points (scattering points T1-1, T1-2,..., T1-N) are scattered X-rays scattered by the object S.
  • the identifying unit 54 calculates the path L1-1 of the scattered X-ray incident on the incident point T1-1 using the above-described equation (8).
  • the calculated path L1-1 is represented by a straight line, and the scattering position of the object S exists at any position on the path L1-1.
  • This path L1-1 is referred to as a scattering point candidate position U1-1. That is, the scattering position on the object S is any point existing on the scattering point candidate position U1-1.
  • the paths of two other scattered X-rays incident in the vicinity of the incident point T 1-1 are illustrated in FIG. 14A by dashed dotted lines. These are paths of X-rays scattered at another scattering position different from the above scattering position in the object S.
  • respective scattering point candidate positions are determined.
  • the identifying unit 54 calculates the path L1-2 of the scattered X-ray incident on the incident point T1-2 using Expression (8).
  • the scattering position of the object S exists at any position on the calculated path L1-2, as in the case of the path L1-1.
  • the scattering position of the object S is any point existing on the scattering point candidate position U1-2.
  • the paths of two other scattered X-rays incident in the vicinity of the incident point T1-2 are illustrated in FIG. 14 (a) by dashed dotted lines.
  • the scattered X-ray is performed by the same procedure as described above in a state where the object S is rotated by a predetermined rotation amount 1 1 (for example, 90 ° clockwise). obtaining scattered point candidate distribution D 1 which path is present. That is, the specifying unit 54 specifies a scattering point candidate position U1-1 ′, a scattering point candidate position U1-2 ′,..., A scattering point candidate position U1-N ′. The specifying unit 54, by superimposing their scattering point candidate position, in the state of FIG. 14 (b), the scattering point candidate distribution D 1 of the scattered X-rays in the object to be measured S is obtained.
  • a predetermined rotation amount 1 1 for example, 90 ° clockwise
  • certain compounds 54, 2 measurement object S ⁇ , ⁇ 3 ⁇ ⁇ N plurality of scattering points in the rotated state candidate distribution D 2, D 3 Request ⁇ ⁇ ⁇ D N.
  • the intersection of these scattering point candidate distributions D 0 to D N is determined.
  • This common part corresponds to the distribution state of scattering inside the object S (for example, the distribution region of the material having a large scattering rate inside the object S).
  • the common part can be determined, for example, by the procedure described below.
  • Scattering point candidate distributions D 1 , D 2 , D 3 ... D N are scattered point candidate distributions D 1 ′, D 2 ′, D 3 ′, ⁇ rotated by the angle obtained by rotating the object S. ⁇ ⁇ Find D N '.
  • the scattering point candidate distributions D 0 , D 1 ′, D 2 ′, D 3 ′,..., D N ′ are added, and the surrounding halo portion is subtracted and removed. It is calculated.
  • the calculated common part is taken as the distribution of scattering within the object S.
  • the specific object 54 can generate image data in which the calculated common part is superimposed on the image data corresponding to the object S. Thereby, it is possible to provide an image that visually represents the internal information of the object S.
  • the X-rays emitted from the X-ray source 2 have an energy distribution. Since the X-ray source emits X-rays by bremsstrahlung, the energy Ei of the emitted X-ray is smaller than the maximum energy Esmax, and the maximum energy Esmax is known.
  • the scattering angle ⁇ when X rays of energy Ei are scattered in the object S is expressed by the following equation (14).
  • the energy E0 can be calculated based on the output of the detector 4A. However, since the energy Ei is unknown, the scattering angle ⁇ can not be determined using equation (14).
  • the scattering angle when the X-ray of the maximum energy Esmax is scattered by the object S by replacing Es in the above equation (13) with Esmax. ⁇ can be calculated. Since the energy Ei is smaller than the maximum energy Esmax, the scattering angle ⁇ of the X-ray of the energy Ei at bremsstrahlung is smaller than the scattering angle ⁇ of the X-ray at the maximum energy Esmax.
  • FIG. 15 shows a case where the X-ray emitted from the X-ray source 2 with the maximum energy Esmax and the X-ray emitted with the energy Ei at the time of bremsstrahlung enter the scattering portion 41A along the path L1.
  • X-rays emitted from the X-ray source 2 with the maximum energy Esmax are scattered at the scattering angle ⁇ at the scattering position P1 of the object S and enter the scattering point (incident position) T1 of the scattering portion 41A at the path L1.
  • X-rays emitted from the X-ray source 2 with energy Ei at the time of bremsstrahlung are scattered at the scattering angle ⁇ at the scattering position P2 of the object S and scatter point (incident position) T1 of the scattering portion 41A at the path L1. Incident to Therefore, it can be estimated that the scattering position P2 is located on a straight line connecting the scattering point (incident position) T1 and the scattering position P1.
  • FIG. 16 is a view schematically showing the positional relationship between the X-ray source 2, the detector 4A, and the object S in the XZ plane as in FIG. 14, and FIG. 16 (b) is a diagram of FIG. The case where the measured object S is rotated by a predetermined rotation amount (90 ° clockwise) from the state of the positional relationship shown in FIG.
  • FIG. 16A X-rays incident on incident positions (scattering points T1-1, T1-2,..., T1-N) are scattered X-rays scattered by the object S.
  • the scattering characteristics of the object S are illustrated as known, but in actuality, both the size of the object S and the scattering characteristics are determined by the following procedure.
  • the identifying unit 54 calculates the path L1-1 of the scattered X-ray incident on the incident point T1-1 using the equation (8).
  • the calculated path L1-1 is represented by a straight line, and the scattering position of the object S exists at any position on the path L1-1. That is, when the path L1-1 is referred to as the scattering point candidate position U1-1, the scattering position on the measuring object S is any point existing on the scattering point candidate position U1-1.
  • the X-ray scattering position P2 of the energy Ei in the measured object S is located in the range from the incident position T1-1 to the scattering position P1 on the path L1-1. it can.
  • the scattering point candidate position U1-1 does not include the region closer to the X-ray source 2 (Z direction ⁇ side) than the object S.
  • the paths of two other scattered X-rays incident in the vicinity of the incident point T 1-1 are illustrated in FIG. 16A by dashed dotted lines. These are paths of X-rays scattered at another scattering position different from the above scattering position in the object S. Also with respect to the paths of these scattered X-rays, respective scattering point candidate positions are determined.
  • the identification unit 54 calculates the path L1-2 of the scattered X-ray incident on the incident point T1-2 using the equation (8).
  • the scattering position of the object S exists at any position on the calculated path L1-2 (ie, the scattering point candidate position U1-2).
  • the scattering point candidate position U1-2 does not include the region closer to the X-ray source 2 (Z direction ⁇ side) than the object S.
  • the paths of two other scattered X-rays incident in the vicinity of the incident point T1-2 are illustrated in FIG. 16 (a) by dashed dotted lines. These are paths of X-rays scattered at another scattering position different from the above scattering position in the object S. Also with respect to the paths of these scattered X-rays, respective scattering point candidate positions are determined.
  • the scattering point candidate positions U1-N are specified for the plurality of incident points T1-N.
  • the specifying unit 54 calculates the coordinates of the scattering point candidate position U1-1 ′, the scattering point candidate position U1-2 ′,..., The scattering point candidate position U1-N ′ on the image data.
  • the course distribution D 1 is obtained.
  • certain compounds 54, 2 measurement object S ⁇ , ⁇ 3 ⁇ ⁇ N plurality of scattering points in the rotated state candidate distribution D 2, D 3 Request ⁇ ⁇ ⁇ D N.
  • the intersection of these scattering point candidate distributions D 0 to D N is determined.
  • This common part corresponds to the distribution state of scattering inside the object S (for example, the distribution region of the material having a large scattering rate inside the object S).
  • the common part can be determined, for example, by the procedure described below.
  • Scattering point candidate distributions D 1 , D 2 , D 3 ... D N are scattered point candidate distributions D 1 ′, D 2 ′, D 3 obtained by rotating the object S by an angle obtained by rotating the object S ' ⁇ Ask for D N '.
  • the scattering point candidate distributions D 0 , D 1 ′, D 2 ′, D 3 ′,..., D N ′ are added, and the surrounding halo portion is subtracted and removed. It is calculated.
  • the intensity of the peripheral halo portion is much smaller than the case where the maximum energy of the X-ray from the X-ray source 2 described with reference to FIG. 14 is unknown.
  • the calculated common part is taken as the distribution of scattering within the object S.
  • the specific object 54 can generate image data in which the calculated common part is superimposed on the image data corresponding to the object S. Thereby, it is possible to provide an image that visually represents the internal information of the object S.
  • FIG. 14 and FIG. 16 an example in which the object S is rotated is shown, but instead, the detector 4A and the X-ray source 2 are rotated around the object S It is also good.
  • the X-ray apparatus 100A of the third embodiment described above identifies the scattered X-rays scattered inside the object S, and represents the internal information and the internal information of the object S based on the scattered X-rays. Image data can be generated.
  • the X-ray apparatus 100A described in the third embodiment described above can be applied to the structure manufacturing system 400 described with reference to FIG. 5 in the first embodiment.
  • the structure manufacturing system 400 has the X-ray apparatus 100A instead of the X-ray apparatus 100 shown in FIG. 5, and the X-ray apparatus 100A is created by the forming apparatus 420 based on the design process of the design apparatus 410.
  • the measurement processing for acquiring the shape information of the structure is performed, and the inspection unit 432 of the control system 430 performs inspection processing for comparing the shape information acquired in the measurement processing with the design information generated in the design processing.
  • the same function and effect as the function and effect obtained by the structure manufacturing system 400 of the first embodiment can be obtained.
  • the mounting table 31 on which the object S is placed is the X axis and Y axis by the X axis moving unit 33, the Y axis moving unit 34, and the Z axis moving unit 35. And is not limited to one that is moved in the Z-axis direction.
  • the mounting table 31 does not move in the X-axis, Y-axis and Z-axis directions, and moves the X-ray source 2 and the detector 4 in the X-axis, Y-axis and Z-axis directions, thereby
  • the X-ray apparatus 100 may have a configuration for moving the radiation source 2 and the detector 4 relative to each other. Further, the X-ray apparatus 100 has a configuration in which the X-ray source 2 and the detector 4 rotate on the rotation axis Yr, instead of the mounting table 31 rotating on the rotation axis Yr, and the mounting table 31 does not rotate. You may have.
  • control device 5 has image generation part 53 and image reconstruction. It is not necessary to have the part 56.
  • the image generation unit 53 and the image reconstruction unit 56 are provided in a processing device or the like separate from the X-ray apparatus 100, and information on X-rays other than the scattered X-ray or the scattered X-ray identified by the identifying unit 54
  • X-ray projection image data or three-dimensional data may be generated by acquiring via, for example, a network or a storage medium.

Landscapes

  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

X線装置は、被測定物を透過した透過X線または被測定物の内部で散乱した散乱X線を含むX線が入射され、入射したX線が散乱した散乱位置を検出する第1検出器と、第1検出器で散乱または透過したX線が入射され、入射したX線の入射位置を検出する第2検出器と、散乱位置および入射位置に基づいて、透過X線および散乱X線の少なくとも一方を特定する特定部と、を備えるX線装置。

Description

X線装置および構造物の製造方法
 本発明は、X線装置および構造物の製造方法に関する。
 従来から、入射した放射線(X線またはγ線)のコンプトン散乱を検出して、その結果に基づいて放射線源の画像を生成するコンプトンカメラが知られている(例えば特許文献1)。しかしながら、コンプトン散乱を検出し、その結果に基づいて被測定物の画像を生成する非破壊検査装置については知られていない。
日本国特許第5991519号
 第1の態様によると、X線装置は、被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射され、入射した前記X線が散乱した散乱位置を検出する第1検出器と、前記第1検出器で散乱または透過した前記X線が入射され、入射した前記X線の入射位置を検出する第2検出器と、前記散乱位置および前記入射位置に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、を備える。
 第2の態様によると、X線装置は、被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射して散乱することにより失う第1エネルギーを検出する第1検出器と、前記第1検出器で散乱または透過した前記X線の第2エネルギーを検出する第2検出器と、前記第1エネルギーおよび前記第2エネルギーに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、を備える。
 第3の態様によると、X線装置は、被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射した後に散乱を起こすことにより生じる電子の進行に関する情報と、前記散乱を起こした後の前記X線の進行に関する情報とを検出する検出器と、前記電子の進行に関する情報と、前記X線の進行に関する情報とに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、を備える。
 第4の態様によると、X線装置は、被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線を検出する検出器と、前記検出器への前記X線の入射方向を推定し、推定した前記入射方向とX線源の位置とに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、を備える。
 第5の態様によると、構造物の製造方法は、構造物の形状に関する設計情報を作成し、 前記設計情報に基づいて前記構造物を作成し、作成された前記構造物の形状を、第1から第4のいずれか一つの態様のX線装置を用いて計測して形状情報を取得し、前記取得された前記形状情報と前記設計情報とを比較する。
第1の実施の形態によるX線装置の構成を模式的に示すブロック図である。 検出器の散乱部と吸収部とに入射するX線の経路を模式的に示す図である。 被測定物の異なる位置を透過したX線が吸収部の同一の位置に入射する場合を模式的に示す図である。 第1の実施の形態によるX線装置が実行する処理を説明するフローチャートである。 実施の形態による構造物製造システムの構成を模式的に示すブロック図である。 実施の形態による構造物製造システムが実行する処理を説明するフローチャートである。 第2の実施の形態によるX線装置の構成を模式的に示すブロック図である。 検出器の散乱部と吸収部とに入射するX線の経路を模式的に示す図である。 検出器の散乱部と吸収部とに入射するX線の経路を模式的に示す図である。 被測定物の異なる位置を透過したX線が吸収部の同一の位置に入射する場合を模式的に示す図である。 第2の実施の形態によるX線装置が実行する処理を説明するフローチャートである。 変形例における検出器の構成を模式的に示す図である。 第3の実施の形態による検出器の散乱部と吸収部とに入射するX線の経路を模式的に示す図の一例である。 第3の実施の形態において、X線の最大エネルギーが未知の場合の例を説明するための、X線源と検出器と被測定物との位置関係を模式的に示す図である。 第3の実施の形態において、X線の最大エネルギーが既知の場合の例を説明するための、検出器の散乱部と吸収部とに入射するX線の経路を模式的に示す図である。 第3の実施の形態において、X線の最大エネルギーが既知の場合の例を説明するための、X線源と検出器と被測定物との位置関係を模式的に示す図である。
-第1の実施の形態-
 図面を参照しながら、第1の実施の形態によるX線装置について説明する。X線装置は、被測定物にX線を照射して、被測定物を透過したX線を検出することにより、被測定物の内部情報(例えば内部構造)等を被測定物を破壊することなく取得する。機械部品や電子部品等の産業用部品を対象とするX線装置は、産業用X線検査装置と呼ばれる。
 図1は本実施の形態によるX線装置100の構成の一例を示す図である。なお、説明の都合上、X軸、Y軸、Z軸からなる座標系を図示の通りに設定する。
 X線装置100は、筐体1、X線源2、載置部3、検出器4および制御装置5を備えている。筐体1は、その下面が工場等の床面に実質的に平行(水平)となるように配置される。筐体1の内部には、X線源2と、載置部3と、検出器4とが収容される。筐体1は、X線が筐体1の外部に漏洩しないようにするため、X線遮蔽材料を含む。なお、X線遮蔽材料として鉛を含む。
 X線源2は、制御装置5による制御に応じて、図1に示す出射点Pを頂点としてZ軸に平行な光軸Zrに沿って、Z軸+方向へ向けてX線を放射する。この出射点Pは後述するX線源2の内部を伝搬する電子線の焦点位置と一致する。すなわち、光軸Zrは、X線源2の電子線の焦点位置である出射点Pと、後述する検出器4の撮像領域の中心とを結ぶ軸である。なお、X線源2から放射するX線は、円錐状に拡がるX線(いわゆるコーンビーム)、扇状のX線(いわゆるファンビーム)、および直線状のX線(いあわゆるペンシルビーム)のいずれでもよい。なお、ファンビームおよびペンシルビームを用いる場合は、被測定物S全体を検査するために、ビームと被測定物Sとを相対的に移動させるスキャン動作を行う必要がある。X線源2からは、短い時間スケールで見た場合、X線光子が1個ずつ放出される。X線源2は、例えば約50eVの超軟X線、約0.1~2keVの軟X線、約2~20keVのX線および約20~数MeVの硬X線の少なくとも1つを照射する。なお、散乱が顕著なのは、硬X線の領域である。
 載置部3は、被測定物Sが載置される載置台31と、回転駆動部32、X軸移動部33、Y軸移動部34およびZ軸移動部35からなるマニピュレータ部36とを備え、X線源2よりもZ軸+側に設けられている。載置台31は、回転駆動部32により回転可能に設けられる。後述するように、回転駆動部32による回転軸YrがX軸、Y軸、Z軸方向に移動する際に、載置台31はともに移動する。回転駆動部32は、例えば電動モータ等によって構成され、後述する制御装置5により制御されて駆動した電動モータが発生する回転力によって、載置台31を回転させる。載置台31の回転軸Yrは、Y軸に平行、かつ、載置台31の中心を通過する。X軸移動部33、Y軸移動部34およびZ軸移動部35は、制御装置5により制御されて、載置台31をX軸方向、Y軸方向およびZ軸方向にそれぞれ移動させる。Z軸移動部35は、制御装置5により制御されて、X線源2から被測定物Sまでの距離が、撮影される画像における被測定物Sの拡大率に応じた距離となるように載置台31をZ軸方向に移動させる。
 図1に示す検出器4は、載置台31よりもZ軸+側に設けられている。すなわち、載置台31は、Z軸方向において、X線源2と検出器4との間に設けられる。検出器4は、散乱部41と吸収部42とを備え、被測定物Sを透過した透過X線(または被測定物Sに照射されずに直進した透過X線)および散乱した散乱X線の少なくとも一方を検出する。散乱部41は、吸収部42よりもZ方向-側、すなわち被測定物Sに近い側に設けられる。散乱部41のZ方向の厚さは、吸収部42のZ方向の厚さよりも小さい。
 散乱部41は、XY平面に平行な入射面411を有する。X線源2から放射され、載置台31上に載置された被測定物Sを透過した、または被測定物S内部で散乱したX線は入射面411に入射する。なお、以後の説明においては、被測定物Sの内部で散乱せずに内部を透過して検出器4に入射するX線を透過X線、被測定物Sの内部で散乱して検出器4に入射するX線を散乱X線と呼ぶ。散乱部41は、例えば公知のシンチレーション物質を含むシンチレータ部と、シンチレータ部から放出された光を増幅して受光する受光部とにより構成される。散乱部41の入射面411に入射したX線は、シンチレータ部において、ある確率でコンプトン散乱する。コンプトン散乱した場合、その反跳電子によりシンチレータ部からは蛍光が放出される。この蛍光は、受光部において光電子増倍管により増幅されて電気信号として制御装置5へ出力される。シンチレータ部においてコンプトン散乱したX線は、進行方向が変化した後、散乱部41を出射し、Z方向+側に設けられた吸収部42に入射する。一方、シンチレータ部においてコンプトン散乱しなかったX線は、散乱部41において進行方向が変化することなく散乱部41を出射し、吸収部42に入射する。
 吸収部42は、XY平面に平行な入射面421を有する。散乱部41から到達したX線は入射面421に入射する。吸収部42は、公知のシンチレーション物質を含むシンチレータ部と、シンチレータ部から放出された光を増幅して受光する受光部等とによって構成される。シンチレータ部の入射面421に入射したX線は吸収部42のシンチレータ部に吸収され蛍光を放出し、放出された蛍光の光エネルギーは上記の受光部で光増倍管により増幅されて電気エネルギーに変換され、電気信号として制御装置5へ出力される。
 検出器4の散乱部41と吸収部42とは共に、シンチレータ部と受光部とが複数の画素が二次元的に配置された構造を有し、シンチレータ部と受光部との複数の画素のそれぞれは、互いに対応するように配列されている。これにより、X線源2から放射され、被測定物Sを透過、散乱または一部が吸収されたX線の強度データを、複数の画素において一括して取得することができ、それぞれの強度データを統合することにより強度分布データを得ることができる。
 制御装置5は、マイクロプロセッサやその周辺回路等を有しており、不図示の記憶媒体(例えばフラッシュメモリ等)に予め記憶されている制御プログラムを読み込んで実行することにより、X線装置100の各部を制御する。制御装置5は、X線源2の動作を制御するX線制御部51、マニピュレータ部36の駆動動作を制御する載置台制御部52、検出器4から出力された電気信号に基づいて被測定物SのX線投影画像データを生成する画像生成部53、および特定部54を有する。特定部54は、散乱部41および吸収部42から出力された電気信号に基づいて、X線源2からのX線が被測定物Sの内部で散乱したか否かを判定する。特定部54は、判定結果に基づいて、被測定物SからのX線のうち、被測定物Sにおいて散乱した散乱X線を特定、または散乱X線以外のX線を特定する。すなわち、特定部54は、検出器4に入射したX線が、被測定物Sで散乱したX線であるか否かを特定する。画像生成部53は、特定部54により特定された散乱X線以外のX線、すなわち、被測定物Sにおいて散乱しなかった透過X線や、一部が吸収されたX線による強度分布データを、X線投影画像データを生成する際の被測定物Sの内部情報として選択する。画像再構成部56は、被測定物Sに対するX線照射方向を相対的に変化させて投影し、それにより得られた複数のX線強度分布データに基づいて、公知の画像再構成処理方法を用いることで、被測定物Sの再構成画像を生成する。画像再構成処理により、被測定物Sの内部構造(断面構造)である断面画像データや3次元データが生成される。なお、断面画像データとは、XZ平面と平行な面内における被測定物Sの構造データを含む。画像再構成処理としては、逆投影法、フィルタ補正逆投影法、逐次近似法等がある。
 以下、特定部54による、被測定物SにおいてX線が散乱したか否かを判定し、透過X線および散乱X線の少なくとも一方を特定する手順について、図2を参照して説明する。図2(a)は、ある時点において、散乱部41の一つの画素である画素41aと吸収部42の一つの画素である画素42aの両方において、電気信号が出力された場合を示している。散乱部41は、画素41aの位置により、入射したX線が検出器4内で散乱した散乱位置を検出し、吸収部42は、画素42aの位置によりX線が入射した入射位置を検出する。この場合、画素41aから電気信号が出力されているので、画素41aに入射したX線は、画素41a内の散乱点T1においてコンプトン散乱し、それにより進行方向が変化して画素42aに入射したと判定される。画素41aから出力された電気信号および画素42aから出力された電気信号に基づいて、特定部54は、それぞれエネルギーE1とエネルギーE2を算出する。エネルギーE1は、画素41aにおけるX線のコンプトン散乱に伴って失われたエネルギー、また、エネルギーE2は、画素41aにおいてコンプトン散乱した後に画素42aに到達したX線のエネルギーに相当する。特定部54は、画素41aにおけるコンプトン散乱によるX線の進行方向の変化角度、すなわちコンプトン散乱角をθ1とすると、コンプトン散乱角θ1は、次の式(1)により算出できる。
Figure JPOXMLDOC01-appb-M000001
 ここで、電子の静止質量me=9.1×10-31(kg)、真空中の光速c=3.0×10(m)として、m=511(keV)である。
 画素41aと画素42aとの位置関係は既知なので、画素41aを出射して画素42aに入射したX線の進路L2を知ることができる。画素41aに入射したX線の進路を進路L1とすると、進路L1と進路L2とは散乱点T1においてコンプトン散乱角θ1をなして交差する。従って、進路L1は、散乱点T1を頂点とし進路L2を中心軸として、中心軸であるL2と母線とのなす角がθ1の円錐面CS上を散乱点T1に向かう直線のいずれかに相当することがわかる。
 次に、検出されたX線の入射位置、円錐面とX線源2との位置とに基づいて透過X線および散乱X線の少なくとも一方を特定する手順、すなわち、上記円錐面CSとX線源2との位置関係に基づいて、図2(a)を参照した上記説明の場合に、X線源2から放射されたX線が被測定物Sで散乱したか否かを判定する手順について説明する。なお、X線装置100においては、図1を用いて上述したようにX線源2と検出器4とが配置されるので、検出器4の位置に対するX線源2の位置は既知である。
 図2(b)に、円錐面CS上にX線源2が存在する場合を示す。この場合には、X線源2が、円錐面CS上を散乱点T1に向かう直線のうちの1つであるL1の上に存在することになる。このことから、特定部54は、X線源2から放射されたX線は、被測定物Sの内部において進行方向を変えることなしに、画素41aに入射し、その進路はL1であると判定できる。なお、円錐面CS上にX線源2が存在しない場合であっても、X線源2と円錐面CSとの距離が所定の長さより小さい場合には、特定部54は、上記と同様に、X線源2から放射されたX線は、被測定物Sの内部で散乱することなしに進路L1に沿って進み、画素41aに入射したと判定できる。すなわち、特定部54は、検出器4に入射したX線が、被測定物Sの内部で散乱せずに透過した透過X線であると特定する。上記の所定の長さは、画素41aおよび画素42aの大きさや、X線装置100における、X線源2、被測定物S、検出器4等の配置等に基づいて、適宜設定することが好ましい。
 次に、図2(c)に、円錐面CS上にX線源2が存在せず、かつ、X線源2と円錐面CSとの距離が所定の長さより大きい場合を示す。この場合には、X線源2が、円錐面CS上を散乱点T1に向かう直線のうちの1つであるL1の上に存在しない。このことから、特定部54は、X線源2から放射されたX線は、被測定物Sの内部において散乱して進行方向を変えて画素41aに入射したと判定できる。すなわち、特定部54は、X線源2から放射されて進路L0を進んだX線は、被測定物Sに入射し、その内部の散乱点T3において散乱して進行方向がL1に変化し、画素41aに入射した散乱X線であると特定する。
 上述した通り、円錐面CSの上、または、円錐面CSの近傍にX線源2が位置すると判定された場合には、特定部54は、X線が被測定物Sで散乱していないと判定し、散乱X線以外のX線であることを特定する。一方、円錐面CSの上および円錐面CSの近傍にX線源2が位置していない場合には、特定部54は、X線が被測定物Sで散乱したと判定し、散乱X線であることを特定する。
 次に、特定部54が行う、被測定物SでX線が散乱したか否かを特定し、散乱X線と散乱X線以外のX線とを特定するための処理(以後、特定処理と呼ぶ)について説明する。
 本実施の形態のX線装置100においては、コンプトン散乱角θ1の検出精度Δθ1は次のようにして求まる。すなわち、散乱部41に入射するX線のエネルギーETを一定と仮定する。エネルギーETをエネルギーE1とエネルギーE2との和とした場合、検出精度Δθ1は、式(1)を微分して求めた次の式(2)の関係を満たす。
Figure JPOXMLDOC01-appb-M000002
 この式(2)と上述した式(1)とから、検出精度Δθ1は、次の式(3)により算出される。
Figure JPOXMLDOC01-appb-M000003
なお、
Figure JPOXMLDOC01-appb-M000004
であるものとする。
 例えば、X線源2から放射されて散乱部41に入射するX線のエネルギーを511keV、エネルギー分解能を3パーセント、コンプトン散乱角θ1が20°と仮定した場合、検出精度Δθ1は5.3°となる。また、エネルギー分解能を1パーセントとした場合には、検出精度Δθ1は1.8°となる。なお、散乱部41と吸収部42との間のZ方向の距離を20mm、検出器4の各画素のサイズを200μm×200μmと仮定した場合、散乱角θ1の分解能は10mradとなる。
 特定部54は、検出器4の散乱部41により得られたエネルギーE1、エネルギーE1の算出に用いた電気信号を出力した画素の位置(XSi、YSi、ZSi)、およびエネルギーE1の算出に用いた電気信号が散乱部41から出力された時刻Tiの各情報を取得する。この画素の位置(XSi、YSi、ZSi)は、図2(a)、図2(b)、および図2(c)のそれぞれにおける散乱点T1に相当する。なお、時刻Tiは、X線が散乱部41に入射してコンプトン散乱した時刻に相当する。
 特定部54は、検出器4の吸収部42により得られたエネルギーE2、エネルギーE2の算出に用いた電気信号を出力した画素の位置(XAj、YAj、ZAj)、およびエネルギーE2の算出に用いた電気信号が吸収部42から出力された時刻Tjの各情報を取得する。この画素の位置(XAj、YAj、ZAj)は、図2(a)、図2(b)、および図2(c)のそれぞれにおける吸収点T2に相当する。なお、時刻Tjは、X線が吸収部42に入射した時刻に相当する。
 特定部54は、エネルギーE1の算出に用いた電気信号が出力された時刻Tiと、エネルギーE2の算出に用いた電気信号が出力された時刻Tjとに基づいて、散乱部41に入射したX線と、吸収部42に入射したX線とが、X線源2からあるタイミングで放射されたX線であるのか否かを判定する。具体的には、特定部54は、時刻TiとTjとの時間差(すなわちTj-Ti)が所定時間T1以内のときに、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とは、同一のX線光子について出力されたものであると判定する。所定時間T1は、X線が散乱部41から吸収部42まで到達するのに必要な時間に基づいて決定される。例えば、所定時間T1は、計算あるいは各種の試験やシミュレーション等により決定され、予めメモリ(不図示)に記憶されている。時刻TiとTjとの時間差であるTijが所定時間T1を超える場合には、特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とは、同一のX線光子について出力されたものではないと判定する。時刻TiとTjとの時間差が所定時間T1を超える場合には、特定部54は、散乱部41から出力されたエネルギーE1と吸収部42から出力されたエネルギーE2とを、散乱X線を特定するための処理に用いないようにする。ここで、前記のように、散乱部41と吸収部42との間のZ方向の距離を20mmとおけば、時刻TiとTjとの時間差Tijは0.1nsec程度の量である。なお、散乱部41と吸収部42との間の距離は、検出器4の各画素サイズと、上述した式(3)で表される検出精度Δθ1とに基づいて、好ましい値に決定することができる。
 時刻TiとTjとの時間差Tijが所定時間T1以内と判定されると、特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とに基づいて、コンプトン散乱角θ1を算出すべきかどうか判定する。コンプトン散乱角θ1を算出すべきでない場合について、次に説明する。
 検出器4の読み出し回路は、所定の時間ごとに、散乱部41および吸収部42からそれぞれエネルギーE1およびエネルギーE2の算出に用いる電気信号を出力する。本明細書においては、上記所定の時間をエネルギー出力間隔Tと呼ぶ。本実施の形態においては、シンチレータの蛍光発光時間を数十nsecとした場合、T=100nsecである。上記の所定時間T1は、エネルギー出力間隔T以下に設定する必要がある。すなわち、
 T1≦T
を満たすように設定する。これにより、エネルギー出力間隔T内に、散乱部41および吸収部42からそれぞれエネルギーE1およびエネルギーE2の算出に用いる電気信号は出力されるようにできる。
 X線源2からは短い時間スケールで見た場合、X線光子が1個ずつ放射される。X線源2からX線光子を放射する平均放射間隔Teと検出器4からエネルギーの算出に用いる電気信号を出力するエネルギー出力間隔Tとの関係について説明する。X線源2から1回の放射により放射されたX線光子が、画素41aおよび画素42aにそれぞれ入射し、検出器4の散乱部41および吸収部42の全ての画素からエネルギーの算出に用いる電気信号の出力を完了するまでを、1イベントと呼ぶ。すなわち、X線源2から放射された単一X線光子に基づく1フレーム分のエネルギー情報を検出器4から得ることを1イベントと呼ぶ。散乱部41におけるコンプトン散乱角θ1を正しく求めるためには、1イベントの間に、X線源2から放射された単一のX線光子の検出値のみを後述する処理に用いるようにする。2回以上の放射によるX線光子が、同時に検出器4に入射した場合には、コンプトン散乱角θ1が正しく求められない可能性があるからである。これをイベントの分離と呼ぶ。従って、放射間隔Teは、1フレーム分のエネルギー情報を得るために要する時間であるエネルギー出力間隔Tより大きいことが必要である。すなわち、
 T≦Te
を満たす必要がある。
 イベントの分離ができずに、コンプトン散乱角θ1が正しく求められない場合について、次に説明する。図3は、X線源2からの放射間隔Teがエネルギー出力間隔TOより小さい場合について示す。なお、図3においては、説明を簡単にするために、X線源2から2つのX線光子が異なる方向にほぼ同時に放射された場合を示す。1つのX線光子(第1のX線光子)は、X線源2から進路L11を進み、被測定物Sを透過して散乱部41の画素41aに入射する。第1のX線光子は、画素41aの散乱点T31においてコンプトン散乱して進路がL12に変化し、吸収部42の画素42aに入射する。一方、X線源2からは、上記の第1のX線光子が放射されてから短時間の後に、別のX線光子(第2のX線光子)が放射される。第2のX線光子は、進路L11とは異なる進路L21を進み、被測定物Sを透過して散乱部41の散乱部41の画素41bに入射する。第2のX線光子は、画素41bの散乱点T32においてコンプトン散乱して進路がL22に変化し、吸収部42の画素42aに入射する。
 すなわち、散乱部41の異なる2つの画素41aと41bでそれぞれコンプトン散乱した異なるX線光子が、同一イベントとして吸収部42の同一画素42aに入射する。このような場合には、散乱部41および吸収部42から得られたエネルギー情報に基づいてコンプトン散乱角θ1を正しく算出できない。このような状態は、例えば、X線源2から放射されるX線光子の数が多く、X線源2からX線光子が放射されたL11を進み始めてすぐに、別のX線光子がX線源2から放射されたような場合に発生することが考えられる。
 このような状態を判別するための手順について次に説明する。特定部54は、単一イベントにおいて、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2との和が、所定値ETLを超えるか否かを判定する。ETLの値は、X線源2から正常状態で放射される単一のX線光子のエネルギーに基づいて設定する。例えば、駆動加速電圧によって決まるX線源2のエネルギーに基づくX線エネルギーの最大値にエネルギー分解能を加えた値に設定する。エネルギーE1とエネルギーE2との和が所定値ETLを超える場合には、特定部54は、図3に示すように、複数のX線光子によるエネルギーE1とエネルギーE2とが算出されたものと判定する。この場合には、特定部54は、散乱部41により得られた2箇所の画素に対応するエネルギーE1と吸収部42により得られたエネルギーE2とを散乱X線の特定に用いない。
 一方、エネルギーE1とエネルギーE2との和が所定値ETL以下で、それぞれが散乱部41および吸収部42におけるそれぞれの単一の画素に対応し、散乱部41および吸収部42のそれぞれにおいて、ほぼ同時刻に複数個の画素に対応するエネルギー値を持たない場合には、これらのエネルギー情報は、単一のX線光子に基づいて正常に出力されたと判定する。この場合には、特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とを用いて、上述した式(1)に基づいて、コンプトン散乱角θ1を算出し、図2を参照して説明したように、X線源2から放射されたX線が被測定物Sの内部において散乱したか否かの判定を行う。一方、エネルギーE1とエネルギーE2との和が所定値ETL以下であっても、散乱部41および吸収部42の少なくとも一方において複数の異なる位置の画素にてほぼ同時刻に電気信号が検出された場合には、特定部54は、エネルギーE1とエネルギーE2とを散乱X線の特定に用いない。
 すなわち、特定部54は、算出したコンプトン散乱角θ1と、取得した位置(XSi、YSi、ZSi)および位置(XAj、YAj、ZAj)とを用いて、図2に示す円錐面CS上またはその近傍にX線源2が存在するか否かを判定する。具体的には、特定部54は、位置(XSi、YSi、ZSi)と、位置(XAj、YAj、ZAj)と、算出したコンプトン散乱角θ1とから、位置(XSi、YSi、ZSi)および位置(XAj、YAj、ZAj)を通る直線を中心軸として、位置(XSi、YSi、ZSi)を頂点とし、母線と中心軸のなす角がθ1の円錐(直円錐)を設定する。特定部54は、この円錐の円錐面CS上または円錐面CSからの距離が所定の距離よりも小さい位置にX線源2が存在するか否か判定する。所定の長さは、上記の通り、画素41aおよび画素42aの大きさや、X線装置100における、X線源2、被測定物S、検出器4等の配置等に基づいて、適宜設定することが好ましい。
 特定部54により、X線源2が円錐面CS上または円錐面CSからの距離が所定の距離よりも小さい位置に存在すると判定された場合、X線源2から放射されたX線は、被測定物Sの内部で散乱することなしに、散乱部41に入射し、散乱部41でコンプトン散乱した後、吸収部42に到達したと判断し、散乱X線以外のX線であることを特定する。この場合には、画像生成部53は、散乱部41と吸収部42とから得られたエネルギー情報(電気信号)に基づいてX線の強度分布データを生成し、X線投影画像データを生成する際の被測定物Sの内部情報として用いる。すなわち、画像生成部53は、被測定物Sの内部情報として用いる。一方、特定部54により、X線源2が円錐面CS上または円錐面CSからの距離が所定の距離よりも小さい位置に存在しないと判定された場合、X線源2から放射されたX線は、被測定物Sの内部で散乱して散乱部41に入射し、散乱部41でコンプトン散乱した後、吸収部42に到達したと判断し、散乱X線であることを特定する。この場合には、画像生成部53は、散乱部41と吸収部42とから得られたエネルギー情報(電気信号)に基づくX線の強度分布データの生成は行わない。すなわち、画像生成部53は、X線投影画像データを生成する際の被測定物Sの内部情報として用いることはない。
 上述した処理を行うことにより、特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とに基づいて、被測定物Sを透過した透過X線(または被測定物Sに照射されずに直進した透過X線)と、被測定物Sで散乱した散乱X線との少なくとも一方を特定する。なお、上述したように、散乱部41により検出器4内でX線が散乱した散乱位置が検出され、吸収部42によりX線が入射した入射位置が検出されるので、特定部54は、散乱位置と入射位置とに基づいて、透過X線と散乱X線との少なくとも一方を特定することができる。
 上記の散乱検出処理が行われることにより、画像生成部53により生成されるX線投影画像データは、散乱を起こした散乱X線に起因する被測定物Sのエッジ等のボケ等が抑制され、その結果、被測定物Sの高画質画像を生成することができる。
 なお、X線源2から放射されて検出器4の散乱部41に入射したX線の一部は、散乱部41でコンプトン散乱せずに、吸収部42に入射することがある。この場合には、単一イベントにおいて、吸収部42からはX線のエネルギー情報は出力されるものの、散乱部41からはX線のエネルギー情報に対応する電気信号は出力されない。このような場合は、画像生成部53は、吸収部42から出力されたエネルギー情報(電気信号)に基づくX線の強度分布データの生成は行わない。すなわち、画像生成部53は、X線投影画像データを生成する際の被測定物Sの内部情報として用いることはない。
 画像再構成部56は、X線源2から被測定物Sに照射するX線の照射角度を変化させながら所定の測定角度ごとに生成された複数のX線投影画像データに対して、公知の画像再構成処理を施して、再構成画像を生成する。これにより、被測定物Sの内部構造(断面構造)である断面画像データや3次元データが生成される。
 図4に示すフローチャートを参照して、制御装置5が行う動作について説明する。図4に示す処理は制御装置5でプログラムを実行して行われる。このプログラムは、メモリ(不図示)に格納されており、制御装置5により起動され、実行される。
 ステップS1では、制御装置5のX線制御部51は、X線源2にX線の照射を開始させてステップS2へ進む。ステップS2では、制御装置5の特定部54は、検出器4の散乱部41の位置(XSi、YSi、ZSi)の画素および吸収部42の位置(XAj、YAj、ZAj)の画素から出力されたX線のエネルギーE1およびE2を取得してステップS3へ進む。
 ステップS3では、特定部54は、散乱部41へX線が入射した時刻Tiと吸収部42へX線が入射した時刻Tjとの時間差ΔTijが所定の値T1以下であるか否かを判定する。時刻TjとTiとの時間差ΔTijが所定の値T1以下の場合、すなわち、散乱部41に入射したX線と、吸収部42に入射したX線とが、同一イベント中のものであると判定された場合には、ステップS3が肯定判定されてステップS4へ進む。時刻TjとTiとの時間差ΔTijが所定の値T1を超える場合、すなわち、散乱部41に入射したX線と、吸収部42に入射したX線とが、同一イベント中のものではないと判定された場合には、ステップS3が否定判定されて、後述するステップS9へ進む。
 ステップS4では、吸収部42に入射したX線が、被測定物Sの異なる位置を透過したものであるか否かを判定する。散乱部41からのエネルギーE1と吸収部42からのエネルギーE2との和ETが、所定値ETLを超える場合、X線は被測定物Sの異なる位置を透過したと判定され、ステップS4が否定判定されて後述するステップS9へ進む。散乱部41からのエネルギーE1と吸収部42からのエネルギーE2との和ETが所定値ETL以下の場合、X線は、X線源2から正常に放射され、単一のX線光子に基づくエネルギー情報(電気信号)が散乱部41および吸収部42から出力されたと判定され、ステップS4が肯定判定されステップS5へ進む。ステップS5では、検出器4の単一箇所(単一画素)からエネルギーE1、E2に対応する電気信号が出力されたか否かを判定する。検出器4の散乱部41および吸収部42のそれぞれにおいて単一個所からエネルギーE1、E2に対応する電気信号が出力された場合には、ステップS5が肯定判定されて、ステップS6へ進む。散乱部41および吸収部42の少なくとも一方において複数個所(複数画素)からエネルギーE1、E2に対応する電気信号が出力された場合には、ステップS5が否定判定されて後述するステップS9へ進む。
 ステップS6では、特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とを用いて、上述した式(1)に基づいて、コンプトン散乱角θ1を算出してステップS7へ進む。ステップS7では、特定部54は、被測定物Sの内部でX線が散乱を起こしたか否かを判定する。具体的には、算出したコンプトン散乱角θ1とから、位置(XSi、YSi、ZSi)および位置(XAj、YAj、ZAj)を通る直線を中心軸として、位置(XSi、YSi、ZSi)を頂点とし、母線と中心軸のなす角がθ1の円錐(直円錐)を設定する。特定部54は、この円錐の円錐面CS上または円錐面CSからの距離が所定の距離よりも小さい位置にX線源2が存在するか否か判定する。被測定物Sの内部でX線が散乱していないと判定された場合には、ステップS7が肯定判定されて、ステップS8へ進む。一方、被測定物Sの内部でX線が散乱したと判定された場合には、ステップS7が否定判定されて、後述するステップS9へ進む。
 ステップS8では、画像生成部53は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2とを、散乱部41に入射したX線の位置(XSi、YSi、ZSi)を入射位置とする内部情報として記憶してステップS9へ進む。すなわち、後述する処理においては、散乱部41に入射したX線の位置(XSi、YSi、ZSi)に基づいて、画像生成部53によりX線投影画像データが生成される。ステップS9では、制御装置5は、検出器4のすべてのイベントについて特定処理を行ったか否かを判定する。全てのイベントについて特定処理が行われた場合には、ステップS9が肯定判定されてステップS10へ進む。ステップS10では、画像生成部53は、生成された内部情報を用いて、被測定物SのX線投影画像データを生成して、後述するステップS12へ進む。
 ステップS9において特定処理が行われていないイベントが存在する場合には、ステップS9が否定判定されてステップS11へ進む。ステップS11では、特定部54は、次のイベントでの散乱部41の位置(XSi、YSi、ZSi)の画素および吸収部42の位置(XAj、YAj、ZAj)の画素から出力された電気信号(すなわち、エネルギーE1およびE2に対応する電気信号)を取得してステップS3へ戻り、以後同様の処理を行う。
 ステップS12では、制御装置5は、全ての測定角度についてX線投影画像データの生成が行われたか否かを判定する。全ての測定角度についてX線投影画像データの生成が行われた場合には、ステップS12が肯定判定されてステップS13へ進む。ステップS13では、画像再構成部56は、複数のX線投影画像データを用いて3次元データを生成して処理を終了する。なお、制御装置5は、この3次元データをモニタ(不図示)に表示したり、メモリ(不図示)に記憶させたりすることができる。
 ステップS12においてX線投影画像データが生成されていない測定角度が存在する場合には、ステップS12が否定判定されてステップS14へ進む。ステップS14においては、載置台制御部52はマニピュレータ部36を駆動させて、載置台31を所定の測定角度に回転駆動させてステップS2へ戻り、以後同様の処理を行う。
 上述した第1の実施の形態によれば、次の作用効果が得られる。
(1)制御装置5の特定部54は、散乱部41により得られたX線の損失エネルギーE1と、吸収部42により得られた散乱X線のエネルギーE2とに基づいて、被測定物Sから検出器4に入射したX線のうち、被測定物Sの内部で散乱を起こした散乱X線を特定する。また、制御装置5の特定部54は、散乱部41により検出された検出器4内でのX線の散乱位置と、吸収部42により検出されたX線の吸収部42への入射位置とに基づいて、透過X線と散乱X線との少なくとも一方を特定することができる。これにより、X線のうち、被測定物Sで散乱した散乱X線と、散乱していないX線とを選別することが可能になる。X線投影画像データを生成する際には、ノイズの要因となる散乱X線を被測定物Sの内部情報から除外することが可能になり、X線投影画像データの画質を向上させることができる。
(2)特定部54は、散乱部41で散乱したX線の散乱角θ1と、散乱部41に入射したX線の位置(散乱位置)と、吸収部42に入射した散乱X線の位置(入射位置)とが満たす関係に基づいて、散乱X線を特定する。これにより、検出器4に入射したX線のうち散乱X線を除外して、X線投影画像データを生成することができる。
(3)特定部54は、散乱角θ1と、散乱部41に入射したX線の位置(すなわち散乱位置)と、吸収部42に入射したX線の位置(すなわち入射位置)とに基づいて、X線源2の存在が推定される位置(円周C2)を算出し、算出された位置にX線源2が存在すると、検出器4に入射したX線を被測定物Sでの散乱X線以外の透過X線として特定する。これにより、検出器4に入射したX線のうち散乱X線を除外できるので、X線投影画像データの画質向上に寄与することができる。
(4)特定部54は、散乱部41からエネルギーE1に対応する電気信号が出力された時刻Tiと、吸収部42からエネルギーE2に対応する電気信号が出力された時刻Tjとの時間差が所定の範囲T1以内のX線を用いて散乱X線を特定する。これにより、異なるタイミングでX線源2から出射され、検出器4に入射したX線を除外して散乱X線を特定できるので、散乱X線を特定する精度を向上させることができる。
(5)特定部54は、X線源2から照射され、被測定物Sの異なる位置を透過して、吸収部42の同一の画素に入射したX線を除外して、散乱X線を特定する。これにより、被測定物Sの内部の異なる位置の情報が、同一の位置の情報として内部情報に含まれることが抑制されるので、X線投影画像データの画質の低下を抑制できる。
(6)特定部54は、散乱部41により得られたエネルギーE1と吸収部42により得られたエネルギーE2との和が、X線源2から照射されたX線のエネルギーE0に基づく値を超えると、被測定物Sの異なる位置を透過して、吸収部42の同一の位置に入射した複数のX線光子であることを検出し、検出したX線光子を除外して散乱X線を特定する。これにより、被測定物Sの異なる位置を透過して、吸収部42の同一の画素に入射したX線を、被測定物Sの内部情報から除外することができる。
 上述した第1の実施の形態を、以下のように変形できる。
(1)上述した第1の実施の形態においては、検出器4の散乱部41と吸収部42とをシンチレーション物質を含むシンチレータ部により構成される場合を例に挙げて説明を行ったが、この例に限定されない。例えば、検出器4の散乱部41と吸収部42とは、入射するX線を光エネルギーに変換することなく、電気信号に直接変換して出力する半導体検出器であってもよい。散乱部41と吸収部42とは、CdTeを有し、入射したX線のエネルギーに比例した数の電子を生成し、電気信号として出力する。なお、この場合も、散乱部41のZ方向の厚さは、吸収部42のZ方向の厚さよりも小さい。
(2)上述した第1の実施の形態においては、画像生成部53は、散乱部41に入射したX線の位置(XSi、YSi、ZSi)に基づいて、X線投影画像データを生成する場合を例に挙げて説明を行った。しかし、画像生成部53は、吸収部42に入射したX線の位置(XAj、YAj、ZAj)に基づいて、X線投影画像データを生成してもよい。この場合、画像生成部53は、X線源2の位置と、散乱部41に入射したX線の位置(XSi、YSi、ZSi)とに基づいて、X線が散乱部41で散乱しないと仮定した場合に入射する吸収部42上の位置(XAj’、YAj’、ZAj’)を算出する。画像生成部53は、この算出した位置(XAj’、YAj’、ZAj’)に基づいて、エネルギーE1とエネルギーE2とからX線投影画像データを生成すればよい。
 図面を参照して、本発明の実施の形態による構造物製造システムを説明する。本実施の形態の構造物製造システムは、例えば自動車のドア部分、エンジン部分、ギア部分および回路基板を備える電子部品等の成型品を作成する。
 図5は本実施の形態による構造物製造システム400の構成の一例を示すブロック図である。構造物製造システム400は、第1の実施の形態にて説明したX線装置100と、設計装置410と、成形装置420と、制御システム430と、リペア装置440とを備える。
 設計装置410は、構造物の形状に関する設計情報を作成する際にユーザが用いる装置であって、設計情報を作成して記憶する設計処理を行う。設計情報は、構造物の各位置の座標を示す情報である。設計情報は成形装置420および後述する制御システム430に出力される。成形装置420は設計装置410により作成された設計情報を用いて構造物を作成、成形する成形処理を行う。この場合、成形装置420は、3Dプリンター技術で代表される積層加工、鋳造加工、鍛造加工および切削加工のうち少なくとも1つを行うものについても本発明の一態様に含まれる。
 X線装置100は、成形装置420により成形された構造物の形状を測定する測定処理を行う。X線装置100は、構造物を測定した測定結果である構造物の座標を示す情報(以後、形状情報と呼ぶ)を制御システム430に出力する。制御システム430は、座標記憶部431と、検査部432とを備える。座標記憶部431は、上述した設計装置410により作成された設計情報を記憶する。
 検査部432は、成形装置420により成形された構造物が設計装置410により作成された設計情報に従って成形されたか否かを判定する。換言すると、検査部432は、成形された構造物が良品か否かを判定する。この場合、検査部432は、座標記憶部431に記憶された設計情報を読み出して、設計情報とX線装置100から入力した形状情報とを比較する検査処理を行う。検査部432は、検査処理として例えば設計情報が示す座標と対応する形状情報が示す座標とを比較し、検査処理の結果、設計情報の座標と形状情報の座標とが一致している場合には設計情報に従って成形された良品であると判定する。設計情報の座標と対応する形状情報の座標とが一致していない場合には、検査部432は、座標の差分が所定範囲内であるか否かを判定し、所定範囲内であれば修復可能な不良品と判定する。
 修復可能な不良品と判定した場合には、検査部432は、不良部位と修復量とを示すリペア情報をリペア装置440へ出力する。不良部位は設計情報の座標と一致していない形状情報の座標を有する部位であり、修復量は不良部位における設計情報の座標と形状情報の座標との差分である。リペア装置440は、入力したリペア情報に基づいて、構造物の不良部位を再加工するリペア処理を行う。リペア装置440は、リペア処理にて成形装置420が行う成形処理と同様の処理を再度行う。
 図6に示すフローチャートを参照しながら、構造物製造システム400が行う処理について説明する。
 ステップS31では、設計装置410はユーザによって構造物の設計を行う際に用いられ、設計処理により構造物の形状に関する設計情報を作成し記憶してステップS32へ進む。なお、設計装置410で作成された設計情報のみに限定されず、既に設計情報がある場合には、その設計情報を入力することで、設計情報を取得するものについても本発明の一態様に含まれる。ステップS32では、成形装置420は成形処理により、設計情報に基づいて構造物を作成、成形してステップS33へ進む。ステップS33においては、X線装置100は測定処理を行って、構造物の形状を計測し、形状情報を出力してステップS34へ進む。
 ステップS34では、検査部432は、設計装置410により作成された設計情報とX線装置100により測定され、出力された形状情報とを比較する検査処理を行って、ステップS35へ進む。ステップS35では、検査処理の結果に基づいて、検査部432は成形装置420により成形された構造物が良品か否かを判定する。構造物が良品である場合、すなわち設計情報の座標と形状情報の座標との差が所定の範囲内の場合には、ステップS35が肯定判定されて処理を終了する。構造物が良品ではない場合、すなわち設計情報の座標と形状情報の座標とが一致しない場合や設計情報には無い座標が検出された場合には、ステップS35が否定判定されてステップS36へ進む。
 ステップS36では、検査部432は構造物の不良部位が修復可能か否かを判定する。不良部位が修復可能ではない場合、すなわち不良部位における設計情報の座標と形状情報の座標との差分が所定範囲を超えている場合には、ステップS36が否定判定されて処理を終了する。不良部位が修復可能な場合、すなわち不良部位における設計情報の座標と形状情報の座標との差分が所定範囲内の場合には、ステップS36が肯定判定されてステップS37へ進む。この場合、検査部432はリペア装置440にリペア情報を出力する。ステップS37においては、リペア装置440は、入力したリペア情報に基づいて、構造物に対してリペア処理を行ってステップS33へ戻る。なお、上述したように、リペア装置440は、リペア処理にて成形装置420が行う成形処理と同様の処理を再度行う。
 上述した実施の形態による構造物製造システムによれば、以下の作用効果が得られる。
(1)構造物製造システム400のX線装置100は、設計装置410の設計処理に基づいて成形装置420により作成された構造物の形状情報を取得する測定処理を行い、制御システム430の検査部432は、測定処理にて取得された形状情報と設計処理にて作成された設計情報とを比較する検査処理を行う。したがって、構造物の欠陥の検査や構造物の内部の情報を非破壊検査によって取得し、構造物が設計情報の通りに作成された良品であるか否かを判定できるので、構造物の品質管理に寄与する。
(2)リペア装置440は、検査処理の比較結果に基づいて、構造物に対して成形処理を再度行うリペア処理を行うようにした。したがって、構造物の不良部分が修復可能な場合には、再度成形処理と同様の処理を構造物に対して施すことができるので、設計情報に近い高品質の構造物の製造に寄与する。
-第2の実施の形態-
 図面を参照しながら、第2の実施の形態によるX線装置について説明する。以下の説明では、第1の実施の形態と同じ構成要素には同じ符号を付し、相違点を主に説明する。特に説明しない点については、第1の実施の形態と同様である。
 図7は第2の実施の形態によるX線装置100Aの構成の一例を示す図である。なお、図7においては、説明の都合上、X軸、Y軸、Z軸からなる座標系を図1と同様に設定する。
 X線装置100Aは、第1の実施の形態と同様の筐体1、X線源2、載置部3および制御装置5と、検出器4Aとを備えている。
 検出器4Aは、載置台31よりもZ軸+側に設けられている。すなわち、載置台31は、Z軸方向において、X線源2と検出器4との間に設けられる。検出器4Aは、散乱部41Aと吸収部42Aとを備える。散乱部41Aは、吸収部42AよりもZ方向-側、すなわち被測定物Sに近い側に設けられる。散乱部41Aには、X線源2から放射され、載置台31上に載置された被測定物Sを透過した、または被測定物S内部で散乱したX線が入射する。なお、以後の説明においても、被測定物Sの内部で散乱せずに内部を透過して検出器4に入射するX線を透過X線、被測定物Sの内部で散乱して検出器4Aに入射するX線を散乱X線と呼ぶ。
 散乱部41Aは、例えば公知の特開2017-67601号と同様の構成を有する。すなわち、図8に示すように、散乱部41Aは、チャンバー411Aとピクセル電極412Aとドリフト電極413Aとを有する。チャンバー411Aの内部には、アルゴンやキセノンなどの希ガスと、エタン、メタンなどの常温でガスのアルカンもしくは二酸化炭素を含む消光作用を有するガスとの混合ガスが封入されている。なお、チャンバー411Aは、混合ガスが封入されるものに限定されず単体のガスが封入されているものでもよい。ピクセル電極412Aはチャンバー411AのZ方向+側の面に設けられ、ドリフト電極413Aはチャンバー411のZ方向-側の面に設けられる。ピクセル電極412Aは、XY平面上に二次元的に配置された複数の画素を有する。チャンバー411A内においは、ピクセル電極412Aとドリフト電極413Aとに印加された電圧により、ピクセル電極412Aとドリフト電極413Aとの間に電場が形成される。
 被測定物SからのX線が散乱部41に入射すると、X線(光子)はチャンバー411A内の気体粒子を構成する電子とある確率で衝突し、コンプトン散乱する。コンプトン散乱したX線光子は、進行方向が変化した状態で散乱部41Aを出射し、Z方向+側に設けられた吸収部42Aに入射する。一方、チャンバー411A内でコンプトン散乱しなかったX線は、チャンバー411A内で進行方向が変化することなく散乱部41Aを出射し、吸収部42Aに入射する。
 チャンバー411A内でX線がコンプトン散乱した場合、X線光子に衝突された電子は反跳電子(荷電粒子)として進行して別の電子に衝突する。衝突された電子は反跳電子として進行し、さらに別の電子に衝突して二次電子を生成する。反跳電子は、電子への衝突を繰り返しながら徐々にエネルギーを失っていき、生成された多数の二次電子は反跳電子の飛跡の近傍に漂う。その結果、チャンバー411A内には、エネルギーを失った電子の分布が生成される。本明細書においては、エネルギーを失った電子の分布を電子雲とよぶ。すなわち、反跳電子の飛跡に沿って電子雲が発生する。電子雲を構成する二次電子は、ピクセル電極412Aとドリフト電極413Aとにより形成される電場によってZ方向+側へ引き寄せられ、ピクセル電極412Aに入射する。ピクセル電極412Aは、X方向に延在するカソード電極と、カソード電極に複数設けられた開口に配置されるアノード電極とからなる複数の画素を有し、電子雲を構成する二次電子が入射した画素は電気信号を制御装置5へ出力する。これにより、散乱部41は、チャンバー411A内に入射したX線の散乱により生じる電子の進行に関する情報を検出することができる。すなわち、電子雲の形状が得られる。電子の進行に関する情報には、反跳電子の初期進行方向(X線光子に衝突された電子の初期進行方向)と多数の二次電子の分布により表される反跳電子が順次進行した進行軌跡(飛跡)との少なくとも一方が含まれる。
 散乱部41Aを出射したX線は吸収部42Aに入射する。吸収部42Aは、公知のシンチレーション物質を含むシンチレータ部と、シンチレータ部から放出された光を増幅して受光する受光部等とによって構成される。シンチレータ部の入射面に入射したX線は吸収部42Aのシンチレータ部に吸収されて蛍光を放出し、放出された蛍光の光エネルギーは上記の受光部で光増倍管により増幅されて電気エネルギーに変換され、電気信号として制御装置5へ出力される。
 検出器4Aの吸収部42Aは、シンチレータ部と受光部とが複数の画素が二次元的に配置された構造を有し、シンチレータ部と受光部との複数の画素のそれぞれは、互いに対応するように配列されている。これにより、X線源2から放射され、被測定物Sを透過、散乱または一部が吸収されたX線の強度データを、複数の画素において一括して取得することができ、それぞれの強度データを統合することにより強度分布データを得ることができる。なお、吸収部42Aの複数の画素のそれぞれは、散乱部41Aのピクセル電極412Aが有する複数の画素のそれぞれと、互いに対応するように配列されている。
 なお、吸収部42Aのシンチレータ部は、複数の画素が二次元的に配置された構造を有していないものでもよい。
 本実施の形態においても、制御装置5の特定部54は、検出器4AへのX線の入射方向を推定し、推定した入射方向とX線源2の位置とに基づいて透過X線(または被測定物Sに照射されずに直進した透過X線)および散乱X線の少なくとも一方を特定する。すなわち、特定部54は、散乱部41Aおよび吸収部42Aから出力された電気信号に基づいて、X線源2からのX線が被測定物Sの内部で散乱したか否かを判定する。特定部54は、判定結果に基づいて、被測定物SからのX線のうち、被測定物Sにおいて散乱した散乱X線を特定、または散乱X線以外のX線を特定する。すなわち、特定部54は、検出器4Aに入射したX線が、被測定物Sで散乱したX線であるか否かを特定する。
 以下、特定部54による、被測定物SにおいてX線が散乱したか否かを判定し、透過X線および散乱X線の少なくとも一方を特定する手順について、図8を参照して説明する。図8は、散乱部41Aのピクセル電極412Aの複数の画素412a~412cにて電気信号を検出し、吸収部42Aの一つの画素42aにて電気信号を検出する場合を示している。なお、ピクセル電極412Aは、画素412c、412b、412aの順序で電気信号を検出したものとする。
 上述したように、ピクセル電極412Aは、散乱部41A内でのX線のコンプトン散乱により生じた反跳電子の飛跡に沿って発生した電子雲を構成する二次電子が、チャンバー411A内のドリフト電場によりZ方向+側へ引き寄せられ、これらの電子がピクセル電極412Aに到達した時点でその信号を検出する。例えば、画素412c、412b、412aにて電子が順次検出される際の時刻と、画素412c、412b、412aの二次元的な位置により、チャンバー411A内における反跳電子の進行に関する情報である電子雲の三次元形状(電子の飛跡)を求めることができる。
 電子雲の三次元形状を求める手順について、具体的に説明する。電子雲の発生に寄与した任意の一つの二次電子がドリフト電場に沿ってZ方向に移動する移動速度Vdはドリフト電場の強さにより決まる。移動速度Vdはチャンバー411A内において実質的に一定であると考えられる。ドリフト電極413Aからピクセル電極412Aまでドリフトする時間をTdとする。電子がピクセル電極412Aに到達した時点で、信号が検出される。
 X線がコンプトン散乱した後、吸収部42Aに達するまでに要する時間、および、それにより電子雲が発生するのに要する時間は、Tdに比べて極めて短いので、コンプトン散乱した時点と、散乱X線が吸収部42Aに到達した時点と、電子雲が形成された時点とは、実質的に同一と考えて問題ない。この点に関して、散乱部41Aのチャンバー411AのZ方向の厚みDを300mmとすると、例えば、X線がドリフト電極に垂直に入射してチャンバー411Aを通過して吸収部42Aに到達するのに要する時間は(0.3[m]/3)×10-8、すなわち1[nsec]である。また、反跳電子の初期速度を高速の1/10とした場合であっても、X線がコンプトン散乱してから電子雲が形成されるまでの時間は、長くても10[nsec]程度より短いと考えられる。従って、コンプトン散乱したX線が吸収部42Aに吸収された時刻をコンプトン散乱した時刻と見做してt=0とする。
 時刻t=0を基準として、ピクセル電極412Aのカソード電極Xi(i=1、m)とアノード電極Yj(j=1、n)から電気信号が読み出された時刻をt=tsigとする。すなわち、時刻t=tsigにて、カソード電極およびアノード電極にて電圧が変化したとする。この場合、ピクセル電極412Aの画素(Xi,Yj)に電圧変化をもたらした電子がドリフト開始時、すなわち、電子雲発生時にチャンバー411Aにおいて存在した位置は、Z方向の基準をピクセル電極412Aとした場合に、(Xi、Yj、Vd×tsig)と考えられる。なお、カソード電極とアノード電極とは時刻tsigにて同時に電位変化したとする。これにより、ピクセル電極412Aの各画素(Xi、Yj)ごとに、(Xi、Yj、Vd×tsig)が算出できる。このようにして、電子雲の発生に寄与した複数の二次電子がドリフトしてピクセル電極412Aに到達することで電圧が検出された時刻tsigと、ピクセル電極412Aにおける画素の位置に基づいて、電子雲の三次元形状(飛跡P1)を求めることができる。
 電子雲の発生に寄与した二次電子のうち、コンプトン散乱した位置に最も近い電子がドリフトしてピクセル電極に到達する時刻が最長(t=tsig0)となる。すなわち、時刻t=tsig0において検出された電子がドリフト開始時に存在した位置(Xi、Yj、Vd×tsig0)がコンプトン散乱位置(散乱点T1)である。また、時刻t=tsig0の次の大きさの時刻であるt=tsig1において検出された電子がドリフト開始時に存在した位置は(Xi、Yj、Vd×tsig1)となる。これらの二つの位置から、散乱点T1からの反跳電子の初期進行方向F1が得られる。
 なお、Z方向の基準をドリフト電極413Aの位置としてもよく、この場合には、電圧変化した画素に対応するチャンバー411Aでの位置は、(Xi、Yj、D(Td-tsig)/Td)となる。この場合、時刻t=tsigが最小のときのZ方向の位置(Vd×tsig)が散乱点T1となる。また、時刻t=tsigの値が所定値より小さな時の複数の位置(Xi、Yj、D(Td-tsig)/Td)に基づいて、反跳電子の初期進行方向F1が得られる。
 換言すると、散乱部41Aは、チャンバー411A内の散乱点T1においてX線がコンプトン散乱することにより生じる反跳電子の初期進行方向F1を検出する第1検出器として機能する。これにより、散乱部41Aは、チャンバー411AでX線が散乱したことにより生じる電子の進行に関する情報として、反跳電子の初期進行方向F1を検出する。
 散乱点T1でコンプトン散乱したX線は、進行方向が変化して吸収部42Aの画素42aに入射したと判定される。すなわち、吸収部42Aの画素42aは、吸収部42Aに入射したX線の入射位置であり、吸収部42Aは散乱部41AからのX線の入射を受け、入射したX線の入射位置を検出する第2検出器として機能する。このようにして、コンプトン散乱したX線の進行方向と、コンプトン散乱による反跳電子により発生した電子雲の形状を把握できる。また、吸収部42Aの画素42aと、散乱点T1に基づいて、コンプトン散乱したX線の進行に関する情報(X線の進行方向)が求められる。
 なお、反跳電子は、多数回のラザフォード散乱により一般に不規則なジグザグ飛跡を描くが、図8においては説明を簡単にするために、コンプトン散乱後の反跳電子の進路である飛跡P1は滑らかな曲線状に示している。
 特定部54は、上記説明に基づいて、散乱部41Aのピクセル電極412Aの画素412a~412cの位置関係と検出時刻から求めた電子雲の形状に基づいて、飛跡P1の長さ(最大長:飛程)を算出する。チャンバー411A内を進む反跳電子はエネルギーを徐々に失って最終的に停止する。従って、飛程から反跳電子のエネルギーを求めることができる。このようにして、特定部54は、飛程から反跳電子のエネルギー、すなわちチャンバー411Aに入射したX線がコンプトン散乱に伴って失われたエネルギーE1(=m)を算出する。ここで、電子の静止質量mは、m=9.1×10-31(kg)である。なお、飛程とエネルギーとの対応関係に関するデータはエネルギーデータとして予めメモリ(不図示)に格納されているものとする。特定部54は、吸収部42Aの画素42aで検出した電気信号に基づいて、チャンバー411A内でコンプトン散乱した後に画素42aに到達したX線のエネルギーE2(=hv’)を算出する。ここで、hはプランク定数である。
 チャンバー411Aに入射したX線と、コンプトン散乱したX線と、反跳電子とは、運動量保存則とエネルギー保存則とを用いて、以下の式(4)~(6)により表される。
Figure JPOXMLDOC01-appb-M000005
 なお、cは真空中の光速、θはコンプトン散乱角、φは反跳電子の出射角(初期進行方向の方位角)である。また、γはローレンツ因子であり以下の式(7)で表される。
Figure JPOXMLDOC01-appb-M000006
 上記の式(4)~(7)により、散乱角θは、以下の式(8)により表される。
Figure JPOXMLDOC01-appb-M000007
 特定部54は、チャンバー411Aに入射したX線のエネルギーE0(=hv)を、算出したエネルギーE1と検出されたエネルギーE2とを加算することにより算出する。特定部54は、算出されたエネルギーE0と検出されたエネルギーE2とを上記式(8)に用いて、チャンバー411Aにおけるコンプトン散乱によるX線の進行方向の変化角度であるコンプトン散乱角θを算出する。特定部54は、ピクセル電極412Aの画素412a~412cのうち、少なくとも散乱点T1のZ方向の投影位置に近い位置に相当する画素412a、412bの位置関係と電子を検出した時間差により、散乱点T1から出射した反跳電子の進行方向(初期進行方向の方位角)を得ることができる。なお、好ましくは、反跳電子が1回目のラザフォード散乱をするまでの位置に相当する直線状に並んだ複数の画素群からの出力に基づいて方位角を得るのがよく、方位角の検出精度を向上させることができる。
 特定部54は、検出器4AへのX線の入射方向を推定し、推定した入射方向とX線源2の位置とに基づいて、透過X線(または被測定物Sに照射されずに直進した透過X線)および散乱X線の少なくとも一方を特定する。以下、具体的に説明する。この場合、散乱点T1および吸収点T2(入射位置)に基づいて、散乱部41Aに入射して散乱したX線の進行に関する情報を取得し、電子の進行に関する情報とX線の進行に関する情報とに基づいて、透過X線(または被測定物Sに照射されずに直進した透過X線)および散乱X線の少なくとも一方を特定する。以下、具体的に説明する。
 既に求められている電子雲の三次元形状により散乱点T1の位置は既知なので、散乱点T1の位置および画素42aの位置から、チャンバー411A内で散乱し吸収部42Aに入射したX線の進路L2が決まる。すなわち、進路L2は、図示の例では、散乱点T1の位置および画素42aの位置を直線で結んだ経路である。進路L2、および、散乱点T1から出射した反跳電子の初期進行方向F1により、散乱X線の進路L2と反跳電子の初期進行方向F1とを含む面(散乱面)が決まる。尚、散乱点T1の位置および画素42aの位置に基づいて、散乱点T1からのX線の進行方向である散乱方向を求め、散乱方向と反跳電子の初期進行方向F1とに基づいて散乱面を決定してもよい。この場合、散乱方向は、進路L2に基づいて算出してもよい。また、反跳電子の初期進行方向F1を用いて散乱面を決定した例を示したが、反跳電子の飛跡P1と、散乱X線の進路L2または散乱方向とを用いて散乱面を決定してもよい。チャンバー411Aに入射したX線の進路L1は、上記の散乱面内に存在するので、進路L1と進路L2とは、散乱面内において、散乱点T1でコンプトン散乱角θをなして交差する。このようにして、進路L1を求めることができ、X線源2が進路L1の延長上、または、その近傍にある場合には、X線源2から出射したX線は、被測定物Sの内部で散乱されることなしにチャンバー411Aに入射したと判定される。この場合を図9(a)に示す。図9(a)では、図8と同様に、散乱部41Aの画素412a~412cと吸収部42Aの画素42aで電気信号が出力された場合を示す。このような場合に、X線源2が、進路L1の延長上、または、その近傍に存在する。一方、X線源2が進路L1の延長上、または、その近傍にない場合には、X線源2から出射したX線は、被測定物Sの内部で散乱されて進行方向が変化した状態でチャンバー411Aに入射したと判定される。この場合を図9(b)に示す。図9(b)は、散乱部41Aの画素412aおよび412fと吸収部42Aの画素42aで電気信号が出力された場合を示す。図9(b)においては、反跳電子の初期進行方向はF2であり、X線は被測定物Sの内部の散乱点T3にて散乱する。このような場合に、X線源2が、進路L1の延長上、または、その近傍に存在しない。
 なお、特定部54は、反跳電子の初期進行の方位角φは、上記の式(4)~(7)から得られた以下の式(9)を用いてを算出することができる。
Figure JPOXMLDOC01-appb-M000008
 式(9)を用いて方位角(出射角)φを算出した場合、上記の進路L1は、散乱点T1を頂点とし反跳電子の初期進行方向F1を中心軸として、この中心軸と母線とのなす角がφの円錐と、散乱点T1を頂点とし進路L2を中心軸とし、この中心軸と母線とのなす角がθの円錐との接線に相当する。
 なお、上記の散乱角θと方位角φの少なくとも一方が算出されれば、チャンバー411Aに入射したX線の進路L1を決定することができるが、散乱角θと方位角φの両方を算出することにより、進路L1をより高精度に決定することが可能となる。
 上述した通り、進路L1上、または進路L1の近傍にX線源2が存在すると判定された場合には、特定部54は、X線が被測定物Sで散乱していないと判定し、散乱X線以外のX線であることを特定する。一方、進路L1上および進路L1の近傍にX線源2が位置していない場合には、特定部54は、X線が被測定物Sで散乱したと判定し、散乱X線であることを特定する。
 次に、特定部54が行う、被測定物SでX線が散乱したか否かを特定し、散乱X線と散乱X線以外のX線とを特定するための処理(以後、特定処理と呼ぶ)について説明する。
 本実施の形態においては、散乱部41に入射するX線のエネルギーET0を一定と仮定した場合、式(8)で算出されるコンプトン散乱角θは次の式(10)を満たす。
Figure JPOXMLDOC01-appb-M000009
 ここで、電子の静止質量m=9.1×10-31(kg)、真空中の光速c=3.0×10(m)として、m=511(keV)である。
 エネルギーET0をエネルギーE1とエネルギーE2との和と仮定した場合、検出精度Δθは、式(10)を微分して求めた式(11)の関係を満たす。
Figure JPOXMLDOC01-appb-M000010
 この式(10)、(11)から、検出精度Δθは、次の式(12)により算出される。
Figure JPOXMLDOC01-appb-M000011
 なお、
Figure JPOXMLDOC01-appb-M000012
であるものとする。
 例えば、X線源2から放射されて散乱部41Aに入射するX線のエネルギーを511keV、エネルギー分解能を3パーセント、コンプトン散乱角θが20°と仮定した場合、検出精度Δθは5.3°となる。また、エネルギー分解能を1パーセントとした場合には、検出精度Δθは1.8°となる。なお、散乱部41Aの射出面(ピクセル電極412側)と吸収部42Aとの間のZ方向の距離を20mm、散乱部41Aのピクセル電極412Aや吸収部42Aの各画素のサイズを200μm×200μmと仮定した場合、散乱角θの分解能は10mradとなる。
 特定部54は、検出器4Aの散乱部41Aの出力に基づいて算出したエネルギーE1、電気信号を出力したピクセル電極412Aの画素の位置、および画素から電気信号が出力された時刻tsigの各情報を取得する。特定部54は、ピクセル電極412Aにおいて電気信号を出力したそれぞれの画素の位置と時刻tsigとに基づいて、電子雲の三次元形状である飛跡を算出する。時刻tsigのうち最も遅い時刻tsig0で電気信号を出力した画素の位置(Xi、Yj)と、位置(Xi、Yj)から電子雲が発生した点までのZ方向の位置Vd×tsig0とにより決まる位置(Xi、Yj、Vd×tsig0)が散乱点T1に相当する。図8、図9(a)、および図9(b)においては、画素412aの位置から、算出した距離Vd×tsigだけZ軸-方向に離れた位置が散乱点T1に相当する。特定部54は、時刻t=tsigの値がtsig0に近い時刻に対応する複数の電子の位置(Xi、Yj、Vd×tsig)に基づいて、反跳電子の初期進行方向F1を算出する。
 特定部54は、検出器4Aの吸収部42Aから得られたエネルギーE2、エネルギーE2に対応する電気信号を出力した画素の位置(XAj、YAj、ZAj)、およびエネルギーE2が吸収部42Aから出力された時刻Tjの各情報を取得する。この画素の位置(XAj、YAj、ZAj)は、図8、図9(a)、および図9(b)のそれぞれにおける吸収点T2に相当する。なお、時刻Tjは、X線が吸収部42Aに入射した時刻に相当する。
 特定部54は、画素412aから電気信号が出力された時刻Tj+tsigと、エネルギーE2に対応する電気信号が出力された時刻Tjとに基づいて、散乱部41Aに入射したX線と、吸収部42Aに入射したX線とが、X線源2からあるタイミングで放射された同一X線光子であるのか否かを判定する。具体的には、特定部54は、時刻Tj+tsigとTjとの時間差(すなわちtsig)が所定時間T1(=Td)以内のときに、散乱部41Aからの電気信号に基づいて算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とは、同一のX線光子について出力されたものであると判定する。例えば、所定時間T1は、計算あるいは各種の試験やシミュレーション等により決定され、予めメモリ(不図示)に記憶されている。時刻Tj+tsigとTjとの時間差であるtsigが所定時間T1を超える場合には、特定部54は、算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とは、同一のX線光子について出力されたものではないと判定する。時刻tsigとTjとの時間差が所定時間T1を超える場合には、特定部54は、算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とを、散乱X線を特定するための処理に用いないようにする。ここで、前記のように、散乱部41Aと吸収部42Aとの間のZ方向の距離を20mmとおけば、時刻Tj+tsigとTjとの時間差tsigは0.1nsec程度の量である。なお、散乱部41Aと吸収部42Aとの間の距離は、検出器4Aの各画素サイズと、上述した式(12)で表される検出精度Δθとに基づいて、好ましい値に決定することができる。
 時刻Tj+tsigとTjとの時間差tsigが所定時間T1以内と判定されると、特定部54は、算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とに基づいて、コンプトン散乱角θを算出すべきかどうか判定する。コンプトン散乱角θを算出すべきでない場合について、次に説明する。
 検出器4Aの読み出し回路は、所定の時間ごとに、散乱部41Aから電気信号を出力し、吸収部42AからエネルギーE2を出力する。本明細書においては、上記所定の時間をエネルギー出力間隔Tと呼ぶ。本実施の形態においては、シンチレータの蛍光発光時間を数十nsecとした場合、T=100nsecである。上記の所定時間T1は、エネルギー出力間隔T以下に設定する必要がある。すなわち、
 T1≦T
を満たすように設定する。これにより、エネルギー出力間隔T内に、散乱部41Aから電気信号が出力され、吸収部42AからエネルギーE2に対応する電気信号が出力されるようにできる。
 X線源2からは短い時間スケールで見た場合、X線光子が1個ずつ放射される。X線源2からX線光子を放射する平均放射間隔Teと検出器4からエネルギー情報(すなわち、ピクセル電極412および吸収部42Aの画素42aからの電気信号)を出力する出力間隔Tとの関係について説明する。X線源2から1回の放射により放射されたX線光子が、チャンバー411Aおよび吸収部42Aにそれぞれ入射し、検出器4Aの散乱部41Aおよび吸収部42Aの全ての画素からエネルギー情報の算出に用いる電気信号の出力を完了するまでを、1イベントと呼ぶ。すなわち、X線源2から放射された単一X線光子に基づく1フレーム分のエネルギー情報を検出器4から得ることを1イベントと呼ぶ。チャンバー411Aにおけるコンプトン散乱角θを正しく求めるためには、1イベントの間に、X線源2から放射された単一のX線光子の検出値のみを後述する処理に用いるようにする。2回以上の放射によるX線光子が、同時に検出器4Aに入射した場合には、コンプトン散乱角θが正しく求められない可能性があるからである。これをイベントの分離と呼ぶ。従って、放射間隔Teは、1フレーム分のエネルギー情報を得るために要する時間であるエネルギー出力間隔Tより大きいことが必要である。すなわち、
 T≦Te
を満たす必要がある。
 イベントの分離ができずに、コンプトン散乱角θが正しく求められない場合について、次に説明する。図10は、X線源2からの放射間隔Teがエネルギー出力間隔Tより小さい場合について示す。なお、図10においては、説明を簡単にするために、X線源2から2つのX線光子が異なる方向にほぼ同時に放射された場合を示す。1つのX線光子(第1のX線光子)は、X線源2から進路L11を進み、被測定物Sを透過して散乱部41Aのチャンバー411Aに入射し、散乱点T31においてコンプトン散乱して進路がL12に変化し、吸収部42Aの画素42aに入射する。散乱点T31では初期進行方向をF3とする反跳電子が発生し、その飛跡に応じて生成した電子雲を構成した二次電子が到達したピクセル電極412Aの画素から電気信号が出力される。図10においては、ピクセル電極412Aの画素412aのZ軸-方向に散乱点T31が存在する。一方、X線源2からは、上記の第1のX線光子の放射直後(ほぼ同時)に、別のX線光子(第2のX線光子)が放射される。第2のX線光子は、進路L11とは異なる進路L21を進み、被測定物Sを透過してチャンバー411Aに入射し、散乱点T32においてコンプトン散乱して進路がL22に変化し、吸収部42Aの画素42aに入射する。散乱点T32では初期進行方向をF4とする反跳電子が発生し、その飛跡に応じて生成した電子雲を構成した二次電子が到達したピクセル電極412Aの画素から電気信号が出力される。図10においては、ピクセル電極412Aの画素412xのZ軸-方向に散乱点T32に存在する。
 すなわち、散乱部41Aの異なる2つの画素412aと412xで検出された散乱点T31とT32のそれぞれにてコンプトン散乱した異なるX線光子が、同一イベントとして吸収部42Aの同一画素42aに入射する。このような場合には、散乱部41Aおよび吸収部42Aから出力されたエネルギー情報に基づいてコンプトン散乱角θを正しく算出できない。このような状態は、例えば、X線源2から放射されるX線光子の数が多く、X線源2からX線光子が放射されたL11を進み始めてすぐに、別のX線光子がX線源2から放射されたような場合に発生することが考えられる。
 このような状態を判別するための手順について次に説明する。特定部54は、単一イベントにおいて、散乱部41Aからの電気信号の出力に基づいて算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2との和が、所定値ETLを超えるか否かを判定する。ETLの値は、X線源2から正常状態で放射される単一のX線光子のエネルギーに基づいて設定する。例えば、駆動加速電圧によって決まるX線源2のエネルギーに基づくX線エネルギーの最大値にエネルギー分解能を加えた値に設定する。エネルギーE1とエネルギーE2との和が所定値ETLを超える場合には、特定部54は、図4に示すように、複数のX線光子によるエネルギーE1とエネルギーE2とが得られたたものと判定する。この場合には、特定部54は、散乱部41Aからの電気信号の出力に基づいて算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とを散乱X線の特定に用いない。
 一方、エネルギーE1とエネルギーE2との和が所定値ETL以下で、散乱部41Aでは単一の飛跡が検出され、吸収部42Aにおいて単一の画素によりエネルギーE2が検出され、吸収部42Aにおいて、ほぼ同時刻に複数個所の画素に対応するエネルギー値を持たない場合には、これらのエネルギー情報は、単一のX線光子に基づいて正常に出力されたと判定する。この場合には、特定部54は、散乱部41Aからの電気信号の出力に基づいて算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2とを用いて、上述した式(8)に基づいて、コンプトン散乱角θを算出し、図8、図9を参照して説明したように、X線源2から放射されたX線が被測定物Sの内部において散乱したか否かの判定を行う。一方、エネルギーE1とエネルギーE2との和が所定値ETL以下であっても、散乱部41Aにおいて複数の異なる飛跡がほぼ同時刻に検出された場合、または吸収部42Aにおいて複数の異なる位置の画素にてほぼ同時刻に電気信号が検出された場合には、特定部54は、エネルギーE1とエネルギーE2とを散乱X線の特定に用いない。
 すなわち、特定部54は、算出したコンプトン散乱角θと、反跳電子の初期進行方向F1と、散乱点T1と、吸収部42Aで検出した画素42aの位置(XAj、YAj、ZAj)とを用いて、図9に示す進路L1上またはその近傍にX線源2が存在するか否かを判定する。具体的には、特定部54は、反跳電子の初期進行方向F1と、散乱点T1の位置と、位置(XAj、YAj、ZAj)と、算出したコンプトン散乱角θとから、反跳電子の初期進行方向F1と、散乱点T1の位置と、位置(XAj、YAj、ZAj)とを含む平面上において、散乱点T1の位置および位置(XAj、YAj、ZAj)を通る直線を進路L1として設定する。特定部54は、この進路L1上または進路L1からの距離が所定の距離よりも小さい位置にX線源2が存在するか否か判定する。所定の距離は、上記の通り、ピクセル電極412Aの各画素および画素42aの大きさや、X線装置100における、X線源2、被測定物S、検出器4A等の配置等に基づいて、適宜設定することが好ましい。
 特定部54により、X線源2が進路L1上または進路L1からの距離が所定の距離よりも小さい位置に存在すると判定された場合、X線源2から放射されたX線は、被測定物Sの内部で散乱することなしに、散乱部41Aに入射し、散乱部41Aでコンプトン散乱した後、吸収部42Aに到達したと判断し、散乱X線以外のX線であることを特定する。この場合には、画像生成部53は、吸収部42Aから出力されたエネルギー情報(電気信号)に基づいてX線の強度分布データを生成し、X線投影画像データを生成する際の被測定物Sの内部情報として用いる。すなわち、画像生成部53は、被測定物Sの内部情報として用いる。一方、特定部54により、X線源2が進路L1上または進路L1からの距離が所定の距離よりも小さい位置に存在しないと判定された場合、X線源2から放射されたX線は、被測定物Sの内部で散乱して散乱部41Aに入射し、散乱部41Aでコンプトン散乱した後、吸収部42Aに到達したと判断し、散乱X線であることを特定する。この場合には、画像生成部53は、吸収部42Aから出力されたエネルギー情報(電気信号)に基づくX線の強度分布データの生成は行わない。すなわち、画像生成部53は、X線投影画像データを生成する際の被測定物Sの内部情報として用いることはない。
 上記の散乱検出処理が行われることにより、画像生成部53により生成されるX線投影画像データは、散乱を起こした散乱X線に起因する被測定物Sのエッジ等のボケ等が抑制され、その結果、被測定物Sの高画質画像を生成することができる。
 なお、X線源2から放射されて検出器4Aの散乱部41Aに入射したX線の一部は、散乱部41Aでコンプトン散乱せずに、吸収部42Aに入射することがある。この場合には、単一イベントにおいて、吸収部42AからはX線のエネルギー情報は出力されるものの、散乱部41Aからは電気信号は出力されない。このような場合は、画像生成部53は、吸収部42Aから出力されたエネルギー情報(電気信号)に基づくX線の強度分布データの生成は行わない。すなわち、画像生成部53は、X線投影画像データを生成する際の被測定物Sの内部情報として用いることはない。
 画像再構成部56は、X線源2から被測定物Sに照射するX線の照射角度を変化させながら所定の測定角度ごとに生成された複数のX線投影画像データに対して、公知の画像再構成処理を施して、再構成画像を生成する。これにより、被測定物Sの内部構造(断面構造)である断面画像データや3次元データが生成される。
 図11に示すフローチャートを参照して、制御装置5が行う動作について説明する。図11に示す処理は制御装置5でプログラムを実行して行われる。このプログラムは、メモリ(不図示)に格納されており、制御装置5により起動され、実行される。
 ステップS101では、制御装置5のX線制御部51は、X線源2にX線の照射を開始させてステップS102へ進む。ステップS102では、制御装置5の特定部54は、検出器4Aの散乱部41Aのピクセル電極412Aから出力された電気信号を取得して、反跳電子の飛跡を検出し、飛跡の最大長と反跳電子の初期進行方向F1と散乱点T1とエネルギーE1とを算出してステップS103へ進む。ステップS103では、特定部54は、吸収部42Aの位置(XAj、YAj、ZAj)の画素から出力された電気信号に基づいて、X線のエネルギーE2を取得(算出)してステップS104へ進む。
 ステップS104では、特定部54は、散乱部41AへX線が入射した時刻Tj+tsigと吸収部42AへX線が入射した時刻Tjとの時間差ΔTijが所定の値T1以下であるか否かを判定する。時刻TjとTj+tsigとの時間差ΔTijが所定の値T1以下の場合、すなわち、散乱部41Aに入射したX線と、吸収部42Aに入射したX線とが、同一イベント中のものであると判定された場合には、ステップS104が肯定判定されてステップS105へ進む。時刻TjとTj+tsigとの時間差ΔTijが所定の値T1を超える場合、すなわち、散乱部41Aに入射したX線と、吸収部42Aに入射したX線とが、同一イベント中のものではないと判定された場合には、ステップS104が否定判定されて、後述するステップS110へ進む。
 ステップS105では、吸収部42Aに入射したX線が、被測定物Sの異なる位置を透過したものであるか否かを判定する。ステップS102で算出されたエネルギーE1とステップS103で取得された吸収部42Aにより得られたエネルギーE2との和ETが、所定値ETLを超える場合、X線は被測定物Sの異なる位置を透過したと判定され、ステップS105が否定判定されて後述するステップS110へ進む。エネルギーE1とエネルギーE2との和ETが所定値ETL以下の場合、X線は、X線源2から正常に放射され、単一のX線光子に基づくエネルギー情報が散乱部41Aおよび吸収部42Aから出力されたと判定され、ステップS105が肯定判定されステップS106へ進む。ステップS106では、検出器4Aで単一の飛跡が検出され、検出器4Aの単一箇所(単一画素)からエネルギーE2に対応する電気信号が出力されたか否かを判定する。検出器4Aの散乱部41Aの出力に基づいて、単一の飛跡が検出され、かつ、吸収部42Aの単一個所からエネルギーE2が出力された場合には、ステップS106が肯定判定されて、ステップS107へ進む。散乱部41Aの出力に基づいて複数の飛跡が検出された場合、または吸収部42Aの複数個所(複数画素)からエネルギーE2に対応する電気信号が出力された場合には、ステップS106が否定判定されて後述するステップS110へ進む。
 ステップS107では、特定部54は、算出されたエネルギーE1と吸収部42から出力されたエネルギーE2とを用いて、上述した式(8)に基づいて、コンプトン散乱角θを算出してステップS108へ進む。ステップS108では、特定部54は、被測定物Sの内部でX線が散乱を起こしたか否かを判定する。具体的には、特定部54は、反跳電子の初期進行方向と、散乱点T1の位置と、位置(XAj、YAj、ZAj)とを含む散乱平面上において、散乱点T1と位置(XAj、YAj、ZAj)とを結ぶ直線に対してコンプトン散乱角θを有して散乱点T1を通る直線を進路L1として設定する。特定部54は、この進路L1上または進路L1からの距離が所定の距離よりも小さい位置にX線源2が存在するか否か判定する。被測定物Sの内部でX線が散乱していないと判定された場合には、ステップS108が肯定判定されて、ステップS109へ進む。一方、被測定物Sの内部でX線が散乱したと判定された場合には、ステップS108が否定判定されて、後述するステップS110へ進む。
 ステップS109では、特定部54は、チャンバー411Aに入射したX線がコンプトン散乱を起こさなかったと仮定した場合に入射することが予想される吸収部42A上の位置(予想入射位置)を算出する。この場合、特定部54は、コンプトン散乱角θと、散乱点T1の位置と、進路L1とに基づいて、上記予想入射位置を算出する。画像生成部53は、吸収部42AからのエネルギーE2と、特定部54により算出された上記の予想入射位置とを内部情報として記憶してステップS110へ進む。ステップS110では、制御装置5は、検出器4Aのすべてのイベントについて特定処理を行ったか否かを判定する。全てのイベントについて特定処理が行われた場合には、ステップS110が肯定判定されてステップS111へ進む。ステップS111では、画像生成部53は、生成された内部情報を用いて、被測定物SのX線投影画像データを生成して、後述するステップS113へ進む。
 ステップS110において特定処理が行われていないイベントが存在する場合には、ステップS110が否定判定されてステップS112へ進む。ステップS112では、特定部54は、次のイベントでの散乱部41Aのピクセル電極412Aからの電気信号および吸収部42Aからの電気信号を取得してステップS104へ戻り、以後同様の処理を行う。
 ステップS113では、制御装置5は、全ての測定角度についてX線投影画像データの生成が行われたか否かを判定する。全ての測定角度についてX線投影画像データの生成が行われた場合には、ステップS113が肯定判定されてステップS114へ進む。ステップS114では、画像再構成部56は、複数のX線投影画像データを用いて3次元データを生成して処理を終了する。なお、制御装置5は、この3次元データをモニタ(不図示)に表示したり、メモリ(不図示)に記憶させたりすることができる。
 ステップS113においてX線投影画像データが生成されていない測定角度が存在する場合には、ステップS113が否定判定されてステップS115へ進む。ステップS115においては、載置台制御部52はマニピュレータ部36を駆動させて、載置台31を所定の測定角度に回転駆動させてステップS102へ戻り、以後同様の処理を行う。
 上述した第2の実施の形態によれば、次の作用効果が得られる。
(1)検出器4Aの散乱部41Aは、被測定物Sを透過した透過X線(または被測定物Sに照射されずに直進した透過X線)または被測定物Sの内部で散乱した散乱X線を含むXが入射して散乱することにより生じる電子の進行に関する情報を検出する。制御装置5の特定部54は、検出器4Aに入射して散乱したX線の散乱位置と、吸収部42Aに入射したX線の入射位置、および電子の進行に関する情報に基づいて、透過X線(または被測定物Sに照射されずに直進した透過X線)および散乱X線の少なくとも一方を特定する。これにより、X線のうち、被測定物Sで散乱した散乱X線と、散乱していないX線とを選別できる。X線投影画像データを生成する際には、ノイズの要因となる散乱X線を被測定物Sの内部情報から除外して、X線投影画像データの画質を向上させることができる。
(2)検出器4Aの散乱部41Aは、被測定物Sを透過した透過X線または被測定物Sの内部で散乱した散乱X線を含むXが入射して散乱することにより生じる反跳電子の進行方向を検出する。制御装置5の特定部54は、反跳電子の進行方向と、吸収部42Aで検出したX線の入射位置と、X線の散乱角θとに基づいて、透過X線および散乱X線の少なくとも一方を特定する。これにより、X線のうち、被測定物Sで散乱した散乱X線と、散乱していないX線とを選別することが可能になる。X線投影画像データを生成する際には、ノイズの要因となる散乱X線を被測定物Sの内部情報から除外することが可能になり、X線投影画像データの画質を向上させることができる。
(3)特定部54は、X線が散乱部41Aでコンプトン散乱することにより失ったエネルギーE1を算出し、反跳電子の進行方向、算出されたエネルギーE1および吸収部42Aにより得られたエネルギーE2に基づいて、透過X線および散乱X線の少なくとも一方を特定する。これにより、X線のうち、被測定物Sで散乱した散乱X線と、散乱していないX線とを選別することが可能になる。
(4)散乱部41は、X線がチャンバー411Aにてコンプトン散乱した散乱点T1を検出し、特定部54は、散乱点T1と、反跳電子の初期進行方向と、チャンバー411Aでコンプトン散乱したX線のコンプトン散乱角θと、吸収部42Aに入射したX線の入射位置とが満たす関係に基づいて、透過X線および散乱X線の少なくとも一方を特定する。これにより、検出器4Aに入射したX線のうち散乱X線を除外して、X線投影画像データを生成することができる。
(5)特定部54は、散乱点T1と、反跳電子の初期進行方向と、コンプトン散乱角θと、吸収部42AへのX線の入射位置とに基づいて、X線源2の存在が推定される位置(進路L1)を算出し、算出された位置にX線源2が存在すると、X線を透過X線として特定する。これにより、検出器4Aに入射したX線のうち散乱X線を除外できるので、X線投影画像データの画質向上に寄与することができる。
(6)特定部54は、散乱部41AからエネルギーE1に対応する電気信号が出力された時刻Tj+tsigと、吸収部42AからエネルギーE2に対応する電気信号が出力された時刻Tjとの時間差が所定の範囲T1以内のX線を用いて散乱X線を特定する。これにより、異なるタイミングでX線源2から出射され、検出器4Aに入射したX線を除外して散乱X線を特定できるので、散乱X線を特定する精度を向上させることができる。
(7)特定部54は、X線源2から照射され、被測定物Sの異なる位置を透過して、吸収部42Aの同一の画素に入射したX線を除外して、散乱X線を特定する。これにより、被測定物Sの内部の異なる位置の情報が、同一の位置の情報として内部情報に含まれることが抑制されるので、X線投影画像データの画質の低下を抑制できる。
(8)特定部54は、算出されたエネルギーE1と吸収部42Aにより得られたエネルギーE2との和が、X線源2から照射されたX線のエネルギーE0に基づく値を超えると、被測定物Sの異なる位置を透過して、吸収部42Aの同一の位置に入射した複数のX線光子であることを検出し、検出したX線光子を除外して散乱X線を特定する。これにより、被測定物Sの異なる位置を透過して、吸収部42Aの同一の画素に入射したX線を、被測定物Sの内部情報から除外することができる。
 上述した第2の実施の形態のX線装置100Aを、以下のように変形できる。
 上述した実施の形態においては、検出器4Aの散乱部41Aが、ガスが封入されたチャンバー411Aを有する場合を例に挙げて説明を行ったが、この例に限定されない。例えば、散乱部41Aは、SOIを用いた公知の三次元量子イメージ検出器により構成されてもよい。図12(a)にこの場合の検出器4の構成を示す。散乱部41Aは、SOIを用いた三次元量子イメージ検出器はZ方向に沿った厚みを有するSiセンサー部500を備え、このSiセンサー部500にて反跳電子の飛跡を検出する。Siセンサー部500は、シリコン層501と、シリコン層501のZ方向+側に集積回路502とを有する。被測定物SからのX線は、シリコン層501を透過またはコンプトン散乱し、集積回路502は、実施の形態と同様にしてコンプトン散乱による反跳電子の飛跡の三次元形状を検出する。
 また、散乱部41Aと吸収部42Aとが、Z方向に沿って複数の三次元量子イメージ検出器を積層した構成を有してもよい。図12(b)、(c)に、この場合の検出器4Aの構成を示す。図12(b)はYZ平面での断面図であり、図12(c)はXY平面での外観図である。散乱部41Aは、Z方向に沿って積層された複数の三次元量子イメージ検出器511a~511d(総称する場合は符号511を付す)と、処理基板512上に設けられた処理回路513とを有する。吸収部42Aは、Z方向に沿って積層された複数の三次元量子イメージ検出器511e~511h(総称する場合は符号511を付す)と、処理基板514上に設けられた処理回路515とを有する。なお、図12(b)、(c)に示す例では、散乱部41Aおよび吸収部42Aは、それぞれ4個の三次元量子イメージ検出器511a~511d、511e~511hを有する場合を示すが、散乱部41Aおよび吸収部42Aが有する三次元量子イメージ検出器511の個数は図12(b)、(c)に示す例に限定されない。また、散乱部41Aと吸収部42Aとが異なる個数の三次元量子イメージ検出器511を有してもよい。
 三次元量子イメージ検出器511は、基板520と、半導体素子521とを有する。半導体素子521は薄板状に形成され、基板520に設けられる。半導体素子521は、その上面に(Z方向-側)に複数のストリップ電極を有する。ストリップ電極はそれぞれワイヤ523により処理基板512、514上に設けられた処理回路513、515に接続される。ストリップ電極は、X線や反跳電子が入射されると電気信号を出力する。散乱部41にX線が入射しコンプトン散乱して反跳電子が出射すると、三次元量子イメージ検出器511a~511dのうち反跳電子が通過したストリップ電極から電気信号が処理回路513に出力される。散乱部41の三次元量子イメージ検出器511a~511dがZ方向に積層されていることにより、反跳電子の飛跡のZ方向の位置が検出される。ストリップ電極の配置により、各三次元量子イメージ検出器511a~511dを通過した反跳電子の軌跡のXY方向の位置が検出される。これにより、反跳電子の飛跡の三次元形状を検出することができる。吸収部42Aの三次元量子イメージ検出器511e~511hはX線が入射するとストリップ電極から電気信号を処理回路515に出力し、X線が吸収された位置が検出される。
 上述した第2の実施の形態および変形例にて説明したX線装置100Aを、第1の実施の形態において図5を参照して説明した構造物製造システム400に適用できる。この場合、構造物製造システム400は、図5に示すX線装置100に代えてX線装置100Aを有し、X線装置100Aが設計装置410の設計処理に基づいて成形装置420により作成された構造物の形状情報を取得する測定処理を行い、制御システム430の検査部432は、測定処理にて取得された形状情報と設計処理にて作成された設計情報とを比較する検査処理を行う。この結果、第1の実施の形態の構造物製造システム400にて得られた作用効果と同様の作用効果が得られる。
-第3の実施の形態-
 図面を参照しながら、第3の実施の形態によるX線装置について説明する。以下の説明では、第2の実施の形態と同じ構成要素には同じ符号を付し、相違点を主に説明する。特に説明しない点については、第2の実施の形態と同様である。第3の実施の形態においては、被測定物の内部で散乱した散乱X線を用いて、被測定物の内部情報を求める点で第2の実施の形態と異なる。以下、詳細に説明する。
 第3の実施の形態においては、図7に示す第2の実施の形態のX線装置100Aと同様の構成を有する。すなわち、第3の実施の形態のX線装置100Aは、第1の実施の形態と同様の筐体1、X線源2、載置部3および制御装置5と、図8に示す第2の実施の形態の検出器4Aを有する。
 本実施の形態においても、制御装置5の特定部54は、第2の実施の形態と同様にして、散乱部41Aおよび吸収部42Aから出力された電気信号に基づいて、X線源2からのX線が被測定物Sの内部で散乱したか否かを判定する。特定部54は、判定結果に基づいて、被測定物SからのX線のうち、被測定物Sにおいて散乱した散乱X線を特定、または散乱X線以外のX線を特定する。すなわち、特定部54は、検出器4Aに入射したX線が、被測定物Sで散乱したX線であるか否かを特定する。
 特定部54は、検出器4Aに入射したX線が散乱X線であることを特定した場合には、散乱X線が被測定物S中で散乱した位置を特定する。以下、被測定物S中でX線が散乱した位置(散乱位置)を特定する方法について詳細に説明する。
 特定部54は、第2の実施の形態の場合と同様にして、検出器4Aの散乱部41Aにてコンプトン散乱によりX線が失ったエネルギーE1と、コンプトン散乱した後に吸収部42Aに到達したX線のエネルギーE2とを加算して、チャンバー411Aに入射したX線のエネルギーE0を算出する。特定部54は、上述した式(8)や式(9)を用いて、散乱角θや方位角φを算出して、チャンバー411Aに入射したX線の進路L1を決定する(図9参照)。
 特定部54は、X線の進路L1を決定すると、被測定物Sの内部でX線が散乱した位置(以下、散乱位置と呼ぶ)を決定する。散乱位置の決定方法について、X線源2から出射されるX線が単色の場合と、X線源2から出射されるX線がエネルギー分布を有し、最大エネルギーが未知の場合と、X線源2から出射されるX線がエネルギー分布を有し、最大エネルギーが既知の場合とに分けて説明を行う。
<X線源2からのX線が単色の場合>
 X線源2が出射する単色のX線のエネルギーをEsとして説明する。この場合、被測定物Sの内部で散乱したX線の散乱角αは、X線源2が出射するX線のエネルギーEsと、チャンバー411Aに入射したX線のエネルギーE0(=E1+E2)とを用いて、以下の式(13)により表すことができる。
Figure JPOXMLDOC01-appb-M000013
 X線の進路L1と、X線源2からのX線とが散乱角αで交差する点が散乱位置である。この関係を図13に示す。なお、図13は、図9(b)と同様に、散乱部41Aの画素412aおよび412fと、吸収部42Aの画素42aとから電気信号が出力された場合を示している。図13に示すように、散乱部41Aと吸収部42Aとの出力に基づいて特性されたX線の経路L1と、X線源2からのX線の経路L11とが、散乱角αにて交差する点P1が散乱位置として決定する。
 具体的に説明する。特定部54は、吸収点T2の座標(吸収部42AのうちエネルギーE2を出力した画素の位置)と、散乱点T1の座標(ピクセル電極412Aのうち最も遅い時刻で電気信号を出力した画素の位置とZ方向の位置とにより決まる位置)と、散乱部41Aでの散乱角θと、X線源2の位置と、散乱角αとに基づいて、XYZ座標系における散乱位置P1の座標を算出する。特定部54は、散乱されたことが特定された複数の散乱X線について上記の処理を行うことにより、複数の散乱位置Pの座標を算出する。その結果、被測定物Sの内部において、散乱位置の度数分布が求まる。この度数分布により、被測定物Sの内部情報、例えば、被測定物Sの内部におけるX線の散乱率が大きい材料の分布状態等を知ることができる。また、特定部54は、算出した複数の散乱位置Pのそれぞれを、被測定物Sに対応する画像データ上に重ねた画像データを生成することができる。これにより、被測定物Sの内部情報を視覚的に表した画像を提供することができる。
 X線が被測定物Sの内部で2回以上散乱(多重散乱)した場合、上記の式(13)で得られる散乱角αと散乱位置P1とは、多重散乱のうち最も検出器4Aに近い側(Z方向+側)で散乱したときの散乱角と散乱位置や、異なる散乱角と散乱位置とを表す場合がある。このため、特定部54により生成された1つの画像データでは、散乱の分布状況の表現精度が劣化する可能性がある。このような劣化を抑制するため、特定部54は、被測定物Sを所定の回転量だけ回転させるごとに、上記と同様の処理を行い、被測定物Sの内部における散乱の分布状況を示す画像データを生成してもよい。これらの複数の画像データを加算して合成画像データを生成することで、多重散乱により実際とは異なる散乱位置のデータ値が平均化され、被測定物Sの内部における散乱の分布状況の表現精度の劣化が抑制される。
<X線源2から出射されるX線がエネルギー分布を有し、最大エネルギーが未知の場合>
 X線源2から出射されるX線がエネルギー分布を有している場合(単色でない場合)であって、かつ、X線の最大エネルギーEsmaxが未知の場合、上記の式(13)を用いて散乱角αを算出することができないため、上記の手順により散乱位置を特定することができない。このため、特定部54は、図14(a)、(b)を参照して説明する以下の方法により、被測定物Sの内部における散乱の分布状況を得る。
 図14は、X線源2と検出器4Aと被測定物SとのXZ平面における位置関係を模式的に示す図であり、図14(b)は図14(a)に示す位置関係の状態から被測定物Sを所定の回転量(時計回りに90°)だけ回転させた場合を示す。なお、特定部54により、第2の実施の形態に示した手順により算出された散乱部41Aにおける散乱点T1を、検出器4Aに入射したX線の入射点として説明する。また、被測定物Sは、一例として内部に散乱特性(散乱率)が異なる物質からなる領域S1を含むものとする。図14では、被測定物Sの散乱特性が既知のように図示されているが、実際は被測定物Sの大きさも散乱特性も以下の手続きにより求めるものである。
 図14(a)において、入射点(散乱点T1-1、T1-2、・・・、T1-N)に入射するX線が被測定物Sで散乱した散乱X線である。
 特定部54は、入射点T1-1に入射した散乱X線の進路L1-1を、上述した式(8)を用いて算出する。算出された進路L1-1は直線で表され、進路L1-1上の何れかの位置に被測定物Sでの散乱位置が存在する。この進路L1-1を散乱点候補位置U1-1と呼ぶ。すなわち、被測定物Sでの散乱位置は、散乱点候補位置U1-1上に存在するいずれかの点である。入射点T1-1近傍に入射した別の2つの散乱X線の進路を一点破線にて図14(a)に図示する。これらは、被測定物S内において、上記の散乱位置とは異なる別の散乱位置にて散乱したX線の進路である。これらの散乱X線の進路についても、それぞれの散乱点候補位置が定まる。同様に、特定部54は、式(8)を用いて、入射点T1-2に入射した散乱X線の進路L1-2を算出する。算出された進路L1-2上の何れかの位置には、進路L1-1と同様に、被測定物Sでの散乱位置が存在する。この進路L1-2を散乱点候補位置U1-2と呼ぶと、被測定物Sでの散乱位置は、散乱点候補位置U1-2上に存在するいずれかの点である。入射点T1-2近傍に入射した別の2つの散乱X線の進路を一点破線にて図14(a)に図示する。これらは、被測定物S内において、上記の散乱位置とは異なる別の散乱位置にて散乱したX線の進路である。これらの散乱X線の進路についても、それぞれの散乱点候補位置が定まる。このようにして、複数の入射点T1-Nについて、散乱点候補位置U1-Nを特定する。以上の手順により、複数の進路(すなわち複数の散乱点候補位置)が算出されるので、それらを重畳することで、図14(a)の状態において、XYZ座標系における被測定物Sでの散乱X線の散乱点候補分布Dが得られる。
 次に、図14(b)に示すように、被測定物Sを所定の回転量Φ(例えば、時計回りに90°)だけ回転させた状態で、上記と同様の手順により、散乱X線の進路が存在する散乱点候補分布Dを求める。すなわち、特定部54は、散乱点候補位置U1-1’、散乱点候補位置U1-2’・・・散乱点候補位置U1-N’を特定する。特定部54は、それらの散乱点候補位置を重畳することにより、図14(b)の状態において、被測定物Sでの散乱X線の散乱点候補分布Dが得られる。
 同様にして、特定物54は、被測定物SをΦ、Φ・・・ΦN回転させた状態において複数の散乱点候補分布D、D・・・Dを求める。これらの散乱点候補分布DからDの共通部分を求める。この共通部分は被測定物Sの内部における散乱の分布状況(例えば、被測定物Sの内部における散乱率が大きい材料の分布領域)に相当する。共通部分は、例えば、次に説明するような手順で求めることができる。散乱点候補分布D、D、D・・・Dのそれぞれを、被測定物Sを回転させた角度だけ回転させた散乱点候補分布D’、D’、D’・・・D’を求める。次に、XYZ座標系において、散乱点候補分布D、D’、D’、D’・・・D’を加算し、周辺のハロー部分を減算して除去すると、共通部分が算出される。算出された共通部分を被測定物Sの内部における散乱の分布状況とする。特定物54は、算出した共通部分を被測定物Sに対応する画像データ上に重ねた画像データを生成することができる。これにより、被測定物Sの内部情報を視覚的に表した画像を提供することができる。
<X線源2から出射されるX線がエネルギー分布を有し最大エネルギーは既知の場合>
 この場合、X線源2から出射されるX線がエネルギー分布を有している。X線源は制動輻射によりX線を出射するため、出射X線のエネルギーEiは最大エネルギーEsmaxより小さく、最大エネルギーEsmaxは既知である。エネルギーEiのX線が被測定物Sにおいて散乱する際の散乱角βは、以下の式(14)により表される。
Figure JPOXMLDOC01-appb-M000014
 ここで、エネルギーE0は、検出器4Aの出力に基づいて算出することができる。しかし、エネルギーEiは未知なので、式(14)を用いて散乱角βを求めることはできない。このような場合、まず、X線の最大エネルギーEsmaxは既知なので、上述した式(13)のEsをEsmaxと置き換えることにより、最大エネルギーEsmaxのX線が被測定物Sで散乱する際の散乱角αを算出することができる。エネルギーEiは最大エネルギーEsmaxより小さいので、制動輻射時のエネルギーEiのX線の散乱角βは、最大エネルギーEsmaxのX線の散乱角αより小さい。
 図15は、X線源2から最大エネルギーEsmaxで出射したX線と、制動輻射時のエネルギーEiで出射したX線がそれぞれ、散乱部41Aに進路L1で入射する場合を示したものである。X線源2から最大エネルギーEsmaxで出射したX線は、被測定物Sの散乱位置P1において散乱角αで散乱して進路L1で散乱部41Aの散乱点(入射位置)T1に入射する。一方、X線源2から制動輻射時のエネルギーEiで出射したX線は、被測定物Sの散乱位置P2において散乱角βで散乱して進路L1で散乱部41Aの散乱点(入射位置)T1に入射する。従って、散乱位置P2は、散乱点(入射位置)T1と散乱位置P1とを結ぶ直線上に位置することが推定できる。
 この場合、特定部54は、以下のようにして被測定物Sの内部での散乱位置を特定する。
 図16は、図14と同様に、X線源2と検出器4Aと被測定物SとのXZ平面における位置関係を模式的に示す図であり、図16(b)は図16(a)に示す位置関係の状態から被測定物Sを所定の回転量(時計回りに90°)だけ回転させた場合を示す。図16(a)において、入射位置(散乱点T1-1、T1-2、・・・、T1-N)に入射するX線が被測定物Sで散乱した散乱X線である。図16では、被測定物Sの散乱特性が既知のように図示されているが、実際には被測定物Sの大きさも散乱特性も以下の手続きにより求めるものである。
 特定部54は、入射点T1-1に入射した散乱X線の進路L1-1を、式(8)を用いて算出する。算出された進路L1-1は直線で表され、進路L1-1上の何れかの位置に被測定物Sの散乱位置が存在する。すなわち、進路L1-1を散乱点候補位置U1-1と呼ぶと、被測定物Sでの散乱位置は、散乱点候補位置U1-1上に存在するいずれかの点である。ただし、図15を用いて説明したように、被測定物SにおけるエネルギーEiのX線の散乱位置P2は、進路L1-1上の入射位置T1-1から散乱位置P1までの範囲に位置すると推定できる。このため、散乱点候補位置U1-1は、被測定物SよりもX線源2側(Z方向-側)の領域を含まない。入射点T1-1近傍に入射した別の2つの散乱X線の進路を一点破線にて図16(a)に図示する。これらは、被測定物S内において、上記の散乱位置とは異なる別の散乱位置にて散乱したX線の進路である。これらの散乱X線の進路についても、それぞれの散乱点候補位置が定まる。同様に、特定部54は、式(8)を用いて入射点T1-2に入射した散乱X線の進路L1-2を算出する。算出された進路L1-2(すなわち、散乱点候補位置U1-2)上の何れかの位置に、被測定物Sでの散乱位置が存在する。この場合も、進路L1-1の場合と同様に、散乱点候補位置U1-2は、被測定物SよりもX線源2側(Z方向-側)の領域を含まない。入射点T1-2近傍に入射した別の2つの散乱X線の進路を一点破線にて図16(a)に図示する。これらは、被測定物S内において、上記の散乱位置とは異なる別の散乱位置にて散乱したX線の進路である。これらの散乱X線の進路についても、それぞれの散乱点候補位置が定まる。このようにして、複数の入射点T1-Nについて、散乱点候補位置U1-Nを特定する。以上の手順により、複数の進路(すなわち複数の散乱点候補位置)が算出されるので、それらを重畳することで、図16(a)の状態において、XYZ座標系における被測定物Sでの散乱X線の散乱点候補分布Dが得られる。
 次に、図16(b)に示すように、被測定物Sを所定の回転量回転量Φ(例えば、時計回りに90°)だけ回転させた状態で、上記と同様の手順により、散乱X線の進路が存在する散乱点候補分布Dを求める。すなわち、特定部54は、散乱点候補位置U1-1’、散乱点候補位置U1-2’・・・散乱点候補位置U1-N’の座標を画像データ上での座標を算出する。これにより、図14(b)の状態において、進路分布Dが得られる。
 同様にして、特定物54は、被測定物SをΦ、Φ・・・ΦN回転させた状態において複数の散乱点候補分布D、D・・・Dを求める。これらの散乱点候補分布DからDの共通部分を求める。この共通部分は被測定物Sの内部における散乱の分布状況(例えば、被測定物Sの内部における散乱率が大きい材料の分布領域)に相当する。共通部分は、例えば、次に説明するような手順で求めることができる。散乱点候補分布D、D、D・・・Dのそれぞれを、被測定物Sの回転を回転させた角度だけ回転させた散乱点候補分布D’、D’、D’・・・D’を求める。次に、XYZ座標系において、散乱点候補分布D、D’、D’、D’・・・D’を加算し、周辺のハロー部分を減算して除去すると、共通部分が算出される。なお、周辺のハロー部分の強度は、図14を用いて説明したX線源2からのX線の最大エネルギーが未知の場合に比べ、はるかに小さい。算出された共通部分を被測定物Sの内部における散乱の分布状況とする。特定物54は、算出した共通部分を被測定物Sに対応する画像データ上に重ねた画像データを生成することができる。これにより、被測定物Sの内部情報を視覚的に表した画像を提供することができる。
 なお、図14及び図16に示す例では、被測定物Sを回転させる例を示したが、これに代えて、検出器4AとX線源2とを被測定物Sの周りに回転させてもよい。
 上述した第3の実施の形態のX線装置100Aは、被測定物Sの内部で散乱した散乱X線を特定し、散乱X線に基づいて、被測定物Sの内部情報および内部情報を表す画像データを生成することができる。
 上述した第3の実施の形態にて説明したX線装置100Aを、第1の実施の形態において図5を参照して説明した構造物製造システム400に適用できる。この場合、構造物製造システム400は、図5に示すX線装置100に代えてX線装置100Aを有し、X線装置100Aが設計装置410の設計処理に基づいて成形装置420により作成された構造物の形状情報を取得する測定処理を行い、制御システム430の検査部432は、測定処理にて取得された形状情報と設計処理にて作成された設計情報とを比較する検査処理を行う。この結果、第1の実施の形態の構造物製造システム400にて得られた作用効果と同様の作用効果が得られる。
 次のような変形も本発明の範囲内であり、変形例の一つ、もしくは複数を上述の実施形態と組み合わせることも可能である。
(1)X線装置100、100Aにおいて、被測定物Sが載置される載置台31がX軸移動部33と、Y軸移動部34と、Z軸移動部35とによってX軸、Y軸およびZ軸方向に移動されるものに限定されない。載置台31はX軸、Y軸およびZ軸方向に移動せず、X線源2および検出器4をX軸、Y軸およびZ軸方向に移動させることにより、被測定物Sに対してX線源2および検出器4を相対移動させる構成をX線装置100が有してもよい。また、載置台31が回転軸Yrにて回転するものに代えて、載置台31は回転せず、X線源2と検出器4とが回転軸Yrにて回転する構成をX線装置100が有してもよい。
(2)上述した各実施の形態と変形例においては、制御装置5は、画像生成部53および画像再構成部56を有するものとして説明したが、制御装置5が画像生成部53と画像再構成部56とを有していなくてもよい。画像生成部53と画像再構成部56とが、X線装置100とは別体の処理装置等に設けられ、特定部54で特定された散乱X線または散乱X線以外のX線に関する情報を、例えばネットワークや記憶媒体等を介して取得して、X線投影画像データや3次元データを生成してもよい。 
 本発明の特徴を損なわない限り、本発明は上記実施の形態に限定されるものではなく、本発明の技術的思想の範囲内で考えられるその他の形態についても、本発明の範囲内に含まれる。
 次の優先権基礎出願の開示内容は引用文としてここに組み込まれる。
 日本国特許出願2017年第239035号(2017年12月13日出願)
 日本国特許出願2017年第245413号(2017年12月21日出願)
2…X線源、4、4A…検出器、5…制御装置、
41、41A…散乱部、42、42A…吸収部、53…画像生成部、
54…特定部、100、100A…X線装置、
400…構造物製造システム、410…計測装置、
420…成形装置、430…制御システム、440…リペア装置
 

Claims (22)

  1.  被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射され、入射した前記X線が散乱した散乱位置を検出する第1検出器と、
     前記第1検出器で散乱または透過した前記X線が入射され、入射した前記X線の入射位置を検出する第2検出器と、
     前記散乱位置および前記入射位置に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、
    を備えるX線装置。
  2.  請求項1に記載のX線装置において、
     前記特定部は、前記第1検出器で散乱した前記X線の散乱角と、前記散乱位置と、前記入射位置とが満たす関係に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  3.  請求項2に記載のX線装置において、
     前記特定部は、前記散乱角と、前記散乱位置と、前記入射位置とに基づいて、X線源の存在が推定される方向を算出し、前記算出された方向にX線源が存在した場合に、前記X線を前記透過X線として特定するX線装置。
  4.  請求項2または3に記載のX線装置において、
     前記第1検出器は、前記X線が前記第1検出器において失った第1エネルギーを検出し、
     前記第2検出器は、前記第1検出器で散乱または透過した前記X線の第2エネルギーを検出し、
     前記特定部は、前記第1エネルギーおよび前記第2エネルギーに基づいて、前記散乱角を特定するX線装置。
  5.  請求項4に記載のX線装置において、
     前記特定部は、前記第1検出器から前記第1エネルギーが出力された時刻と、前記第2検出器から前記第2エネルギーが出力された時刻との差が所定の時間差以内の前記X線に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  6.  請求項4または5に記載のX線装置において、
     前記特定部は、前記第1検出器からの前記第1エネルギーが出力された時刻と、前記被測定物の異なる位置を透過して、前記第1検出器の別の位置において前記第1エネルギーが出力された時刻との差が所定の時間差以内の、前記検出したX線を除外して、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  7.  請求項1から6のいずれか一項に記載のX線装置であって、
     前記第1検出器は、入射した前記X線の散乱により生じる電子の進行に関する情報を検出し、
     前記特定部は、前記散乱位置、前記入射位置、および、前記第1検出器で検出された前記電子の進行に関する情報に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  8.  請求項7に記載のX線装置であって、
     前記特定部は、前記散乱位置および前記入射位置に基づいて、前記第1検出器に入射して散乱した前記X線の進行に関する情報を取得し、前記電子の進行に関する情報と、前記X線の進行に関する情報とに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  9.  請求項7または8に記載のX線装置において、
     前記特定部は、前記散乱位置、前記入射位置、および、前記電子の進行に関する情報に基づいて、X線源の存在が推定される方向を算出し、前記算出された方向にX線源が存在した場合、前記X線を前記透過X線として特定する、X線装置。
  10.  請求項7から9のいずれか一項に記載のX線装置において、
     前記特定部は、前記第1検出器で散乱した前記X線の散乱角および前記電子の出射角の少なくとも一方に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  11.  請求項10に記載のX線装置において、
     前記第1検出器は、前記X線が前記第1検出器において失った第1エネルギーを検出し、
     前記第2検出器は、前記第1検出器で散乱または透過して入射した前記X線の第2エネルギーを検出し、
     前記特定部は、前記第1エネルギーおよび前記第2エネルギーに基づいて、前記散乱角および前記出射角の少なくとも一方を算出するX線装置。
  12.  請求項11に記載のX線装置において、
     前記特定部は、前記第1検出器で前記X線が散乱した時刻と、前記第2検出器から前記第1エネルギーが出力された時刻との差が所定の時間差以内の前記X線に基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する、X線装置。
  13.  請求項3に従属する請求項4、5、11、12のいずれか一項に記載のX線装置において、
     前記特定部は、前記第1エネルギーと前記第2エネルギーとの和が、前記X線源から照射された前記X線のエネルギーに基づいて設定された所定の値を超えた前記検出したX線を除外して前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  14.  請求項7から12のいずれか一項に記載のX線装置において、
     前記特定部は、前記第1検出器で前記X線が散乱された時刻と、前記被測定物の異なる位置を透過して、前記第1検出器の別の位置において前記電子が検出された時刻との差が所定の時間差以内の、前記検出したX線を除外して、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  15.  請求項7から12のいずれか一項に記載のX線装置において、
     前記進行に関する情報は、前記電子または前記X線の進路および進行方向の前記少なくとも一方を含むX線装置。
  16.  請求項1から15の何れか一項に記載のX線装置において、
     前記特定部は、X線源から照射され、前記被測定物の異なる位置を透過して、前記第2検出器の同一の画素に所定の時間内に入射した前記X線を除外して、前記透過X線および前記散乱X線の少なくとも一方を特定するX線装置。
  17.  請求項1から16のいずれか一項に記載のX線装置において、
     前記特定部により特定された前記透過X線に基づいて前記被測定物の内部情報を取得し、前記被測定物の内部情報に基づいて、X線投影画像を生成する画像生成部をさらに備えるX線装置。
  18.  請求項1から16のいずれか一項に記載のX線装置において、
     前記特定部により特定された前記散乱X線に基づいて前記被測定物の内部情報を取得し、前記被測定物の内部情報に基づいて、X線散乱点分布像を再構成する画像生成部をさらに備えるX線装置。
  19.  被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射して散乱することにより失う第1エネルギーを検出する第1検出器と、
     前記第1検出器で散乱または透過した前記X線の第2エネルギーを検出する第2検出器と、
     前記第1エネルギーおよび前記第2エネルギーに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、
    を備えるX線装置。
  20.  被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線が入射した後に散乱を起こすことにより生じる電子の進行に関する情報と、前記散乱を起こした後の前記X線の進行に関する情報とを検出する検出器と、
     前記電子の進行に関する情報と、前記X線の進行に関する情報とに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、
     を備えるX線装置。
  21.  被測定物を透過した透過X線または前記被測定物の内部で散乱した散乱X線を含むX線を検出する検出器と、
     前記検出器への前記X線の入射方向を推定し、推定した前記入射方向とX線源の位置とに基づいて、前記透過X線および前記散乱X線の少なくとも一方を特定する特定部と、
    を備えるX線装置。
  22.  構造物の形状に関する設計情報を作成し、
     前記設計情報に基づいて前記構造物を作成し、
     作成された前記構造物の形状を、請求項1から21のいずれか一項に記載のX線装置を用いて計測して形状情報を取得し、
     前記取得された前記形状情報と前記設計情報とを比較する構造物の製造方法。
PCT/JP2018/045516 2017-12-13 2018-12-11 X線装置および構造物の製造方法 WO2019117147A1 (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2017-239035 2017-12-13
JP2017239035A JP2021028575A (ja) 2017-12-13 2017-12-13 X線装置および構造物の製造方法
JP2017-245413 2017-12-21
JP2017245413A JP2021028576A (ja) 2017-12-21 2017-12-21 X線装置および構造物の製造方法

Publications (1)

Publication Number Publication Date
WO2019117147A1 true WO2019117147A1 (ja) 2019-06-20

Family

ID=66819694

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/045516 WO2019117147A1 (ja) 2017-12-13 2018-12-11 X線装置および構造物の製造方法

Country Status (1)

Country Link
WO (1) WO2019117147A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7437337B2 (ja) 2021-03-02 2024-02-22 株式会社日立製作所 内部状態画像化装置および内部状態画像化方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080224061A1 (en) * 2007-03-13 2008-09-18 Smith Bruce D Compton Camera Configuration and Imaging Method
JP2015138027A (ja) * 2014-01-23 2015-07-30 公益財団法人高輝度光科学研究センター コンプトン散乱を用いた元素濃度の決定方法
JP5991519B2 (ja) * 2012-07-13 2016-09-14 学校法人早稲田大学 コンプトンカメラ
JP2017146102A (ja) * 2016-02-15 2017-08-24 独立行政法人石油天然ガス・金属鉱物資源機構 放射線位置特定装置、放射線位置特定方法、及び放射線位置特定用プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080224061A1 (en) * 2007-03-13 2008-09-18 Smith Bruce D Compton Camera Configuration and Imaging Method
JP5991519B2 (ja) * 2012-07-13 2016-09-14 学校法人早稲田大学 コンプトンカメラ
JP2015138027A (ja) * 2014-01-23 2015-07-30 公益財団法人高輝度光科学研究センター コンプトン散乱を用いた元素濃度の決定方法
JP2017146102A (ja) * 2016-02-15 2017-08-24 独立行政法人石油天然ガス・金属鉱物資源機構 放射線位置特定装置、放射線位置特定方法、及び放射線位置特定用プログラム

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7437337B2 (ja) 2021-03-02 2024-02-22 株式会社日立製作所 内部状態画像化装置および内部状態画像化方法

Similar Documents

Publication Publication Date Title
JP6167182B2 (ja) エックス線を利用した薄膜層の厚さ測定方法及び装置
US7130375B1 (en) High resolution direct-projection type x-ray microtomography system using synchrotron or laboratory-based x-ray source
US7141790B2 (en) Defect inspection instrument and positron beam apparatus
US20200041429A1 (en) Nondestructive inspection apparatus and method for micro defect inspection of semiconductor packaging using a plurality of miniature x-ray tubes
US20080240344A1 (en) X-ray tomosynthesis device
JP5583993B2 (ja) 小型多焦点x線源、x線回折イメージングシステム、及び小型多焦点X線源を製作するための方法
EP3109625B1 (en) Combination of cad data and computer tomography
US10285252B2 (en) Dual-energy ray scanning system, scanning method and inspecting system
JP2007309685A (ja) 検査装置及び検査方法
US11448603B1 (en) Methods and apparatuses for microscopy and spectroscopy in semiconductor systems
Syryamkin et al. Digital X-ray Tomography: Edited by VI Syryamkin
JP2015121547A (ja) 単一エネルギガンマ線源による、同位体識別、分析およびイメージングのための二重同位体ノッチ観測機
WO2019117147A1 (ja) X線装置および構造物の製造方法
KR20140059012A (ko) 비파괴 검사 장치
Nachtrab et al. Laboratory X-ray microscopy with a nano-focus X-ray source
JP7340476B2 (ja) 放射線計測装置および放射線計測方法
JP2021028576A (ja) X線装置および構造物の製造方法
JP2021026926A (ja) 画像生成方法、非一時的コンピューター可読媒体、及びシステム
US20140341342A1 (en) Equipment for the radiography of a load for performing detection, and associated method
JP2021028575A (ja) X線装置および構造物の製造方法
JP6729801B2 (ja) 検査装置、検査方法および検査対象物の製造方法
Tran et al. Monte Carlo simulation of a cone-beam CT system for lightweight casts
JP6281229B2 (ja) X線源、x線装置、構造物の製造方法、及び構造物製造システム
WO2023189135A1 (ja) 検査装置及び検査方法
JP2009008560A (ja) 非破壊検査方法及び装置

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: 18888025

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18888025

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP