PHASE CONTRAST X-RAY IMAGING SYSTEM

Abstract

A phase contrast X-ray imaging system includes an X-ray source; a plurality of gratings; a detector; a grating movement mechanism; and an image processor that generates a phase contrast image. The image processor generates the phase contrast image by using a pitch of an intensity change and a function which has the pitch as a variable and expresses the intensity change in a pixel value as a grating moves.

Claims

1-8. (canceled)

9. A phase contrast X-ray imaging system comprising: an X-ray source; a plurality of gratings including a first grating that is irradiated with an X-ray from the X-ray source, and a second grating that is irradiated with the X-ray from the first grating; a detector that detects the X-ray irradiated from the X-ray source; a grating movement mechanism that moves at least one of the plurality of gratings; and an image processor that generates a phase contrast image based on an intensity change indicating a change in a pixel value of each pixel detected by the detector, wherein the image processor generates the phase contrast image by using a pitch of the intensity change and a function which has the pitch as a variable and expresses the intensity change in the pixel value as the grating moves.

10. The phase contrast X-ray imaging system according to claim 9, wherein the pitch of the intensity change is a pitch including at least one shift of a shift in a pitch of the grating and a shift in a movement amount of the grating by the grating movement mechanism.

11. The phase contrast X-ray imaging system according to claim 9, wherein the image processor determines the pitch of the intensity change based on data on the intensity change which is acquired by moving at least one of the plurality of gratings, and the function.

12. The phase contrast X-ray imaging system according to claim 9, wherein the image processor determines the pitch of the intensity change by fitting a waveform shape of the intensity change and a waveform shape of the function to each other.

13. The phase contrast X-ray imaging system according to claim 12, wherein the image processor obtains a deviation between the waveform shape of the intensity change and the waveform shape of the function using, as the pitch, at least both of a value which is larger than a design value for a pitch of the grating moved by the grating movement mechanism and a value which is smaller than the design value.

14. The phase contrast X-ray imaging system according to claim 12, wherein the image processor acquires a coefficient of the function which corresponds to a predetermined pitch, and acquires a deviation between the waveform shape of the intensity change and the waveform shape of the function based on the coefficient and the predetermined pitch.

15. The phase contrast X-ray imaging system according to claim 9, wherein the image processor acquires a coefficient of the function from a pitch of the intensity change which is acquired in advance, and generates the phase contrast image based on the acquired coefficient and the pitch of the intensity change which is acquired in advance.

16. The phase contrast X-ray imaging system according to claim 9, wherein the plurality of gratings further include a third grating that is disposed between the X-ray source and the first grating.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0019] FIG. 1 is a schematic diagram of a phase contrast X-ray imaging system according to a first embodiment of the present invention when seen from a direction Y.

[0020] FIG. 2 is a schematic view for describing a process of acquiring an intensity change in the pixel value of a pixel.

[0021] FIG. 3(A) is a graph for describing a shift in the pitch of the intensity change in the pixel value of each pixel, and FIG. 3(B) is a schematic view for describing an artifact occurring in a phase contrast image to be generated.

[0022] FIGS. 4(A) to 4(C) are schematic views for describing a change in the artifact when the pitch is changed.

[0023] FIG. 5 is a flowchart for describing a method for generating a phase contrast image according to the first embodiment of the present invention.

[0024] FIG. 6 is a flowchart for describing a method for generating a phase contrast image according to a second embodiment of the present invention.

[0025] FIG. 7 is a schematic diagram of a phase contrast X-ray imaging system according to a modification example of the first embodiment of the present invention when seen from the direction Y.

DESCRIPTION OF EMBODIMENTS

[0026] Hereinafter, embodiments in which the present invention is materialized will be described based on the drawings.

First Embodiment

[0027] A configuration of a phase contrast X-ray imaging system 100 according to a first embodiment of the present invention will be described with reference to FIGS. 1 to 5.

[0028] (Configuration of Phase Contrast X-Ray Imaging System)

[0029] Firstly, the configuration of the phase contrast X-ray imaging system 100 according to the first embodiment of the present invention will be described with reference to FIG. 1.

[0030] FIG. 1 is a diagram of the phase contrast X-ray imaging system 100 when seen from a direction Y. As illustrated in FIG. 1, the phase contrast X-ray imaging system 100 includes an X-ray source 1; a first grating 2; a second grating 3; a detector 4; an image processor 5; a device control unit 9; and a grating movement mechanism 10. Incidentally, in this specification, a direction from the X-ray source 1 toward the first grating 2 is a direction Z2, and a direction reverse thereto is a direction Z1. In addition, a rightward and leftward direction in a plane perpendicular to a direction Z is the direction Y, a direction toward the back of the drawing sheet is a direction Y2, and a direction toward the front of the drawing sheet is a direction Y1. In addition, an upward and downward direction in the plane perpendicular to the direction Z is a direction X, an upward direction is a direction X1, and a downward direction is a direction X2.

[0031] The X-ray source 1 is configured to generate an X-ray and to irradiate the detector 4 (direction Z2) with the generated X-ray if a high voltage is applied to the X-ray source 1 based on a signal from the device control unit 9.

[0032] The first grating 2 has a plurality of slits 2a and X-ray phase change portions 2b which are arranged at a predetermined pitch pi in the direction X. Each of the slit 2a and the X-ray phase change portion 2b extends straightly. In addition, the slit 2a and the X-ray phase change portion 2b extend in parallel to each other. The first grating 2 is a so-called phase grating.

