WO2010089319A1 - Low dose single step grating based x-ray phase contrast imaging - Google Patents

Low dose single step grating based x-ray phase contrast imaging Download PDF

Info

Publication number
WO2010089319A1
WO2010089319A1 PCT/EP2010/051291 EP2010051291W WO2010089319A1 WO 2010089319 A1 WO2010089319 A1 WO 2010089319A1 EP 2010051291 W EP2010051291 W EP 2010051291W WO 2010089319 A1 WO2010089319 A1 WO 2010089319A1
Authority
WO
WIPO (PCT)
Prior art keywords
grating
ray
phase
pixel
detector
Prior art date
Application number
PCT/EP2010/051291
Other languages
French (fr)
Inventor
Peiping Zhu
Ziyu Wu
Marco Stampanoni
Original Assignee
Institute Of High Energy Physics
Paul Scherrer Institut
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute Of High Energy Physics, Paul Scherrer Institut filed Critical Institute Of High Energy Physics
Priority to EP10703260.9A priority Critical patent/EP2400891B1/en
Priority to CA2751442A priority patent/CA2751442C/en
Priority to AU2010210169A priority patent/AU2010210169B2/en
Priority to JP2011548673A priority patent/JP5606455B2/en
Priority to US13/148,198 priority patent/US8972191B2/en
Priority to CN2010800066507A priority patent/CN102325498B/en
Publication of WO2010089319A1 publication Critical patent/WO2010089319A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/484Diagnostic techniques involving phase contrast X-ray imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/041Phase-contrast imaging, e.g. using grating interferometers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/42Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating

Abstract

Phase sensitive X-ray imaging methods can provide substantially increased contrast over conventional absorption based imaging, and therefore new and otherwise inaccessible information. The use of gratings as optical elements in hard X-ray phase imaging overcomes some of the problems that have impaired the wider use of phase contrast in X-ray radiography and tomography. So far, to separate the phase information from other contributions detected with a grating interferometer, a phase-stepping approach has been considered, which implies the acquisition of multiple radiographic projections. Here, an innovative, highly sensitive X-ray tomographic phase contrast imaging approach is presented based on grating interferometry, which extracts the phase contrast signal without the need of phase stepping. Compared to the existing phase step approach, the main advantage of this new method dubbed "reverse projection" is the significantly reduced delivered dose, without degradation of the image quality. The new technique sets the pre-requisites for future fast and low dose phase contrast imaging methods, fundamental for imaging biological specimens and in-vivo studies.

Description

