WO2012029005A1

WO2012029005A1 - Differential phase-contrast imaging with improved sampling

Info

Publication number: WO2012029005A1
Application number: PCT/IB2011/053750
Authority: WO
Inventors: Thomas Koehler
Original assignee: Koninklijke Philips Electronics N.V.; Philips Intellectual Property & Standards Gmbh
Priority date: 2010-09-03
Filing date: 2011-08-26
Publication date: 2012-03-08
Also published as: US9105369B2; JP2013541699A; CN103079469A; JP5961614B2; US20130170618A1; EP2611364B1; EP2611364A1; CN103079469B

Abstract

The present invention relates to differential phase-contrast imaging of an object. For increasing spatial resolution of an X-ray imaging system (2) the size of a detector pixel element (8) may be considered a limiting factor. Accordingly, it may be beneficial to increase the resolution of an apparatus (38) for phase-contrast imaging without further reducing the area of an individual pixel element (8). Accordingly, an apparatus (38) for phase-contrast imaging with improved sampling is provided, comprising an X-ray source (4), a first grating element G₁ (24), a second grating element G₂ (26) and an X-ray detector element (6) comprising a plurality of detector pixel elements (8), each detector pixel element (8) having a pixel area A. An object to be imagined (14) is arrangeable between the X-ray source (4) and the X-ray detector element (6). The first grating element G₁ (24) as well as the second grating element G₂ (26) are arrangeable between the X-ray source (4) and the X-ray detector element (6). The X-ray source (4), the first grating element G₁ (24), the second grating element G₂ (26) and the X-ray detector (6) are operatively coupled for acquisition of a phase-contrast image of the object (14). At least one of the first grating element G₁ (24) and the second grating element G₂ (26) comprise a first area A ₁ having a first grating pitch p ₁ and a second area A ₂ having a second grating pitch p ₂ different from the first grating pitch.

Description

DIFFERENTIAL PHASE-CONTRAST IMAGING WITH IMPROVED SAMPLING

FIELD OF THE INVENTION

The present invention relates to X-ray imaging technology in general.

More particularly, the present invention relates to differential phase-contrast imaging.

In particular, the present invention relates to an apparatus for phase-contrast imaging with improved sampling, an X-ray system comprising an apparatus according to the present invention as well as the use of an apparatus according to the present invention in one of an X-ray system and a CT system.

BACKGROUND OF THE INVENTION

When acquiring an X-ray image, an object to be examined, e.g. a patient, is arranged between an X-ray generating device and an X-ray detector. X-ray radiation emanating from the X-ray generating device is penetrating the object to be examined, subsequently arriving at the X-ray detector. The object to be examined, situated in the path of the X-ray radiation is spatially attenuating the X-ray beam, depending on the specific tissue structure within the object. The X-ray detector is subsequently detecting the spatially attenuated X-ray radiation by determining an intensity distribution of the X-ray radiation, which image information is employed for generating, further processing, and subsequently displaying an X-ray image of the object to be examined.

However, an object to be examined may provide only minor differences when attenuating the X-ray radiation, resulting in a relatively uniformly attenuated X-ray image having low contrast, thus lacking detail of the imaged inner structure of the object.

While certain objects or regions within an object may comprise similar attenuation properties, a phase of X-ray radiation penetrating the object may be influenced to a larger extent by the structure of the object.

In phase-contrast imaging, at least partly spatially coherent X-ray radiation is employed, e.g., generated by a source grating arranged adjacent to, in the vicinity of an X-ray source, e.g. an X-ray tube. Coherent X-rays penetrating the object may allow for subsequent retrieval of phase information.

However, a phase of a wave may not be measured directly, rather a phase-shift may be required to be converted to an intensity modulation, e.g., by interfering two or more waves. For generating an according interference pattern, a so-called phase grating is employed, arranged between the object to be examined and an X-ray detector. However, an interference pattern generated by only employing a phase grating may be too small to be detectable with a current X-ray detector, due to a lack of spatial resolution of the X-ray detector.