[0033] The first grating 2 is disposed between the X-ray source 1 and the second grating 3, and is irradiated with an X-ray from the X-ray source 1. The first grating 2 is provided so as to form a self-image 20 (refer to FIG. 2) of the first grating 2 using the Talbot effect. If an X-ray having coherence passes through a grating provided with slits, an image (self-image 20) of the grating is formed at a position which is a predetermined distance (Talbot distance) apart from the grating. This refers to the Talbot effect.

[0034] The second grating 3 has a plurality of radioparent portions 3a and X-ray absorbent portions 3b which are arranged at a predetermined pitch p.sub.2 in the direction X. Each of the radioparent portion 3a and the X-ray absorbent portion 3b extends straightly. In addition, the radioparent portion 3a and the X-ray absorbent portion 3b extend in parallel to each other. The second grating 3 is a so-called absorbent grating. The first grating 2 and the second grating 3 are gratings having different roles; however, each of the slit 2a and the radioparent portion 3a transmits an X-ray. In addition, the X-ray absorbent portion 3b plays a role in shielding an X-ray. The X-ray phase change portion 2b changes an X-ray phase due to a difference in refractive index between the slit 2a and the X-ray phase change portion 2b.

[0035] The second grating 3 is disposed between the first grating 2 and the detector 4, and is irradiated with an X-ray passing through the first grating 2. In addition, the second grating 3 is disposed at a position which is the Talbot distance apart from the first grating 2. The second grating 3 interferes with the self-image 20 of the first grating 2 to form moire fringes (not illustrated) on a front detection surface of the detector 4.

[0036] The detector 4 is configured to detect an X-ray, to transduce the detected X-ray into an electric signal, and to read the converted electric signal as an image signal. The detector 4 is, for example, a flat panel detector (FPD). The detector 4 is formed of a plurality of transduction elements (not illustrated) and pixel electrodes (not illustrated) which are disposed on the plurality of transduction elements. The plurality of transduction elements and the pixel electrodes are arranged in an array at a predetermined pitch (pixel pitch) in the direction X and the direction Y. In addition, the detector 4 is configured to output the acquired image signal to the image processor 5.

[0037] The image processor 5 includes an image generation unit 6; a control unit 7; and a storage unit 8. The image generation unit 6 is configured to generate an X-ray image (not illustrated) based on the image signal output from the detector 4. In addition, the image generation unit 6 is configured to generate a phase contrast image 11 (refer to FIG. 3(B)) based on the generated X-ray image. The image generation unit 6 includes a processor such as a graphics processing unit (GPU) or a field-programmable gate array (FPGA) configured for image processing.

[0038] In addition, the control unit 7 is configured to acquire a pitch (period) p of an intensity change 30 (refer to FIG. 2) of the X-ray detected by the detector 4, and the like. A detailed configuration where the control unit 7 acquires the pitch p of the intensity change 30 will be described later. In addition, the control unit 7 includes a processor such as a central processing unit (CPU).

[0039] The storage unit 8 is configured to store the X-ray image generated by the image generation unit 6, a program for generating the phase contrast image 11, and the like. The storage unit 8 includes a hard disk drive (HDD), a non-volatile memory, or the like.

[0040] The device control unit 9 is configured to stepwise move, via the grating movement mechanism 10, the second grating 3 in a direction (direction X) perpendicular to a direction (direction Y) of a grating pattern in a grating plane. In addition, the device control unit 9 includes a processor such as a CPU.

[0041] The grating movement mechanism 10 is configured to stepwise move the second grating 3 in the direction (direction X) perpendicular to the direction (direction Y) of the grating pattern of the grating in the grating plane (in a plane XY), based on a signal from the device control unit 9. Specifically, the grating movement mechanism 10 divides the pitch p.sub.2 of the second grating 3 into N steps, and stepwise moves the second grating 3 in a step of p.sub.2/N. Incidentally, N is a positive integer, and in the first embodiment, for example, N=4. In addition, the grating movement mechanism 10 includes a stepping motor, a piezo actuator, or the like.

[0042] (Configuration for Acquiring Intensity Change in Pixel Value of Each Pixel)

[0043] Subsequently, a configuration where the image processor 5 acquires the intensity change 30 in the pixel value of each pixel will be described with reference to FIG. 2.

[0044] FIGS. 2(A) to 2(D) are schematic views showing a positional relationship between the self-image 20 of the first grating 2 and the second grating 3 when the grating movement mechanism 10 translationally moves the second grating 3 and the image generation unit 6 according to the first embodiment generates the phase contrast image 11 by a fringe scanning method, and FIG. 2 is a graph showing the intensity change 30 in the pixel value of each pixel of the X-ray image captured while translationally moving the second grating 3.

[0045] In the first embodiment, the device control unit 9 is configured to capture images while translationally moving, via the grating movement mechanism 10, the second grating 3 four times in a step of p.sub.2/4. Namely, in the first embodiment, as illustrated in FIGS. 2(A) to 2(D), the phase contrast X-ray imaging system 100 is configured to capture images while moving the second grating 3 in a step of p.sub.2/4 in the direction X2. In addition, the control unit 7 acquires, as shown in the graph in FIG. 2, the intensity change 30 as an intensity change in the pixel value of each pixel when capturing images while moving the second grating 3 in a step of p.sub.2/4 in the direction X2.

[0046] (Description of Comparative Example)

[0047] Here, in a general fringe scanning method, the phase contrast image 11 is generated based on the intensity change 30 in an X-ray image captured in a state where an object Q is not disposed and the intensity change 30 in an X-ray image captured in a state where the object Q is disposed. For example, a case where images are captured while the second grating 3 is translated N times can be considered.

[0048] The intensity of an X-ray detected by the detector 4 when an image is captured in a state where the object Q is disposed is I.sub.k (x, y), the intensity of an X-ray when an image is captured in a state where the object Q is not disposed is I.sub.0k (x, y), and S (x, y) and S.sub.0 (x, y) are defined as follows.

