[go: up one dir, main page]

WO2010050483A1 - Ｘ線撮像装置およびｘ線撮像方法 - Google Patents

Ｘ線撮像装置およびｘ線撮像方法 Download PDF

Info

Publication number
WO2010050483A1
WO2010050483A1 PCT/JP2009/068434 JP2009068434W WO2010050483A1 WO 2010050483 A1 WO2010050483 A1 WO 2010050483A1 JP 2009068434 W JP2009068434 W JP 2009068434W WO 2010050483 A1 WO2010050483 A1 WO 2010050483A1
Authority
WO
WIPO (PCT)
Prior art keywords
phase
moire
intensity distribution
grating
ray imaging
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.)
Ceased
Application number
PCT/JP2009/068434
Other languages
English (en)
French (fr)
Inventor
大内　千種
英之助 伊藤
長井　健太郎
田　透
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.)
Canon Inc
Original Assignee
Canon Inc
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 to DE112009002606.0T priority Critical patent/DE112009002606B4/de
Priority to EP09823593.0A priority patent/EP2343537B1/en
Priority to KR1020117011759A priority patent/KR101258927B1/ko
Priority to JP2010535804A priority patent/JP5174180B2/ja
Priority to CN2009801428377A priority patent/CN102197303A/zh
Application filed by Canon Inc filed Critical Canon Inc
Publication of WO2010050483A1 publication Critical patent/WO2010050483A1/ja
Priority to US12/842,937 priority patent/US8009797B2/en
Priority to US13/085,199 priority patent/US8559594B2/en
Anticipated expiration legal-status Critical
Priority to US13/190,770 priority patent/US8340243B2/en
Priority to US13/682,445 priority patent/US8537966B2/en
Ceased 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/041—Phase-contrast imaging, e.g. using grating interferometers
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • 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
    • 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/20—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 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/20075—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 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 interferences of X-rays, e.g. Borrmann effect
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K1/00—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating
    • G21K1/06—Arrangements for handling particles or ionising radiation, e.g. focusing or moderating using diffraction, refraction or reflection, e.g. monochromators
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/484—Diagnostic techniques involving phase contrast X-ray imaging
    • 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
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K2201/00—Arrangements for handling radiation or particles
    • G21K2201/06—Arrangements for handling radiation or particles using diffractive, refractive or reflecting elements
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
    • G21K2201/00—Arrangements for handling radiation or particles
    • G21K2201/06—Arrangements for handling radiation or particles using diffractive, refractive or reflecting elements
    • G21K2201/067—Construction details

