WO2020090168A1 - X線位相差撮影システム - Google Patents
X線位相差撮影システム Download PDFInfo
- Publication number
- WO2020090168A1 WO2020090168A1 PCT/JP2019/029232 JP2019029232W WO2020090168A1 WO 2020090168 A1 WO2020090168 A1 WO 2020090168A1 JP 2019029232 W JP2019029232 W JP 2019029232W WO 2020090168 A1 WO2020090168 A1 WO 2020090168A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- ray
- subject
- imaging
- rotation angle
- condition
- Prior art date
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
- 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-contrast imaging system, and more particularly to an X-ray phase-contrast imaging system that captures an image while relatively rotating a subject and an imaging system.
- the X-ray phase contrast imaging system of Non-Patent Document 1 includes an X-ray source, a phase grating, a detector, and an image processing device.
- the X-ray from the X-ray source is scattered by the subject, passes through the phase grating, and is applied to the detector.
- the dark field image is generated by the image processing device based on the interference intensity of the X-ray detected by the detector. Further, a dark field image is generated by photographing at each rotation angle when the subject is rotated by 360 degrees.
- Non-Patent Document 1 in the conventional X-ray phase contrast imaging system as described in Non-Patent Document 1, when taking an image while rotating the subject, the rotation angle is changed.
- the path length through which the X-rays pass may differ depending on the case. For example, when a carbon fiber reinforced plastic (CFRP) having a plate-like shape is rotated and photographed as an object, X-rays are transmitted at a specific range of rotation angle (imaging angle) along the longitudinal direction of the plate-like shape. As the path length increases, the degree of X-ray scattering increases, and the visibility (definition) may decrease to near zero. When the visibility is reduced to near zero, it becomes difficult to accurately extract the X-ray interference intensity, and there is a problem that the image quality of the dark field image deteriorates.
- CFRP carbon fiber reinforced plastic
- the present invention has been made to solve the above problems, and one object of the present invention is to increase the degree of scattering of X-rays based on the rotation angle, resulting in a dark object.
- An object of the present invention is to provide an X-ray phase contrast imaging system capable of suppressing deterioration of the image quality of a field image.
- an X-ray phase contrast imaging system includes an X-ray source, a detector that detects X-rays emitted from the X-ray source, an X-ray source and a detector.
- a plurality of gratings arranged between the two, a subject, and an imaging system composed of an X-ray source, a detector, and a plurality of gratings, and a rotating mechanism that relatively rotates the rotating mechanism.
- the image processing unit that generates at least a dark field image caused by X-ray scattering based on the intensity distribution of the X-ray detected by the detector, and And a control unit that controls switching of imaging conditions based on the degree of X-ray scattering depending on the rotation angle.
- a control unit that performs control for switching imaging conditions based on the magnitude of the degree of X-ray scattering depending on the rotation angle when imaging a subject.
- a control unit that performs control for switching imaging conditions based on the magnitude of the degree of X-ray scattering depending on the rotation angle when imaging a subject.
- the control unit performs imaging under a first imaging condition in a first rotation angle range among the plurality of rotation angles, and the imaging unit is more than the first rotation angle range.
- control is performed to switch to the second imaging condition in which the effective energy of X-rays is relatively larger than the first imaging condition.
- the control unit sets the rotation angle for switching the imaging conditions to a preset value based on the shape of the subject that contributes to the degree of X-ray scattering.
- the rotation angle for switching the imaging conditions to a preset value based on the shape of the subject that contributes to the degree of X-ray scattering.
- the control unit sets the rotation angle for switching the imaging conditions to a preset value based on the shape of the subject that contributes to the degree of X-ray scattering.
- the image processing unit is configured to acquire the feature amount based on the degree of X-ray scattering or the pixel value of the dark field image, and the control unit sets the imaging condition based on the feature amount. It is configured to determine the rotation angle to be switched.
- the degree of X-ray scattering for each detection pixel of the detector at each of the plurality of rotation angles or the pixel value for each pixel of the dark field image at each of the plurality of rotation angles is compared.
- the shooting condition is switched by determining the rotation angle for switching the shooting condition according to the feature amount acquired based on the X-ray scattering degree or the pixel value in a plurality of pixels.
- the amount of information used when determining the rotation angle can be reduced. As a result, it is possible to reduce the calculation load when making the determination of switching the shooting conditions.
- the feature amount is the minimum value or average value of the X-ray scattering degree in the detector, or the minimum value or average value of the pixel values of the region in which the subject appears in the dark field image.
- the control unit is configured to determine the rotation angle for switching the shooting condition by comparing the feature amount with a threshold value. ing. According to this structure, by comparing the feature amount at each rotation angle with the threshold value, it is possible to easily grasp the degree of decrease in the X-ray visibility at each rotation angle. As a result, it is possible to easily determine the rotation angle for switching the shooting conditions.
- the control unit controls the imaging of the subject at a plurality of rotation angles according to the first imaging condition, and after the imaging, the effective energy of the X-ray is relatively larger than that of the first imaging condition. It is configured to switch to the two shooting conditions and perform control to shoot the subject again in each of the rotation angle ranges in which the feature amount is smaller than the first threshold value. According to this structure, it is possible to perform the shooting under the first shooting condition and the shooting under the second shooting condition in two steps. As a result, since it is possible to change the number of times the photographing conditions are changed to one during photographing, the photographing is performed in real time while switching the first photographing condition and the second photographing condition based on the magnitude of the degree of X-ray scattering. Comparing with the configuration described above, it is possible to suppress the control from becoming complicated.
- the control unit performs the shooting under the first shooting condition while the feature amount is smaller than the first threshold value.
- the characteristic amount is greater than the first threshold value while imaging under the second imaging condition.
- the control unit performs at least one of increasing the voltage applied to the X-ray source more than the first imaging condition and changing the X-ray filter provided in the X-ray source. By doing so, the control is performed such that the effective energy of the X-rays emitted toward the detector is relatively increased to switch to the second imaging condition.
- control unit further includes a lattice moving mechanism that translates any one of the plurality of lattices, and the control unit applies the X-ray source to the X-ray source more than the first imaging condition when photographing the subject under the second imaging condition.
- At least one of increasing the current to be applied, increasing the storage time in the detector, and increasing the number of steps for translating any one of the plurality of grids by the grid moving mechanism. Is configured to increase the X-ray dose reaching the detector.
- the S / N ratio signal-to-noise ratio
- the image processing unit when the subject is imaged at a plurality of rotation angles, the degree of scattering or darkness of X-rays discontinuous due to switching of the imaging conditions.
- the correction is performed so as to smoothly connect the pixel values of the visual field image.
- the present invention it is possible to prevent the image quality of the dark-field image of the subject from being deteriorated due to the increase in the degree of X-ray scattering based on the rotation angle as described above.
- FIG. 1 It is a schematic diagram showing a configuration of an X-ray phase contrast imaging system according to the first embodiment.
- FIG. 9 is a graph for explaining the difference in the degree of reduction in visibility due to the difference in effective energy of X-rays. It is a schematic diagram for demonstrating a 1st rotation angle range and a 2nd rotation angle range.
- the schematic diagram (A) showing the graph of the characteristic amount when the X-ray phase contrast imaging system according to the second embodiment images the subject while switching between the first imaging condition and the second imaging condition and the graph of the corrected feature amount. It is a schematic diagram (B) which shows. 9 is a flowchart for explaining a process of acquiring a three-dimensional dark field image by the X-ray phase contrast imaging system according to the second embodiment.
- the X-ray phase contrast imaging system 100 uses the diffusion (scattering) of X-rays that have passed through the subject T to generate a dark field image 22 (see FIG. 2) of the subject T. Specifically, the X-ray phase contrast imaging system 100 uses the Talbot effect to generate at least the dark field image 22 of the subject T.
- the X-ray phase contrast imaging system 100 can be used to image the interior of an object, for example, in non-destructive inspection applications.
- CFRP carbon fiber reinforced plastic
- FIG. 1 is a view of the X-ray phase contrast imaging system 100 viewed from the X direction.
- an X-ray phase contrast imaging system 100 includes an X-ray source 1, a plurality of gratings including a first grating 2, a second grating 3, and a third grating 4, a detector 5, an image.
- the processing unit 6, the control unit 7, the rotation mechanism 8, and the lattice moving mechanism 9 are provided.
- the direction from the X-ray source 1 to the first grating 2 is the Z2 direction
- the opposite direction is the Z1 direction.
- one of the two directions orthogonal to each other is defined as the X direction
- the direction toward the back of the paper is the X2 direction
- the direction toward the front of the paper is the X1 direction
- the other of the two directions orthogonal to each other is the Y direction
- the upward direction of the paper is the Y1 direction
- the downward direction is the Y2 direction.
- the X-ray source 1 is configured to generate X-rays by applying a high voltage and irradiate the generated X-rays in the Z2 direction. Further, the X-ray source 1 is provided with an X-ray filter 10.
- the X-ray filter 10 includes a metal plate and is configured to transmit only X-rays in a predetermined wavelength range.
- the first grating 2 has a plurality of slits 2a arranged at a predetermined period (pitch) p1 in the Y direction, and an X-ray phase changing portion 2b.
- Each of the slits 2a and the X-ray phase changing portion 2b is formed so as to extend linearly. Further, each slit 2a and X-ray phase changing portion 2b are formed so as to extend in parallel.
- the first grating 2 is a so-called phase grating.
- the first grating 2 is arranged 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 to form a self-image (not shown) of the first grating 2 by the Talbot effect.
- an image (self-image) of the grating is formed at a position separated from the grating by a predetermined distance (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 in the Y direction at a predetermined cycle (pitch) p2.
- the X-ray absorbing portion 3b extends along the direction in which the X-ray phase changing portion 2b extends.
- Each X-ray transmitting portion 3a and X-ray absorbing portion 3b is formed so as to extend linearly. Further, each X-ray transmitting portion 3a and X-ray absorbing portion 3b are formed so as to extend in parallel.
- the second grating 3 is a so-called absorption grating.
- the slit 2a and the X-ray transmitting portion 3a respectively transmit X-rays.
- the X-ray absorbing section 3b plays 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 the refractive index with the slit 2a.
- the second grating 3 is arranged between the first grating 2 and the detector 5 and is irradiated with the 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 the Talbot distance.
- the second grating 3 interferes with the self-image of the first grating 2 to form moire 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 in the Y direction at a predetermined cycle (pitch) p3. Each X-ray transmitting portion 4a and X-ray absorbing portion 4b is formed so as to extend linearly. Further, each X-ray transmitting portion 4a and X-ray absorbing portion 4b are 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.
- lattice 4 is comprised so that the X-ray which passed each X-ray transmission part 4a may be used as a line light source.
- the pitches 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. This is called the low effect. Thereby, the interference intensity can be maintained even if the tube size of the X-ray source 1 is large.
- the X-ray phase contrast imaging system 100 in the first embodiment is configured by the so-called Talbot-Lau interferometer.
- the detector 5 is configured to detect X-rays, convert the detected X-rays into an electric signal, and read the converted electric signal as an image signal.
- 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). Further, the detector 5 is configured to output the acquired image signal to the image processing unit 6.
- the image processing unit 6 is configured to generate an absorption image 20 (see FIG. 2) based on the X-ray intensity distribution detected by the detector 5. Further, the image processing unit 6 is configured to generate the phase differential image 21 (see FIG. 2) based on the intensity distribution of the X-ray detected by the detector 5. The image processing unit 6 is also configured to generate a dark-field image 22 (see FIG. 2) due to X-ray scattering, based on the X-ray intensity distribution detected by the detector 5.
- the absorption image 20 is an image of a contrast caused by a difference in absorption of X-rays by the subject T. That is, the absorption image 20 is an image of the X-ray transmittance of the subject T.
- the phase differential image 21 is an image formed based on the phase shift of X-rays generated when the X-rays pass through the subject T. That is, the phase differential image 21 is an image when the change in the phase of the X-ray is used as contrast.
- the dark field image 22 is an image in which the degree of refraction (scattering) of X-rays due to the fine structure inside the subject T is used as contrast. In other words, the dark field image 22 is an image of the decrease in the visibility V in the detector 5, and the change in the visibility V depends on the degree of scattering of the subject T.
- the image processing unit 6 generates a plurality of absorption images 20, a plurality of phase differential images 21, and a plurality of dark field images 22 which are captured while rotating the subject T by the rotating mechanism 8 (at each of a plurality of rotation angles).
- the three-dimensional absorption image 30 (see FIG. 2), the three-dimensional phase image 31 (see FIG. 2), and the three-dimensional dark field image 32 (see FIG. 2) are generated by reconstructing each.
- the image processing unit 6 includes a processor such as a GPU (Graphics Processing Unit) or an FPGA (Field-Programmable Gate Array) configured for image processing.
- the control unit 7 is configured of the subject T, the X-ray source 1, the detector 5, and a plurality of gratings (the first grating 2, the second grating 3, and the third grating 4) by the rotation mechanism 8. It is configured to relatively rotate the imaging system 11. Further, the control unit 7 is configured to move the first grating 2 stepwise in the lattice plane in a direction orthogonal to the lattice direction by the lattice moving mechanism 9.
- the X-ray phase contrast imaging system 100 uses a method (a fringe scanning method) of acquiring a reconstructed image from a plurality of moire fringes (images) obtained by scanning the first grating 2 at regular intervals. ..
- control unit 7 is configured to perform control for switching imaging conditions based on the magnitude of the degree of X-ray scattering. Details of the configuration in which the control unit 7 switches the shooting conditions will be described later.
- the control unit 7 includes a processor such as a CPU (Central Processing Unit), for example.
- the rotation mechanism 8 is configured to relatively rotate the subject T and the imaging system 11 based on a signal from the control unit 7. Specifically, the rotation mechanism 8 is configured to rotate the subject T about an axis AR extending along the Y direction, thereby rotating the subject T relative to the imaging system 11.
- the rotating mechanism 8 includes a rotating stage 8a driven by, for example, a motor.
- the image processing unit 6 reconstructs the plurality of absorption images 20, the plurality of phase differential images 21, and the plurality of dark field images 22 obtained by photographing the subject T at each of the plurality of rotation angles.
- a three-dimensional absorption image 30, a three-dimensional phase image 31, and a three-dimensional dark field image 33 are generated.
- the grating moving mechanism 9 is configured to move the first grating 2 stepwise in a direction orthogonal to the grating direction (Y direction in FIG. 1) in the grating plane (XY plane) based on the signal from the control unit 7. Has been done. Specifically, the grating moving mechanism 9 divides the period p1 of the first grating 2 into n, and moves the first grating 2 stepwise by p1 / n. The grating moving mechanism 9 is configured to move the first grating 2 stepwise by at least one cycle p1 of the first grating 2. The image processing unit 6 acquires the step curve SC (see FIG.
- n is a positive integer, such as 4, for example.
- the grid moving mechanism 9 also includes, for example, a stepping motor, a piezo actuator, or the like.
- the image processing unit 6 determines the ratio of the visibility V acquired by performing the image capturing in the state where the subject T is not arranged and the visibility V obtained by performing the image capturing in the state where the subject T is disposed. Based on this, the dark field image 22 is generated.
- photographing is performed while the subject T and the imaging system 11 are relatively rotated by the rotation mechanism 8.
- the path angle of X-rays passing through the subject T changes due to a change in the rotation angle when the subject T is taken.
- the degree of X-ray scattering by the subject T becomes greater.
- the example shown in FIG. 3 is an example in which the path length of the X-ray passing through the subject T is short.
- a subject T is a CFRP having a plate-like shape
- the rotation angle at which the X-rays are emitted from the shorter direction of the subject T is reached.
- the example shown in FIG. 4 is an example in which the path length of the X-ray passing through the subject T is long.
- the rotation angle is such that X-rays are emitted from the longitudinal direction of the subject T.
- the path length through which the X-rays emitted along the optical axis XR of the X-rays penetrate the subject T is defined as a distance d1. ..
- the path length through which the X-rays emitted along the optical axis XR of the X-rays pass through the subject T is defined as a distance d2. .. Since the subject T has a plate-like shape, the path length distance d1 through the subject T is shorter than the path length distance d2 through the subject T.
- the magnitude of the amplitude A1 of the step curve SC1 is the magnitude of the amplitude A2 of the step curve SC2 shown in the graph G2 of FIG. 4, as shown in the graph G1 of FIG. Will be bigger than In other words, as the degree of X-ray scattering increases, the amplitude A of the obtained step curve SC decreases.
- the magnitude of the amplitude A of the step curve SC becomes small, it becomes difficult to accurately extract the phase differential image 21 and the dark field image 22. Since the amplitude A of the step curve SC is the difference between Imax and Imin, it can be regarded as equal to the visibility V when the amount of X-ray absorption by the subject T is small. That is, when the X-ray is emitted at a rotation angle at which the path length of the X-ray passing through the subject T becomes long, the visibility V decreases.
- Graphs G3a and G3b shown in FIG. 5 are graphs showing the relationship between the path length and the visibility V through which the X-rays emitted from the X-ray source 1 pass through the subject T.
- the horizontal axis represents the path length through which X-rays pass
- the vertical axis represents the visibility.
- a graph G3a is a graph showing the relationship between the visibility V and the path length of the X-rays that pass through the subject T when the X-rays are irradiated with a predetermined effective energy (design energy). As shown in the graph G3a, the visibility V decreases as the X-ray path length increases. It is considered that this is because the degree of scattering of X-rays by the subject T increases as the path length of X-rays increases.
- the graph G3b shown in FIG. 5 is a graph when an X-ray having a larger effective energy than that of the graph G3a is irradiated.
- the X-rays are less likely to be scattered because the effective energy of the irradiated X-rays is large. Therefore, the visibility V when the path length of the X-rays passing through the subject T is 0 (zero) is smaller than the visibility V when the path length of the X-rays passing through the subject T in the graph G3a is 0 (zero). Become.
- the value of visibility V of the graph G3b becomes larger than the value of visibility V of the graph G3a. It is considered that this is because the increase in the effective energy of the X-rays makes it difficult for the X-rays to scatter, and thus makes it less likely to be affected by the degree of scattering of the X-rays. That is, when the path length of X-rays passing through the subject T is long, increasing the effective energy of the X-rays makes it less susceptible to scattering by the subject T, and thus the amplitude A of the step curve SC becomes smaller. Can be suppressed.
- the control unit 7 is configured to perform switching control of the imaging conditions based on the magnitude of the degree of X-ray scattering depending on the rotation angle when imaging the subject T.
- the imaging conditions include the voltage (tube voltage) applied to the X-ray source 1, the X-ray filter 10 provided in the X-ray source 1, the current (tube current) applied to the X-ray source 1, and the detector 5. Is a condition based on the accumulation time, the number of steps of moving the first grating 2 by the grating moving mechanism 9, and the like.
- FIG. 6 is a schematic view of the subject T placed on the rotation mechanism 8 as viewed from the Y1 direction.
- the control unit 7 sets the rotation angle range in which the degree of X-ray scattering by the subject T is small as the first rotation angle range FR.
- the rotation mechanism 8 is omitted for convenience.
- the control unit 7 performs imaging in the first rotation angle range FR of the plurality of rotation angles under the first imaging condition, and scatters X-rays by the subject T more than in the first rotation angle range FR.
- the second rotation angle range SR in which the degree is relatively large, control is performed to switch to the second imaging condition in which the effective energy of X-rays is relatively larger than the first imaging condition.
- the first imaging condition is determined based on the design energy of the X-ray phase contrast imaging system 100 (imaging system 11) such as the voltage applied to the X-ray source 1 and the X-ray filter 10 provided in the X-ray source 1. These are the taken shooting conditions.
- the control unit 7 sets the voltage (tube voltage) applied to the X-ray source 1 larger than that in the first imaging condition, so that the effective energy of the X-rays emitted toward the detector 5 is reduced. It is configured to be relatively large. Thereby, the control unit 7 is configured to perform control to switch to the second imaging condition in which the effective energy of X-ray is relatively larger than the first imaging condition.
- control unit 7 acquires the rotation angle for switching the imaging conditions at the rotation angle determined based on the pixel value of the dark field image 22 caused by the scattering of the X-rays generated at each of the plurality of rotation angles.
- the image processing unit 6 is configured to acquire the feature amount based on the pixel value of the dark field image 22.
- the feature amount for example, the minimum value or the average value of the pixel values of the region in which the subject T is shown in the dark field image 22 can be adopted.
- the image processing unit 6 is configured to acquire, as the feature amount, the minimum pixel value of the region in the dark field image 22 in which the subject T is captured.
- the control unit 7 is configured to determine the rotation angle for switching the shooting conditions based on the feature amount. Specifically, the control unit 7 is configured to determine the rotation angle for switching the shooting condition by comparing the feature amount with the threshold value. More specifically, as shown in FIG. 7, the control unit 7 controls the imaging of the subject T at a plurality of rotation angles according to the first imaging condition, and after the imaging is performed, the X-rays are acquired more than the first imaging condition. Is switched to the second shooting condition in which the effective energy is relatively large, and control is performed to shoot the subject T again in each of the rotation angle ranges (second rotation angle range SR) in which the feature amount is smaller than the first threshold ThrA. Is configured.
- the first threshold ThrA may be set to any value as long as it is within a range in which the dark field image 22 can be imaged.
- the first threshold ThrA is set to, for example, a value that is 10% of the pixel value when the subject T is photographed without being placed.
- the control unit 7 shoots while rotating the subject T 360 degrees according to the first shooting condition.
- the graph G4 is a graph in which the horizontal axis represents the rotation angle of the subject T and the vertical axis represents the feature amount. As shown in the graph G4, when the rotation angle of the subject T is near 90 degrees and around 270 degrees, the feature amount becomes smaller than the first threshold ThrA. Therefore, the control unit 7 controls to switch the shooting condition between the second shooting condition after the shooting of the subject T under the first shooting condition is completed, and at the rotation angle at which the feature amount is smaller than the first threshold value ThrA. , The subject T is photographed again.
- the feature amount is less than or equal to the first threshold ThrA in the range of 80 degrees or more and 100 degrees or less and the range of 260 degrees or more and 270 degrees or less. Therefore, the control unit 7 sets the range of 80 degrees or more and 100 degrees or less and the range of 260 degrees or more and 270 degrees or less as the second rotation angle range SR. As shown in a graph G5 shown in FIG. 7B, the control unit 7 switches to the second shooting condition and shoots the subject T in the second rotation angle range SR. In the example shown in the graph G5 of FIG. 7B, the feature amount is a value larger than the first threshold value ThrA in the entire second rotation angle range SR.
- the feature amount is larger than the first threshold value ThrA in the entire range of the second rotation angle range SR. It doesn't have to be. That is, in the second rotation angle range SR, it is sufficient that the minimum value of the feature amount at the time of shooting under the second shooting condition is larger than the minimum value of the feature amount at the time of shooting under the first shooting condition.
- the image processing unit 6 obtains, for example, the sensitivity of the pixel value of the dark field image 22 under the first shooting condition and the sensitivity of the pixel value of the dark field image 22 under the second shooting condition in advance, and the sensitivity becomes equal. It may be configured to correct. Further, the image processing unit 6 may be configured to perform correction by aligning the pixel values of the dark field image 22 at each rotation angle.
- the control unit 7 suppresses deterioration of the image quality of the dark-field image 22 due to noise (quantum noise) depending on the number of X-ray photons and noise from the detector 5. In order to do so, a correction process is performed. Specifically, when the subject T is imaged under the second imaging condition, the control unit 7 makes the current applied to the X-ray source 1 larger than that in the first imaging condition, and increases the accumulation time in the detector 5. And increasing the number of steps for translationally moving the first grating 2 by the grating moving mechanism 9 to control at least one of increasing the X-ray dose reaching the detector 5. Is configured to do.
- control unit 7 increases the current applied to the X-ray source 1, increases the accumulation time in the detector 5, and the number of steps when translating the first grating 2 by the grating moving mechanism 9.
- the control may be performed by combining some controls or by increasing all the controls.
- FIG. 8 a process in which the X-ray phase contrast imaging system 100 according to the first embodiment generates the three-dimensional dark field image 32 will be described. Note that the processing illustrated in FIG. 8 is triggered by the fact that the control unit 7 has acquired a signal to start shooting input by a user operation.
- step S1 the control unit 7 sets the shooting condition to the first shooting condition. If the shooting condition is already set to the first shooting condition, the process of step S1 can be omitted.
- step S2 the control unit 7 controls the rotation mechanism 8 to take an image while rotating the subject T by 360 degrees. That is, the control unit 7 performs the photographing with the rotation range of 360 degrees as the first rotation angle range FR. Thereafter, the process proceeds to step S3.
- step S3 the image processing unit 6 acquires a feature amount from the dark field image 22 in each rotation.
- the control unit 7 compares the feature amount of the dark field image 22 at each rotation angle with the first threshold ThrA.
- the control unit 7 stores the rotation angle at which the feature amount is equal to or less than the first threshold value ThrA in a storage unit (not shown) or the like. If there is a rotation angle at which the feature amount is less than or equal to the first threshold ThrA, the process proceeds to step S4. If there is no rotation angle at which the feature amount is equal to or less than the first threshold ThrA, the process proceeds to step S6.
- step S4 the control unit 7 switches the shooting condition to the second shooting condition.
- step S5 the control unit 7 rotates the subject T at the rotation angles stored in the storage unit or the like, and shoots at each rotation angle. That is, in step S2, the control unit 7 sets, as the second rotation angle range SR, an angle range in which the feature amount is equal to or less than the first threshold ThrA among the rotation angles photographed as the first rotation angle range FR, and the second photography. Take the picture again depending on the conditions.
- the control unit 7 also performs control to increase the X-ray dose reaching the detector 5 described above.
- step S6 the image processing unit 6 generates the dark field image 22 at each rotation angle. Further, the image processing unit 6 performs a correction for smoothly connecting the pixel values of the dark field image 22 described above.
- step S7 the image processing unit 6 reconstructs the dark field image 22 at each rotation angle to generate the three-dimensional dark field image 32, and ends the process.
- the image processing unit 6 generates the three-dimensional dark field image 32, the dark field image obtained by re-imaging under the second imaging condition at the rotation angle included in the second rotation angle range SR. 22 is used.
- step S6 the image processing unit 6 reconstructs the plurality of absorption images 20 and the plurality of phase differential images 21 captured at each rotation angle, respectively, to thereby obtain the three-dimensional absorption image 30 and the three-dimensional phase image 31.
- the image processing unit 6 captures the image again at the rotation angle included in the second rotation angle range SR under the second capturing condition, similarly to the generation of the three-dimensional dark field image 32.
- the phase differential image 21 obtained by doing is used.
- the X-ray phase contrast imaging system 100 includes the X-ray source 1, the detector 5 that detects the X-rays emitted from the X-ray source 1, and the X-ray source 1. It is composed of a plurality of gratings (first grating 2, second grating 3 and third grating 4) arranged between the detector 5 and the object T, the X-ray source 1, the detector 5 and a plurality of gratings. Based on the intensity distribution of the X-ray detected by the detector 5 at each of the rotation mechanism 8 that relatively rotates the imaging system 11 that rotates and the plurality of rotation angles when rotated by the rotation mechanism 8.
- the photographing conditions so as to suppress deterioration of the visibility V of X-rays and perform photographing.
- the visibility V from decreasing to near zero, and the CT image (three-dimensional dark field) of the reconstructed subject T is caused by the increase in the degree of X-ray scattering. It is possible to suppress deterioration of the image quality of the image 32).
- the control unit 7 performs shooting under the first shooting condition in the first rotation angle range FR of the plurality of rotation angles, and the shooting speed is lower than that in the first rotation angle range FR.
- control is performed to switch to the second imaging condition in which the effective X-ray energy is relatively larger than the first imaging condition.
- the imaging is performed under the second imaging condition in which the effective energy of X-rays is relatively larger than in the first imaging condition. This makes it difficult to scatter X-rays.
- the control unit 7 sets the rotation angle for switching the imaging conditions to the pixel value of the dark field image 22 caused by the scattering of the X-rays generated at each of the plurality of rotation angles. Is configured to be acquired at the rotation angle determined based on.
- the subject T has a structure in which the degree of X-ray scattering increases in a predetermined angle range, it is possible to prevent the visibility V from decreasing in the predetermined angle range.
- the dark-field image 22 in the angle range can be accurately generated.
- the image processing unit 6 is configured to acquire the feature amount based on the pixel value of the dark field image 22, and the control unit 7 is based on the feature amount.
- the rotation angle for switching the shooting condition is determined. Accordingly, the X-ray scattering degree for each detection pixel of the detector 5 at each of the plurality of rotation angles, or the pixel value for each pixel of the dark-field image 22 at each of the plurality of rotation angles is compared to obtain the imaging condition.
- the rotation angle for switching the imaging conditions is determined by determining the rotation angle for switching the imaging conditions according to the feature amount acquired based on the X-ray scattering degree or the pixel value in a plurality of pixels, as compared with the configuration for performing the determination for switching the imaging conditions. It is possible to reduce the amount of information used when determining As a result, it is possible to reduce the calculation load when making the determination of switching the shooting conditions.
- control unit 7 is configured to determine the rotation angle for switching the shooting condition by comparing the feature amount with the threshold value. Accordingly, by comparing the feature amount at each rotation angle with the threshold value, the degree of decrease in the visibility V of the X-ray at each rotation angle can be easily grasped. As a result, it is possible to easily determine the rotation angle for switching the shooting conditions.
- the control unit 7 controls the imaging of the subject T at a plurality of rotation angles according to the first imaging condition, and after the imaging, the X-rays are taken more than the first imaging condition. Is switched to the second shooting condition in which the effective energy is relatively large, and the control is performed again to shoot the subject T in each of the rotation angle ranges in which the feature amount is smaller than the first threshold ThrA.
- the shooting under the first condition and the shooting under the second shooting condition can be performed twice.
- the photographing is performed in real time while switching the first photographing condition and the second photographing condition based on the magnitude of the degree of X-ray scattering. Comparing with the configuration described above, it is possible to suppress the control from becoming complicated.
- the control unit 7 increases the voltage applied to the X-ray source 1 compared to the first imaging condition, so that the X-rays emitted toward the detector 5 are controlled. It is configured to control to switch to the second imaging condition by relatively increasing the effective energy.
- the effective energy of the X-rays emitted can be increased without replacing the X-ray source 1.
- the control unit 7 makes the current applied to the X-ray source 1 larger than that in the first imaging condition when imaging the subject T under the second imaging condition.
- the detector 5 is reached by performing at least one of increasing the accumulation time in the detector 5 and increasing the number of steps for translationally moving the first grating 2 by the grating moving mechanism 9.
- the X-ray dose is controlled to increase.
- the S / N ratio signal-to-noise ratio
- the image processing unit 6 sets the pixel values of the dark-field image 22 that are discontinuous due to switching of shooting conditions when shooting the subject T at a plurality of rotation angles. It is configured to perform a smooth connection correction. As a result, it is possible to suppress the occurrence of artifacts (noise) or the like due to the discontinuity of the pixel values of the dark field image 22 caused by the switching of the shooting conditions. As a result, for example, when the CT image (three-dimensional dark field image 32) is acquired by reconstructing the dark field image 22 at each rotation angle, the CT image (3 It is possible to suppress deterioration of the image quality of the three-dimensional dark field image 32).
- an X-ray phase contrast imaging system 200 (see FIG. 1) according to the second embodiment will be described with reference to FIGS. 1, 9 and 10.
- the control unit 70 shoots while rotating the subject T.
- the feature amount and the threshold value are compared with each other, and control is performed to perform shooting while switching between the first shooting condition and the second shooting condition.
- the same components as those in the first embodiment are designated by the same reference numerals and the description thereof will be omitted.
- the X-ray phase contrast imaging system 200 has the same configuration as the X-ray phase contrast imaging system 100 according to the first embodiment except that a control unit 70 is provided.
- the control unit 70 performs the second imaging condition in which the effective energy of the X-ray is relatively larger than the first imaging condition when the feature amount becomes smaller than the first threshold ThrA while performing the imaging under the first imaging condition.
- Control for switching to the first shooting condition and shooting when the feature amount becomes larger than the second threshold ThrB that is larger than the first threshold ThrA while shooting under the second shooting condition Is configured to do.
- the second threshold ThrB may be set to any value as long as it is within the range in which the dark field image 22 can be imaged and is larger than the first threshold ThrA.
- the second threshold ThrB is set to, for example, 20% of the pixel value when the subject T is photographed without being placed.
- the control unit 70 captures the first rotation angle range FR under the first capturing condition while the feature amount of the dark field image 22 (see FIG. 2) is the first threshold ThrA.
- the shooting condition is switched to the second shooting condition without changing the rotation angle.
- the subject T is photographed again under the second photographing condition.
- the subject T is photographed in the second rotation angle range SR under the second photographing condition.
- the control unit 70 determines whether or not the feature amount exceeds the second threshold value ThrB each time the image capturing is performed at each rotation angle.
- the control unit 70 switches the shooting condition to the first shooting condition and shoots again. That is, the control unit 70 is configured to perform control while performing shooting while switching between the first shooting condition and the second shooting condition in the first rotation angle range FR and the second rotation angle range SR.
- the imaging conditions are switched by changing the voltage applied to the X-ray source 1, it is preferable to provide a predetermined time interval until the X-ray source 1 stabilizes.
- the graph G7 of the feature amount is discontinuous. Therefore, like the control unit 7 in the first embodiment, the image processing unit 6 changes the pixel values of the dark field image 22 discontinuous due to the switching of the shooting conditions, as shown in the graph G8 in FIG. 9B. Perform a correction to smoothly connect the.
- FIG. 10 a process in which the X-ray phase contrast imaging system 200 according to the second embodiment generates the three-dimensional dark field image 32 will be described. Note that the processing shown in FIG. 10 is also started by the acquisition of a signal to start shooting input by the user's operation by the control unit 70, as in the first embodiment. Further, detailed description of steps that perform the same processing as the processing in the first embodiment will be omitted.
- step S1 the control unit 7 sets the shooting condition to the first shooting condition.
- the process of step S1 can be omitted.
- step S8 the subject T is photographed. Then, a process progresses to step S9.
- step S9 the control unit 70 determines whether or not it is necessary to switch the shooting conditions based on the feature amount of the dark field image 22. Specifically, when shooting under the first shooting condition, the control unit 70 compares the feature amount with the first threshold ThrA to determine whether the shooting condition needs to be switched. Further, the control unit 70 determines whether or not the switching of the shooting conditions is necessary by comparing the feature amount with the second threshold value ThrB while shooting under the second shooting conditions. If it is necessary to switch the shooting conditions, the process proceeds to step S10. If it is not necessary to switch the shooting conditions, the process proceeds to step S12.
- step S10 the control unit 70 switches shooting conditions. Specifically, the control unit 70 controls to switch to the second shooting condition when shooting is performed under the first shooting condition. In addition, the control unit 70 controls to switch to the first shooting condition when shooting is performed under the second shooting condition.
- step S11 the control unit 7 photographs the subject T without changing the rotation angle. Although the subject T is photographed without changing the rotation angle in step S11, the subject T may be photographed after rotating the subject T by a predetermined angle. In this case, the process of step S11 can be omitted. Then, a process progresses to step S12.
- step S12 the control unit 70 determines whether or not the subject T is photographed by rotating 360 degrees.
- the process proceeds to steps S6 and S7, and the image processing unit 6 generates the three-dimensional dark field image 32 and ends the process. If the subject T has not been photographed by rotating 360 degrees, the process proceeds to step S13.
- step S13 the control unit 70 controls the rotation mechanism 8 to rotate the subject T by a predetermined rotation angle. Thereafter, the processing proceeds to step S8.
- step S6 when the subject T is photographed under the second photographing condition, the control unit 70 may perform control to increase the X-ray dose reaching the detector 5. Good. Further, when the dark field image 22 is generated in step S7, the image processing unit 6 may perform correction for smoothly connecting the pixel values of the dark field image 22.
- the control unit 70 when the feature amount becomes smaller than the first threshold value ThrA while performing the imaging under the first imaging condition, the control unit 70 emits X-rays more than the first imaging condition.
- the characteristic amount becomes larger than the second threshold value ThrB which is larger than the first threshold value ThrA while performing the photographing under the second shooting condition while switching to the second shooting condition where the effective energy is relatively large.
- control unit 7 (70) determines the rotation angle for switching the shooting conditions based on the pixel value of the dark field image 22 . It is not limited to this.
- the control unit 7 (70) acquires the rotation angle for switching the imaging conditions at the rotation angle determined based on the preset value set based on the shape of the subject T contributing to the degree of X-ray scattering. May be configured. For example, when the subject T has a plate shape, as shown in FIG. 6, a path of X-rays that penetrates the subject T in a predetermined angle range from 90 degrees and 180 degrees of the rotation angle of the subject T. You can see that the length becomes longer.
- the user sets a predetermined angle range from 90 degrees and 180 degrees as the second rotation angle range SR.
- the angle range of 80 degrees to 100 degrees and the angle range of 260 degrees to 280 degrees are set as the second rotation angle range SR.
- the user also sets the other rotation angle range as the first rotation angle range FR.
- the control unit 7 (70) may be configured to perform control for switching the shooting conditions according to the rotation angle range set by the user.
- control unit 7 (70) determines the rotation angle for switching the photographing conditions based on the pixel value of the dark field image 22 in the first and second embodiments
- the present invention is not limited to this. It is not limited to this.
- the control unit 7 (70) is configured to acquire the rotation angle for switching the imaging condition at the rotation angle determined based on the degree of X-ray scattering detected by the detector 5 at each of the plurality of rotation angles. It may have been done. According to this structure, even when the rotation angle at which the degree of scattering of X-rays becomes large cannot be grasped in advance from the shape of the subject T, once the subject T is photographed, based on the measured value of the degree of scattering of X-rays.
- control unit 7 (70) determines the rotation angle for switching the photographing conditions based on the pixel value of the dark field image 22 in the first and second embodiments
- the present invention is not limited to this. It is not limited to this.
- the rotation angle for switching the shooting condition may be determined by the user designating a range in which the pixel value of the dark field image 22 is low.
- the image processing unit 6 acquires the feature amount based on the pixel value of the dark field image 22 is shown, but the present invention is not limited to this.
- the image processing unit 6 may be configured to acquire the feature amount based on the degree of X-ray scattering. Any value may be used as the feature amount as long as the magnitude of the amplitude A of the step curve SC can be grasped.
- the image processing unit 6 may acquire the sum of all pixel values of the dark field image 22 or the sum of logarithmic values of all pixel values of the dark field image 22 as the feature amount.
- the control unit 7 (70) shows an example of a configuration in which the voltage applied to the X-ray source 1 is increased to switch to the second imaging condition. Is not limited to this.
- the control unit 7 (70) may be configured to switch to the second imaging condition by replacing the X-ray filter 10 provided in the X-ray source 1.
- the control unit 7 (70) controls the X-ray wavelength that can be transmitted by the X-ray filter 10 to a lower wavelength side. By exchanging with a filter, the effective energy of X-rays can be increased.
- control unit 7 (70) shows an example of the configuration for performing the control for reducing the noise of the dark field image 22, but the present invention is not limited to this.
- the control unit 7 (70) does not have to perform the control for reducing the noise of the dark field image 22.
- the noise generated in the dark field image 22 increases, the image quality of the three-dimensional dark field image 32 deteriorates. Therefore, it is preferable to perform the control for reducing the noise of the dark field image 22.
- the image processing unit 6 shows an example of the configuration for performing correction for smoothly connecting the pixel values of the dark field image 22, but the present invention is not limited to this.
- the image processing unit 6 may be configured to perform correction for smoothly connecting the X-ray scattering degree.
- CT imaging may be performed by rotating the imaging system 11.
- the third lattice 4 is provided, but the present invention is not limited to this.
- the X-ray source 1 that emits X-rays having high coherence is used, the third grating 4 may not be provided.
- the dark field image 22 is generated by the fringe scanning method
- the present invention is not limited to this.
- a method of rotating any of the first grating 2, the second grating 3, or the third grating 4 in a plane orthogonal to the X-ray optical axis XR direction (Z direction) (so-called moire 1
- the dark-field image 22 may be generated by a single-shot method.
- 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.
- first grating 2 is step-moved in the lattice plane
- present invention is not limited to this. Any of the plurality of grids may be moved in steps.
- the present invention is not limited to this. Any subject T may be used as long as it is a subject T whose path length through which X-rays pass changes depending on the rotation angle of the subject T.
- CFRP carbon fiber reinforced plastic
- control unit 7 (70) shows an example of a configuration for performing control for switching between the first shooting condition and the second shooting condition, but the present invention is not limited to this. Absent.
- control unit 7 (70) may be configured to control switching of three or more shooting conditions.
- control unit 7 (70) has been described as an example of a configuration in which the shooting condition is controlled to suppress the deterioration of the image quality of the dark field image 22.
- the present invention is not limited to this.
- the control unit 7 (70) may be configured to perform control to switch the shooting conditions in order to suppress deterioration of the image quality of the phase differential image 21.
- the image processing unit 6 determines the average of the pixel values of the phase differential image 21 and the minimum of the pixel values of the phase differential image 21.
- a value, or a phase shift between the step curve SC when the subject T is photographed without arranging and the step curve SC when the subject T is photographed, may be acquired as a feature amount. .. Further, the control unit 7 (70) may be configured to perform control for switching the shooting conditions based on the feature amount of the phase differential image 21 acquired by the image processing unit 6.
- control unit 7 (70) is described using the “flow driven type” flowchart for convenience of description, but the present invention is not limited to this.
- the process of the control unit 7 (70) may be performed by an “event-driven type” in which the process is performed in event units.
- the event driving may be performed completely, or the event driving and the flow driving may be combined.
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
このX線位相差撮影システム(100)は、X線源(1)と、検出器(5)と、複数の格子と、被写体(T)と撮像系(11)とを相対的に回転させる回転機構(8)と、回転機構により回転された場合における複数の回転角度の各々において、少なくとも、X線の散乱に起因する暗視野像(22)を生成する画像処理部(6)と、被写体を撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行う制御部(7)とを備える。
Description
本発明は、X線位相差撮影システムに関し、特に、被写体と撮像系とを相対回転させながら撮影するX線位相差撮影システムに関する。
従来、被写体と撮像系とを相対回転させながら撮影するX線位相差撮影システムが知られている。このようなX線位相差撮影システムは、たとえば、M Bech, et al., Quantitative x-ray dark-field computed tomography, PHYSICS IN MEDICINE AND BIOLOGY, 55(2010), P.5529- 5539(以下、非特許文献1という)に開示されている。
上記非特許文献1のX線位相差撮影システムは、X線源と、位相格子と、検出器と、画像処理装置とを備えている。X線源からのX線は、被写体により散乱され、位相格子を通過して検出器に照射される。検出器により検出されるX線の干渉強度に基づいて、画像処理装置により暗視野像が生成される。また、被写体が360度回転された場合における各々の回転角度において撮影することにより暗視野像の生成が行われる。
M Bech, et al., Quantitative x-ray dark-field computed tomography, PHYSICS IN MEDICINE AND BIOLOGY, 55(2010), P.5529- 5539
しかしながら、上記非特許文献1には明記されていないが、上記非特許文献1に記載されているような従来のX線位相差撮影システムにおいて、被写体を回転させながら撮影する際に、回転角度に応じてX線が透過する経路長が異なる場合がある。たとえば、被写体として、板状形状を有する炭素繊維強化プラスチック(CFRP)を回転させながら撮影する際に、板状形状の長手方向に沿った特定範囲の回転角度(撮影角度)においてX線が透過する経路長が長くなることによってX線の散乱度合いが大きくなり、ビジビリティ(鮮明度)がゼロ付近まで低下する場合がある。ビジビリティがゼロ付近まで低下した場合、X線の干渉強度を正確に抽出することが困難になり、暗視野像の画質が劣化するという問題点がある。
この発明は、上記のような課題を解決するためになされたものであり、この発明の1つの目的は、回転角度に基づいてX線の散乱度合いが大きくなることに起因して、被写体の暗視野像の画質が劣化することを抑制することが可能なX線位相差撮影システムを提供することである。
上記目的を達成するために、この発明の一の局面におけるX線位相差撮影システムは、X線源と、X線源から照射されたX線を検出する検出器と、X線源と検出器との間に配置された複数の格子と、被写体と、X線源と検出器と複数の格子とによって構成される撮像系とを相対的に回転させる回転機構と、回転機構により回転された場合における複数の回転角度の各々において、検出器により検出されたX線の強度分布に基づいて、少なくとも、X線の散乱に起因する暗視野像を生成する画像処理部と、被写体を撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行う制御部とを備える。
この発明の一の局面におけるX線位相差撮影システムでは、上記のように、被写体を撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行う制御部を備える。これにより、たとえば、回転角度に基づいてX線の散乱度合いが大きくなった場合でも、X線のビジビリティが低下することを抑制することが可能な撮影条件に変更して撮影することができる。その結果、ビジビリティがゼロ付近まで低下することを抑制することが可能となるので、回転角度に基づいてX線の散乱度合いが大きくなることに起因して、被写体の暗視野像の画質が劣化することを抑制することができる。
上記一の局面におけるX線位相差撮影システムにおいて、好ましくは、制御部は、複数の回転角度のうち、第1回転角度範囲において、第1撮影条件によって撮影を行い、第1回転角度範囲よりも被写体によるX線の散乱度合いが相対的に大きくなる第2回転角度範囲では、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える制御を行うように構成されている。このように構成すれば、X線の散乱度合いが相対的に大きくなる第2回転角度範囲において撮影する場合に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件によって撮影することにより、X線を散乱しにくくすることができる。その結果、X線の散乱度合いが大きい被写体を撮影する場合でも、第1撮影条件によって撮影する場合と比較して、ビジビリティが低下することを抑制することが可能となるので、暗視野像を精度よく生成することができる。
上記一の局面におけるX線位相差撮影システムにおいて、好ましくは、制御部は、撮影条件を切り替える回転角度を、X線の散乱度合いに寄与する被写体の形状に基づいて予め設定された設定値、複数の回転角度の各々において検出器により検出されたX線の散乱度合い、または、複数の回転角度の各々において生成されたX線の散乱に起因する暗視野像の画素値のいずれかに基づいて決定した回転角度において取得するように構成されている。このように構成すれば、板状形状を有する被写体を撮影する場合など、X線の散乱度合いが大きくなる回転角度をユーザが予め把握できる場合には、被写体を撮影することなく撮影条件を切り替える角度を設定することができる。その結果、撮影条件を切り替える角度を取得するために事前に撮影を行う必要がなくなるので、撮影回数が増加することを抑制することができる。また、被写体の形状からではX線の散乱度合いが大きくなる回転角度が予め把握できない場合でも、一度被写体を撮影することによって、X線の散乱度合いまたは暗視野像の画素値の実測値に基づいて、正確に撮影条件を切り替える角度を決定することができる。その結果、被写体が所定の角度範囲でX線の散乱度合いが大きくなる構造を有している場合でも、所定の角度範囲においてビジビリティが低下することを抑制することが可能となるので、所定の角度範囲における暗視野像を精度よく生成することができる。また、X線の散乱度合いが大きくなる回転角度を予め把握できる場合でも、ユーザが入力することなく撮影条件を切り替えることができる。その結果、ユーザビリティ(ユーザの利便性)を向上させることができる。
この場合、好ましくは、画像処理部は、X線の散乱度合いまたは暗視野像の画素値に基づく特徴量を取得するように構成されており、制御部は、特徴量に基づいて、撮影条件を切り替える回転角度を決定するように構成されている。このように構成すれば、複数の回転角度の各々における検出器の検出画素ごとのX線の散乱度合い、または、複数の回転角度の各々における暗視野像の画素ごとの画素値を比較することにより撮影条件を切り替える判定を行う構成と比較して、複数の画素におけるX線の散乱度合いまたは画素値に基づいて取得された特徴量によって撮影条件を切り替える回転角度を決定することにより、撮影条件を切り替える回転角度を決定する際に用いる情報量を低減することができる。その結果、撮影条件を切り替える判定を行う際の計算負荷を軽減することができる。なお、特徴量とは、検出器におけるX線の散乱度合いの大きさの最小値または平均値、または、暗視野像において被写体が写っている領域の画素値の最小値または平均値などである。
上記特徴量に基づいて、撮影条件を切り替える回転角度を決定する構成において、好ましくは、制御部は、特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定するように構成されている。このように構成すれば、各回転角度における特徴量と閾値とを比較することにより、各回転角度におけるX線のビジビリティの低下度合いを容易に把握することができる。その結果、撮影条件を切り替える回転角度を容易に決定することができる。
この場合、好ましくは、制御部は、第1撮影条件によって複数の回転角度において被写体を撮影する制御を行うとともに、撮影した後に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて、特徴量が第1閾値よりも小さい回転角度範囲のそれぞれにおいて、再び被写体を撮影する制御を行うように構成されている。このように構成すれば、第1条件による撮影と、第2撮影条件による撮影とを2回に分けて行うことができる。その結果、撮影中において撮影条件を切り替える回数を1回にすることが可能となるので、X線の散乱度合いの大きさに基づいて第1撮影条件と第2撮影条件とを切り替えながらリアルタイムに撮影する構成と比較して、制御が複雑化することを抑制することができる。
上記特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定する構成において、好ましくは、制御部は、第1撮影条件による撮影をしながら、特徴量が第1閾値よりも小さくなった場合に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて撮影するとともに、第2撮影条件による撮影をしながら、特徴量が、第1閾値よりも大きい第2閾値よりも大きくなった場合に、第1撮影条件に切り替えて撮影する制御を行うように構成されている。このように構成すれば、1度の撮影中に第1撮影条件と第2撮影条件とを切り替えながらリアルタイムに撮像することができる。その結果、第1撮影条件によって撮影した後に、特徴量が第1閾値よりも小さくなった撮影角度範囲において第2撮影条件によって再び撮影する場合と比較して、撮影回数を低減することが可能となるので、撮影に要する時間を短縮することができる。
上記第1回転角度範囲よりも被写体によるX線の散乱度合いが相対的に大きくなる第2回転角度範囲において、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える構成おいて、好ましくは、制御部は、第1撮影条件よりもX線源に印加する電圧を大きくすること、および、X線源に設けるX線フィルタを変更することのうち、少なくともいずれかを行うことにより、検出器に向けて照射されるX線の実効エネルギーを相対的に大きくすることによって、第2撮影条件に切り替える制御を行うように構成されている。このように構成すれば、X線源に印加する電圧を大きくすることまたはX線源に設けるX線フィルタを変更することにより、X線源を交換することなく照射されるX線の実効エネルギーを大きくすることができる。その結果、たとえば、互いに異なるX線の実効エネルギーが予め設定されている複数のX線源を交換することにより撮影条件を切り替える構成と比較して、システム構成が複雑化することを抑制することができる。
この場合、好ましくは、複数の格子のいずれかを並進移動させる格子移動機構をさらに備え、制御部は、第2撮影条件において被写体を撮影する際に、第1撮影条件よりもX線源に印加する電流を大きくすること、検出器における蓄積時間を大きくすること、および、格子移動機構による複数の格子のいずれかを並進移動させる際のステップ数を大きくすることのうち、少なくともいずれかを行うことによって、検出器に到達するX線量を増加させる制御を行うように構成されている。このように構成すれば、各回転角度における暗視野像のS/N比(信号対雑音比(signal-to-noise ratio))を大きくすることができる。その結果、たとえば、各回転角度における暗視野像を再構成することによってCT画像を取得する場合、取得されるCT画像のS/N比を大きくすることが可能となるので、ノイズに起因してCT画像の画質が劣化することを抑制することができる。
上記一の局面におけるX線位相差撮影システムにおいて、好ましくは、画像処理部は、複数の回転角度において被写体を撮影する際に、撮影条件の切り替えによって不連続になったX線の散乱度合いまたは暗視野像の画素値を滑らかに連結する補正を行うように構成されている。このように構成すれば、撮影条件の切り替えにより生じたX線の散乱度合いの不連続または暗視野像の画素値の不連続に起因してアーチファクト(ノイズ)などが発生することを抑制することができる。その結果、たとえば、各回転角度における暗視野像を再構成することによってCT画像を取得する場合、暗視野像に生じるアーチファクトに起因して、CT画像の画質が劣化することを抑制することができる。
本発明によれば、上記のように、回転角度に基づいてX線の散乱度合いが大きくなることに起因して、被写体の暗視野像の画質が劣化することを抑制することができる。
以下、本発明を具体化した実施形態を図面に基づいて説明する。
[第1実施形態]
図1~図7を参照して、第1実施形態によるX線位相差撮影システム100の構成について説明する。
図1~図7を参照して、第1実施形態によるX線位相差撮影システム100の構成について説明する。
(X線位相差撮影システムの構成)
図1に示すように、X線位相差撮影システム100は、被写体Tを通過したX線の拡散(散乱)を利用して、被写体Tの暗視野像22(図2参照)を生成する。具体的には、X線位相差撮影システム100は、タルボ(Talbot)効果を利用して、少なくとも被写体Tの暗視野像22を生成する。X線位相差撮影システム100は、たとえば、非破壊検査用途では、物体の内部の画像化に用いることが可能である。第1実施形態では、被写体Tとして、たとえば、板状形状を有する炭素繊維強化プラスチック(CFRP)を撮影する例を示す。
図1に示すように、X線位相差撮影システム100は、被写体Tを通過したX線の拡散(散乱)を利用して、被写体Tの暗視野像22(図2参照)を生成する。具体的には、X線位相差撮影システム100は、タルボ(Talbot)効果を利用して、少なくとも被写体Tの暗視野像22を生成する。X線位相差撮影システム100は、たとえば、非破壊検査用途では、物体の内部の画像化に用いることが可能である。第1実施形態では、被写体Tとして、たとえば、板状形状を有する炭素繊維強化プラスチック(CFRP)を撮影する例を示す。
図1は、X線位相差撮影システム100をX方向から見た図である。図1に示すように、X線位相差撮影システム100は、X線源1と、第1格子2と第2格子3と第3格子4とを含む複数の格子と、検出器5と、画像処理部6と、制御部7と、回転機構8と、格子移動機構9とを備えている。なお、本明細書において、X線源1から第1格子2に向かう方向をZ2方向、その逆向きの方向をZ1方向とする。また、Z方向と直交する面内において、互いに直交する2方向のうち、一方方向をX方向とし、紙面の奥に向かう方向をX2方向、紙面の手前側に向かう方向をX1方向とする。また、Z方向と直交する面内において、互いに直交する2方向のうち、他方方向をY方向とし、紙面の上方向をY1方向、下方向をY2方向とする。
X線源1は、高電圧が印加されることにより、X線を発生させるとともに、発生されたX線をZ2方向に向けて照射するように構成されている。また、X線源1には、X線フィルタ10が設けられている。X線フィルタ10は、金属板を含み、所定の波長範囲のX線のみを透過するように構成されている。
第1格子2は、Y方向に所定の周期(ピッチ)p1で配列される複数のスリット2a、および、X線位相変化部2bを有している。各スリット2aおよびX線位相変化部2bはそれぞれ、直線状に延びるように形成されている。また、各スリット2aおよびX線位相変化部2bはそれぞれ、平行に延びるように形成されている。第1格子2は、いわゆる位相格子である。
第1格子2は、X線源1と、第2格子3との間に配置されており、X線源1からX線が照射される。第1格子2は、タルボ効果により、第1格子2の自己像(図示せず)を形成するために設けられている。なお、可干渉性を有するX線が、スリットが形成された格子を通過すると、格子から所定の距離(タルボ距離)離れた位置に、格子の像(自己像)が形成される。これをタルボ効果という。
第2格子3は、Y方向に所定の周期(ピッチ)p2で配列される複数の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方向に所定の周期(ピッチ)p3で配列される複数の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線を線光源とするように構成されている。3枚の格子(第1格子2、第2格子3、および、第3格子4)のピッチと格子間の距離とが一定の条件を満たすことにより、X線源1から照射されるX線の可干渉性を高めることが可能である。これを、ロー効果という。これにより、X線源1の管球の焦点サイズが大きくても干渉強度を保持できる。このように、第1実施形態におけるX線位相差撮影システム100は、いわゆるタルボ・ロー干渉計により構成される。
検出器5は、X線を検出するとともに、検出されたX線を電気信号に変換し、変換された電気信号を画像信号として読み取るように構成されている。検出器5は、たとえば、FPD(Flat Panel Detector)である。検出器5は、複数の変換素子(図示せず)と複数の変換素子上に配置された画素電極(図示せず)とにより構成されている。複数の変換素子および画素電極は、所定の周期(画素ピッチ)で、X方向およびY方向にアレイ状に配列されている。また、検出器5は、取得した画像信号を、画像処理部6に出力するように構成されている。
画像処理部6は、検出器5により検出されたX線の強度分布に基づいて、吸収像20(図2参照)を生成するように構成されている。また、画像処理部6は、検出器5により検出されたX線の強度分布に基づいて、位相微分像21(図2参照)を生成するように構成されている。また、画像処理部6は、検出器5により検出されたX線の強度分布に基づいて、X線の散乱に起因する暗視野像22(図2参照)を生成するように構成されている。ここで、吸収像20とは、被写体TによるX線の吸収の差によって生じるコントラストを画像化したものである。すなわち、吸収像20とは、被写体TによるX線の透過率を画像化したものである。また、位相微分像21とは、X線が被写体Tを通過した際に発生するX線の位相のずれをもとに画像化した像である。すなわち、位相微分像21とは、X線の位相の変化をコントラストとして画像化してものである。また、暗視野像22とは、被写体Tの内部にある微細構造によるX線の屈折(散乱)度合いをコントラストとして画像化したものである。言い換えると、暗視野像22は、検出器5におけるビジビリティVの低下を画像化したものであり、ビジビリティVの変化は被写体Tの散乱の程度に依存する。
また、画像処理部6は、回転機構8によって被写体Tを回転させながら(複数の回転角度の各々において)撮像された複数の吸収像20、複数の位相微分像21および複数の暗視野像22をそれぞれ再構成することにより、3次元吸収像30(図2参照)、3次元位相像31(図2参照)および3次元暗視野像32(図2参照)を生成する。また、画像処理部6は、たとえば、GPU(Graphics Processing Unit)や画像処理用に構成されたFPGA(Field-Programmable Gate Array)などのプロセッサを含む。
制御部7は、回転機構8により、被写体Tと、X線源1、検出器5、および、複数の格子(第1格子2、第2格子3、および、第3格子4)によって構成される撮像系11とを相対的に回転させるように構成されている。また、制御部7は、格子移動機構9により、第1格子2を格子面内において格子方向と直交する方向にステップ移動させるように構成されている。X線位相差撮影システム100では、第1格子2を一定周期間隔に走査することにより得られた複数のモアレ縞(画像)から再構成画像を取得する手法(縞走査法)が用いられている。また、制御部7は、X線の散乱度合いの大きさに基づいて、撮影条件を切り替える制御を行うように構成されている。制御部7が撮影条件を切り替える構成の詳細については後述する。制御部7は、たとえば、CPU(Central Processing Unit)などのプロセッサを含む。
回転機構8は、制御部7からの信号に基づいて、被写体Tと撮像系11とを相対的に回転させるように構成されている。具体的には、回転機構8は、被写体TをY方向に沿って延びる軸線AR周りに回転させることにより、撮像系11に対して被写体Tを相対的に回転させるように構成されている。回転機構8は、たとえば、モータなどによって駆動される回転ステージ8aを含む。
画像処理部6は、複数の回転角度のそれぞれにおいて被写体Tを撮影することにより取得した複数の吸収像20、複数の位相微分像21、および、複数の暗視野像22をそれぞれ再構成することにより、3次元吸収像30、3次元位相像31、および、3次元暗視野像33を生成する。
格子移動機構9は、制御部7からの信号に基づいて、第1格子2を格子面内(XY面内)において格子方向と直交する方向(図1ではY方向)にステップ移動させるように構成されている。具体的には、格子移動機構9は、第1格子2の周期p1をn分割し、p1/nずつ第1格子2をステップ移動させる。格子移動機構9は、少なくとも第1格子2の1周期p1分、第1格子2をステップ移動させるように構成されている。画像処理部6は、格子移動機構9によって第1格子2をステップ移動させながら撮影した際の各ステップにおけるX線の信号強度に基づいて、ステップカーブSC(図3参照)を取得する。なお、nは正の整数であり、たとえば、4などである。また、格子移動機構9は、たとえば、ステッピングモータやピエゾアクチュエータなどを含む。
(暗視野像の生成方法)
ここで、画像処理部6による暗視野像22の生成の方法について説明する。画像処理部6は、被写体Tが配置されていない状態で撮影を行うことにより取得されたビジビリティVと、被写体Tが配置されている状態で撮影を行うことにより取得されたビジビリティVとの比に基づいて、暗視野像22を生成する。なお、ビジビリティVは、鮮明度を意味する。また、ステップカーブSC(図3参照)における輝度値の最大値および最小値をそれぞれImaxおよびIminとすると、ビジビリティVは、V=(Imax-Imin)/(Imax+Imin)で表される。
ここで、画像処理部6による暗視野像22の生成の方法について説明する。画像処理部6は、被写体Tが配置されていない状態で撮影を行うことにより取得されたビジビリティVと、被写体Tが配置されている状態で撮影を行うことにより取得されたビジビリティVとの比に基づいて、暗視野像22を生成する。なお、ビジビリティVは、鮮明度を意味する。また、ステップカーブSC(図3参照)における輝度値の最大値および最小値をそれぞれImaxおよびIminとすると、ビジビリティVは、V=(Imax-Imin)/(Imax+Imin)で表される。
(被写体の撮影角度とビジビリティとの関係)
第1実施形態では、回転機構8によって被写体Tと撮像系11とを相対回転させながら撮影を行う。板状形状を有する被写体Tを撮影する場合、被写体Tを撮影する際の回転角度が変化することによって、被写体Tを透過するX線の経路長が変化する。被写体Tを透過する経路長が長くなると、被写体TによるX線の散乱度合いが大きくなる。
第1実施形態では、回転機構8によって被写体Tと撮像系11とを相対回転させながら撮影を行う。板状形状を有する被写体Tを撮影する場合、被写体Tを撮影する際の回転角度が変化することによって、被写体Tを透過するX線の経路長が変化する。被写体Tを透過する経路長が長くなると、被写体TによるX線の散乱度合いが大きくなる。
図3に示す例は、被写体Tを透過するX線の経路長が短い場合の例である。たとえば、板状形状を有するCFRPなどを被写体Tとする場合、被写体Tの短手方向からX線が照射される回転角度となった際の例である。図4に示す例は、被写体Tを透過するX線の経路長が長い場合の例である。たとえば、板状形状を有する被写体Tを撮影する場合、被写体Tの長手方向からX線が照射される回転角度となった際の例である。
図3に示すように、被写体Tの短手方向からX線が照射された場合に、X線の光軸XRに沿って照射されるX線が被写体Tを透過する経路長を距離d1とする。また、図4に示すように被写体Tの長手方向からX線が照射された場合に、X線の光軸XRに沿って照射されるX線が被写体Tを透過する経路長を距離d2とする。被写体Tは板状形状を有するため、被写体Tを透過する経路長の距離d1は、被写体Tを透過する経路長の距離d2よりも短くなる。
被写体Tを透過する経路長の距離d1が短い場合、図3のグラフG1に示すように、ステップカーブSC1の振幅A1の大きさは、図4のグラフG2に示すステップカーブSC2の振幅A2の大きさよりも大きくなる。言い換えると、X線の散乱度合いが大きくなると、得られるステップカーブSCの振幅Aの大きさが小さくなる。ステップカーブSCの振幅Aの大きさが小さくなると、位相微分像21および暗視野像22を正確に抽出することが困難になる。なお、ステップカーブSCの振幅Aは、ImaxとIminとの差であるため、被写体TによるX線吸収量が小さい場合は、ビジビリティVと等しいと見なすことができる。すなわち、被写体Tを透過するX線の経路長が長くなる回転角度においてX線を照射した場合、ビジビリティVが低下する。
(被写体を透過するX線の経路長とビジビリティとの関係)
図5を参照して、X線の経路長とビジビリティVとの関係について説明する。図5に示すグラフG3aおよびグラフG3bは、X線源1から照射されるX線が被写体Tを透過する経路長とビジビリティVとの関係を示したグラフである。グラフG3aおよびグラフG3bは、それぞれ、横軸が、X線が透過する経路長であり、縦軸がビジビリティである。所定の実効エネルギー(設計エネルギー)によってX線を照射した場合の被写体Tを透過するX線の経路長とビジビリティVとの関係を示すグラフが、グラフG3aである。グラフG3aに示すように、X線の経路長が増加するにつれて、ビジビリティVが低下する。これは、X線の経路長が長くなることにより、被写体TによるX線の散乱度合いが大きくなることが原因であると考えられる。
図5を参照して、X線の経路長とビジビリティVとの関係について説明する。図5に示すグラフG3aおよびグラフG3bは、X線源1から照射されるX線が被写体Tを透過する経路長とビジビリティVとの関係を示したグラフである。グラフG3aおよびグラフG3bは、それぞれ、横軸が、X線が透過する経路長であり、縦軸がビジビリティである。所定の実効エネルギー(設計エネルギー)によってX線を照射した場合の被写体Tを透過するX線の経路長とビジビリティVとの関係を示すグラフが、グラフG3aである。グラフG3aに示すように、X線の経路長が増加するにつれて、ビジビリティVが低下する。これは、X線の経路長が長くなることにより、被写体TによるX線の散乱度合いが大きくなることが原因であると考えられる。
図5に示すグラフG3bは、グラフG3aよりも実効エネルギーが大きいX線を照射した場合のグラフである。グラフG3bでは、照射されるX線の実効エネルギーが大きいため、X線が散乱されにくくなる。そのため、被写体Tを透過するX線の経路長が0(ゼロ)の場合のビジビリティVは、グラフG3aにおける被写体Tを透過するX線の経路長が0(ゼロ)の場合のビジビリティVよりも小さくなる。一方、被写体Tを透過するX線の経路長が増加するにつれて、グラフG3bのビジビリティVの値は、グラフG3aのビジビリティVの値よりも大きくなる。これは、X線の実効エネルギーが大きくなったことにより、X線が散乱しにくくなるため、X線の散乱度合いによる影響を受けにくくなるためであると考えられる。すなわち、被写体Tを透過するX線の経路長が長い場合に、X線の実効エネルギーを大きくすることにより、被写体Tによる散乱の影響を受けにくくなるため、ステップカーブSCの振幅Aが小さくなることを抑制することができる。
(撮影条件の切り替え)
そこで、第1実施形態では、制御部7は、被写体Tを撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行うように構成されている。なお、撮影条件とは、X線源1に印加される電圧(管電圧)、X線源1に設けられるX線フィルタ10、X線源1に印加される電流(管電流)、検出器5による蓄積時間、格子移動機構9による第1格子2をステップ移動させる回数などに基づく条件である。
そこで、第1実施形態では、制御部7は、被写体Tを撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行うように構成されている。なお、撮影条件とは、X線源1に印加される電圧(管電圧)、X線源1に設けられるX線フィルタ10、X線源1に印加される電流(管電流)、検出器5による蓄積時間、格子移動機構9による第1格子2をステップ移動させる回数などに基づく条件である。
図6は、回転機構8に配置された被写体TをY1方向から見た模式図である。第1実施形態では、制御部7は、被写体TによるX線の散乱度合いが小さい回転角度範囲を、第1回転角度範囲FRと設定する。なお、図6に示す例は、便宜上、回転機構8の図示を省略している。
第1実施形態では、制御部7は、複数の回転角度のうち、第1回転角度範囲FRにおいて、第1撮影条件によって撮影を行い、第1回転角度範囲FRよりも被写体TによるX線の散乱度合いが相対的に大きくなる第2回転角度範囲SRでは、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える制御を行うように構成されている。なお、第1撮影条件とは、X線位相差撮影システム100(撮像系11)の設計エネルギーに基づいて、X線源1に印加する電圧やX線源1に設けるX線フィルタ10などが決定された撮影条件である。
第1実施形態では、制御部7は、第1撮影条件よりもX線源1に印加する電圧(管電圧)を大きくすることにより、検出器5に向けて照射されるX線の実効エネルギーを相対的に大きくするように構成されている。これにより、制御部7は、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える制御を行うように構成されている。
また、制御部7は、撮影条件を切り替える回転角度を、複数の回転角度の各々において生成されたX線の散乱に起因する暗視野像22の画素値に基づいて決定した回転角度において取得するように構成されている。具体的には、画像処理部6は、暗視野像22の画素値に基づく特徴量を取得するように構成されている。特徴量としては、たとえば、暗視野像22において被写体Tが写っている領域の画素値の最小値または平均値を採用できる。第1実施形態では、画像処理部6は、特徴量として、暗視野像22において被写体Tが写っている領域の画素値の最小値を取得するように構成されている。
(撮影条件を切り替える回転角度の決定)
制御部7は、特徴量に基づいて、撮影条件を切り替える回転角度を決定するように構成されている。具体的には、制御部7は、特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定するように構成されている。より具体的には、図7に示すように、制御部7は、第1撮影条件によって複数の回転角度において被写体Tを撮影する制御を行うとともに、撮影した後に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて、特徴量が第1閾値ThrAよりも小さい回転角度範囲(第2回転角度範囲SR)のそれぞれにおいて、再び被写体Tを撮影する制御を行うように構成されている。なお、第1閾値ThrAは、暗視野像22を画像化可能な範囲であれば、任意の値に設定すればよい。第1実施形態では、第1閾値ThrAは、たとえば、被写体Tを配置せずに撮影した際の画素値の10%の値に設定する。
制御部7は、特徴量に基づいて、撮影条件を切り替える回転角度を決定するように構成されている。具体的には、制御部7は、特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定するように構成されている。より具体的には、図7に示すように、制御部7は、第1撮影条件によって複数の回転角度において被写体Tを撮影する制御を行うとともに、撮影した後に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて、特徴量が第1閾値ThrAよりも小さい回転角度範囲(第2回転角度範囲SR)のそれぞれにおいて、再び被写体Tを撮影する制御を行うように構成されている。なお、第1閾値ThrAは、暗視野像22を画像化可能な範囲であれば、任意の値に設定すればよい。第1実施形態では、第1閾値ThrAは、たとえば、被写体Tを配置せずに撮影した際の画素値の10%の値に設定する。
図7(A)に示すように、制御部7は、第1撮影条件によって被写体Tを360度回転させながら撮影を行う。グラフG4は、被写体Tの回転角度を横軸に、特徴量を縦軸にしたグラフである。グラフG4に示すように、被写体Tの回転角度が90度付近および270度付近になった場合、特徴量が第1閾値ThrAよりも小さくなる。そこで、制御部7は、第1撮影条件による被写体Tの撮影が終了した後に、撮影条件を、第2撮影条件を切り替える制御を行い、特徴量が第1閾値ThrAよりも小さくなった回転角度において、再び被写体Tを撮影する。
図7(A)に示す例では、たとえば、80度以上100度以下の範囲と、260度以上270度以下の範囲とにおいて、特徴量が第1閾値ThrA以下となっている。そこで、制御部7は、80度以上100度以下の範囲と、260度以上270度以下の範囲とを、それぞれ、第2回転角度範囲SRとする。図7(B)に示すグラフG5のように、制御部7は、第2回転角度範囲SRにおいて、第2撮影条件に切り替えて被写体Tを撮影する。なお、図7(B)のグラフG5に示す例では、第2回転角度範囲SRの全範囲において、特徴量が第1閾値ThrAよりも大きい値となっている。しかし、第2撮影条件における撮影では、第1撮影条件における撮影時よりも特徴量が大きくなっていればよいので、第2回転角度範囲SRの全範囲において特徴量が第1閾値ThrAよりも大きくなっていなくてもよい。すなわち、第2回転角度範囲SRにおいて、第2撮影条件における撮影時の特徴量の最小値が、第1撮影条件における撮影時の特徴量の最小値よりも大きくなっていればよい。
図7(A)および図7(B)に示すように、撮影条件を切り替えて被写体Tを撮影したため、X線の実効エネルギーが変化し、特徴量(暗視野像22の画素値)のグラフが不連続となる。したがって、各回転角度における暗視野像22を再構成することにより3次元暗視野像32を生成した場合、3次元暗視野像32において、アーチファクトが生じる原因となる。そこで、図7(C)に示すグラフG6のように、画像処理部6は、複数の回転角度において被写体Tを撮影する際に、撮影条件の切り替えによって不連続になった暗視野像22の画素値を滑らかに連結する補正を行うように構成されている。画像処理部6は、たとえば、第1撮影条件における暗視野像22の画素値の感度と、第2撮影条件における暗視野像22の画素値の感度とを予め取得しておき、感度が等しくなるように補正するように構成されていてもよい。また、画像処理部6は、各回転角度における暗視野像22の画素値を揃えることにより補正するように構成されていてもよい。
また、第1実施形態では、制御部7は、X線のフォトン数に依存するノイズ(量子ノイズ)や、検出器5由来のノイズに起因して暗視野像22の画質が劣化することを抑制するために、補正処理を行うように構成されている。具体的には、制御部7は、第2撮影条件において被写体Tを撮影する際に、第1撮影条件よりもX線源1に印加する電流を大きくすること、検出器5における蓄積時間を大きくすること、および、格子移動機構9による第1格子2を並進移動させる際のステップ数を大きくすることのうち、少なくともいずれかを行うことによって、検出器5に到達するX線量を増加させる制御を行うように構成されている。なお、制御部7は、X線源1に印加する電流を大きくすること、検出器5における蓄積時間を大きくすること、および、格子移動機構9による第1格子2を並進移動させる際のステップ数を大きくすることのうち、いくつかの制御を組み合わせて行ってもよいし、全ての制御を行ってもよい。
次に、図8を参照して、第1実施形態によるX線位相差撮影システム100が3次元暗視野像32を生成する処理について説明する。なお、図8に示す処理は、ユーザの操作によって入力された撮影開始の信号を制御部7が取得したことをトリガーとして開始される。
ステップS1において、制御部7は、撮影条件を第1撮影条件に設定する。なお、撮影条件がすでに第1撮影条件に設定されている場合には、ステップS1の処理は省略可能である。
次に、ステップS2において、制御部7は、回転機構8を制御することにより、被写体Tを360度回転させながら撮影を行う。すなわち、制御部7は、360度の回転範囲を、第1回転角度範囲FRとして撮影を行う。その後、処理は、ステップS3へ進む。
ステップS3において、画像処理部6は、各回転における暗視野像22から、特徴量を取得する。制御部7は、各回転角度における暗視野像22の特徴量と第1閾値ThrAとを比較する。制御部7は、特徴量が第1閾値ThrA以下となる回転角度を記憶部(図示せず)などに記憶する。特徴量が第1閾値ThrA以下となる回転角度がある場合には、処理は、ステップS4へ進む。特徴量が第1閾値ThrA以下となる回転角度がない場合には、処理は、ステップS6へ進む。
ステップS4において、制御部7は、撮影条件を第2撮影条件に切り替える。次に、ステップS5において、制御部7は、記憶部などに記憶しておいた回転角度に被写体Tを回転させ、各回転角度において撮影を行う。すなわち、制御部7は、ステップS2において、第1回転角度範囲FRとして撮影した回転角度のうち、特徴量が第1閾値ThrA以下となった角度範囲を第2回転角度範囲SRとし、第2撮影条件によって再度撮影を行う。また、制御部7は、上記した検出器5に到達するX線量を増加させる制御を行う。
次に、ステップS6において、画像処理部6は、各回転角度における暗視野像22を生成する。また、画像処理部6は、上記した暗視野像22の画素値を滑らかに連結する補正を行う。次に、ステップS7において、画像処理部6は、各回転角度における暗視野像22を再構成することにより、3次元暗視野像32を生成し、処理を終了する。なお、画像処理部6は、3次元暗視野像32を生成する際には、第2回転角度範囲SRに含まれる回転角度において、第2撮影条件によって再度撮影することにより得られた暗視野像22を用いる。
なお、ステップS6において、画像処理部6は、各回転角度において撮影された複数の吸収像20および複数の位相微分像21をそれぞれ再構成することにより、3次元吸収像30および3次元位相像31を生成する。画像処理部6は、3次元位相像31を生成する際には、3次元暗視野像32の生成と同等に、第2回転角度範囲SRに含まれる回転角度において、第2撮影条件によって再度撮影することにより得られた位相微分像21を用いる。
(第1実施形態の効果)
第1実施形態では、以下のような効果を得ることができる。
第1実施形態では、以下のような効果を得ることができる。
第1実施形態では、上記のように、X線位相差撮影システム100は、X線源1と、X線源1から照射されたX線を検出する検出器5と、X線源1と検出器5との間に配置された複数の格子(第1格子2、第2格子3および第3格子4)と、被写体Tと、X線源1と検出器5と複数の格子とによって構成される撮像系11とを相対的に回転させる回転機構8と、回転機構8により回転された場合における複数の回転角度の各々において、検出器5により検出されたX線の強度分布に基づいて、少なくとも、X線の散乱に起因する暗視野像22を生成する画像処理部6と、被写体Tを撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行う制御部7とを備える。これにより、X線の散乱度合いが大きくなった場合でも、X線のビジビリティVが低下することを抑制することが可能な撮影条件に変更して撮影することができる。その結果、ビジビリティVがゼロ付近まで低下することを抑制することが可能となるので、X線の散乱度合いが大きくなることに起因して、再構成された被写体TのCT画像(3次元暗視野像32)の画質が劣化することを抑制することができる。
また、第1実施形態では、上記のように、制御部7は、複数の回転角度のうち、第1回転角度範囲FRにおいて、第1撮影条件によって撮影を行い、第1回転角度範囲FRよりも被写体TによるX線の散乱度合いが相対的に大きくなる第2回転角度範囲SRでは、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える制御を行うように構成されている。これにより、X線の散乱度合いが相対的に大きくなる第2回転角度範囲SRにおいて撮影する場合に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件によって撮影することにより、X線を散乱しにくくすることができる。その結果、X線の散乱度合いが大きい被写体Tを撮影する場合でも、第1撮影条件によって撮影する場合と比較して、ビジビリティVが低下することを抑制することが可能となるので、暗視野像22を精度よく生成することができる。また、X線の散乱度合いが大きくなる回転角度を予め把握できる場合でも、ユーザが入力することなく撮影条件を切り替えることができる。その結果、ユーザビリティ(ユーザの利便性)を向上させることができる。
また、第1実施形態では、上記のように、制御部7は、撮影条件を切り替える回転角度を、複数の回転角度の各々において生成されたX線の散乱に起因する暗視野像22の画素値に基づいて決定した回転角度において取得するように構成されている。これにより、被写体Tの形状からではX線の散乱度合いが大きくなる回転角度が予め把握できない場合は、一度被写体Tを撮影することによって、暗視野像22の画素値の実測値に基づいて、正確に撮影条件を切り替える角度を決定することができる。その結果、被写体Tが所定の角度範囲でX線の散乱度合いが大きくなる構造を有している場合でも、所定の角度範囲においてビジビリティVが低下することを抑制することが可能となるので、所定の角度範囲における暗視野像22を精度よく生成することができる。
また、第1実施形態では、上記のように、画像処理部6は、暗視野像22の画素値に基づく特徴量を取得するように構成されており、制御部7は、特徴量に基づいて、撮影条件を切り替える回転角度を決定するように構成されている。これにより、複数の回転角度の各々における検出器5の検出画素ごとのX線の散乱度合い、または、複数の回転角度の各々における暗視野像22の画素ごとの画素値を比較することにより撮影条件を切り替える判定を行う構成と比較して、複数の画素におけるX線の散乱度合いまたは画素値に基づいて取得された特徴量によって撮影条件を切り替える回転角度を決定することにより、撮影条件を切り替える回転角度を決定する際に用いる情報量を低減することができる。その結果、撮影条件を切り替える判定を行う際の計算負荷を軽減することができる。
また、第1実施形態では、上記のように、制御部7は、特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定するように構成されている。これにより、各回転角度における特徴量と閾値とを比較することにより、各回転角度におけるX線のビジビリティVの低下度合いを容易に把握することができる。その結果、撮影条件を切り替える回転角度を容易に決定することができる。
また、第1実施形態では、上記のように、制御部7は、第1撮影条件によって複数の回転角度において被写体Tを撮影する制御を行うとともに、撮影した後に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて、特徴量が第1閾値ThrAよりも小さい回転角度範囲のそれぞれにおいて、再び被写体Tを撮影する制御を行うように構成されている。これにより、第1条件による撮影と、第2撮影条件による撮影とを2回に分けて行うことができる。その結果、撮影中において撮影条件を切り替える回数を1回にすることが可能となるので、X線の散乱度合いの大きさに基づいて第1撮影条件と第2撮影条件とを切り替えながらリアルタイムに撮影する構成と比較して、制御が複雑化することを抑制することができる。
また、第1実施形態では、上記のように、制御部7は、第1撮影条件よりもX線源1に印加する電圧を大きくすることにより、検出器5に向けて照射されるX線の実効エネルギーを相対的に大きくすることによって、第2撮影条件に切り替える制御を行うように構成されている。これにより、X線源1に印加する電圧を大きくすることによって、X線源1を交換することなく照射されるX線の実効エネルギーを大きくすることができる。その結果、たとえば、互いに異なるX線の実効エネルギーが予め設定されている複数のX線源1を交換することにより撮影条件を切り替える構成と比較して、システム構成が複雑化することを抑制することができる。
また、第1実施形態では、上記のように、制御部7は、第2撮影条件において被写体Tを撮影する際に、第1撮影条件よりもX線源1に印加する電流を大きくすること、検出器5における蓄積時間を大きくすること、および、格子移動機構9による第1格子2を並進移動させる際のステップ数を大きくすることのうち、少なくともいずれかを行うことによって、検出器5に到達するX線量を増加させる制御を行うように構成されている。これにより、各回転角度における暗視野像22のS/N比(信号対雑音比(signal-to-noise ratio))を大きくすることができる。その結果、たとえば、各回転角度における暗視野像22を再構成することによってCT画像(3次元暗視野像32)を取得する場合、取得されるCT画像(3次元暗視野像32)のS/N比を大きくすることが可能となるので、ノイズに起因してCT画像(3次元暗視野像32)の画質が劣化することを抑制することができる。
また、第1実施形態では、上記のように、画像処理部6は、複数の回転角度において被写体Tを撮影する際に、撮影条件の切り替えによって不連続になった暗視野像22の画素値を滑らかに連結する補正を行うように構成されている。これにより、撮影条件の切り替えにより生じた暗視野像22の画素値の不連続に起因してアーチファクト(ノイズ)などが発生することを抑制することができる。その結果、たとえば、各回転角度における暗視野像22を再構成することによってCT画像(3次元暗視野像32)を取得する場合、暗視野像22に生じるアーチファクトに起因して、CT画像(3次元暗視野像32)の画質が劣化することを抑制することができる。
[第2実施形態]
次に、図1、図9および図10を参照して、第2実施形態によるX線位相差撮影システム200(図1参照)について説明する。第1撮影条件によって撮影した後に、第2撮影条件に切り替えて撮影する第1実施形態とは異なり、第2実施形態では、制御部70(図1参照)は、被写体Tを回転させながら撮影する際に、特徴量と閾値とを比較し、第1撮影条件と第2撮影条件とを切り替えながら撮影する制御を行うように構成されている。なお、上記第1実施形態と同様の構成については同様の符号を付し、説明を省略する。
次に、図1、図9および図10を参照して、第2実施形態によるX線位相差撮影システム200(図1参照)について説明する。第1撮影条件によって撮影した後に、第2撮影条件に切り替えて撮影する第1実施形態とは異なり、第2実施形態では、制御部70(図1参照)は、被写体Tを回転させながら撮影する際に、特徴量と閾値とを比較し、第1撮影条件と第2撮影条件とを切り替えながら撮影する制御を行うように構成されている。なお、上記第1実施形態と同様の構成については同様の符号を付し、説明を省略する。
図1に示すように、X線位相差撮影システム200は、制御部70を備える点を除いて、上記第1実施形態によるX線位相差撮影システム100と同様の構成を有している。
制御部70は、第1撮影条件による撮影をしながら、特徴量が第1閾値ThrAよりも小さくなった場合に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて撮影するとともに、第2撮影条件による撮影をしながら、特徴量が、第1閾値ThrAよりも大きい第2閾値ThrBよりも大きくなった場合に、第1撮影条件に切り替えて撮影する制御を行うように構成されている。なお、第2閾値ThrBは、暗視野像22を画像化可能な範囲で、かつ、第1閾値ThrAのよりも大きければ、任意の値に設定すればよい。第2実施形態では、第2閾値ThrBは、たとえば、被写体Tを配置せずに撮影した際の画素値の20%の値に設定する。
図9(A)に示すように、制御部70は、第1回転角度範囲FRにおいて、第1撮影条件によって撮影をしながら、暗視野像22(図2参照)の特徴量が第1閾値ThrAよりも小さくなった場合、回転角度を変更せずに撮影条件を第2撮影条件に切り替える。そして、暗視野像22の特徴量が第1閾値ThrAよりも小さくなった回転角度において、第2撮影条件において再び被写体Tの撮影を行う。その後、第2回転角度範囲SRにおいて、第2撮影条件によって被写体Tの撮影を行う。制御部70は、各回転角度において撮影を行う度に、特徴量が第2閾値ThrBを超えているか否かの判定を行う。特徴量が第2閾値ThrBを超えた場合には、制御部70は、撮影条件を第1撮影条件に切り替えて、再び撮影を行う。すなわち、制御部70は、第1回転角度範囲FRと第2回転角度範囲SRとにおいて、第1撮影条件と第2撮影条件とを切り替えながら撮影を行う制御を行うように構成されている。なお、X線源1に印加する電圧を変更することにより撮影条件を切り替えた場合には、X線源1が安定するまでに、所定の時間間隔を設けることが好ましい。
第2実施形態においても、撮影条件を切り替えて撮影するため、特徴量のグラフG7は、不連続となる。そこで、画像処理部6は、第1実施形態における制御部7と同様に、図9(B)に示すグラフG8のように、撮影条件の切り替えによって不連続になった暗視野像22の画素値を滑らかに連結する補正を行う。
次に、図10を参照して、第2実施形態によるX線位相差撮影システム200が3次元暗視野像32を生成する処理について説明する。なお、図10に示す処理も、上記第1実施形態と同様に、ユーザの操作によって入力された撮影開始の信号を制御部70が取得したことをトリガーとして開始される。また、上記第1実施形態における処理と同様の処理を行うステップについては、詳細な説明は省略する。
ステップS1において、制御部7は、撮影条件を第1撮影条件に設定する。なお、上記第1実施形態と同様に、撮影条件がすでに第1撮影条件に設定されている場合には、ステップS1の処理は省略可能である。
次に、ステップS8において、被写体Tの撮影を行う。その後、処理は、ステップS9へ進む。
ステップS9において、制御部70は、暗視野像22の特徴量に基づいて、撮影条件の切り替えが必要か否かを判定する。具体的には、制御部70は、第1撮影条件において撮影している際には、特徴量と第1閾値ThrAとを比較することにより、撮影条件の切り替えが必要か否かを判定する。また、制御部70は、第2撮影条件において撮影している際には、特徴量と第2閾値ThrBとを比較することにより、撮影条件の切り替えが必要か否かを判定する。撮影条件の切り替えが必要な場合には、処理は、ステップS10へ進む。撮影条件の切り替えが必要ない場合には、処理は、ステップS12へ進む。
ステップS10において、制御部70は、撮影条件を切り替える。具体的には、制御部70は、第1撮影条件によって撮影していた場合には、第2撮影条件に切り替える制御を行う。また、制御部70は、第2撮影条件によって撮影していた場合には、第1撮影条件に切り替える制御を行う。次に、ステップS11において、制御部7は、回転角度を変更することなく、被写体Tを撮影する。ステップS11では回転角度を変更することなく、被写体Tを撮影しているが、所定角度分回転させてから、被写体Tを撮影しても良い。この場合は、ステップS11の処理は省略可能である。その後、処理は、ステップS12へ進む。
次に、ステップS12において、制御部70は、360度回転させて被写体Tを撮影したか否かの判定を行う。360度回転させて被写体Tを撮影した場合には、処理は、ステップS6、ステップS7へと進み、画像処理部6は、3次元暗視野像32を生成して処理を終了する。360度回転させて被写体Tを撮影していない場合には、処理は、ステップS13へ進む。
ステップS13において、制御部70は、回転機構8を制御することにより、被写体Tを所定の回転角度分だけ回転させる。その後、処理は、ステップS8へ進む。
なお、上記第1実施形態と同様に、ステップS6において、第2撮影条件によって被写体Tを撮影する際には、制御部70は、検出器5に到達するX線量を増加させる制御を行ってもよい。また、ステップS7において暗視野像22を生成する際には、画像処理部6は、暗視野像22の画素値を滑らかに連結する補正を行ってもよい。
なお、第2実施形態のその他の構成は、上記第1実施形態と同様である。
(第2実施形態の効果)
第2実施形態では、以下のような効果を得ることができる。
第2実施形態では、以下のような効果を得ることができる。
第2実施形態では、上記のように、制御部70は、第1撮影条件による撮影をしながら、特徴量が第1閾値ThrAよりも小さくなった場合に、第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて撮影するとともに、第2撮影条件による撮影をしながら、特徴量が、第1閾値ThrAよりも大きい第2閾値ThrBよりも大きくなった場合に、第1撮影条件に切り替えて撮影する制御を行うように構成されている。これにより、1度の撮影中に第1撮影条件と第2撮影条件とを切り替えながらリアルタイムに撮像することができる。その結果、第1撮影条件によって撮影した後に、特徴量が第1閾値ThrAよりも小さくなった撮影角度範囲において第2撮影条件によって再び撮影する場合と比較して、撮影回数を低減することが可能となるので、撮影に要する時間を短縮することができる。
なお、第2実施形態のその他の効果は、上記第1実施形態と同様である。
(変形例)
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく、請求の範囲によって示され、さらに請求の範囲と均等の意味および範囲内でのすべての変更(変形例)が含まれる。
なお、今回開示された実施形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施形態の説明ではなく、請求の範囲によって示され、さらに請求の範囲と均等の意味および範囲内でのすべての変更(変形例)が含まれる。
たとえば、上記第1および第2実施形態では、制御部7(70)が暗視野像22の画素値に基づいて、撮影条件を切り替える回転角度を決定する構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、撮影条件を切り替える回転角度を、X線の散乱度合いに寄与する被写体Tの形状に基づいて予め設定された設定値に基づいて決定した回転角度において取得するように構成されていてもよい。たとえば、被写体Tが板状の形状を有している場合、図6に示すように、被写体Tの回転角度が90度および180度から所定の角度範囲において、被写体Tを透過するX線の経路長が長くなることがわかる。そこで、ユーザは、90度および180度から所定の角度範囲を、第2回転角度範囲SRと設定する。図6に示す例では、たとえば、80度から100度の角度範囲、および、260度から280度の角度範囲を、第2回転角度範囲SRと設定する。また、ユーザは、それ以外の回転角度範囲を、第1回転角度範囲FRと設定する。制御部7(70)は、ユーザが設定した回転角度範囲に応じて、撮影条件を切り替える制御を行うように構成されればよい。
上記のように構成すれば、板状形状を有する被写体Tを撮影する場合など、X線の散乱度合いが大きくなる回転角度をユーザが予め把握できる場合には、被写体Tを撮影することなく撮影条件を切り替える角度を設定することができる。その結果、撮影条件を切り替える角度を取得するために事前に撮影を行う必要がなくなるので、撮影回数が増加することを抑制することができる。
また、上記第1および第2実施形態では、制御部7(70)が暗視野像22の画素値に基づいて、撮影条件を切り替える回転角度を決定する構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、撮影条件を切り替える回転角度を、複数の回転角度の各々において検出器5により検出されたX線の散乱度合いに基づいて決定した回転角度において取得するように構成されていてもよい。このように構成すれば、被写体Tの形状からではX線の散乱度合いが大きくなる回転角度が予め把握できない場合でも、一度被写体Tを撮影することによって、X線の散乱度合いの実測値に基づいて、正確に撮影条件を切り替える角度を決定することができる。その結果、被写体Tが所定の角度範囲でX線の散乱度合いが大きくなる構造を有している場合でも、ビジビリティVが低下することを抑制することが可能となるので、暗視野像22を精度よく生成することができる。
また、上記第1および第2実施形態では、制御部7(70)が暗視野像22の画素値に基づいて、撮影条件を切り替える回転角度を決定する構成の例を示したが、本発明はこれに限られない。たとえば、撮影条件を切り替える回転角度は、暗視野像22の画素値が低い範囲をユーザが指定することにより決定されてもよい。
また、上記第1および第2実施形態では、画像処理部6が、暗視野像22の画素値に基づいて特徴量を取得する構成の例を示したが、本発明はこれに限られない。たとえば、画像処理部6は、X線の散乱度合いに基づいて、特徴量を取得するように構成されていてもよい。ステップカーブSCの振幅Aの大きさを把握することができれば、どのような値を特徴量としてもよい。たとえば、画像処理部6は、暗視野像22の全画素値の和や、暗視野像22の全画素値の対数値の和などを特徴量として取得してもよい。
また、上記第1および第2実施形態では、制御部7(70)が、X線源1に印加する電圧を大きくすることにより、第2撮影条件に切り替える構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、X線源1に設けられるX線フィルタ10を交換することにより、第2撮影条件に切り替えるように構成されていてもよい。X線フィルタ10を交換することにより第2撮影条件に切り替える構成の場合、制御部7(70)は、X線フィルタ10が透過可能なX線の波長をより低波長側に設計されたX線フィルタに交換することにより、X線の実効エネルギーを大きくすることができる。
また、上記第1および第2実施形態では、制御部7(70)が、暗視野像22のノイズを低減する制御を行う構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、暗視野像22のノイズを低減する制御を行わなくてもよい。しかしながら、暗視野像22に生じるノイズが増加すると、3次元暗視野像32の画質が劣化するため、暗視野像22のノイズを低減する制御を行う方が好ましい。
また、上記第1および第2実施形態では、画像処理部6が、暗視野像22の画素値を滑らかに連結する補正を行う構成の例を示したが、本発明はこれに限られない。たとえば、画像処理部6は、X線の散乱度合いが滑らかに連結する補正を行うように構成されていてもよい。
また、上記第1および第2実施形態では、被写体Tを回転機構8により回転させてCT撮影を行う例を示したが、本発明はこれに限られない。たとえば、撮像系11を回転させてCT撮影を行ってもよい。
また、上記第1および第2実施形態では、第3格子4が設けられている例を示したが、本発明はこれに限られない。可干渉性が高いX線を照射するX線源1を用いる場合、第3格子4が設けられていなくてもよい。
また、上記第1および第2実施形態では、縞走査法によって暗視野像22を生成する例を示したが、本発明はこれに限られない。たとえば、第1格子2、第2格子3、または、第3格子4のうちのいずれかを、X線の光軸XR方向(Z方向)に直交する平面内において回転させる手法(いわゆる、モアレ1枚撮り手法)によって、暗視野像22を生成してもよい。
また、上記第1および第2実施形態では、第1格子2が位相格子である例を示したが、本発明はこれに限られない。たとえば、第1格子2は吸収格子であってもよい。
また、上記第1および第2実施形態では、第1格子2を格子面内においてステップ移動させる例を示したが、本発明はこれに限られない。複数の格子のうち、いずれの格子をステップ移動させてもよい。
また、上記第1および第2実施形態では、被写体Tとして、板状形状を有する炭素繊維強化プラスチック(CFRP)を撮像する例を示したが、本発明はこれに限られない。被写体Tの回転角度に応じてX線が透過する経路長が変化する被写体Tであれば、どのようなものを被写体Tとしてもよい。
また、上記第1および第2実施形態では、制御部7(70)が、第1撮影条件と第2撮影条件とを切り替える制御を行う構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、3つ以上の撮影条件を切り替える制御を行うように構成されていてもよい。
また、上記第1および第2実施形態では、制御部7(70)が、暗視野像22の画質が劣化することを抑制するために撮影条件を切り替える制御を行う構成の例を示したが、本発明はこれに限られない。たとえば、制御部7(70)は、位相微分像21の画質が劣化することを抑制するために撮影条件を切り替える制御を行うように構成されていてもよい。位相微分像21の画質が劣化することを抑制するために撮影条件を切り替える制御を行う場合、画像処理部6は、位相微分像21の画素値の平均値、位相微分像21の画素値の最小値、または被写体Tを配置せずに撮影した際のステップカーブSCと、被写体Tを配置して撮影した際のステップカーブSCとの位相のずれなどを特徴量として取得するように構成すればよい。また、制御部7(70)は、画像処理部6が取得した位相微分像21の特徴量に基づいて、撮影条件を切り替える制御を行うように構成すればよい。
また、上記第1および第2実施形態では、説明の便宜上、制御部7(70)の処理を「フロー駆動型」のフローチャートを用いて説明したが、本発明はこれに限られない。制御部7(70)の処理をイベント単位で実行する「イベント駆動型」により行ってもよい。この場合、完全なイベント駆動型で行ってもよいし、イベント駆動およびフロー駆動を組み合わせて行ってもよい。
1 X線源
2 第1格子(複数の格子)
3 第2格子(複数の格子)
4 第3格子(複数の格子)
5 検出器
6 画像処理部
7、70 制御部
8 回転機構
9 格子移動機構
10 X線フィルタ
11 撮像系
22 暗視野像
100、200 X線位相差撮影システム
FR 第1回転角度範囲
SR 第2回転角度範囲
T 被写体
ThrA 第1閾値
ThrB 第2閾値
2 第1格子(複数の格子)
3 第2格子(複数の格子)
4 第3格子(複数の格子)
5 検出器
6 画像処理部
7、70 制御部
8 回転機構
9 格子移動機構
10 X線フィルタ
11 撮像系
22 暗視野像
100、200 X線位相差撮影システム
FR 第1回転角度範囲
SR 第2回転角度範囲
T 被写体
ThrA 第1閾値
ThrB 第2閾値
Claims (10)
- X線源と、
前記X線源から照射されたX線を検出する検出器と、
前記X線源と前記検出器との間に配置された複数の格子と、
被写体と、前記X線源と前記検出器と前記複数の格子とによって構成される撮像系とを相対的に回転させる回転機構と、
前記回転機構により回転された場合における複数の回転角度の各々において、前記検出器により検出されたX線の強度分布に基づいて、少なくとも、X線の散乱に起因する暗視野像を生成する画像処理部と、
被写体を撮影する際の回転角度によるX線の散乱度合いの大きさに基づいて撮影条件を切り替える制御を行う制御部とを備える、X線位相差撮影システム。 - 前記制御部は、複数の回転角度のうち、第1回転角度範囲において、第1撮影条件によって撮影を行い、第1回転角度範囲よりも被写体によるX線の散乱度合いが相対的に大きくなる第2回転角度範囲では、前記第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替える制御を行うように構成されている、請求項1に記載のX線位相差撮影システム。
- 前記制御部は、撮影条件を切り替える回転角度を、X線の散乱度合いに寄与する被写体の形状に基づいて予め設定された設定値、複数の回転角度の各々において前記検出器により検出されたX線の散乱度合い、または、複数の回転角度の各々において生成されたX線の散乱に起因する前記暗視野像の画素値のいずれかに基づいて決定した回転角度において取得するように構成されている、請求項1に記載のX線位相差撮影システム。
- 前記画像処理部は、前記X線の散乱度合いまたは前記暗視野像の画素値に基づく特徴量を取得するように構成されており、
前記制御部は、前記特徴量に基づいて、撮影条件を切り替える回転角度を決定するように構成されている、請求項3に記載のX線位相差撮影システム。 - 前記制御部は、前記特徴量と閾値とを比較することにより、撮影条件を切り替える回転角度を決定するように構成されている、請求項4に記載のX線位相差撮影システム。
- 前記制御部は、第1撮影条件によって複数の回転角度において被写体を撮影する制御を行うとともに、撮影した後に、前記第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて、前記特徴量が第1閾値よりも小さい回転角度範囲のそれぞれにおいて、再び被写体を撮影する制御を行うように構成されている、請求項5に記載のX線位相差撮影システム。
- 前記制御部は、第1撮影条件による撮影をしながら、前記特徴量が第1閾値よりも小さくなった場合に、前記第1撮影条件よりもX線の実効エネルギーが相対的に大きい第2撮影条件に切り替えて撮影するとともに、前記第2撮影条件による撮影をしながら、前記特徴量が、前記第1閾値よりも大きい第2閾値よりも大きくなった場合に、前記第1撮影条件に切り替えて撮影する制御を行うように構成されている、請求項5に記載のX線位相差撮影システム。
- 前記制御部は、前記第1撮影条件よりも前記X線源に印加する電圧を大きくすること、および、前記X線源に設けるX線フィルタを変更することのうち、少なくともいずれかを行うことにより、前記検出器に向けて照射されるX線の実効エネルギーを相対的に大きくすることによって、前記第2撮影条件に切り替える制御を行うように構成されている、請求項2に記載のX線位相差撮影システム。
- 前記複数の格子のいずれかを並進移動させる格子移動機構をさらに備え、
前記制御部は、前記第2撮影条件において被写体を撮影する際に、前記第1撮影条件よりも前記X線源に印加する電流を大きくすること、前記検出器における蓄積時間を大きくすること、および、前記格子移動機構による前記複数の格子のいずれかを並進移動させる際のステップ数を大きくすることのうち、少なくともいずれかを行うことによって、前記検出器に到達するX線量を増加させる制御を行うように構成されている、請求項8に記載のX線位相差撮影システム。 - 前記画像処理部は、複数の回転角度において被写体を撮影する際に、撮影条件の切り替えによって不連続になった前記X線の散乱度合いまたは前記暗視野像の画素値を滑らかに連結する補正を行うように構成されている、請求項1に記載のX線位相差撮影システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020554766A JP7163969B2 (ja) | 2018-10-31 | 2019-07-25 | X線位相差撮影システム |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018-205454 | 2018-10-31 | ||
JP2018205454 | 2018-10-31 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020090168A1 true WO2020090168A1 (ja) | 2020-05-07 |
Family
ID=70462561
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2019/029232 WO2020090168A1 (ja) | 2018-10-31 | 2019-07-25 | X線位相差撮影システム |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7163969B2 (ja) |
WO (1) | WO2020090168A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114569145A (zh) * | 2022-02-28 | 2022-06-03 | 西安大医集团股份有限公司 | 图像校正方法、成像系统、电子设备和存储介质 |
WO2023149113A1 (ja) * | 2022-02-03 | 2023-08-10 | 国立大学法人東北大学 | 弾性率算出方法及び弾性率算出装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07204189A (ja) * | 1994-01-26 | 1995-08-08 | Hitachi Medical Corp | X線ct装置 |
JP2003079611A (ja) * | 2001-07-04 | 2003-03-18 | Toshiba Corp | X線コンピュータ断層診断装置 |
JP2003290214A (ja) * | 2002-03-27 | 2003-10-14 | Ge Medical Systems Global Technology Co Llc | 透過x線データ獲得装置およびx線断層像撮影装置 |
JP2006214879A (ja) * | 2005-02-03 | 2006-08-17 | Sony Corp | X線断層撮像方法及び装置 |
JP2009279289A (ja) * | 2008-05-26 | 2009-12-03 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2013524897A (ja) * | 2010-04-19 | 2013-06-20 | 清華大学 | X線源回折格子のステップ撮像システムおよび撮像方法 |
US20170227476A1 (en) * | 2014-11-04 | 2017-08-10 | Tsinghua University | X-ray phase-contrast imaging system and imaging method |
JP2017524139A (ja) * | 2014-08-13 | 2017-08-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | トモグラフィにおける暗視野イメージング |
-
2019
- 2019-07-25 JP JP2020554766A patent/JP7163969B2/ja active Active
- 2019-07-25 WO PCT/JP2019/029232 patent/WO2020090168A1/ja active Application Filing
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07204189A (ja) * | 1994-01-26 | 1995-08-08 | Hitachi Medical Corp | X線ct装置 |
JP2003079611A (ja) * | 2001-07-04 | 2003-03-18 | Toshiba Corp | X線コンピュータ断層診断装置 |
JP2003290214A (ja) * | 2002-03-27 | 2003-10-14 | Ge Medical Systems Global Technology Co Llc | 透過x線データ獲得装置およびx線断層像撮影装置 |
JP2006214879A (ja) * | 2005-02-03 | 2006-08-17 | Sony Corp | X線断層撮像方法及び装置 |
JP2009279289A (ja) * | 2008-05-26 | 2009-12-03 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2013524897A (ja) * | 2010-04-19 | 2013-06-20 | 清華大学 | X線源回折格子のステップ撮像システムおよび撮像方法 |
JP2017524139A (ja) * | 2014-08-13 | 2017-08-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | トモグラフィにおける暗視野イメージング |
US20170227476A1 (en) * | 2014-11-04 | 2017-08-10 | Tsinghua University | X-ray phase-contrast imaging system and imaging method |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023149113A1 (ja) * | 2022-02-03 | 2023-08-10 | 国立大学法人東北大学 | 弾性率算出方法及び弾性率算出装置 |
CN114569145A (zh) * | 2022-02-28 | 2022-06-03 | 西安大医集团股份有限公司 | 图像校正方法、成像系统、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2020090168A1 (ja) | 2021-09-09 |
JP7163969B2 (ja) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7044764B2 (ja) | X線撮像のための線源格子 | |
JP5331940B2 (ja) | 放射線撮影システム及び放射線画像生成方法 | |
US8451975B2 (en) | Radiographic system, radiographic method and computer readable medium | |
JP5475737B2 (ja) | 放射線撮影装置及び画像処理方法 | |
JP5238786B2 (ja) | 放射線撮影装置及び放射線撮影システム | |
JP5343065B2 (ja) | 放射線撮影システム | |
JP2012090945A (ja) | 放射線検出装置、放射線撮影装置、放射線撮影システム | |
JP2012115576A (ja) | 放射線画像検出装置、放射線撮影装置、放射線撮影システム | |
JP2012200567A (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP2012120653A (ja) | 放射線撮影装置、及び放射線撮影システム | |
WO2020090168A1 (ja) | X線位相差撮影システム | |
JP2011206490A (ja) | 放射線撮影システム及び放射線撮影方法 | |
WO2013065502A1 (ja) | 放射線撮影方法及び装置 | |
JP2012125423A (ja) | 放射線画像検出装置、放射線撮影装置、放射線撮影システム | |
WO2012057278A1 (ja) | 放射線撮影システム及び放射線撮影方法 | |
JP6642676B2 (ja) | 放射線撮影システム及び画像処理装置 | |
WO2013038881A1 (ja) | 放射線撮影装置及び画像処理方法 | |
JP2013042788A (ja) | 放射線撮影装置及びアンラップ処理方法 | |
WO2012147528A1 (ja) | 放射線撮影装置及びその制御方法 | |
JP2012115621A (ja) | 放射線画像検出装置、放射線撮影装置、放射線撮影システム | |
WO2013027519A1 (ja) | 放射線撮影装置及びアンラップ処理方法 | |
WO2013027536A1 (ja) | 放射線撮影装置及び放射線撮影方法 | |
JP7136033B2 (ja) | X線位相差撮影装置 | |
JP2013063098A (ja) | 放射線撮影装置及び画像処理方法 | |
JP6992902B2 (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: 19877812 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2020554766 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19877812 Country of ref document: EP Kind code of ref document: A1 |