Robust reconstruction for dark-field and phase contrast CT
A system and related method for signal processing. Interferometric projection data reconstructed into one or more images for a spatial distribution of a physical property of an imaged object. The interferometric projection data is derived from signals acquired by an X-ray detector (D), said signals caused by X-ray radiation after interaction of said X-ray radiation with an interferometer and with the object (OB) to be imaged, said interferometer (IF) having a reference phase. A reconstructor (RECON) reconstructs for the image(s) by fitting said data to a signal model by adapting fitting variables, said fitting variables including i) one or more imaging variables for the one or more images and ii), in addition to said one or more imaging variables, a dedicated phase variable for a fluctuation of said reference phase.
1. A signal processing system, comprising:
at least one processor configured to:
receive interferometric projection data derived from signals acquired by an X-ray detector, said signals caused by X-ray radiation after interaction of said X-ray radiation with an interferometer and with an object to be imaged, said interferometer having a reference phase;
reconstruct one or more images of a spatial distribution of one or more physical properties of said object and fit said interferometric projection data to a signal model by adapting a plurality of fitting variables, said fitting variables including i) one or more imaging variables for the one or more images, and ii) a dedicated phase variable for a fluctuation of said reference phase; and
output said one or more images.
2. The signal processing system of claim 1 , wherein the projection data has been acquired from different projection directions.
3. The signal processing system of claim 2 , wherein the reference phase fluctuation is modelled by said phase variable as a constant offset independent of said different projection directions.
4. The signal processing system of claim 2 , wherein the reference phase fluctuation is modeled by said phase variable as a non-constant offset that depends on said different projection directions.
5. The signal processing system of claim 1 , wherein the reference phase fluctuation is modeled by said phase variable to depend on a position of a detector element of said X-ray detector.
6. The signal processing system of claim 5 , wherein the detector element is a single detector pixel or a group of detector pixels.
7. The signal processing system of claim 1 , wherein the phase variable includes a temporal dependency to model a change over time of said fluctuation.
8. An imaging system, comprising:
an X-ray imager comprising an interferometer and an X-ray detector configured to supply interferometric projection data; and
a signal processing system comprising:
at least one processor configured to:
receive the interferometric projection data derived from signals acquired by an X-ray detector, said signals caused by X-ray radiation after interaction of said X-ray radiation with an interferometer and with an object to be imaged, said interferometer having a reference phase;
reconstruct one or more images of a spatial distribution of one or more physical properties of said object and fit said interferometric projection data to a signal model by adapting a plurality of fitting variables, said fitting variables including i) one or more imaging variables for the one or more images, and ii) a dedicated phase variable for a fluctuation of said reference phase; and
output said one or more images.
9. A signal processing method, comprising:
receiving interferometric projection data derived from signals acquired by a detector, said signals caused by radiation after interaction of said radiation with an interferometer and an object to be imaged, said interferometer having a reference phase;
reconstructing one or more images of a spatial distribution of one or more physical properties of said object by fitting, said interferometric projection data to a signal model by adapting a plurality of fitting variables, said fitting variables including i) one or more imaging variables for the one or more images, and ii) a dedicated phase variable for a fluctuation of said reference phase; and
outputting said one or more images.
10. The method according to claim 9 , wherein said physical property includes at least one of: i) attenuation, ii) refraction or electron density distribution, and iii) small angle scattering.
11. A non-transitory computer-readable medium having one or more executable instructions stored thereon which, when executed by at least one processor, cause the at least one processor to perform a signal processing method, comprising:
receiving interferometric projection data derived from signals acquired by a detector, said signals caused by radiation after interaction of said radiation with an interferometer and an object to be imaged, said interferometer having a reference phase;
reconstructing one or more images of a spatial distribution of one or more physical properties of said object by fitting said interferometric projection data to a signal model by adapting a plurality of fitting variables, said fitting variables including i) one or more imaging variables for the one or more images, and ii) a dedicated phase variable for a fluctuation of said reference phase; and
outputting said one or more images.