Definitions

  • the present invention relates to an X-ray imaging apparatus and an X-ray imaging method.
  • X-rays Since X-rays have high substance permeability and high spatial resolution imaging is possible, they are used for non-destructive inspection of objects for industrial use and radiography for medical use. In these methods, a contrast image is formed by utilizing a difference in absorption at the time of X-ray transmission due to a constituent element or density difference in an object or a living body, which is called an X-ray absorption contrast method.
  • an X-ray absorption contrast method a difference in absorption at the time of X-ray transmission due to a constituent element or density difference in an object or a living body.
  • X-ray absorption contrast method since X-ray absorption is very small in light elements, it is difficult to image a biological soft tissue or soft material composed of carbon, hydrogen, oxygen, etc., which are constituent elements of a living body, by the X-ray absorption contrast method.
  • Patent Document 1 X-ray phase contrast method using Talbot interference as a method that can use a normal X-ray tube
  • the method using Talbot interference described in Patent Document 1 is an X-ray tube that generates X-rays, a phase grating that modulates the phase of the X-rays to generate an interference intensity distribution, and the interference intensity distribution is an intensity of moire. It comprises an absorption grating that converts to a distribution and an X-ray detector that detects the interference intensity distribution.
  • imaging is performed by scanning the absorption grating in the grating period direction. By this scanning, the detected moire moves, and when the scanning amount reaches one period of the absorption grating, the moire image returns to the original state.
  • a differential phase image is acquired by performing arithmetic processing using at least three or more pieces of imaging data being scanned.
  • Patent Document 1 is a method of obtaining a differential phase image by performing imaging of at least three sheets and calculating a phase image from the differential phase image.
  • Patent Document 1 since the method described in Patent Document 1 requires three or more images, there is a problem that the image quality deteriorates when the subject moves during the imaging.
  • the longer imaging time causes an increase in the X-ray dose to the subject, which is not preferable in consideration of medical use.
  • an object of the present invention is to provide an X-ray imaging apparatus and an X-ray imaging method capable of acquiring a differential phase image or a phase image of a subject from at least one imaging.
  • An X-ray imaging apparatus includes an X-ray source, a phase grating that forms an interference intensity distribution due to the Talbot effect by transmitting X-rays from the X-ray source, and interference formed by the phase grating.
  • An absorption grating that generates moiré by shielding a part of the intensity distribution, a detector that detects the intensity distribution of the moiré generated by the absorption grating, and the intensity distribution of the moiré detected by the detector
  • An arithmetic device that images and outputs information, and the arithmetic device performs a Fourier transform on the moire intensity distribution obtained by the detector to obtain a spatial frequency spectrum, and the Fourier transform step.
  • the phase that separates the spectrum corresponding to the carrier frequency from the spatial frequency spectrum obtained by, and acquires the differential phase image using the inverse Fourier transform Characterized in that it is configured to perform a recovery step.
  • an X-ray imaging apparatus and an X-ray imaging method capable of acquiring a differential phase image or phase image of a subject from at least one imaging.
  • FIG. 1 An example of the configuration of an X-ray imaging apparatus using Talbot interference is shown in FIG. A process until obtaining a phase image using an X-ray imaging apparatus will be specifically described.
  • X-ray source X-rays 111 generated from the X-ray source 110 pass through the subject 120.
  • X-ray 111 passes through the subject 120, a change in phase and absorption according to the composition, shape, and the like of the subject 120 occur.
  • the X-ray either a continuous X-ray or a characteristic X-ray may be used.
  • the wavelength is appropriately selected from about 0.1 to 5 mm.
  • a wavelength selection filter or a grating for the radiation source may be provided as appropriate.
  • phase grating When the X-ray 111 transmitted through the subject 120 passes through the phase grating 130, an interference intensity distribution 140 is formed by the Talbot effect.
  • the phase grating 130 is arranged on the upstream side or the downstream side of the subject 120.
  • the phase grating 130 includes a phase advance portion 131 and a phase delay portion 132 that are configured by periodically changing the thickness of the X-ray transmitting member.
  • the phase progression unit 131 and the phase delay unit 132 may be formed so as to have a phase difference with respect to the transmitted X-ray.
  • the configuration is such that the phase of the X-ray transmitted through the phase advancer 131 advances by ⁇ with respect to the phase of the X-ray transmitted through the phase delay unit 132.
  • the amount of change in thickness is determined by the wavelength and member of the X-ray used.
  • the phase grating 130 modulates the phase of the X-ray transmitted through the phase delay unit 132 by ⁇ or ⁇ / 2 with respect to the phase of the X-ray transmitted through the phase delay unit 132.
  • the former may be referred to as a ⁇ phase grating and the latter as a ⁇ / 2 phase grating.
  • the phase modulation amount may be periodic, and may be, for example, ⁇ / 3 modulated.
  • the phase grating 130 may have a one-dimensional linear shape, or may be configured with a two-dimensional checkered pattern as shown in FIG. Moreover, you may be comprised by the screen door-like pattern as shown to FIG. 2 (B).
  • d is a period
  • 201 is a two-dimensional phase grating
  • 210 is a phase progression unit
  • 220 is a phase delay unit.
  • phase advance part 210 or the phase delay part 220 is a square in FIG. 2 (A) and FIG. 2 (B)
  • an outer edge may deform
  • it can be used as a phase grating.
  • phase grating 130 When the phase grating 130 has a one-dimensional period, only the phase gradient information of the subject 120 in the one-dimensional direction can be acquired. However, when the phase grating 130 has a two-dimensional period, it is advantageous in that phase gradient information in a two-dimensional direction can be acquired.
  • the material constituting the phase grating 130 is preferably a substance that transmits X-rays.
  • silicon or the like can be used.
  • the interference intensity distribution formed after transmission through the phase grating 130 is the distance Z 0 between the X-ray source and the phase grating 130, the distance Z 1 from the absorption grating 150 satisfies the following formula (1). Appears most clearly.
  • the “interference intensity distribution” is a periodic intensity distribution reflecting the grating period of the phase grating 130.
  • ⁇ is the wavelength of the X-ray
  • d is the grating period of the phase grating 130.
  • N is a different value depending on the form of the phase grating, and is a real number that can be expressed as follows. Note that n is a natural number.
  • One-dimensional array of ⁇ phase grating: N n / 4 ⁇ 1 / 8
  • Two-dimensional checkered pattern ⁇ phase grating: N n / 4 ⁇ 1 / 8
  • Two-dimensional checkered pattern ⁇ / 2 phase grating: N n / 2 ⁇ 1 / 4
  • Absorption grid Since the period of the interference intensity distribution is generally smaller than the pixel size of the X-ray detector 170, the interference intensity distribution cannot be detected as it is. Therefore, by using the absorption grating 150, moire having a period larger than the pixel size of the X-ray detector 170 is generated, and the moire intensity distribution is detected by the X-ray detector 170.
  • Absorption grating 150 is preferably provided apart from the phase grating 130 by the distance Z 1 position.
  • the absorption grating 150 includes a transmission portion 151 and a light shielding portion 152 that are periodically arranged, and is arranged so as to shield a part of a bright portion of the interference intensity distribution 140 formed by the phase grating 130.
  • the transmission part 151 should just be comprised so that a part of X-rays can permeate
  • the material constituting the absorption grating 150 may be any material that absorbs X-rays well, and for example, gold can be used.
  • the period of the absorption grating 150 is the same as or slightly different from the interference intensity distribution.
  • the transmission portions 151 may be arranged one-dimensionally or two-dimensionally.
  • the transmission part 351 and the light shielding part 352 are two-dimensionally arranged as shown in FIG.
  • a screen door-shaped absorption grating 300 is used.
  • the checkered ⁇ / 2 phase grating shown in FIG. 2A is used, a checkered pattern in which the transmission part 351 and the light shielding part 352 are two-dimensionally arranged as shown in FIG. 3B.
  • the absorption grating 300 is used.
  • phase grating and the absorption grating are possible.
  • the X-ray detector 170 is an element that can detect information on the X-ray interference intensity distribution.
  • an FPD Full Panel Detector
  • FPD Full Panel Detector
  • the arithmetic device 180 has, for example, a CPU (Central Processing Unit).
  • the interference intensity distribution When the interference intensity distribution is formed, a large number of diffracted lights interfere and interfere with each other, and therefore include a fundamental frequency (hereinafter referred to as a carrier frequency) and a number of harmonic components thereof.
  • a carrier frequency a fundamental frequency
  • the moire since the moire has a shape in which the carrier frequency component of the interference intensity distribution is spatially enlarged, the moire can be expressed by equation (2) when a one-dimensional phase grating whose score line is orthogonal to the x axis is used. it can.
  • Equation (2) indicates that moire is represented by the sum of the first term of the background and the second term having periodicity.
  • a (x, y) represents the background
  • b (x, y) represents the amplitude of the carrier frequency component.
  • F 0 represents the carrier frequency of the interference fringes
  • ⁇ (x, y) represents the phase of the carrier frequency component.
  • the carrier frequency component is generated by the interference between the 0th-order diffracted light and the + 1st-order diffracted light, and the interference between the 0th-order diffracted light and the ⁇ 1st-order diffracted light. Further, when a checkered ⁇ phase grating is used as the phase grating 130, a carrier frequency component is generated due to interference between the + 1st order diffracted light and the ⁇ 1st order diffracted light.
  • these interferences are shearing interference in which a shear amount s is Nd in the case of a ⁇ / 2 phase grating and a shear amount s is 2Nd in the case of a ⁇ phase grating.
  • phase image of the subject 120 at the position of the phase grating 130 is W (x, y)
  • the phase ⁇ (x, y) and the phase image W (x, y) have the following relationship.
  • phase ⁇ (x, y) is information obtained by differentiating the phase image W (x, y) of the subject 120. Therefore, the phase image W (x, y) of the subject 120 can be obtained by integrating ⁇ (x, y).
  • FIG. 4 is a spectrum pattern of the interference intensity distribution when a one-dimensional grating is used. Usually, three peaks occur as shown in FIG. Central peak is the peak derived from predominantly A (f x, f y) to. On the other hand, the both sides of each peak C (f x, f y) and C * (f x, f y ) is the carrier frequency peak derived from, these peaks occur at the position of ⁇ f 0.
  • the peaks derived from f y ) or C * (f x , f y ) are separated.
  • phase ⁇ (x, y) that is, differential phase information is obtained from complex number information obtained by inverse Fourier transform.
  • 5A shows a checkered ⁇ / 2 phase grating (FIG. 2A) and a screened absorption grating (FIG. 3A), or a checkered absorption grating (FIG. 3A).
  • B)) is an example of the intensity distribution of moiré
  • 510 is a bright portion of moiré
  • 520 is a dark portion of moiré.
  • a checkered ⁇ phase grating (FIG. 2A) and a checkered absorption grating (FIG. 3B) are used, a moire intensity distribution is generated in such an oblique direction.
  • FIG. 5B shows the intensity distribution of moiré when a checkered ⁇ phase grating (FIG. 2A) and a screened absorption grating (FIG. 3A) are used.
  • Moire light part 540 is a moire dark part. In this case, a moire intensity distribution is generated in the vertical and horizontal directions.
  • moiré intensity distribution as described above can also be obtained using a screen door-shaped phase grating (FIG. 2B).
  • FIGS. 5C and 5D are spatial frequencies obtained by performing FFT (Fast Fourier Transform) processing, which is a kind of Fourier transform, on the moire intensity distribution shown in FIGS. 5A and 5B, respectively. It is a spectrum.
  • the maximum spatial frequency that can be calculated by FFT is 1 / (2P), where P is the pixel period of the X-ray detector 170.
  • Each of the two peaks 570 and 571 and 580 and 581 at orthogonal positions is cut out in the same manner as in the above-described one-dimensional case, and moved to the origin to perform inverse Fourier transform.
  • a broken line indicates a cutout region. Then, differential phase information in two orthogonal directions is obtained from the complex number information obtained by the inverse Fourier transform.
  • the differential phase information is in the direction of ⁇ 45 degrees
  • the differential phase information is in the X and Y directions.
  • the differential phase information obtained in this way is often convoluted (wrapped) between 2 ⁇ regions. That is, if the true phase at any point on a screen is ⁇ (x, y) and the convolved phase is ⁇ wrap (x, y)
  • n is an integer
  • ⁇ wrap (x, y) is determined so as to be in a region having a width of 2 ⁇ , for example, between 0 to 2 ⁇ and ⁇ to + ⁇ .
  • phase connection (unwrapping) of ⁇ (x, y) wrap is performed to restore the original ⁇ (x, y).
  • phase image W (x, y) of the subject can be obtained by integrating ⁇ (x, y) restored by the equation (8).
  • the integration direction can only be a direction orthogonal to the grid engraving direction. Therefore, in order to correctly measure the phase image W (x, y), the phase image is set so that one side parallel to the engraving direction of the X-ray detector 170 is irradiated with X-rays that do not pass through the subject 120. It is necessary to provide a portion where W (x, y) is known.
  • phase image W (x, y) can be correctly measured even when the entire surface of the X-ray detector 170 is irradiated with X-rays passing through the subject 120.
  • a peak separation step of cutting out a spectrum corresponding to the carrier frequency (a spectrum carrying the phase information) is performed (S631).
  • a peak separation step of cutting out a spectrum corresponding to the carrier frequency (a spectrum carrying the phase information) is performed (S631).
  • phase ⁇ (x, y) that is a differential phase image is obtained from the complex number information obtained in S632 (S633).
  • S631, S632, and S633 may be referred to as a phase recovery step (S630).
  • ⁇ (x, y) is wrapped, unwrap processing is performed to obtain true ⁇ (x, y) (S640). This is sometimes referred to as a phase connection process. If ⁇ (x, y) is not wrapped, this S640 can be omitted.
  • ⁇ (x, y) is differential phase information (differential phase image).
  • phase image W (x, y) is acquired by integrating ⁇ (x, y) (S650).
  • an X-ray imaging apparatus and an X-ray imaging method capable of acquiring a phase image of a subject from at least one imaging. It is also possible to provide a program that causes a computer to execute the above steps.
  • FIG. 7B shows the spatial frequency spectrum described in the present embodiment.
  • the basic period of the two-dimensional moire generated in the interference intensity distribution and the absorption grating is set to the pixel period of the X-ray detector.
  • the moiré direction is adjusted to be inclined 45 degrees with respect to the pixel array.
  • FIG. 7A is a diagram showing the moire intensity distribution on the X-ray detector in such a state.
  • 710 is a light receiving surface of the X-ray detector
  • 720 is a bright portion of moire
  • d is a moire cycle
  • P is a pixel cycle of the X-ray detector.
  • a checkered ⁇ / 2 phase grating (FIG. 2A) and a checkered absorption grating (FIG. 3B) are used.
  • other phase gratings and absorption gratings may be used.
  • FIG. 7B is a spatial frequency spectrum obtained by performing FFT processing on the moiré intensity distribution of FIG. 7A. If the number of pixel arrays is n both vertically and horizontally, the spectrum space obtained by FFT also becomes n ⁇ n discrete data. The maximum frequency that can be expressed is 1 / (2P) where P is the pixel period of the X-ray detector 170.
  • a carrier peak 711 that is a peak corresponding to the carrier frequency of the moiré intensity distribution is generated at the position.
  • the phase image of the subject can be restored.
  • the spatial resolution can be increased by cutting out the spectral region as wide as possible.
  • an unnecessary peak 721 that is the peak of the high frequency component and the DC component exists at the position of the sum or difference of the peak coordinates of the carrier frequency component.
  • the cutout spectrum region is a cutout region 731 inside the intermediate line between the carrier frequency peak and the unnecessary peak 721.
  • the spatial frequency of the phase image that can be restored in the present embodiment is 1 ⁇ 2 of the size of the cutout region 731. Therefore, as can be seen from FIG. 7B, the maximum frequency in the pixel array direction is 1 / (4P), and the maximum frequency in the 45 degree direction is
  • the cutout region in FIG. 7B is wider than the cutout region in FIG. Therefore, according to the present embodiment, it is possible to improve the spatial resolution as compared with the above embodiment.
  • FIG. 3 A third embodiment of the X-ray imaging apparatus according to the present invention will be described with reference to FIG. This embodiment is also devised to improve the spatial resolution compared to the spatial frequency spectrum of FIG. 5D described in the first embodiment.
  • FIG. 8B shows the spatial frequency spectrum described in the present embodiment.
  • the basic period of the two-dimensional moire generated in the interference intensity distribution and the absorption grating is set to three times the pixel period of the X-ray detector, and the direction of moire coincides with the pixel array.
  • FIG. 8 (A) is a diagram showing the moire intensity distribution on the X-ray detector in such a state.
  • 810 is a light receiving surface of the X-ray detector
  • 820 is a bright portion of moire
  • d is a moire cycle
  • P is a pixel cycle of the X-ray detector.
  • a checkered ⁇ phase grating (FIG. 2A) and a screen door absorption grating (FIG. 3A) are used.
  • other phase gratings and absorption gratings may be used.
  • FIG. 8A is a spatial frequency spectrum obtained by performing FFT processing on the moire intensity distribution shown in FIG. 8B.
  • the absolute value of the carrier frequency is 1 / (3P). Therefore, in spectral space
  • a carrier peak 811 that is a peak corresponding to the carrier frequency in the intensity distribution of the moire occurs.
  • the process described in the first embodiment is performed. The phase image of the subject can be restored.
  • the cutout spectral region is a cutout region 831 inside the intermediate line between the carrier frequency peak and the unnecessary peak 821.
  • the spatial frequency of the phase image that can be restored in this embodiment is 1 ⁇ 2 of the size of the cutout region 831. Therefore, from FIG. 8B, the maximum frequency in the pixel array direction is 1 / (6P), and the maximum frequency in the 45 degree direction is
  • the present embodiment is superior to the second embodiment in the spatial resolution in the direction of 45 degrees with respect to the pixel array.
  • the X-ray imaging apparatus according to the present embodiment is an X-ray imaging apparatus according to any one of Embodiments 1 to 3 equipped with an inspection object moving device 900.
  • the subject moving device 900 can move the subject 920 in the X-ray optical axis direction.
  • the imaging magnification of the subject 920 in the X-ray detector is L1 / L2 if the distance between the X-ray source 910 and the absorption grating 940 is L1, and the distance between the X-ray source 910 and the absorption grating 940 is L2.
  • L2 becomes large, so that imaging can be performed at a low magnification.
  • L2 becomes small, so that imaging can be performed at a high magnification.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Optics & Photonics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measurement Of Radiation (AREA)

Abstract

　少なくとも１回の撮像から、被検体の微分位相像または位相像を取得することができるＸ線撮像装置およびＸ線撮像方法を提供する。 　本発明に係るＸ線撮像装置は、位相格子１３０と吸収格子１５０と検出器１７０と演算装置１８０を有する。演算装置１８０は、検出器で取得したモアレの強度分布をフーリエ変換することにより、空間周波数スペクトルを取得するフーリエ変換工程を実行する。また、フーリエ変換工程により得られた空間周波数スペクトルから、キャリア周波数に対応したスペクトルを分離し、逆フーリエ変換を用いて微分位相像を取得する位相回復工程を実行する。

Description

