Open Access
20 February 2020 Phase-contrast x-ray tomography of neuronal tissue at laboratory sources with submicron resolution
Marina Eckermann, Mareike Töpperwien, Anna-Lena Robisch, Franziska van der Meer, Christine Stadelmann, Tim Salditt
Author Affiliations +

Purpose: Recently, progress has been achieved in implementing phase-contrast tomography of soft biological tissues at laboratory sources. This opens up opportunities for three-dimensional (3-D) histology based on x-ray computed tomography (μ- and nanoCT) in the direct vicinity of hospitals and biomedical research institutions. Combining advanced x-ray generation and detection techniques with phase reconstruction algorithms, 3-D histology can be obtained even of unstained tissue of the central nervous system, as shown, for example, for biopsies and autopsies of human cerebellum. Depending on the setup, i.e., source, detector, and geometric parameters, laboratory-based tomography can be implemented at very different sizes and length scales.

We investigate the extent to which 3-D histology of neuronal tissue can exploit the cone-beam geometry at high magnification M using a nanofocus transmission x-ray tube (nanotube) with a 300 nm minimal spot size (Excillum), combined with a single-photon counting camera. Tightly approaching the source spot with the biopsy punch, we achieve high M  ≈  101  −  102, high flux density, and exploit the superior efficiency of this detector technology.

Approach: Different nanotube configurations such as spot size and flux, M, as well as exposure time, Fresnel number, and coherence are varied and selected in view of resolution, field of view, and/or phase-contrast requirements.

Results: The data show that the information content for the cytoarchitecture is enhanced by the phase effect. Comparison of results to those obtained at a microfocus rotating-anode x-ray tomography setup with a high-resolution detector, i.e., in low-M geometry, reveals similar to slightly superior data quality for the nanotube setup. In addition to its compactness, reduced power consumption by a factor of 103, and shorter scan duration, the particular advantage of the nanotube setup also lies in its suitability for pixel detector technology, enabling an increased range of opportunities for applications in laboratory phase-contrast x-ray tomography.

Conclusions: The phase retrieval scheme utilized mixes amplitude and phase contrast, with results being robust with respect to reconstruction parameters. Structural information content is comparable to slightly superior to previous results achieved with a microfocus rotating-anode setup but can be obtained in shorter scan time. Beyond advantages as compactness, lowered power consumption, and flexibility, the nanotube setup’s scalability in view of the progress in pixel detector technology is particularly beneficial. Further progress is thus likely to bring 3-D virtual histology to the performance in scan time and throughput required for clinical practice in neuropathology.



Phase-contrast x-ray tomography offers a unique potential to realize 3-D virtual histology, with cellular and even subcellular resolution, and for 3-D volumes, which are inaccessible by more established techniques. While volume and throughput in conventional histology are limited by slicing and staining, volume penetration for both light- and electron-based microscopy techniques is unsuitable for larger tissue volumes. Several propagation-based x-ray phase-contrast tomography studies with synchrotron radiation have demonstrated this potential, both for stained13 and unstained 4,5 soft tissues, and including mouse models6 as well as human tissue.5,7 Many different types of tissues have been imaged, including heart and cardiovascular systems,810 skin,11 cancerous tissue in particular for mamma,12 the peripheral nervous system,1,13 and tissues of the central nervous system.7,1416

Toward broader accessibility and use of the technique in the clinical setting, translation from synchrotron to laboratory sources is an important goal of ongoing technique and instrumentation development. Laboratory-based phase-contrast computed tomography (μCT) has already been implemented at different sources: transmission microfocus sealed tubes,17 liquid metal-jet anodes,16,1821 as well as microfocus rotating-anode sources.22,23 The detection technology has been equally diverse. Nanotubes offer yet another opportunity to implement phase-contrast tomography with possible advantages, in particular, in view of increased spatial coherence length. For stained tissue, high-resolution tomography of different tissues has been demonstrated in Refs. 24 and 25. A particular challenge is to reach sufficient image quality for unstained soft tissues with laboratory radiation. To this end, we have recently demonstrated the capabilities of optimized image acquisition (by geometry and detection) and reconstruction, using liquid-anode metal-jet sources as well as microfocus rotating-anode sources.23,26

In this work, we want to investigate the suitability of a home-built setup installed at a nanotube (NanoTube N1 60 kV, Excillum AB, Stockholm, Sweden) for 3-D histology of unstained “postmortem” human brain tissue. As in our earlier study,5 we chose cerebellar tissue as a reference structure for phase-contrast tomography data evaluation, due to its well-known anatomy and features covering various relevant length scales and e densities. Figure 1 illustrates its basic anatomical features in order to place the tomographic results obtained on biopsy punches from paraffin-embedded tissue blocks with typical diameters of a millimeter into proper anatomical perspective. The “white matter” (WM) consists of myelinated neuronal structures as “axons” in “fiber tracts” (FT). From there, fibers traverse the “granular layer” (GL), consisting of “granular cells” (GC) and “Golgi cells” (GgC). Bounding the GL, a single layer of “Purkinje cells” (PC, PCL) with sparse occurrences of “basket cells” (BC) is situated. PC’s long, planar, and widely branched dendritic trees reach into the “molecular layer” (ML), which distinguishes itself from the GL by the “molecular cells” (MC), having more sparsely spread neurons and having larger nuclei.