[00001]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.1\right].Math.& \\ S\left(x,y\right)=\underset{k=1}{\overset{N}{.Math.}}.Math.I_k\left(x,y\right).Math.\exp \left(-\frac{2.Math.i.Math..Math.k}{N}\right)& \left(1\right)\\ S_0\left(x,y\right)=\underset{k=1}{\overset{N}{.Math.}}.Math.I_{0.Math.k}\left(x,y\right).Math.\exp \left(-\frac{2.Math.i.Math..Math.k}{N}\right)& \left(2\right)\end{array}$

[0049] Here, k is a positive integer from 1 to N. In addition, x and y are an x coordinate and a y coordinate of each pixel.

[0050] The pixel value of each of pixels representing a phase differential image 11a can be calculated by the following equation (3) using the equations (1) and (2). Incidentally, a pixel value (x, y) of each pixel of the phase differential image 11a is a difference in phase value between when the object Q is disposed and when the object Q is not disposed.

[00002]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.2\right].Math.& \\ {}_x\left(x,y\right)=\frac{p_2}{2.Math..Math.Z_T}.Math.\arg \left[\frac{S\left(x,y\right)}{S_0\left(x,y\right)}\right]& \left(3\right)\end{array}$

[0051] Here, Z.sub.T is a distance between the first grating 2 and the second grating 3. In addition, a pixel value T (x, y) of each pixel representing an absorption image 11b is expressed in the following equation (4).

[00003]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.3\right].Math.& \\ T\left(x,y\right)=\frac{\underset{k=1}{\overset{N}{.Math.}}.Math.I_k\left(x,y\right)}{\underset{k=1}{\overset{N}{.Math.}}.Math.I_{0.Math.k}\left(x,y\right)}& \left(4\right)\end{array}$

[0052] V (x, y) which is a visibility when the object Q is disposed and V.sub.0 (x, y) which is a visibility when the object Q is not disposed are expressed in the following equations (5) and (6).

[00004]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.4\right].Math.& \\ V\left(x,y\right)=\frac{2.Math..Math.S\left(x,y\right).Math.}{\underset{k=1}{\overset{N}{.Math.}}.Math.I_k\left(x,y\right)}& \left(5\right)\\ V_0\left(x,y\right)=\frac{2.Math..Math.S_0\left(x,y\right).Math.}{\underset{k=1}{\overset{N}{.Math.}}.Math.I_{0.Math.k}\left(x,y\right)}& \left(6\right)\end{array}$

[0053] Since a dark field image 11c is a visibility ratio between when the object Q is disposed and the object Q is not disposed, the dark field image 11c is expressed in the following equation (7).

[00005]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.5\right].Math.& \\ D\left(x,y\right)=\frac{V\left(x,y\right)}{V_0\left(x,y\right)}=\frac{.Math.S\left(x,y\right).Math.}{.Math.S_0\left(x,y\right).Math.}/T\left(x,y\right)& \left(7\right)\end{array}$

[0054] As shown in the equations (1) and (2), in the general fringe scanning method, the phase contrast image 11 is imaged based on a value calculated on the assumption that the second grating 3 is moved only by a value, which is obtained by dividing the pitch p.sub.2 of the second grating 3 by the number of steps, per step. In other words, since the equations (1) and (2) are on the assumption that an image to be acquired in each step is an image in each of phases , 2/4, , and 4/4 of the pitch p.sub.2 of the second grating 3, it is necessary to reliably acquire an image in each of the phases to 4/4 by capturing images while translationally moving the second grating 3. Therefore, if a shift occurs between the pitch p of the intensity change 30 to be acquired and the pitch p.sub.2 of the second grating 3, since a calculation is performed using an image (pixel value) in a phase which is shifted from each of the phases to 4/4, an artifact 12 is formed in the phase contrast image 11 to be generated.

[0055] (Shift in Pitch of Intensity Change)

[0056] Subsequently, the artifact 12 occurring in the phase contrast image 11 when a shift occurs between the pitch p of the intensity change 30 and the pitch p.sub.2 of the second grating 3 will be described with reference to FIG. 3.

[0057] The graphs illustrated in FIG. 3(A) schematically show an intensity change 30a (solid line graph) when the pitch p of the intensity change 30 coincides with the pitch p.sub.2 of the second grating 3, and an intensity change 30b (dotted line graph) when a shift occurs between the pitch p of the intensity change 30 and the pitch p.sub.2 of the second grating 3. Incidentally, the intensity change 30b is actually N points data; however, in the example shown in FIG. 3(A), each of the N points data is regarded as data for a pitch k/N. Here, in the phase contrast X-ray imaging system 100, the second grating 3 may be deformed due to heat generated from the X-ray source 1. If the second grating 3 is deformed, the pitch p.sub.2 of the second grating 3 changes. In addition, the grating movement mechanism 10 is configured to move the second grating 3 in parallel to the first grating 2; however, due to a positional shift during installation, or the like, a movement direction of the second grating 3 may not be parallel to the first grating 2. As described above, if the pitch p.sub.2 of the second grating 3 changes and/or a shift occurs in the movement direction of the second grating 3, a shift occurs between a total movement amount of the second grating 3 in the direction (direction X) perpendicular to the direction (direction Y) of the grating pattern and the pitch p.sub.2 of the second grating 3. If a shift occurs between the total movement amount of the second grating 3 in the direction (direction X) perpendicular to the direction (direction Y) of the grating pattern and the pitch p.sub.2 of the second grating 3, the pitch p of the intensity change 30 acquired by the control unit 7 is shifted from the pitch p.sub.2 of the second grating 3. Namely, the pitch p of the intensity change 30 is the pitch p including at least one shift of a shift in the pitch p.sub.2 of the second grating 3 and a shift in the movement amount of the second grating 3 by the grating movement mechanism 10.