Ｘ線撮像装置およびＸ線撮像方法
　本発明は、Ｘ線撮像装置およびＸ線撮像方法に関する。
　Ｘ線は高い物質透過性を持ち、かつ高空間分解能イメージングが可能であることから、工業的用途として物体の非破壊検査、医療用途としてレントゲン撮影、等に用いられている。これらは、物体や生体内の構成元素や密度差によりＸ線透過時の吸収の違いを利用してコントラスト画像を形成するものであり、Ｘ線吸収コントラスト法と言われる。しかし、軽元素ではＸ線吸収が非常に小さいため、生体の構成元素である炭素、水素、酸素などからなる生体軟組織、あるいはソフトマテリアルをＸ線吸収コントラスト法により画像化することは困難である。
　これに対して、軽元素で構成される組織でも明瞭な画像化を可能とする方法として、Ｘ線の位相差を用いたＸ線位相コントラスト法の研究が、１９９０年代より行なわれている。
　数多く開発されたＸ線位相コントラスト法のなかでも、通常のＸ線管を用いることが可能な方法として、タルボ干渉を利用したＸ線位相コントラスト法がある（特許文献１）。
　特許文献１に記載のタルボ干渉を利用した方式は、Ｘ線を発生するＸ線管、Ｘ線の位相を変調し干渉強度分布を発生させるための位相格子、および前記干渉強度分布をモアレの強度分布に変換する吸収格子、前記干渉強度分布を検出するＸ線検出器から構成される。
　特許文献１に記載の方式では、吸収格子を格子周期方向に走査することにより撮像を行なう。この走査により、検出されるモアレは移動し、走査量が吸収格子の１周期に達すると、モアレ画像は元の状態へ戻る。前記走査中の少なくとも３枚以上の撮像データを用いて演算処理することにより、微分位相像を取得する。
米国特許７１８０９７９号明細書
　上記特許文献１に記載の方法は、少なくとも３枚以上の撮像を行なうことにより微分位相像を取得し、その微分位相像から位相像を演算する方法である。
　しかしながら、特許文献１に記載の方法は３枚以上の撮像が必要になることから、撮像中に被検体が移動すると、画質が低下するという問題がある。
　また、撮像時間の長時間化は被検体へのＸ線線量の増加を招き、医療用の用途を考慮すると好ましくない。
　そこで、本発明は、少なくとも１回の撮像から、被検体の微分位相像または位相像を取得することができるＸ線撮像装置およびＸ線撮像方法の提供を目的とする。
　本発明に係るＸ線撮像装置は、Ｘ線源と、前記Ｘ線源からのＸ線を透過させることにより、タルボ効果による干渉強度分布を形成する位相格子と、前記位相格子により形成された干渉強度分布の一部を遮光することによりモアレを生じさせる吸収格子と、前記吸収格子により生じたモアレの強度分布を検出する検出器と、前記検出器により検出されたモアレの強度分布から被検体の情報を画像化し、出力する演算装置を有し、前記演算装置は、前記検出器で取得したモアレの強度分布をフーリエ変換することにより、空間周波数スペクトルを取得するフーリエ変換工程と、前記フーリエ変換工程により得られた空間周波数スペクトルから、キャリア周波数に対応したスペクトルを分離し、逆フーリエ変換を用いて微分位相像を取得する位相回復工程と、を実行するように構成されていることを特徴とする。
　本発明によれば、少なくとも１回の撮像から、被検体の微分位相像または位相像を取得することができるＸ線撮像装置およびＸ線撮像方法を提供することができる。
本発明の実施形態１に係るＸ線撮像装置を説明する図である。 本発明の実施形態２、および３に係る２次元位相格子を説明する図である。 本発明の実施形態１、および２に係る２次元吸収格子を説明する図である。 干渉強度分布のスペクトルパターンを表している図である。 ２次元位相格子を用いた場合のモアレの強度分布およびスペクトルパターンを表している図である。 本発明の演算装置で実施する解析方法のフローチャートを説明する図である。 本発明の実施形態２に係るモアレの強度分布および空間周波数スペクトルを説明する図である。 本発明の実施形態３に係るモアレの強度分布および空間周波数スペクトルを説明する図である。 本発明の実施形態４に係るズーム機構を説明する図である。
　（実施形態１）
　タルボ干渉を用いたＸ線撮像装置の構成の一例を図１に示す。Ｘ線撮像装置を用いて位相像を得るまでの工程を具体的に説明する。
　（Ｘ線源）
　Ｘ線源１１０から発生したＸ線１１１が被検体１２０を透過する。Ｘ線１１１は、被検体１２０を透過する際に、被検体１２０の組成や形状等に応じた位相の変化及び吸収を生じる。
　Ｘ線としては、連続Ｘ線を用いても特性Ｘ線を用いても良い。また、波長としては、０．１Åから５Å程度から適宜選択される。また、Ｘ線源１１０の後に、波長選択フィルタや線源用の格子を適宜設けても良い。
　（位相格子）
　前記被検体１２０を透過したＸ線１１１は、位相格子１３０を透過すると、タルボ効果により干渉強度分布１４０を形成する。
　位相格子１３０は、被検体１２０よりも上流側または下流側に配置される。
　位相格子１３０は、Ｘ線透過部材の厚みを周期的に変化させることにより構成された位相進行部分１３１と位相遅延部１３２からなる。位相進行部１３１と位相遅延部１３２は、透過したＸ線に対して位相差がつくように形成されていればよい。たとえば、位相進行部１３１を透過したＸ線の位相が位相遅延部１３２を透過したＸ線の位相に対してπだけ進行するように構成されている。厚みの変化量は、用いるＸ線の波長および部材により決定される。
　位相格子１３０は、位相遅延部１３２を透過したＸ線の位相に対して、位相進行部１３１を透過したＸ線の位相がπ、またはπ／２変調するものが一般的である。前者をπ位相格子、後者をπ／２位相格子と称する場合もある。ただし、位相の変調量は周期的であればよく、例えばπ／３変調するものであっても良い。
　位相格子１３０は、１次元の線状形状でもいいし、図２（Ａ）に示すような２次元の市松模様で構成されていてもよい。また、図２（Ｂ）に示すような網戸状の模様で構成されていてもよい。図２において、ｄは周期であり、２０１は２次元位相格子、２１０は位相進行部、２２０は位相遅延部を表している。
　なお、位相進行部２１０または位相遅延部２２０の形状は図２（Ａ）および図２（Ｂ）では正方形であるが、作製上、外縁が円形状に変形することもある。このように円形状に変形した場合においても、位相格子として利用することが可能である。
　位相格子１３０が１次元の周期を有する場合には、被検体１２０の１次元方向の位相勾配情報しか取得できない。しかし、位相格子１３０が２次元の周期を有する場合には、２次元方向の位相勾配情報を取得することができる点で有利である。
　なお、位相格子１３０を構成する材料はＸ線を透過する物質が好ましく、例えば、シリコン等を用いることができる。
　位相格子１３０の透過後に形成された干渉強度分布は、Ｘ線源と位相格子１３０との距離Ｚ０とすると、吸収格子１５０からの距離Ｚ１が下記の式（１）を満たしている位置に最も明瞭に現れる。ここで「干渉強度分布」とは、前記位相格子１３０の格子周期を反映した周期的強度分布のことである。
　式（１）において、λはＸ線の波長、ｄは位相格子１３０の格子周期である。
Figure JPOXMLDOC01-appb-M000005
 　Ｎは位相格子の形態により異なる値であり、以下のように表現できる実数である。なお、ｎは自然数である。
  １次元配列のπ位相格子：Ｎ＝ｎ／４－１／８
  １次元配列のπ／２位相格子：Ｎ＝ｎ－１／２
  ２次元配列の市松模様π位相格子：Ｎ＝ｎ／４－１／８
  ２次元で市松模様π／２位相格子：Ｎ＝ｎ／２－１／４
　（吸収格子）
　干渉強度分布の周期は、一般にＸ線検出器１７０の画素サイズより小さいため、そのままでは干渉強度分布を検出できない。そこで、吸収格子１５０を用いることによって、Ｘ線検出器１７０の画素サイズより大きな周期のモアレを発生させ、モアレの強度分布をＸ線検出器１７０で検出する。吸収格子１５０は、位相格子１３０から距離Ｚ１だけ離れた位置に設けることが好ましい。
　吸収格子１５０は、周期的に配列した透過部１５１と遮光部１５２を有し、位相格子１３０で形成された干渉強度分布１４０の明部の一部を遮光するように配置される。なお、透過部１５１はＸ線の一部が透過可能に構成されていればよく、開口が貫通している必要はない。吸収格子１５０を構成する材料はＸ線をよく吸収する物質であればよく、例えば、金を用いることができる。
　ここで、吸収格子１５０の周期は、干渉強度分布と同一、または僅かに異なる。
　干渉強度分布と同一周期を持つ吸収格子を用いた場合、吸収格子の面内回転によってモアレが発生する。干渉強度分布の周期をＤ、干渉強度分布の明暗の方位と吸収格子の方位のなす角をθ（ただし、θ＜＜１）とすれば、モアレの周期Ｄｍは、Ｄ／θとなる。
　一方、干渉強度分布と僅かに異なる周期を持つ吸収格子を用いた場合、吸収格子の面内回転を行なうことなくモアレが発生する。吸収格子の周期をＤａ＝Ｄ＋δ（ただしδ＜＜Ｄ）とすれば、モアレの周期Ｄｍは、Ｄ２／δとなる。
　吸収格子１５０において、透過部１５１は、１次元的に配列されていてもよいし、また２次元的に配列されていてもよい。
　例えば、図２（Ａ）に示した市松状の位相格子であって、π位相格子を用いた場合には、図３（Ａ）のように透過部３５１と遮光部３５２が２次元的に配列されている網戸状の吸収格子３００を用いる。また、図２（Ａ）に示した市松状のπ／２位相格子を用いた場合には、図３（Ｂ）のように透過部３５１と遮光部３５２が２次元的に配列されている市松状の吸収格子３００を用いる。
　なお、上記の組合せは一例であり、位相格子と吸収格子は種々の組合せが可能である。
　（検出器）
　前記吸収格子１４０を透過したＸ線の干渉強度分布の情報はモアレの強度分布として、Ｘ線検出器１７０によって検出される。Ｘ線検出器１７０は、Ｘ線の干渉強度分布の情報を検出することのできる素子である。例えば、デジタル信号への変換が可能な、ＦＰＤ（Ｆｌａｔ　Ｐａｎｅｌ　Ｄｅｔｅｃｔｏｒ）等を用いることができる。
　（演算装置）
　前記Ｘ線検出器１７０で検出されたモアレの強度分布の情報は、後述の解析方法による演算装置１８０によって解析され、微分位相像や位相像を画像化する。取得された微分位相像や位相像は、表示部１９０に表示するための出力画像である。演算装置１８０は例えばＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）を有する。
　以下、検出器で取得されるモアレの強度分布情報から位相像を取得する解析方法について説明し、その後に演算装置で行なう処理工程について説明する。
　（解析方法）
　干渉強度分布が形成される際には多数の回折光が重なって干渉するため、基本周波数（以後キャリア周波数と呼ぶ）と数多くのその高調波成分を含んでいる。しかし、モアレは干渉強度分布のキャリア周波数成分を空間的に拡大した形状を持つので、刻線がｘ軸に直交する１次元の位相格子を使用した場合にモアレは（２）式で表すことができる。
Figure JPOXMLDOC01-appb-M000006
　また、２次元の位相格子を使用した場合には、（２）式にｙ方向のキャリア周波数成分が重畳した形状となる。
　（２）式は、モアレが、バックグラウンドの第一項と、周期性を持った第二項との和で表されることを示している。ここで、ａ（ｘ，ｙ）はバックグラウンドを表し、ｂ（ｘ，ｙ）はキャリア周波数成分の振幅を表す。また、ｆ０は干渉縞のキャリア周波数を表し、φ（ｘ，ｙ）はキャリア周波数成分の位相を表している。
　キャリア周波数成分は、位相格子１３０として市松状のπ／２位相格子を用いた場合は、０次回折光と＋１次回折光の干渉、および０次回折光と－１次回折光の干渉によって生じる。また、位相格子１３０として市松状のπ位相格子を用いた場合は＋１次回折光と－１次回折光の干渉によってキャリア周波数成分が生じる。
　０次回折光と１次回折光は位相格子１３０で互いにＮｄだけ離れた位置にある光線が重なり、＋１次回折光と－１次回折光は位相格子１３０で互いに２Ｎｄだけ離れた位置にある光線が重なる。すなわち、これらの干渉はπ／２位相格子の場合はシア量ｓがＮｄ、π位相格子の場合はシア量ｓが２Ｎｄのシアリング干渉である。
　したがって、位相格子１３０の位置における被検体１２０の位相像をＷ（ｘ，ｙ）とすると、位相φ（ｘ，ｙ）と位相像Ｗ（ｘ，ｙ）は次式の関係がある。
Figure JPOXMLDOC01-appb-M000007
　一般にｓは微小なので、
Figure JPOXMLDOC01-appb-M000008
 となる。