Fig. 1

Depiction of the cerebellum anatomy, giving a rather macroscopic perspective on the left and a zoom-in to the cellular scale on the right. Notations correspond to the ones defined in the text.


The paper is organized as follows. After this introduction, the setup realization and data processing of the nanotube x-ray source are outlined in Sec. 2.1, complemented by a brief presentation of a well-established rotating-anode microfocus setup in Sec. 2.2. The results (Sec. 3) start with the inspection of commissioning test measurements regarding stability and spectrum in Sec. 3.1. However, the analysis of the tomographic data occupies the main part of this work, being devoted to the variation of physical parameters as source spot size, lateral coherence length, dose, or Fresnel number in Secs. 3.2.1 and 3.2.2 and also to the comparison with a well-established microfocus setup in Sec. 3.2.3. This paper closes with a summary and outlook in Sec. 4.




NanoTube Implementation


Experimental setup

We have designed a nanofocus x-ray setup for propagation-based x-ray phase-contrast tomography, as depicted in Fig. 2. Its central component is the NanoTube N1 60 kV (Excillum AB) x-ray source, with two-dimensional (2-D) spatial resolution down to 150 nm according to the manufacturer (lines-and-spaces, metal test objects in absorption contrast). The source is operated at 60 keV with a power between 0.2 and 1.2 W, depending on the spot size. The NanoTube system was calibrated for three different spot sizes s [full width at half maximum, (FWHM)]: “big spot, high flux” at 1  μm, “middle spot, middle flux” at 0.5  μm, and “small spot, low flux” at 0.3  μm (FWHM). Figure 2(a) shows the schematic of the transmission-anode target, which consists of a 0.50  μm W-film providing a rather high e stopping power S (S=54.4  MeV/cm27) on a 100-μm layer of diamond, serving as carrier layer and for heat mitigation at fewer e interactions (S=18.1  MeV/cm). As shown in the zoomed-in optical micrograph in Fig. 2(c), samples can be positioned in direct proximity to the target, since it also serves as a vacuum exit window. The cone-shaped front end of the NanoTube allows to position the sample tower underneath the x-ray target and to achieve small source-to-sample distances down to z01100  μm. Owing to the small submicron spot size s, high geometrical magnification M=z02/z011 can be realized without source blurring, making it possible to use direct (photon-counting) pixel detectors, which are available only with relatively large pixel size px. Accordingly, the effective pixel size in the sample plane is reduced to28

Eq. (1)


Fig. 2

Illustration of the nanotube setup: (a) schematic of the transmission target x-ray source; (b) recording geometry: the cone-shaped front end of the nanotube allows for small distances z01 between source s and the sample. The detector is aligned at z02 from the source. The axis camera allows for visual inspection of the setup in general, while the Manta camera gives a microscopic control of the tight z01 environment; see exemplary view in (c).


At the same time, phase contrast is still feasible, since small s also assures the spatial coherence length ζ29 to be of the same order of magnitude as pxeff

Eq. (2)

where λ denotes the x-ray wavelength. The sample tower is equipped with three translational motors at the bottom (x,y,z) to position the tomographic rotation axis. The field of view (FOV) on the sample is selected by two further translational motors (x,z) on top of the rotation, enabling automated alignment routines, as described in Ref. 23. For high-resolution optical monitoring of the sample environment, a Manta camera (Allied Vision Technologies GmbH, Stadtroda, Germany) is installed with a view along the x axis, as shown in Figs. 2(b) and 2(c). Broad overviews of the scene are provided by an Axis camera (Axis Communications AB, Lund, Sweden). A single-photon counting Timepix Hexa H05-W0154 detector (XIE, Freiburg, Germany) was used with a 500-μm Si sensor, px=55  μm and 768×512  pixels (w×h) distributed over six modules to record the 2-D projections.


Data recording settings

For tomographic acquisition, 2-D x-ray projections were recorded at 1201 rotation angles, equally distributed over 192 deg (due to the cone-beam angle of 12 deg, an increased angular sampling range was required), with sets of 25 flat-field images before and after each scan. The z0220  cm was kept fixed, varying z01 and s. The detector was operated with a lower-limit cut-off energy of 4 keV. The acquisition time was adjusted depending on the spot size s to avoid overexposure, which is at maximum counts of 11810 ph. Most scans were split into four consecutive tomographic scans with 1/4 exposure time each, and then recombined via cross correlation in Fourier space.30 Thus, for these cases, four different exposures were recorded and averaged for each projection. All scan parameters are tabulated in Tables 1Table 23. This acquisition scheme was found to reduce ring artifacts arising from a statistically varying response of the modules.