[0058] In the general fringe scanning method, the phase contrast image 11 is generated, as shown in the equations (1) to (7), by applying the acquired intensity change 30b to a function which expresses the intensity change 30a having the same pitch as the pitch p.sub.2 of the second grating 3. Therefore, if in a state where a shift occurs between the pitch p of the intensity change 30 and the pitch p.sub.2 of the second grating 3, the phase contrast image 11 is generated by the general fringe scanning method using the equations (1) to (8), as illustrated in FIG. 3(B), the artifact 12 with a fringe shape is formed in each image.

[0059] Then, in the first embodiment, the image processor 5 is configured to generate the phase contrast image 11 by using the pitch p of the intensity change 30b and a function which has the pitch as a variable and expresses the intensity change 30b in pixel value as the grating moves. Specifically, the image generation unit 6 generates the phase contrast image 11 by performing the following process such that even in a state where a shift occurs between the pitch p of the intensity change 30 and the pitch p.sub.2 of the second grating 3, the artifact 12 can be prevented from occurring in the phase contrast image 11.

[0060] (Process of Generating Phase Contrast Image)

[0061] In the first embodiment, the control unit 7 is configured to determine the pitch p of the intensity change 30b based on data on the intensity change 30b acquired by moving the second grating 3, and the function. Specifically, the control unit 7 is configured to determine the pitch p of the intensity change 30b by fitting a waveform shape of the intensity change 30b and a waveform shape of the function to each other. The control unit 7 is configured to determine coefficients a to c of the function, and to acquire the pitch p where a sum Z of squares of deviation of the function becomes smallest by performing the least square fitting on N data points on the acquired intensity change 30b.

[0062] X coordinates x.sub.k of the N data points can be expressed in the following equation (8).

[00006]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.6\right].Math.& \\ x_k=p_2.Math.\frac{k}{N}& \left(8\right)\end{array}$

[0063] A function expressed by the following equation (9) is used in performing fitting.

[00007]$\begin{array}{cc}\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.7\right].Math.& \end{array}& \\ y=a.Math..Math.\sin .Math..Math.\frac{2.Math.}{p}.Math.x+b.Math..Math.\cos .Math..Math.\frac{2.Math.}{p}.Math.x+c& \left(9\right)\end{array}$

[0064] Here, p is the pitch of the intensity change 30b. In addition, a, b, and c are the coefficients of the function which are determined if the pitch p is determined.

[0065] The image processor 5 is configured to acquire the coefficients a to c of the equation (9) which correspond to a predetermined pitch. In addition, the control unit 7 is configured to acquire the sum Z of squares of deviation between the waveform shape of the intensity change 30b based on the acquired coefficients a to c and the predetermined pitch p, and a waveform shape of the equation (9). Specifically, the control unit 7 is configured to acquire the pitch p of the equation (9) by calculating the sum Z of squares of deviation between the function of the equation (9) and the N data points. The sum Z of squares of deviation can be expressed as shown in the following equation (10).

[00008]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.8\right].Math.& \\ Z=\underset{k=1}{\overset{N}{.Math.}}.Math.{\left(a.Math..Math.\sin .Math..Math.{\mathrm{mx}}_k+b.Math..Math.\cos .Math..Math.{\mathrm{mx}}_k+c-y_k\right)}^2& \left(10\right)\end{array}$

[0066] Here, m is a frequency of the intensity change 30b, and m=2/p.

[0067] In the first embodiment, the control unit 7 partially differentiates the coefficients a, b, and c of the equation (10) so as to perform fitting using the least square method. The partial differentiation of the coefficients a, b, and c is expressed in the following equations (11) to (13).

[00009]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.9\right].Math.& \\ \frac{Z}{a}=2.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left[\sin .Math.m.Math.x_k.Math.\left(a.Math..Math.\sin .Math.m.Math.x_k+b.Math..Math.\cos .Math..Math.{\mathrm{mx}}_k+c-y_k\right)\right]& \left(11\right)\\ \frac{Z}{b}=2.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left[\cos .Math..Math.{\mathrm{mx}}_k\left(a.Math..Math.\sin .Math..Math.{\mathrm{mx}}_k+b.Math..Math.\cos .Math..Math.{\mathrm{mx}}_k+c-y_k\right)\right]& \left(12\right)\\ \frac{Z}{c}=2.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left[a.Math..Math.\sin .Math..Math.{\mathrm{mx}}_k+b.Math..Math.\cos .Math..Math.{\mathrm{mx}}_k+c-y_k\right]& \left(13\right)\end{array}$

[0068] Here, y.sub.k is the pixel value of each pixel.

[0069] According to the definition of the least square method, since a, b, and c which cause each of the equations (11) to (13) to become zero may be obtained, the equations (11) to (13) are transformed into the following equations (14) to (16).