　ここで、（３）式より位相φ（ｘ，ｙ）は被検体１２０の位相像Ｗ（ｘ，ｙ）を微分した情報となっていることがわかる。したがって、被検体１２０の位相像Ｗ（ｘ，ｙ）はφ（ｘ，ｙ）を積分することで求めることができる。
　ところで、位相φ（ｘ，ｙ）は、（２）式からフーリエ変換法で求めることができる。すなわち、（２）式は
Figure JPOXMLDOC01-appb-M000009
 と書き表すことができる。ここで
Figure JPOXMLDOC01-appb-M000010
 である。
　したがって、干渉縞からｃ（ｘ，ｙ）の成分またはｃ＊（ｘ，ｙ）の成分を取り出すことにより、位相φ（ｘ，ｙ）の情報を得ることが可能となる。
　ここで、フーリエ変換によって（４）式は
Figure JPOXMLDOC01-appb-M000011
 となる。
　ここでＧ（ｆｘ，ｆｙ）、Ａ（ｆｘ，ｆｙ）、Ｃ（ｆｘ，ｆｙ）はそれぞれｇ（ｘ，ｙ）、ａ（ｘ，ｙ）、ｃ（ｘ，ｙ）に対する二次元のフーリエ変換である。
　図４は、１次元格子を使用した場合の干渉強度分布のスペクトルパターンである。通常、図４のように３つのピークが生じる。中央のピークは主としてＡ（ｆｘ，ｆｙ）に由来するピークである。一方、その両脇のピークはそれぞれＣ（ｆｘ，ｆｙ）とＣ＊（ｆｘ，ｆｙ）に由来するキャリア周波数ピークであり、これらのピークは±ｆ０の位置に生じる。
　次に、Ｃ（ｆｘ，ｆｙ）あるいはＣ＊（ｆｘ，ｆｙ）に由来するピークのある領域を切り出す。例えば、Ａ（ｆｘ，ｆｙ）に由来するピークと、Ｃ（ｆｘ，ｆｙ）あるいはＣ＊（ｆｘ，ｆｙ）に由来するピークの周囲を切り出すことにより、Ｃ（ｆｘ，ｆｙ）あるいはＣ＊（ｆｘ，ｆｙ）に由来するピークを分離する。
　次に、分離したＣ（ｆｘ，ｆｙ）あるいはＣ＊（ｆｘ，ｆｙ）に由来するピークを周波数空間の原点に移動させ、逆フーリエ変換を行なう。そして逆フーリエ変換により得られた複素数情報から位相φ（ｘ，ｙ）、すなわち微分位相情報を得る。
　また、図５（Ａ）は、市松模様のπ／２位相格子（図２（Ａ））と、網戸模様の吸収格子（図３（Ａ））、または、市松模様の吸収格子（図３（Ｂ））を用いた場合のモアレの強度分布の例であり、５１０はモアレの明部、５２０はモアレの暗部である。なお、市松模様のπ位相格子（図２（Ａ））と、市松模様の吸収格子（図３（Ｂ））を用いた場合にも、このような斜め方向にモアレの強度分布が生じる。
　また、図５（Ｂ）は、市松模様のπ位相格子（図２（Ａ））と、網戸模様の吸収格子（図３（Ａ））を用いた場合のモアレの強度分布であり、５３０はモアレの明部、５４０はモアレの暗部である。この場合、縦横方向にモアレの強度分布が生じる。
　なお、網戸模様の位相格子（図２（Ｂ））を用いて、上記のようなモアレの強度分布を得ることもできる。
　図５（Ｃ）、（Ｄ）はそれぞれ、図５（Ａ）、（Ｂ）で示したモアレの強度分布をフーリエ変換の一種であるＦＦＴ（Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ）処理することによって得られた空間周波数スペクトルである。ＦＦＴで算出可能な最大空間周波数は、Ｘ線検出器１７０の画素周期をＰとすれば１／（２Ｐ）である。
　それぞれ、直交する位置にある２つのピ－ク５７０と５７１および５８０と５８１の周囲を、上述の１次元の場合と同様に切り出して、それぞれ原点に移動させ逆フーリエ変換を行なう。破線は切り出し領域を示している。そして逆フーリエ変換により得られた複素数情報から、直交する２方向の微分位相情報を得る。
　ここで、図５（Ｃ）場合は±４５度方向の微分位相情報となり、図５（Ｄ）場合はＸおよびＹ方向の微分位相情報となる。
　このようにして得られた微分位相情報は、２πの領域の間に畳み込まれている（ラップされている）ことが多い。すなわち、ある画面上の任意の点における真の位相をφ（ｘ，ｙ）とし畳み込まれた位相をφｗｒａｐ（ｘ，ｙ）とすると
Figure JPOXMLDOC01-appb-M000012
 の関係となる。ここでｎは整数であり、φｗｒａｐ（ｘ，ｙ）がある２πの幅を持った領域、例えば０から２πや－πから＋πの間に入るように決定される。
　これらの情報を元にして、φ（ｘ，ｙ）ｗｒａｐの位相接続（アンラップ）を行ない元のφ（ｘ，ｙ）に復元する。
　被検体の位相像Ｗ（ｘ，ｙ）は式（８）により復元したφ（ｘ，ｙ）を積分して求めることができる。
Figure JPOXMLDOC01-appb-M000013
　１次元格子を使用した場合は積分方向が格子刻線方向と直交する方向にしかできない。そのため、位相像Ｗ（ｘ，ｙ）を正しく計測するには、Ｘ線検出器１７０の前記刻線方向と平行な１辺に被検体１２０を通過しないＸ線が照射されるようにして位相像Ｗ（ｘ，ｙ）が既知の部分を設ける必要がある。
　一方、２次元格子を使用した場合は２方向の積分ができるのでＸ線検出器１７０全面に被検体１２０を通過したＸ線が照射されても正しく位相像Ｗ（ｘ，ｙ）を計測できる。
　（演算装置で行なう処理工程）
　上記を踏まえて、前記演算装置１８０において実行する処理フローの一例を図６に示す。
　まず、前記Ｘ線検出器１７０からモアレの強度分布の情報を取得する（Ｓ６１０）。
　次に、Ｓ６１０で得られたモアレの強度分布の情報をフーリエ変換し、空間周波数スペクトルを取得するフーリエ変換工程を行なう（Ｓ６２０）。
　次に、Ｓ６２０で得られた周波数空間において、キャリア周波数に対応したスペクトル（位相情報を担っているスペクトル）を切り出すピーク分離工程を行なう（Ｓ６３１）。なお、キャリア周波数に対応したスペクトル自体を切り出すのが難しい場合には、スペクトルの周辺領域の情報を切り出す。
　次に、Ｓ６３１で切り出されたスペクトルを周波数空間の原点に移動させ、逆フーリエ変換を行う（Ｓ６３２）。これにより、位相情報を有した複素数情報を得ることができる。
　次に、Ｓ６３２で得られた複素数情報から微分位相像である位相φ（ｘ，ｙ）を得る（Ｓ６３３）。なお、Ｓ６３１、Ｓ６３２、Ｓ６３３を位相回復工程ということもある（Ｓ６３０）。
　次に、φ（ｘ，ｙ）がラップされている場合には、アンラップ処理を行ない、真のφ（ｘ，ｙ）を得る（Ｓ６４０）。このことを位相接続工程ということもある。なお、φ（ｘ，ｙ）がラップされていない場合にはこのＳ６４０を省略することも可能である。ここで、φ（ｘ，ｙ）は微分位相情報（微分位相像）である。
　次に、φ（ｘ，ｙ）を積分することにより、位相像Ｗ（ｘ，ｙ）を取得する（Ｓ６５０）。
　以上のような構成によれば、少なくとも１枚の撮像から、被検体の位相像を取得することができるＸ線撮像装置およびＸ線撮像方法を提供することが可能となる。また、上記の工程をコンピュータに実行させるプログラムを提供することも可能となる。
　（実施形態２）
　本発明によるＸ線撮像装置の第２の実施形態を図７を用いて説明する。本実施形態では、実施形態１で説明した図５（Ｃ）の空間周波数スペクトルよりも空間分解能を向上させる工夫がなされている。
　図７（Ｂ）は、本実施形態で説明する空間周波数スペクトルを示したものである。このような周波数スペクトルを得るためには干渉強度分布と吸収格子で生じる２次元のモアレの基本周期をＸ線線検出器の画素周期の
Figure JPOXMLDOC01-appb-M000014
 とし、モアレの方位を画素配列に対し４５度傾けるように調整する。
　図７（Ａ）は、このような状態におけるＸ線検出器上のモアレの強度分布を示している図である。７１０はＸ線検出器の受光面、７２０はモアレの明部、ｄはモアレの周期、ＰはＸ線検出器の画素周期である。なお、本実施形態では、市松模様のπ／２位相格子（図２（Ａ））と、市松模様の吸収格子（図３（Ｂ））を用いている。ただし、発生するモアレの強度分布が同一ならば、他の位相格子および吸収格子を用いてもよい。
　図７（Ｂ）は図７（Ａ）のモアレの強度分布をＦＦＴ処理することによって得られた空間周波数スペクトルである。画素配列数を縦横ともにｎとすれば、ＦＦＴにより得られるスペクトル空間もｎ×ｎの離散的なデータとなる。また、表現できる最大の周波数は、Ｘ線検出器１７０の画素周期をＰとして１／（２Ｐ）である。
　本実施形態ではモアレの基本周期は
Figure JPOXMLDOC01-appb-M000015
 であるからその周波数であるキャリア周波数の絶対値は