Table 1

Parameters and quality measured for the tomographic scans at fixed F. Exposure times also indicate whether the scan was performed as a single or multiple ones. FSC is based on volumes of 4003  voxels, with a Kaiser–Bessel window of 7 pixels.31 The SNR was calculated as (μft−μbg)/σbg on five GC nuclei, four PC bodies, and their nuclei, respectively. Scan time refers to the full scan, including readout and motor positioning.

Experimental parameters
s (μm)Fz01 (mm)pxeff (μm)ζ (μm)Exp. time (s)Power (W)Scan time (h)Dose (kGy)
Quality measures
s (μm)FαBACγBACFSC (μm)SNR—GC nucleusSNR—PC nucleusSNR—PC body

Table 2

Analysis of image quality as a function of s and F. Exposure times also indicate whether the scan was performed as a single or multiple ones. FSC evaluation is based on volumes of K3 voxels, with a Kaiser–Bessel window of 7 pixels.31 SNR was calculated as (μft−μbg)/σbg on five GC nuclei, four PC bodies, and their nuclei, respectively.

Experimental parameters
s (μm)Fz01 (mm)pxeff (μm)ζ (μm)Exp. time (s)Power (W)Scan time (h)Dose (kGy)
Image quality parameters
s (μm)FαBACγBACFSC (μm)SNR—GC nucleusSNR—PC nucleusSNR—PC body

Table 3

Analysis of the tomographic scans for nanofocus setup with pxeff≈1.07  μm, in comparison to a rotating-anode scan with the same pxeff. Exposure times also indicate whether the scan was performed as a single or multiple ones. FSC analysis is based on volumes of 4003   voxels, with a Kaiser–Bessel window of 7 pixels.31 SNR was calculated as (μft−μbg)/σbg on five GC nuclei, four PC bodies, and their nuclei, respectively.

Experimental parameters
s (μm)Fz01 (mm)pxeff (μm)ζ (μm)Exp. time (s)Power (W)Scan time (h)
Image quality parameters
s (μm)FαBACγBACFSC (μm)SNR—GC nucleusSNR—PC nucleusSNR—PC body


Phase retrieval

The data are recorded in the “direct contrast” or “edge enhancement regime” at Fresnel numbers:

Eq. (3)

where pxeff=px/M for the effective pixel size and zeff=z12/M have been used for the effective propagation distance (Fresnel scaling theorem28). As demonstrated in Refs. 20 and 32, in this regime, even data from low coherence sources can be successfully reconstructed based on the “Bronnikov-aided correction” (BAC) scheme. To this end, the “transport of intensity equation” (TIE) serves as a starting point to describe the propagation of a paraxial wave along z with intensity I(r) and phase distribution ϕ(r): [I(r)·ϕ(r)]=kzI(r), with k being the wavenumber. Under the assumption of small propagation distances zeff and a purely phase-shifting object, the TIE can be linearized, and an approximate phase ϕ˜ can be computed as33,34

Eq. (4)

where F denotes the Fourier transform, I0 denotes the (uniform) intensity distribution of the illumination, and k is the spatial frequency in units of px1. The parameter α is introduced to regularize the singularity at zero spatial frequencies, and in practice, is chosen such that edge enhancement is cancelled. In a second step, the approximate phase ϕ˜ is used to compute a (corrected) sharp intensity distribution in the object exit plane, according to35

Eq. (5)


In practice, the relevance of using the BAC is to choose the two regularization parameters α and γ such that the resulting somewhat hypothetical intensity distribution contains contributions from intensity and phase contrast, even though expressed only in terms of transmitted intensity. As for reconstructions assuming a homogeneous object, the two contributions cannot be separated. For high-resolution laboratory μ- and nano-CTs, this scheme yields unparalleled image quality. Afterward, the phase-retrieved 2-D projection data are processed for ring artifact mitigation according to Ref. 36 and then recombined to a 3-D volume using the ASTRA-toolbox.37,38 Visualization was done with Avizo (Thermo Fisher Scientific, Waltham, Massachusetts).


Rotating-Anode Setup

For comparison, the samples were also scanned at a (home-built) laboratory μCT setup installed at a rotating Cu-anode x-ray source (Rigaku, Tokyo, Japan) with main line 8.048 keV (Kα) and source size s=70  μm.23 It was operated at 40 keV and 30 mA. A high resolution, lens-coupled single crystal scintillator CCD detector (Xsight, Rigaku, Prague, Czech Republic), resulting in px=0.54  μm was used. Images were recorded at M=500  mm/494  mm1 to achieve a spatial coherence ζ1  μm reasonable for micron resolution (after 2×2  pixel-binning).