[00010]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.10\right]& \\ a.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\sin .Math..Math.{\mathrm{mx}}_k.Math..Math.\sin .Math..Math.{\mathrm{mx}}_k\right)+b.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\sin .Math..Math.{\mathrm{mx}}_k.Math.\cos .Math..Math.{\mathrm{mx}}_k\right)+c.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\sin .Math..Math.{\mathrm{mx}}_k\right)=\underset{k=1}{\overset{N}{.Math.}}.Math.\left(y_k.Math.\sin .Math..Math.{\mathrm{mx}}_k\right)& \left(14\right)\\ a.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\sin .Math..Math.{\mathrm{mx}}_k.Math.\cos .Math..Math.{\mathrm{mx}}_k\right)+b.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\cos .Math..Math.{\mathrm{mx}}_k.Math.\cos .Math..Math.{\mathrm{mx}}_k\right)+c.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\cos .Math..Math.{\mathrm{mx}}_k\right)=\underset{k=1}{\overset{N}{.Math.}}.Math.\left(y_k.Math.\cos .Math..Math.{\mathrm{mx}}_k\right)& \left(15\right)\\ .Math.a.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\sin .Math..Math.{\mathrm{mx}}_k\right)+b.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(\cos .Math..Math.{\mathrm{mx}}_k\right)+c.Math.\underset{k=1}{\overset{N}{.Math.}}.Math.\left(1\right)=\underset{k=1}{\overset{N}{.Math.}}.Math.\left(y_k\right)& \left(16\right)\end{array}$

[0070] The coefficients a, b, and c can be expressed in the following equations (17) to (19) by solving simultaneous equations including the equations (14) to (16).

[00011]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.11\right].Math.& \\ a=\frac{a_1}{d}& \left(17\right)\\ b=\frac{b_1}{d}& \left(18\right)\\ c=\frac{c_1}{d}& \left(19\right)\end{array}$

[0071] At the time, each of a.sub.1, b.sub.1, c.sub.1, and d can be expressed in the following equations (20) to (23).

[00012]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.12\right]& \\ a_1=\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}\right)+\left(.Math.y\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)+\left(.Math.1\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.1\right).Math.(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right)-\left(.Math.y\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.& \left(20\right)\\ b_1=\left(.Math.1\right).Math.\left(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\sin .Math..Math.\mathrm{mx}\right)+\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right)+\left(.Math.y\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.y\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\sin .Math..Math.\mathrm{mx}\right)-\left(.Math.1\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.& \left(21\right)\\ c_1=\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}\right)+\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right)+\left(.Math.y\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.y\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-.Math.\sin .Math..Math.\mathrm{mx}).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.y.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.& \left(22\right)\\ d=\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)+\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)+\left(.Math.1\right).Math.\left(.Math.\sin .Math..Math.m.Math..Math.x.Math..Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)-\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\sin .Math..Math.\mathrm{mx}\right)-\left(.Math.1\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right).Math.\left(.Math.\sin .Math..Math.\mathrm{mx}.Math..Math.\cos .Math..Math.\mathrm{mx}\right)& \left(23\right)\end{array}$

[0072] The coefficients a, b, and c acquired from the equations (17) to (19) and the pitch p (frequency m) included in the equation (9) are unknown numbers. Therefore, in the first embodiment, the image processor 5 is configured to acquire the pitch p where the sum Z of squares of deviation becomes smallest. In the example shown in FIG. 3(A), the pitch p of the intensity change 30b is shifted to become smaller than the pitch p.sub.2 of the second grating 3. For this reason, the control unit 7 may acquire the sum Z of squares of deviation using a plurality of the pitches p, the values of which are smaller than the pitch p.sub.2 of the second grating 3. However, actually, the pitch p of the intensity change 30b is an unknown number. For this reason, the control unit 7 is configured to obtain the sum Z of squares of deviation between the waveform shape of the intensity change 30b and the waveform shape of the equation (9) using, as the pitch, at least both of a value which is larger than a design value for the pitch p.sub.2 of the second grating 3 moved by the grating movement mechanism 10 and a value which is smaller than the design value.

[0073] For example, the control unit 7 acquires a first deviation using the value larger than the pitch p.sub.2 of the second grating 3 as a pitch, and a second deviation using the value smaller than the pitch p.sub.2 of the second grating 3 as a pitch. The control unit 7 compares the acquired first deviation with the acquired second deviation. If the first deviation is smaller than the second deviation, the control unit 7 acquires a plurality of the sums Z of squares of deviation using a plurality of the pitches p larger than the pitch p.sub.2 of the second grating 3. On the other hand, if the second deviation is smaller than the first deviation, the control unit 7 acquires a plurality of the sums Z of squares of deviation using a plurality of the pitches p smaller than the pitch p.sub.2 of the second grating 3. The control unit 7 is configured to determine, as the pitch p of the intensity change 10b, the pitch p where the acquired plurality of sums Z of squares of deviation become smallest.

[0074] FIG. 4 shows an example where the phase differential image 11a, the absorption image 11b, and the dark field image 11c are generated by using the plurality of pitches, namely, the pitches p of the intensity change 30b which are smaller than the pitch p.sub.2 of the second grating 3. An analysis pitch in FIG. 4 is a pitch obtained by multiplying the pitch p.sub.2 of the second grating 3 by a predetermined coefficient. For example, an analysis pitch of 0.95 is a pitch (0.95p.sub.2) obtained by multiplying the pitch p.sub.2 of the second grating 3 by 0.95. It can be seen that as illustrated in FIG. 4, as the analysis pitch decreases from 1.00 to 0.95, the artifact 12 occurring in each image decreases. In addition, it can also be seen that if the analysis pitch further decreases from 0.95, the artifact 12 increases. Therefore, in the example shown in FIG. 4, a pitch where the sum Z of squares of deviation becomes smallest is an analysis pitch of 0.95 (0.95p.sub.2). Incidentally, the example shown in FIG. 4 is an example where a simulation is performed using 0.95p.sub.2 as the pitch p of the intensity change 30b. Therefore, actually, in order to determine the pitch p of the intensity change 30b, the control unit 7 is configured to set a plurality of the pitches p including the point of inflection at which the pitch p where the magnitude of the sum Z of squares of deviation decreases changes to the pitch p where the magnitude increases, and to determine, as the pitch p of the intensity change 30b, the pitch p where the sum Z of squares of deviation becomes the minimum value.

