An X-ray machine with gratings takes not one but three pictures at once: a regular image, a map of beam deflections invisible to the eye, and an image of tiny structures that scatter light. Together, they help spot, for example, emphysema or pulmonary fibrosis. But there's a catch: to separate these three layers, you need a clever mathematical trick. For decades, the same one was used — fitting the data to a sine wave, as if laying a wave with the right rhythm over the image. The problem is that any scratch on the grating or a weak signal throws off that rhythm, and stripes appear on the image, like ripples from a stone thrown into water.
The authors of the new study took as their starting point an equation that describes how light, passing through matter, slightly changes direction and scatters — much like a flashlight beam in fog. This is the Fokker–Planck equation, and it can account for the random jolts each photon receives. Instead of fitting a sine wave, the algorithm looks at how smoothly brightness changes from pixel to pixel and figures out on its own where the real shadow is and where the junk from grating defects is. It's as if the image is wiped with a soft cloth, removing dust but leaving important details untouched.
This isn't just a cosmetic improvement. When an X-ray is taken quickly (for example, to avoid asking the patient to hold their breath for too long) the signal weakens, and noise can hide a tiny nodule. The Fokker–Planck-based algorithm extracts more useful information from such noise than the usual sine wave. It's like an experienced restorer who, from a few surviving brushstrokes, reconstructs the entire painting without inventing extra details.
as written by the authors arXiv:2605.01265
Grating interferometry is a promising diagnostic technique that enables simultaneous acquisition of three complementary, synergistic X-ray images: transmission, differential phase, and dark-field. Its key advantage over other setups is its ability to use large pixels and, hence, large-area detectors, as well as its compatibility with low-coherence, compact X-ray sources, both of which are key factors for human-scale imaging. It has already demonstrated strong potential for chest imaging applications, including the diagnosis of pulmonary emphysema, fibrosis, and cancer. To retrieve transmission, differential phase, and dark-field images from data, an algorithm is required to separate the distinct mechanisms contributing to measured contrast. Since its realization, this image-retrieval step has remained fundamentally unchanged. In this work, we develop a novel transmission- and dark-field retrieval algorithm for grating-interferometry derived from the X-ray Fokker-Planck equation. To demonstrate and validate our Fokker-Planck algorithm, we apply it to experimental measurements of a test sample and to data from a mouse chest acquired with varying exposure times and added Poisson noise. The retrieved images were qualitatively and quantitatively compared with those retrieved using a conventional sinusoidal-fitting approach. Across both samples, the Fokker--Planck method produced images consistent with conventional retrieval, with a comparable signal-to-noise ratio. Notably, our Fokker-Planck method suppresses artefacts arising in the conventional approach under grating perturbations (e.g., structural defects like scratches) and reduced flux or visibility, yielding smoother and more reproducible images. Additionally, we demonstrate that our Fokker-Planck method has an advantage over the conventional dark-field retrieval method for fast sample imaging with short exposure times and high noise.
The Fokker–Planck equation, which underlies the method, was originally devised to describe Brownian motion: those chaotic jumps of pollen in water that botanist Robert Brown once observed under a microscope.
In science fiction stories about future medical scanners, the idea of a device that sees through a person without any distortion often appears: this algorithm takes a step in exactly that direction.