The number of projection angles was identical to that of the nanotube scans, distributed over an angular range of 180 deg from the cone-beam angle of 0  deg. The exposure time, however, needed to be increased significantly due to the very different detector technology, i.e., to 50 s. In addition to the flat fields, 10 dark images were recorded prior to each scan to account for the dark current of the CCD. With F=1.24, BAC also applies here.




Commissioning Tests

Prior to performing high-resolution tomography scans, the stability of the source and its environment was verified. To this end, the source was operated at spot size s=0.56  μm, and images of a custom-fabricated JIMA target (1.5  μmW; ZonePlates Ltd, London), positioned at z01=19.6  mm, were recorded at z02=900  mm, using an flat panel detectors camera (Photonic Science, UK) with a 5μm Gadox scintillator 2208×2744  pixels (w×h) and px=4.54  μm. As indicated in Fig. 3(a), images of 600 nm lines-and-spaces were acquired every 12 min for 5 min of exposure over the course of 6.6 h. Figure 3(b) shows the modulations of the lines-and-spaces corresponding to the marked area in Fig. 3(a) after 0, 3, and 6 h (from top to bottom). Only a minor shift on the order of 30  nm/h is observed. At the same time, the variations of the total intensity (integrated flat fields) shown in Fig. 3(c) are about 0.6% (RMS value).

Fig. 3

Pre-characterization of the imaging setup. (a)–(c) Test of stability by recording a time series for a test pattern, recorded at s=0.56  μm and pxeff=0.099  μm. (a) Image of the JIMA target (1.5  μmW). The annotation “600” refers to the structure measured, i.e., 600 nm lines and 600 nm spaces. The red rectangle indicates the area for which line profiles are shown in (b) at three time points: (top) 0 h, (center) 3 h, and (bottom) 6 h. (c) Integrated intensity over time. (d) Source spectrum: (red) measured data (detector output), (yellow) corrected for absorption in air, and (blue) the spectrum taking into account the energy dependence in the quantum efficiency of the detector used for the tomographic scans.


Next, the spectrum of the x-ray source was assessed using an energy-resolving XR-100CdTe detector (Amptek, Bedford) with a sensor thickness of 1 mm, positioned at z02=1.2  m to avoid saturation. Spectral bins were calibrated based on the Kα and Kβ fluorescence signals of Mo, Ni, and Ag foils. The spectrum is plotted in red in Fig. 3(d). The counts-weighted mean energy of this spectrum Ew is 21.40 keV. However, at the given distance, absorption in air is already substantial and has to be corrected for using the tabulated values (Henke tables), accessed through the CXRO data base39 and matched to the detector bins with a MATLAB-implementation of shape-preserving piecewise cubic extrapolation. The resulted corrected curve is shown in yellow, representing the source emission spectrum, with a correspondingly lower mean energy Ew=11.46  keV. This spectrum was further corrected for the X-ray absorption in air corresponding to the detector position z0220  cm and for the detector sensitivity: the sensor material was 500  μm of Si, hence exhibiting a stronger sensitivity for low-energetic photons. The final, weighted spectrum is given in blue. It is characterized by Ew=9.11  keV and therefore, represents the mean photon energy for the tomography data (without taking beam hardening in the sample into account).

The photon flux is dependent on s and the actual, specific source calibration. For s=0.95  μm, it was found to be on the order of 2·1011  ph/s in 2π space, 4·1010  ph/s for s=0.51  μm, and 1·1010  ph/s for s=0.30  μm.


Phase-Contrast Tomography of Unstained Human Brain Tissue

In this work, paraffin-embedded human cerebellum was used to investigate the suitability of the setup for 3-D virtual histology and neuropathology. From the tissue block, a 0.5-mm biopsy punch was extracted and transferred into a 0.5-mm polyimide tube (Professional Plastics, Fullerton, California) on a custom-fabricated Huber brass pin (Huber Diffraktionstechnik GmbH & Co. KG, Rimsting, Germany), as illustrated in Fig. 2.


Variation of s/ζ at constant F