Low dose single step grating based X-ray phase contrast imaging
The present invention relates of a method and a system for low dose single step grating based X-ray phase contrast imaging .
It is well known that, differently from conventional visible light optics, the refractive index in X-ray optics is very close to and smaller than unity since the X-ray photon energy is often much larger than the atomic resonance energies. In first approximation, for small and negligible anisotropy in the medium, the index of refraction characterizing the optical properties of a tissue can be expressed - including X-ray absorption - with its complex form: n=l-δ-iβ where δ is the decrement of the real part of the refractive index, characterizing the phase shifting property, while the imaginary part β describes the absorption property of the sample. In conventional absorption-based radiography, the X-ray phase shift information is usually not directly utilized for image reconstruction. However, at photon energies greater than 10 keV and for light materials (made up of low-Z elements) , the phase shift term plays a more prominent role than the attenuation term because δ is typically three orders of magnitude larger than β. As a consequence, phase-contrast modalities can generate significantly greater image contrast compared to conventional, absorption-based imaging. Furthermore, far from absorption edges, δ is inversely proportional to the square of the X-ray energy whilst β decreases as the fourth power of energy. A significant consequence of this mechanism is that phase signals can be obtained with much lower dose deposition than absorption, a very important issue when radiation damage has to be taken into account such as in biological samples or in living systems. Several approaches have been developed in order to record the phase signal. They can be classified as interferometric methods (with crystals) , phase propagation methods, techniques based on an analyzer crystal or on grating interferometry.
In the prior art the feasibility of two-grating interferometry in the hard X-ray region using a pair of transmission gratings made by gold stripes on glass plates has been demonstrated. This work has since been extended to achieve a three-dimensional tomographic phase reconstruction using a hard X-ray two-gratings interferometer. Recently, three-grating interferometry in the hard X-ray region with low-brilliance tube-based X-ray sources has been demonstrated. This laboratory-based instrument is of great interest for applications in biology, medicine and for non-destructive testing. A grating interferometer setup is mechanically robust, is easy to align, has low sensitivity to mechanical drift and its requirements on temporal coherence (ΔE/E ~ 0.1-0.2) and spatial coherence (few microns) are moderate: as a consequence the instrument can be easily scaled up to large fields of view, an important asset when used in combination with a conventional X-ray tube.
These characteristics make grating interferometry superior to other phase contrast approaches and set the prerequisites for a broad use of phase contrast X-ray radiography and tomography.
To separate the phase information from other contributions, a phase-stepping approach is normally adopted. One of the gratings is displaced transversely to the incident beam whilst acquiring multiple projections. The intensity signal at each pixel in the detector plane oscillates as a function of the displacement and the phase of this intensity oscillation can be directly linked to the wave- front phase profile and to the decrement of the real part δ of the object's refractive index.
Obviously, this approach is loaded with the limitation of both (long) data acquisition time and severe dose released to specimen.
It is therefore the objective of the present invention to provide a method and a system for the extraction of the phase information, which does not require a stepping procedure, thus overcoming limitations of both data acquisition time and dose imparted to the specimen.
This objective is achieved according to the present invention with respect to the system by an interferometer for x-rays, in particular hard x-rays, for obtaining quantitative x-ray images from a sample including: a) an X-ray source; b) a diffractive optical element, hereafter referred to as beam splitter grating, other than a Bragg crystal, preferably in transmission geometry; c) a position-sensitive detector with spatially modulated detection sensitivity having a number of individual pixels; d) means for recording the images of the detector; e) means for evaluating the intensities for each pixel in a series of images in order to identify the characteristic of the object for each individual pixel as an absorption dominated pixel and/or a differential phase contrast dominated pixel and/or an x-ray scattering dominated pixel; wherein the series of images is collected by continuously or stepwise rotating from 0 to π or 2π either the sample or the interferometer and the source relative to the sample.
With respect to the method the objective is achieved according to the present invention by a method for obtaining quantitative x-ray images from a sample, comprising the steps of: a) providing an X-ray source; b) providing a diffractive optical element, hereafter referred to as beam splitter grating, other than a Bragg crystal, preferably in transmission geometry; c) providing a position-sensitive detector with spatially modulated detection sensitivity having a number of individual pixels; d) applying shots of the X-ray source to the probe and recording the images of the detector; e) evaluating the intensities for each pixel in a series of images in order to identify the characteristic of the object for each individual pixel as an absorption dominated pixel and/or a differential phase contrast dominated pixel and/or an x-ray scattering dominated pixel; wherein the series of images is collected by continuously or stepwise rotating from 0 to 2π either the sample or the interferometer relative to the sample.
The invented system and method therefore present an innovative, highly sensitive X-ray tomographic phase contrast imaging approach based on grating interferometry, which extracts the phase contrast signal without the need of phase stepping (PS) . Compared to the existing phase step approach, the main advantage of this invention dubbed "reverse projection (RP)" is the significantly reduced delivered dose, without degradation of the image quality. The new technique sets the pre-requisites for future fast and low dose phase contrast imaging methods, fundamental for imaging biological specimens and in-vivo studies.
Typically, the beam splitter grating may be a line grating, preferably a phase grating; that is, a grating with low X- ray absorption, but considerable X-ray phase shift (Φ), the latter preferably of either Φe f (2/ -I)- -arcsinθ.8, (2/ -I)- + arcsinθ.8 1 or
Φe ((2/ -l)π - arcsin0.8, (2/ -l)π + arcsinθ.8) , where / = 1,2,3 A further preferred embodiment of the present invention may provide the phase grating that acts as the beam splitter is made by deep etching into silicon, a polymer or similar material.
A further preferred embodiment of the present invention may provide the analyzer grating with one-dimensional grating structure being integrated into the detector, the pixel of the detector is in range of 2 to 10 times the size of the period of the grating, half lines with sensor in a pixel are sensitive to X-ray and half lines without sensor let X- ray go through. In this way the analyzer grating with 100% absorption can be achieved without needing to make heavy metal absorption gratings with high aspect ratio, in particular it is possible to avoid gold gratings.
A further preferred embodiment of the present invention may provide an analyzer grating having a one-dimensional grating structure with high X-ray absorption contrast, its period is the same as that of the self image of the phase grating, placed closely in front of the detector, with its lines parallel to those of the phase grating; preferably this analyzer grating serves as an anti-scatter grid, or an anti-scatter grid is used as a modulation mask.
Dimensioning the interferometer is fundamental for the present invention. Advantageously, the distance between the beam splitter grating and the analyzer grating is chosen to be an odd fractional Talbot distance, given by the equation
_ L - Dn L - n - A 2 /2η 2λ
Dn,spk = -r-p- = - 2 * 2. • whe re « = U, 5 , and
L - Dn L -n - P1 /2η λ ' where
Figure imgf000007_0001
/ = 1,2,3 , Dn is an odd fractional Talbot distance when the parallel X-ray beam is used, while Dnsph is that when the fan or cone X-ray beam is used, L is the distance between the source and the phase grating.
Further, the position of half slope on the shifting curve may be achieved by positioning at least one of the beam splitter grating and the analyzer grating relative to the probe in a direction substantially perpendicular to the orientation of the lines in at least one of the two gratings .
In order to establish a rather simple set-up of the present interferometer, a mechanism can be comprised to place the sample to be investigated between the source and the beam splitter grating or between the beam splitter grating and the analyzer grating being rotated from 0 to n or to 2n.
A further preferred embodiment of the present invention may provide a collimator placed between the source and the beam splitter grating limiting the spatial extent of the illuminating X-rays to a fan beam; a line-array detector is used, and a mechanism is comprised that allows to rotate (either stepwise or continuously) the sample relative to the rest of the apparatus, the rotational axis being perpendicular to the opening angle of the fan, and preferably at same time allows to translate (either stepwise or continuously) the sample relative to the rest of the apparatus along the direction parallel to the rotational axis. Alternatively, a collimator placed between the source and the beam splitter grating may limit the spatial extent of the illuminating X-rays to a cone beam, a pixel-array detector is used, and a mechanism is comprised that allows to rotate the sample relative to the rest of the apparatus, perpendicular to the opening angle of the fan.
Excellent results with respect to the quality of the image can be achieved when an analysis procedure is implemented for reverse-projection data that comprises the steps of calculating, for each element of the detector, the absorption signal M and the refraction angle θr according to the following equations (8) and (9) resp.:
Figure imgf000008_0001
Preferred embodiments of the present invention are hereinafter described in more detail thereby referring to the attached drawings .
Brief description of the drawings:
Figure 1 (a) shows the working principle of the grating interferometer: through the Talbot effect, a periodic interference pattern (known as self image) is formed behind the phase grating (Gl), in the plane of the analyzer grating (G2) . Figure 1 (b) is a plot of the intensity oscillation (shifting curve) as a function of the grating position x for a detector pixel over one period of the analyzer grating. The dots corresponds to the measured values (normalized to unit) while the gray line shows a sinusoidal fit.
Figure 2 illustrates phase contrast tomographic reconstructions of a demineralised mouse joint, acquired at a voxel size of 3.5 X 3.5 X 3.5 μm3. Sub-Figures al to a3 show the data obtained with the classical phase stepping (PS) protocol, while Sub-Figures bl to b3 the reconstruction using the reverse projection (RP) method, al and bl shows an axial slice: bl is sharper than al and there are no ring artifacts (see text below) . a2 and b2 depict a coronal slice through the joint, clearly showing that the RP protocol is less sensitive to typical horizontal stripes artifacts observed with the PS method (see enlarged inset) . a3 and b3 show a sagittal view through the joint. The dotted lines mark the locations where the axial views (al and bl) have been taken. Scale bar is 500 microns.
Figure 3 represents a phase contrast reconstructed coronal slice of a rat brain, obtained after tomographic reconstruction using the PS- (a) and the RP-protocol (b) . Qualitatively, both reconstructions are very similar. In (b) the effects of the grating imperfection (ring artifacts), as expected, are more evident. Figure (c) shows a quantitative comparison of two line profiles extracted at the position marked by colour bars (hippocampus region) . Scale bar is 1 mm.
Figure 4 (a) shows differential phase contrast radiography of a rat paw (7 stacks, RP-protocol) . Figures 4(bl-2) and (cl-2) show axial and coronal slices through the paw acquired with the PS and RP protocol, respectively. Structural details of both soft (muscles, fat) and hard tissue (bone) are well visible. Scale bars are 2 mm in (a) and 1 mm for (bl-2) and (cl-2) . Figure 5(al-2) is a tomographic reconstruction of a rat brain - (al-2) obtained with the PS protocol, (bl-2) obtained with the RP-protocol using Eq. 11 to calculate the map of the index of refraction. Scale bar is 1 mm in (al,bl) and 2 mm in (a2,b2) .
Table 1 summarizes the experimental parameters for the tomographic scans of the three investigated samples: a rat brain (4% PFA, paraffin embedded), a (demineralized) mouse joint in PBS and a rat paw (4% PFA) . All experiments have been performed at 25 keV and at the 3rd Talbot distance. Visibility of the interferometer was ~ 30%.
With reference to the above-mentioned figures, an innovative approach for the extraction of the phase information is presented which does not require a stepping procedure, thus overcoming limitations of both data acquisition time and dose released to specimens.
This novel approach relies on the physical similarities between a crystal analyzer based system and a grating interferometer. Both techniques record refraction angle signals and, analogously to the rocking curve of a crystal analyzer, the properties of the shifting curve (see Figure 1) can be exploited to fully describe the performance of a grating interferometer. The refraction angle, i.e., the phase information of the sample, can be extracted by setting the grating interferometer in the central position where the intensity curve follows a linear behavior.
According to the aforementioned analogy, the intensity I recorded by a detector positioned after the grating interferometer can be expressed as: x
/ = I0 exp -^ μ(x,y,z)dyr
D [1]
where μ is the linear absorption coefficient, xg denotes the relative displacement between the phase grating and the analyzer grating along the direction perpendicular to both the incoming beam and the line of gratings, θr is the refraction angle, D is the distance between the phase and the analyzer grating, S is the shifting curve. For the
Figure imgf000011_0001
sake of simplicity, the scattering contribution - which would induce a weak increment of the background noise - is neglected. (xr,yr,z) are the coordinates of the reference frame associated to the X-ray beam and (x,y,z) those associated with the sample. The two frames are linked by the rotation matrix
Figure imgf000011_0002
being φ the rotation angle between the xr and the x-axis around the z-axis.
With a good approximation, the behaviour of the shifting curve near its half slope may be considered linear so that, if P2 is the period of the analyzer grating, θr<-^- can be replaced by a first-order Taylor expansion. Further we can write :
Figure imgf000011_0003
where Q is a constant .
Figure imgf000011_0004
The refraction angle in the X-Y plane (Fig. 1 (a) ) is determined by the line integral of the first-order derivative of the refractive index along the X-ray path and it may be written as:
Figure imgf000012_0001
where δ corresponds to the decrement of the real part of the refractive index of the sample as mentioned in the introduction. Substituting Eqs . 3 and 4 into Eq. 1, the projected image for a grating interferometer can be described by:
Figure imgf000012_0002
μ is a scalar and therefore rotational-invariant, while
3δ strongly depends on the direction along which it is dxr measured.
The projected image at the rotation angles φ and its corresponding reverse image at φ+π can be written as:
/(xr,φ,z) = /o exp -j μ(x,j,z)dyr \s U M- i1--cC?j -9δ(w) '-ά-y [6] dxr
I(-xr ,φ +π , z) = I0 exp -jμteΛI», ls(iϊ1+c/»|ϋ≤φ. [7]
The absorption signal can be obtained from the two projected images by the sum of Eq. 6 and 7 and solving the Beer-Lambert relationship, i.e.,
Figure imgf000012_0003
In the same way, the angle of refraction can be obtained by a proper combination of Eq. 6 and 7, as shown in the following expression:
Figure imgf000013_0001
According to fundamentals of computed tomography reconstruction, i.e., the Fourier Slice Theorem, the absorption coefficient as well as the refractive index can be easily obtained by the inverse Fourier Transform and a Hubert filter: μ(x,y,z) = Jc*|> J [M(xr,<b,z)*F-1 (|p|)] -δ(xcosφ +
Figure imgf000013_0002
-xr)dxr [10]
0 δ (x,y,z) = + jvsinφ -xr)dxr [11]
Figure imgf000013_0003
where p is the spatial frequency and F"1 denotes the inverse Fourier transform.
Based on Eq. 10 and 11, we introduce here a novel acquisition protocol, dubbed "reverse projection" (RP) method.
It can be described in five steps:
(i) without sample, scan the phase grating or the analyzer grating along the transverse direction xg over one period of the analyzer grating and record the normalized intensity — on the
detector versus the angle — , i.e., get the
shifting curve S — , (ii) set the grating interferometer at the center of the linear region of the shifting curve by positioning the phase grating or the analyzer grating at xg = p214 or xg =—p2/ 4 , (iϋ) put the sample in front of or behind the phase grating, collect m angular projections of the sample over a rotation of 360°, (iv) extract M and θr according to Eq. 8 and Eq. 9 and finally (v) reconstruct either the absorption coefficient or the refractive index using the filtered back- projection. Therefore, the total number of acquired projection images is m.
On the contrary, the Phase Stepping (PS) acquisition protocol can be described in four steps:
(i) put the sample in front of or behind the phase grating, scan one of the two gratings along the transverse direction xg (k points over one period of the analyzer grating) and record one projection image at each point, (ii) repeat step (i) for a total of m/2 times over a sample rotation of 180°,
(iii) extract the gradient signal via the Fourier analysis of the intensity signal and
(iv) reconstruct the phase via a filtered back- projection .
For this second method, the total number of acquired projections images is k*m/2.
As a consequence, the total number of projections required by the RP protocol is reduced by a factor of k/2 compared to the PS. The method was validated by performing both phase stepping (PS) and reverse projection (RP) experiments using the grating interferometer installed at the TOMCAT beamline of the Swiss Light Source at the Paul Scherrer Institute, Villigen, Switzerland. The energy was tuned at 25 keV and the interferometer was operated in the 3rd Talbot distance. In this configuration, the visibility has been measured to be 30%. Additional details on the grating interferometer installed at TOMCAT can be found in public documentation related to this installation with the Paul Scherrer Institute .
In a first case study, we investigated two different samples: a rat brain first fixed in 4% paraformaldehyde (PFA) and then embedded in paraffin and a demineralized mouse joint, fixed only in a phosphate buffer solution (PBS) (no embedding) . We used the mouse joint and the rat brain to test the reconstruction method both on small (< 4 mm) and large (> 10 mm) samples. Both specimens are weakly absorbing objects and therefore ideal candidates for phase contrast imaging. Reconstructions based on the two methods (PS and RP) are shown in Fig. 2 and 3 while experimental parameters are summarized in Table 1.
Fig. 2 shows axial, sagittal and coronal views of a mouse joint obtained with both PS and RP protocols (Table 1) . The joint was immersed and fixed in an Eppendorf vial containing PBS to avoid any movements during the acquisition. A qualitative comparison of the images clearly shows that RP-reconstructions are comparable to those obtained with the PS approach. Moreover, looking at the inset shown in Figure (2,a2) and (2,b2) the RP-slice appears to be sharper than the PS-reconstructions. This can be explained by the fact that the shifting curve is directly proportional to the refraction angle and that this - in the RP protocol - is obtained by simple subtraction of a reference image (with no sample) from the paired images described in Eq. 9. In addition, since with the RP method no phase-stepping is required, the system is less sensitive to mechanical instabilities.
The largest investigated sample, a rat brain, was mounted vertically on the flat surface of the sample support to match the horizontal field of view of the detector system. The vertical sample arrangement also enabled a direct reconstruction of coronal slices through the sample, an approach very useful when trying to identify anatomical brain regions (Fig. 3) . The height of the sample was larger than the vertical height of the beam and therefore four scans have been collected along the vertical direction to image the whole brain. To achieve phase matching between sample and surroundings, we used an aquarium bath filled with room temperature liquid paraffin (chemical formula CriΑ2ri+2 where n = 5-17, density ~ 0.7 g cm"3) . For large samples too, a qualitative comparison of the images clearly shows that the RP-reconstruction is as good as the one obtained with the PS-approach. In addition, a line profile taken at the level of the hippocampus, see Fig. 3c, shows a quantitative good agreement between RP and PS approaches.
In the second case, the novel method has been validated using a more realistic sample, namely a specimen containing both soft and hard tissue. For this purpose, we investigated a rat paw (containing both bone and muscles) which was only fixed in 4% PFA. This fixation procedure is frequently used to maintain biological samples in a status as close as possible to their natural, original conditions. The rat paw was also mounted vertically in order to best match the horizontal field of view of the detector. Seven stacked scans were necessary to image the full sample volume . The measurement of the rat paw was the most challenging experiment since the sample has been measured in air. This usually causes large phase jumps at the air-specimen interface and explains the "star" artifacts visible in Fig. 4b-l and, less serious, in 4c-l. This is because the shifting curve is saturated when θr≥ — Pl— and, as a consequence, the RP method is not very sensitive to large refraction angles. This is not the case for the PS method, which has to cope with angles as large as θr < —— .
Our invention introduces a novel approach for fast and low dose extraction of both the absorption coefficient and the refractive index of a sample using a grating interferometer is introduced. It is demonstrated that this new approach yields comparable information to the established phase stepping technique but with a reduction factor of k/2 in the total dose delivered to the sample. Moreover, the reverse projection approach makes high-sensitivity phase contrast computed tomography (CT) as straightforward as conventional, absorption based CT. It is first shown that this new method works well with parallel beam geometries but it is not difficult to generalize it to either cone or fan beam setups, making it accessible also to X-ray tube- based applications.
In particular, the significant decrease of the dose and the straight forward acquisition protocol does no affect image quality, while representing a major advancement in hard X- ray phase contrast tomography for synchrotron radiation and laboratory X-ray sources, enabling experiments impossible so far.
The next and probably most challenging application of the RP-protocol will be in-vivo phase contrast imaging. With the advent of new, high efficient and high speed detectors it will be possible to acquire the same amount of data within a fraction of a second. We estimate that it will be realistic to obtain a full tomographic data set with the RP protocol with a total exposure time of 2-3 seconds. This, together with the ongoing efforts regarding robust and reliable iterative reconstruction algorithms, requiring a significant smaller amount of projections, can push the total acquisition time below Is and hence opening up the possibility of phase contrast tomographic microscopy of small living animals.
Another very challenging application of the RP-protocol will be the quantitative 3D description of the scattering signal. This image contrast is generated by small-angle scattering within the sample and it provides complementary and otherwise inaccessible structural information at micrometer and sub-micrometer length scale. However, the signal is not rotational-invariant and therefore it will be very challenging to quantitatively obtain such information in 3D.
Further developments will concern the manufacturing of optimized gratings for high X-ray energies leading to the implementation of the RP-protocol in new medical X-ray CT scanners that would offer a significant increase in soft tissue sensitivity, a characteristic now provided (at much lower resolutions however) only by much more expensive techniques such as magnetic resonance imaging.
Finally, we would like to point out that this approach is not limited to X-ray imaging and may be easily generalized to other methods such as grating based neutron phase imaging and visible light differential interference contrast (DIC) microscopy where a similar shifting curve is considered and a quantitative phase description appears possible . Table 1
Figure imgf000019_0001
1 The optical system always produces field of views with a square shape However, along the vertical direction the size of the field of view is reduced and limited by the vertical size of the beam (approximately 3 6 mm at 25 keV and at 25 m from the source for the TOMCAT beamhne of the Swiss Light Source) SUPPLEMENTARY MATERIAL I
Ring-like artefacts due to grating imperfections are discussed hereinafter.
A careful study of the reconstructed images reveals that the reverse projections (RP) protocol imposes more stringent requirements on gratings, i.e., it needs better full-field uniformity and lower local grating imperfections than gratings suitable for the phase stepping (PS) method. In particular, imperfections may induce small ring-like artifacts in the RP-reconstructed slices, which are less evident or missing in PS-reconstructions.
Ring artifacts are clearly visible in Figure 5, for both coronal (bl) and sagittal (b2) cuts. Due to the averaging effect associated to the phase stepping extraction, the PS- protocol is significantly less sensitive to grating defects and therefore the rings artifact are less pronounced, see Figure 5 al-2.

Claims

Claims
1. An interferometer for x-rays, in particular hard x-rays, for obtaining quantitative x-ray images from a sample including: a) a X-ray source (x-ray) b) a diffractive optical element, hereafter referred to as beam splitter grating (Gl), other than a Bragg crystal, preferably in transmission geometry; c) a position-sensitive detector (PSD) with spatially modulated detection sensitivity having a number of individual pixels; d) means for recording the images of the detector (PSD) ; e) means for evaluating the intensities for each pixel in a series of images in order to identify the characteristic of the object for each individual pixel as an absorption dominated pixel and/or a differential phase contrast dominated pixel and/or an x-ray scattering dominated pixel; wherein the series of images is collected by continuously or stepwise rotating from 0 to π or 2π either the sample or the interferometer and the source relative to the sample.
2. The interferometer according to claim 1, wherein the beam splitter grating is a line grating (Gl), preferably a phase grating; that is, a grating with low X-ray absorption, but considerable X-ray phase shift (Φ), the latter preferably of either Φe | (2/ -l)--arcsin0.8, (2/ -l)- + arcsin0.8 ] or
Φe ((2/ -l)π - arcsin0.8, (2/ -l)π + arcsinθ.8) , where / = 1,2,3
3. The interferometer according to claim 1 or 2, wherein the phase grating that acts as the beam splitter grating is made by deep etching into silicon, a polymer or similar material.
4. The interferometer according to claim 1, wherein an analyzer grating having a one-dimensional grating structure (G2) with high X-ray absorption contrast, its period is the same as that of the self image of the phase grating, is placed closely in front of the detector (PSD) with its lines parallel to those of the phase grating; preferably this grating structure serves as an anti-scatter grid, or an anti-scatter grid is used as the modulation mask.
5. The interferometer according to claim 1, wherein the analyzer grating with one-dimensional grating structure is integrated into the detector, the pixel of the detector is from 2 to 10 times the size of the period of the grating, half lines with sensor in a pixel are sensitive to X-ray and half lines without sensor let X- ray go through.
6. The interferometer according to any of the preceding claims, wherein the distance (D) between the beam splitter grating (Gl) and the analyzer grating (G2) is chosen to be an odd fractional Talbot distance, given by the equation
n,sph where « = 1,3,5 , and
Figure imgf000022_0001
' where
Figure imgf000022_0002
/ = 1,2,3 , Dn is an odd fractional Talbot distance when the parallel X-ray beam is used, while Dnsph is that when the fan or cone X-ray beam is used, L is the distance between the source and the phase grating.
7. The interferometer according to any of the preceding claims, wherein a mechanism is comprised to place the sample to be investigated between the source (X-ray) and the beam splitter grating (Gl) or between the beam splitter grating (Gl) and the analyzer grating (G2) being rotated from 0 to π or to 2π.
8. The interferometer according to claim 1, wherein a collimator placed between the source and the beam splitter grating (Gl) limits the spatial extent of the illuminating X-rays to a fan beam, a line-array detector is used, and a mechanism is comprised that allows to rotate (either stepwise or continuously) the sample relative to the rest of the apparatus, the rotational axis being perpendicular to the opening angle of the fan, and preferably at the same time allows to translate (either stepwise or continuously) the sample relative to the rest of the apparatus along the direction parallel to the rotational axis.
9. The interferometer according to claim 1, wherein a collimator placed between the source and the beam splitter grating (Gl) limits the spatial extent of the illuminating X-rays to a cone beam, a pixel-array detector is used, and a mechanism is comprised that allows to rotate the sample relative to the rest of the apparatus, perpendicular to the opening angle of the fan .
10. The interferometer according to any of the preceding claims, wherein an analysis procedure is implemented for reverse-projection data that comprises the steps of calculating, for each element of the detector, the absorption signal M and the refraction angle θr according to the following equations (8) and (9) resp.: )
Figure imgf000024_0001
1 7(xr,φ,z) -/(-xr,φ +π ,z) p θδ(x,j;,z) o:
C /(χ ,φ,z) + /(-x ,φ +π ,z) r Ψ J 3x
11. A method for obtaining quantitative x-ray images from an sample, comprising the steps of: a) providing an X-ray source (x-ray) ; b) providing a diffractive optical element, hereafter referred to as beam splitter grating (Gl) , other than a Bragg crystal, preferably in transmission geometry; c) providing a position-sensitive detector (PSD) with spatially modulated detection sensitivity having a number of individual pixels; d) applying shots of the X-ray source to the probe and recording the images of the detector (PSD) ; e) evaluating the intensities for each pixel in a series of images in order to identify the characteristic of the object for each individual pixel as an absorption dominated pixel and/or a differential phase contrast dominated pixel and/or an x-ray scattering dominated pixel; wherein collecting the series of images by continuously or stepwise rotating from 0 to 2n either the sample or the interferometer relative to the sample.
12. The method according to claim 11, wherein the beam splitter grating is a line grating (Gl) , preferably a phase grating; that is, a grating with low X-ray absorption, but considerable X-ray phase shift (Φ), the latter preferably of either Φef (2/-I)--arcsinθ.8, (2/-I)-+ arcsinθ.81 or
Φe((2/-l)π-arcsin0.8, (2/-l)π +arcsinθ.8) , where / = 1,2,3
13. The method according to claim 11 or 12, wherein the phase grating that acts as the beam splitter grating is made by deep etching into silicon, a polymer or similar material .
14. The method according to claim 11, wherein an analyzer grating having a one-dimensional grating structure (G2) with high X-ray absorption contrast, its period is the same as that of the self image of the phase grating, is placed closely in front of the detector (PSD) with its lines parallel to those of the phase grating; preferably this grating structure serves as an anti-scatter grid, or an anti-scatter grid is used as the modulation mask.
15. The interferometer according to any of the preceding claims, wherein the distance between the beam splitter and the analyzer is chosen to be an odd fractional Talbot distance, given by the equation
, where , = 1,3,5 and
Figure imgf000025_0001
π L + D n,sph if the phase shift of G1 is (2/ -I)- , p2 = η = i L+D^ph Pl ' where if the phase shift of G1 is (2/ -l)π , p2 =
L 2
/ = 1,2,3 , Dn is an odd fractional Talbot distance when the parallel X-ray beam is used, while Dnsph is that when the fan or cone X-ray beam is used, L is the distance between the source and the phase grating.
16. The method according to any of the preceding claims 11 to 15, wherein the position of half slope on the shifting curve is achieved by positioning at least one of the gratings (Gl, G2) relative to the probe in a direction (xg) substantially perpendicular to both the incoming beam and the orientation of the lines in at least one grating (Gl, G2) .
17. The method according to any of the preceding claims 11 to 16, wherein a mechanism is comprised to place the sample to be investigated between the source (X-ray) and the beam splitter grating (Gl) or between the beam splitter grating (Gl) and the analyzer grating (G2) being rotated from 0 to π or to 2π.
18. The method according to claim 11, wherein a collimator placed between the source and the beam splitter grating
(Gl) limits the spatial extent of the illuminating X- rays to a fan beam, a line-array detector is used, and a mechanism is comprised that allows to rotate (either stepwise or continuously) the sample relative to the rest of the apparatus, the rotational axis being perpendicular to the opening angle of the fan, and preferably at the same time allows to translate (either stepwise or continuously) the sample relative to the rest of the apparatus along the direction parallel to the rotational axis.
19. The method according to claim 11, wherein a collimator placed between the source and the beam splitter grating
(Gl) limits the spatial extent of the illuminating X- rays to a cone beam, a pixel-array detector is used, and a mechanism is comprised that allows to rotate the sample relative to the rest of the apparatus, perpendicular to the opening angle of the fan.
20. The method according to any of the preceding claims 11 to 19, wherein an analysis procedure is implemented for reverse-projection data that comprises the steps of calculating, for each element of the detector, the absorption signal M and the refraction angle θr according to the following equations (8) and (9) resp.:
Figure imgf000027_0001
21. The method according to any of the preceding claims 11 to 20, wherein the analyzer grating with one-dimensional grating structure is integrated into the detector, the pixel of the detector is from 2 to 10 times the size of the period of the grating, half lines with sensor in a pixel are sensitive to X-ray and half lines without sensor let X-ray go through.
PCT/EP2010/051291 2009-02-05 2010-02-03 Low dose single step grating based x-ray phase contrast imaging WO2010089319A1 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
EP10703260.9A EP2400891B1 (en) 2009-02-05 2010-02-03 Low dose single step grating based x-ray phase contrast imaging
CA2751442A CA2751442C (en) 2009-02-05 2010-02-03 Low dose single step grating based x-ray phase contrast imaging
AU2010210169A AU2010210169B2 (en) 2009-02-05 2010-02-03 Low dose single step grating based X-ray phase contrast imaging
JP2011548673A JP5606455B2 (en) 2009-02-05 2010-02-03 Imaging apparatus for backprojection and method of operating the same
US13/148,198 US8972191B2 (en) 2009-02-05 2010-02-03 Low dose single step grating based X-ray phase contrast imaging
CN2010800066507A CN102325498B (en) 2009-02-05 2010-02-03 Low dose single step grating based X-ray phase contrast imaging

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP09100099 2009-02-05
EP09100099.2 2009-02-05

