WO2020066135A1 - X線位相イメージング装置 - Google Patents
X線位相イメージング装置 Download PDFInfo
- Publication number
- WO2020066135A1 WO2020066135A1 PCT/JP2019/022384 JP2019022384W WO2020066135A1 WO 2020066135 A1 WO2020066135 A1 WO 2020066135A1 JP 2019022384 W JP2019022384 W JP 2019022384W WO 2020066135 A1 WO2020066135 A1 WO 2020066135A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- information
- relative angle
- ray
- phase contrast
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 107
- 238000012545 processing Methods 0.000 claims description 59
- 238000010191 image analysis Methods 0.000 claims description 37
- 230000003287 optical effect Effects 0.000 claims description 14
- 230000000153 supplemental effect Effects 0.000 claims 1
- 239000000835 fiber Substances 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000010521 absorption reaction Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 9
- 238000012805 post-processing Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 238000000034 method Methods 0.000 description 6
- 238000000333 X-ray scattering Methods 0.000 description 5
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 239000004918 carbon fiber reinforced polymer Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 239000011152 fibreglass Substances 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
-
- 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/06—Diaphragms
-
- 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
- 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/046—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 tomography, e.g. computed tomography [CT]
Definitions
- the present invention relates to an X-ray phase imaging apparatus, and more particularly, to an X-ray phase imaging apparatus including a rotation mechanism that changes a relative angle between a subject and a plurality of gratings around an X-ray optical axis.
- an X-ray phase imaging apparatus imaging is performed while changing the relative angle between a subject and a plurality of gratings around the optical axis of X-rays, and a plurality of imagings having different relative angles between the obtained subject and the plurality of gratings.
- a technique for performing image analysis from a result is known.
- Such an X-ray phase imaging apparatus is disclosed in, for example, Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems”, PHYSICAL REVIEW B 82, 214103 (2010).
- the Jensen T H, et al., “Directional x-ray dark-field imaging of strong ordered systems” PHYSICAL REVIEW B 82, 214103 (2010) has an X-ray source, a plurality of gratings, It has a detector and an image processing device. X-rays from the X-ray source pass through a plurality of gratings via the subject and irradiate the detector. An absorption image, a phase differential image, and a dark-field image are generated by the image processing device based on the signal detected by the detector.
- the X-ray phase imaging apparatus focuses on the X-ray optical axis.
- an imaging method for changing an angle of a subject with respect to a plurality of grids is disclosed.
- a dark-field image or the like (hereinafter simply referred to as an image) is generated by the image processing device.
- the image processing apparatus is configured to perform image analysis such as the directivity and angle of minute angle scattering based on each of the plurality of generated images and the corresponding relative angle.
- an advantage of some aspects of the invention is that image analysis is performed based on an incorrect relative angle between a subject and a plurality of grids. Accordingly, an object of the present invention is to provide an X-ray phase imaging apparatus capable of suppressing information of analysis results from becoming inappropriate due to the above.
- an X-ray phase imaging apparatus includes an X-ray source, a detector that detects X-rays emitted from the X-ray source, an X-ray source and a detector. And a plurality of gratings including a first grating for forming a self-image by X-rays emitted from the X-ray source, and a second grating for causing interference with the self-image of the first grating.
- Rotating at least one of a subject disposed between the X-ray source and the detector and the plurality of gratings about the optical axis of the X-ray to change the relative angle between the subject and the plurality of gratings A first rotation mechanism; and an image processing unit that generates a phase contrast image based on a signal detected by the detector at each of the plurality of relative angles changed by the rotation of the first rotation mechanism.
- the relative angle information is stored in the corresponding position. It is configured to store in association with the contrast image.
- the X-ray phase imaging apparatus by storing information of each of a plurality of relative angles in association with a corresponding phase contrast image, image analysis of the plurality of phase contrast images is performed. At this time, image analysis can be performed using the relative angles stored in association with each phase contrast image. As a result, unlike the case where the corresponding relative angle is calculated from each phase contrast image, it is possible to suppress the use of the wrong relative angle in the image analysis. Thus, it is possible to prevent the information of the analysis result from being inappropriate due to the image analysis being performed based on the erroneous relative angle between the subject and the plurality of grids.
- the relative angle information stored in association with each of the plurality of phase contrast images is read, and the image is formed based on the read relative angle information. It is configured to perform the analysis.
- image analysis can be performed only by reading out the relative angle. The control load on the phase imaging device can be reduced.
- the phase contrast image and the information on the relative angle are associated with each other by including the information on the relative angle in each file name of the plurality of phase contrast images. It is configured to memorize. With this configuration, unlike the case where the corresponding phase contrast image and the information on the relative angle are stored in a common folder to associate them with each other, it is not necessary to newly generate a folder. Further, unlike the case where the information on the relative angles of the corresponding phase contrast images is described in a file to associate them with each other, it is not necessary to newly generate a file.
- phase contrast image with the information on the relative angle only by giving the information on the relative angle to the file name without generating a new folder and file, so that the data used in the memory etc. An increase in capacity can be suppressed.
- the supplementary data section included in each of the plurality of phase contrast images includes information on the corresponding relative angle, thereby obtaining information on the phase contrast image and the relative angle.
- the corresponding phase contrast image and the information on the relative angle are stored in a common folder to associate them with each other, it is not necessary to newly generate a folder.
- the information on the relative angles of the corresponding phase contrast images is described in a file to associate them with each other, it is not necessary to newly generate a file.
- phase contrast image with the information on the relative angle only by providing the information on the relative angle in the supplementary data section without generating a new folder and file, and thus it is used in a memory or the like. An increase in data capacity can be suppressed.
- each of the plurality of phase contrast images is stored in a separate folder from each other, and information on the relative angle corresponding to each of the plurality of phase contrast images is shared.
- the phase contrast image and the information on the relative angle are stored in association with each other.
- a management file different from the plurality of phase contrast images which indicates a relationship between each of the plurality of phase contrast images and information on a corresponding relative angle, is created.
- the phase contrast image and the information on the relative angle are stored in association with each other.
- the data capacity of the phase contrast image can be reduced as compared with the case where the information of the relative angle is provided in the phase contrast image itself.
- by associating the phase contrast image with the information on the relative angle based on the management file it is possible to collectively associate each of the plurality of dark field images with the information on the corresponding relative angle by the management file. it can. As a result, it is not necessary to generate a file for describing information on the corresponding relative angle for each dark-field image, so that the number of files stored in a memory or the like can be reduced.
- the X-ray phase imaging apparatus further includes a control unit that controls the first rotation mechanism so as to change a relative angle between the subject and the plurality of gratings, and the plurality of the plurality of grids generated by the image processing unit. Is stored in association with the information on the relative angle transmitted from the control unit.
- the control load on the image processing unit can be reduced as compared with the case where both the generation of the phase contrast image and the acquisition of the information on the relative angle are performed by the image processing unit.
- the subject and an imaging system including an X-ray source, a detector, and a plurality of gratings are relatively rotated about a direction orthogonal to the optical axis.
- the image processing unit further includes a second rotation mechanism that causes the first rotation mechanism to rotate the subject and causes the second rotation mechanism to relatively move the subject and the imaging system at each of the plurality of relative angles changed.
- a three-dimensional phase contrast image is generated at each of the plurality of relative angles, and information of each of the plurality of relative angles is stored in association with a corresponding three-dimensional phase contrast image. It is configured. With this configuration, when a three-dimensional phase contrast image is generated using the second rotation mechanism, image analysis is performed based on an incorrect relative angle between the subject and the plurality of grids. Thus, it is possible to prevent the information of the analysis result from becoming inappropriate.
- the relative angle information stored in association with each of the plurality of three-dimensional phase contrast images is read, and the corresponding three-dimensional phase contrast image is rotated based on the read relative angle information.
- the plurality of three-dimensional phase contrast images are configured to be oriented in the same direction.
- the three-dimensional phase contrast image generated at each of the plurality of relative angles is displaced by an amount corresponding to the difference in the relative angles, so that image analysis cannot be performed properly.
- the corresponding three-dimensional phase contrast images are rotated about the optical axis, and the orientations of the plurality of three-dimensional phase contrast images are aligned to thereby obtain a plurality of rotated three-dimensional phase contrast images.
- Image analysis can be appropriately performed using the two-dimensional phase contrast image.
- the image processing unit is configured to generate a dark field image at each of the plurality of relative angles
- the plurality of dark field images generated by the image processing unit Is configured to calculate the relationship between the orientation angle of the subject and the scattering intensity by performing image analysis based on each of the above and information on the corresponding relative angle.
- FIG. 1 is a diagram illustrating a configuration of an X-ray phase imaging apparatus according to a first embodiment.
- FIG. 3 is a diagram showing a subject imaged by the X-ray phase imaging apparatus according to the first and second embodiments. 5 is a dark-field image taken at each of a plurality of relative angles by the X-ray phase imaging apparatus according to the first embodiment.
- FIG. 3 is a diagram illustrating data of a plurality of dark-field images stored in a memory in the X-ray phase imaging apparatus according to the first embodiment.
- FIG. 3 is a diagram illustrating a relationship between pixel values and a relative angle of a dark-field image captured by the X-ray phase imaging apparatus according to the first embodiment.
- FIG. 2 is a diagram illustrating images acquired from a plurality of dark-field images captured by the X-ray phase imaging apparatus according to the first embodiment.
- FIG. 6A is a diagram illustrating an image of the orientation angle of the subject
- FIG. 6B is a diagram illustrating an image of the anisotropy of X-ray scattering in the subject.
- FIG. 5 is a flowchart illustrating control in the image processing unit according to the first embodiment. It is a figure showing the composition of the X-ray phase imaging device by a 2nd embodiment.
- FIG. 11 is a flowchart illustrating control in an image processing unit according to a second embodiment.
- FIG. 11 is a diagram illustrating a plurality of dark-field images stored in a memory in the X-ray phase imaging apparatus according to the first modification.
- FIG. 13 is a diagram illustrating a file in which data of a dark field image and information of a relative angle stored in a memory in an X-ray phase imaging apparatus according to a second modified example.
- FIG. 13 is a diagram showing dark field image data and a management file stored in a memory in an X-ray phase imaging apparatus according to a third modification.
- the X-ray phase imaging apparatus 100 generates a dark-field image of the subject T using the diffusion (scattering) of the X-rays that have passed through the subject T. Specifically, the X-ray phase imaging apparatus 100 generates a dark-field image of the subject T using the Talbot effect.
- the X-ray phase imaging apparatus 100 can be used for imaging the inside of an object in, for example, non-destructive inspection applications.
- the dark field image is an example of the “phase contrast image” in the claims.
- the subject T includes the fiber bundle 11 (see FIG. 2) inside.
- the subject T is, for example, a carbon fiber reinforced plastic (CFRP) in which carbon fibers are used as the fiber bundle 11.
- CFRP carbon fiber reinforced plastic
- the fiber bundle is a bundle of a large number of fibers.
- the fiber bundle 11 is formed in a plate shape by a large number of fibers.
- the subject T is arranged between an X-ray source 1 and a detector 5 described later. Specifically, the subject T is arranged between a first grid 2 and a third grid 4 described later.
- FIG. 1 is a diagram of the X-ray phase imaging apparatus 100 viewed from the X direction.
- the X-ray phase imaging apparatus 100 includes an X-ray source 1, a first grating 2, a second grating 3, a third grating 4, a detector 5, a PC 6, a control unit 7, , A grid rotating mechanism 8 and a grid moving mechanism 9.
- a direction from the X-ray source 1 toward the first grating 2 is defined as a Z2 direction
- a direction opposite thereto is defined as a Z1 direction.
- a left-right direction in a plane orthogonal to the Z direction is defined as an X direction
- a direction toward the back of the paper is defined as an X2 direction
- a direction toward the near side of the paper is defined as an X1 direction
- a vertical direction in a plane orthogonal to the Z direction is defined as a Y direction
- an upward direction is defined as a Y1 direction
- a downward direction is defined as a Y2 direction.
- the lattice rotating mechanism 8 is an example of the “first rotating mechanism” in the claims.
- the X-ray source 1 is configured to generate X-rays when a high voltage is applied, and to irradiate the generated X-rays in the Z2 direction.
- the first grating 2 has a plurality of slits 2a arranged at a predetermined period (pitch) d1 in the Y direction and an X-ray phase changing section 2b. Each slit 2a and X-ray phase changing portion 2b are formed so as to extend linearly. The slits 2a and the X-ray phase changing portions 2b are formed so as to extend in parallel.
- the first grating 2 is a so-called phase grating.
- the first grating 2 is disposed between the X-ray source 1 and the second grating 3, and is irradiated with X-rays from the X-ray source 1.
- the first grating 2 is provided for forming a self-image (indicated by a broken line in FIG. 1) of the first grating 2 by the Talbot effect.
- a self-image indicated by a broken line in FIG. 1
- Talbot distance This is called the Talbot effect.
- the second grating 3 has a plurality of X-ray transmitting portions 3a and X-ray absorbing portions 3b arranged at a predetermined period (pitch) d2 in the Y direction.
- the X-ray absorbing section 3b extends along the direction in which the X-ray phase changing section 2b extends.
- Each of the X-ray transmitting portions 3a and the X-ray absorbing portions 3b is formed to extend linearly. Further, each of the X-ray transmitting portions 3a and the X-ray absorbing portions 3b is formed so as to extend in parallel.
- the second grating 3 is a so-called absorption grating.
- the first grating 2 and the second grating 3 are gratings having different roles, respectively, but the slit 2a and the X-ray transmitting portion 3a respectively transmit X-rays.
- the X-ray absorbing section 3b has a role of blocking X-rays, and the X-ray phase changing section 2b changes the phase of X-rays depending on the difference in refractive index from the slit 2a.
- the second grating 3 is arranged between the first grating 2 and the detector 5, and is irradiated with X-rays that have passed through the first grating 2.
- the second grating 3 is arranged at a position away from the first grating 2 by a Talbot distance.
- the second grating 3 interferes with the self-image of the first grating 2 and forms moiré fringes (not shown) on the detection surface of the detector 5.
- the third grating 4 has a plurality of X-ray transmitting portions 4a and X-ray absorbing portions 4b arranged at a predetermined period (pitch) d3 in the Y direction. Each of the X-ray transmitting portions 4a and the X-ray absorbing portions 4b is formed so as to extend linearly. Further, each of the X-ray transmitting portions 4a and the X-ray absorbing portions 4b is formed so as to extend in parallel.
- the third grating 4 is a so-called multi-slit.
- the third grating 4 is arranged between the X-ray source 1 and the first grating 2.
- the third grating 4 is configured to convert the X-rays from the X-ray source 1 into a multi-point light source by using the X-rays that have passed through the respective X-ray transmitting sections 4a as a line light source.
- the pitch of the three gratings (the first grating 2, the second grating 3, and the third grating 4) and the distance between the gratings satisfy certain conditions, the X-rays emitted from the X-ray source 1 It is possible to increase coherence. Thereby, the interference intensity can be maintained even if the focal size of the tube of the X-ray source 1 is large.
- the detector 5 is configured to detect X-rays, convert the detected X-rays into electric signals, and read the converted electric signals as image signals.
- the detector 5 is, for example, an FPD (Flat @ Panel @ Detector).
- the detector 5 includes a plurality of conversion elements (not shown) and pixel electrodes (not shown) arranged on the plurality of conversion elements.
- the plurality of conversion elements and the pixel electrodes are arranged in an array in the X direction and the Y direction at a predetermined cycle (pixel pitch).
- the detector 5 is configured to output the acquired image signal to the PC 6.
- the PC 6 includes an image processing unit 60 and a memory 61.
- the image processing unit 60 has an image generation unit 60a and a post-processing unit 60b.
- the image processing unit 60 (image generation unit 60a) is configured to generate an absorption image (not shown) based on the image signal output from the detector 5.
- the image processing unit 60 (image generation unit 60a) generates a dark-field image (see FIG. 3) due to X-ray scattering based on the X-ray intensity distribution detected by the detector 5.
- the absorption image is an image of a contrast generated by a difference in X-ray absorption by the subject T.
- the dark field image is an image of a contrast generated by minute angle scattering of X-rays by a fine structure inside the subject T.
- the image processing unit 60 includes, for example, a CPU (Central Processing Unit) and a processor such as a GPU (Graphics Processing Unit) and an FPGA (Field-Programmable Gate Array) configured for image processing.
- the control unit 7 calculates the relative angle ⁇ between the subject T and the first grid 2, the second grid 3, and the third grid 4 (hereinafter, three grids are collectively referred to as a plurality of grids).
- the lattice rotating mechanism 8 is controlled to change. Specifically, the relative angle ⁇ is changed by rotating a plurality of gratings around the optical axis ⁇ of the X-ray by the grating rotating mechanism 8. A plurality of gratings may be rotated simultaneously or sequentially.
- the grating rotating mechanism 8 includes a rotating mechanism 8a for rotating the first grating 2, a rotating mechanism 8b for rotating the second grating 3, and a rotating mechanism 8c for rotating the third grating 4.
- the X-ray phase imaging apparatus 100 includes a rotating grating unit 20 including the first grating 2, the rotating mechanism 8a, and the grating moving mechanism 9. Further, the X-ray phase imaging apparatus 100 includes a rotating grating unit 30 including the second grating 3 and the rotating mechanism 8b. Further, the X-ray phase imaging apparatus 100 includes a rotating grating unit 40 including the third grating 4 and the rotating mechanism 8c.
- the controller 7 is configured to control the grid moving mechanism 9 to move the first grid 2 stepwise in the grid plane in a direction orthogonal to the grid direction.
- the X-ray phase imaging apparatus 100 uses a technique (fringe scanning method) of acquiring a reconstructed image from a plurality of moiré fringes (images) obtained by scanning the first grating 2 at regular intervals.
- Control unit 7 includes a processor such as a CPU (Central Processing Unit), for example.
- the lattice direction is a direction in which the lattice pattern of the lattice extends.
- the grating pattern is a slit 2a, an X-ray phase changing unit 2b, an X-ray transmitting unit 3a, an X-ray absorbing unit 3b, and the like of each grating.
- the lattice moving mechanism 9 is configured to move the first lattice 2 stepwise in a direction orthogonal to the lattice direction in the lattice plane (XY plane) (Y direction in FIG. 1) based on a signal from the control unit 7. Have been. Specifically, the grating moving mechanism 9 divides the period d1 of the first grating 2 by n, and moves the first grating 2 stepwise by d1 / n. The grating moving mechanism 9 is configured to move the first grating 2 stepwise at least for one period d1 of the first grating 2. Note that n is a positive integer, for example, 9 or the like.
- the grid moving mechanism 9 includes, for example, a stepping motor, a piezo actuator, and the like. In the case where the first grating 2 is rotated by the grating rotating mechanism 8, the direction of the step movement changes along with the direction orthogonal to the grating direction.
- the image processing unit 60 performs darkening based on the signal detected by the detector 5 at each of the plurality of relative angles ⁇ changed by the rotation of the grid rotation mechanism 8. Generate a field image. Specifically, the image processing unit 60 (image generation unit 60a) generates a dark-field image in each of the states where the relative angle ⁇ is 0 degree, 45 degrees, 90 degrees, and 135 degrees. Note that the relative angle ⁇ is an angle between the direction of the plurality of gratings (for example, the direction in which the X-ray phase changing unit 2b extends) and the X direction. Further, as shown in FIG. 3, in each dark-field image, the scattering intensity of the fiber bundle 11 extending in parallel with the directions of the plurality of lattices is displayed strong (darker black).
- the X-ray phase imaging apparatus 100 is configured to store information of each of the plurality of relative angles ⁇ in association with the corresponding dark field image. Specifically, the data of each dark field image associated with the information of the corresponding relative angle ⁇ is stored (stored) in the memory 61 of the PC 6.
- the X-ray phase imaging apparatus 100 includes the information of the relative angle ⁇ in each of the file names of the plurality of dark-field images, so that the The image and the information on the relative angle are stored in association with each other.
- the user of the X-ray phase imaging apparatus 100 arbitrarily specifies the file name (eg, test_0deg).
- the X-ray phase imaging apparatus 100 reads information on the relative angle ⁇ stored in association with each of the plurality of dark field images, and performs image analysis based on the read information on the plurality of relative angles ⁇ . It is configured as follows. Specifically, the image processing unit 60 (post-processing unit 60b) searches for and reads out the file names of the data of the plurality of dark field images stored in the memory 61, and associates the data with the data of each dark field image. Of the obtained relative angle ⁇ .
- the X-ray phase imaging apparatus 100 performs image analysis based on each of the plurality of dark-field images generated by the image processing unit 60 (image generation unit 60a) and information on the corresponding relative angle ⁇ . , And is configured to calculate the relationship between the orientation angle of the subject T and the scattering intensity.
- the image processing unit 60 (post-processing unit 60b) plots the value of the pixel value (the value based on the scattering intensity) with respect to the relative angle ⁇ in each pixel of the dark field image. . Then, the image processing unit 60 (post-processing unit 60b) approximates the waveform created by plotting by the expression A + Bcos [2 ( ⁇ ) + ⁇ ].
- A is offset information in the waveform.
- B is information on the amplitude in the waveform.
- ⁇ which is the orientation angle of the fiber bundle 11 of the subject T, is calculated based on the values of A, B, and ⁇ .
- FIG. 6A shows an image obtained by imaging the calculated orientation angle ⁇ .
- FIG. 6B shows an image obtained by imaging B / A, which is the anisotropy of X-ray scattering. Then, information indicating the direction of the fiber bundle 11 (see FIG. 7) can be obtained by coloring the pixel positions of only the portion having a high numerical value of B / A with ⁇ , which is information on the orientation angle, and displaying them.
- FIG. 7 shows that the color of the image of the fiber bundle 11 changes in accordance with the orientation angle of the fiber bundle 11. Thereby, it is possible to analyze the orientation of the fiber bundle 11 visually and quantitatively.
- step S1 a dark-field image is generated at each of the plurality of relative angles ⁇ .
- step S2 information of the relative angle ⁇ is included in the file names of the plurality of dark field images and stored (stored) in the memory 61.
- step S3 information on the relative angle ⁇ is read from each file name of the plurality of dark field images stored (stored) in the memory 61, and based on the dark field image and the read relative angle ⁇ . Perform image analysis.
- the X-ray phase imaging apparatus 100 is configured to store information on each of the plurality of relative angles ⁇ in association with the corresponding dark-field image. Accordingly, when performing image analysis on a plurality of dark field images, image analysis can be performed using the relative angle ⁇ stored in association with each dark field image. As a result, unlike the case where the corresponding relative angle ⁇ is calculated from each dark field image, it is possible to suppress the use of the wrong relative angle in the image processing. As a result, it is possible to prevent the information of the analysis result from becoming inappropriate due to the image analysis being performed based on the erroneous relative angle ⁇ between the subject T and the plurality of grids.
- the information on the relative angle ⁇ stored in association with each of the plurality of dark field images is read, and the image based on the read information on the plurality of relative angles ⁇ is read.
- the X-ray phase imaging apparatus 100 is configured to perform the analysis.
- image analysis can be performed only by reading out the relative angle ⁇ .
- the control load on the imaging device 100 can be reduced.
- the file name of each of the plurality of dark field images includes the information of the corresponding relative angle ⁇ , thereby associating the dark field image with the information of the relative angle.
- the X-ray phase imaging apparatus 100 is configured to store the information. Accordingly, unlike the case where the corresponding dark-field image and the information of the relative angle ⁇ are stored in a common folder and are associated with each other, it is not necessary to newly generate a folder. Also, unlike the case where the information on the relative angle ⁇ of the corresponding dark-field image is described in a file to associate them with each other, it is not necessary to newly generate a file.
- the plurality of dark-field images generated by the image processing unit 60 and the information on the corresponding relative angles ⁇ transmitted from the control unit 7 are stored in association with each other.
- the X-ray phase imaging apparatus 100 The control load on the image processing unit 60 can be reduced as compared with the case where both the generation of the dark field image and the acquisition of the information on the relative angle ⁇ are performed by the image processing unit 60.
- the subject T is analyzed by performing image analysis based on each of the plurality of dark-field images generated by the image processing unit 60 and information on the corresponding relative angle ⁇ .
- the X-ray phase imaging apparatus 100 is configured so as to calculate the relationship between the orientation angle and the scattering intensity. Accordingly, when performing image analysis on a plurality of dark field images, image analysis can be performed using the relative angle ⁇ stored in association with each dark field image. As a result, unlike the case where the corresponding relative angle ⁇ is calculated from each dark field image, it is possible to suppress the use of the wrong relative angle ⁇ in the image analysis.
- the X-ray phase imaging apparatus 200 of the second embodiment is different from the X-ray phase imaging apparatus 100 of the first embodiment in that it is configured to generate a three-dimensional dark field image at each of a plurality of relative angles. I have.
- the same components as those in the first embodiment are denoted by the same reference numerals as those in the first embodiment, and are not described.
- the X-ray phase imaging apparatus 200 includes a PC 16, a control unit 17, a subject rotation mechanism 18, and a CT rotation mechanism.
- the PC 16 includes an image processing unit 160 and a memory 161.
- the image processing section 160 has an image generation section 160a and a post-processing section 160b.
- the CT rotation mechanism 28 includes, for example, a rotation stage 28a driven by a motor or the like.
- the subject rotation mechanism 18 and the CT rotation mechanism 28 are examples of the “first rotation mechanism” and the “second rotation mechanism” in the claims, respectively.
- the subject rotating mechanism 18 changes the relative angle ⁇ between the subject T and the plurality of grids by rotating the subject T about the optical axis ⁇ of the X-ray.
- the CT rotation mechanism 28 relatively rotates the subject T, the X-ray source 1, the detector 5, and the imaging system 300 including a plurality of gratings based on a signal from the control unit 17. It is configured. Specifically, the CT rotation mechanism 28 is configured to rotate the subject T about an axis ⁇ orthogonal to the optical axis ⁇ , thereby rotating the subject T relative to the imaging system 300.
- the X-ray phase imaging apparatus 200 includes a grating unit 120 including the first grating 2 and the grating moving mechanism 9. Further, the X-ray phase imaging apparatus 200 includes a grating unit 130 including the second grating 3. Further, the X-ray phase imaging apparatus 200 includes a grating unit 140 including the third grating 4. Further, the X-ray phase imaging apparatus 200 includes a subject holding unit 150 including the subject rotation mechanism 18 and the CT rotation mechanism 28.
- the image processing unit 16 reconstructs a plurality of absorption images and a plurality of dark-field images captured at each of the plurality of rotation angles while rotating the CT rotation mechanism 28. Generates a three-dimensional absorption image (not shown) and a three-dimensional dark-field image (not shown). Note that the three-dimensional dark-field image is an example of the “three-dimensional phase contrast image” in the claims.
- the image processing unit 16 (the image generation unit 160 a) performs the CT rotation at each of the plurality of relative angles ⁇ changed by rotating the subject T by the subject rotation mechanism 18.
- a three-dimensional dark-field image is generated at each of the plurality of relative angles ⁇ .
- the image processing unit 16 (image generation unit 160a) generates a three-dimensional dark field image in each of the cases where the relative angle ⁇ is 45 °, 0 °, and ⁇ 45 °.
- the image of the subject T is shown in a simplified manner.
- the X-ray phase imaging apparatus 200 is configured to store information of each of the plurality of relative angles ⁇ in association with a corresponding three-dimensional dark field image. Specifically, the data of each three-dimensional dark field image associated with the information of the corresponding relative angle ⁇ is stored (stored) in the memory 161 of the PC 16.
- the X-ray phase imaging apparatus 200 reads the information on the relative angle ⁇ stored in association with each of the plurality of three-dimensional dark field images, and also reads the information on the read relative angle ⁇ .
- the corresponding three-dimensional dark-field images are rotated about the optical axis ⁇ to align the directions of the plurality of three-dimensional dark-field images.
- the image processing unit 160 (post-processing unit 160b) is configured to perform image analysis using each slice image of the plurality of three-dimensional dark field images in a state where the directions of the plurality of three-dimensional dark field images are aligned. ing.
- the slice image means a cross-sectional image in a three-dimensional dark field image.
- a plurality of slice images are acquired in a three-dimensional dark field image.
- the orientation of the fiber bundle 11 can be visually and quantitatively analyzed for each of the plurality of slice images.
- step S11 a three-dimensional dark field image is generated at each of the plurality of relative angles ⁇ .
- step S12 the file names of the plurality of three-dimensional dark field images are stored (stored) in the memory 161 with the corresponding information of the relative angle ⁇ included.
- step S13 information on the relative angle ⁇ is read from each file name of the plurality of three-dimensional dark-field images, and at least two of the three-dimensional dark-field images are rotated to align their directions.
- step S14 image analysis is performed for each slice image of the three-dimensional dark field image based on the information of the associated relative angle ⁇ .
- the X-ray phase imaging apparatus 200 is configured to store the information of each of the plurality of relative angles ⁇ in association with the corresponding three-dimensional dark field image. Thereby, when a three-dimensional dark-field image is generated using the CT rotation mechanism 28, image analysis is performed based on an incorrect relative angle ⁇ between the subject T and a plurality of grids. It is possible to prevent the information of the analysis result from becoming inappropriate.
- the information of the relative angle ⁇ stored in association with each of the plurality of three-dimensional dark field images is read, and based on the read information of the relative angle ⁇ ,
- the X-ray phase imaging apparatus 200 is configured to rotate the corresponding three-dimensional dark-field image so that the directions of the plurality of three-dimensional dark-field images are aligned.
- the three-dimensional dark field image generated at each of the plurality of relative angles ⁇ is displaced by an amount corresponding to the difference between the relative angles ⁇ , so that image analysis cannot be appropriately performed. Then, based on the read information on the relative angle ⁇ , the corresponding three-dimensional dark-field images are rotated to align the directions of the plurality of three-dimensional dark-field images. Image processing can be appropriately performed using the image.
- the X-ray phase imaging apparatus 400 includes a header included in each of the plurality of dark-field images, by including information of the corresponding relative angle ⁇ , thereby obtaining a
- the information on the relative angle ⁇ is stored in association with the information.
- the dark field image including the information of the relative angle ⁇ in the header is stored (stored) in the memory 261.
- the header is not actually displayed in the dark field image, it is illustrated in FIG. 13 as being displayed for explanation.
- the header is an example of the “supplementary data section” in the claims.
- the X-ray phase imaging apparatus 500 stores each of the plurality of dark-field images in a separate folder from each other, and stores information on each of the plurality of dark-field images and the corresponding relative angle ⁇ .
- a folder 361a, a folder 361b, a folder 361c, and a folder 361d are created in the memory 361.
- the file 500c
- ° and a file 500b are stored.
- the files 500b to 500d are the same as the file 500a, and a detailed description will be omitted.
- the X-ray phase imaging apparatus 600 includes a management file 600a that indicates the relationship between each of the plurality of dark-field images and the information of the corresponding relative angle ⁇ and is different from the plurality of dark-field images.
- the dark field image and the information on the relative angle ⁇ are stored in association with each other.
- the memory 461 stores a management file 600a (for example, a text file or a binary file) in which the correspondence between each of the plurality of dark-field images and the relative angle ⁇ is described.
- data of a plurality of dark field images may be stored in separate folders, and the correspondence between the folder name and the relative angle ⁇ may be described in the management file 600a.
- the management file 600a collectively associates each of the plurality of dark field images with the corresponding information of the relative angle ⁇ . It can be performed. As a result, it is not necessary to generate a file for describing the information of the corresponding relative angle ⁇ for each dark-field image, so that the number of files stored in the memory 461 can be reduced.
- information on a plurality of relative angles ⁇ is stored in association with a corresponding dark-field image
- the present invention is not limited to this.
- information on a plurality of relative angles ⁇ may be stored in association with an absorption image or a phase differential image.
- the relative angle ⁇ is changed by rotating a plurality of gratings, but the present invention is not limited to this.
- the relative angle ⁇ may be changed by rotating the subject T.
- both the plurality of grids and the subject T may be rotated.
- the present invention is not limited to this.
- a plurality of gratings may be rotated to change the relative angle ⁇ .
- the example in which the subject T is arranged between the first grid 2 and the third grid 4 has been described. May be arranged between them.
- CT imaging may be performed by rotating the imaging system 300.
- the present invention is not limited to this.
- the third grating 4 may not be provided.
- a dark-field image is generated by the fringe scanning method
- the present invention is not limited to this.
- a method of rotating any one of the first grating 2, the second grating 3, and the third grating 4 on a plane orthogonal to the direction in which the optical axis ⁇ of the X-ray extends (so-called single moire shooting) Method), a dark field image may be generated.
- the first grating 2 is a phase grating
- the present invention is not limited to this.
- the first grating 2 may be an absorption grating.
- the present invention is not limited to this. Any of the plurality of grids may be moved stepwise.
- CFRP carbon fiber reinforced plastic
- GFRP glass fiber reinforced plastic
- the relative angle ⁇ between the subject T and the plurality of grids is included in the file name of the data of the dark field image, but the present invention is not limited to this.
- both the subject T and the rotation angles of each of the plurality of grids may be included in the file name.
- only the rotation angle of the rotation may be included in the file name.
- the image processing unit 60 (160) associates the corresponding dark-field image with the relative angle ⁇
- the present invention is not limited to this.
- the control unit 7 (17) may associate the corresponding dark-field image with the relative angle ⁇ .
- the processing of the image processing unit 60 (160) has been described with reference to the flow-driven type flowchart for convenience of description, but the present invention is not limited to this.
- the processing of the image processing unit 60 (160) may be performed by an “event driven type” in which the processing is performed in event units. In this case, it may be performed in a completely event-driven manner, or may be performed in a combination of event-driven and flow-driven.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Surgery (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Chemical & Material Sciences (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- Animal Behavior & Ethology (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
このX線位相イメージング装置(100)は、被写体(T)と複数の格子との相対角度(θ)の各々の情報を、対応する暗視野像に関連付けて記憶するように構成されている。
Description
本発明は、X線位相イメージング装置に関し、特に、X線の光軸を中心に被写体と複数の格子との相対角度を変化させる回転機構を備えるX線位相イメージング装置に関する。
従来、X線位相イメージング装置において、X線の光軸を中心に被写体と複数の格子との相対角度を変化させて撮影し、得られた被写体と複数の格子との相対角度が異なる複数の撮影結果から画像解析する手法が知られている。このようなX線位相イメージング装置は、たとえば、Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)に開示されている。
上記Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)のX線位相イメージング装置は、X線源と、複数の格子と、検出器と、画像処理装置とを備えている。X線源からのX線は、被写体を介して、複数の格子を通過して検出器に照射される。検出器により検出される信号に基づいて、画像処理装置により吸収像、位相微分像、および、暗視野像が生成される。
また、上記Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)のX線位相イメージング装置では、X線の光軸を中心に、被写体の複数の格子に対する角度を変化させる撮影手法が開示されている。上記撮影手法における回転により変化された被写体と複数の格子との複数の相対角度の各々において、画像処理装置により暗視野像等(以下、単に画像とする)が生成される。そして、画像処理装置は、生成した複数の画像の各々と、対応する相対角度とに基づいて、微少角散乱の指向性や角度など画像解析を行うように構成されている。
Jensen T H, et al., "Directional x-ray dark-field imaging of strongly ordered systems" PHYSICAL REVIEW B 82, 214103 (2010)
ここで、複数の画像を画像解析する際には、各々の画像に対応する被写体と複数の格子との相対角度の情報が必要である。そして、上記Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)には明記されていないが、上記Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)に記載されているような従来のX線位相イメージング装置では、複数の画像の各々を解析して、各々の画像が生成された際の被写体と複数の格子との相対角度が算出される場合がある。上記Jensen T H, et al., “Directional x-ray dark-field imaging of strongly ordered systems” PHYSICAL REVIEW B 82, 214103 (2010)に記載されているように被写体を光軸回転させる場合、被写体がたとえば円状や球状のような回転させても形状が変わらないような形状である場合、相対角度の算出(予測)は困難である。また被写体ではなく、格子を回転させる場合、得られる画像の被写体の形状は変化せず、相対角度を算出(予測)することは困難である。このため、被写体と複数の格子との相対角度が誤って算出される場合もある。この場合、誤った相対角度に基づいて画像解析が行われるので、解析結果の情報が不適切になるという問題点がある。
この発明は、上記のような課題を解決するためになされたものであり、この発明の1つの目的は、被写体と複数の格子との間の誤った相対角度に基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することが可能なX線位相イメージング装置を提供することである。
上記目的を達成するために、この発明の一の局面におけるX線位相イメージング装置は、X線源と、X線源から照射されたX線を検出する検出器と、X線源と検出器との間に配置され、X線源から照射されるX線により自己像を形成するための第1格子と、第1格子の自己像と干渉させるための第2格子と、を含む複数の格子と、X線源と検出器との間に配置される被写体、および、複数の格子の少なくともいずれか一方をX線の光軸を中心に回転させ、被写体と複数の格子との相対角度を変化させる第1回転機構と、第1回転機構による回転により変化された複数の相対角度の各々において、検出器によって検出された信号に基づいて位相コントラスト画像を生成する画像処理部と、を備え、複数の相対角度の各々の情報を、対応する位相コントラスト画像に関連付けて記憶するように構成されている。
この発明の一の局面におけるX線位相イメージング装置では、上記のように、複数の相対角度の各々の情報を、対応する位相コントラスト画像に関連付けて記憶することによって、複数の位相コントラスト画像を画像解析する際に、各々の位相コントラスト画像に関連付けて記憶された相対角度を用いて画像解析を行うことができる。その結果、各々の位相コントラスト画像から対応する相対角度を算出する場合と異なり、画像解析において誤った相対角度が用いられることを抑制することができる。これにより、被写体と複数の格子との間の誤った相対角度に基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、複数の位相コントラスト画像の各々に関連付けて記憶された相対角度の情報を読み出すとともに、読み出された複数の相対角度の情報に基づいて画像解析を行うように構成されている。このように構成すれば、複数の位相コントラスト画像の各々から対応する被写体と複数の格子との相対角度を算出する場合と異なり、相対角度を読み出すだけで画像解析を行うことができるので、X線位相イメージング装置の制御負荷を軽減することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、複数の位相コントラスト画像の各々のファイル名に、対応する相対角度の情報を含ませることにより、位相コントラスト画像と相対角度の情報とを関連付けて記憶するように構成されている。このように構成すれば、対応する位相コントラスト画像と相対角度の情報とを共通のフォルダに保存することにより互いを関連付ける場合と異なり、フォルダを新しく生成する必要がない。また、対応する位相コントラスト画像の相対角度の情報をファイルに記載することにより互いを関連付ける場合と異なり、ファイルを新しく生成する必要がない。これらの結果、新しくフォルダおよびファイルを生成することなくファイル名に相対角度の情報を持たせるだけで位相コントラスト画像と相対角度の情報との関連付けを行うことができるので、メモリ等において使用されるデータ容量が増大するのを抑制することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、複数の位相コントラスト画像の各々に含まれる補足データ部に、対応する相対角度の情報を含ませることにより、位相コントラスト画像と相対角度の情報とを関連付けて記憶するように構成されている。このように構成すれば、対応する位相コントラスト画像と相対角度の情報とを共通のフォルダに保存することにより互いを関連付ける場合と異なり、フォルダを新しく生成する必要がない。また、対応する位相コントラスト画像の相対角度の情報をファイルに記載することにより互いを関連付ける場合と異なり、ファイルを新しく生成する必要がない。これらの結果、新しくフォルダおよびファイルを生成することなく補足データ部に相対角度の情報を持たせるだけで位相コントラスト画像と相対角度の情報との関連付けを行うことができるので、メモリ等において使用されるデータ容量が増大するのを抑制することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、複数の位相コントラスト画像の各々を互いに別個のフォルダに格納するとともに、複数の位相コントラスト画像の各々と、対応する相対角度の情報とを共通のフォルダに格納することにより、位相コントラスト画像と相対角度の情報とを関連付けて記憶するように構成されている。このように構成すれば、位相コントラスト画像自体に相対角度の情報を持たせる場合に比べて、位相コントラスト画像のデータ容量を低減することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、複数の位相コントラスト画像の各々と、対応する相対角度の情報との関係を示す、複数の位相コントラスト画像とは異なる管理ファイルを作成することにより、位相コントラスト画像と相対角度の情報とを関連付けて記憶するように構成されている。このように構成すれば、位相コントラスト画像自体に相対角度の情報を持たせる場合に比べて、位相コントラスト画像のデータ容量を低減することができる。さらに、管理ファイルに基づいて位相コントラスト画像と相対角度の情報との関連付けを行うことによって、管理ファイルにより一括で複数の暗視野像の各々と、対応する相対角度の情報との関連付けを行うことができる。その結果、各々の暗視野像毎に、対応する相対角度の情報を記載するためのファイルを生成する必要がないので、メモリ等に格納するファイルの数を低減することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、被写体と複数の格子との相対角度を変化させるように第1回転機構を制御する制御部をさらに備え、画像処理部により生成された複数の位相コントラスト画像と、制御部から送信された対応する相対角度の情報とを関連付けて記憶するように構成されている。このように構成すれば、画像処理部により位相コントラスト画像の生成および相対角度の情報の取得の両方を行う場合に比べて、画像処理部の制御負荷を軽減することができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、被写体と、X線源と検出器と複数の格子とによって構成される撮像系とを光軸に直交する方向を中心に相対的に回転させる第2回転機構をさらに備え、画像処理部は、第1回転機構により被写体を回転させることによって変化された複数の相対角度の各々において、第2回転機構により被写体と撮像系とが相対的に回転された場合に、複数の相対角度の各々において3次元位相コントラスト画像を生成するように構成され、複数の相対角度の各々の情報を、対応する3次元位相コントラスト画像に関連付けて記憶するように構成されている。このように構成すれば、第2回転機構を用いて3次元位相コントラスト画像を生成する場合において、被写体と複数の格子との間の誤った相対角度に基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することができる。
この場合、好ましくは、複数の3次元位相コントラスト画像の各々に関連付けて記憶された相対角度の情報を読み出すとともに、読み出された相対角度の情報に基づいて、対応する3次元位相コントラスト画像を回転させて、複数の3次元位相コントラスト画像の向きを揃えるように構成されている。ここで、3次元位相コントラスト画像を生成する場合においては、複数の格子を回転させるよりも被写体を回転させることにより相対角度を変化させるのが好ましいことが一般的に知られている。被写体を回転させた場合、複数の相対角度の各々において生成された3次元位相コントラスト画像は、相対角度の差の分、向きがずれるので、画像解析を適切に行うことができない。そこで、読み出された相対角度の情報に基づいて、対応する3次元位相コントラスト画像を光軸中心に回転させて、複数の3次元位相コントラスト画像の向きを揃えることによって、回転後の複数の3次元位相コントラスト画像を用いて適切に画像解析を行うことができる。
上記一の局面におけるX線位相イメージング装置において、好ましくは、画像処理部は、複数の相対角度の各々において暗視野像を生成するように構成され、画像処理部により生成された複数の暗視野像の各々と、対応する相対角度の情報とに基づいて画像解析を行うことにより、被写体の配向角度と散乱強度との関係を算出するように構成されている。このように構成すれば、複数の暗視野像を画像解析する際に、各々の暗視野像に関連付けて記憶された相対角度を用いて画像解析を行うことができる。その結果、各々の暗視野像から対応する相対角度を算出する場合と異なり、画像解析において誤った相対角度が用いられることを抑制することができる。これにより、被写体と複数の格子との間の誤った相対角度に基づいて暗視野像の画像解析が行われることに起因して、算出された被写体の配向角度と散乱強度との関係が不適切になるのを抑制することができる。
本発明によれば、上記のように、被写体と複数の格子との間の誤った相対角度に基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することができる。
以下、本発明を具体化した実施形態を図面に基づいて説明する。
[第1実施形態]
図1~図8を参照して、第1実施形態によるX線位相イメージング装置100の構成について説明する。
図1~図8を参照して、第1実施形態によるX線位相イメージング装置100の構成について説明する。
(X線位相イメージング装置の構成)
図1に示すように、X線位相イメージング装置100は、被写体Tを通過したX線の拡散(散乱)を利用して、被写体Tの暗視野像を生成する。具体的には、X線位相イメージング装置100は、タルボ(Talbot)効果を利用して、被写体Tの暗視野像を生成する。X線位相イメージング装置100は、たとえば、非破壊検査用途では、物体の内部の画像化に用いることが可能である。なお、暗視野像は、請求の範囲の「位相コントラスト画像」の一例である。
図1に示すように、X線位相イメージング装置100は、被写体Tを通過したX線の拡散(散乱)を利用して、被写体Tの暗視野像を生成する。具体的には、X線位相イメージング装置100は、タルボ(Talbot)効果を利用して、被写体Tの暗視野像を生成する。X線位相イメージング装置100は、たとえば、非破壊検査用途では、物体の内部の画像化に用いることが可能である。なお、暗視野像は、請求の範囲の「位相コントラスト画像」の一例である。
被写体Tは、内部に繊維束11(図2参照)を含む。被写体Tは、たとえば、繊維束11として炭素繊維が用いられる炭素繊維強化プラスチック(CFRP)である。なお、繊維束とは、繊維が多数集まって束状になったものである。第1実施形態では、繊維束11は、多数の繊維によって板状に形成されたものである。なお、被写体Tは、後述するX線源1と検出器5との間に配置されている。詳細には、被写体Tは、後述する第1格子2と第3格子4との間に配置されている。
図1は、X線位相イメージング装置100をX方向から見た図である。図1に示すように、X線位相イメージング装置100は、X線源1と、第1格子2と、第2格子3と、第3格子4と、検出器5と、PC6と、制御部7と、格子回転機構8と、格子移動機構9とを備えている。なお、本明細書において、X線源1から第1格子2に向かう方向をZ2方向、その逆向きの方向をZ1方向とする。また、Z方向と直交する面内の左右方向をX方向とし、紙面の奥に向かう方向をX2方向、紙面の手前側に向かう方向をX1方向とする。また、Z方向と直交する面内の上下方向をY方向とし、上方向をY1方向、下方向をY2方向とする。また、格子回転機構8は、請求の範囲の「第1回転機構」の一例である。
X線源1は、高電圧が印加されることにより、X線を発生させるとともに、発生されたX線をZ2方向に向けて照射するように構成されている。
第1格子2は、Y方向に所定の周期(ピッチ)d1で配列される複数のスリット2a、および、X線位相変化部2bを有している。各スリット2aおよびX線位相変化部2bはそれぞれ、直線状に延びるように形成されている。また、各スリット2aおよびX線位相変化部2bはそれぞれ、平行に延びるように形成されている。第1格子2は、いわゆる位相格子である。
第1格子2は、X線源1と、第2格子3との間に配置されており、X線源1からX線が照射される。第1格子2は、タルボ効果により、第1格子2の自己像(図1の破線で表示)を形成するために設けられている。なお、可干渉性を有するX線が、スリットが形成された格子を通過すると、格子から所定の距離(タルボ距離)離れた位置に、格子の像(自己像)が形成される。これをタルボ効果という。
第2格子3は、Y方向に所定の周期(ピッチ)d2で配列される複数のX線透過部3aおよびX線吸収部3bを有する。X線吸収部3bは、X線位相変化部2bが延びる方向に沿って延びている。各X線透過部3aおよびX線吸収部3bはそれぞれ、直線状に延びるように形成されている。また、各X線透過部3aおよびX線吸収部3bはそれぞれ、平行に延びるように形成されている。第2格子3は、いわゆる、吸収格子である。第1格子2および第2格子3はそれぞれ異なる役割を持つ格子であるが、スリット2aおよびX線透過部3aはそれぞれX線を透過させる。また、X線吸収部3bはX線を遮蔽する役割を担っており、X線位相変化部2bはスリット2aとの屈折率の違いによってX線の位相を変化させる。
第2格子3は、第1格子2と検出器5との間に配置されており、第1格子2を通過したX線が照射される。また、第2格子3は、第1格子2からタルボ距離離れた位置に配置される。第2格子3は、第1格子2の自己像と干渉して、検出器5の検出表面上にモアレ縞(図示せず)を形成する。
第3格子4は、Y方向に所定の周期(ピッチ)d3で配列される複数のX線透過部4aおよびX線吸収部4bを有する。各X線透過部4aおよびX線吸収部4bはそれぞれ、直線状に延びるように形成されている。また、各X線透過部4aおよびX線吸収部4bはそれぞれ、平行に延びるように形成されている。第3格子4は、いわゆる、マルチスリットである。
第3格子4は、X線源1と第1格子2との間に配置されている。第3格子4は、各X線透過部4aを通過したX線を線光源とすることにより、X線源1からのX線を多点光源化するように構成されている。3枚の格子(第1格子2、第2格子3、および、第3格子4)のピッチと格子間の距離とが一定の条件を満たすことにより、X線源1から照射されるX線の可干渉性を高めることが可能である。これにより、X線源1の管球の焦点サイズが大きくても干渉強度を保持できる。
検出器5は、X線を検出するとともに、検出されたX線を電気信号に変換し、変換された電気信号を画像信号として読み取るように構成されている。検出器5は、たとえば、FPD(Flat Panel Detector)である。検出器5は、複数の変換素子(図示せず)と複数の変換素子上に配置された画素電極(図示せず)とにより構成されている。複数の変換素子および画素電極は、所定の周期(画素ピッチ)で、X方向およびY方向にアレイ状に配列されている。また、検出器5は、取得した画像信号を、PC6に出力するように構成されている。
PC6は、画像処理部60と、メモリ61とを含む。画像処理部60は、画像生成部60aと、後処理部60bとを有している。
画像処理部60(画像生成部60a)は、検出器5から出力された画像信号に基づいて、吸収像(図示せず)を生成するように構成されている。また、画像処理部60(画像生成部60a)は、検出器5により検出されたX線の強度分布に基づいて、X線の散乱に起因する暗視野像(図3参照)を生成するように構成されている。ここで、吸収像とは、被写体TによるX線の吸収の差によって生じるコントラストを画像化したものである。また、暗視野像とは、被写体Tの内部にある微細構造によるX線の微少角散乱によって生じるコントラストを画像化したものである。また、画像処理部60は、たとえば、CPU(Central Processing Unit)およびGPU(Graphics Processing Unit)や画像処理用に構成されたFPGA(Field-Programmable Gate Array)などのプロセッサを含む。
制御部7は、被写体Tと、第1格子2、第2格子3、および、第3格子4(以下では、3枚の格子をまとめて複数の格子と記載する)との相対角度θを、格子回転機構8を制御して変化させる。具体的には、X線の光軸αを中心に、複数の格子を格子回転機構8により回転させることにより相対角度θを変化させる。複数の格子を同時に回転させても、順に回転させてもよい。
詳細には、格子回転機構8は、第1格子2を回転させる回転機構8aと、第2格子3を回転させる回転機構8bと、第3格子4を回転させる回転機構8cとを含む。なお、X線位相イメージング装置100は、第1格子2と回転機構8aと格子移動機構9とを含む回転格子ユニット20を備える。また、X線位相イメージング装置100は、第2格子3と回転機構8bとを含む回転格子ユニット30を備える。また、X線位相イメージング装置100は、第3格子4と回転機構8cとを含む回転格子ユニット40を備える。
また、制御部7は、格子移動機構9を制御して、第1格子2を格子面内において格子方向と直交する方向にステップ移動させるように構成されている。X線位相イメージング装置100では、第1格子2を一定周期間隔に走査することにより得られた複数のモアレ縞(画像)から再構成画像を取得する手法(縞走査法)が用いられている。また、制御部7は、たとえば、CPU(Central Processing Unit)などのプロセッサを含む。なお、格子方向とは、格子の格子パターンが延びる方向である。また、格子パターンとは、各格子のスリット2a、X線位相変化部2b、X線透過部3a、および、X線吸収部3bなどのことである。
格子移動機構9は、制御部7からの信号に基づいて、第1格子2を格子面内(XY面内)において格子方向と直交する方向(図1ではY方向)にステップ移動させるように構成されている。具体的には、格子移動機構9は、第1格子2の周期d1をn分割し、d1/nずつ第1格子2をステップ移動させる。格子移動機構9は、少なくとも第1格子2の1周期d1分、第1格子2をステップ移動させるように構成されている。なお、nは正の整数であり、たとえば、9などである。また、格子移動機構9は、たとえば、ステッピングモータやピエゾアクチュエータなどを含む。なお、格子回転機構8により第1格子2を回転させる場合は、格子方向と直交する方向と共にステップ移動する方向も変わる。
図3に示すように、画像処理部60(画像生成部60a)は、格子回転機構8による回転により変化された複数の相対角度θの各々において、検出器5によって検出された信号に基づいて暗視野像を生成する。具体的には、画像処理部60(画像生成部60a)は、相対角度θが0度、45度、90度、および、135度の各々の状態において、暗視野像を生成する。なお、相対角度θとは、複数の格子の向き(たとえばX線位相変化部2bが延びる方向)と、X方向とのなす角度である。また、図3に示すように、各々の暗視野像において、複数の格子の向きと平行に延びる繊維束11の散乱強度が強く(より濃い黒色に)表示されている。
ここで、第1実施形態では、X線位相イメージング装置100は、複数の相対角度θの各々の情報を、対応する暗視野像に関連付けて記憶するように構成されている。具体的には、対応する相対角度θの情報と関連付けられた各々の暗視野像のデータは、PC6のメモリ61に格納(記憶)される。
詳細には、X線位相イメージング装置100は、画像処理部60(画像生成部60a)により生成された複数の暗視野像と、制御部7から送信された対応する相対角度θの情報とを関連付けて記憶するように構成されている。たとえば、相対角度θ=0°の状態で暗視野像が生成された場合、制御部7から画像処理部60(画像生成部60a)に相対角度θ=0°という情報が送信される。そして、画像処理部60(画像生成部60a)は、生成した暗視野像と、制御部7から送信された相対角度θ=0°という情報とを関連付けて、メモリ61に格納(記憶)する。
また、第1実施形態では、図4に示すように、X線位相イメージング装置100は、複数の暗視野像の各々のファイル名に、対応する相対角度θの情報を含ませることにより、暗視野像と相対角度の情報とを関連付けて記憶するように構成されている。具体的には、θ=0°、45°、90°、および、135°の各々の暗視野像のデータは、それぞれ、test_0deg、test_45deg、test_90deg、および、test_135degというファイル名でメモリ61に格納(記憶)される。なお、ファイル名(上記のtest_0degなど)は、X線位相イメージング装置100の使用者が任意に指定する。
また、X線位相イメージング装置100は、複数の暗視野像の各々に関連付けて記憶された相対角度θの情報を読み出すとともに、読み出された複数の相対角度θの情報に基づいて画像解析を行うように構成されている。具体的には、画像処理部60(後処理部60b)は、メモリ61に格納された複数の暗視野像のデータのファイル名を検索して読み出すことにより、各々の暗視野像のデータに関連付けられた相対角度θの情報を取得する。
また、X線位相イメージング装置100は、画像処理部60(画像生成部60a)により生成された複数の暗視野像の各々と、対応する相対角度θの情報とに基づいて画像解析を行うことにより、被写体Tの配向角度と散乱強度との関係を算出するように構成されている。
具体的には、図5に示すように、画像処理部60(後処理部60b)は、暗視野像の各画素において、相対角度θに対する画素値の値(散乱強度に基づく値)をプロットする。そして、画像処理部60(後処理部60b)は、プロットして作成された波形を、A+Bcos[2(θ―φ)+π]という式により近似する。Aとは、波形におけるオフセット情報である。また、Bとは、波形における振幅の情報である。A、B、および、θの値に基づいて、被写体Tの繊維束11の配向角度であるφが算出される。算出された配向角度であるφを画像化した画像を図6(A)に示す。また、B/Aは、X線の散乱の異方性度を示す。X線の散乱の異方性度であるB/Aを画像化した画像を図6(B)に示す。そして、B/Aの数値が高い部分のみの画素位置を配向角度の情報であるφにより色付けして表示することによって、繊維束11の方向を示す情報(図7参照)が得られる。図7では、繊維束11の配向角度に応じて繊維束11の画像の色が変化するように表示されている。これにより、繊維束11の配向を視覚的、定量的に解析することが可能である。
次に、図8を参照して、画像処理部60による制御フローを説明する。
まず、ステップS1において、複数の相対角度θの各々において暗視野像を生成する。次に、ステップS2において、複数の暗視野像のファイル名に、対応する相対角度θの情報を含ませてメモリ61に格納(記憶)する。そして、ステップS3において、メモリ61に格納(記憶)された複数の暗視野像の各々のファイル名から相対角度θの情報を読み出すとともに、暗視野像と読み出された相対角度θとに基づいて画像解析を行う。
(第1実施形態の効果)
第1実施形態では、以下のような効果を得ることができる。
第1実施形態では、以下のような効果を得ることができる。
第1実施形態では、上記のように、複数の相対角度θの各々の情報を、対応する暗視野像に関連付けて記憶するように、X線位相イメージング装置100を構成する。これにより、複数の暗視野像を画像解析する際に、各々の暗視野像に関連付けて記憶された相対角度θを用いて画像解析を行うことができる。その結果、各々の暗視野像から対応する相対角度θを算出する場合と異なり、画像処理において誤った相対角度が用いられることを抑制することができる。これにより、被写体Tと複数の格子との間の誤った相対角度θに基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することができる。
また、第1実施形態では、上記のように、複数の暗視野像の各々に関連付けて記憶された相対角度θの情報を読み出すとともに、読み出された複数の相対角度θの情報に基づいて画像解析を行うように、X線位相イメージング装置100を構成する。これにより、複数の暗視野像の各々から対応する被写体Tと複数の格子との相対角度θを算出する場合と異なり、相対角度θを読み出すだけで画像解析を行うことができるので、X線位相イメージング装置100の制御負荷を軽減することができる。
また、第1実施形態では、上記のように、複数の暗視野像の各々のファイル名に、対応する相対角度θの情報を含ませることにより、暗視野像と相対角度の情報とを関連付けて記憶するように、X線位相イメージング装置100を構成する。これにより、対応する暗視野像と相対角度θの情報とを共通のフォルダに保存することにより互いを関連付ける場合と異なり、フォルダを新しく生成する必要がない。また、対応する暗視野像の相対角度θの情報をファイルに記載することにより互いを関連付ける場合と異なり、ファイルを新しく生成する必要がない。これらの結果、新しくフォルダおよびファイルを生成することなくファイル名に相対角度θの情報を持たせるだけで暗視野像と相対角度θの情報との関連付けを行うことができるので、メモリ等において使用されるデータ容量が増大するのを抑制することができる。
また、第1実施形態では、上記のように、画像処理部60により生成された複数の暗視野像と、制御部7から送信された対応する相対角度θの情報とを関連付けて記憶するように、X線位相イメージング装置100を構成する。これにより、画像処理部60により暗視野像の生成および相対角度θの情報の取得の両方を行う場合に比べて、画像処理部60の制御負荷を軽減することができる。
また、第1実施形態では、上記のように、画像処理部60により生成された複数の暗視野像の各々と、対応する相対角度θの情報とに基づいて画像解析を行うことにより、被写体Tの配向角度と散乱強度との関係を算出するように、X線位相イメージング装置100を構成する。これにより、複数の暗視野像を画像解析する際に、各々の暗視野像に関連付けて記憶された相対角度θを用いて画像解析を行うことができる。その結果、各々の暗視野像から対応する相対角度θを算出する場合と異なり、画像解析において誤った相対角度θが用いられることを抑制することができる。これにより、被写体Tと複数の格子との間の誤った相対角度θに基づいて暗視野像の画像解析が行われることに起因して、算出された被写体Tの配向角度と散乱強度との関係が不適切になるのを抑制することができる。
[第2実施形態]
次に、図9~図12を参照して、第2実施形態によるX線位相イメージング装置200の構成について説明する。この第2実施形態のX線位相イメージング装置200は、上記第1実施形態のX線位相イメージング装置100とは異なり、複数の相対角度の各々において3次元暗視野像を生成するように構成されている。なお、上記第1実施形態と同様の構成は、第1実施形態と同じ符号を付して図示するとともに説明を省略する。
次に、図9~図12を参照して、第2実施形態によるX線位相イメージング装置200の構成について説明する。この第2実施形態のX線位相イメージング装置200は、上記第1実施形態のX線位相イメージング装置100とは異なり、複数の相対角度の各々において3次元暗視野像を生成するように構成されている。なお、上記第1実施形態と同様の構成は、第1実施形態と同じ符号を付して図示するとともに説明を省略する。
(X線位相イメージング装置の構成)
図9に示すように、X線位相イメージング装置200は、PC16と、制御部17と、被写体回転機構18と、CT回転機構28とを備えている。PC16は、画像処理部160と、メモリ161とを含む。画像処理部160は、画像生成部160aと、後処理部160bとを有している。また、CT回転機構28は、たとえば、モータなどによって駆動される回転ステージ28aを含む。なお、被写体回転機構18およびCT回転機構28は、それぞれ、請求の範囲の「第1回転機構」および「第2回転機構」の一例である。
図9に示すように、X線位相イメージング装置200は、PC16と、制御部17と、被写体回転機構18と、CT回転機構28とを備えている。PC16は、画像処理部160と、メモリ161とを含む。画像処理部160は、画像生成部160aと、後処理部160bとを有している。また、CT回転機構28は、たとえば、モータなどによって駆動される回転ステージ28aを含む。なお、被写体回転機構18およびCT回転機構28は、それぞれ、請求の範囲の「第1回転機構」および「第2回転機構」の一例である。
被写体回転機構18は、X線の光軸αを中心に、被写体Tを回転させることにより、被写体Tと複数の格子との相対角度θを変化させる。
CT回転機構28は、制御部17からの信号に基づいて、被写体Tと、X線源1、検出器5、および、複数の格子によって構成される撮像系300とを相対的に回転させるように構成されている。具体的には、CT回転機構28は、被写体Tを光軸αと直交する軸線β周りに回転させることにより、撮像系300に対して被写体Tを相対的に回転させるように構成されている。
なお、X線位相イメージング装置200は、第1格子2と格子移動機構9を含む格子ユニット120を備える。また、X線位相イメージング装置200は、第2格子3を含む格子ユニット130を備える。また、X線位相イメージング装置200は、第3格子4を含む格子ユニット140を備える。また、X線位相イメージング装置200は、被写体回転機構18とCT回転機構28とを含む被写体保持ユニット150を備える。
また、画像処理部16(画像生成部160a)は、CT回転機構28を回転させながら(複数の回転角度の各々において)撮像された複数の吸収像および複数の暗視野像をそれぞれ再構成することにより、3次元吸収像(図示せず)および3次元暗視野像(図示せず)を生成する。なお、3次元暗視野像は、請求の範囲の「3次元位相コントラスト画像」の一例である。
具体的には、図10に示すように、画像処理部16(画像生成部160a)は、被写体回転機構18により被写体Tを回転させることによって変化された複数の相対角度θの各々において、CT回転機構28により被写体Tと撮像系300とが相対的に回転された場合に、複数の相対角度θの各々において3次元暗視野像を生成するように構成されている。詳細には、画像処理部16(画像生成部160a)は、相対角度θ=45°、0°、-45°の各々の場合において、3次元暗視野像を生成する。なお、図10では、被写体Tの画像を簡略化して図示している。
ここで、第2実施形態では、X線位相イメージング装置200は、複数の相対角度θの各々の情報を、対応する3次元暗視野像に関連付けて記憶するように構成されている。具体的には、対応する相対角度θの情報と関連付けられた各々の3次元暗視野像のデータは、PC16のメモリ161に格納(記憶)される。
また、図11に示すように、X線位相イメージング装置200は、複数の3次元暗視野像の各々に関連付けて記憶された相対角度θの情報を読み出すとともに、読み出された相対角度θの情報に基づいて、対応する3次元暗視野像を、光軸αを中心に回転させて、複数の3次元暗視野像の向きを揃えるように構成されている。たとえば、画像処理部160(後処理部160b、図9参照)は、θ=0°の3次元暗視野像と被写体Tの像の向きを揃えるように、θ=45°の3次元暗視野像を-45°回転させる。また、画像処理部160(後処理部160b)は、θ=0°の3次元暗視野像と被写体Tの像の向きを揃えるように、θ=-45°の3次元暗視野像を45°回転させる。なお、θ=0°の3次元暗視野像を回転させてもよい。
複数の3次元暗視野像の向きが揃った状態で、画像処理部160(後処理部160b)は、複数の3次元暗視野像の各々のスライス画像を用いて画像解析を行うように構成されている。スライス画像とは、3次元暗視野像における断面画像を意味する。スライス画像は、3次元暗視野像において複数枚取得される。これにより、複数のスライス画像毎に、繊維束11の配向を視覚的、定量的に解析することが可能である。
次に、図12を参照して、画像処理部160による制御フローを説明する。
まず、ステップS11において、複数の相対角度θの各々において3次元暗視野像を生成する。次に、ステップS12において、複数の3次元暗視野像のファイル名に、対応する相対角度θの情報を含ませてメモリ161に格納(記憶)する。次に、ステップS13において、複数の3次元暗視野像の各々のファイル名から相対角度θの情報を読み出すとともに、3次元暗視野像の少なくとも2つを回転させて互いの向きを揃える。そして、ステップS14において、3次元暗視野像のスライス画像毎に、関連付けられた相対角度θの情報に基づいて画像解析を行う。
第2実施形態のその他の構成は、上記第1実施形態と同様である。
(第2実施形態の効果)
第2実施形態では、以下のような効果を得ることができる。
第2実施形態では、以下のような効果を得ることができる。
第2実施形態では、上記のように、複数の相対角度θの各々の情報を、対応する3次元暗視野像に関連付けて記憶するように、X線位相イメージング装置200を構成する。これにより、CT回転機構28を用いて3次元暗視野像を生成する場合において、被写体Tと複数の格子との間の誤った相対角度θに基づいて画像解析が行われることに起因して、解析結果の情報が不適切になるのを抑制することができる。
また、第2実施形態では、上記のように、複数の3次元暗視野像の各々に関連付けて記憶された相対角度θの情報を読み出すとともに、読み出された相対角度θの情報に基づいて、対応する3次元暗視野像を回転させて、複数の3次元暗視野像の向きを揃えるように、X線位相イメージング装置200を構成する。ここで、3次元暗視野像を生成する場合においては、複数の格子を回転させるよりも被写体Tを回転させることにより相対角度θを変化させるのが好ましいことが一般的に知られている。被写体Tを回転させた場合、複数の相対角度θの各々において生成された3次元位暗視野像は、相対角度θの差の分、向きがずれるので、画像解析を適切に行うことができない。そこで、読み出された相対角度θの情報に基づいて、対応する3次元暗視野像を回転させて、複数の3次元暗視野像の向きを揃えることによって、回転後の複数の3次元暗視野像を用いて適切に画像処理を行うことができる。
なお、第2実施形態のその他の効果は、上記第1実施形態と同様である。
(変形例)
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく、請求の範囲によって示され、さらに請求の範囲と均等の意味および範囲内でのすべての変更(変形例)が含まれる。
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく、請求の範囲によって示され、さらに請求の範囲と均等の意味および範囲内でのすべての変更(変形例)が含まれる。
たとえば、上記第1および第2実施形態では、暗視野像(位相コントラスト画像)のファイル名に相対角度θの情報を含ませる例を示したが、本発明はこれに限られない。ファイル名以外に相対角度θの情報を含ませてもよい。
具体的には、図13に示すように、X線位相イメージング装置400は、複数の暗視野像の各々に含まれるヘッダに、対応する相対角度θの情報を含ませることにより、暗視野像と相対角度θの情報とを関連付けて記憶するように構成されている。たとえば、θ=0°、45°、90°、および、135°の暗視野像のヘッダには、それぞれ、0deg、45deg、90deg、135degと記載されている。ヘッダに相対角度θの情報が含まれる暗視野像は、メモリ261に格納(記憶)される。ヘッダは実際には暗視野像には表示されないが、図13では説明のために表示されているように図示する。なお、ヘッダは、請求の範囲の「補足データ部」の一例である。
これにより、対応する暗視野像と相対角度θの情報とを共通のフォルダに保存することにより互いを関連付ける場合と異なり、フォルダを新しく生成する必要がない。また、対応する暗視野像の相対角度θの情報をファイルに記載することにより互いを関連付ける場合と異なり、ファイルを新しく生成する必要がない。これらの結果、新しくフォルダおよびファイルを生成することなくヘッダに相対角度θの情報を持たせるだけで暗視野像と相対角度θの情報との関連付けを行うことができるので、メモリ261において使用されるデータ容量が増大するのを抑制することができる。
また、図14に示すように、X線位相イメージング装置500は、複数の暗視野像の各々を互いに別個のフォルダに格納するとともに、複数の暗視野像の各々と、対応する相対角度θの情報とを共通のフォルダに格納することにより、暗視野像と相対角度θの情報とを関連付けて記憶するように構成されている。具体的には、メモリ361には、フォルダ361a、フォルダ361b、フォルダ361c、フォルダ361dが作成されている。フォルダ361aには、θ=0°の暗視野像と、フォルダ361aに格納されている暗視野像はθ=0°で生成された暗視野像であることを示すファイル500a(たとえばテキストファイルまたはバイナリファイル)が格納されている。フォルダ361bには、θ=45°の暗視野像とファイル500bとが格納され、フォルダ361cには、θ=90°の暗視野像とファイル500cとが格納され、フォルダ361dには、θ=135°の暗視野像とファイル500bとが格納されている。ファイル500b~500dについては、ファイル500aと同様であるので、詳細な説明は省略する。
これにより、暗視野像自体に相対角度の情報を持たせる場合に比べて、暗視野像のデータ容量を低減することができる。
また、図15に示すように、X線位相イメージング装置600は、複数の暗視野像の各々と、対応する相対角度θの情報との関係を示す、複数の暗視野像とは異なる管理ファイル600aを作成することにより、暗視野像と相対角度θの情報とを関連付けて記憶するように構成されている。具体的には、メモリ461には、複数の暗視野像の各々と、相対角度θとの対応関係が記載されている管理ファイル600a(たとえばテキストファイルまたはバイナリファイル)が格納されている。
また、複数の暗視野像のデータが別個のフォルダに格納されており、管理ファイル600aに、フォルダ名と相対角度θとの対応関係が記載されていてもよい。
これにより、暗視野像自体に相対角度θの情報を持たせる場合に比べて、暗視野像のデータ容量を低減することができる。さらに、管理ファイル600aに基づいて暗視野像と相対角度θの情報との関連付けを行うことによって、管理ファイル600aにより一括で複数の暗視野像の各々と、対応する相対角度θの情報との関連付けを行うことができる。その結果、各々の暗視野像毎に、対応する相対角度θの情報を記載するためのファイルを生成する必要がないので、メモリ461に格納するファイルの数を低減することができる。
なお、上記第1および第2実施形態の構成と、図13~15の構成とを組み合わせてもよい。
また、上記第1および第2実施形態では、複数の相対角度θの情報を、対応する暗視野像に関連付けて記憶する例を示したが、本発明はこれに限られない。たとえば、複数の相対角度θの情報を、吸収像または位相微分像に関連付けて記憶してもよい。
また、上記第1実施形態では、複数の格子を回転させて相対角度θを変化させる例を示したが、本発明はこれに限られない。たとえば、被写体Tを回転させて相対角度θを変化させてもよい。また、複数の格子および被写体Tの両方を回転させてもよい。
また、上記第2実施形態では、被写体Tを回転させて相対角度θを変化させる例を示したが、本発明はこれに限られない。たとえば、第1実施形態と同様に、複数の格子を回転させて相対角度θを変化させてもよい。
また、上記第1(第2)実施形態では、相対角度θを4つ(3つ)の角度に変化させる例を示したが、本発明はこれに限られない。3つ以上であれば、相対角度θの角度の数は限定されない。
また、上記第1(第2)実施形態では、被写体Tは、第1格子2と第3格子4との間に配置されている例を示したが、第1格子2と第2格子3との間に配置されていてもよい。
また、上記第2実施形態では、被写体TをCT回転機構28(第2回転機構)により回転させてCT撮影を行う例を示したが、本発明はこれに限られない。たとえば、撮像系300を回転させてCT撮影を行ってもよい。
また、上記第1および第2実施形態では、第3格子4が設けられている例を示したが、本発明はこれに限られない。第3格子4が設けられていなくてもよい。
また、上記第1および第2実施形態では、縞走査法によって暗視野像を生成する例を示したが、本発明はこれに限られない。たとえば、第1格子2、第2格子3、または、第3格子4のうちのいずれかを、X線の光軸αが延びる方向に直交する平面上において回転させる手法(いわゆる、モアレ1枚撮り手法)によって、暗視野像を生成してもよい。
また、上記第1および第2実施形態では、第1格子2が位相格子である例を示したが、本発明はこれに限られない。たとえば、第1格子2は吸収格子であってもよい。
また、上記第1および第2実施形態では、第1格子2を格子面内においてステップ移動させる例を示したが、本発明はこれに限られない。複数の格子のうち、いずれの格子をステップ移動させてもよい。
また、上記第1および第2実施形態では、被写体Tとして、炭素繊維強化プラスチック(CFRP)を撮像する例を示したが、本発明はこれに限られない。たとえば、ガラス繊維強化プラスチック(GFRP)などを被写体Tとして用いてもよい。
また、上記第1および第2実施形態では、被写体Tと複数の格子との相対角度θを暗視野像のデータのファイル名に含ませる例を示したが、本発明はこれに限られない。たとえば、被写体Tおよび複数の格子の各々の回転角度の両方をファイル名に含ませてもよい。また、被写体Tおよび複数の格子のいずれか一方だけを回転させる場合は、回転させる方の回転角度のみをファイル名に含ませておいてもよい。
また、上記第1および第2実施形態では、画像処理部60(160)において、対応する暗視野像と相対角度θとの関連付けを行う例を示したが、本発明はこれに限られない。たとえば、制御部7(17)において、対応する暗視野像と相対角度θとの関連付けを行ってもよい。
また、上記第1および第2実施形態では、説明の便宜上、画像処理部60(160)の処理を「フロー駆動型」のフローチャートを用いて説明したが、本発明はこれに限られない。画像処理部60(160)の処理をイベント単位で実行する「イベント駆動型」により行ってもよい。この場合、完全なイベント駆動型で行ってもよいし、イベント駆動およびフロー駆動を組み合わせて行ってもよい。
1 X線源
2 第1格子
3 第2格子
5 検出器
7、17 制御部
8 格子回転機構(第1回転機構)
18 被写体回転機構(第1回転機構)
28 CT回転機構(第2回転機構)
60、160 画像処理部
100、200、400、500、600 X線位相イメージング装置
300 撮像系
361a、361b、361c、361d フォルダ
600a 管理ファイル
T 被写体
α 光軸
θ 相対角度
2 第1格子
3 第2格子
5 検出器
7、17 制御部
8 格子回転機構(第1回転機構)
18 被写体回転機構(第1回転機構)
28 CT回転機構(第2回転機構)
60、160 画像処理部
100、200、400、500、600 X線位相イメージング装置
300 撮像系
361a、361b、361c、361d フォルダ
600a 管理ファイル
T 被写体
α 光軸
θ 相対角度
Claims (10)
- X線源と、
前記X線源から照射されたX線を検出する検出器と、
前記X線源と前記検出器との間に配置され、前記X線源から照射される前記X線により自己像を形成するための第1格子と、前記第1格子の自己像と干渉させるための第2格子と、を含む複数の格子と、
前記X線源と前記検出器との間に配置される被写体、および、前記複数の格子の少なくともいずれか一方をX線の光軸を中心に回転させ、前記被写体と前記複数の格子との相対角度を変化させる第1回転機構と、
前記第1回転機構による回転により変化された複数の前記相対角度の各々において、前記検出器によって検出された信号に基づいて位相コントラスト画像を生成する画像処理部と、を備え、
前記複数の相対角度の各々の情報を、対応する前記位相コントラスト画像に関連付けて記憶するように構成されている、X線位相イメージング装置。 - 複数の前記位相コントラスト画像の各々に関連付けて記憶された前記相対角度の情報を読み出すとともに、読み出された前記複数の相対角度の情報に基づいて画像解析を行うように構成されている、請求項1に記載のX線位相イメージング装置。
- 複数の前記位相コントラスト画像の各々のファイル名に、対応する前記相対角度の情報を含ませることにより、前記位相コントラスト画像と前記相対角度の情報とを関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。
- 複数の前記位相コントラスト画像の各々に含まれる補足データ部に、対応する前記相対角度の情報を含ませることにより、前記位相コントラスト画像と前記相対角度の情報とを関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。
- 複数の前記位相コントラスト画像の各々を互いに別個のフォルダに格納するとともに、前記複数の位相コントラスト画像の各々と、対応する前記相対角度の情報とを共通の前記フォルダに格納することにより、前記位相コントラスト画像と前記相対角度の情報とを関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。
- 複数の前記位相コントラスト画像の各々と、対応する前記相対角度の情報との関係を示す、前記複数の位相コントラスト画像とは異なる管理ファイルを作成することにより、前記位相コントラスト画像と前記相対角度の情報とを関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。
- 前記被写体と前記複数の格子との前記相対角度を変化させるように前記第1回転機構を制御する制御部をさらに備え、
前記画像処理部により生成された複数の前記位相コントラスト画像と、前記制御部から送信された対応する前記相対角度の情報とを関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。 - 前記被写体と、前記X線源と前記検出器と前記複数の格子とによって構成される撮像系とを前記光軸に直交する方向を中心に相対的に回転させる第2回転機構をさらに備え、
前記画像処理部は、前記第1回転機構により前記被写体を回転させることによって変化された前記複数の相対角度の各々において、前記第2回転機構により前記被写体と前記撮像系とが相対的に回転された場合に、前記複数の相対角度の各々において3次元位相コントラスト画像を生成するように構成され、
前記複数の相対角度の各々の情報を、対応する前記3次元位相コントラスト画像に関連付けて記憶するように構成されている、請求項1に記載のX線位相イメージング装置。 - 複数の前記3次元位相コントラスト画像の各々に関連付けて記憶された前記相対角度の情報を読み出すとともに、読み出された前記相対角度の情報に基づいて、対応する前記3次元位相コントラスト画像を回転させて、前記複数の3次元位相コントラスト画像の向きを揃えるように構成されている、請求項8に記載のX線位相イメージング装置。
- 前記画像処理部は、前記複数の相対角度の各々において暗視野像を生成するように構成され、
前記画像処理部により生成された複数の前記暗視野像の各々と、対応する前記相対角度の情報とに基づいて画像解析を行うことにより、前記被写体の配向角度と散乱強度との関係を算出するように構成されている、請求項1に記載のX線位相イメージング装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018179357 | 2018-09-25 | ||
JP2018-179357 | 2018-09-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020066135A1 true WO2020066135A1 (ja) | 2020-04-02 |
Family
ID=69951840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2019/022384 WO2020066135A1 (ja) | 2018-09-25 | 2019-06-05 | X線位相イメージング装置 |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2020066135A1 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7460577B2 (ja) | 2020-06-03 | 2024-04-02 | 株式会社リガク | X線画像生成装置 |
WO2024070776A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
WO2024070777A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012170618A (ja) * | 2011-02-22 | 2012-09-10 | Konica Minolta Medical & Graphic Inc | X線撮影システム |
WO2012128335A1 (ja) * | 2011-03-23 | 2012-09-27 | コニカミノルタエムジー株式会社 | 医用画像表示システム |
JP2012189446A (ja) * | 2011-03-10 | 2012-10-04 | Canon Inc | 撮像装置及び撮像方法 |
JP2017072399A (ja) * | 2015-10-05 | 2017-04-13 | 株式会社日立ハイテクサイエンス | X線検査装置及びx線検査方法 |
JP2017147532A (ja) * | 2016-02-16 | 2017-08-24 | 大日本印刷株式会社 | 撮像装置、表示装置、撮像装置のプログラム、撮像システム、表示装置のプログラム及び携帯端末 |
WO2017216178A1 (en) * | 2016-06-13 | 2017-12-21 | Technische Universität München | X-ray tensor tomography system |
WO2018096759A1 (ja) * | 2016-11-22 | 2018-05-31 | 株式会社島津製作所 | X線位相イメージング装置 |
-
2019
- 2019-06-05 WO PCT/JP2019/022384 patent/WO2020066135A1/ja active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012170618A (ja) * | 2011-02-22 | 2012-09-10 | Konica Minolta Medical & Graphic Inc | X線撮影システム |
JP2012189446A (ja) * | 2011-03-10 | 2012-10-04 | Canon Inc | 撮像装置及び撮像方法 |
WO2012128335A1 (ja) * | 2011-03-23 | 2012-09-27 | コニカミノルタエムジー株式会社 | 医用画像表示システム |
JP2017072399A (ja) * | 2015-10-05 | 2017-04-13 | 株式会社日立ハイテクサイエンス | X線検査装置及びx線検査方法 |
JP2017147532A (ja) * | 2016-02-16 | 2017-08-24 | 大日本印刷株式会社 | 撮像装置、表示装置、撮像装置のプログラム、撮像システム、表示装置のプログラム及び携帯端末 |
WO2017216178A1 (en) * | 2016-06-13 | 2017-12-21 | Technische Universität München | X-ray tensor tomography system |
WO2018096759A1 (ja) * | 2016-11-22 | 2018-05-31 | 株式会社島津製作所 | X線位相イメージング装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7460577B2 (ja) | 2020-06-03 | 2024-04-02 | 株式会社リガク | X線画像生成装置 |
WO2024070776A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
WO2024070777A1 (ja) * | 2022-09-30 | 2024-04-04 | コニカミノルタ株式会社 | 画像表示装置、画像表示システム、画像表示方法、画像処理装置及びプログラム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6422123B2 (ja) | 放射線画像生成装置 | |
JP5777360B2 (ja) | X線撮像装置 | |
WO2020066135A1 (ja) | X線位相イメージング装置 | |
JP6741080B2 (ja) | X線位相イメージング装置 | |
US20140286475A1 (en) | Interferometer and object information acquisition system | |
JP6058860B2 (ja) | X線撮像装置及びx線撮像方法 | |
JP2019045394A (ja) | X線イメージング装置 | |
JP6897801B2 (ja) | X線イメージング装置 | |
JP7006784B2 (ja) | X線イメージング装置 | |
WO2019220689A1 (ja) | X線イメージング装置 | |
JP7180566B2 (ja) | X線イメージング装置およびx線イメージング方法 | |
EP3395253B1 (en) | X-ray phase imaging apparatus | |
WO2017159229A1 (ja) | 放射線画像生成装置 | |
JP7276210B2 (ja) | X線位相撮像システム、x線画像処理装置、および、x線画像処理方法 | |
JP7226171B2 (ja) | X線位相イメージング装置 | |
WO2019123758A1 (ja) | X線位相差撮像システム | |
WO2018168621A1 (ja) | 放射線画像生成装置 | |
JP6797762B2 (ja) | 放射線画像生成装置及び放射線画像生成方法 | |
US20230000453A1 (en) | Estimation of full-field scattering for dax imaging | |
JP6673189B2 (ja) | X線位相コントラスト撮影装置 | |
WO2020039660A1 (ja) | X線位相イメージング装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19867881 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19867881 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: JP |