Thus, a further analyzer grating may be employed arranged between the phase grating and the X-ray detector, subsequently providing an interference pattern, which is large enough to be detectable by current X-ray detectors.

To obtain appropriate image information, a so-called phase stepping is performed. In phase stepping, one of the source grating, the phase grating, and the analyzer grating is displaced laterally with respect to the other grating and the X-ray detector element by a fraction of its grating pitch, e.g., a fourth, sixth, eighth of the grating pitch, e.g. of the phase grating. If the phase stepping is performed using a particular grating, then the phase stepping shall cover a full period of this particular grating.

The spatial resolution of a grating based differential phase-contrast imaging system is typically limited by the size of the focal spot of the X-ray source and the size of a single detector pixel element.

While differential phase-contrast imaging may provide enhanced image information over transmission X-ray imaging taking into account attenuation only, the spatial resolution of the so obtained image information may be considered to still be limited by the individual detector pixel elements and their respective size.

Thus, for obtaining image information with increased detail, an improved or enhanced spatial resolution of acquired image information of a differential phase-contrast projection may be beneficial.

X-ray phase-contrast imaging is described in both Weitkamp T., Diaz A.,

David C. et al.: "X-ray phase imaging with a grating interferometer"; Optics Express 6296, 8.

August 2005, Vol. 13, No. 16 as well as Bartl P., Durst J., Haas W. et al. "Simulation of X- ray phase-contrast computed tomography of a medical phantom comprising particle and wave contributions", Proc of SPIE, Volume 7622, 76220Q-1. SUMMARY OF THE INVENTION

An aspect of the present invention may be seen as providing means for overcoming the limit of spatial resolution directly depending on the physical size of a detector pixel element, thus increasing spatial resolution beyond the actual size of the detector pixel element. This may allow improving the spatial resolution of the differential phase-contrast projection while employing a current X-ray detector or while maintaining the spatial resolution of the differential phase-contrast projection relaxing the requirements of the spatial resolution of the detector, e.g., by increasing the individual pixel area of an individual detector pixel element.

Accordingly, an apparatus for phase-contrast imaging with improved sampling, an X-ray system comprising an apparatus according to the present invention as well as the use of an apparatus according to the present invention in at least one of an X-ray system and a CT system according to the independent claims is provided.

Preferred embodiments may be taken from the dependent claims.

In particular, it may be beneficial to increase the resolution of an apparatus for phase-contrast imaging without further reducing the area of an individual pixel element.

These and other aspects of the present invention will become apparent from and elucidated with reference to the embodiments described hereinafter.

Exemplary embodiments of the present invention will be described below with reference to the following drawings.

The illustration in the drawings is schematic. In different drawings, similar or identical elements are provided with similar or identical reference numerals.

The figures are not drawn to scale, however may depict qualitative

proportions.

BRIEF DESCRIPTION OF THE DRAWINGS

Fig 1 shows an exemplary embodiment of an X-ray system according to the present invention;

Fig 2 shows an exemplary detail of an arrangement for phase contrast

imaging according to the present invention; and

Fig 3 shows an exemplary detail of a grating arrangement according to the present invention. DETAILED DESCRIPTION OF EMBODIMENTS

One aspect of the present invention may be seen as employing a phase grating and an analyzer grating having a non-uniform or varying pitch structure with regard to each detector pixel. In particular, each of the phase grating and the analyzer grating may be seen as being divided in two individual areas with each individual area corresponding to the area and/or size of a single detector pixel element of the X-ray detector.

Within each area corresponding to a single detector pixel element, the grating structure of the phase grating and the analyzer grating may be a non-uniform grating structure. The non-uniform grating structure may be seen as employing at least two individual grating pitches for each detector pixel element.