[0075] In the first embodiment, the image generation unit 6 calculates each pixel value of the phase contrast image 11 using the coefficients a, b, and c calculated from the equations (17) to (23) and the determined pitch p. The following equation (24) rewritten from the equation (9) is used when calculating each pixel value of the phase contrast image 11.

[00013]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.13\right].Math.& \\ y=A.Math.\sin \left(\frac{2.Math.}{p}.Math.x+B\right)+C& \left(24\right)\end{array}$

[0076] Here, there is a relationship between coefficients A, B, and C and the coefficients a, b, and c which is expressed in the following equations (25) to (27). Incidentally, the unit in the following equation (26) is rad.

[00014]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.14\right].Math.& \\ A=\sqrt{a^2+b^2}& \left(25\right)\\ B={\tan }^{-1}\left(\frac{b}{a}\right)& \left(26\right)\\ C=c& \left(27\right)\end{array}$

[0077] If the coefficients A, B, and C of the equation (24) are used, the phase differential image 11a, the absorption image 11b, and the dark field image 11c can be expressed in the following equations (28) to (30).

[00015]$\begin{array}{cc}\left[\mathrm{Numeral}.Math..Math.15\right].Math.& \\ =\frac{p_2}{2.Math..Math.Z_T}.Math.\left(B_s-B_r\right)& \left(28\right)\\ T=\frac{C_s}{C_r}& \left(29\right)\\ D=\frac{A_s/C_s}{A_r/C_r}& \left(30\right)\end{array}$

[0078] Here, is the phase differential image. In addition, T is the absorption image. In addition, D is the dark field image. In addition, subscripts s and r are a coefficient when an image is captured in a state where the object Q is present, and a coefficient when an image is captured in a state where the object Q is not present.

[0079] As described above, in the first embodiment, using the least square method, the control unit 7 is configured to fit the waveform shape of the intensity change 30b of the actually captured X-ray image and the waveform shape of the equation (9) to each other. Therefore, even though the pitch p of the intensity change 30b is shifted from the pitch p.sub.2 of the second grating 3, it is possible to generate the phase contrast image 11 by obtaining the function (equation (9)) corresponding to the pitch p of the intensity change 30b.

[0080] (Method for Generating Phase Contrast Image)

[0081] Subsequently, a method for generating the phase contrast image 11 according to the first embodiment will be described with reference to FIG. 5.

[0082] In Step S1, the image generation unit 6 acquires an X-ray image which is a captured image of the object Q. Thereafter, the control unit 7 acquires data on the intensity change 30b in the pixel value of each pixel of the X-ray image which is a captured image of the object Q. Subsequently, in Step S2, the control unit 7 acquires a sum of squares of the first deviation by substituting a first pitch having a value, which is larger than the pitch p.sub.2 of the second grating 3, into the function. In addition, the control unit 7 acquires a sum of squares of the second deviation by substituting a second pitch having a value, which is smaller than the pitch p.sub.2 of the second grating 3, into the function. Incidentally, in Step S2, in order to acquire the sum of squares of the first deviation and the sum of squares of the second deviation, an equation including only the pitch p as a variable, which is obtained by solving the equation (9) using the equation (9), the equations (17) to (23), and a relational expression of m=2/p, is used. Thereafter, the process proceeds to Step S3.

[0083] In Step S3, the control unit 7 compares the magnitude of the value of the sum of square of the first deviation with the magnitude of the value of the sum of squares of the second deviation. If the sum of squares of the first deviation is smaller than the sum of squares of the second deviation, the process proceeds to Step S4. If the sum of squares of the first deviation is larger than the sum of squares of the second deviation, the process proceeds to Step S5.

[0084] In Step S4, the control unit 7 sets a plurality of the pitches p, the values of which are smaller than the pitch p.sub.2 of the second grating 3. In addition, in Step S5, the control unit 7 sets a plurality of the pitches p, the values of which are larger than the pitch p.sub.2 of the second grating 3. Thereafter, the process proceeds to Step S6.

[0085] In Step S6, the control unit 7 acquires a plurality of the sums Z of squares of deviation using the plurality of pitches p which are set in Step S4 (Step S5). Thereafter, the control unit 7 acquires the pitch p where the plurality of sums Z of squares of deviation become smallest, and sets the pitch p as the pitch p of the function (equation (9)). Thereafter, in Step S7, the control unit 7 acquires the coefficients a to c of the function (equation (9)). Thereafter, in Step S8, the image generation unit 6 generates the phase contrast image 11 by using the function (equation (9)) where the pitch p acquired in Step S6 is determined, and the coefficients a to c acquired in Step S7, and ends the process.

[0086] (Effects of First Embodiment)

[0087] In the first embodiment, it is possible to obtain the following effects.

[0088] In the first embodiment, as described above, the phase contrast X-ray imaging system 100 includes the X-ray source 1; a plurality of gratings including the first grating 2 irradiated with an X-ray from the X-ray source 1 and the second grating 3 irradiated with the X-ray from the first grating 2; the detector 4 that detects the X-ray irradiated from the X-ray source 1; the grating movement mechanism 10 that moves the second grating 3; and the image processor 5 that generates the phase contrast image 11 based on the intensity change 30b indicating a change in the pixel value of each pixel detected by the detector 4. The image processor 5 is configured to generate the phase contrast image 11 by using the pitch p of the intensity change 30b and the function (equation (9)) which has the pitch as a variable and expresses the intensity change 30b in pixel value as the grating moves. Therefore, since the pitch p of the function (equation (9)) is a variable, it is possible to approximate the pitch p of the intensity change 30b, which is actually acquired while moving the second grating 3, by obtaining an optimal solution for the pitch p of the function (equation (9)). As a result, even though a shift occurs between the pitch p of the intensity change 30b in the pixel value of each pixel and the pitch p.sub.2 of the second grating 3, since it is possible to adapt the pitch p of the function (equation (9)) to the pitch p of the intensity change 30b, it is possible to prevent the artifact 12 from occurring in the phase contrast image 11 to be obtained.