Publications (1)

Publication Number Publication Date
WO2010089319A1 true WO2010089319A1 (en) 2010-08-12

Family

ID=42078937

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2010/051291 WO2010089319A1 (en) 2009-02-05 2010-02-03 Low dose single step grating based x-ray phase contrast imaging

Country Status (7)

Country Link
US (1) US8972191B2 (en)
EP (1) EP2400891B1 (en)
JP (1) JP5606455B2 (en)
CN (1) CN102325498B (en)
AU (1) AU2010210169B2 (en)
CA (1) CA2751442C (en)
WO (1) WO2010089319A1 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012080125A1 (en) 2010-12-13 2012-06-21 Paul Scherrer Institut A method and a system for image integration using constrained optimization for phase contrast imaging with an arrangement of gratings
JP2012130586A (en) * 2010-12-22 2012-07-12 Fujifilm Corp X-ray image detecting apparatus, radiographing apparatus, and radiographing system
WO2013014083A1 (en) 2011-07-28 2013-01-31 Paul Scherrer Institut Method for image fusion based on principal component analysis
WO2013124164A1 (en) * 2012-02-24 2013-08-29 Paul Scherrer Institut A system for non-invasively classification of different types of micro-calcifications in human tissue
WO2014148828A1 (en) * 2013-03-20 2014-09-25 삼성전자 주식회사 Method for projection and back-projection for image processing, and image processing apparatus therefor
CN110095481A (en) * 2019-05-24 2019-08-06 清华大学 X-ray grating imaging system and imaging method
CN112568923A (en) * 2020-12-10 2021-03-30 中国科学院深圳先进技术研究院 X-ray phase contrast image extraction method, device, terminal and storage medium
CN114137002B (en) * 2021-11-18 2023-07-14 北京航空航天大学 Low-dose X-ray differential phase contrast imaging method based on contrast enhancement

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012157690A (en) * 2011-01-14 2012-08-23 Fujifilm Corp Radiation image capturing apparatus and radiation image detecting device
RU2585801C2 (en) 2011-06-01 2016-06-10 Тоталь Са Device for x-ray tomography
WO2012164092A1 (en) * 2011-06-01 2012-12-06 Total Sa An x-ray tomography device
US20150117599A1 (en) 2013-10-31 2015-04-30 Sigray, Inc. X-ray interferometric imaging system
WO2013111050A1 (en) * 2012-01-24 2013-08-01 Koninklijke Philips N.V. Multi-directional phase contrast x-ray imaging
US9357975B2 (en) * 2013-12-30 2016-06-07 Carestream Health, Inc. Large FOV phase contrast imaging based on detuned configuration including acquisition and reconstruction techniques
US10096098B2 (en) 2013-12-30 2018-10-09 Carestream Health, Inc. Phase retrieval from differential phase contrast imaging
US10578563B2 (en) 2012-12-21 2020-03-03 Carestream Health, Inc. Phase contrast imaging computed tomography scanner
WO2014154188A1 (en) * 2013-03-26 2014-10-02 Institute Of Experimental And Applied Physics Method of phase gradient radiography and arrangement of an imaging system for application of the method
CN103164863B (en) * 2013-04-02 2016-03-02 中国科学院高能物理研究所 For rebuilding the method for positron emission computerized tomography image
WO2014180901A1 (en) 2013-05-10 2014-11-13 Siemens Aktiengesellschaft Phase contrast x-ray imaging device and refracting grid therefor
WO2014187885A1 (en) 2013-05-22 2014-11-27 Siemens Aktiengesellschaft Phase-contrast x-ray imaging device
US10297359B2 (en) 2013-09-19 2019-05-21 Sigray, Inc. X-ray illumination system with multiple target microstructures
US10269528B2 (en) 2013-09-19 2019-04-23 Sigray, Inc. Diverging X-ray sources using linear accumulation
US10295485B2 (en) 2013-12-05 2019-05-21 Sigray, Inc. X-ray transmission spectrometer system
US9719947B2 (en) * 2013-10-31 2017-08-01 Sigray, Inc. X-ray interferometric imaging system
USRE48612E1 (en) * 2013-10-31 2021-06-29 Sigray, Inc. X-ray interferometric imaging system
US10304580B2 (en) 2013-10-31 2019-05-28 Sigray, Inc. Talbot X-ray microscope
JP2015118081A (en) * 2013-11-12 2015-06-25 キヤノン株式会社 Radiation detection system and radiation imaging device
WO2015122542A1 (en) * 2014-02-14 2015-08-20 Canon Kabushiki Kaisha X-ray talbot interferometer and x-ray talbot interferometer system
CN106104318B (en) * 2014-04-09 2019-05-03 拉姆伯斯公司 Low-power image modification detector
CN103900502B (en) * 2014-04-16 2017-01-04 中国科学技术大学 Precise displacement measurement device based on X-ray geometric projection Moire fringe and method
US10401309B2 (en) 2014-05-15 2019-09-03 Sigray, Inc. X-ray techniques using structured illumination
CN107076682B (en) * 2014-05-15 2021-02-02 斯格瑞公司 X-ray method for measuring, characterizing and analyzing periodic structures
JP2016106721A (en) * 2014-12-03 2016-06-20 キヤノン株式会社 Image processing device and image processing method
US10352880B2 (en) 2015-04-29 2019-07-16 Sigray, Inc. Method and apparatus for x-ray microscopy
EP3314576B1 (en) * 2015-06-26 2019-11-27 Koninklijke Philips N.V. Robust reconstruction for dark-field and phase contrast ct
US10295486B2 (en) 2015-08-18 2019-05-21 Sigray, Inc. Detector for X-rays with high spatial and high spectral resolution
CN105675631A (en) * 2016-01-05 2016-06-15 合肥泰禾光电科技股份有限公司 Fast fan-beam geometric phase contrast CT imaging device and method
CN107807139B (en) * 2016-09-05 2020-04-24 天津工业大学 Dual-energy X-ray phase contrast imaging system without stepping device and implementation method thereof
US10247683B2 (en) 2016-12-03 2019-04-02 Sigray, Inc. Material measurement techniques using multiple X-ray micro-beams
CN106618623B (en) * 2017-01-11 2019-08-30 合肥工业大学 The imaging method of the hard X ray grating interferometer of single exposure
US10598612B2 (en) 2017-02-01 2020-03-24 Washington University Single-shot method for edge illumination X-ray phase-contrast tomography
WO2018175570A1 (en) 2017-03-22 2018-09-27 Sigray, Inc. Method of performing x-ray spectroscopy and x-ray absorption spectrometer system
JP6943090B2 (en) * 2017-09-05 2021-09-29 株式会社島津製作所 X-ray imaging device
EP3494885A1 (en) * 2017-12-07 2019-06-12 Koninklijke Philips N.V. Apparatus for presentation of dark field x-ray image information
US10578566B2 (en) 2018-04-03 2020-03-03 Sigray, Inc. X-ray emission spectrometer system
US10845491B2 (en) 2018-06-04 2020-11-24 Sigray, Inc. Energy-resolving x-ray detection system
GB2591630B (en) 2018-07-26 2023-05-24 Sigray Inc High brightness x-ray reflection source
US10656105B2 (en) 2018-08-06 2020-05-19 Sigray, Inc. Talbot-lau x-ray source and interferometric system
DE112019004433T5 (en) 2018-09-04 2021-05-20 Sigray, Inc. SYSTEM AND PROCEDURE FOR X-RAY FLUORESCENCE WITH FILTERING
CN112823280A (en) 2018-09-07 2021-05-18 斯格瑞公司 System and method for depth-selectable X-ray analysis
WO2021046059A1 (en) 2019-09-03 2021-03-11 Sigray, Inc. System and method for computed laminography x-ray fluorescence imaging
CN110916713B (en) * 2019-11-29 2022-04-29 清华大学 Grating imaging system and scanning method thereof
CN110916712B (en) * 2019-11-29 2022-04-29 清华大学 Grating imaging system and scanning method thereof
NL2024368B1 (en) * 2019-12-03 2021-08-31 Xeikon Prepress Nv Method and system for processing a raster image file
US11175243B1 (en) 2020-02-06 2021-11-16 Sigray, Inc. X-ray dark-field in-line inspection for semiconductor samples
JP7395775B2 (en) 2020-05-18 2023-12-11 シグレイ、インコーポレイテッド Systems and methods for X-ray absorption spectroscopy using a crystal analyzer and multiple detector elements
JP2023542674A (en) 2020-09-17 2023-10-11 シグレイ、インコーポレイテッド System and method for depth-resolved measurement and analysis using X-rays
WO2022126071A1 (en) 2020-12-07 2022-06-16 Sigray, Inc. High throughput 3d x-ray imaging system using a transmission x-ray source
CN113237901A (en) * 2021-05-07 2021-08-10 中国科学院上海应用物理研究所 Biological feature recognition system and biological feature recognition method
US11885755B2 (en) 2022-05-02 2024-01-30 Sigray, Inc. X-ray sequential array wavelength dispersive spectrometer
CN115684222B (en) * 2022-12-21 2023-04-11 济南汉江光电科技有限公司 Rapid low-dose X-ray multi-mode CT system and imaging method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1731099A1 (en) * 2005-06-06 2006-12-13 Paul Scherrer Institut Interferometer for quantitative phase contrast imaging and tomography with an incoherent polychromatic x-ray source
US20070183579A1 (en) * 2006-02-01 2007-08-09 Joachim Baumann X-ray optical transmission grating of a focus-detector arrangement of an X-ray apparatus for generating projective or tomographic phase contrast recordings of a subject
US20070183581A1 (en) * 2006-02-01 2007-08-09 Bjorn Heismann Focus/detector system of an X-ray apparatus for generating phase contrast recordings

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000098449A (en) 1998-09-21 2000-04-07 Canon Inc Exposure controller, exposure control method and electronic camera
JP2003135438A (en) 2001-11-02 2003-05-13 Fuji Photo Film Co Ltd Method, device, and program for image generation
JP2004203066A (en) 2002-12-24 2004-07-22 Kobe Steel Ltd Bumper device of automobile
US8841046B2 (en) 2004-10-22 2014-09-23 Eulitha Ag System and a method for generating periodic and/or quasi-periodic pattern on a sample
DE102006063048B3 (en) 2006-02-01 2018-03-29 Siemens Healthcare Gmbh Focus / detector system of an X-ray apparatus for producing phase-contrast images
CN101011256A (en) 2006-02-01 2007-08-08 西门子公司 Method and measuring arrangement for nondestructive analysis of an examination object by means of X-radiation
EP1879020A1 (en) 2006-07-12 2008-01-16 Paul Scherrer Institut X-ray interferometer for phase contrast imaging
JP5273955B2 (en) * 2007-06-26 2013-08-28 株式会社日立製作所 X-ray imaging apparatus and X-ray imaging method
US20100327175A1 (en) * 2007-12-14 2010-12-30 Yakov Nesterets Phase-contrast imaging method and apparatus
US7949095B2 (en) * 2009-03-02 2011-05-24 University Of Rochester Methods and apparatus for differential phase-contrast fan beam CT, cone-beam CT and hybrid cone-beam CT
US8760817B2 (en) * 2009-05-22 2014-06-24 HGST Netherlands B.V. Three-terminal design for spin accumulation magnetic sensor

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1731099A1 (en) * 2005-06-06 2006-12-13 Paul Scherrer Institut Interferometer for quantitative phase contrast imaging and tomography with an incoherent polychromatic x-ray source
US20070183579A1 (en) * 2006-02-01 2007-08-09 Joachim Baumann X-ray optical transmission grating of a focus-detector arrangement of an X-ray apparatus for generating projective or tomographic phase contrast recordings of a subject
US20070183581A1 (en) * 2006-02-01 2007-08-09 Bjorn Heismann Focus/detector system of an X-ray apparatus for generating phase contrast recordings

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012080125A1 (en) 2010-12-13 2012-06-21 Paul Scherrer Institut A method and a system for image integration using constrained optimization for phase contrast imaging with an arrangement of gratings
JP2014506352A (en) * 2010-12-13 2014-03-13 パウル・シェラー・インスティトゥート Method and system for image integration using constrained optimization for phase contrast imaging with a grating device
JP2012130586A (en) * 2010-12-22 2012-07-12 Fujifilm Corp X-ray image detecting apparatus, radiographing apparatus, and radiographing system
US9117296B2 (en) 2011-07-28 2015-08-25 Paul Scherrer Institut Method for image fusion based on principal component analysis
WO2013014083A1 (en) 2011-07-28 2013-01-31 Paul Scherrer Institut Method for image fusion based on principal component analysis
CN103858002A (en) * 2011-07-28 2014-06-11 保罗·谢勒学院 Method for image fusion based on principal component analysis
JP2014521945A (en) * 2011-07-28 2014-08-28 パウル・シェラー・インスティトゥート Image fusion method based on principal component analysis
WO2013124164A1 (en) * 2012-02-24 2013-08-29 Paul Scherrer Institut A system for non-invasively classification of different types of micro-calcifications in human tissue
WO2014148828A1 (en) * 2013-03-20 2014-09-25 삼성전자 주식회사 Method for projection and back-projection for image processing, and image processing apparatus therefor
US9993218B2 (en) 2013-03-20 2018-06-12 Samsung Electronics Co., Ltd. Method for projection and back-projection for image processing, and image processing apparatus therefor
CN110095481A (en) * 2019-05-24 2019-08-06 清华大学 X-ray grating imaging system and imaging method
CN112568923A (en) * 2020-12-10 2021-03-30 中国科学院深圳先进技术研究院 X-ray phase contrast image extraction method, device, terminal and storage medium
CN114137002B (en) * 2021-11-18 2023-07-14 北京航空航天大学 Low-dose X-ray differential phase contrast imaging method based on contrast enhancement