Figure JPOXMLDOC01-appb-M000016
となる。またモアレの方位が４５度傾いているので、スペクトル空間上で
Figure JPOXMLDOC01-appb-M000017
 の位置にモアレの強度分布のキャリア周波数に対応するピークであるキャリアピーク７１１が生じる。
　４個のキャリアピーク７１１のうち隣り合う２個のピークについて、１辺が
Figure JPOXMLDOC01-appb-M000018
 である４５度傾いた正方形領域を切り出して実施形態１に記載の処理をすれば被検体の位相像を復元することができる。
　この際、スペクトル領域をなるべく広く切り出すことで空間分解能を大きくすることができる。しかしながら、スペクトル空間にはキャリア周波数のピークの他に、キャリア周波数成分のピーク座標の和または差の位置に高周波成分とＤＣ成分のピークである不要ピーク７２１が存在する。
　切り出し領域を広げすぎることにより、これらの不要ピーク７２１付近の領域まで含めてしまうと、不要ピーク７２１の影響で正確な位相像が得られなくなる。従って、切り出すスペクトル領域は、キャリア周波数のピークと不要ピーク７２１の中間線より内側の切り出し領域７３１としている。
　本実施形態で復元できる位相像の空間周波数は、切り出し領域７３１の大きさの１／２である。従って、図７（Ｂ）を見てわかるように画素配列方向の最大周波数は１／（４Ｐ）、４５度方向の最大周波数は
Figure JPOXMLDOC01-appb-M000019
 である。これを分解能として復元できる最小周期を画素単位で表現すれば、最大周波数の逆数になるので画素配列方向で４画素、４５度方向で
Figure JPOXMLDOC01-appb-M000020
 である。図５（Ｃ）の切り出し領域と比較すると、図７（Ｂ）の切り出し領域の方が、図５（Ｃ）の切り出し領域よりも広いため、復元できる空間周波数が大きい。したがって、本実施形態によれば、上記実施形態と比較して空間分解能を向上させることが可能となる。
　（実施形態３）
　本発明によるＸ線撮像装置の第３の実施形態を図８を用いて説明する。本実施形態も、実施形態１で説明した図５（Ｄ）の空間周波数スペクトルよりも空間分解能を向上させる工夫がなされている。
　図８（Ｂ）は、本実施形態で説明する空間周波数スペクトルを示したものである。このような周波数スペクトルを得るためには干渉強度分布と吸収格子で生じる２次元のモアレの基本周期をＸ線線検出器の画素周期の３倍とし、モアレの方位を画素配列に一致させる。
　図８（Ａ）は、このような状態におけるＸ線検出器上のモアレの強度分布を示している図である。８１０はＸ線検出器の受光面、８２０はモアレの明部、ｄはモアレの周期、ＰはＸ線検出器の画素周期である。なお、本実施形態では、市松模様のπ位相格子（図２（Ａ））と、網戸模様の吸収格子（図３（Ａ））を用いている。ただし、発生するモアレの強度分布が同一ならば他の位相格子および吸収格子を用いてもよい。
　図８（Ａ）は図８（Ｂ）に示したモアレの強度分布をＦＦＴ処理することによって得られた空間周波数スペクトルである。本実施形態ではモアレの基本周期は３Ｐであるから、そのキャリア周波数の絶対値は１／（３Ｐ）となる。従って、スペクトル空間上で
Figure JPOXMLDOC01-appb-M000021
 のモアレの強度分布のキャリア周波数に対応するピークであるキャリアピーク８１１が生じる。実施形態２と同様に４個のキャリアピーク８１１のうち隣り合う２個のピークのうち、１辺が１／（３Ｐ）の正立した正方形領域を切り出して実施形態１に記載の処理をすれば被検体の位相像を復元することができる。
　本実施形態でも、スペクトル空間にはキャリア周波数のピークの他に、キャリア周波数成分のピーク座標の和または差の位置に高周波成分とＤＣ成分のピークである不要ピーク８２１が存在する。そのため、切り出すスペクトル領域は、キャリア周波数のピークと不要ピーク８２１の中間線より内側の切り出し領域８３１としている。
　本実施形態で復元できる位相像の空間周波数は、切り出し領域８３１の大きさの１／２である。従って、図８（Ｂ）より、画素配列方向の最大周波数は１／（６Ｐ）、４５度方向の最大周波数は
Figure JPOXMLDOC01-appb-M000022
 である。これを分解能として復元できる最小周期を画素単位で表現すれば、最大周波数の逆数になるので画素配列方向で６画素、４５度方向で
Figure JPOXMLDOC01-appb-M000023
 である。したがって、画素配列に対して４５度方向の空間分解能において、本実施形態は実施形態２よりも優れている。
　（実施形態４）
　本発明によるＸ線撮像装置の第３の実施形態を図９を用いて説明する。本実施形態のＸ線撮像装置は、実施形態１から３記載のＸ線撮像装置に検査物移動装置９００を装備したものである。検査物移動装置９００によって被検体９２０をＸ線の光軸方向に移動させることができる。
　Ｘ線検出器における被検体９２０の撮像倍率は、Ｘ線源９１０と吸収格子９４０の距離をＬ１、Ｘ線源９１０と吸収格子９４０の距離をＬ２とすればＬ１／Ｌ２となる。
　そのため、被検体９２０を位相格子９３０に近づければＬ２が大きくなるため低倍率で撮像することができる。一方、Ｘ線源９１０に近づければＬ２が小さくなるので高倍率で撮像することができる。
　１１０　Ｘ線源
　１１１　Ｘ線
　１２０　被検体
　１３０　位相格子
　１５０　吸収格子
　１５１　透過部
　１５２　遮光部
　１７０　Ｘ線検出器
　１８０　演算装置

Claims (13)

  1. 　Ｘ線源と、
    　前記Ｘ線源からのＸ線を透過させることにより、タルボ効果による干渉強度分布を形成する位相格子と、
    　前記位相格子により形成された干渉強度分布の一部を遮光することによりモアレを生じさせる吸収格子と、
    　前記吸収格子により生じたモアレの強度分布を検出する検出器と、
    　前記検出器により検出されたモアレの強度分布から被検体の情報を画像化し、出力する演算装置を有し、
    　前記演算装置は、
    　前記検出器で取得したモアレの強度分布をフーリエ変換することにより、空間周波数スペクトルを取得するフーリエ変換工程と、
    　前記フーリエ変換工程により得られた空間周波数スペクトルから、キャリア周波数に対応したスペクトルを分離し、逆フーリエ変換を用いて微分位相像を取得する位相回復工程と、を実行するように構成されていることを特徴とするＸ線撮像装置。
  2. 　前記演算装置は、前記位相回復工程により得られた微分位相像を積分することにより位相像を取得する工程を実行するように構成されている請求項１に記載のＸ線撮像装置。
  3. 　前記演算装置は、前記位相回復工程により得られた微分位相像に対してアンラップ処理を行なう位相接続工程を実行するように構成されている請求項１に記載のＸ線撮像装置。
  4. 　前記位相格子は、位相進行部と位相遅延部とが２次元的に周期的に配されている請求項１に記載のＸ線撮像装置。
  5. 　前記位相格子は、前記位相進行部と前記位相遅延部とが市松状に配されている請求項４に記載のＸ線撮像装置。
  6. 　前記位相格子は、前記位相進行部を透過したＸ線と、前記位相遅延部を透過したＸ線との位相差がπ／２またはπとなるように構成されている請求項４に記載のＸ線撮像装置。
  7. 　前記フーリエ変換工程で得られたスペクトル空間上において、前記検出器の画素周期をＰとしたときに、
    Figure JPOXMLDOC01-appb-M000001
     の位置に前記キャリア周波数に対応したスペクトルが生じるように、前記位相格子、前記吸収格子、前記検出器が調整されていることを特徴とする請求項１に記載のＸ線撮像装置。
  8. 　前記吸収格子は、前記モアレの周期が
    Figure JPOXMLDOC01-appb-M000002
     でモアレの方位が前記検出器の画素配列に対し４５度となるように配置され、
    前記演算装置は、前記キャリア周波数に対応したスペクトルを分離するために、検出器で得られたモアレ画像の周波数空間から１辺が
    Figure JPOXMLDOC01-appb-M000003
     で、画素配列方向に対し４５度傾いた２つの正方形領域を切り出す工程を実行する請求項７に記載のＸ線撮像装置。
  9. 　前記フーリエ変換工程で得られたスペクトル空間上において、前記検出器の画素周期をＰとしたときに、
    Figure JPOXMLDOC01-appb-M000004
     の位置に前記キャリア周波数に対応したスペクトルが生じるように、前記位相格子、前記吸収格子、前記検出器が調整されていることを特徴とする請求項１に記載のＸ線撮像装置。
  10. 　前記吸収格子は、前記モアレの周期が３Ｐでモアレの方位が検出器の画素配列に一致するように配置され、
    　前記演算装置は、前記キャリア周波数に対応したスペクトルを分離するために、検出器で得られたモアレ画像の周波数空間から１辺が１／（３Ｐ）で、画素配列方向に対し正立した２つの正方形領域を切り出す工程を実行する請求項９に記載のＸ線撮像装置。
  11. 　前記被検体をＸ線の光軸方向に移動させることが可能な検査物移動装置を有する請求項１に記載のＸ線撮像装置。
  12. 　Ｘ線撮像装置において用いるＸ線撮像方法であって、
    　Ｘ線を透過させることにより、タルボ効果による干渉強度分布を形成する工程と、
    　前記干渉強度分布の一部を遮光することによりモアレを生じさせる工程と、
    　前記モアレの強度分布を検出する工程と、
    　前記モアレの強度分布をフーリエ変換することにより、空間周波数スペクトルを取得する工程と、
    　前記空間周波数スペクトルから、キャリア周波数に対応したスペクトルを分離し、逆フーリエ変換を用いて微分位相像を取得する工程と、
    　前記微分位相像を積分することにより、位相像を得る工程とを有することを特徴とするＸ線撮像方法。
  13. 　コンピュータに、
    　Ｘ線を透過させることにより、タルボ効果による干渉強度分布を形成する工程と、
    　前記干渉強度分布の一部を遮光することによりモアレを生じさせる工程と、
    　前記モアレの強度分布を検出する工程と、
    　前記モアレの強度分布をフーリエ変換することにより、空間周波数スペクトルを取得する工程と、
    　前記空間周波数スペクトルから、キャリア周波数に対応したスペクトルを分離し、逆フーリエ変換を用いて微分位相像を得る工程と、
    　前記微分位相像を積分することにより、位相像を得る工程と、
    を実行させることを特徴とするプログラム。
PCT/JP2009/068434 2008-10-29 2009-10-27 Ｘ線撮像装置およびｘ線撮像方法 Ceased WO2010050483A1 (ja)

Priority Applications (9)