The grating structure of each grating element may be seen as comprising individual barrier elements, each forming a barrier region, spaced apart from one another, thus forming a trench region between the barrier elements. Preferably, the trench region and the barrier region both comprise the same width, thus the trench region and the barrier region or barrier element are substantially of the same dimension.

The distance between two barrier elements arranged adjacent to each other may be referred to as the pitch of the grating. Thus, the pitch of a grating structure is either the width of a trench region plus the width of a barrier region or, since barrier region and trench region preferably comprise the same width, the pitch of the grating structure equals also two times either the width of a trench region or a barrier region.

A pitch of the grating may also be referred to as the periodicity of the grating.

Preferably, the area or size of the gratings corresponding to the area or size of a single detector pixel element may be halved with each half of the area corresponding to a single detector pixel element comprising an individual grating structure having an individual grating pitch. In other words, a first half of the area of one detector pixel element is covered by a grating structure having a first grating pitch and a second half of the area of one detector pixel element is covered by a grating structure having a second grating pitch, different from the first grating pitch.

In a particular embodiment, two different grating pitches p_\ and /¾ for phase grating Gi may be employed each covering the respective areas

i.e., half of the pixel area A of a single detector pixel element. A distance d between a first grating element, a phase grating Gi, and a second grating element, an analyzer grating G₂, may be exemplarily selected as the first fractional Talbot distance related to a pitch, e.g., pitch p_\ . The n-t fractional Talbot distance may be determined by employing the equation

Equation 1

where λ is the wavelength at a defined design energy of the X-ray source and n is a positive integer. If a spherical wave is employed, a further correction may be required, which is known in the art. Thus, the first fractional Talbot distance equals

8 - A

Equation 2

In case the precise wavelength λ is employed, a sharply defined interference pattern without blurring may be obtained at a fractional Talbot distance.

A second grating pitch p₂ may be selected to be

_i

p₂ = (2n - 1) 2 . _Pl

Equation 3

Employing such a second grating pitch results in a further fractional Talbot distance relating to pitch p₂ to equal distance d as well. In other words, distance <i is a fractional Talbot distance for pitch p_\ as well as a further, different fractional Talbot distance for/¾.

In case of equation 2 and equation 3 with n=2, distance d constitutes the first fractional Talbot distance for pitch p_\ as well as the third fractional Talbot distance for pitch Ρ2 at the same design energy.

Analyzer grating G₂ also comprises areas A i and A₂ having individual pitches q\ and q₂. The respective pitches of p\ and q\ as well as p₂ and q₂ may be preferably related by equation 4.

I + d

<li = Pi - ^~2j- ⁾ i = l,2

Equation 4

Each individual area A_\, A₂ may be considered to create an individual virtual pixel on the physical detector pixel element.

If now phase stepping is conducted, the intensity variation on the first virtual pixel may be considered to comprise the period of pitch p\, whereas the intensity variation of the second virtual pixel comprise the period of pitch p₂. Thus, the phase gradients of the respective two virtual pixels may be seen as being coded at different frequencies and thus they may both be retrieved from phase stepping. However, the minimum number of projections acquired during the phase stepping phase may have to be increased for allowing obtaining the independent parameters required for both virtual pixels.

E.g., with two different grating pitches p\ and /¾ the minimum number of projections may to be required to be increased from three different phase stepping states to five phase stepping states. In other words, when employing a single pitch only, three individual alignment positions of the phase grating Gi and the analyzer grating G₂ relative to one another may be considered sufficient for obtaining image information, while, when employing two different grating pitches five independent projections may have to be acquired.

Exemplarily, the intensity modulation in area A _\ may be modeled by the following equation.

Equation 5