We have first investigated the influence of the lateral coherence length ζ on the tomographic image quality by setting two different spot sizes (I) s=0.95  μm and (II) s=0.51  μm, respectively, while keeping F2.0 constant, and the dose also approximately at the same level. Hence the larger source size resulted in a significantly reduced total scan time. All experimental and phase reconstruction parameters are listed in Table 1. In this work, the radiation dose was calculated combining the photon counts from the projection data with the spectra in Fig. 3. Figure 4 illustrates the data reconstruction steps from projection to orthoslices through the reconstructed volume. The projections in Fig. 4(a) reveal slight edge enhancement in both cases, especially for the polyimide–air interface, as plotted in Fig. 4(b). As expected, the higher coherence for the 0.501μm source spot yields a more pronounced edge enhancement compared to 0.934  μm. In the same plot, the blue curves show the profiles after phase retrieval, and the respective 2-D images are given in Fig. 4(c). Virtual slices through the same position in the xz plane and xy plane are shown in Figs. 4(d) and 4(f), respectively, with a corresponding zoom in shown in Fig. 4(e). Data for (I) were recorded as a single tomographic scan, acquisition for (II) was split into four as described in Sec. 2.1. By fractionating the dose over four scans, ring artifacts were found to be reduced. Based on visual inspection, both data sets appear to be of very similar quality. Different cerebellum-specific regions, as outlined in Fig. 1, are clearly identifiable: as for the zoom ins (I.e) and (II.e), the cell-dense GL on the right and the ML on the left are separated by the sparsely distributed, bold PCL cells (red arrows). The resolution was quantified exploiting Fourier Shell Correlation (FSC).40,41 In this analysis scheme, two independently recorded Kaiser–Bessel-filtered data are evaluated for their consistency, defined via the intersection of the cross correlation in Fourier space with a threshold criteria of 1/2-bit in this case.42 Therefore, FSC values are not only governed by spatial resolution but also by the overall noise level. In this work, two independent data sets are reached performing the tomographic reconstruction twice, using half of the angular projections each. Hence, FSC-based resolution evaluation serves as an upper limit estimate. For both scans inspected here, the analysis results in similar values of around 1.8  μm. However, case (I) for s=0.51  μm has superior feature contrast (see Table 1).

Fig. 4

Data from tomography scans at fixed Fresnel number of 2.0, at (I) s=0.95  μm (ζ=0.501  μm) and (II) 0.51  μm (ζ=0.934  μm). (a) Respective flat-field-corrected projections and (c) BAC-reconstructed projections. Colored bars indicate the position of the profiles shown in (b). The same virtual slice (d) in the xz plane with a zoom-in in (e) and (f) in the xy plane. Red arrows mark PCs. Data from quantitative analysis are summarized in Table 1. Scale bars: 100  μm.



Variations of F at constant source spot size s

Next, we investigate the influence of F on tomographic image quality, as controlled by the source-to-sample distance z01, while keeping z02 fixed. Hence pxeff [Eq. (3)] and ζ [Eq. (2)] vary accordingly. Most scans were recorded at s=0.51  μm=const. In addition, we include a scan at the minimally achievable s=0.30  μm, as well as at s=0.95  μm. For the latter, the 0.95  μm scan from Table 1 was repeated four times to reach an equal overall scan time of ζ. All parameters are detailed in Table 2. Figure 5 shows the virtual slices along the xy plane through the reconstruction volume in similar positions. By repeating the 2  h-0.95  μm-scan from Table 1 four times, and increasing the dose accordingly, the 3-D resolution (FSC) is increased to 1.00  μm, i.e., close to the voxel size. The signal-to-noise ratio (SNR) also increased for GCs and PC nuclei, but only marginally for PC bodies. To increase the SNR of PC bodies substantially, it was necessary to double ζ (via reduction of s) at constant F (constant z01). Despite the significantly lower dose, which compromised the resolution, the contrast for the rather large cell bodies is higher. As an overall trend, reducing F (via z01) results in increased resolution, which is expected based on higher dose (smaller pxeff). At the same time, despite the rise in dose, the SNRs are lower, indicating that the increase in ζ is more important, in particular for intrinsically low-contrasted features such as PC bodies. This conclusion is also confirmed by comparing the results in Table 2 for s=0.95  μm and s=0.30  μm. In these scans, a roughly similar coherence length ζ0.5 to 0.7  μm results in a constant SNR0.83 to 0.86 for PC bodies, even though the radiation dose deviates by more than a factor of 4.

Fig. 5

Virtual xy slices through the tomographic reconstruction at similar positions in the biopsy punch for spot sizes of (a) 0.95  μm, (b)–(d) 0.51  μm, and (e) 0.30  μm and Fresnel numbers of (a) 1.97, (b) 1.96, (c) 1.40, (d) 0.83, and (e) 0.60. Note the rise in geometrical magnification from (a) and (b) to (e). The respective analysis parameters are given in Table 2. Scale bars: 50  μm.



Comparison with data from microfocus laboratory setup

