US20180024079A1 - Image reconstruction method for x-ray measuring device, structure manufacturing method, image reconstruction program for x-ray measuring device, and x-ray measuring device - Google Patents

Image reconstruction method for x-ray measuring device, structure manufacturing method, image reconstruction program for x-ray measuring device, and x-ray measuring device Download PDF

Info

Publication number
US20180024079A1
US20180024079A1 US15/549,010 US201515549010A US2018024079A1 US 20180024079 A1 US20180024079 A1 US 20180024079A1 US 201515549010 A US201515549010 A US 201515549010A US 2018024079 A1 US2018024079 A1 US 2018024079A1
Authority
US
United States
Prior art keywords
measurement object
image
estimated
ray
data
Prior art date
Legal status (The legal status 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 status listed.)
Abandoned
Application number
US15/549,010
Other languages
English (en)
Inventor
Atsushi Yamada
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nikon Corp
Original Assignee
Nikon Corp
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
Application filed by Nikon Corp filed Critical Nikon Corp
Assigned to NIKON CORPORATION reassignment NIKON CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: YAMADA, ATSUSHI
Publication of US20180024079A1 publication Critical patent/US20180024079A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/044—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using laminography or tomosynthesis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/401—Imaging image processing

Definitions

  • the present invention relates to an image reconstruction method for an X-ray measuring device, a structure manufacturing method, an image reconstruction program for an X-ray measuring device, and an X-ray measuring device.
  • an X-ray device that irradiates an X-ray to an object and detects the X-ray that has passed through the object is known.
  • This X-ray device has an X-ray source from which an X-ray is irradiated and detects the X-ray that has passed through the object to enable observation of the inside thereof (cf., PTL 1). In this manner, the information about the inside of the object is acquired.
  • an object of the present invention is to provide an image reconstruction method for an X-ray measuring device, a structure manufacturing method, an image reconstruction program for an X-ray measuring device, and an X-ray measuring device, which can prevent a decrease in the accuracy of inspection.
  • an image reconstruction method for an X-ray measuring device comprises: causing a measurement object to be irradiated with X-rays from multiple different irradiation directions by rotating the measurement object and/or changing the irradiation directions of the X-rays, detecting X-rays that passed through the measurement object in each of the multiple different irradiation directions, and generating multiple sets of detection data relating to the intensity of the X-rays that passed through the measurement object; estimating the intensity of the X-rays that are assumed to have been passed through an estimated structure, whose shape has been estimated based on shape information about the measurement object, in each of the multiple different irradiation directions if the estimated structure is assumed to be irradiated with the X-rays from the multiple different irradiation directions, and generating multiple sets of estimated data; and, extracting differential data indicating a difference between the estimated data and the detection data from the detection data using, for each of the multiple sets of detection data
  • an image reconstruction method for an X-ray measuring device comprises: causing a measurement object to be irradiated with X-rays from multiple different irradiation directions by rotating the measurement object and/or changing the irradiation directions of the X-rays, detecting X-rays that have been passed through the measurement object in each of the multiple different irradiation directions, and generating multiple sets of detection data relating to the intensity of the X-rays passed through the measurement object; estimating the intensity of the X-rays that are assumed to be passed through an estimated structure in each of the multiple different irradiation directions of the X-rays if the X-rays are assumed to be irradiated to the estimated structure of which a shape has been estimated based on shape information of the measurement object along each of the multiple different irradiation directions, and generating a plurality of estimated data; extracting differential data indicating a difference between the estimated data and the detection data from the detection data using
  • a structure manufacturing method comprises: generating design information relating to the shape of a structure; forming the structure based on the design information; measuring the shape of the formed structure by performing the image reconstruction method for an X-ray measuring device according to any one of the aspects described above to acquire shape information; and comparing the acquired shape information with the design information.
  • an X-ray measuring device comprises: an acquisition unit that acquires multiple sets of detection data relating to the intensity of an X-rays passed through a measurement object, the detection data being generated by causing the measurement object to be irradiated with X-rays from multiple different irradiation directions by rotating the measurement object and/or changing the irradiation directions of the X-rays, and detecting X-rays that have been passed through the measurement object in each of the multiple different irradiation directions; and an image generation unit that estimates the intensity of the X-rays that are assumed to be passed through an estimated structure, whose shape has been estimated based on shape information about the measurement object, in each of the multiple different irradiation directions if the estimated structure is assumed to be irradiated with the X-rays from the multiple different irradiation directions, and generates multiple sets of estimated data and extracting differential data indicating a difference between the estimated data and the detection data from the detection data using, for
  • FIG. 1 is a diagram explaining the configuration of an X-ray device according to an embodiment of the present invention.
  • FIG. 2 is a diagram schematically showing a concept for generating a three-dimensional image
  • FIG. 3 is a diagram schematically showing a concept for generating a three-dimensional image using various types of data
  • FIG. 4 is a diagram schematically showing a concept of correcting the shape of an estimated structure
  • FIG. 5 is a flowchart explaining a process for generating a three-dimensional image according to the first and second embodiments
  • FIG. 6 is a flowchart explaining a process for generating a three-dimensional image according to the third embodiment
  • FIG. 7 is a block diagram showing the configuration of a structure manufacturing system according to the fourth embodiment.
  • FIG. 8 is a flowchart explaining the operation of a structure manufacturing system according to the fourth embodiment.
  • FIG. 9 is a diagram explaining the overall structure of a device that is used for providing a program product.
  • the X-ray device irradiates an X-rays to a measurement object and detects the X-rays that passed through the measurement object to acquire information about the inside (for instance, inside structure) of the measurement object or the like non-destructively.
  • the target of measurement i.e., the measurement object is a part for industrial use, for instance, a machine part or an electronic component
  • the X-ray device is called an industrial X-ray CT (Computed Tomography) inspection device that inspects the part for industrial use.
  • the inspection target of the measurement object may be an animate being such as a human body, an animal, a plant, etc.
  • the target may be one of tissues of a living organism.
  • the target may be a combination of an animate being and a non-living material such as a part for industrial use.
  • the present embodiment is presented in order to concretely explain the gist of the present invention and unless indicated specifically, it does not limit the present invention.
  • FIG. 1 is a diagram showing an example of the configuration of an X-ray device 100 according to the present embodiment. Note that for convenience of explanation, a coordinate system consisting of an X-axis, a Y-axis and a Z-axis is set as shown in FIG. 1 .
  • the X-ray device 100 includes a housing 1 , an X-ray source 2 , a mounting unit 3 , a detector 4 , a control unit 5 , a display monitor 6 , and a frame 8 .
  • the housing 1 is disposed on a floor surface of factory or the like so as to be substantially parallel (horizontal) to an XZ-plane and in the inside thereof are housed the X-ray source 2 , the mounting unit 3 , the detector 4 , and the frame 8 . To prevent the X-rays from leaking to the outside, the housing 1 is made of a material that contains lead.
  • the X-ray source 2 radiates X-rays extending along an optical axis Zr that is parallel to the Z axis, spreading in the form of a cone (so-called cone beam) toward the positive direction of the Z axis with an output point Q shown in FIG. 1 being at the top of the cone.
  • the output point Q corresponds to a focal spot of the X-ray source 2 . That is, the optical axis Zr connects the output point Q, which is the focal spot of the X-ray source 2 , with the center of an imaging area of the detector 4 described later.
  • the X-ray source 2 that radiates X-rays in a fan-like form (so-called fan beam) and one that radiates X-rays in a form of a line (so-called pencil beam) instead of one that emits X-rays in a cone form are also included as the aspects of the present invention.
  • the X-ray source 2 can radiate at least one of, for instance, an ultrasoft X-ray of about 50 eV, a soft X-ray of about 0.1 to 2 keV, an X-ray of about 2 to 20 keV, and a hard X-ray of about 20 to 100 keV.
  • the X-ray source 2 may radiate an X-ray of, for instance, 1 to 10 Mev.
  • the mounting unit 3 includes a mount stage 30 on which a measurement object S is mounted and a manipulator unit 36 that includes a rotary drive unit 32 , a Y-axis movement unit 33 , an X-axis movement unit 34 and a Z-axis movement unit 35 and is provided on the positive (+) side along the Z-axis relative to the X-ray source 2 .
  • the mount stage 30 is provided rotatable by the rotary drive unit 32 and moves together with the manipulator unit 36 when the manipulator unit 36 moves in the direction along the X axis, Y axis, and Z axis respectively.
  • the rotary drive unit 32 is structured so as to include, for instance, an electric motor that rotates the mount stage 30 by rotative force generated by the electric motor being driven under control of the control unit 5 as described later such that the mount stage 30 rotates about an axis that is parallel to the Y-axis and that passes the center of the mount stage 30 as a rotation axis Yr. That is, the rotary drive unit 32 changes relative direction of the mount stage 30 and the measurement object S thereon, with respect to the X-rays radiated from the X-ray source 2 by rotating the mount stage 30 .
  • the Y-axis movement unit 33 , the X-axis movement unit 34 , and the Z-axis movement unit 35 under control by the control unit 5 , move the mount stage 30 along the X-axis direction, along the Y-axis direction, and along the Z-axis direction, respectively, so that the measurement object S can be positioned within the area in which the X-rays outputted from the X-ray source 2 is irradiated.
  • the Z-axis movement unit 35 under control by the control unit 5 , moves the mount stage 30 along the Z-axis direction to a position at which the distance of the measurement object S from the X-ray source 2 is a distance such that a projection image of the measurement object S has a desired magnification ratio.
  • AY position detector 331 an X position detector 341 , and a Z position detector 351 are encoders that detect the positional movements of the mount stage 30 along the X-axis direction, along the Y-axis direction, and along the Z-axis direction, respectively, by the Y-axis movement unit 33 , the X-axis movement unit 34 , and the Z-axis movement unit 35 , respectively, and output signals that indicate the detected position (hereafter, referred to as detected positional movements) to the control unit 5 .
  • the rotational position detector 321 is an encoder that detects the rotational position of the mount stage 30 that is rotated about the rotation axis Yr by the rotary drive unit 32 and outputs a signal that indicates the detected rotational position (hereafter, referred to as detected rotational position) to the control unit 5 . That is, the detected rotational position indicates a relative direction of the measurement object S on the mount stage 30 with respect to the X-rays radiated from the X-ray source 2 .
  • the detector 4 is provided on the positive (+) side of Z-axis relative to the X-ray source 2 and the mount stage 30 . That is, the mount stage 30 is provided between the X-ray source 2 and the detector 4 along the Z-axis direction.
  • the detector 4 has an incident surface 41 that is parallel to the XY-plane. At the incident surface 41 enter the X-rays which include the X-rays that have been irradiated from the X-ray source 2 and passed through the measurement object S mounted on the mount stage 30 .
  • the detector 4 is constituted by a scintillator unit that contains a scintillation material that is known in the art, a photomultiplier tube, a light receiving unit, and so on.
  • the detector 4 converts energy of the X-rays incident to the incident surface 41 of the scintillator unit into light energy of visible light or ultraviolet light, amplifies the light energy with the photomultiplier tube, converts the amplified light energy into electric energy with the light receiving unit, and outputs the electric energy as electric signals to the control unit 5 .
  • the scintillator unit, the photomultiplier tube and the light receiving unit are respectively have structures in which divided into a plurality of pixels. These pixels are arranged in a two-dimensional array. Due to this, it is possible to obtain, at one time, the intensity distribution of the X-rays that has passed through the measurement object S after radiating from the X-ray source 2 .
  • the detector 4 may be one that converts the energy of the incident X-rays into electric energy and then outputs the converted electric energy in the form of electric signal without converting the energy of the incident X rays into light energy.
  • the detector 4 is not limited to one in which the pixels are arranged in a two-dimensional array.
  • the detector 4 has an incident surface 41 that extends on a plane parallel to the XY-plane, for instance, extends to the X direction.
  • the incident surface 41 may be constituted by a line sensor that includes each of pixels disposed one after another in the Y direction.
  • the direction along which the pixels of the line sensor are disposed is not limited to the Y direction but they may be arranged in the X direction or in the Z direction.
  • the detector 4 may have a structure in which no photomultiplier tube is provided and the scintillator unit is directly formed on the light receiving unit (photoelectric conversion unit).
  • the frame 8 supports the X-ray source 2 , the mounting unit 3 , and the detector 4 .
  • the frame 8 is produced to have sufficient rigidity. Due to this, it is possible to stably support the X-ray source 2 , the mounting unit 3 , and the detector 4 while a projection image of the measurement object S is being acquired.
  • the frame 8 is supported by a vibration isolation mechanism 81 to prevent transmitting the vibration that is generated in the outside to the frame 8 directly without being attenuated.
  • the control unit 5 which includes a microprocessor and its peripheral circuits and so on, reads in a control program that is stored in advance in a storage medium (for instance, a flash memory or the like) (not shown) and executes the program to thereby control each of units of the X-ray device 100 .
  • the control unit 5 includes an X-ray control unit 51 , a movement control unit 52 , an image generation unit 53 , and a work memory 55 .
  • the X-ray control unit 51 controls the operation of the X-ray source 2 .
  • the movement control unit 52 controls the movement operation of the manipulator unit 36 .
  • the image generation unit 53 performs image reconstruction processing based on electric signals that are outputted from the detector 4 in correspondence to the intensity distribution of the X-rays that passed through the measurement object S at every predetermined angle as the measurement object S is being rotated according to the rotation of the mount stage 30 to generate a three-dimensional image of the measurement object S.
  • the image generation unit 53 performs, as the image reconstruction processing, processing of generating X-ray projection image data (detection data) of the measurement object S based on the electric signals outputted from the detector 4 , and processing a three-dimensional image showing an inside structure (cross-section structure) of the measurement object S by generating a back projection image based on the X-ray projection image data of the measurement object S under the different directions of projection, that is, different detected rotational positions.
  • Examples of the processing of generating a back projection image include a back projection method, a filtered back projection method, a successive approximation method, and so on.
  • the work memory 55 is constituted by, for instance, a volatile storage medium, at which the X-ray projection image data generated by the image generation unit 53 is temporarily stored.
  • FIG. 2 is a conceptual diagram schematically showing the processing of generating a three-dimensional image according to the present embodiment.
  • FIG. 2 ( a ) schematically shows a stereoimage of the measurement object S generated based on the detection data, which is generated based on the electric signals outputted from the detector 4 .
  • the stereoimage of the measurement object S can be obtained using plurality of sets of detection data which are obtained based on the X-rays that passed through the measurement object S by irradiation the X-rays from different directions as the mount stage is rotating.
  • the stereoimage of the measurement object S is based on the transmission intensity distribution of the X-rays that passed through the measurement object S.
  • FIG. 2 ( a ) shows a stereoimage of a defect S 1 in the inside of the measurement object S as an example of the error.
  • FIG. 3 ( a ) schematically shows the transmission intensity of X-ray that is used for generating the stereoimage of the measurement object S, that is, a specified single set of detection data D 1 among plurality of sets of detection data.
  • FIG. 3 ( a ) shows the transmission intensity of the X-ray at a cross-sectional surface (a plane P indicated in broken line in FIG. 2 ( a ) ) that passes through the defect S 1 in the inside of the measurement object S when the mount stage 30 is in one of a moved and rotated position for detection.
  • the intensity in the vertical axis of FIG. 3 ( a ) shows values calculated by dividing the transmission intensity of the X-ray by the transmission length (distance) of the X-ray. The same is true for FIG. 3 ( b ) and FIG. 3 ( c ) .
  • FIG. 3 ( a ) shows an example in which the transmission intensity of the X-ray that passed through the defect S 1 is higher than the intensity of the X-ray that did not passed through the measurement object S, because a ratio of X-ray absorbed at the defect S 1 is smaller than the ratio of X-ray absorbed at the surrounding of the defect S 1 .
  • FIG. 2 ( b ) schematically shows a stereoimage of an imaginary estimated structure S 2 that is based on a shape information such as CAD or the like as design information and a material information of the measurement object S.
  • the stereoimage of the estimated structure S 2 corresponds to an ideal state in which the measurement object S is manufactured according to the design values. Therefore, the estimated structure S 2 does not contain the defect S 1 .
  • FIG. 3 ( b ) schematically shows the transmission intensity of the X-ray which is estimated to be detected if the X-rays passes through the estimated structure S 2 , that is, estimated data D 2 .
  • the estimated data D 2 is an information about the intensity of the X-ray that is assumed to be passed through the estimated structure S 2 that is acquired by simulation for the case in which the X-rays are irradiated to the estimated structure S 2 , which corresponds to that the measurement object S is manufactured according to design values (in an ideal state). Since the estimated structure S 2 does not contain the defect S 1 , the transmission intensity of the X-ray which is estimated to be passed through the estimated structure S 2 is influenced only by the absorption coefficient of the estimated structure S 2 as shown in FIG. 3 ( b ) . Note that the estimated data D 2 is generated under the same condition as the condition under which the detection data D 1 is obtained.
  • the estimated data D 2 is generated on the presumption that X-rays are irradiated to the estimated structure S 2 at the moved and rotated position for detection of the mount stage 30 as that are the same as the moved and rotated position for detection of the mount stage 30 when the detection data D 1 is generated. Moreover, the estimated data D 2 is generated on the presumption that X-rays are irradiated to the estimated structure S 2 at an output that is the same as the output of the X-rays that are outputted from the X-ray source 2 when the detection data D 1 is generated. Therefore, each of the estimated data D 2 corresponds to each of the detection data D 1 , that are detected for respective different irradiation directions of the X-ray according to the measurement object S.
  • the differential data D 3 which is the difference between the detection data D 1 and the estimated data D 2 .
  • the differential data D 3 indicates an absolute value of intensity.
  • the differential data D 3 contains the defect S 1 , such as the cavity, or the like, inside the measurement object S and an error in shape between the measurement object S and the estimated structure S 2 .
  • the differential data D 3 is obtained by extracting the difference between each of the detection data D 1 and each of the estimated data D 2 , which are in correspondence to each other. Therefore, each of the differential data D 3 that correspond to the detection data D 1 , respectively, detected for respective irradiation directions of the X-ray according to the measurement object S.
  • the extracted stereoimage of the defect S 1 and the extracted stereoimage of the estimated structure S 2 are combined with each other.
  • the stereoimage of the estimated structure S 2 corresponds to the shape manufactured according to the design values and to the ideal state that does not contain any defect or the like.
  • a stereoimage which is a combination of the stereoimage of the estimated structure S 2 in which generation of artifacts is prevented with the defect S 1 that is extracted from the stereoimage of the measurement object S. That is, the stereoimage of the defect S 1 and the stereoimage that has substantially the same shape as the measurement object S and of which generation off artifacts is prevented or reduced are generated.
  • plurality of calculated differential data D 3 are subjected to back projection to generate a back projection image Im 1 relative to the difference between the measurement object S and the estimated structure S 2 .
  • the back projection image Im 1 and an image Im 2 that corresponds to the estimated structure S 2 are combined with each other to generate an image Im 3 .
  • the image Im 3 is a three-dimensional image in which the defect or the like that is present inside the actual measurement object S is reproduced in the image Im 2 of the estimated structure S 2 , which has substantially the same shape as the measurement object S and of which generation of artifacts is prevented or reduced. Therefore, in the image Im 3 , generation of artifacts by beam hardening is prevented or reduced.
  • the estimated data D 2 is calculated as supposing that the image obtained as the image Im 3 is treated as the image of the estimated structure S 2 .
  • differential data D 3 is extracted. Using the extracted differential data D 3 , the stereoimage of the defect S 1 is generated.
  • the number of times of repetition may be determined based on the differential data D 3 that is extracted. Note that in case a line sensor is used as the detector 4 , back projection may be performed using X-ray intensity data groups in correspondence to different positions of the measurement object S along the Y direction to generate a three-dimensional image showing the inside structure of the overall measurement object S.
  • the image generation unit 53 generates the estimated data D 2 using the expressions (1) and (1)′ and extracts the differential data D 3 from the detection data D 1 based on the output from the detector 4 .
  • ⁇ 0(X,Y,Z, ⁇ ) is an attenuation coefficient of the X-ray at a position (X,Y,Z) when supposing the X-ray having a photon energy ⁇ passes through the inside of the imaginary estimated structure S 2 that has been estimated based on the design information of the measurement object S.
  • ⁇ (X,Y,Z, ⁇ ) in the expression (1) is a attenuation coefficient of the X-ray having a photon energy ⁇ at a position (X,Y,Z) when the X-ray passes through the inside of the measurement object S.
  • expression (2) The relationship between these attenuation coefficients is represented by expression (2).
  • ⁇ (X,Y,Z) is an difference between both attenuation coefficients.
  • the expressions (1) and (1)′ are derived from expression (2) as follows.
  • the attenuation amount of the X-ray when the X-ray travels over a minute portion having a distance ⁇ l in the measurement object S is represented by expression (3) as below.
  • the X-ray having an intensity Ii ( ⁇ ) that has entered the minute portion having a distance ⁇ l has an intensity represented by expression (4) at a point that the X-ray has been passed over a minute portion having a distance ⁇ l.
  • the intensity yi of the X-ray that enters the detector 4 after passed through the inside of the measurement object S along its traveling direction is represented by expression (5).
  • i indicates the position of one detection pixel among a plurality of detection pixels included in the detector 4 .
  • yi represents the intensity of the X-ray that enters an i-th detection pixel from a beginning.
  • ri represents a noise component such as dark current contained in the outputted from the detection pixel arranged at i-th position in the detector 4 . That is, according to expression (5), intensities yi in numbers in correspondence to the number of detection pixels in the detector 4 are generated for each detected rotational position.
  • the attenuation coefficient ⁇ 0(X,Y,Z, ⁇ ) is a attenuation coefficient of the X-ray at the position (X,Y,Z) of the inside of the estimated structure S 2 and it depends on the design information of the measurement object S.
  • the material information about the material that constitutes the measurement object S may be used.
  • the material information of the measurement object S may be, for instance, information about the ratio of materials contained in the measurement object S.
  • the material information of the measurement object S may be, for instance, information as to whether the material contained in the measurement object S is metal or nonmetal.
  • the material information of the measurement object S may be, for instance, information about elements or compounds contained in respective members that constitute the measurement object S.
  • As the design information for instance, information about an external shape and/or an internal shape of the measurement object S may be used.
  • the intensity Ii ( ⁇ ) depends on the intensities and spectra of the X-rays irradiated to the measurement object S, which are measuring conditions under which the measuring device is operated. Therefore, Bi represented by the expression (1)′ corresponds to the estimated detection intensity estimated to be detected by i-th detection pixel of the detector 4 that is calculated based on the design information, the material information, and the measuring conditions.
  • Bi corresponds to the transmission intensity of the X-ray that has passed through the estimated structure S 2 when the X-rays are irradiated to the estimated structure S 2 according to the measuring conditions.
  • ⁇ (X,Y,Z) is a value that is based on a difference between a distance along which the X-ray that has passed through the inside of the measurement object S travels and a distance along which the X-ray that has passed through the inside of the estimated structure S 2 travels at each of the plurality of detection pixels of the detector 4 . That is, ⁇ (X,Y,Z) corresponds to a difference between the detection data D 1 and the estimated data D 2 . Therefore, the image generation unit 53 calculates ⁇ (X,Y,Z) for a detected intensity that is outputted from each of the plurality of detection pixels arrayed in the detector 4 at each irradiation position of the X-ray irradiated to the measurement object S. By so doing, the image generation unit 53 extracts a plurality of sets of differential data D 3 at each irradiation position of the X-ray irradiated to the measurement object S.
  • the image generation unit 53 evaluates each of the ⁇ (X,Y,Z) values, whose number corresponds to the number of the detection pixels of the detector 4 , calculated by expression (1) for any irradiation position of the X-ray The image generation unit 53 , when all the calculated ⁇ (X,Y,Z) values do not exceed a first predetermined value, judges that the differential data D 3 corresponds to the difference caused by a cavity or the like in the inside of the measurement object S (for instance, the defect S 1 in FIG. 2 ).
  • the image generation unit 53 when at least one of the calculated ⁇ (X,Y,Z) values is equal to or exceeds the first predetermined value, judges that the error in shape between the measurement object S 1 and the estimated structure S 2 is considerably large, that is, the size of the measurement object S differs from the design value considerably.
  • the image generation unit 53 When all the calculated ⁇ (X,Y,Z) values do not exceed the first predetermined value, the image generation unit 53 performs back projection of the differential data D 3 to generate a back projection image Im 1 and combines it with the image Im 2 of the estimated structure S 2 to generate the image Im 3 . This enables one to evaluate location, shape, size and so on of defects, if any, in the inside of the measurement object S.
  • the image generation unit 53 performs, based on the ⁇ (X,Y,Z) value, correction of the estimated structure S 2 so that it becomes closer to the shape of the measurement object S.
  • the image generation unit 53 based on the corrected size of the estimated structure S 2 , obtains ⁇ (X,Y,Z) in the same procedure as described above and evaluates its value. This procedure is repeated until all the calculated ⁇ (X,Y,Z) values will be below the first predetermined value.
  • FIG. 4 schematically illustrates the concept of correction of the estimated structure S 2 .
  • a solid line indicates a detected intensity distribution L 1 of the X-ray detected by the detector 4 , that is, a detected intensity distribution of the X-ray that corresponds to the detection data D 1 , when the X-rays, that are outputted from the X-ray source 2 , is irradiated to the measurement object S having the defect S 1 in the inside thereof.
  • L 2 indicates the estimated intensity distribution of the X-ray to which the first processing was performed but the correction of the shape of the estimated structure S 2 was not performed yet.
  • L 1 is higher than L 2 in the transmission intensity distribution of the X-ray. This indicates that the transmission distance of the X-ray that transmits through the measurement object S is shorter than that of the designed size.
  • L 3 and L 4 indicate each an estimated intensity distribution of the X-ray that is estimated after the correction of the size of the estimated structure S 2 in the second and subsequent processing are performed.
  • the estimated transmission intensity of the X-ray for detection pixels included in the area ir among the plurality of the detection pixels included in the detector 4 is higher than the estimated transmission intensity of the X-ray for any one of detection pixels of the detector 4 other than those included in the area ir. This indicates that the defect S 1 in the inside of the measurement object S gives some influence on the estimated intensity of the X-ray for pixels in the area ir.
  • Each of the differences between the intensity distribution of L 1 and the intensity distributions of L 2 , L 3 , or L 4 shown in FIG. 4 correspond to ⁇ (X,Y,Z).
  • the image generation unit 53 changes the estimated transmission intensity distribution of the X-ray from L 2 in the direction indicated by an arrow A, that is, in the direction in which the intensity increases in FIG. 4 by performing correction of the shape of the estimated structure S 2 .
  • Correcting the shape of the estimated structure S 2 corresponds to changing the distance in which the X-ray transmits through the inside of the estimated structure S.
  • the image generation unit 53 changes the shape of the estimated structure S 2 in a direction in which it becomes closer to the shape of the measurement object S In the case of this example, the intensities in correspondence to the area ir for L 3 and L 4 become higher than the intensities in correspondence to a remaining area other than the area ir.
  • the transmission distance of the X-ray that transmits through the measurement object is larger than that of the designed size
  • the size of the measurement object is larger than the design value
  • the relationship of the magnitude of L 2 to that of L 1 is in reverse to the above relationship.
  • new ⁇ (X,Y,Z) is calculated as a difference in transmission intensity distributions of the X-ray between L 1 and L 3 .
  • this ⁇ (X,Y,Z) value is larger than the first predetermined value, the shape of the estimated structure S 2 is further corrected based on ⁇ (X,Y,Z) value.
  • the intensity distribution of the X-ray transmits through the measurement object S and the intensity distribution of the X-ray that is estimated to be passed through the estimated structure S 2 after the correction of its shape are deemed identical to each other. That is, the value of ⁇ (X,Y,Z) that is newly calculated as a difference in transmission intensity between L 1 and L 4 is equal to or lower than the first predetermined value. This means that the estimated data D 2 reproduces the transmission intensity distribution of the X-ray that has been passed through the measurement object S.
  • the image generation unit 53 may be configured to align on data the measurement object S represented by the detection data D 1 and the estimated structure S 2 represented by the estimated data D 2 . For instance, the image generation unit 53 detects the position of a boundary (edge) (pixel position i) between the measurement object S represented by the detection data D 1 and the background. The image generation unit 53 detects the pixel position i+a in correspondence to the boundary on the estimated data D 2 . The image generation unit 53 performs calculation by expression (1) using the detection data D 1 at pixel i (that is, yi) and the estimated data D 2 at pixel i+a (that is, Bi+a). As a result, at time of the first processing, the value of ⁇ (X,Y,Z) can be calculated as a relatively small value and hence the number of times of processing can be decreased even when subsequently correction of ⁇ (X,Y,Z) may become necessary.
  • the image generation unit 53 performs the above-mentioned processing on each of the plurality sets of the detection data D 1 and each of the plurality of sets of the estimated data D 2 , each are based on different projection directions of X-ray toward the measurement object S, that is, each are based on different detected rotational position to effect alignment and subsequently synthesizes the back projection image Im 1 with the image Im 2 relating to the estimated structure S 2 to generate the image Im 3 .
  • the image generation unit 53 commands the display monitor 6 to display the generated image Im 3 in a three-dimensional image.
  • Step S 1 the shape information and the material information, which constitute design information of the measurement object S, and the measuring conditions for measuring the measurement object S, which is measuring device information, are acquired, and the procedure proceeds to Step S 2 .
  • Step S 2 Bi for each of the pixels is calculated based on the shape information, the material information, and the measuring information using the above-mentioned expression (1)′ (that is, estimated data D 2 is generated) and the procedure proceeds to Step S 3 .
  • Step S 3 transmission images for all the irradiation directions are acquired and the procedure proceeds to Step S 4 .
  • the transmission image in all measuring directions is acquired for using in image reconfiguration. For instance, in FIG. 1 , the transmission images that are acquired when the measurement object S is rotated by one turn of 360° with the stage 30 .
  • Step S 4 a detected value detected at each of the pixels by causing X-rays to be actually passed through the measurement object S and calculated Bi are substituted into expression (1) and ⁇ (X,Y,Z) value is calculated according to the back projection method (that is, differential data D 3 is extracted) and the procedure proceeds to Step S 5 .
  • Step S 5 it is determined whether all the position-dependent ⁇ (X,Y,Z) values that are calculated for the detection pixels, respectively, are equal to or smaller than the first predetermined value.
  • Step S 5 is determined to be affirmative and the procedure proceeds to Step S 7 .
  • Step S 5 is determined to be negative and the procedure proceeds to Step S 6 .
  • Step S 6 the shape of the estimated structure S 2 is corrected and the procedure proceeds to Step S 2 .
  • the shape of the estimated structure S 2 that has been corrected in Step S 6 is used as design information (shape information) of the measurement object S.
  • Step S 7 since the back projection image Im 1 of the differential data D 3 is 0, the shape of the estimated structure S 2 , correction of which has thus far been continued, is deemed to match the shape of the measurement object S (structure).
  • the image of the estimated structure S 2 is displayed on the display monitor 6 as a three-dimensional image of the measurement object S and the processing is terminated. Note that the image to be displayed on the display monitor 6 is not limited to the image Im 3 of the estimated structure S 2 .
  • the image generation unit 53 generates a plurality of sets of detection data D 1 relating to the transmission intensity of the X-ray that has been passed through the measurement object S and a plurality of sets of estimated data D 2 relating to the estimated transmission intensity of the X-ray on the assumption that the X-rays have been passed through an imaginary estimated structure S 2 that is configured based on design information under the same irradiation condition as that under which the X-rays have been passed through the measurement object S.
  • the image generation unit 53 extracts differential data D 3 , which indicates the difference between the detection data D 1 and the estimated data D 2 .
  • the image generation unit 53 performs back projection of the extracted differential data D 3 to generate a back projection image Im 1 relating to the difference and combines the back projection image Im 1 with the estimated structure S 2 to generate the image Im 3 relating to the inside structure of the measurement object S. Therefore, the image Im 3 , of which generation of artifacts due to beam hardening is prevented or more reduced as compared with the back projection image obtained by performing back projection of the detection data D 1 , is generated. As a result, it becomes easier to grasp the inside defect of the measurement object S by the image Im 3 .
  • the image generation unit 53 generates estimated data D 2 by estimating the transmission intensity of the X-ray when the X-rays are irradiated to the estimated structure S 2 based on the material information of the measurement object S and spectrum information of the X-ray as the measuring condition of the measurement object S. Therefore, the estimated data D 2 can be generated without requiring a large load of processing.
  • the attenuation coefficient ⁇ (X,Y,Z, ⁇ ) of the X-ray is expressed using the attenuation coefficient ⁇ 0 (X,Y,Z, ⁇ ), which is a first component that depends on the photon energy ⁇ , and the error term ⁇ (X,Y,Z), which is a second component.
  • the image generation unit 53 calculates the value of ⁇ (X,Y,Z) using the detection data D 1 and the estimated data D 2 , thereby extracting differential data D 3 .
  • the image generation unit 53 performs back projection of the differential data D 3 to generate a back projection image Im 1 when all the values of ⁇ (X,Y,Z) that are calculated with respect to a certain irradiation direction of the X-ray do not exceed the first predetermined value. Therefore, a back projection image Im 1 can be generated so that a small cavity or the like included in the measurement object S can be observed.
  • the image generation unit 53 when a plurality of ⁇ (X,Y,Z) values calculated for a certain irradiation direction of an X-ray exceeds the first predetermined value, corrects the estimated data D 2 based on the ⁇ (X,Y,Z) values to generate new estimated data D 2 . Therefore, the defect such as a small cavity or the like included in the measurement object S can be correctly evaluated in condition that the error in shape between the measurement object S and the estimated structure S 2 is sufficiently small.
  • the image generation unit 53 when it generates new estimated data D 2 , extracts new differential data D 3 using the detection data D 1 and the new estimated data D 2 .
  • the image generation unit 53 when the new differential data D 3 does not exceed the first predetermined value, performs back projection of the new differential data D 3 to generate a back projection image Im 1 . Thereafter, the image generation unit 53 combines the back projection image Im 1 with the image Im 2 relating to the estimated structure S 2 to generate an image Im 3 that relates to the inside structure of the measurement object S.
  • the image of the estimated structure S 2 whose shape is substantially identical to that of the measurement object S, is combined with the back projection image Im 1 that shows a defect or the like, the shape that may be deemed to be identical to that of the measurement object S that is actually measured can be reproduced on an image of which occurrence of artifacts due to beam hardening is reduced. That is, a defect or the like in the inside of the actually measured measurement object S can be understood with ease.
  • the image generation unit 53 generates estimated data D 2 based on the concept that is explained in the first embodiment using the following expressions (7) and (7)′ and extracts differential data D 3 from the detection data D 1 that is based on the output from the detector 4 .
  • the image generation unit 53 generates estimated data D 2 by calculating Ai using expression (7)′ and calculates ⁇ (X,Y,Z) by substituting the calculated Ai into expression (7), thereby extracting the differential data D 3 .
  • the expressions (7) and (7)′ are derived as follows.
  • the attenuation coefficient ⁇ (X,Y,Z, ⁇ ) of the X-ray having a photon energy ⁇ at a position (X,Y,Z) where the X-rays pass through the measurement object S is assumed to be represented by expression (8) as follows.
  • ⁇ ( X,Y,Z , ⁇ ) ⁇ 1 ( X,Y,Z ) f PE ( ⁇ )+ ⁇ 2 ( X,Y,Z ) f KN ( ⁇ )+ ⁇ ( X,Y,Z ) (8)
  • ⁇ 1 (X,Y,Z) indicates a spatial distribution of reduction of the X-ray by photoelectric absorption in the measurement object S
  • ⁇ 2 (X,Y,Z) indicates a spatial distribution of reduction of the X-ray due to the Compton effect in the measurement object S.
  • ⁇ 1 (X,Y,Z) and ⁇ 2 (X,Y,Z) are known values that depend on the material of the measurement object S and are stored in a memory (not shown) in advance.
  • f PE ( ⁇ ) and f KN ( ⁇ ) which are functions that depend on the photon energy ⁇ , are represented by the following expressions (9) to (11), respectively.
  • ⁇ f PE ⁇ ( ⁇ ) 1 ⁇ 1 ( 9 )
  • 510.975 keV in expression (11) corresponds to static energy of an electron and ⁇ 0 represents photon energy of the X-ray that is normalized with the static energy of the electron.
  • ⁇ 1 (X,Y,Z)f PE ( ⁇ )+ ⁇ 2 (X,Y,Z)f KN ( ⁇ ) in expression (8) is information relating to a attenuation of the X-ray when the X-ray having a photon energy ⁇ passes through the inside of the estimated structure S 2 , which is estimated as the measurement object S in an ideal state based on the material information of the measurement object S.
  • ⁇ (X,Y,Z) is a term indicating an error between the attenuation coefficient ⁇ (X,Y,Z, ⁇ ) of the X-ray having a photon energy ⁇ at a position (X,Y,Z) when the X-ray passes through the inside of the measurement object S and the information relating to the attenuation described above.
  • the attenuation coefficient ⁇ (X,Y,Z, ⁇ ) the functions f PE ( ⁇ ) and f KN ( ⁇ ) are elements that depend on the photon energy ⁇ .
  • ⁇ (X,Y,Z) is an element that is assumed not to depend on the photon energy ⁇ .
  • the attenuation coefficient of the X-ray is represented by the following expression (12) using expression (8).
  • ⁇ ( X,Y,Z , ⁇ ) ⁇ l ⁇ 1 ( X,Y,Z ) f PE ( ⁇ )+ ⁇ 2 ( X,Y,Z ) f KN ( ⁇ ) ⁇ l + ⁇ ( X,Y,Z ) ⁇ l (12)
  • the expression (12) is changed using the spectrum Ii ( ⁇ ) of the X-ray that enters the pixels of the detector 4 in the same manner as the case in the first embodiment in which using expression (4) and then derives expression (5), the intensity yi of the X-ray that enters the detector 4 according to the present embodiment is substituted to obtain expression (13).
  • i indicates the position of one pixel among a plurality of pixels included in the detector 4 .
  • ri represents a noise component such as dark current contained in the output from the i-th pixel arranged in the detector 4 .
  • expression (7)′ is derived.
  • ⁇ 1 (X,Y,Z) and ⁇ 2 (X,Y,Z) are values that depend on the design information of the measurement object S.
  • the intensity Ii ( ⁇ ) depends on the intensity and spectrum of the X-ray, as the measuring conditions of the measuring device, which is irradiated to the measurement object S. Therefore, Ai represented by expression (7)′ corresponds to the calculated detection intensity assuming to be detected by the i-th detection pixel of the detector 4 based on the design information and the measuring conditions.
  • Ai corresponds to the transmission intensity of the X-ray that is assumed to being passed through the estimated structure S 2 when the X-ray assumed to be irradiated to the estimated structure S 2 according to the measuring conditions. That is, it corresponds to the estimated data D 2 . Therefore, the image generation unit 53 generates the estimated data D 2 by applying the design information and the measuring conditions to expression (7)′.
  • ⁇ (X,Y,Z) represents an error between the sum of the attenuation coefficient due to the photoelectric effect and the attenuation coefficient due to the Compton effect and the attenuation coefficient ⁇ (X,Y,Z, ⁇ ), and ⁇ (X,Y,Z) does not depend on the photon energy ⁇ , so that the error term ⁇ (X,Y,Z) can be calculated.
  • the image generation unit 53 performs the same processing as that in the first embodiment. That is, the image generation unit 53 evaluates the magnitude relationship between the calculated ⁇ (X,Y,Z) value and the first predetermined value and depending on the result of the evaluation, performs a generation of a back projection image and displaying a three-dimensional image, or performs a generation of a back projection image and displaying a three-dimensional image after correcting the estimated data D 2 .
  • the image generation unit 53 performs series of processing shown in the flowchart of FIG. 5 in the first embodiment to generate a three-dimensional image.
  • the image generation unit 53 calculates Ai using expression (7)′ and generates the estimated data D 2 .
  • the image generation unit 53 calculates ⁇ (X,Y,Z) using expression (7).
  • the image generation unit 53 uses ⁇ 1 (X,Y,Z), which is information relating to the photoelectric absorption, and ⁇ 2 (X,Y,Z), which is information relating to the Compton effect for expressing the attenuation coefficient ⁇ (X,Y,Z, ⁇ ).
  • the attenuation coefficient ⁇ (X,Y,Z, ⁇ ) may be expressed using either one of ⁇ 1 (X,Y,Z) or ⁇ 2 (X,Y,Z).
  • the image generation unit 53 uses the information relating to the photoelectric absorption, ⁇ 1 (X,Y,Z), and when the photon energy ⁇ is relatively large, the image generation unit 53 uses the information relating to the Compton effect, ⁇ 2 (X,Y,Z).
  • the image generation unit 53 generates the estimated data D 2 based on the spectrum information of the X-ray that enters the measurement object S and at least one of the information of, ⁇ 1(X,Y,Z), which is information relating to the photoelectric absorption of the X-ray that is passed through the measurement object S, and ⁇ 2(X,Y,Z), which is information relating to the Compton effect of the X-ray that is passed through the measurement object S. That is, the image generation unit 53 calculates Ai using the expression (7)′. Therefore, the estimated data D 2 can be generated without requiring a large load of processing.
  • the image generation unit 53 calculates ⁇ (X,Y,Z), that is, extracts differential data D 3 using expression (1) in the same manner as that in the first embodiment.
  • the image generation unit 53 determines the magnitude relationship between the calculated ⁇ (X,Y,Z) and the second predetermined value.
  • the second predetermined value may be, for instance, an acceptable tolerance value of the measurement object S.
  • ⁇ (X,Y,Z) is equal to or smaller than the second predetermined value, it is determined that the shape of the measurement object S is within an acceptable tolerance range with respect to the shape of the estimated structure S 2 . That is, in this case the image generation unit 53 determines that the measurement object S does not have a large error in shape with respect to the design dimensions and determines that this measurement object S is a conforming product.
  • the image generation unit 53 performs generation of a back projection image Im 1 , combination of the back projection image Im 1 with the image Im 2 relating to the estimated structure S 2 , and generation of the image Im 3 in the same manner as that in the first embodiment.
  • the image generation unit 53 generates a three-dimensional image using the image Im 3 and displays it on the monitor 6 . That is, the image generation unit 53 deems that the estimated structure S 2 obtained by repeating corrections until the back projection image Im 1 of the differential data D 3 becomes 0 matches the shape of the actual measurement object S and the image generation unit 53 generates a three-dimensional image using the image Im 3 .
  • the image Im 3 may be displayed in comparison with the design image, or only a differential image between the image Im 3 and the design image may be displayed.
  • the image generation unit 53 determines that the shape of the measurement object S exceeds an acceptable tolerance range with respect to the shape of the estimated structure S 2 . That is, it determines that the measurement object S has a significantly large error in shape with respect to the design dimensions. In this case the image generation unit 53 determines that the measurement object S is a defective product and does not perform subsequent processing. Note that the image generation unit 53 may not display the three-dimensional image of the measurement object S on the display monitor 6 and may, for instance, display a warning notifying that the measurement object S is a defective product on the monitor 6 .
  • the image generation unit 53 may generate a three-dimensional image of the measurement object S and displays it on the display monitor 6 . In this case, the image generation unit 53 may display a warning notifying that the measurement object S is a defective product in superposition to the three-dimensional image of the measurement object S displayed on the display monitor 6 . Note that although in the above explanation, the image generation unit 53 calculates ⁇ (X,Y,Z) in the same manner as that in the first embodiment, it may be configured to calculate ⁇ (X,Y,Z) in the same manner as that in the second embodiment.
  • Step S 11 acquisition of design information and measuring conditions
  • Step S 14 extraction of differential data
  • Step S 4 extraction of differential data
  • Step S 15 it is determined whether all the position-dependent ⁇ (X,Y,Z) values that are calculated for respective pixels are equal to or smaller than the second predetermined value. When all the ⁇ (X,Y,Z) values are equal to or smaller than the second predetermined value, Step S 15 is determined to be affirmative and the procedure proceeds to Step S 17 . When at least one of the ⁇ (X,Y,Z) values exceeds the second predetermined value, Step S 15 is determined to be negative and the procedure proceeds to Step S 16 . In Step S 16 , the measurement object S is determined to be a defective product and the processing is terminated. Note that in Step S 16 , a message that the measurement object S is a defective product or the like may be displayed on the display monitor 6 .
  • Step S 17 determines whether the shape of the corrected estimated structure is the shape of the actual measurement object (structure) is similar to each processing of from Step S 5 (determination of magnitude relationships between each of the ⁇ (X,Y,Z) values and the first predetermined value) to Step S 7 (deeming the shape of the corrected estimated structure as the shape of the actual measurement object (structure) in FIG. 5 .
  • Step S 13 determines whether the processing in Step S 13 is performed again via Step S 18 .
  • the image generation unit 53 when ⁇ (X,Y,Z) values that correspond to the differential data D 3 do not exceed the second predetermined value, determines the measurement object S to be a conforming product. Therefore, it becomes possible to use ⁇ (X,Y,Z) values in processing other than the generation of images, so that convenience can be increased.
  • the image generation unit 53 when it determines the measurement object S as a conforming product based on the differential data D 3 , performs back projection of the differential data D 3 to generate a back projection image Im 1 and combines it with the image Im 2 relating to the estimated structure S 2 to generate an image Im 3 of the inside structure of the measurement object S. Therefore, in case where the measurement object S does not have a large error in shape, the Im 3 relating to the inside structure of the measurement object S can be generated. This enables prevention of any possible increase in load to be caused by generation of images of low necessity for the measurement object S that has been determined to be a defective product can be prevented.
  • the structure manufacturing system according to the present embodiment manufactures molding products, for instance, door parts, engine parts, or gear parts for automotive use, and electronic components provided with circuit boards, and the like.
  • FIG. 7 is a block diagram showing an example of the configuration of a structure manufacturing system 600 according to the present embodiment.
  • the structure manufacturing system 600 includes the X-ray device 100 explained in any one of the first to third embodiments and variation examples thereof, a design device 610 , a molding device 620 , a control system 630 , and a repairing device 640 .
  • the design device 610 which is a device that is used by the user for generating design information relating to the shape of a structure, performs design processing to generate design information and store it.
  • the design information is information that indicates respective positions of the structure on the coordinate.
  • the design information is outputted to the molding device 620 and the control system 630 described later.
  • the molding device 620 performs molding processing to manufacture the structure based on the design information generated by the design device 610 .
  • the molding device 620 that performs at least one of a lamination processing represented by 3D printing processing, a casting processing, a forging processing, and a cutting processing is included as an aspect of the present invention.
  • the X-ray device 100 performs measurement processing in which the shape of the structure molded by the molding device 620 is measured.
  • the X-ray device 100 outputs information that indicates coordinate information (hereafter, referred to as “shape information”) of the structure, which information is a result of measurement of the structure to the control system 630 .
  • the control system 630 includes a coordinate storage unit 631 and an inspection unit 632 .
  • the coordinate storage unit 631 stores the design information generated by the design device 610 .
  • the inspection unit 632 determines whether the structure molded by the molding device 620 has been manufactured according to the design information generated by the design device 610 . In other words, the inspection unit 632 determines whether the manufactured structure is a conforming product. In this case, the inspection unit 632 reads out the design information stored at the coordinate storage unit 631 and performs inspection processing, in which the design information is compared with the shape information that is inputted from the image measuring device 100 .
  • the inspection unit 632 performs, as inspection processing, comparison of, for instance, coordinates indicated by the design information with corresponding coordinates indicated by the shape information and determines, when the coordinates of the design information matches the coordinates of the shape information as a result of the inspection processing, that the structure is a conforming product that has been manufactured according to the design information.
  • the inspection unit 632 determines whether a difference in coordinates is within a predetermined range. If the difference is within the predetermined range, the inspection unit 632 determines that the structure is a repairable defective product.
  • the inspection unit 632 When the structure is determined to be a repairable defective product, the inspection unit 632 outputs repair information that indicates a defective region and an amount of correction to a repairing device 640 .
  • the defective region indicates the region in which the coordinates of the shape information does not match the coordinates of the design information, and the amount of correction is a difference between the coordinates of the design information and the coordinates of the shape information in the defective region.
  • the repairing device 640 performs repair processing in which the defective region of the structure is reprocessed based on the inputted repair information.
  • the repairing device 640 in repair processing performs again processing similar to the molding processing performed by the molding device 620 .
  • Step S 111 the design device 610 , which is used by the user for designing a structure, performs design processing to generate design information relating to the shape of the structure and stores it. Then the procedure proceeds to Step S 112 .
  • the design information is not limited to the one that is generated by the design device 610 .
  • the design device 610 may, when there is the design information that has already been prepared, acquire such design information by inputting it therein.
  • Step S 112 the molding device 620 performs molding processing to manufacture a structure based on the design information and the procedure proceeds to Step S 113 .
  • Step S 113 the image measuring device 100 performs measurement processing to measure the shape of the structure and outputs shape information, and then the procedure proceeds to Step S 114 .
  • Step S 114 the inspection unit 632 performs inspection processing in which the design information generated by the design device 610 is compared with the shape information acquired by measurement and outputted by the image measuring device 100 , and then the procedure proceeds to Step S 115 .
  • Step S 115 the inspection unit 632 determines whether the structure manufactured by the molding device 620 is a conforming product based on the result of the inspection processing. When the structure is a conforming product, that is, when the coordinates of the design information matches the coordinates of the shape information, Step S 115 is determined to be affirmative and the processing is terminated.
  • Step S 115 is determined to be negative, and then the procedure proceeds to Step S 116 .
  • Step S 116 the inspection unit 632 determines whether the defective region of the structure is repairable.
  • Step S 116 is determined to be negative, and then the processing is terminated.
  • Step S 116 is determined to be affirmative, and then the procedure proceeds to Step S 117 .
  • the inspection unit 632 outputs repair information to the repairing device 640 .
  • Step S 117 the repairing device 640 performs repair processing to the structure based on the inputted repair information, and then the procedure returns back to Step S 113 . Note that as mentioned above, the repairing device 640 performs again processing similar to the molding processing performed by the molding device 620 in the repair processing.
  • the X-ray device 100 in the structure manufacturing system 600 performs measurement processing, in which the shape information of the structure generated by the molding device 620 by the design processing performed by the design device 610 is acquired and the inspection unit 632 in the control system 630 performs inspection processing, in which the shape information acquired by the measurement processing is compared with the design information generated in the design processing. Therefore, it is possible to inspect the defect of the structure or acquire the information of the inside of the structure by a non-destructive inspection and determine whether the structure is a conforming product that has been manufactured according to the design information. This contributes quality control of the structure.
  • the repairing device 640 is configured to perform repair processing in which the molding processing is performed again to the structure depending on the result of comparison in the inspection processing. Therefore, when the defective region of the structure is repairable, processing similar to the molding processing can be performed to the structure again. This contributes to manufacture of a high quality structure that is approximate to the design information.
  • data showing a difference between the measurement object and an estimated structure obtained based on the shape information of the measurement object can be extracted.
  • variation examples as described below are within the scope of the present invention and one or more variation examples may be combined with the above-mentioned embodiments.
  • the image generation unit 53 has been explained as estimating the shape of the estimated structure S 2 using the shape information contained in the design information.
  • the X-ray device 100 may be one that generates the estimated data D 2 based on the detection data D 1 of the measurement object S.
  • the image generation unit 53 performs back projection of the detection data D 1 of the measurement object S using, for instance, filtered back projection (FBP) to generate an image. Since this image is expressed by a contrasting density in correspondence to the intensity of the X-ray that has been passed through the measurement object S, the image generation unit 53 estimates a substance that constitutes the measurement object S, that is, a material based on this contrasting density.
  • FBP filtered back projection
  • the image generation unit 53 is configured to have data in which the density of the image is in relation to the material and estimate the material based on the density of the generated image, that is, the intensity of the X-ray.
  • the image generation unit 53 may use the estimated material as the material information and calculate Bi according to expression (1)′ to generate the estimated data D 2 .
  • the material information of the measurement object S is acquired from projection image of the measurement object S in Step S 1 of FIG. 5 or Step S 11 of FIG. 6 . Therefore, when no material information is obtained as the design information of the measurement object S, the estimated data D 2 can be generated.
  • the image generation unit 53 may be configured to extract shape information of the measurement object S by performing edge detection or the like using the back projection image of the measurement object S.
  • the shape of the measurement object S can be obtained with high precision.
  • the image generation unit 53 may be configured to extract the profile shape of the measurement object S using an image of the measurement object S captured by an image-capturing device such as a camera and to use the extracted external contour as the shape of the estimated structure S 2 .
  • the X-ray device 100 includes an image-capturing device (not shown) provided with image sensors that are constituted by CMOS, CCD or the like.
  • the image-capturing device which is provided on a part of the ceiling of the housing 1 (on the internal wall surface on the positive side of Y-axis), captures an image of the external appearance of the measurement object S mounted on the mount stage 30 along the Y-axis direction that is substantially orthogonal to the direction of projection of the X-ray (Z-axis) and the image-capturing device outputs the generated image signal to the control unit 5 .
  • the image generation unit 53 performs known (edge) detection processing and the like to the inputted image signal to extract a contour of the measurement object S on the image signal.
  • the image-capturing device is preferably one that is capable of capturing images over a wide range so that the measurement object S in whole can be captured regardless of the position of the measurement object S that may be varied by the Y-axis movement unit 33 and the X-axis movement unit 34 .
  • the image-capturing device may be provided so that it can be moved in synchronization with the movements of the Y-axis movement unit 33 and the X-axis movement unit 34 .
  • the shape information is not limited to one that is obtained by the image-capturing device.
  • the shape information may include information obtained by measuring the measurement object S using a profile projector that projects the optical image of the measurement object S on a screen, or a shape measuring device, which is a contact type three-dimensional measuring device using a touch probe, or which is a non-contact type three-dimensional measuring device, such as a scanning laser probe type or an optical type.
  • the mount stage 30 on which the measurement object S is mounted is not limited to one that is moved along X-axis direction, Y-axis direction, and Z-axis direction by the X-axis movement unit 34 , the Y-axis movement unit 33 , and the Z-axis movement unit 35 , respectively.
  • the mount stage 30 may be one that does not move along the X-axis, the e Y-axis, and the Z-axis direction but instead causes the X-ray source 2 and the detector 4 move along the X-axis direction, the Y-axis direction, and the Z-axis direction, respectively, relative to the measurement object S.
  • the X-ray device 100 may be structured so that the mount stage 30 does not rotate but the X-ray source 2 and the detector 4 rotate about the rotation axis Yr instead of the mount stage 30 that rotates about the rotation axis Yr.
  • the differential data D 3 may be generated using only a part of the transmission images.
  • the function of an interface through which the detection data D 1 , the design information of the measurement object S, and the measuring conditions are inputted and the function of the image generation unit 53 may be implemented by a computer.
  • the implementation may be achieved by storing a program for implementing the function of the image generation processing into a computer-readable recording medium and causing a computer system to read in the above-mentioned program relating to image generation that is stored in the recording medium and execute it.
  • the “computer system” as used herein includes hardware such as an OS (Operating System) and peripheral devices.
  • the “computer-readable recording medium” as used herein includes portable recording media such as a flexible disk, a magneto-optical disk, an optical disk, and a memory card and so on, and a storage device such as a hard disk incorporated in the computer system. Furthermore, the “computer-readable recording medium” as used herein may include a thing that can hold a program for a short time and dynamically such as a communication wire in case a program is transmitted through a network such as an internet or a telephone line, and a thing that holds the program for a certain period of time such as a volatile memory in the inside of the computer system serving as a server or a client in such a case.
  • the program may be one that implements a part of the above-mentioned function or that implements the above-mentioned function in combination with the program that has already been stored in the computer system.
  • FIG. 9 is a diagram illustrating this situation.
  • a personal computer 950 receives a program via a CD-ROM 953 .
  • the personal computer 950 has a connection function to be connected with a communication line 951 .
  • a computer 952 is a server-computer that stores above described program into a recording medium such as a hard disk.
  • a communication line 951 is a communication line such as the internet, a personal computer communication, or the like, or a dedicated communication line.
  • a computer 952 reads out the program from the hard disk and transmits the read out program to the personal computer 950 through the communication line 951 . That is, it carries the program as data signal with a carrier wave to transmit it through the communication line 951 .
  • the program can be provided as a computer-readable computer program product in one of various forms such as a recording medium having it stored therein or a carrier wave.
  • the present invention is not limited to the above-mentioned embodiments and/or variation examples and other forms conceivable within the technical idea of the present invention are included within the scope of the present invention.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
US15/549,010 2015-02-09 2015-02-09 Image reconstruction method for x-ray measuring device, structure manufacturing method, image reconstruction program for x-ray measuring device, and x-ray measuring device Abandoned US20180024079A1 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/053518 WO2016129035A1 (fr) 2015-02-09 2015-02-09 Procédé de reconstruction d'image pour dispositif de mesure par rayons x, procédé de fabrication de la structure, programme de reconstruction d'image pour dispositif de mesure par rayons x et dispositif de mesure par rayons x

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/053518 A-371-Of-International WO2016129035A1 (fr) 2015-02-09 2015-02-09 Procédé de reconstruction d'image pour dispositif de mesure par rayons x, procédé de fabrication de la structure, programme de reconstruction d'image pour dispositif de mesure par rayons x et dispositif de mesure par rayons x

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US16/879,389 Continuation US11860111B2 (en) 2015-02-09 2020-05-20 Image reconstruction method for X-ray measuring device, structure manufacturing method, image reconstruction program for X-ray measuring device, and X-ray measuring device

Publications (1)

Publication Number Publication Date
US20180024079A1 true US20180024079A1 (en) 2018-01-25

Family

ID=56615118

Family Applications (2)

Application Number Title Priority Date Filing Date
US15/549,010 Abandoned US20180024079A1 (en) 2015-02-09 2015-02-09 Image reconstruction method for x-ray measuring device, structure manufacturing method, image reconstruction program for x-ray measuring device, and x-ray measuring device
US16/879,389 Active 2036-07-19 US11860111B2 (en) 2015-02-09 2020-05-20 Image reconstruction method for X-ray measuring device, structure manufacturing method, image reconstruction program for X-ray measuring device, and X-ray measuring device

Family Applications After (1)

Application Number Title Priority Date Filing Date
US16/879,389 Active 2036-07-19 US11860111B2 (en) 2015-02-09 2020-05-20 Image reconstruction method for X-ray measuring device, structure manufacturing method, image reconstruction program for X-ray measuring device, and X-ray measuring device

Country Status (5)

Country Link
US (2) US20180024079A1 (fr)
EP (1) EP3258251B1 (fr)
JP (1) JP6658552B2 (fr)
CN (1) CN107209130B (fr)
WO (1) WO2016129035A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180365847A1 (en) * 2017-06-20 2018-12-20 Mitutoyo Corporation Three-dimensional geometry measurement apparatus and three-dimensional geometry measurement method

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6729801B2 (ja) * 2017-04-26 2020-07-22 株式会社ニコン 検査装置、検査方法および検査対象物の製造方法
US11446861B2 (en) * 2017-08-30 2022-09-20 Ecole Polytechnique Federale De Lausanne (Epfl) Methods and apparatus for three-dimensional fabrication by tomographic back projections
EP4251981A1 (fr) * 2020-11-25 2023-10-04 Shenzhen Xpectvision Technology Co., Ltd. Appareil et procédé pour l'étude d'analytes biologiques par fluorescence de rayons x

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4852131A (en) * 1988-05-13 1989-07-25 Advanced Research & Applications Corporation Computed tomography inspection of electronic devices
US5436829A (en) 1992-11-05 1995-07-25 General Electric Company Method of achieving reduced dose X-ray fluoroscopy by employing transform-based estimation of Poisson noise
JP2002248344A (ja) * 2001-02-26 2002-09-03 Nikon Corp 極端紫外光発生装置並びにそれを用いた露光装置及び半導体製造方法
US6823044B2 (en) * 2001-11-21 2004-11-23 Agilent Technologies, Inc. System for collecting multiple x-ray image exposures of a sample using a sparse configuration
JP3909048B2 (ja) * 2003-09-05 2007-04-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置およびｘ線管
US7023951B2 (en) * 2003-12-09 2006-04-04 General Electric Company Method and apparatus for reduction of artifacts in computed tomography images
JP4709162B2 (ja) * 2003-12-16 2011-06-22 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ ヒール効果によって引き起こされるアーチファクトの補正
WO2006070316A1 (fr) * 2004-12-29 2006-07-06 Philips Intellectual Property & Standards Gmbh Appareil et procede de correction d'artefacts de projections de rayons x
US7310403B2 (en) * 2006-03-20 2007-12-18 Moshe Ein-Gal Target reconstruction
EP2126931B1 (fr) 2007-01-26 2011-09-21 Koninklijke Philips Electronics N.V. Filtre de compensation de l'effet talon qui conserve le spectre, réalisé dans le même matériau qu'une plaque anodique
CN101266216A (zh) * 2007-03-14 2008-09-17 清华大学 标定双能ct系统的方法和图像重建方法
DE112008002137T5 (de) * 2007-08-07 2010-09-02 Shimadzu Corp. Fluoroskop
JP5177633B2 (ja) * 2007-11-16 2013-04-03 株式会社Ｉｈｉ 材質識別検査装置および方法
CN101726503B (zh) * 2008-10-17 2012-08-29 清华大学 用于x射线相衬层析成像的系统和方法
US8180138B2 (en) * 2009-03-23 2012-05-15 Morpho Detection, Inc. Method and system for inspection of containers
US20120093287A1 (en) * 2010-10-18 2012-04-19 Satpal Singh Two unit portable x-ray scanner
JP5834584B2 (ja) 2011-07-25 2015-12-24 ソニー株式会社 情報処理装置、情報処理方法、プログラム及び蛍光スペクトルの強度補正方法
JP2013217775A (ja) * 2012-04-09 2013-10-24 Nikon Corp Ｘ線装置、方法、及び構造物の製造方法
JP2013217797A (ja) * 2012-04-10 2013-10-24 Nikon Corp 装置、判定方法、及び構造物の製造方法
JP2014016239A (ja) * 2012-07-09 2014-01-30 Tokyo Electron Ltd Ｘ線検査方法及びｘ線検査装置
WO2014192889A1 (fr) * 2013-05-29 2014-12-04 地方独立行政法人東京都立産業技術研究センター Dispositif et procédé de reconstruction d'image à différentes énergies de rayons x, et dispositif et procédé de mesure tridimensionnelle aux rayons x
CN104240270B (zh) * 2013-06-14 2017-12-05 同方威视技术股份有限公司 Ct成像方法和系统
CN103489223A (zh) * 2013-09-29 2014-01-01 华南理工大学 一种面向3d封装ic的3d断层扫描成像方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180365847A1 (en) * 2017-06-20 2018-12-20 Mitutoyo Corporation Three-dimensional geometry measurement apparatus and three-dimensional geometry measurement method
US10529082B2 (en) * 2017-06-20 2020-01-07 Mitutoyo Corporation Three-dimensional geometry measurement apparatus and three-dimensional geometry measurement method

Also Published As

Publication number Publication date
JP6658552B2 (ja) 2020-03-04
CN107209130B (zh) 2021-07-09
EP3258251A1 (fr) 2017-12-20
EP3258251B1 (fr) 2021-03-24
WO2016129035A1 (fr) 2016-08-18
CN107209130A (zh) 2017-09-26
US20200284735A1 (en) 2020-09-10
JPWO2016129035A1 (ja) 2017-11-24
EP3258251A4 (fr) 2018-10-17
US11860111B2 (en) 2024-01-02

Similar Documents

Publication Publication Date Title
US11860111B2 (en) Image reconstruction method for X-ray measuring device, structure manufacturing method, image reconstruction program for X-ray measuring device, and X-ray measuring device
JP6631624B2 (ja) Ｘ線検査装置、ｘ線検査方法および構造物の製造方法
US7801264B2 (en) Method for calibrating a dual -spectral computed tomography (CT) system
JP6280544B2 (ja) Ｘ線エネルギー別画像再構成装置及び方法並びにｘ線三次元測定装置及び方法
US9031195B2 (en) Imaging detector and methods for image detection
US10809209B2 (en) Measurement processing device, x-ray inspection device, measurement processing method, measurement processing program, and structure manufacturing method
JP6455516B2 (ja) Ｘ線装置および構造物の製造方法
US20070114428A1 (en) Radiation Image Capturing Condition Correction Apparatus and Radiation Image Capturing Condition Correction Method
JP2015075336A (ja) 再構成画像生成装置、形状測定装置、構造物製造システム、再構成画像生成方法及び再構成画像生成プログラム
Rogers et al. Reduction of wobble artefacts in images from mobile transmission x-ray vehicle scanners
JP5884351B2 (ja) Ｘ線検査装置、ｘ線検査装置の制御方法、ｘ線検査装置を制御するためのプログラム、および、当該プログラムを格納したコンピュータ読み取り可能な記録媒体
WO2019163960A1 (fr) Procédé de mesure à rayons x, dispositif à rayons x et procédé de fabrication de structure
US20240057955A1 (en) System and method for incorporating lidar-based techniques with a computed tomography system
JP2013215472A (ja) 放射線断層画像生成装置、放射線断層撮影装置および放射線断層画像生成プログラム
US11937970B2 (en) System and method for calibrating a camera feature detection system of an x-ray system
JP5492634B2 (ja) 放射線断層撮影方法および放射線断層撮影装置
KR20220127726A (ko) Ct 장치
KR20210136108A (ko) 바디 변형의 시분할 계산을 위한 방법 및 컴퓨터 프로그램
JP2013205267A (ja) Ｘ線断層撮影方法およびｘ線断層撮影装置
JP2021028575A (ja) Ｘ線装置および構造物の製造方法
JP2009109416A (ja) 組立品の検査方法及び検査システム
WO2016056107A1 (fr) Générateur de données de projection, dispositif de mesure, et procédé de fabrication de structure

Legal Events

Date Code Title Description
AS Assignment

Owner name: NIKON CORPORATION, JAPAN

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:YAMADA, ATSUSHI;REEL/FRAME:043315/0613

Effective date: 20170711

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO PAY ISSUE FEE