Application Number Priority Date Filing Date Title
DE112009002606.0T DE112009002606B4 (de) 2008-10-29 2009-10-27 Röntgenstrahlabbildungsgerät und Röntgenstrahlabbildungsverfahren
EP09823593.0A EP2343537B1 (en) 2008-10-29 2009-10-27 X-ray imaging device and x-ray imaging method
KR1020117011759A KR101258927B1 (ko) 2008-10-29 2009-10-27 Ｘ선 촬상장치 및 ｘ선 촬상방법
JP2010535804A JP5174180B2 (ja) 2008-10-29 2009-10-27 Ｘ線撮像装置およびｘ線撮像方法
CN2009801428377A CN102197303A (zh) 2008-10-29 2009-10-27 X射线成像装置和x射线成像方法
US12/842,937 US8009797B2 (en) 2008-10-29 2010-07-23 X-ray imaging apparatus, X-ray imaging method, and X-ray imaging program
US13/085,199 US8559594B2 (en) 2008-10-29 2011-04-12 Imaging apparatus and imaging method
US13/190,770 US8340243B2 (en) 2008-10-29 2011-07-26 X-ray imaging apparatus, X-ray imaging method, and X-ray imaging program
US13/682,445 US8537966B2 (en) 2008-10-29 2012-11-20 X-ray imaging apparatus, X-ray imaging method, and X-ray imaging program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2008-278425 2008-10-29
JP2008278425 2008-10-29

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/842,937 Continuation US8009797B2 (en) 2008-10-29 2010-07-23 X-ray imaging apparatus, X-ray imaging method, and X-ray imaging program

Publications (1)

Publication Number Publication Date
WO2010050483A1 true WO2010050483A1 (ja) 2010-05-06

Family

ID=41664696

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/JP2009/068434 Ceased WO2010050483A1 (ja) 2008-10-29 2009-10-27 Ｘ線撮像装置およびｘ線撮像方法
PCT/JP2009/068863 Ceased WO2010050611A1 (en) 2008-10-29 2009-10-28 Analysis method, radiation imaging apparatus using analysis method, and analysis program for executing analysis method

