Differential phase contrast imaging with energy sensitive detection

1. A method for correcting differential phase image data ( 52 ), the method comprising the steps of

Receiving differential phase image data ( 52 ) acquired with radiation at different energy levels (E −1 , E 0 , E 1 ), wherein the differential phase image data ( 52 ) comprises pixels ( 60 ), each pixel ( 60 ) having a phase gradient value ( 62 a , 62 b , 62 c ) for each energy level (E −1 , E 0 , E I );

Determining an energy dependent behavior ( 44 a to 44 g ) of phase gradient values ( 62 a , 62 b , 62 c ) of a pixel ( 60 );

Determining a corrected phase gradient value ( 68 ) for the pixel ( 60 ) from the phase gradient values ( 62 a , 62 b , 62 c ) of the pixel ( 60 ) and a model for the energy dependence of the phase gradient values ( 62 a , 62 b , 62 c ).

2. The method of claim 1 ,

wherein the energy dependent behavior is determined by fitting the phase gradient values ( 62 a , 62 b , 62 c ) of the pixel ( 60 ) to a plurality of phase wrapping curves ( 44 a to 44 g ), each wrapping curve being associated with a wrapping number;

wherein the wrapping number of the pixel ( 60 ) is determined by determining a wrapping curve ( 44 a to 44 g ) with a best fit of the fitted phase gradient values ( 62 a , 62 b , 62 c ) of the pixel ( 60 ) and by selecting the wrapping number associated with the determined wrapping curve ( 44 a to 44 g ).

3. The method of claim 2 ,

wherein a wrapping curve ( 44 a to 44 g ) is represented by phase gradient values ( 42 , 46 ) associated with energy levels (E −1 , E 0 , E 1 );

wherein a phase gradient value ( 42 , 46 ) of an energy level (E −1 , E 0 , E 1 ) of a wrapping curve ( 44 a to 44 g ) is fitted with a phase gradient value ( 62 a , 62 b , 62 c ) of the pixel ( 60 ) at this energy level.

4. The method of claim 3 ,

wherein the phase gradient values ( 42 , 46 ) of the wrapping curves ( 44 a to 44 g ) are determined such that phase gradient values ( 42 ) at a reference energy level (E 0 ) are equal.

5. The method of claim 1 ,

wherein the pixel ( 60 ) has an attenuation value ( 64 );

wherein the method comprises the steps of:

Determining the gradient of the attenuation value ( 64 ) at the pixel ( 60 );

Choosing an estimated region of wrapping numbers;

Fitting the phase gradient values ( 62 a , 62 b , 62 c ) of a pixel ( 60 ) to a plurality of wrapping curves ( 44 a to 44 g ) associated with wrapping numbers in the estimated region.

6. The method claim 1 ,

wherein the corrected phase gradient value ( 68 ) is determined by shifting the phase gradient value ( 62 b ) of the pixel ( 60 ) associated with a reference energy level (E 0 ) with the determined wrapping number of the pixel ( 60 ).

7. The method of claim 1 ,

wherein a pixel has phase gradient values ( 62 a , 62 b , 62 c ) corresponding to at least three energy levels (E −1 , E 0 , E 1 ).

8. The method of claim 1 ,

wherein the radiation is electromagnetic radiation.

9. The method of claim 1 ,

wherein the energy levels E −1 , E 0 , E 1 of the radiation comprise a reference energy level (E 0 ) and two neighboring energy levels (E −1 , E 1 ) differing from 8% to 12% from the reference energy level (E 0 ).

10. A method for generating corrected differential phase image data ( 66 ), the method comprising the steps:

Generating radiation at different energy levels (E −1 , E 0 , E 1 );

Detecting the generated radiation penetrating an object of interest ( 28 );

Acquiring differential phase image data ( 52 ) from the detected radiation;

Executing the steps of claim 1 with the generated differential phase image data ( 52 ).

11. The method of claim 10 ,

wherein the differential phase image data is acquired in different directions with respect to the object of interest;

wherein the method comprises the step of:

Generating tomography image data from the corrected differential phase image data.

12. A non-transitory computer-readable storage medium having stored therein machine readable instructions configured to be executed by a processor for correcting differential phase image data ( 52 ), which, when executed by the processor ( 30 ), is adapted to carry out the steps of claim 1 .

13. A differential phase imaging system ( 10 ), comprising:

a radiation source ( 12 );

a detector ( 14 ); and

a controller ( 16 );

wherein the radiation source ( 12 ) is adapted to generate radiation of different energy levels;

wherein the detector ( 14 ) is adapted to detect differential phase image data ( 52 ) of an object of interest ( 28 ) penetrated by the radiation;

wherein the controller ( 16 ) is adapted to carry out the method according to claim 1 .

14. The differential phase imaging system ( 10 ) of claim 13 ,

wherein the differential phase imaging system is an X-ray CT system ( 10 ).