The present invention relates to an X-ray imaging apparatus.
An X-ray has high material transparency, and can achieve imaging with high spatial resolution. For these reasons, the X-ray is used for a nondestructive inspection of an object as industrial utilization, X-raying as medical utilization, and the like.
That is, by the X-ray in the above utilization, a contrast image is formed by using a difference of absorption in a case where the X-ray transmits through an object or a living body, due to constituent elements and density differences of the object or the living body. It should be noted that such a process is called an X-ray absorption contrast method. However, since an X-ray absorption capability of a light element is very small, it is difficult by the X-ray absorption contrast method to take an image of living soft tissue which consists of carbon, hydrogen, oxygen and the like being constituent elements of the living body, or a soft material.
On the other hand, in order to provide a method which can clearly take an image of even tissue consisting of light elements, a research for a phase contrast method using a phase difference of X-rays has been performed since 1990's.
Among numerous phase contrast methods developed, a method capable of obtaining an image by a single X-ray imaging is desirable for somatoscopy, and particularly for medical somatoscopy. For example, a method described in PTL 1 and a method described in NPL 1 have been proposed.
More specifically, in the method described in PTL 1, an X-ray is divided by periodically arranged slits so that the divided X-ray corresponds to a pixel of a detector. Further, an X-ray shielding mask is arranged at the edge of the pixel of the detector, and slight deviation in an X-ray direction which arises when the X-ray transmits through an object to be detected is detected as an output change of the pixel, whereby a phase image is obtained from a first imaging.
Further, in a Scanning Double Grating method (hereinafter, called an SDG method) in NPL 1, imaging is performed while a diffraction grating which forms a self-image and an absorption grating which is disposed at the position of the self-image are simultaneously scanned in the same direction. Thus, an inclination of wavefront of an X-ray is detected as an intensity change through a single imaging, whereby a phase image is obtained.
In each of the methods described in PTL 1 and NPL 1, in a case where absorption of the X-rays existing in the object to be detected is non-uniformly distributed, the non-uniformly distributed absorption becomes an error factor, whereby quality of the obtained phase image deteriorates.
In consideration of such a problem, the present invention aims to provide an X-ray imaging apparatus which can reduce an influence of a transmissivity distribution of the X-ray in the object to be detected even in the case where there is the non-uniformly distributed absorption in the object to be detected, and thus take the phase image having exceptional quality.
In one aspect of the present invention, an X-ray imaging apparatus, which takes an image of an object to be detected, is characterized by comprising; an X-ray source, a first grating configured to form a periodic bright-dark pattern by a Talbot effect, on the basis of an X-ray from the X-ray source, a second grating, disposed at a position where the bright-dark pattern is formed, configured to block a part of the bright-dark pattern, a detector configured to detect an X-ray intensity distribution of the X-ray which passed through the second grating; and a calculator configured to calculate phase information of the X-ray on the basis of the X-ray intensity distribution detected by the detector, wherein the second grating includes a first region in which a first blocking pattern has been formed and a second region in which a second blocking pattern has been formed, and a direction in which the first blocking pattern blocks a bright section of the bright-dark pattern and a direction in which the second blocking pattern blocks the bright section of the bright-dark pattern are different from each other. Other aspects of the present invention will be clarified in the following exemplary embodiments of the present invention.
Further features of the present invention will become apparent from the following description of exemplary embodiments with reference to the attached drawings.
Hereinafter, a constructive example of an X-ray imaging apparatus which takes a phase image by a Talbot effect using an X-ray, according to an exemplary embodiment of the present invention will be described.
The X-ray imaging apparatus according to the present embodiment comprises:
an X-ray source which radiates an X-ray;
a first grating which is disposed at the position to which the X-ray is incident from the X-ray source and forms a periodic bright-dark pattern by a Talbot effect using the transmitting X-ray;
a second grating which is disposed at the position where the bright-dark pattern is formed by the first grating, and blocks a part of the bright-dark pattern; and
an X-ray intensity detector which detects X-ray intensity information based on an X-ray intensity distribution of the X-ray transmitted through an object to be detected disposed between the first grating and the second grating or between the first grating and the X-ray source.
Thus, the X-ray imaging apparatus obtains, by a calculator, phase information of the X-ray from the X-ray intensity information detected by the X-ray intensity detector.
In the X-ray imaging apparatus thus constructed, the X-ray necessary for measurement is radiated by the X-ray source toward the object to be detected.
The object to be detected is disposed between the first grating and the second grating or between the first grating and the X-ray source, and the periodic bright-dark pattern is formed by the Talbot effect using the X-ray transmitting through the first grating.
The first grating is constructed by an absorption member on which openings are periodically arranged. Concretely, the first grating is constructed by, for example, a phase grating composed of an X-ray transmitting member of which the thickness periodically changes, or a slit which has periodically arranged openings.
The second grating is disposed so as to block a part of the bright section of the bright-dark pattern formed by the first grating.
That is, in order to block a part of the bright-dark pattern, the second grating has a first region in which a first blocking pattern is formed and a second region in which a second blocking pattern is formed, as the regions on which the blocking patterns extending mutually opposite directions respectively are formed in regard to the periodicity direction of the bright-dark pattern.
Here, concretely, the respective blocking patterns block, for example, right and left halves, upper and lower halves, or oblique-45° halves of the pattern.
Further, with respect to the bright-dark patterns which are incident to the same pixel of the X-ray intensity detector, the blocking directions thereof are set to be the same.
Therefore, if the bright-dark pattern moves in the direction perpendicular to the blocking direction, an amount of received light of the corresponding pixel changes.
If an inclination of a wavefront changes according to transmission of the X-ray through the object to be detected, an X-ray direction changes, and the bright-dark pattern on the second grating thus moves.
Therefore, an inclination of a transmitted wavefront (hereinafter, called a differential wavefront) of the object to be detected can be obtained based on an output change of the relevant pixel. However, since the object to be detected generally includes a uniform transmissivity distribution caused by absorption of the X-ray, the amount of received light of the relevant pixel also changes by the uniform transmissivity distribution.
In the present embodiment, in order to reduce an influence by the transmissivity distribution in the object to be detected, the second gratings respectively provided on the adjacent pixels are set so that the blocking directions thereof are mutually opposite to each other.
In the output changes of these two pixels, the changes caused by the absorption occur in phase. However, since the changes caused by the movement of the bright-dark pattern occur in reversed phases, the influence of the transmissivity distribution can be reduced by separating these changes, whereby the accurate differential wavefront can be obtained.
In the calculator of the X-ray imaging apparatus according to the present embodiment, calculation is performed as follows.
Namely, a ratio of a difference between the outputs of a pair of the pixels respectively blocked in the opposite directions and a sum of these outputs is obtained, and the obtained ratio is multiplied by a proportionality constant determined a distance between the first grating and the second grating and a period of the first grating, whereby the inclination of the wavefront at a relevant portion is obtained. Then, the above-described calculation is performed to the overall effective region of the X-ray intensity detector, whereby the differential wavefront is calculated.
Also, the calculator can calculate the transmitted wavefront of the object to be detected by integrating the differential wavefronts concerning the perpendicular two directions.
Hereinafter, examples of the present invention will be described.
As an example 1, a construction of the X-ray imaging apparatus to which the present invention is applied will be described with reference to
This drawing illustrates an X-ray source 10 which radiates an X-ray, an X-ray 11 which is radiated by the X-ray source 10, and an object 20 to be detected which is to be measured by the X-ray imaging apparatus. Further, this drawing illustrates a first grating 30 which is used to form a periodic bright-dark pattern by the transmitting X-ray, a bright section 40 of the bright-dark pattern formed by the first grating 30, and a second grating 50 which blocks a part of the bright section 40 of the bright-dark pattern.
Furthermore, this drawing illustrates an X-ray intensity detector 70 which takes and detects an X-ray image received, and pixels Pi (i=1, 2, 3, . . . ) of the X-ray intensity detector 70.
Furthermore, this drawing illustrates a calculator 80 which calculates a differential wavefront and a transmitted wavefront from the image which has been taken by the X-ray intensity detector 70.
In the above construction, the X-ray 11 radiated from the X-ray source 10 passes through the object 20 to be detected and the first grating 30, and then reaches the second grating 50.
The first grating 30 is manufactured by silicone of which X-ray transmissivity is large and which has excellent workability.
A phase modulation grating of a phase difference π or π/2 (hereinafter, called a one-dimensional π phase grating or a one-dimensional π/2 phase grating) that a portion where a phase of the transmitting X-ray is relatively different by π or π/2 is one-dimensionally and periodically arranged is formed by periodically changing the thickness of the silicone.
That is, the thickness of a hatched portion 32 is made different from the thickness of a non-hatched portion 31 so that a transmission phase of the hatched portion 32 is different from a transmission phase of the non-hatched portion 31 by n or n/2.
The X-ray transmitting section and the blocking section are periodically arranged on the second grating 50, and a blocking section 51 is manufactured by gold which is a member for satisfactorily blocking the X-ray.
The second grating 50 is disposed at a position where the Talbot effect by the first grating 30 arises based on the X-ray radiated from the X-ray source 10 (hereinafter, called a Talbot position).
That is, if it is assumed that a distance between the X-ray source and the first grating is Z0 and a distance between the first grating 30 and the second grating 50 is Z1, the second grating is disposed so as to satisfy the following expression (1).
In the expression (1), λ indicates a wavelength of the X-ray, d indicates a grating period of the first grating 30, n indicates a natural number, and N indicates a real number which is expressed as n/4−⅛ in case of the one-dimensional π phase grating and as n− 1/2 in case of the one-dimensional π/2 phase grating.
1/Z0+1/Z1=1/N×λ/d2 (1)
The periodic bright-dark pattern is generated by the Talbot effect on the second grating 50.
The blocking section 51 of the second grating 50 is disposed at the position of blocking a part of the bright section 40 of the bright-dark pattern.
When passing through the object 20 to be detected, a wavefront 12 of the X-ray 11 deforms according to its refractive-index distribution, whereby a transmitted wavefront 13 is obtained.
The bright section 40 moves in proportion to an inclination change of the corresponding section of the transmitted wavefront 13. An amount of light of the X-ray which reaches the X-ray intensity detector 70 through the second grating 50 changes according to an amount of the movement.
The size of the pixel is set to be equal to an integral multiple of the period of the bright-dark pattern so as to make the number of the bright sections 40 for light reception of each pixel constant for all the pixels.
If the intensity of the bright section 40 radiated to the second grating 50 is identical, an inclination of the transmitted wavefront 13 can be calculated because it is proportional to a change amount of the amount of received light.
However, since the intensity of the bright section 40 of the bright-dark pattern depends on X-ray transmissivity of the object 20 to be detected, an amount of movement of the bright section 40 cannot be calculated only by the change of the amount of received light of each pixel if the X-ray transmissivity is uneven.
For this reason, in this example, in order to be able to calculate the amount of movement of the bright section 40 even if the X-ray transmissivity is uneven, the direction in which the blocking section 51 of the second grating 50 blocks the bright section 40 is set to be identical on the same pixel of the X-ray intensity detector 70. On the other hand, the direction in which the blocking section 51 of the second grating 50 blocks the bright section 40 on one pixel is set to be opposite to the direction in which the blocking section 51 of the second grating 50 blocks the bright section 40 on the adjacent pixel.
That is, the blocking by the first blocking pattern of the first region on the pixel P2n−1 aims at the left half, while the blocking by the second blocking pattern of the second region on the pixel P2n aims at the right half.
Here, it is assumed that n=1, 2, 3, . . . . Therefore, for example, the amount of received light in a case where the bright section 40 moves in the x direction being the direction perpendicular to the gratings (lines) of the first grating 30 due to an inclination of the wavefront increases in the pixels P1 and P3, but decreases in the pixels P2 and P4. Further, since the size of each pixel is generally small, the magnitude of intensity of the bright section 40 in one pixel can be considered to be approximately the same as that in its adjacent pixel, and the magnitude of movement of the bright section 40 in one pixel can be considered to be approximately the same as that in its adjacent pixel. Therefore, a normalized intensity change which is obtained by dividing a difference between the outputs of the adjacent pixels by a sum of these outputs does not depend on the intensity of the bright sections 40, but is proportional only to an average value of the amounts of movements of the bright sections 40 of these two pixels.
That is, the normalized intensity change A1 at an intermediate position P1′ of the pixels P1 and P2 is defined by the following expression (2), and the normalized intensity change A2 at the intermediate position P2′ of the pixels P2 and P3 is defined by the following expression (3). Hereinafter, if the normalized intensity change Ai at the intermediate position Pi′ of the pixels Pi and Pi+1 is likewise defined by the following expression (4), in general, the normalized intensity change Ai (i=1, 2, 3, . . . ) is proportional only to the change of the bright section in the x direction without being influenced by the transmissivity distribution of the object 20 to be detected.
A
1=−(I2−I1)/(I2+I1) (2)
A
2=(I3−I2)/(I3+I2) (3)
Generally, the following expression is established.
A
i=(−1)i×(Ii+1−Ii)/(Ii+1+Ii) (4)
Here, it should be noted that I1 is the output of the pixel P1, I2 is the output of the pixel P2, and Ii is the output of the pixel Pi. Further, if it is assumed that the period of the bright-dark pattern is D, the width of the bright section 40 is D/2. Further, since the half of the bright section is blocked, the width of the X-ray received by the pixel is D/4.
If it is assumed that an average of the amount of the movements of the bright sections 40 in the pixels Pi and Pi+1 is ΔLi, as each of Ii and Ii+1 is proportional to an area of the not-blocked portion of the brightness section, Ii is proportional to D/4+ΔLi if i is an odd number, while Ii is proportional to D/4−ΔLi if i is an even number.
Consequently, the following expression (5) is obtained from the expression (4).
Here, the period D of the bright-dark pattern is expressed by the following expressions (6) and (7) using the grating period d of the first grating 30, the distance Z0 between the X-ray source and the first grating, and the distance Z1 between the first grating and the second grating 50. That is, in the case where the one-dimensional π phase grating is used, the period D is expressed by the following expression (6).
D=(1+Z1/Z0)×d/2 (6)
Further, in the case where the one-dimensional r/2 phase grating is used, the period D is expressed by the following expression (7).
D=(1+Z1/Z0)×d (7)
The inclination Wxi of the wavefront in the periodic direction (x direction in
Wx
i
=ΔL
i
/Z
1 (8)
Therefore, in the case where the one-dimensional π phase grating is used, the following expression (9) is obtained from the expressions (8), (5) and (6).
Wx
i
=A
i/8×(1/Z0+1/Z1)×d (9)
Likewise, in the case where the one-dimensional π/2 phase grating is used, the following expression (10) is obtained from the expressions (8), (5) and (7).
Wx
i
=A
i/4×(1/Z0+1/Z1)×d (10)
Furthermore, in the case where the one-dimensional π phase grating is used, the following expression (11) is obtained from the expressions (9) and (10) on the basis of the expression (1).
Wx
i=λ/8Nd×Ai (11)
Likewise, in the case where the one-dimensional r/2 phase grating is used, the following expression (12) is obtained.
Wx
i=λ/4Nd×Ai (12)
The above-described calculation is performed for the overall pixels in the light reception region of the X-ray intensity detector 70, whereby the calculator 80 calculates the differential wavefront Wx of the transmitted wavefront 13 with a high degree of accuracy even if an uneven distribution is included in the transmissivity of the object 20 to be detected.
At one edge (right edge in
Since the X-ray wavefront in the region A is not influenced by the object 20 to be detected, the transmitted wavefront 13 is obtained by integrating the differential wavefront Wx from the region A toward the other edge (x direction in
Since the shape of the transmitted wavefront indicates a path length in the path through which the X-ray of the object 20 to be detected is transmitted, a spatial change of a refractive index in the object 20 to be detected can be known by measuring the transmitted wavefront 13, and further an aspect of a spatial change of the composition can be known based on the spatial change of the refractive index.
Incidentally, the object 20 to be detected is disposed above the first grating 30 (+ side on the z axis in
An example 2 will be described with reference to
The first grating in this example is a phase modulation grating of a phase difference π which is two-dimensionally and periodically arranged (hereinafter, called a two-dimensional n phase grating).
Therefore, also the second grating has a two-dimensional periodic structure. Since the constituent elements other than the first grating and the second grating are the same as those in the example 1, these elements respectively have the same corresponding numerals as those illustrated in
The distance Z0 between the X-ray source and the first grating and the distance Z1 between the first grating and the second grating satisfy the above expression (1)
Here, d indicates the grating period of the first grating illustrated in
The distance Z0 between the X-ray source and the first grating, the distance Z1 between the first grating and the second grating, and the grating period d are selected so that the period D is equivalent to 1/integer number of a pixel period p of the X-ray intensity detector 70.
More specifically, respective portions Pij (i, j=1, 2, 3, . . . ) correspond to the pixels of the X-ray intensity detector 70, in which the borders of the pixels are indicated by the dotted lines, and a portion 52 corresponds to the second grating. A hatched portion 53 corresponds to the blocking section of the second grating 52, and a non-hatched portion corresponds to the light transmitting section. The blocking sections and the light transmitting sections are periodically arranged like 45°-inclined checkers. A period D2 of the second grating in the pixel arrangement direction is twice the pixel period p.
The direction in which the blocking section of the second grating 52 blocks a bright section 41 is set to be identical on the same pixel, the directions in which the blocking sections respectively block the bright sections on the adjacent pixels are set to be perpendicular to each other, and the directions in which the blocking sections respectively block the bright sections on the oblique pixels are set to be opposite to each other.
That is, the first blocking pattern of the first region on the pixel P2i−1,2j−1 blocks the lower right half. Further, the second blocking pattern of the second region on the pixel P2i,2j blocks the upper left half. Furthermore, the third blocking pattern of the third region on the pixel P2i−1,2j blocks the lower left half. Furthermore, the fourth blocking pattern of the fourth region on the pixel P2i,2j−1 blocks the upper right half.
Here, the direction in which the blocking section 53 blocks the bright section 41 is the lower right of the pixel P11, the lower left of the pixel P12, the upper right of the pixel P21, and the upper left of the pixel P22. Hereinafter, the bright sections are likewise blocked every two pixel periods.
In the above state, if the bright section moves in the x direction due to an inclination of the wavefront, the amount of received light increases in each of the pixels P12 and P22, and the amount of received light decreases in each of the pixels P11 and P21. On the other hand, if the bright section moves in the y direction, the amount of received light increases in each of the pixels P11 and P12, and the amount of received light decreases in each of the pixels P21 and P22. Here, the x and y directions coincide with arrangement directions of pixels 61 (this is the same also in the following). Further, since the size of each pixel is generally small, the magnitude of intensity of the bright section 41 in one pixel can be considered to be approximately the same as that in its adjacent pixel, and also the magnitude of movement of the bright section 41 in one pixel can be considered to be approximately the same as that in its adjacent pixel. Therefore, a normalized intensity change which is obtained by dividing a difference between the outputs of the adjacent pixels by a sum of these outputs does not depend on the intensity of the bright sections 41, but is proportional only to an average value of the amounts of movements of the bright sections 41 of these four pixels.
That is, the normalized intensity changes A11 and B11 respectively concerning the movement components in the x and y directions at an intermediate position P11′ of the pixels P11, P12, P21 and P22 are respectively defined by the following expressions (13) and (14).
Further, the normalized intensity changes A21 and B21 respectively concerning the movement components in the x and y directions at an intermediate position P21′ of the pixels P21, P22, P31 and P32 are respectively defined by the following expressions (15) and (16).
In general, if the normalized intensity changes Aij and Bij respectively concerning the movement components in the x and y directions at an intermediate position of the pixels Pi,j, Pi,j+1, Pi+1,j and Pi+1,j+1 are respectively defined by the following expressions (17) and (18), the normalized intensity changes Aij and Bij (i, j=1, 2, 3, . . . ) are proportional only to the changes of the bright section 41 in the x and y directions without being influenced by the transmissivity distribution of the object 20 to be detected.
A
11=(I12+I22−I11−I21)/(I12+I22+I11+I21) (13)
B
11=(I11+I12−I21−I22)/(I11+I12+I21+I22) (14)
A
21=(I22+I32−I21−I31)/(I22+I32+I21+I31) (15)
B
21=(I32+I31−I21−I22)/(I32+I31+I21+I22) (16)
A
ij=(−1)j×(Ii,j+Ii+1,j−Ii,j+1−Ii+1,j+1)/(Ii,j+Ii+1,j+Ii,j+1+Ii+1,j+1) (17)
B
ij=(−1)i×(Ii+1,j+Ii+1,j+1−Ii,j−Ii,j+1)/(Ii+1,j+Ii+1,j+1+Ii,j+Ii,j+1) (18)
The inclination Wxij of the wavefront in the x direction is obtained from the following expression (19), in the same manner as that for obtaining Wxi by the expression (11) in the example 1.
Wx
ij=λ/4Nd×Aij (19)
Likewise, the inclination Wyij of the wavefront in the y direction is obtained from the following expression (20).
Wy
ij=λ/4Nd×Bij (20)
Since the inclinations of the wavefront in the x and y directions are obtained as described above, the transmitted wavefront of the object 20 to be detected can be obtained by integrating the obtained inclinations.
In this example, since the inclination can be integrated in both the x and y directions, the transmitted wavefront of the object 20 to be detected can accurately be calculated even if there is no region which is provided on one edge of the X-ray intensity detector and to which the X-ray not passing through the object 20 to be detected is radiated.
Although the transmitted wavefront of the object 20 to be detected is calculated based on the movement components in the x and y directions of the bright section 41 in this example, the transmitted wavefront can also be calculated based on the movement components in 45° and −45° directions.
If the bright section 41 moves in the 450 direction in
That is, the normalized intensity changes A11 and B11 respectively concerning the movement components in the 45° and −450 directions at the intermediate position P11′ of the pixels P11, P12, P21 and P22 are respectively defined by the following expressions (21) and (22).
Further, the normalized intensity changes A21 and B21 respectively concerning the movement components in the 45° and −45° directions at the intermediate position P21′ of the pixels P21, P22, P31 and P32 are respectively defined by the following expressions (23) and (24).
In general, if the normalized intensity changes Aij and Bij respectively concerning the movement components in the 45° and −45° directions at the intermediate position Pi,j′ of the pixels Pi,j, Pi,j+1, Pi+1,j and Pi+1,j+1 are respectively defined by the following expressions (25) and (26), the normalized intensity changes Aij and Bij (i, j=1, 2, 3, . . . ) are proportional only to the changes of the bright section 41 in the 45° and −45° directions without being influenced by the transmissivity distribution of the object 20 to be detected.
A
11=(I12−I21)/(I12+I21) (21)
B
11=(I22−I11)/(I22+I11) (22)
A
21=(I32−I21)/(I32+I21) (23)
B
21=(I22−I31)/(I22+I31) (24)
A
ij=(−1)j×(Ii,j−Ii+1,j+1)/(Ii,j+Ii+1,j+1) (i+j: even number)
A
ij=(−1)j×(Ii+1,j−Ii,j+1)/(Ii+1,j+Ii,j+1) (i+j: odd number) (25)
B
ij=(−1)j×(Ii+1,j−Ii,j+1)/(Ii+1,j+Ii,j+1) (i+j: even number)
B
ij=(−1)j×(Ii+1,j−Ii,j+1)/(Ii+1,j−Ii,j+1) (i+j: odd number) (26)
The inclination Waij of the wavefront in the 450 direction is obtained from the following expression (27), in the same manner as that for obtaining Wxi by the expression (11) in the example 1.
Wa
ij=λ/4Nd×Aij (27)
Likewise, the inclination Wbij of the wavefront in the −45° direction is obtained from the following expression (28).
Wb
ij=λ/4Nd×Bij (28)
Since the inclinations of the wavefront in the 45° and −45° directions are obtained as described above, the transmitted wavefront of the object 20 to be detected can be obtained by integrating the obtained inclinations.
As just described, the case where the period D of the bright-dark pattern is equivalent to the pixel period p of the X-ray intensity detector 70 is explained. The wavefront inclinations Wxij and Wyij can be likewise obtained even in a case where the period D is equivalent to 1/integer number of the pixel period p. To decrease the period D implies to decrease the grading period d of the first grating. Hence, since the distance Z0 between the X-ray source and the first grating and the distance Z1 between the first grating and the second grating can be made small according to the above expression (1), the overall apparatus can be made compact in size.
More specifically, a portion 54 corresponds to the second grating, a hatched portion 55 corresponds to the blocking section of the second grating, and the dotted lines indicate the borders of the pixels.
An example 3 will be described with reference to
Here, the blocking direction of the second grating in this example is different from that in the example 2, but other constituent elements are the same as those in the example 2.
Here, respective portions Pij (i, j=1, 2, 3, . . . ) correspond to the pixels of the X-ray intensity detector 70, in which the borders of the pixels are indicated by the dotted lines, and a portion 56 corresponds to the second grating.
A hatched portion 57 corresponds to the blocking section of the second grating 56, and a non-hatched portion corresponds to the light transmitting section. The blocking sections and the light transmitting sections are periodically arranged so as to alternately block the upper, lower, left and right halves of the pixels. The period D2 of the second grating 56 in the pixel arrangement direction is twice the pixel period p. As well as the example 2, the direction in which the blocking section of the second grating 56 blocks the bright section is set to be identical on the same pixel, the directions in which the blocking sections respectively block the bright sections on the adjacent pixels are set to be perpendicular to each other, and the directions in which the blocking sections respectively block the bright sections on the oblique pixels are set to be opposite to each other.
That is, the first blocking pattern of the first region on the pixel P2i−1,2j−1 blocks the right half. Further, the second blocking pattern of the second region on the pixel P2i,2j blocks the left half. Furthermore, the third blocking pattern of the third region on the pixel P2i−1,2j blocks the lower half. Furthermore, the fourth blocking pattern of the fourth region on the pixel P2i,2j−1 blocks the upper half.
Here, the blocking direction in this example is different from those of the second gratings 52 and 54 in the example 2. Namely, the direction in which the blocking section blocks the bright section is the right part of the pixel P11, the lower part of the pixel P12, the upper part of the pixel P21, and the left part of the pixel P22. Hereinafter, the bright sections are likewise blocked every two pixel periods.
In the above state, if the bright section 43 moves in the x direction due to an inclination of the wavefront, the amount of received light increases in the pixel P22, the amount of received light decreases in the pixel P11, and the amount of received light does not change in each of the pixels P22 and P21. On the other hand, if the bright section moves in the y direction, the amount of received light increases in the pixel P12, the amount of received light decreases in the pixel P21, and the amount of received light does not change in each of the pixels P11 and P22.
As well as the example 2, the normalized intensity change is introduced.
That is, the normalized intensity changes A11 and B11 respectively concerning the movement components in the x and y directions at the intermediate position P11′ of the pixels P11, P12, P21 and P22 are respectively defined by the following expressions (29) and (30).
Further, the normalized intensity changes A21 and B21 respectively concerning the movement components in the x and y directions at the intermediate position P21′ of the pixels P21, P22, P31 and P32 are respectively defined by the following expressions (31) and (32).
In general, if the normalized intensity changes Aij and Bij respectively concerning the movement components in the x and y directions at the intermediate position Pi,j′ of the pixels Pi,j, Pi,j+1, Pi+1,j and Pi+1,j+1 are respectively defined by the following expressions (33) and (34), the normalized intensity changes Aij and Bij (i, j=1, 2, 3, . . . ) are proportional only to the changes of the bright section 43 in the x and y directions without being influenced by the transmissivity distribution of the object 20 to be detected.
A
11=(I22−I11)/(I22+I11) (29)
B
11=(I12−I21)/(I12+I21) (30)
A
21=(I22−I31)/(I22+I31) (31)
B
21=(I32−I21)/(I32+I21) (32)
A
ij=(−1)j×(Ii,j−Ii+1,j+1)/(Ii,j+Ii+1,j+1) (i+j: even number)
A
ij=(−1)j×(Ii,j+1−Ii+1,j)/(Ii+1,j+Ii,j+1) (i+j: odd number) (33)
B
ij=(−1)j×(Ii+1,j−Ii,j+1)/(Ii+1,j+Ii,j+1) (i+j: even number)
B
ij=(−1)j×(Ii,j−Ii+1,j+1)/(Ii,j+Ii+1,j+1) (i+j: odd number) (34)
After then, the inclinations of the wavefront in the x and y directions are obtained from the expressions (19) and (20) as well as the example 2, whereby the transmitted wavefront of the object 20 to be detected can be obtained by integrating the obtained inclinations.
As just described, the case where the period D of the bright-dark pattern is equivalent to the pixel period p of the X-ray intensity detector 70 is explained. The wavefront inclinations Wxij and Wyij can be likewise obtained even in the case where the period D is equivalent to 1/integer number of the pixel period p.
Here, the direction in which the blocking section 59 blocks the bright section 44 is set to be identical on the same pixel, the directions in which the blocking sections respectively block the bright sections on the adjacent pixels are set to be perpendicular to each other, and the directions in which the blocking sections respectively block the bright sections on the oblique pixels are set to be opposite to each other. Namely, the overlapping state is the same as that illustrated in
Therefore, the wavefront inclinations Wxij and Wyij can be calculated by using the expressions (33), (34), (19) and (20).
An example 4 will be described with reference to FIGS. 7A, 7B, 8A, 8B, 9A and 9B.
Here, the first grating in this example is a phase modulation grating of a phase difference n/2 which is two-dimensionally and periodically arranged (hereinafter, called a two-dimensional n/2 phase grating).
Therefore, also the second grating has a two-dimensional periodic structure of which the pattern is different from that in each of the examples 2 and 3.
Since the constituent elements other than the first grating and the second grating are the same as those in the example 1, these elements respectively have the same corresponding numerals as those illustrated in
The distance Z0 between the X-ray source and the first grating and the distance Z1 between the first grating and the second grating satisfy the above expression (1). In this case, d indicates the grating period of the first grating illustrated in
The period D of the bright-dark pattern is expressed by the expression (7) as well as the example 1. The distance Z0 between the X-ray source and the first grating, the distance Z1 between the first grating and the second grating, and the grating period d are selected so that the period D is equivalent to 1/integer number of the pixel period p of the X-ray intensity detector 70.
Here, respective portions Pij (i, j=1, 2, 3, . . . ) correspond to the pixels of the X-ray intensity detector 70, and the borders of the pixels are indicated by the dotted lines, and a portion 60 corresponds to the second grating. A hatched portion 61 corresponds to the blocking section of the second grating 60, and a non-hatched portion corresponds to the light transmitting section. The blocking sections and the light transmitting sections are two-dimensionally and periodically arranged. The period D2 of the second grating 60 in the pixel arrangement direction is twice the pixel period p. As well as the example 2, the direction in which the blocking section of the second grating 60 blocks the bright section is set to be identical on the same pixel, the directions in which the blocking sections respectively block the bright sections on the adjacent pixels are set to be perpendicular to each other, and the directions in which the blocking sections respectively block the bright sections on the oblique pixels are set to be opposite to each other.
That is, the first blocking pattern of the first region on the pixel P2i−1,2j−1 blocks the upper left half. Further, the second blocking pattern of the second region on the pixel P2i,2j blocks the lower right half. Furthermore, the third blocking pattern of the third region on the pixel P2i−1,2j blocks the lower left half. Furthermore, the fourth blocking pattern of the fourth region on the pixel P2i−1,2j−1 blocks the upper right half.
Here, the direction in which the blocking section 61 blocks the bright section 43 is the upper left part of the pixel P11, the lower left part of the pixel P12, the upper right part of the pixel P21, and the lower right part of the pixel P22. Hereinafter, the bright sections are likewise blocked every two pixel periods.
In the above state, if the bright section 45 moves in the x direction due to an inclination of the wavefront, the amount of received light increases in each of the pixels P11 and P12, and the amount of received light decreases in each of the pixels P21 and P22. On the other hand, if the bright section moves in the y direction, the amount of received light increases in each of the pixels P12 and P22, and the amount of received light decreases in each of the pixels P11 and P21. As well as the example 2, the normalized intensity change is introduced.
That is, the normalized intensity changes A11 and B11 respectively concerning the movement components in the x and y directions at the intermediate position P11′ of the pixels P11, P12, P21 and P22 are respectively defined by the following expressions (35) and (36).
Further, the normalized intensity changes A21 and B21 respectively concerning the movement components in the x and y directions at the intermediate position P21′ of the pixels P21, P22, P31 and P32 are respectively defined by the following expressions (37) and (38).
In general, if the normalized intensity changes Aij and Bij respectively concerning the movement components in the x and y directions at the intermediate position Pi,j′ of the pixels Pi,j, Pi,j+1, Pi+1,j and Pi+1,j+1 are respectively defined by the following expressions (39) and (40), the normalized intensity changes Aij and Bij (i, j=1, 2, 3, . . . ) are proportional only to the changes of the bright section 41 in the x and y directions without being influenced by the transmissivity distribution of the object 20 to be detected.
A
11=(I11+I12−I21−I22)/(I11+I12+I21+I22) (35)
B
11=(I12+I22−I11−I21)/(I12+I22+I11+I21) (36)
A
21=(I31+I32−I21−I22)/(I31+I32+I21+I22) (37)
B
21=(I22+I32−I21−I31)/(I22+I32+I21+I31) (38)
A
ij=(−1)i×(Ii+1,j+Ii+1,j+1−Ii,j−Ii,j+1)/(Ii+1,j+Ii+1,j+1−Ii,j−Ii,j+1) (39)
B
ij=(−1)j×(Ii,j+Ii+1,j−Ii,j+1−Ii+1,j+1)/(Ii,j+Ii+1,j−Ii,j+1−Ii+1,j+1) (40)
The inclination Wxij of the wavefront in the x direction is obtained from the following expression (41), in the same manner as that for obtaining Wxi by the expression (12) in the example 1.
Wx
ij=λ/2Nd×Aij (41)
Likewise, the inclination Wyij of the wavefront in the y direction is obtained from the following expression (42).
Wy
ij=λ/2Nd×Bij (42)
Since the inclinations of the wavefront in the x and y directions are obtained as described above, the transmitted wavefront of the object to be detected can be obtained by integrating the obtained inclinations.
Although the transmitted wavefront of the object 20 to be detected is calculated based on the movement components in the x and y directions of the bright section 45, it is apparent that the transmitted wavefront can also be calculated based on the movement components in 45° and −45° directions as well as the example 2. Therefore, the detail of such calculation will be omitted.
As just described, the case where the period D of the bright-dark pattern is equivalent to the pixel period p of the X-ray intensity detector 70 is explained. The wavefront inclinations Wxij and Wyij can be likewise obtained even in a case where the period D is equivalent to 1/integer number of the pixel period p.
Here, a portion 62 corresponds to the second grating, a hatched portion 63 corresponds to the blocking section, and the dotted lines indicate the borders of the pixels.
Here, the directions in which the blocking sections 63 respectively block the bright sections 46 are set to be the same on the same pixel, the directions in which the blocking sections 63 respectively block the bright sections 46 on the adjacent pixels are set to be perpendicular to each other, and the directions in which the blocking sections 63 respectively block the bright sections 46 on the oblique pixels are set to be opposite to each other. Namely, the above overlapping state is the same as that illustrated in
Therefore, the wavefront inclinations Wxij and Wyij can be calculated by using the expressions (39) to (42).
An example 5 will be described with reference to
Here, the first grating in this example is a two-dimensional amplitude grating of which the transmitting section is two-dimensionally and periodically arranged.
In the drawing, a first grating 80 includes an absorption section 81 through which an X-ray does not transmit, and transmitting sections 82 through which X-rays transmit and are two-dimensionally arranged with periods d.
If the distance Z1 between the first grating and the second grating is made small so that N in the expression (1) has a value sufficiently smaller than 1, the X-ray which transmitted through the transmitting section 82 approximately goes straight, thereby generating the bright-dark pattern on the second grating.
The bright-dark pattern is similar to the arrangement pattern of the transmitting sections 82 of the first grating 80, and the period D thereof is expressed by the following expression (43).
In the expression, each of the distance Z0 between the X-ray source and the first grating, the distance Z1 between the first grating and the second grating, and the grating period d is selected to become 1/integer number of the pixel period p of the X-ray intensity detector 70.
D=(1+Z1/Z0)×d (43)
The second grating has the periodic-pattern blocking section. Thus, as well as the above examples, the blocking section blocks the bright sections 47 in the same direction in the pixel of the arranged X-ray intensity detector 70, and in the different directions respectively in the adjacent pixels.
In this example, since the bright-dark pattern is the same as those in the examples 2 and 3, the second grating has the pattern of the second grating illustrated in
Although a calculation expression of the transmitted wavefront of the object 20 to be detected depends on the pattern of the second grating, it will be easily led with reference to the calculation expressions as described in the above examples.
As described above, according to the present invention, even if there is the non-uniformly distributed absorption in the object to be detected, it is possible to reduce the influence of the transmissivity distribution of the X-ray in the object to be detected, and it is thus possible to achieve the X-ray imaging apparatus capable of taking the phase image having exceptional quality.
While the present invention has been described with reference to the exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
This application claims the benefit of Japanese Patent Application No. 2010-036909, filed Feb. 23, 2010, which is hereby incorporated by reference herein in its entirety.
Number | Date | Country | Kind |
---|---|---|---|
2010-036909 | Feb 2010 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP2011/053585 | 1/14/2011 | WO | 00 | 6/21/2012 |