[0089] In addition, in the first embodiment, as described above, the pitch p of the intensity change 30b is a pitch including at least one shift of a shift in the pitch p.sub.2 of the second grating 3 and a shift in the movement amount of the second grating 3 by the grating movement mechanism 10. Therefore, even though a shift occurs in either of the pitch p.sub.2 of the second grating 3 and the movement amount of the second grating 3 (or, even though a shift occurs in both), it is possible to prevent the artifact 12 from occurring in the phase contrast image 11 to be generated by adapting the pitch p (variable) of the function (equation (9)) to the pitch p of the intensity change 30b which includes the shift.

[0090] In addition, in the first embodiment, as described above, the image processor 5 is configured to determine the pitch p of the intensity change 30b based on the data on the intensity change 30b acquired while moving the second grating 3, and the function (equation (9)). Therefore, it is possible to acquire the pitch p of the intensity change 30b from the actually measured data on the intensity change 30b. As a result, even though a shift occurs between the pitch p of the intensity change 30b and the pitch p.sub.2 of the second grating 3, it is possible to acquire the pitch p of the intensity change 30b which includes the shift.

[0091] In addition, in the first embodiment, as described above, the image processor 5 is configured to determine the pitch p of the intensity change 30b by fitting the waveform shape of the intensity change 30b and the waveform shape of the function (equation (9)) to each other. Therefore, it is possible to determine the pitch p of the intensity change 30b based on the sum Z of squares of deviation between the waveform shape of the intensity change 30b and the waveform shape of the function (equation (9)). As a result, since it is possible to determine the pitch p of the intensity change 30b by comparing together the sums Z of squares of deviation acquired using a plurality of the pitches, for example, compared to the case where the pitch p of the intensity change 30b is determined by performing image processing on images acquired using a plurality of the pitches p, it is possible to further simplify the process of determining the pitch p of the intensity change 30b.

[0092] In addition, in the first embodiment, as described above, the sum Z of squares of deviation between the waveform shape of the intensity change 30b and the waveform shape of the function (equation (9)) is obtained using, as the pitch p, at least both of a value which is larger than the design value for the pitch p.sub.2 of the second grating 3 moved by the grating movement mechanism 10 and a value which is smaller than the design value. Therefore, whether the pitch p where the sum Z of squares of deviation becomes smallest is a value larger than the design value for the pitch p.sub.2 of the second grating 3 or a value smaller than the design value can be determined by comparing the first deviation calculated using the value larger than the design value for the pitch p.sub.2 of the second grating 3 with the second deviation calculated using the value smaller than the design value for the pitch p.sub.2 of the second grating 3. As a result, it is possible to easily determine the pitch p where the sum Z of squares of deviation becomes smallest.

[0093] In addition, in the first embodiment, as described above, the image processor 5 is configured to acquire the coefficients a to c of the function (equation (9)) which correspond to the predetermined pitch p, and to acquire the sum Z of squares of deviation between the waveform shape of the intensity change 30b and the waveform shape of the function (equation (9)) based on the coefficients a to c and the predetermined pitch p. Therefore, it is possible to acquire the sum Z of squares of deviation which corresponds to the predetermined pitch p. As a result, it is possible to easily determine the pitch p where the sum Z of squares of deviation becomes smallest by comparing together the sums Z of squares of deviation which correspond to a plurality of the pitches p.

Second Embodiment

[0094] Subsequently, a phase contrast X-ray imaging system 200 according to a second embodiment of the present invention will be described with reference to FIGS. 1 and 6. Unlike the first embodiment where the pitch p of the function is determined by comparing together a plurality of the sums Z of squares of deviation acquired using a plurality of the pitches p, in the second embodiment, the image processor 5 is configured to acquire the coefficients a to c of the function (equation (9)) from the pitch p of the intensity change 30b which is acquired in advance, and to generate the phase contrast image 11 based on the acquired coefficients a to c and the pitch p of the intensity change 30b which is acquired in advance. Incidentally, the same reference signs are assigned to the same configurations as those in the first embodiment, and the descriptions thereof are omitted.

[0095] Here, since the pitch p.sub.2 of the second grating 3 is a very small value, the pitch p.sub.2 may be slightly shifted due to manufacturing errors. In this case, even though a pitch (hereinafter, referred to as p.sub.x) which is slightly shifted from the pitch p.sub.2 of the second grating 3 is already known, in the general fringe scanning method using the equations (1) and (2), it is necessary to acquire data for the pitch k/N. However, if the accuracy of the grating movement mechanism 10 is not sufficient, since it is not possible to translationally move the second grating 3 by a movement amount (p.sub.x(k/N)) obtained by dividing the already-known pitch p.sub.x by N, it is not possible to acquire the data for the pitch k/N, and the artifact 12 occurs. Then, in the second embodiment, the control unit 7 is configured to premeasure a shift in the pitch p.sub.2 of the second grating 3. In addition, the control unit 7 is configured to determine the coefficients a to c by fitting an actual measurement value and the equation (9), into which the pitch p.sub.2 of the intensity change 30b which is acquired in advance is substituted, to each other.

[0096] (Method for Generating Phase Contrast Image)

[0097] Subsequently, a method for generating the phase contrast image 11 according to the second embodiment will be described with reference to FIG. 6.