To put the nanotube results into perspective with earlier implementations of 3-D histology with laboratory x-ray phase-contrast tomography, Fig. 6 and Table 3 present a comparison to a reconstruction obtained at the microfocus rotating-anode x-ray source with instrumentation described in Sec. 2.2. For the comparison, we selected the nanotube data set recorded for s=0.51  μm and F=1.96 in order to achieve similar effective pixel size pxeff1.07  μm and coherence length ζ0.93 to 1.09  μm. The nanotube setup achieves slightly higher resolution (FSC-based) and increased SNRs (for the features considered), but requires only half of the scan time. However, as is directly apparent from Fig. 6(a), this comes at the cost of reduced FOV, reflecting the different detector technologies. For sufficiently narrow samples, this could be compensated by consecutively scanning two volumes and stacking them, which would result in similar overall scan time for both setups. The plots in Fig. 6(b) show profiles across the capillary edges indicated in Fig. 6(a). The two setups give very similar intensity profiles; absorption and edge enhancement are slightly more pronounced for the microfocus rotating-anode data. Finally, Figs. 6(c) and 6(d) show virtual slices through the xz and xy planes to judge image quality by visual inspection.

Fig. 6

Comparison of (I) nanotube tomography with (II) microfocus rotating-anode results. Settings were chosen such that pxeff1.07  μm same for both data sets. (a) Flat-field-corrected projections. The colored bars show the position of the intensity profile plots in (b). Note the identical ranges of the y axis. Virtual slices through the same position in the reconstructed sample volume are shown in (c) for the xz plane and (d) the xy plane. Scale bars: 100  μm.



Summary and Outlook

In summary, we have successfully demonstrated phase-contrast tomography of unstained neuronal tissue using a home-built laboratory setup with a nanofocus x-ray transmission tube (nanotube) and a photon-counting pixel detector. Sufficient image quality for the detection of neurons and hence representation of the cytoarchitecture was achieved. In particular, FSC analysis indicated a resolution of 0.90  μm for a source setting of s=0.30  μm. The phase retrieval and reconstruction scheme presented here mixes amplitude and phase information but is very suitable to visualize the small electron density differences in unstained tissue and hence the cytoarchitecture, for example, of neuronal tissue. Final results can be represented as β or δ up to a factor of k, expressing the amplitude (real-valued transmission function) or the phase shift, respectively. Importantly, the information content for the cytoarchitecture is enhanced by the phase effect, as evidenced by the variation of F (see also Fig. 7). Experimental determination of αBAC, as illustrated in Fig. 7, is unproblematic and results are robust with respect to variations up to about ±40%. Since magnification and cone-beam angle are high, cone-beam tomographic reconstruction according to Ref. 43 is required. Effects of the cone-beam geometry can easily appear in the form of distortions, such as at the edges of the source-facing side of the reconstructed volume [cf. Figs. 6(I.d) and 6(II.d)]. Compared to a microfocus rotating-anode setup, data quality and information content appear to be on a similar to slightly superior level. However, the particular advantage of the nanotube setup (apart from its more compact size and much reduced power consumption) is in its scalability with respect to progress in pixel detector technology. While already performing on the same level for a relatively thin silicon sensor, future replacement by 1-mm-thick Si sensors or even GaAs sensors, along with an increase in the detection panels, may result in significantly reduced scan times. We also note that new pixel detector technology allows for registration of counting with subpixel registration and for photon energy determination,44 opening up an entirely new opportunity for laboratory phase-contrast tomography. Further progress is hence likely to bring virtual 3-D histology to the performance in scan time and throughput required for clinical practice in neuropathology. This, however, would be useful only in combination with automated evaluation of cytoarchitecture as in Ref. 5 and more knowledge on alterations associated with pathologies. Importantly, the demonstrated laboratory-based image quality would already be sufficient for this task.

Fig. 7

Edge enhancement and choice of reconstruction parameters. The plots correspond to the data shown in Fig. 5 and give the image profiles before reconstruction (red) exhibiting pronounced edge enhancement between air and Kapton tubes as a hallmark of phase contrast. In addition, profiles are shown after BAC reconstruction, carried out for three different values of αBAC, to illustrate the process of finding the best-suited parameter value.



An earlier version of this work was presented at SPIE Optics and Photonics 2019 as contribution 11113-21 “Phase contrast x-ray tomography of neuronal tissue at laboratory sources with submicron resolution.” Presentation date: August 12, 2019. It is published in the corresponding conference proceedings (Optics and Photonics 2019: X-Ray, Gamma-Ray, and Particle Technologies). Apart from that, there are no relevant financial interests in the paper and no other potential conflicts of interest to disclose. Human cerebellar tissue was retrieved at routine autopsy in agreement with the ethics committee of the University Medical Center Göttingen. Furthermore, the authors declare no conflict of interest.


We thank Excillum AB for their support with the NanoTube x-ray source, Peter Luley and Matthias Meister for engineering support, and Jasper Frohn for assistance with the microfocus rotating-anode x-ray source. This work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—EXC 2067/1-390729940.