where x is the displacement of phase grating Gi and the phase ( i is related to the mean gradient

of the x-ray wave front in the area A i according to the following equation. Furthermore, C_\ is related to the mean intensity during the phase stepping, which quantity relates to a conventional X-ray transmission image. V\ is related to the so-called visibility function, which is affected by the bandwidth of the X-ray source, the spatial coherence of the X-ray beam, the scattering in the object and other effects.

Equation 6

Correspondingly, the intensity modulation in the area A₂ can be modeled by the following equation

2 · Tc ' x

l₂ = C₂ + V₂ - cos( + ψ₂)

P

with

Equation 7

Subsequently, the total intensity measured by the individual detector pixel element may be taken from the following equation. / = I + I₂ = (Ci + C₂) + V₁ · cos( ^{2 X} + + V · cos ^{2 X} + φ₂)

Equation 8

The two phase gradients of interest d<i>/dx(A_\), d<i>/dx(A₂) may thus be obtained during phase stepping because they are coded in different frequencies l/pi and l//¾. In order to retrieve the two phase gradients from the phase stepping, the grating Gi should be moved now at least over the larger length of p_\ and /¾.

The concept of the present invention may not be restricted to the partitioning of the physical pixel element into two virtual pixels only. Rather by employing further different grating pitches, e.g., three different grating pitches within the area of a single detector pixel element, e.g., three virtual pixels may be generated. In such a case, seven individual projections may be required to be acquired. Accordingly, it may be beneficial that distance d equals a distance of a certain fractional Talbot distance or all grating pitches p_\, pi, P3, . .. p_n employed, while the number of individual images required to be acquired for reconstructing the plurality of virtual pixels it may be 2n+l .

Furthermore, it is also conceivable to select different fractional Talbot distances for the virtual pixels. In particular, distance d may also be for example the third fractional Talbot distance for pitch p_\ and e.g., the fifth fractional Talbot distance for pitch /¾. However, the sensitivity of the system to phase gradients may be seen as increasing linearly with the fractional Talbot order. In order to achieve a similar sensitivity for the different virtual pixels, it may be beneficial to employ similar fractional Talbot orders for different virtual pixels.

Now referring to Fig. 1 , an exemplary embodiment of an X-ray system according to the present invention is depicted.

In Fig. 1 , X-ray system 2 exemplarily is depicted as a CT system 2 comprising an X-ray source 4 and an X-ray detector 6 mounted on a gantry 10 to allow rotation of both the X-ray detector 6 and the X-ray source 4 about an axis.

X-ray source 4 generates X-ray radiation 20 emanating from X-ray source 4 and subsequently arriving at X-ray detector 6.

A support 12 is arranged within the gantry 10 of the CT system 2 for supporting an object 14 to be examined. Object 14 is arrangeable in the path of X-ray radiation 20 to allow acquisition of X-ray image information by X-ray detector 6 by its individual detector pixel elements 8. X-ray detector 6 is displayed exemplarily as a two- dimensional array of detector pixel elements 8. A control system 16 comprising a display 18 is communicatively coupled to X-ray system 2 for controlling the image acquisition process and also for providing means for a subsequent display of acquired and reconstructed image information on display 18.

Individual grating elements for differential phase-contrast imaging acquisition are not depicted in Fig. 1. However, a detailed illustration of the path of X-ray radiation 20 from X-ray tube 4 to detector element 6 may be taken from Fig. 2.

Now referring to Fig. 2, an exemplary embodiment of an arrangement for phase contrast imaging according to the present invention is depicted.

In arrangement 38, X-ray source 4 is arranged in the vicinity of a source grating 22 for generating at least partially spatial coherent X-ray radiation 20. Source grating 22 is spaced apart from phase grating Gi 24 by distance /. Object 14 is arranged in the path of X-ray radiation 20 between X-ray source 4 and phase grating Gi 24. Wave front 28a having a uniform phase is depicted as arriving at object 14 while a further phase front 28b having a changed phase relationship within the wave front due to a phase shift imposed on the wave front while penetrating the object 14 is depicted.

Subsequently, wave fronts arrive at phase grating Gi 24, which is arranged spaced apart by distance d from analyzer grating 26. Distance d may be a fractional Talbot distance for grating pitch p.

Analyzer grating G₂26 is displaceable 32 relative to phase grating Gi 24 for acquisition of phase-contrast images. However it is also conceivable to displace the source grating Go 22 or phase grating Gi 24 instead of analyzer grating G₂ 26.

X-ray radiation 20 passing through phase grating Gi 24 and analyzer grating G₂ 26 subsequently generates an interference pattern on X-ray detector element 6 which is detectably by detector pixel elements 8.

For sake of clarity in Fig. 2, phase grating Gi 24 is depicted having a uniform pitch p and analyzer grating G₂ 26 is depicted having a uniform pitch q. However, the detailed illustration regarding the individual pitch arrangements of both the phase grating Gi 24 as well as the analyzer grating G₂ 26 may be taken from Fig. 3.

Now referring to Fig. 3, an exemplary detail of a grating arrangement according to the present invention is depicted.

In Fig. 3 in grating arrangement 40, area A of a single detector pixel element 8 as well as correspondingly sized phase grating Gi 24 as well as analyzer grating G₂ 26 are depicted. It is to be understood that the depicted arrangement of Fig. 3 may in particular be repeated for every single detector pixel element 8 of X-ray detector 6. Area A of detector pixel element 8 is divided into two similarly sized areas A i and A₂, corresponding to areas of different pitch with regard to phase grating Gi 24 as well as analyzer grating G₂ 26. In particular, A = A_\ + A₂.

In area A_\, phase grating Gi 24 comprises pitch p_\ and analyzer grating G₂ 26 comprises pitch q_\ while in area A₂ phase grating Gi 24 comprises pitch p₂ and analyzer grating G₂ 26 comprises pitch q₂.

The size and dimension of the respective pitches p_\, p₂, q_\, q₂ may be determined and set according to the foregoing description. E.g., p_\ and p₂ may be related by

_ 1 _{j + £j}

p₂ = 3~2 · p while p_£and q_i may be related by q_i = p_t ·—.

Each grating comprises barrier elements or barrier regions 34 with intermediate trench regions 36, which may be considered to constitute gaps between the barrier elements or barrier regions 34. Each barrier region 24 and trench region 36 comprises half of the respective pitch p_\, p₂, q_\, q₂, depending on which area the respective barrier region 24 and trench region 36 is arranged in.

While with regard to phase grating Gi 24, the trench regions 36 may be substantially a space or gap between adjacent barrier regions 34, the trench regions 36 of the analyzer grating G₂ 26 may be filled with a material, e.g., gold by electroplating.

Barrier region 34 as well as trench region 36 may be of equal width, both constituting and adding up to the respective pitch of the respective area A_\, A₂ of the phase grating Gi 24 and the analyzer grating G₂ 26.

Exemplarily, the analyzer grating 26 may be displaceable 32 parallel to the phase grating Gi 24 for conducting phase stepping. Depending on the current phase stepping state, an interference pattern may be projected onto detector pixel element 8 either in the area of virtual pixel 1 30a or virtual pixel 2 30b. Since both virtual pixels 30a,b are arranged on the same detector pixel element 8, both virtual pixels are also read out via the same physical detector pixel element 8. The detector element 6 may thus not be required to be altered.

Subsequently, the acquired image information may be sorted with regard to an individual phase stepping state, i.e., with regard to whether an acquired image information is considered to belong to virtual pixel 1 30a or virtual pixel 2 30b. Thus, image information that is to be reconstructed by the plurality of readout values of the individual detector pixel elements 8 taken into account virtual pixels 30a,b may have a spatial resolution, at least in the direction perpendicular to the grating structure, i.e., perpendicular to the barrier regions 34 and trench regions 36, that is double the spatial resolution of the detector element 6 itself. It should be noted that the term "comprising" does not exclude other elements or steps and that "a" or "an" does not exclude a plurality. Also, elements described in association with different embodiments may be combined.

It should also be noted, that reference numerals in the claims shall not be construed as limiting the scope of the claims.

LIST OF REFERENCE SIGNS:

2 X-ray system

4 X-ray source

6 X-ray detector element

8 Detector pixel element

10 Gantry

12 Support

14 Object

16 Control system

18 Display

20 X-ray radiation

22 Source grating

24 Phase grating Gi

26 Analyzer grating G₂

28a,b Wave front

30a,b Virtual pixel

32 Displacement/phase stepping

34 Barrier element/barrier region

36 Trench region

38 Apparatus for phase contrast imaging

40 Grating arrangement

Claims

CLAIMS:

1. Grating arrangement (40) for phase contrast imaging, comprising

a first grating element Gi (24);

a second grating element G₂ (26); and

wherein at least one of the first grating element (24) and the second grating element (26) comprises a first area A i having a first grating pitch p_\,q_\ and a second area A₂ having a second grating pitch p qi different from the first grating pitch p_\,q_\ .

2. Grating arrangement according to the preceding claim, further comprising an X-ray detector element (6), comprising a plurality of detector pixel elements (8), each detector pixel element (8) having a pixel area A;

wherein the first area A 1 and the second area A₂ are arranged adjacently; and wherein the size of the first area A 1 and the second area A₂ correspond to the pixel area A of a single detector pixel element (8).

3. Grating arrangement according to at least one of the preceding claims, wherein the first area and the second area are of similar size.

4. Grating arrangement according to at least one of the preceding claims, wherein at least one of the first grating element (24) and the second grating element (26) comprises the first area having the first grating pitch p_\,q_\ and the second area having the second grating pitch /¾, ₂ different from the first grating pitch p_\,q_\ for each pixel element.

5. Apparatus (38) for phase contrast imaging, comprising

an X-ray source (4); and

a grating arrangement (40) according to one of claims 1 to 4;

wherein an object (14) to be imaged is arrangeable between the X-ray source (4) and the X-ray detector element (6);

wherein the first grating element (24) and the second grating element (26) are arrangeable between the X-ray source (4) and the X-ray detector element (6); wherein the X-ray source (4), the first grating element (24), the second grating element (26) and the X-ray detector element (6) are operatively coupled for acquisition of a phase contrast image of the object (14).

6. Apparatus according to the preceding claim,

wherein the X-ray source (4) is adapted to emit X-ray radiation (20) having a design energy level with a defined wave length;

wherein the first grating element (24) and the second grating element (26) are spaced apart by a distance d;

wherein distance d corresponds to a fractional Talbot distance of one of the first grating pitch p\ and the second grating pitch p₂ of the defined wave length.

7. Apparatus according to the preceding claim;

wherein the distance d corresponds to a fractional Talbot distance for both the first grating pitch p\ and the second grating pitch p₂.

8. Apparatus according to at least one of claims 5 to 7;

wherein the first grating pitch p_\ and the second grating pitch /¾ relate by

wherein distance d is the (2m-l)th fractional Talbot distance for the pitch p_\, and

wherein distance d is the (2/?-l)th fractional Talbot distance for the pitch /¾.

9. Apparatus according to at least one of claims 5 to 8;

wherein the distance d is one of the first and third fractional Talbot distance of the first grating pitch p_\ and the third and firth fractional Talbot distance of the second grating pitch

10. Apparatus according to at least one of claims 5 to 9,

wherein the first grating element (24) and the second grating element (26) are arranged substantially parallel; and

wherein the first grating element (24) and the second grating element (26) are adapted to be movable relative to one another for providing phase stepping.

11. X-ray system, comprising an apparatus (38) according to one of claims 5 to 10.

12. Use of an apparatus (38) according to one of claims 5 to 10 in at least one of an X-ray system and a CT system.