Family Applications After (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/068863 Ceased WO2010050611A1 (en) 2008-10-29 2009-10-28 Analysis method, radiation imaging apparatus using analysis method, and analysis program for executing analysis method

Country Status (8)

Country Link
US (5) US8520799B2 (ja)
EP (2) EP2343537B1 (ja)
JP (4) JP5174180B2 (ja)
KR (1) KR101258927B1 (ja)
CN (3) CN102197303A (ja)
DE (1) DE112009002606B4 (ja)
RU (1) RU2519663C2 (ja)
WO (2) WO2010050483A1 (ja)

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012005820A (ja) * 2010-05-27 2012-01-12 Canon Inc Ｘ線撮像装置
WO2012023356A1 (ja) * 2010-08-19 2012-02-23 富士フイルム株式会社 放射線撮影システム及びその画像処理方法
WO2012038857A1 (en) * 2010-09-20 2012-03-29 Koninklijke Philips Electronics N.V. Phase gradient unwrapping in differential phase contrast imaging
WO2012056724A1 (ja) * 2010-10-29 2012-05-03 富士フイルム株式会社 放射線位相画像撮影装置
JP2012103237A (ja) * 2010-10-14 2012-05-31 Canon Inc 撮像装置
JP2012108098A (ja) * 2010-10-20 2012-06-07 Canon Inc トールボット干渉を用いた撮像装置および撮像装置の調整方法
JP2012143553A (ja) * 2010-12-24 2012-08-02 Fujifilm Corp 放射線画像撮影装置および放射線画像検出器
WO2012128335A1 (ja) * 2011-03-23 2012-09-27 コニカミノルタエムジー株式会社 医用画像表示システム
WO2012133553A1 (ja) * 2011-03-29 2012-10-04 富士フイルム株式会社 放射線撮影システム及び放射線撮影方法
JP2012187341A (ja) * 2011-03-14 2012-10-04 Canon Inc Ｘ線撮像装置
US20120275564A1 (en) * 2011-04-26 2012-11-01 Fujifilm Corporation Radiation imaging apparatus
WO2012147749A1 (ja) * 2011-04-25 2012-11-01 富士フイルム株式会社 放射線撮影システム及び放射線撮影方法
JP2012228371A (ja) * 2011-04-26 2012-11-22 Canon Inc 撮像装置
WO2012169427A1 (ja) * 2011-06-10 2012-12-13 富士フイルム株式会社 放射線撮影システム
WO2012169426A1 (ja) * 2011-06-08 2012-12-13 富士フイルム株式会社 放射線撮影システム
JP2013002845A (ja) * 2011-06-13 2013-01-07 Canon Inc 撮像装置、干渉縞解析プログラム及び干渉縞解析方法
WO2013027519A1 (ja) * 2011-08-22 2013-02-28 富士フイルム株式会社 放射線撮影装置及びアンラップ処理方法
JP2013050441A (ja) * 2011-08-03 2013-03-14 Canon Inc 波面測定装置、波面測定方法、及びプログラム並びにｘ線撮像装置
WO2013038881A1 (ja) * 2011-09-12 2013-03-21 富士フイルム株式会社 放射線撮影装置及び画像処理方法
WO2013047011A1 (ja) * 2011-09-30 2013-04-04 富士フイルム株式会社 放射線画像検出器及びその製造方法、並びに放射線画像検出器を用いた放射線撮影システム
EP2578156A1 (en) 2011-10-04 2013-04-10 Fujifilm Corporation Radiation imaging apparatus and image processing method
JP2013102951A (ja) * 2011-11-14 2013-05-30 Canon Inc 撮像装置および画像処理方法
JP2013536723A (ja) * 2010-09-03 2013-09-26 コーニンクレッカ フィリップス エヌ ヴェ 微分位相差イメージングにおける正規化位相回復
JP2013540031A (ja) * 2010-10-19 2013-10-31 コーニンクレッカ フィリップス エヌ ヴェ 微分位相コントラスト画像形成
JP2014108358A (ja) * 2012-11-30 2014-06-12 Canon Inc 逆リース変換による差分画像の合成
WO2014092206A1 (en) 2012-12-13 2014-06-19 Canon Kabushiki Kaisha Object information obtaining apparatus, program, and imaging system
US8767916B2 (en) 2011-04-20 2014-07-01 Fujifilm Corporation Radiation imaging apparatus and image processing method
JP2014121614A (ja) * 2012-12-24 2014-07-03 Canon Inc 位相画像を再構成する方法、コンピュータ可読記憶媒体、及び装置
EP2827339A1 (en) 2013-07-16 2015-01-21 Canon Kabushiki Kaisha Source grating, interferometer, and object information acquisition system
EP2924972A1 (en) 2014-03-27 2015-09-30 Canon Kabushiki Kaisha Image processing apparatus and imaging system
US9855018B2 (en) 2013-04-08 2018-01-02 Konica Minolta, Inc. Diagnostic medical image system and method of introducing Talbot capturing device to diagnostic medical image system used for general capturing
EP3530189A1 (en) 2018-02-23 2019-08-28 Konica Minolta, Inc. X-ray imaging system
WO2019240165A1 (ja) 2018-06-12 2019-12-19 国立大学法人筑波大学 位相画像撮影方法とそれを利用した位相画像撮影装置
US10725186B2 (en) 2016-02-23 2020-07-28 Canon Kabushiki Kaisha Scintillator plate, radiation detector, and radiation measurement system

Families Citing this family (76)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5339975B2 (ja) * 2008-03-13 2013-11-13 キヤノン株式会社 Ｘ線位相イメージングに用いられる位相格子、該位相格子を用いたｘ線位相コントラスト像の撮像装置、ｘ線コンピューター断層撮影システム
WO2010050483A1 (ja) * 2008-10-29 2010-05-06 キヤノン株式会社 Ｘ線撮像装置およびｘ線撮像方法
CN101943668B (zh) * 2009-07-07 2013-03-27 清华大学 X射线暗场成像系统和方法
JP5586899B2 (ja) * 2009-08-26 2014-09-10 キヤノン株式会社 Ｘ線用位相格子及びその製造方法
US8532252B2 (en) * 2010-01-27 2013-09-10 Canon Kabushiki Kaisha X-ray shield grating, manufacturing method therefor, and X-ray imaging apparatus
JP5631013B2 (ja) * 2010-01-28 2014-11-26 キヤノン株式会社 Ｘ線撮像装置
JP5725870B2 (ja) * 2010-02-22 2015-05-27 キヤノン株式会社 Ｘ線撮像装置およびｘ線撮像方法
JP5935693B2 (ja) * 2010-09-29 2016-06-15 コニカミノルタ株式会社 医用画像表示方法
RU2573114C2 (ru) * 2010-10-19 2016-01-20 Конинклейке Филипс Электроникс Н.В. Формирование изображений методом дифференциального фазового контраста
JP2012095865A (ja) * 2010-11-02 2012-05-24 Fujifilm Corp 放射線撮影装置、放射線撮影システム
JP2012200567A (ja) * 2011-03-28 2012-10-22 Fujifilm Corp 放射線撮影システム及び放射線撮影方法
JP2013024731A (ja) 2011-07-21 2013-02-04 Canon Inc 放射線検出装置
EP2586373B1 (en) * 2011-10-28 2014-12-03 CSEM Centre Suisse D'electronique Et De Microtechnique SA X-ray interferometer
US20150117599A1 (en) 2013-10-31 2015-04-30 Sigray, Inc. X-ray interferometric imaging system
CN104066375B (zh) * 2012-01-24 2017-08-11 皇家飞利浦有限公司 多方向相衬x射线成像
WO2013151342A1 (ko) * 2012-04-05 2013-10-10 단국대학교 산학협력단 방사선 위상차 영상 장치
JP2014006247A (ja) * 2012-05-28 2014-01-16 Canon Inc 被検体情報取得装置、被検体情報取得方法及びプログラム
JP6079204B2 (ja) * 2012-12-18 2017-02-15 コニカミノルタ株式会社 医用画像システム
US9357975B2 (en) 2013-12-30 2016-06-07 Carestream Health, Inc. Large FOV phase contrast imaging based on detuned configuration including acquisition and reconstruction techniques
US10578563B2 (en) 2012-12-21 2020-03-03 Carestream Health, Inc. Phase contrast imaging computed tomography scanner
US10096098B2 (en) 2013-12-30 2018-10-09 Carestream Health, Inc. Phase retrieval from differential phase contrast imaging
US9364191B2 (en) * 2013-02-11 2016-06-14 University Of Rochester Method and apparatus of spectral differential phase-contrast cone-beam CT and hybrid cone-beam CT
JP2014171799A (ja) * 2013-03-12 2014-09-22 Canon Inc Ｘ線撮像装置及びｘ線撮像システム
KR20140111818A (ko) 2013-03-12 2014-09-22 삼성전자주식회사 엑스선 영상 장치 및 그 제어 방법
JP2014178130A (ja) * 2013-03-13 2014-09-25 Canon Inc Ｘ線撮像装置及びｘ線撮像システム
US10295485B2 (en) 2013-12-05 2019-05-21 Sigray, Inc. X-ray transmission spectrometer system
DE102013221818A1 (de) * 2013-10-28 2015-04-30 Siemens Aktiengesellschaft Bildgebendes System und Verfahren zur Bildgebung
US9726622B2 (en) * 2013-10-31 2017-08-08 Tohoku University Non-destructive inspection device
USRE48612E1 (en) 2013-10-31 2021-06-29 Sigray, Inc. X-ray interferometric imaging system
CN105828716B (zh) 2013-12-17 2019-05-10 皇家飞利浦有限公司 用于扫描差分相位衬度系统的相位恢复
AU2013273822A1 (en) * 2013-12-23 2015-07-09 Canon Kabushiki Kaisha Modulation guided phase unwrapping
WO2015122542A1 (en) 2014-02-14 2015-08-20 Canon Kabushiki Kaisha X-ray talbot interferometer and x-ray talbot interferometer system
JP2015166676A (ja) * 2014-03-03 2015-09-24 キヤノン株式会社 Ｘ線撮像システム
JP6362914B2 (ja) * 2014-04-30 2018-07-25 キヤノンメディカルシステムズ株式会社 Ｘ線診断装置及び画像処理装置
US10401309B2 (en) 2014-05-15 2019-09-03 Sigray, Inc. X-ray techniques using structured illumination
CN104111120B (zh) * 2014-07-25 2017-05-31 中国科学院上海光学精密机械研究所 基于朗奇剪切干涉仪的相位提取方法
JP2016032573A (ja) * 2014-07-31 2016-03-10 キヤノン株式会社 トールボット干渉計、トールボット干渉システム、及び縞走査法
JP6688795B2 (ja) * 2014-11-24 2020-04-28 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｘ線位相コントラストトモシンセシス撮像に対する検出器及び撮像システム
US10117629B2 (en) 2014-12-03 2018-11-06 Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College High energy grating techniques
KR101636438B1 (ko) * 2015-03-18 2016-07-05 제이피아이헬스케어 주식회사 단일 그리드를 이용한 ｐｃｉ 기반의 엑스선 영상 생성 방법 및 그 장치
JP6604772B2 (ja) * 2015-08-05 2019-11-13 キヤノン株式会社 Ｘ線トールボット干渉計
JP6608246B2 (ja) * 2015-10-30 2019-11-20 キヤノン株式会社 Ｘ線回折格子及びｘ線トールボット干渉計
DE102015226571B4 (de) * 2015-12-22 2019-10-24 Carl Zeiss Smt Gmbh Vorrichtung und Verfahren zur Wellenfrontanalyse
JP6613988B2 (ja) * 2016-03-30 2019-12-04 コニカミノルタ株式会社 放射線撮影システム
DE102016206559B3 (de) * 2016-04-19 2017-06-08 Siemens Healthcare Gmbh Verfahren zur Korrektur eines Röntgenbilds auf Effekte eines Streustrahlenrasters, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
US10247683B2 (en) 2016-12-03 2019-04-02 Sigray, Inc. Material measurement techniques using multiple X-ray micro-beams
JP6937380B2 (ja) 2017-03-22 2021-09-22 シグレイ、インコーポレイテッド Ｘ線分光を実施するための方法およびｘ線吸収分光システム
CN109087348B (zh) * 2017-06-14 2022-04-29 北京航空航天大学 一种基于自适应区域投射的单像素成像方法
JP6838531B2 (ja) * 2017-09-06 2021-03-03 株式会社島津製作所 放射線位相差撮影装置
US10863958B2 (en) * 2017-10-11 2020-12-15 Shimadzu Corporation X-ray phase difference imaging system and phase contrast image correction method
WO2019111505A1 (ja) * 2017-12-06 2019-06-13 株式会社島津製作所 Ｘ線位相差撮像システム
US10578566B2 (en) 2018-04-03 2020-03-03 Sigray, Inc. X-ray emission spectrometer system
US10845491B2 (en) 2018-06-04 2020-11-24 Sigray, Inc. Energy-resolving x-ray detection system
CN112470245B (zh) 2018-07-26 2025-03-18 斯格瑞公司 高亮度x射线反射源
US10656105B2 (en) 2018-08-06 2020-05-19 Sigray, Inc. Talbot-lau x-ray source and interferometric system
CN112638261B (zh) 2018-09-04 2025-06-27 斯格瑞公司 利用滤波的x射线荧光的系统和方法
WO2020051221A2 (en) 2018-09-07 2020-03-12 Sigray, Inc. System and method for depth-selectable x-ray analysis
DE102019111463A1 (de) * 2019-05-03 2020-11-05 Wipotec Gmbh Röntgenstrahlungsdetektorvorrichtung und Vorrichtung zur Röntgeninspektion von Produkten, insbesondere von Lebensmitteln
BR112021019318A2 (pt) * 2019-05-30 2021-12-14 Becton Dickinson Co Correção de fase de sinais multiplexados de frequência de rádio
DE112020004169T5 (de) 2019-09-03 2022-05-25 Sigray, Inc. System und verfahren zur computergestützten laminografieröntgenfluoreszenz-bildgebung
CN111089871B (zh) * 2019-12-12 2022-12-09 中国科学院苏州生物医学工程技术研究所 X射线光栅相衬图像的相位信息分离方法及系统、储存介质、设备
US11175243B1 (en) 2020-02-06 2021-11-16 Sigray, Inc. X-ray dark-field in-line inspection for semiconductor samples
US11215572B2 (en) 2020-05-18 2022-01-04 Sigray, Inc. System and method for x-ray absorption spectroscopy using a crystal analyzer and a plurality of detector elements
US11549895B2 (en) 2020-09-17 2023-01-10 Sigray, Inc. System and method using x-rays for depth-resolving metrology and analysis
US12480892B2 (en) 2020-12-07 2025-11-25 Sigray, Inc. High throughput 3D x-ray imaging system using a transmission x-ray source
WO2022126071A1 (en) 2020-12-07 2022-06-16 Sigray, Inc. High throughput 3d x-ray imaging system using a transmission x-ray source
US12360067B2 (en) 2022-03-02 2025-07-15 Sigray, Inc. X-ray fluorescence system and x-ray source with electrically insulative target material
DE112023001408T5 (de) 2022-03-15 2025-02-13 Sigray, Inc. System und verfahren für die kompakte laminographie unter verwendung einer mikrofokus-transmissionsröntgenquelle und eines röntgendetektors mit variabler vergrösserung
CN115265524B (zh) * 2022-03-24 2024-10-25 北京控制工程研究所 一种临近空间x射线角秒级星跟踪器
JP2025515085A (ja) 2022-05-02 2025-05-13 シグレイ、インコーポレイテッド Ｘ線シーケンシャルアレイ波長分散型分光計
CN121013975A (zh) 2023-02-16 2025-11-25 斯格瑞公司 具有至少两个堆叠的平面布拉格衍射器的x射线探测器系统
US12181423B1 (en) 2023-09-07 2024-12-31 Sigray, Inc. Secondary image removal using high resolution x-ray transmission sources
WO2025101530A1 (en) 2023-11-07 2025-05-15 Sigray, Inc. System and method for x-ray absorption spectroscopy using spectral information from two orthogonal planes
CN117475172B (zh) * 2023-12-28 2024-03-26 湖北工业大学 一种基于深度学习的高噪声环境相位图解包裹方法和系统
US12429436B2 (en) 2024-01-08 2025-09-30 Sigray, Inc. X-ray analysis system with focused x-ray beam and non-x-ray microscope
WO2025174966A1 (en) 2024-02-15 2025-08-21 Sigray, Inc. System and method for generating a focused x‑ray beam

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003090807A (ja) * 2001-09-19 2003-03-28 Hitachi Ltd Ｘ線撮像法
WO2004058070A1 (ja) * 2002-12-26 2004-07-15 Atsushi Momose X線撮像装置および撮像方法
US20050117699A1 (en) 2003-11-28 2005-06-02 Akio Yoneyama X-ray imaging apparatus and x-ray imaging method
DE102006037257A1 (de) 2006-02-01 2007-08-02 Siemens Ag Verfahren und Messanordnung zur zerstörungsfreien Analyse eines Untersuchungsobjektes durch Röntgenstrahlung
JP2008026098A (ja) * 2006-07-20 2008-02-07 Hitachi Ltd Ｘ線撮像装置及び撮像方法
JP2008200359A (ja) * 2007-02-21 2008-09-04 Konica Minolta Medical & Graphic Inc Ｘ線撮影システム
JP2008200360A (ja) * 2007-02-21 2008-09-04 Konica Minolta Medical & Graphic Inc Ｘ線撮影システム
JP2009244260A (ja) * 2008-03-13 2009-10-22 Canon Inc Ｘ線位相イメージングに用いられる位相格子、該位相格子を用いたｘ線位相コントラスト像の撮像装置、ｘ線コンピューター断層撮影システム

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1673933A1 (ru) * 1988-12-28 1991-08-30 Ереванский политехнический институт им.К.Маркса Рентгеноинтерферометрический способ исследовани кристаллов
SU1748030A1 (ru) * 1990-06-07 1992-07-15 Ереванский политехнический институт им.К.Маркса Способ получени рентгеновских проекционных топограмм
US5424743A (en) * 1994-06-01 1995-06-13 U.S. Department Of Energy 2-D weighted least-squares phase unwrapping
US5812629A (en) * 1997-04-30 1998-09-22 Clauser; John F. Ultrahigh resolution interferometric x-ray imaging
JP2000088772A (ja) * 1998-09-11 2000-03-31 Hitachi Ltd Ｘ線撮像装置
US7268891B2 (en) * 2003-01-15 2007-09-11 Asml Holding N.V. Transmission shear grating in checkerboard configuration for EUV wavefront sensor
JP2006263180A (ja) * 2005-03-24 2006-10-05 Fuji Photo Film Co Ltd 画像処理装置およびこれを用いた放射線撮影システム
EP1731099A1 (en) * 2005-06-06 2006-12-13 Paul Scherrer Institut Interferometer for quantitative phase contrast imaging and tomography with an incoherent polychromatic x-ray source
DE102006063048B3 (de) * 2006-02-01 2018-03-29 Siemens Healthcare Gmbh Fokus/Detektor-System einer Röntgenapparatur zur Erzeugung von Phasenkontrastaufnahmen
CN101011260A (zh) * 2006-02-01 2007-08-08 西门子公司 用于识别和区分患者血管结构中的斑块的方法和ct系统
JP2008197593A (ja) * 2007-02-16 2008-08-28 Konica Minolta Medical & Graphic Inc Ｘ線用透過型回折格子、ｘ線タルボ干渉計およびｘ線撮像装置
WO2008102598A1 (ja) * 2007-02-21 2008-08-28 Konica Minolta Medical & Graphic, Inc. 放射線画像撮影装置及び放射線画像撮影システム
WO2010050483A1 (ja) * 2008-10-29 2010-05-06 キヤノン株式会社 Ｘ線撮像装置およびｘ線撮像方法
JP2010253157A (ja) * 2009-04-28 2010-11-11 Konica Minolta Medical & Graphic Inc Ｘ線干渉計撮影装置及びｘ線干渉計撮影方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003090807A (ja) * 2001-09-19 2003-03-28 Hitachi Ltd Ｘ線撮像法
WO2004058070A1 (ja) * 2002-12-26 2004-07-15 Atsushi Momose X線撮像装置および撮像方法
US20050286680A1 (en) 2002-12-26 2005-12-29 Atsushi Momose X-ray imaging system and imaging method
US20050117699A1 (en) 2003-11-28 2005-06-02 Akio Yoneyama X-ray imaging apparatus and x-ray imaging method
JP2005152500A (ja) * 2003-11-28 2005-06-16 Hitachi Ltd Ｘ線撮像装置及び撮像方法
DE102006037257A1 (de) 2006-02-01 2007-08-02 Siemens Ag Verfahren und Messanordnung zur zerstörungsfreien Analyse eines Untersuchungsobjektes durch Röntgenstrahlung
JP2008026098A (ja) * 2006-07-20 2008-02-07 Hitachi Ltd Ｘ線撮像装置及び撮像方法
JP2008200359A (ja) * 2007-02-21 2008-09-04 Konica Minolta Medical & Graphic Inc Ｘ線撮影システム
JP2008200360A (ja) * 2007-02-21 2008-09-04 Konica Minolta Medical & Graphic Inc Ｘ線撮影システム
JP2009244260A (ja) * 2008-03-13 2009-10-22 Canon Inc Ｘ線位相イメージングに用いられる位相格子、該位相格子を用いたｘ線位相コントラスト像の撮像装置、ｘ線コンピューター断層撮影システム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
M. JIANG ET AL.: "X-Ray Phase-Contrast Imaging with Three 2D Gratings", INTERNATIONAL JOURNAL OF BIOMEDICAL IMAGING, vol. 2008, 24 March 2008 (2008-03-24), pages 1 - 8, XP002671455, DOI: doi:10.1155/2008/827152
See also references of EP2343537A4

Cited By (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012005820A (ja) * 2010-05-27 2012-01-12 Canon Inc Ｘ線撮像装置
US8824629B2 (en) 2010-08-19 2014-09-02 Fujifilm Corporation Radiation imaging system and image processing method
JP2012061300A (ja) * 2010-08-19 2012-03-29 Fujifilm Corp 放射線撮影システム及びその画像処理方法
WO2012023356A1 (ja) * 2010-08-19 2012-02-23 富士フイルム株式会社 放射線撮影システム及びその画像処理方法
CN103068311A (zh) * 2010-08-19 2013-04-24 富士胶片株式会社 放射线摄影系统及其图像处理方法
JP2013536723A (ja) * 2010-09-03 2013-09-26 コーニンクレッカ フィリップス エヌ ヴェ 微分位相差イメージングにおける正規化位相回復
WO2012038857A1 (en) * 2010-09-20 2012-03-29 Koninklijke Philips Electronics N.V. Phase gradient unwrapping in differential phase contrast imaging
JP2012103237A (ja) * 2010-10-14 2012-05-31 Canon Inc 撮像装置
JP2013540031A (ja) * 2010-10-19 2013-10-31 コーニンクレッカ フィリップス エヌ ヴェ 微分位相コントラスト画像形成
JP2012108098A (ja) * 2010-10-20 2012-06-07 Canon Inc トールボット干渉を用いた撮像装置および撮像装置の調整方法
US9194825B2 (en) 2010-10-20 2015-11-24 Canon Kabushiki Kaisha Imaging apparatus using talbot interference and adjusting method for imaging apparatus
US8781069B2 (en) 2010-10-29 2014-07-15 Fujifilm Corporation Radiographic phase-contrast imaging apparatus
CN103188996A (zh) * 2010-10-29 2013-07-03 富士胶片株式会社 放射线照相相衬成像设备
WO2012056724A1 (ja) * 2010-10-29 2012-05-03 富士フイルム株式会社 放射線位相画像撮影装置
JP5796908B2 (ja) * 2010-10-29 2015-10-21 富士フイルム株式会社 放射線位相画像撮影装置
CN103188996B (zh) * 2010-10-29 2015-06-24 富士胶片株式会社 放射线照相相衬成像设备
JP2012143553A (ja) * 2010-12-24 2012-08-02 Fujifilm Corp 放射線画像撮影装置および放射線画像検出器
JP2012187341A (ja) * 2011-03-14 2012-10-04 Canon Inc Ｘ線撮像装置
JP5915645B2 (ja) * 2011-03-23 2016-05-11 コニカミノルタ株式会社 医用画像表示システム
WO2012128335A1 (ja) * 2011-03-23 2012-09-27 コニカミノルタエムジー株式会社 医用画像表示システム
WO2012133553A1 (ja) * 2011-03-29 2012-10-04 富士フイルム株式会社 放射線撮影システム及び放射線撮影方法
US8767916B2 (en) 2011-04-20 2014-07-01 Fujifilm Corporation Radiation imaging apparatus and image processing method
WO2012147749A1 (ja) * 2011-04-25 2012-11-01 富士フイルム株式会社 放射線撮影システム及び放射線撮影方法
JP2012228371A (ja) * 2011-04-26 2012-11-22 Canon Inc 撮像装置
US20120275564A1 (en) * 2011-04-26 2012-11-01 Fujifilm Corporation Radiation imaging apparatus
WO2012169426A1 (ja) * 2011-06-08 2012-12-13 富士フイルム株式会社 放射線撮影システム
WO2012169427A1 (ja) * 2011-06-10 2012-12-13 富士フイルム株式会社 放射線撮影システム
JP2013002845A (ja) * 2011-06-13 2013-01-07 Canon Inc 撮像装置、干渉縞解析プログラム及び干渉縞解析方法
JP2013050441A (ja) * 2011-08-03 2013-03-14 Canon Inc 波面測定装置、波面測定方法、及びプログラム並びにｘ線撮像装置
WO2013027519A1 (ja) * 2011-08-22 2013-02-28 富士フイルム株式会社 放射線撮影装置及びアンラップ処理方法
WO2013038881A1 (ja) * 2011-09-12 2013-03-21 富士フイルム株式会社 放射線撮影装置及び画像処理方法
WO2013047011A1 (ja) * 2011-09-30 2013-04-04 富士フイルム株式会社 放射線画像検出器及びその製造方法、並びに放射線画像検出器を用いた放射線撮影システム
EP2578156A1 (en) 2011-10-04 2013-04-10 Fujifilm Corporation Radiation imaging apparatus and image processing method
JP2013102951A (ja) * 2011-11-14 2013-05-30 Canon Inc 撮像装置および画像処理方法
JP2014108358A (ja) * 2012-11-30 2014-06-12 Canon Inc 逆リース変換による差分画像の合成
JP2014117336A (ja) * 2012-12-13 2014-06-30 Canon Inc 演算装置、プログラム及び撮像システム
US20150308967A1 (en) * 2012-12-13 2015-10-29 Canon Kabushiki Kaisha Object information obtaining apparatus, program, and imaging system
WO2014092206A1 (en) 2012-12-13 2014-06-19 Canon Kabushiki Kaisha Object information obtaining apparatus, program, and imaging system
US9702831B2 (en) 2012-12-13 2017-07-11 Canon Kabushiki Kaisha Object information obtaining apparatus, program, and imaging system
JP2014121614A (ja) * 2012-12-24 2014-07-03 Canon Inc 位相画像を再構成する方法、コンピュータ可読記憶媒体、及び装置
US9855018B2 (en) 2013-04-08 2018-01-02 Konica Minolta, Inc. Diagnostic medical image system and method of introducing Talbot capturing device to diagnostic medical image system used for general capturing
EP2827339A1 (en) 2013-07-16 2015-01-21 Canon Kabushiki Kaisha Source grating, interferometer, and object information acquisition system
EP2924972A1 (en) 2014-03-27 2015-09-30 Canon Kabushiki Kaisha Image processing apparatus and imaging system
US10725186B2 (en) 2016-02-23 2020-07-28 Canon Kabushiki Kaisha Scintillator plate, radiation detector, and radiation measurement system
EP3530189A1 (en) 2018-02-23 2019-08-28 Konica Minolta, Inc. X-ray imaging system
WO2019240165A1 (ja) 2018-06-12 2019-12-19 国立大学法人筑波大学 位相画像撮影方法とそれを利用した位相画像撮影装置
US11181487B2 (en) 2018-06-12 2021-11-23 University Of Tsukuba Phase imaging method and phase imaging apparatus using phase imaging method

Also Published As

Publication number Publication date
EP2343537A1 (en) 2011-07-13
US8009797B2 (en) 2011-08-30
EP2343537A4 (en) 2012-05-02
US20110158493A1 (en) 2011-06-30
RU2519663C2 (ru) 2014-06-20
US20100290590A1 (en) 2010-11-18
JP5174180B2 (ja) 2013-04-03
US20130070893A1 (en) 2013-03-21
JP2012507690A (ja) 2012-03-29
KR20110079742A (ko) 2011-07-07
CN102197302B (zh) 2013-06-26
CN103876761B (zh) 2016-04-27
CN102197303A (zh) 2011-09-21
US8681934B2 (en) 2014-03-25
US8340243B2 (en) 2012-12-25
JP2013106962A (ja) 2013-06-06
KR101258927B1 (ko) 2013-04-29
US20130301795A1 (en) 2013-11-14
JP5456032B2 (ja) 2014-03-26
JP5595473B2 (ja) 2014-09-24
US8520799B2 (en) 2013-08-27
JPWO2010050483A1 (ja) 2012-03-29
DE112009002606B4 (de) 2024-02-01
CN103876761A (zh) 2014-06-25
CN102197302A (zh) 2011-09-21
DE112009002606T5 (de) 2012-08-02
EP2343537B1 (en) 2019-04-10
US20110280368A1 (en) 2011-11-17
US8537966B2 (en) 2013-09-17
RU2012137875A (ru) 2014-03-27
EP2342551A1 (en) 2011-07-13
WO2010050611A1 (en) 2010-05-06
JP2014205079A (ja) 2014-10-30

Similar Documents

Publication Publication Date Title
JP5174180B2 (ja) Ｘ線撮像装置およびｘ線撮像方法
JP5777360B2 (ja) Ｘ線撮像装置
JP5538936B2 (ja) 解析方法、プログラム、記憶媒体、ｘ線位相イメージング装置
JP5631013B2 (ja) Ｘ線撮像装置
JP5586986B2 (ja) Ｘ線撮像装置
JP2013063099A (ja) Ｘ線撮像装置
JP2012005820A (ja) Ｘ線撮像装置
JP6537293B2 (ja) Ｘ線トールボット干渉計及びｘ線トールボット干渉計システム
JP2013042983A (ja) トモシンセシス撮像装置及びトモシンセシス画像の撮像方法
JP6673188B2 (ja) Ｘ線位相撮影装置
JP2013164339A (ja) Ｘ線撮像装置
WO2015156379A1 (en) Image processing unit and control method for image processing unit
JP7180566B2 (ja) Ｘ線イメージング装置およびｘ線イメージング方法
JP2017090414A (ja) 二次元干渉パターン撮像装置
RU2472137C1 (ru) Устройство формирования рентгеновских изображений и способ формирования рентгеновских изображений
JPWO2018168621A1 (ja) 放射線画像生成装置
JP2011237773A (ja) 撮像装置及び撮像方法
WO2013051647A1 (ja) 放射線撮影装置及び画像処理方法
JP2020038153A (ja) 放射線画像生成装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980142837.7

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09823593

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2010535804

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2009823593

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 1120090026060

Country of ref document: DE

ENP Entry into the national phase

Ref document number: 20117011759

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 3616/CHENP/2011

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 2011121670

Country of ref document: RU