M. Bartels et al., “Myelinated mouse nerves studied by x-ray phase contrast zoom tomography,” J. Struct. Biol., 192 561 –568 (2015). JSBIEM 1047-8477 Google Scholar


M. C. Strotton et al., “Optimising complementary soft tissue synchrotron x-ray microtomography for reversibly-stained central nervous system samples,” Sci. Rep., 8 (1), 12017 (2018). SRCEC3 2045-2322 Google Scholar


M. Töpperwien et al., “Phase-contrast tomography of neuronal tissues: from laboratory-to high resolution synchrotron CT,” Proc. SPIE, 9967 99670T (2016). PSISDG 0277-786X Google Scholar


M. Saccomano et al., “Synchrotron inline phase contrast μCT enables detailed virtual histology of embedded soft-tissue samples with and without staining,” J. Synchrotron Radiat., 25 1153 –1161 (2018). JSYRES 0909-0495 Google Scholar


M. Töpperwien et al., “Three-dimensional virtual histology of human cerebellum by x-ray phase-contrast tomography,” Proc. Natl. Acad. Sci. U. S. A., 115 (27), 6940 –6945 (2018). Google Scholar


M. Töpperwien et al., “Contrast enhancement for visualizing neuronal cytoarchitecture by propagation-based x-ray phase-contrast tomography,” NeuroImage, 199 70 –80 (2019). NEIMEF 1053-8119 Google Scholar


A. Khimchenko et al., “Hard x-ray nanoholotomography: large-scale, label-free, 3D neuroimaging beyond optical limit,” Adv. Sci., 5 1700694 (2018). Google Scholar


P. Garcia-Canadilla et al., “Complex congenital heart disease associated with disordered myocardial architecture in a midtrimester human fetus,” Circ.: Cardiovasc. Imaging, 11 e007753 (2018). Google Scholar


W. Vågberg et al., “Cellular-resolution 3D virtual histology of human coronary arteries using x-ray phase tomography,” Sci. Rep., 8 (1), 11014 (2018). SRCEC3 2045-2322 Google Scholar


H. Dejea et al., “Microstructural analysis of cardiac endomyocardial biopsies with synchrotron radiation-based x-ray phase contrast imaging,” Functional Imaging and Modelling of the Heart, 23 –31 Springer, Cham, Switzerland (2017). Google Scholar


L. A. Walton et al., “Morphological characterisation of unstained and intact tissue micro-architecture by x-ray computed micro- and nano-tomography,” Sci. Rep., 5 10074 (2015). SRCEC3 2045-2322 Google Scholar


Y. Zhao et al., “High-resolution, low-dose phase contrast x-ray tomography for 3D diagnosis of human breast cancers,” Proc. Natl. Acad. Sci. U. S. A., 109 18290 –18294 (2012). Google Scholar


J. F. Scopel et al., “Are human peripheral nerves sensitive to x-ray imaging?,” PLoS One, 10 e0116831 (2015). POLNCL 1932-6203 Google Scholar


S. E. Hieber et al., “Tomographic brain imaging with nucleolar detail and automatic cell counting,” Sci. Rep., 6 32156 (2016). SRCEC3 2045-2322 Google Scholar


A. Khimchenko et al., “Extending two-dimensional histology into the third dimension through conventional micro computed tomography,” NeuroImage, 139 26 –36 (2016). NEIMEF 1053-8119 Google Scholar


M. Krenkel et al., “X-ray phase contrast tomography from whole organ down to single cells,” Proc. SPIE, 9212 92120R (2014). PSISDG 0277-786X Google Scholar


F. Pfeiffer et al., “Phase retrieval and differential phase-contrast imaging with low-brilliance x-ray sources,” Nat. Phys., 2 258 –261 (2006). NPAHAX 1745-2473 Google Scholar


T. Tuohimaa, M. Otendal and H. M. Hertz, “Phase-contrast x-ray imaging with a liquid-metal-jet-anode microfocus source,” Appl. Phys. Lett., 91 (7), 074104 (2007). APPLAB 0003-6951 Google Scholar


T. Zhou et al., “Comparison of two x-ray phase-contrast imaging methods with a microfocus source,” Opt. Express, 21 30183 –30195 (2013). OPEXFF 1094-4087 Google Scholar


M. Bartels et al., “Phase contrast tomography of the mouse cochlea at microfocus x-ray sources,” Appl. Phys. Lett., 103 (8), 083703 (2013). APPLAB 0003-6951 Google Scholar


W. Vågberg et al., “X-ray phase-contrast tomography for high-spatial-resolution zebrafish muscle imaging,” Sci. Rep., 5 (16625), 1 –6 (2015). SRCEC3 2045-2322 Google Scholar


A. Tkachuk et al., “X-ray computed tomography in Zernike phase contrast mode at 8 keV with 50-nm resolution using Cu rotating anode x-ray source,” Z. Kristallogr. Cryst. Mater., 222 650 –655 (2007). Google Scholar