[0098] In Step S1, the control unit 7 acquires data on the intensity change 30b in the pixel value of each pixel of an X-ray image which is a captured image of the object Q. Subsequently, in Step S10, the control unit 7 acquires the pitch p of the intensity change 30b. Thereafter, in Step S11, the control unit 7 acquires the coefficients a to c of the equation (9) from the pitch p of the intensity change 30b which is acquired in advance. Thereafter, the process proceeds to Step S8.

[0099] In Step S8, the image generation unit 6 generates the phase contrast image 11 based on the acquired coefficients a to c and the pitch p of the intensity change 30b which is acquired in advance.

[0100] Incidentally, other configurations of the second embodiment are the same as those of the first embodiment.

[0101] (Effects of Second embodiment)

[0102] In the second embodiment, it is possible to obtain the following effects.

[0103] In the second embodiment, as described above, the image processor 5 is configured to acquire the coefficients a to c of the function (equation (9)) from the pitch p of the intensity change 30b which is acquired in advance, and to generate the phase contrast image 11 based on the acquired coefficients a to c and the pitch p of the intensity change 30b which is acquired in advance. Therefore, it is possible to acquire the coefficients a to c of the function (equation (9)) by fitting an actual measurement value and the function (equation (9)), into which the pitch p of the intensity change 30b which is acquired in advance is substituted, to each other. As a result, even though the movement amount of the second grating 3 is shifted, it is possible to determine the function (equation (9)) which accurately approximates actual data, and it is possible to prevent the artifact 12 from occurring in the phase contrast image 11 to be generated.

[0104] Incidentally, other effects of the second embodiment are the same as those of the first embodiment.

[0105] (Modification Example)

[0106] Incidentally, it should be considered that the embodiments disclosed this time are examples in all aspects and the present invention is not limited thereto. The scope of the present invention is not determined by the descriptions of the embodiments, but is as described in the claims, and further includes meanings equivalent to the claims and all changes (modification examples) within the scope.

[0107] For example, the first to third embodiments show an example where the first grating 2 and the second grating 3 are provided as a plurality of gratings; however, the present invention is not limited thereto. For example, as in a phase contrast X-ray imaging system 300 illustrated in FIG. 7, a third grating 40 may be provided between the X-ray source 1 and the first grating 2. The third grating 40 has a plurality of slits 40a and radioparent portions 40b which are arranged at a predetermined pitch p.sub.3 in the direction X. Each of the slit 40a and the radioparent portion 40b extends straightly. In addition, the slit 40a and the radioparent portion 40b extend in parallel to each other. In addition, the third grating 40 is disposed between the X-ray source 1 and the first grating 2, and is irradiated with an X-ray from the X-ray source 1. The third grating 40 is configured such that an X-ray passing through each of the slits 40a serves as a line light source corresponding to the position of each of the slits 40a. Therefore, it is possible to enhance the coherence of the X-ray irradiated from the X-ray source 1 due to the third grating 40. As a result, since it is possible to form the self-image 20 of the first grating 2 without depending on a focal diameter of the X-ray source 1, it is possible to improve the degree of freedom in selecting the X-ray source 1.

[0108] In addition, the first and second embodiments show an example where a phase grating is used as the first grating 2; however, the present invention is not limited thereto. For example, an absorbent grating may be used as the first grating 2. If an absorbent grating is used as the first grating 2, the image processor 5 (image generation unit 6) generates the phase contrast image 11 by using a fringe pattern of an X-ray transmitting through the first grating 2, and the second grating 3. Therefore, since it is possible to acquire the phase contrast image 11 without using the self-image 20 of the first grating 2, it is possible to improve the degree of freedom in the position of disposition of the first grating 2. However, if an absorbent grating is used as the first grating 2, the image quality of the phase contrast image 11 to be obtained decreases, and thus if the phase contrast image 11 with a high image quality is desired to be obtained, a phase grating is preferably used as the first grating 2.

[0109] In addition, the first embodiment shows an example where values smaller than the pitch p.sub.2 of the second grating 3 are used as a plurality of the pitches p which are used to acquire the sum Z of squares of deviation; however, the present invention is not limited thereto. For example, if the pitch p of the intensity change 30b is larger than the pitch p.sub.2 of the second grating 3, values larger than the pitch p.sub.2 of the second grating 3 may be used as the plurality of pitches p.

[0110] In addition, the first embodiment shows an example where the least square method is used to fit the waveform shape of the intensity change 30b and the waveform shape of the function to each other; however, the present invention is not limited thereto. Any fitting method may be used as long as it is possible to acquire the coefficients a to c of the function (equation (9)).

[0111] In addition, the first and second embodiments show an example where the grating movement mechanism 10 translationally moves the second grating 3; however, the present invention is not limited thereto. The first grating 2 may be translationally moved. A grating to be moved may be any grating.

[0112] In addition, the first and second embodiments show an example where the control unit 7 determines the sum Z of squares of deviation of the function (equation (9)); however, the present invention is not limited thereto. For example, a display unit may be provided to display, as illustrated in FIG. 4, a plurality of the phase contrast images 11 where the pitch p is changed, and the sum Z of squares of deviation of the function (equation (9)) may be determined by user's selection of the phase contrast images 11 displayed on the display unit.

REFERENCE SIGNS LIST

[0113] 1: X-ray source

[0114] 2: first grating

[0115] 3: second grating

[0116] 4: detector

[0117] 5: image processor

[0118] 10: grating movement mechanism

[0119] 12: phase contrast image

[0120] 30, 30a, 30b: intensity change

[0121] 40: third grating

[0122] 100, 200, 300: phase contrast X-ray imaging system

[0123] Q: object

[0124] p: pitch of intensity change

[0125] p.sub.2: pitch of grating

[0126] Z: sum of squares of deviation