Also Published As

Publication number Publication date
AU2010210169A1 (en) 2011-08-25
EP2400891B1 (en) 2017-07-19
CA2751442A1 (en) 2010-08-12
JP2012516739A (en) 2012-07-26
CA2751442C (en) 2018-06-12
US8972191B2 (en) 2015-03-03
EP2400891A1 (en) 2012-01-04
CN102325498B (en) 2013-07-10
CN102325498A (en) 2012-01-18
AU2010210169B2 (en) 2015-04-09
JP5606455B2 (en) 2014-10-15
US20120041679A1 (en) 2012-02-16

Similar Documents

Publication Publication Date Title
CA2751442C (en) Low dose single step grating based x-ray phase contrast imaging
EP2652708B1 (en) A method and a system for image integration using constrained optimization for phase contrast imaging with an arrangement of gratings
US10058300B2 (en) Large FOV phase contrast imaging based on detuned configuration including acquisition and reconstruction techniques
US9795350B2 (en) Material differentiation with phase contrast imaging
US8855265B2 (en) Correction method for differential phase contrast imaging
Mokso et al. Advantages of phase retrieval for fast x-ray tomographic microscopy
Vila-Comamala et al. High sensitivity X-ray phase contrast imaging by laboratory grating-based interferometry at high Talbot order geometry
JP2013541699A (en) Differential phase contrast imaging with improved sampling
JP2016535642A (en) X-ray grating phase contrast imaging apparatus and method
Stampanoni et al. Tomographic hard X‐ray phase contrast micro‐and nano‐imaging at TOMCAT
EP3140641A1 (en) Tilted-grating approach for scanning-mode x-ray grating interferometry
JP5143147B2 (en) X-ray imaging apparatus and X-ray imaging method
Thüring et al. Compact hard x-ray grating interferometry for table top phase contrast micro CT
JP6788521B2 (en) X-ray imaging device
Pyakurel Phase and dark field radiography and CT with mesh-based structured illumination and polycapillary optics
Zdora et al. Principles and State of the Art of X-ray Speckle-Based Imaging
Bopp et al. X-ray Phase Contrast: Research on a Future Imaging Modality
Pyakurel et al. Optimization of Signal and Noise in X-Ray Phase and Dark Field Imaging with a Wire Mesh
Wolf et al. Effect of object size, position, and detector pixel size on X-ray absorption, differential phase-contrast and dark-field signal
Han et al. Developments of x-ray grating imaging and trying of multiple information fusion
Li et al. Developments of X-Ray Grating Imaging Technology
Li et al. Measurements of imaging parameters of a phase contrast x-ray imaging prototype

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201080006650.7

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 10703260

Country of ref document: EP

Kind code of ref document: A1

DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)
REEP Request for entry into the european phase

Ref document number: 2010703260

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2010703260

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2751442

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2011548673

Country of ref document: JP

ENP Entry into the national phase

Ref document number: 2010210169

Country of ref document: AU

Date of ref document: 20100203

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 13148198

Country of ref document: US