M. Reichardt et al., “Nanoscale holographic tomography of heart tissue with x-ray waveguide optics,” Proc. SPIE, 10391 1039105 (2017). PSISDG 0277-786X Google Scholar


M. Müller et al., “Myoanatomy of the velvet worm leg revealed by laboratory-based nanofocus x-ray source tomography,” Proc. Natl. Acad. Sci. U. S. A., 114 12378 –12383 (2017). Google Scholar


M. Busse et al., “Three-dimensional virtual histology enabled through cytoplasm-specific x-ray stain for microscopic and nanoscopic computed tomography,” Proc. Natl. Acad. Sci. U. S. A., 115 (10), 2293 –2298 (2018). Google Scholar


M. Töpperwien et al., “Multiscale x-ray phase-contrast tomography in a mouse model of transient focal cerebral ischemia,” Biomed. Opt. Express, 10 92 –103 (2019). BOEICL 2156-7085 Google Scholar


“NIST chemistry webbook, NIST standard reference database number 69,” National Institute of Standards and Technology, Gaithersburg, Maryland (1998). Google Scholar


D. M. Paganin, Coherent X-Ray Optics, Oxford University, New York (2006). Google Scholar


H. H. Barrett and K. J. Myers, Foundations of Image Science, Wiley(2003). Google Scholar


M. Guizar-Sicairos, S. T. Thurman and J. R. Fienup, “Efficient subpixel image registration algorithms,” Opt. Lett., 33 156 –158 (2008). OPLEDP 0146-9592 Google Scholar


T. Butz, Fouriertransformation für Fußgänger, Vieweg+Teubner Verlag / Springer Fachmedien Wiesbaden GmbH, Wiesbaden (2012). Google Scholar


M. Töpperwien et al., “Three-dimensional mouse brain cytoarchitecture revealed by laboratory-based x-ray phase-contrast tomography,” Sci. Rep., 7 42847 (2017). SRCEC3 2045-2322 Google Scholar


A. Groso, R. Abela and M. Stampanoni, “Implementation of a fast method for high resolution phase contrast tomography,” Opt. Express, 14 8103 –8110 (2006). OPEXFF 1094-4087 Google Scholar


A. Groso et al., “Phase contrast tomography: an alternative approach,” Appl. Phys. Lett., 88 214104 (2006). APPLAB 0003-6951 Google Scholar


Y. D. Witte et al., “Bronnikov-aided correction for x-ray computed tomography,” J. Opt. Soc. Am. A, 26 890 –894 (2009). JOAOD6 0740-3232 Google Scholar


B. Münch et al., “Stripe and ring artifact removal with combined wavelet-Fourier filtering,” Opt. Express, 17 (10), 8567 –8591 (2009). OPEXFF 1094-4087 Google Scholar


W. van Aarle et al., “The ASTRA toolbox: a platform for advanced algorithm development in electron tomography,” Ultramicroscopy, 157 35 –47 (2015). ULTRD6 0304-3991 Google Scholar


W. J. Palenstijn, K. J. Batenburg and J. Sijbers, “The ASTRA tomography toolbox,” in 13th Int. Conf. Comput. and Math. Methods in Sci. and Eng. CMMSE, (2013). Google Scholar


“CXRO—the center for x-ray optics,” (2019) Google Scholar


G. Harauz and M. van Heel, “Exact filters for general geometry three dimensional reconstruction,” Optik, 73 146 –156 (1986). Google Scholar


M. van Heel, “Similarity measures between images,” Ultramicroscopy, 21 (1), 95 –100 (1987). ULTRD6 0304-3991 Google Scholar


M. van Heel and M. Schatz, “Fourier shell correlation threshold criteria,” J. Struct. Biol., 151 (3), 250 –262 (2005). JSBIEM 1047-8477 Google Scholar


L. A. Feldkamp, L. C. Davis and J. W. Kress, “Practical cone-beam algorithm,” J. Opt. Soc. Am. A, 1 612 –619 (1984). JOAOD6 0740-3232 Google Scholar


S. Cartier et al., “Micron resolution of MÖNCH and GOTTHARD, small pitch charge integrating detectors with single photon sensitivity,” J. Instrum., 9 C05027 (2014). Google Scholar

Biographies of the authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Marina Eckermann, Mareike Töpperwien, Anna-Lena Robisch, Franziska van der Meer, Christine Stadelmann, and Tim Salditt "Phase-contrast x-ray tomography of neuronal tissue at laboratory sources with submicron resolution," Journal of Medical Imaging 7(1), 013502 (20 February 2020).
Received: 11 August 2019; Accepted: 21 January 2020; Published: 20 February 2020 Logo
Cited by 17 scholarly publications.




Image quality


Signal to noise ratio

Back to Top