# Towards routine 3D characterization of intact mesoscale samples by multi-scale and multimodal scanning X-ray tomography

Ruiqiao Guo, Andrea Somogyi, Dominique Bazin, Elise Bouderlique, Emmanuel Letavernier, Catherine Curie, Marie-Pierre Isaure, Kadda Medjoubi

<figure><img src="image_1.png" /></figure>

## To cite this version:

Ruiqiao Guo, Andrea Somogyi, Dominique Bazin, Elise Bouderlique, Emmanuel Letavernier, et al. Towards routine 3D characterization of intact mesoscale samples by multi-scale and multimodal scanning X-ray tomography. Scientific Reports, 2022, 12 (1), pp.16924. 10.1038/s41598-022-21368-0. hal-03842035

HAL Id: hal-03842035

https://univ-pau.hal.science/hal-03842035

Submitted on 7 Nov 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Towards routine 3D characterization of intact mesoscale samples by multi-scale and multimodal scanning X-ray tomography

Ruiqiao Guo<fnref n="1" />,², Andrea Somogyi<fnref n="1" />✉, Dominique Bazin³, Elise Bouderlique⁴,⁵, Emmanuel Letavernier⁴,⁵,⁶, Catherine Curie⁷, Marie-Pierre Isaure⁸ & Kadda Medjoubi<fnref n="1" />✉

Non-invasive multi-scale and multimodal 3D characterization of heterogeneous or hierarchically structured intact mesoscale samples is of paramount importance in tackling challenging scientific problems. Scanning hard X-ray tomography techniques providing simultaneous complementary 3D information are ideally suited to such studies. However, the implementation of a robust on-site workflow remains the bottleneck for the widespread application of these powerful multimodal tomography methods. In this paper, we describe the development and implementation of such a robust, holistic workflow, including semi-automatic data reconstruction. Due to its flexibility, our approach is especially well suited for on-the-fly tuning of the experiments to study features of interest progressively at different length scales. To demonstrate the performance of the method, we studied, across multiple length scales, the elemental abundances and morphology of two complex biological systems, Arabidopsis plant seeds and mouse renal papilla samples. The proposed approach opens the way towards routine multimodal 3D characterization of intact samples by providing relevant information from pertinent sample regions in a wide range of scientific fields such as biology, geology, and material sciences.

Complex scientific problems in biology, earth-, environmental, and material sciences are inherently multi-scale. This requires the investigation of nanoscale features and functionalities within system-representative mesoscale samples to link those to emergent properties and functionalities at larger scales. This triggers an ever-increasing demand for new analytical tools capable of providing spatially resolved multi-scale information on intact, highly heterogeneous, or hierarchically structured samples in situ or *in operando*. Scanning hard X-ray imaging and tomography techniques are ideally suited to tackle this challenge due to their large penetration depth and inherently multimodal nature, where complementary information on the elemental distribution, morphology, crystal-line structure, and chemical speciation can be obtained simultaneously. Moreover, these non-invasive scanning techniques provide straightforward access to multiple-length scale experiments. Recent developments in fast continuous scanning, data acquisition<fnref n="1" />⁻⁴, and the high flux obtained at modern synchrotron-based hard X-ray nanoprobes have opened routine access to scanning 2D multimodal imaging. Amongst the scanning techniques, the high analytical sensitivity of X-ray Fluorescence (XRF) imaging provides unique possibilities in several scientific fields to study the role and fate of trace elements⁵⁻¹³. However, the unambiguous interpretation of 2D elemental distribution maps is not always straightforward or feasible and is especially problematic in the case of thick, complex samples. As such, the scientific community is highly demanding the extension of XRF imaging and other complementary scanning techniques to 3D tomography. Scanning XRF and multimodal tomography, where measuring projection images at different projection angles permits the reconstruction of the internal

<footnote n="1">¹Synchrotron SOLEIL, 91190 Saint-Aubin, France. ²Université Paris-Saclay, 91190 Gif-sur-Yvette, France. ³Institut de Chimie Physique, Université Paris-Saclay, CNRS, 91405 Orsay, France. ⁴UMR S 1155, Sorbonne Université, 75020 Paris, France. ⁵Institut National de la Santé et de la Recherche Médicale (INSERM), Sorbonne Université, 75013 Paris, France. ⁶Physiology Unit, Hôpital Tenon, AP-HP, 75020 Paris, France. ⁷IPSiM, Univ Montpellier, CNRS, INRAE, Institut Agro, 34000 Montpellier, France. ⁸IPREM, Université de Pau Et Pays de l'Adour, E2S UPPA, CNRS, 64053 Pau, France. ✉email: andrea.somogyi@synchrotron-soleil.fr; kadda.medjoubi@synchrotron-soleil.fr</footnote>
features by adapted reconstruction methods, provides unambiguous information about the internal elemental distributions in the context of sample morphology, crystalline structure, chemical speciation, etc.<sup>2,14–21</sup>. Up to now, the lengthy acquisition time necessary for these experiments<sup>2</sup> at 3rd generation synchrotrons has been one of the practical difficulties to its comprehensive utilization. Indeed, even if emerging sparse tomography techniques permit boosting the measurement throughput<sup>15,22,23</sup> by compromising spatial resolution, overcoming the time constraint remains a challenge at 3rd generation synchrotrons.

Meanwhile, at the dawn of 4th generation synchrotron sources, routine 3D scanning X-ray tomography is becoming within reach. Indeed, the two orders of magnitude larger flux available at 4th generation hard X-ray nanoprobes boosts the speed of scanning tomography measurements proportionally, paving the way towards high-throughput scanning X-ray tomography of mesoscale samples<sup>24</sup>. Hence, the implementation of robust, user-friendly scanning tomography workflow is crucial for the routine application of these techniques, similarly as has been published recently for high-throughput electron tomography<sup>25</sup>. Such workflow imposes flexible and fast on-site data processing and tomographic reconstruction adapted to the different number of projections (sparse and high-resolution tomography), diverse data quality (e.g., high and low-level elemental abundances, missing wedge), different imaging modalities (e.g., XRF, X-ray absorption, X-ray diffraction), and the possibility of adapting the field of view and spatial resolution to the examined phenomenon by using multi-scale or local tomography. Some recent developments tackle part of these requirements by sample-type specific processing of multimodal tomography data sets<sup>15,18</sup>, also in a semi-automatic way in the case of similar processing requirements<sup>17,26</sup>. A sparse sampling approach followed by sophisticated data processing has also been reported<sup>15,27</sup>. However, according to our knowledge, a robust, holistic approach addressing all requirements of flexible multi-scale and multimodal scanning 2D/3D X-ray tomography does not exist yet.

In the present paper, we introduce such a robust workflow for scanning multi-length scale XRF-tomography and complementary modalities. The presented workflow has been developed and implemented at the Nanoskopium beamline<sup>21</sup> of SOLEIL Synchrotron and includes semi-automatic data reconstruction. The proposed reconstruction algorithm yields good reconstruction data quality for diverse scientific fields with no need for parameter readjustment depending on the sample type. As the first step of this approach, sparse tomography provides a 3D overview of the entire meso- or microscale sample. The reconstructed sparse tomograms, containing relevant information to the investigated scientific problem, are used to choose pertinent regions for high spatial resolution single slice tomography, projection imaging, and local tomography. This approach permits optimizing the scanning tomography experiments and obtaining relevant information from pertinent sample regions in 2D or 3D during a user project. Thanks to recent developments<sup>3,22,28,29</sup>, this method paves the way towards statistically significant 3D studies, similar to those already available in full-field X-ray tomography<sup>30</sup>, and electron tomography<sup>25</sup>.

The performance of the workflow is demonstrated through the study of two highly heterogeneous mesoscale samples. A 700 µm thick wild type *Arabidopsis thaliana* seed, widely used as a model organism for plant biology studies, has been measured by sparse XRF tomography followed by high-resolution 2D tomography of some virtual slices and 2D projection imaging. The multi-length scale study of mesoscale renal papilla samples is crucial to investigating pathological calcification<sup>31</sup>. This study highly profited from high-resolution local XRF tomography of a micron-sized calcification sphere chosen from the reconstructed sparse tomograms.

## Results

**Multi-length scale scanning X-ray imaging/tomography workflow.** The workflow implemented at Nanoscopium for multi-length scale and multimodal scanning X-ray imaging and tomography experiments is presented in Fig. 1: the sample mounting and alignment are followed by sparse tomography of the whole sample, and then the visualization of the volume rendering of the reconstructed tomograms. The strategy of high-resolution (HR) measurements is based on these medium-resolution results. This workflow permits even users who are new to X-ray imaging and tomography, to accomplish all measurement and data reconstruction steps during their experiment. At Nanoscopium all user projects apply for XRF multi-scale imaging or tomography, which is often complemented by other modalities (absorption- or phase-contrast imaging/tomography, XANES, or XRD) to best tackle the actual scientific question. As such, in this paper, we have chosen scanning XRF tomography to demonstrate our approach.

As a first step, a sparse XRF/multimodal tomography measurement with 20 projection angles (see details in "Data acquisition" section) is followed by semi-automatic on-site tomography reconstruction. Direct visualization of a specimen's reconstructed medium-resolution 3D elemental distribution and morphology enables immediate identification of information pertinent to the research project. This permits users to choose the strategy for succeeding in high-resolution (HR) measurements for studying the smallest sample features crucial to tackling the scientific problem<sup>32</sup>. This HR measurement can be 2D projection imaging, 2D single slice tomography, local scanning 3D tomography, or any combination of these. Notably, a reasonable trade-off must be made during a tomography experiment regarding acquisition time, spatial resolution, and the number of samples wished to be investigated. After the high-resolution experiment, an estimate of the spatial resolution achievable by different projection numbers can be obtained by the Fourier ring correlation (FRC) calculation (detailed in the "Methods" section). If the resolution of the first sparse tomography results does not meet the user's requirements, then a second sparse tomography measurement can be performed. This will start with an angular offset equal to the half angular step of the first sparse tomography experiment. The second optional sparse tomography, having the same projection numbers as the first, results in doubling the total number of projections and hence improving the spatial resolution. As demonstrated in Fig. 1, all these modalities can be reconstructed semi-automatically on-site during the experiment using predefined and pre-parameterized reconstruction algorithms with no need for interaction from the users. Hence, users can guide their experiments and redefine the scientific objectives on-the-fly depending on the results obtained on-site.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Semi-automatic multi-scale XRF and multimodal scanning imaging and tomography workflow implemented at Nanoscopium. After sample mounting and alignment, sparse scanning tomography is performed on the whole sample. Some minutes after the experiment, medium resolution volume rendering is available for data interpretation by the proposed semi-automatic on-site data reconstruction. This permits optimizing the strategy of successive high-resolution (HR) measurements. The on-site reconstruction algorithm provides HR tomograms for data interpretation and further data processing some minutes after the HR experiments. The measurement and data reconstruction steps are included in orange and grey rectangles, respectively.</figcaption></figure>

Moreover, the HDF5 data format of the reconstructed tomograms is fully compatible with widely used cross-platform freewares and open-source data analysis tools, such as PyMCA (XRF data processing), ImageJ, Chimera (imaging/tomography data processing), etc. As such, further on-site or post-experiment processing, such as extraction of mean XRF spectra, quantification, combined treatment of different imaging modalities, and multivariate statistical analysis (Principal Component Analysis, Cluster Analysis, etc.)³³, can be obtained straightforwardly in any 2D/3D feature identified from the projection images or the reconstructed 2D/3D tomograms.

**Proof of principle test measurements: Wild-type *Arabidopsis thaliana* seeds.** *Arabidopsis thaliana* is a weed native to Eurasia and Africa with a short life cycle (≥ 6 weeks). This annual plant is a popular model organism in plant biology due to the knowledge of its genome and the availability of numerous mutants. The wild type (Col-0) is commonly used in plant biological laboratory experiments and even in Space experiments to study genetics, evolution, and development of flowering plants. As a model plant, *Arabidopsis thaliana* is also a powerful tool to investigate metal homeostasis and nutrient distribution, relevant questions about the world's food production and the agriculture industry. In this context, essential elements Fe and Mn in the plant seed are crucial for plant germination. Mn, involved in the photosystem II in chloroplasts, is needed for the vigour of germinating plants, and Fe is involved in various metabolic processes (respiration, photosynthesis...). These two elements have different physiological functions and different transporter pathways. Thus, differences are expected in their distribution and concentration within the seed. 2D XRF imaging is a powerful tool for locating trace elements in plants. However, assigning the measured metal distribution unambiguously to the grain's ultrastructure can be challenging due to the several hundreds of microns information depth of Mn and Fe. Moreover, sample sectioning, required to study the internal elemental distribution by 2D scanning XRF imaging, can be intricate for small and hard samples such as *A. thaliana* seeds and can induce tissue alteration and artifacts.

XRF and multimodal tomography are useful complementary tools for determining metal concentrations and distribution in intact seeds with minimal sample preparation¹⁶,³⁴,³⁵.

**Sparse scanning X-ray tomography of a whole seed.** We tested the performance of our tomography approach on wild-type *Arabidopsis thaliana* seeds. One advantage of using seeds, in general, is their low water content, resulting in lower amounts of radiation damage than more hydrated tissues. In order to obtain an overview of the elemental distribution of a whole mesoscale seed of 700 µm dimensions, we performed simultaneous sparse XRF (see Fig. 2) and absorption (Fig. S1) micro-tomography. The figure of merit of measuring a limited number of projections is to identify, within acceptable measurement times, the features of interest, in our case the distribution of Mn and Fe within different seed compartments. Since the total measurement time of a tomography experiment scales proportionally with the number of projections, it can be drastically reduced by measuring only a few angular projections. This also reduces the eventual radiation damage of the sample, which might be important in the case of scanning imaging of larger samples. However, the reduced number of projections should be chosen without significantly compromising the quality of the obtained 3D mesoscale tomograms. The reconstructed sparse 3D tomograms of Mn and Fe can be seen in Fig. 2A, B. It comprises 223 virtual slices
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Reconstructed sparse 3D Fe and Mn tomograms and a virtual slice showing the internal Fe and Mn distributions. (A,B) volume rendering of the Fe and Mn tomograms reconstructed by MLEM from 20 (A) and 5 angular projections (B). The white dashed line indicates the altitude of the virtual slice shown in (C). The red line marks the altitude of the high-resolution single slice tomography shown in Fig. 3. (C) Comparison of the reconstructed results obtained by the FBP and MLEM algorithms in the function of the measured number of angular projections. The columns correspond to two different algorithms: filtered back-projection (FBP) and maximum likelihood expectation maximization (MLEM). The lines correspond to 20, 10, and 5 angular projections, respectively. (D) Comparison of the intensity profiles of Fe obtained by FBP and MLEM. The Fe intensity profiles were extracted along the white straight line indicated in (C).</figcaption></figure>

and is calculated from 20 measured angular projections by the maximum-likelihood expectation-maximization (MLEM) algorithm (described in “Methods” section). It is clear from the reconstructed sparse tomograms that Fe is preferentially located in the provascular systems of the seed, and that Mn is mainly distributed in the abaxial area of both cotyledons, at the subepidermal level, as well as in the cortex area of the hypocotyl. These observations are in agreement with the results published in the literature³⁴,³⁵.

To optimize the measurement conditions and the corresponding reconstruction parameters of sparse tomography, we compared two different reconstruction algorithms, Filtered back-projection (FBP) and MLEM, in the function of the number of measured angular projections (5, 10, and 20). The comparison is demonstrated in Fig. 2C, D using the virtual slice marked by a white dashed line in Fig. 2A.

FBP algorithm is the standard solution for 3D tomography reconstruction due to its fast reconstruction time and easy implementation. However, because of the intrinsic nature of the algorithm, FBP has severe limitations
<figure><img src="image_6.png" /><figcaption>**Figure 3.** High-resolution 2D scanning XRF imaging and tomography of the Arabidopsis seed. (A) High-resolution 2D projection image of the Mn distribution at an appropriate angle chosen from the 3D sparse tomograms. (B,C) Reconstructed HR single slice of the Mn tomogram obtained by the FBP algorithm (B) and by the MLEM algorithm (C) measured at the altitude marked by the red line in Fig. 2A. (D) Mn sparse tomography reconstruction at the same slice obtained by MLEM.</figcaption></figure>

in the case of noisy datasets and highly under-determined measurement conditions with a small number of projections³⁶. This can be clearly seen in the first column in Fig. 2C, where the signal-to-noise ratio of all sparse tomograms obtained from 5, 10, and 20 projections is very low. Indeed, the strong streak artifacts of FBP, caused by the small number of projections, shadow the elemental distributions: the Mn and Fe distributions are hardly visible even in the case of 20 angular projections. In effect, as shown by the Fe intensity profiles (shown in Fig. 2D by the blue curves), these strong artifacts, resulting also in negative intensity values, are hiding the information on the Fe variation within the sample. As such, since FBP cannot significantly reduce streak artifacts at low number of projections, it is not adapted to sparse tomography reconstruction.

MLEM algorithm (second column of Fig. 2C) is superior to FBP in handling noisy datasets. Moreover, it includes the non-negativity constraint assumption. In our workflow, we included a simple automatic stopping criterion for the MLEM algorithm with the smallest possible noise as a figure of merit (as detailed in “Meth-ods” section). This choice results in terminating the reconstruction process at small number of iterations (thus with short reconstruction time) introducing only a few, weak artifacts in the reconstructed tomograms. The reconstructed tomograms in Fig. 2C illustrate the robustness of this algorithm in treating limited number of projections. Moreover, with increasing number of projections the reconstruction artifacts of MLEM are becom-ing weaker. This is illustrated by the purple Fe profiles in Fig. 2D, where the strongest reconstruction artifacts marked by the black circles in the 5-projection tomogram (last line in D) diminished in the 20-projection tomogram (first line of D).

**High-resolution single slice tomography.** The above-described 3D sparse tomography permits choos-ing the best angular position for high-resolution 2D projection imaging, as shown in Fig. 3A for Mn distribution. It also allows choosing the altitude(s) for high-resolution single slice tomography (Fig. 3B, C). This provides insight into the internal Mn and Fe distributions of subcellular features within intact seeds. We performed HR single slice XRF tomography at the altitude shown by the red line in Fig. 2A. The FBP and MLEM reconstruc-tions of the internal 2D Mn distribution can be seen in Fig. 3B, C. We included for comparison the sparse tomogram of Mn obtained by MLEM in the very same virtual slice (Fig. 3D). The reconstructed features are in
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Spatial resolution estimation by Fourier ring correlation (FRC) method. Estimation of spatial resolution for tomograms reconstructed by MLEM from different numbers of projections chosen from the measured 360 angular projection dataset. The fixed 0.5 threshold was used for FRC analysis. The Nyquist frequency is 0.25 µm⁻¹.</figcaption></figure>

<table><caption>**Table 1.** Estimation of the spatial resolution of sparse tomograms (reconstructed by the MLEM algorithm) by FRC analysis in the function of the measured number of projections. The spatial resolution defined by the Nyquist sampling condition is included in the 3rd column for comparison.</caption><thead><tr><td>Number of projections</td><td>FRC analysis: Resolution (µm)</td><td>Nyquist sampling condition: resolution (µm)</td></tr></thead><tbody><tr><td>5</td><td>20.5</td><td>184.7</td></tr><tr><td>10</td><td>12.9</td><td>92.4</td></tr><tr><td>20</td><td>8.4</td><td>46.2</td></tr><tr><td>40</td><td>5.6</td><td>23.1</td></tr></tbody></table>

good agreement between the sparse (Fig. 3D) and high-resolution tomography (Fig. 3B, C) results. However, as expected, high-resolution tomography reveals fine details with improved spatial resolution, which are non-detectable or hardly identifiable by sparse tomography. For such HR tomography, the use of FBP (Fig. 3B) is straightforward and is the fastest reconstruction algorithm, where the computation time is proportional to the registered number of projections. However, even with 360 angular positions, there are non-negligible streak artifacts due to insufficient angular projections compared to the number of scanned pixels in this mesoscale sample (see Nyquist sampling condition in Eq. 3). The reconstruction result of the MLEM algorithm (Fig. 3C) provides better contrast. Moreover, the computation time of < 2 s/sinogram is significantly faster than the measurement time of some hours (and will be comparable with the measurement time of some minutes/sparse tomography at a 4th generation synchrotron). As such, the MLEM algorithm proved to be the best compromise for semi-automatic image reconstruction for both 3D sparse tomography and 2D high-resolution single slice tomography experiments.

**Spatial resolution of the tomograms reconstructed by MLEM.** To estimate the spatial resolution of the tomograms reconstructed by the MLEM algorithm, we extracted two subsets with uniform angular sampling from the 2D high-resolution dataset. Taking the 2D high-resolution sinogram of 360 projections as an example, the full projection dataset was divided into two subsets, both having the same projection numbers. These two subsets were reconstructed independently by MLEM with the automatically determined iteration numbers. The two independent tomograms were then used for FRC analysis.

Figure 4 shows the FRC curves between the two tomograms reconstructed using a different number of projections from the same 2D high-resolution sinogram. The spatial resolution was determined at the intersection of the FRC curve and the fixed 0.5 threshold³⁷. Table 1 shows the spatial resolution determined by the 0.5 threshold criterion in the function of the number of projections. The spatial resolution improves with the number of projections as expected (see Figs. S1, C and E). A resolution estimate of 8.4 µm was obtained for the 2D projection tomogram. The MLEM iterative method has resulted in a significant resolution improvement compared to the spatial resolution given by the Nyquist sampling condition.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Sparse and local scanning X-ray tomography of a renal papilla sample. (A) Optical microscope image of the mouse papilla sample, mimicking human kidney calcifications. The scanned region is marked in red. (B) Transmission image of the tip of the sample. The small red rectangle in the middle marks the sample region chosen for local tomography. The zoom-in image of this region is inserted into the upper-right corner. (C) Sparse tomography reconstruction of Zn and Ca. The small red rectangle in the middle marks the sample region chosen for local tomography. The zoom-in image of this region is inserted into the upper-right corner. (D) Reconstructed 3D local tomogram of Ca and Zn of the calcification micro-sphere marked by the red rectangles in (B,C). (E) Cut-off view at the middle of the calcification micro-sphere presented in (D).</figcaption></figure>

**Renal papilla sample.** In order to test the performance of our workflow in missing wedge measurement conditions and for HR 3D local tomography, we investigated mouse renal papilla samples. Nowadays, the formation of renal stones affects 10% of the global population³⁸. Most of the stones develop on Randall's Plaque, a mineral deposit at the tip of renal papillae³⁹. The major components of Randall's Plaque are calcium phosphate apatite and amorphous carbonated calcium phosphate⁴⁰. Therefore, the distribution of Ca indicates the position of calcification in the renal papilla sample. Studying the correlation and colocalization between Ca and other major and trace components permits to reveal those trace elements, that are involved in early-stage calcification processes. Understanding their role in the pathological process opens the way toward efficient prevention and treatment of the renal stone formation. The sample we used in the present study is a mouse papilla affected by calcium phosphate deposits similar to the human Randall's plaque (see⁴¹ for more details on the sample preparation).

**Sparse and high-resolution XRF tomography of renal papilla.** The semi-automatic workflow has been used to study the elemental abundancies within 50 µm thick slices of renal papilla samples (Fig. 5A, B). In this case, the fixation of the thin biological sample of ~500 µm lateral dimensions on a Si₃N₄ membrane caused a missing wedge of 2×28° in the tomography measurement. This poses a specific challenge to tomographic reconstruction⁴². This was successfully tackled by the MLEM algorithm, which proved to be the best compromise also for missing wedge tomography. 3D sparse tomography performed on the top of the sample showed substantial colocalizations between Zn and Ca (Fig. 5C) within dense sample regions revealed by simultaneous transmission imaging (Fig. 5B). However, the ~8.5 µm medium spatial resolution of sparse tomography does not permit to reveal the details of micrometer-sized Ca- and Zn-rich features, that are related to early-stage calcification process. Gaining insight, with high spatial resolution, into the elemental distribution of these micrometer-sized Ca and Zn containing spherical structures is crucial. We performed 3D local tomography around a chosen micro-sphere (red rectangles in Fig. 5B, C). The angular projection images were measured with 500 nm pixel size. Figure 5D, E show the volume rendering of Zn and Ca within the measured micrometer-scale calcification sphere. In the cut-off view of the 3D volume rendering in Fig. 5E the internal distribution of Zn and Ca is shown within the ~10 µm dimension calcification sphere. These results reveal that Zn is enriched within a few micrometers thick rim on the surface of the calcification micro-sphere. Since Zn is considered to be a marker of inflammation⁴³, this result also indicates that the calcium phosphate deposition in the medullary interstitial is a pathological process.

**Mean XRF voxel-spectrum.** Next to the elemental distribution maps, complete local XRF spectral information is necessary to obtain detailed information on the chemical composition of local features. This might reveal rare characteristics appearing only in particular sample locations/voxels ("needle in the haystack problem").
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Extraction of the mean XRF spectra of a chosen area in a virtual slice of the Arabidopsis seed. (A) Reconstructed high-resolution virtual slice of the sum XRF spectra of the Arabidopsis seed. (B) Comparison of the mean XRF spectra of the selected area, marked by a red dashed circle in (A), obtained by FBP (red curve) and MLEM (blue curve) algorithm. (C) The mean XRF spectra of the selected area obtained by MLEM from 360 projections (blue curve) and from 20 projections (fuchsine curve).</figcaption></figure>

lem"). Since we collect the full XRF spectrum in each measured pixel during data collection, next to the tomograms of the predefined elements, it is also possible to reconstruct the full XRF spectrum in each voxel of the tomogram. In other words, we can include a spectral dimension to the reconstructed dataset, which can be used to extract the mean XRF spectrum of any chosen sample area or volume. Such complete hyperspectral tomography reconstruction was realized for each energy channel of the measured XRF spectra. At 10 keV excitation energy, this results in 1000 sinograms for each virtual slice, and the reconstruction process, even with the MLEM iterative method and multi-core, is becoming quite time-consuming. For example, the reconstruction of the high-resolution 3D hyperspectral XRF tomogram (with the energy as the 3rd dimension) took ~ 1 h by MLEM for the Arabidopsis seed.

Figure 6A shows the virtual slice of the Arabidopsis seed reconstructed from the sum-XRF spectra of the dataset containing 360 projections. The mean XRF spectra extracted from the area marked by a red circle in Fig. 6A by FBP and MLEM can be seen in Fig. 6B, C. The mean XRF spectra of the chosen region obtained by FBP and MLEM do not show a significant difference; the total intensities of the XRF spectra agree within 5%. As such, in the case of high-resolution tomography, FBP is a good compromise to obtain local spectral information within a limited reconstruction time.

As the next step, we selected equiangularly 20 projections out of the 360 ones to simulate a sparse tomography in identical experimental conditions. This allows comparing the mean XRF spectra obtained by sparse and high-resolution tomographies. In the case of sparse tomography, FBP cannot be used for data reconstruction due to the low number of projections. As such, we reconstructed the sparse tomogram for each energy channel by MLEM. The mean XRF spectrum of the same region as before (red circle in Fig. 6A) was calculated from this sparse tomography dataset. Figure 6C shows that the two mean XRF spectra agree well, which illustrates that the MLEM algorithm provides reliable hyperspectral tomograms and mean XRF spectra even for a severely limited number of projections.

The mean XRF spectrum within a selected volume of interest of the calcification sphere (shown in Fig. 5D) of the renal papilla is demonstrated in Fig. 7. The reconstruction process for the 4D hyperspectral tomography dataset (with the energy as the 4th dimension) took around 1 h by the MLEM algorithm. To obtain the mean XRF spectra within a volume of interest, we selected the 3D volumes of interest by ImageJ, a freeware frequently used by our users. We applied this mask to each energy channel in the reconstructed 4D dataset. The extracted mean XRF spectra shown in Fig. 7B reveal that Ca is the dominant element in the core of the calcification microsphere (red XRF spectrum), which also contains Zn. Within the surface rim (indicated by green in Fig. 7A and green XRF spectrum in Fig. 7B), the increased Zn to Ca ratio, due to the threefold larger Zn and ~ 20% smaller Ca content, provides evidence of the association of Zn to the pathological process.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Extraction of the mean XRF spectra within the volume of interest of the reconstructed calcification micro-sphere shown in Fig. 5D,E. (A) Reconstructed result by MLEM algorithm from the total intensity of the XRF spectra, grey volume: reconstructed total sphere volume, red and green volumes: core and surface rim regions, respectively. (B) Comparison of the mean XRF spectra of the core (red) and surface rim (green) volumes.</figcaption></figure>

Discussion

The proposed holistic multi-scale and multimodal scanning X-ray tomography workflow, implemented at Nanoscopium, was tested on mesoscale samples. To tackle the scientific problems presented in the recent paper, multi-scale XRF- and scanning absorption contrast tomography were the best-adapted imaging modalities.

The best compromise for semi-automatic on-site tomography reconstruction for sparse-, local- and high-resolution scanning X-ray tomography is proved to be the MLEM algorithm if it is benefited from the apt early stopping strategy. This algorithm also efficiently handles the missing wedge sampling conditions. The FBP algorithm is, in general, not adapted for the reconstruction of sparse and missing wedge tomographies since the measured dataset does not fulfill the optimal Nyquist sampling criterion. In the case of high-resolution hyperspectral tomography, which results in large 3D/4D datasets, the figure of merit is a compromise between the reconstructed data quality and the reconstruction time. Here, FBP is the best compromise to obtain, within a limited reconstruction time, reliable mean XRF-spectra from regions/volumes chosen from the reconstructed XRF hyper-spectral tomograms. The semi-automatic use of these two algorithms according to the above-detailed conditions permits obtaining a flexible semi-automatic workflow providing good quality on-site reconstruction for various samples in diverse experimental conditions. Due to its flexibility, our approach is especially well suited for on-the-fly tuning of the experiments to study features of interest progressively at different length scales.

We must note that scanning XRF imaging is generally a semi-quantitative approach. However, in the present paper, we presented the 2D and 3D elemental intensity distributions (number of characteristic X-ray Fluorescence photons per given dwell time). As the next step, we intend to fully integrate the quantification method in the semi-automatic user-friendly data-processing pipeline and to provide it for all scanning XRF tomography experiments as an option. For this, the measurement conditions will be calibrated by an adequate reference sample, and simultaneous scanning X-ray absorption tomography will be used for self-absorption correction in each voxel of the tomogram for each reconstructed element.

The robust scanning tomography method implemented at the Nanoscopium beamline opens the way for non-expert users towards routine non-destructive multi-length-scale characterization of complex samples during a standard beam-time of 3–5 days. The reconstructed tomograms can be treated by widely used freewares such as ImageJ, Chimera, etc. This allows straightforward data handling for users during and after their experiments. The workflow provides the possibility to study relevant 3D micro-features of several mesoscale samples during a routine user experiment, even at 3rd generation synchrotrons. Next to XRF and scanning absorption tomography, we have also extensively tested the workflow for scanning phase contrast- and X-ray diffraction tomography studying perovskite samples and bio-mineralization. These results will be presented in a separate paper. Moreover, the application of our approach for XANES tomography, which can be considered as a variant of hyperspectral tomography⁴⁴, is straightforward.

Such a robust semi-automatic flexible scanning multimodal tomography workflow will be a scientific game-changer at emerging 4th generation synchrotron sources, where data throughput of scanning hard X-ray techniques is boosted by ~ two orders of magnitude. The new possibility to explore a multitude of sample characteristics simultaneously, with high analytical sensitivity, at hierarchical length-scales in 3D, in a statistically significant manner in meso- and micro-scale samples, will revolutionize a wide range of scientific fields in ways that we can currently only dream of and will provide a unique complement to already existing state-of-the-art multi-scale and fast full-field X-ray tomography techniques⁴⁵⁻⁴⁷. These high through-put scanning X-ray imaging techniques will also complement emerging state-of-the-art laboratory scanning charged-particle microscopy and tomography and other conventional laboratory microscopy techniques providing 2D/3D elemental and morphology information with nanometer resolution. For example, scanning electron microscopy⁴⁸ has limited analytical sensitivity (~ 0.1–1%) for the analysis of heavy elements, and its small depth of information (a few microns) limits its non-invasive (without sample sectioning) application to surface studies of large samples or the 3D study of
light, major, and minor components of tiny, micrometer-sized samples. Due to sample radiation damage caused by charged-particle microscopies, multiple measurement on the very same sample region is often problematic. Scanning hard XRF imaging and tomography has high analytical sensitivity (with trace, ≤ ppm detection limit) for transition metals and heavier elements. Moreover, their large information depth of several tens/hundreds of microns permits the non-invasive multilength-scale 3D study of mesoscale samples. Multiple measurements and hence, multi-scale and local tomography is readily available by scanning X-ray imaging techniques, where radiation damage is smaller than by charged particle microscopies.

## Methods

**Data acquisition.** The experiments were performed at the KB-based nanoprobe station of the Nanoscope beamline of Synchrotron Soleil²¹. Two energy-dispersive silicon drift detectors, placed at 120° to the incident beam path, collected the X-ray fluorescence (XRF) spectra. The intensity of the incident and the transmitted beam has been collected in each pixel by two Si diodes placed before and behind the sample, respectively. An RT100 air-bearing rotation stage has been used for tomography measurements. The experiments have been performed by the Flyscan architecture¹.

In the case of scanning X-ray single slice tomography, the number of angular projections has been chosen according to the Nyquist sampling condition (see Eq. 3) in order to preserve the spatial resolution along the horizontal direction in the reconstructed tomograms.

The medium resolution 3D elemental distribution of the Arabidopsis seed has been obtained by sparse XRF tomography at 11 keV excitation X-ray energy. We collected 20 projections at 18 degrees angular intervals over 360°. These scanning parameters are proposed by default for sparse tomography by our workflow since, according to our experience with multiple samples and imaging modalities, they are proved to be the best-balanced parameters for almost all experimental conditions. At each projection angle, a total area of 581 × 445 µm² has been scanned with 2 µm pixel size and 20 ms dwell time in continuous scanning mode. The collection of the full XRF sparse tomography dataset of this mesoscale sample took 8 h 40 min. For high-resolution single slice tomography, we measured 360 angular projections over 360° with a lateral step size of 2 µm and a dwell time of 20 ms/pixel, with a total measurement time of 35 min. For high-resolution 2D projection imaging, a field of view of 375 × 434 µm² was scanned with a step of 500 nm and a dwell time of 40 ms per pixel. The total acquisition time was 7 h 30 min.

In the case of the mouse renal papilla sample, the sparse tomography was performed at 12 keV. A 586 × 272 µm² region of the sample was mapped at 22 angular positions over 360° with a lateral step size of 2 µm and an exposure time of 20 ms/pixel, with a total measurement time of 5 h. This is followed by a high-resolution local tomography of a local sample volume with 65 angular positions over 360° with a dwell time of 40 ms/pixel and a uniform lateral step size of 500 nm, covering a total area of 67 × 16 µm². The total acquisition time was 3 h.

**Data processing.** The elements present in the investigated samples were identified from the sum-spectra of all angular projections. The sum XRF spectra were fitted by the PyMCA software⁴⁹ for the identification of the elements present in the sample. The elemental distribution maps and sinograms were then extracted from the raw dataset by integrating pre-selected spectral regions of interest corresponding to the detected elements. The transmitted and incident beam intensity ratio at each pixel gives access to the sample absorption maps and sinograms. The sinograms were then reconstructed either by the MLEM or FBP algorithms. The overall reduction process is performed automatically with a robust in-house MATLAB code. The reconstructed volumetric data were exported from MATLAB as 16-bit z-stacks and imported either to ImageJ⁵⁰ for analysis or to Chimera⁵¹ for 3D visualization.

The overall data processing is performed on a workstation with an Intel® Xeon® Processor E5-2630 v3 @ 2.40 GHz × 32 CPU with 125.8 GB of system RAM.

**Reconstruction algorithms.** Two different tomography reconstruction algorithms have been integrated into the workflow: the FBP⁴²,⁵²,⁵³ and MLEM⁵⁴ methods. The first, which is an analytical filtering inversion technique, is the most commonly used method in routine tomography reconstruction. It is the fastest method, exploiting fast Fourier transforms. However, FBP suffers from a lack of robustness when the measurements are sparse, of low contrast, or noisy⁵⁵,⁵⁶, which is often the case in XRF tomography. MLEM is an iterative algorithm that explicitly accounts for the noise affecting the data and imposes positivity on the estimated pixels. MLEM was initially developed for the analysis of positron emission tomography (PET) data. It assumes a Poisson distribution of the acquired photon statistics and thus might be more noise-tolerant than FBP²². MLEM belongs to the class of majorization-minimization algorithms⁵⁷. By construction, it decreases the negative log-likelihood of the estimated image of the given data monotonically. Moreover, it involves only simple multiplicative updates in such a way that the estimated image has positive pixel values, if all the pixels of the starting assumption image are positive. The MLEM algorithm provides good performance with sparse datasets, and the streak artifacts can be reduced as long as it is combined with an appropriate early stopping strategy, as is described hereafter.

**Determination of the number of iterations: automatic early stopping strategy for MLEM.** Although the loss of the MLEM algorithm decreases along with iterations, there is also an increase in the noise amplitude, which significantly impacts the final reconstruction results. This is typical for non-regularized reconstruction methods. Therefore a so-called early stopping strategy with a well-tuned criterion is crucial and must be integrated into the workflow⁵⁸. We implemented the figure of merit about the Normalized Root Mean Square Error Deviation (NRMSED) as an indicator:
$$NRMSED = \sqrt{\frac{\sum_{i=1}^{N} (x(i) - \tilde{x}(i))^2}{\sum_{i=1}^{N} x(i)^2}} \quad (1)$$

where $x(i)$ is the sinogram of the measured data, $\tilde{x}(i)$ is the estimated sinogram by MLEM and N is the pixel number in the sinogram. NRMSED indicates the accordance between the reconstructed and measured data. In addition to this, information about the derivative of NRMSED between two successive iterations should also be considered, and the ratio R is therefore defined as:

$$R = \frac{\Delta NRMSED}{NRMSED} \quad (2)$$

where $\Delta NRMSED = NRMSED^{(i)} - NRMSED^{(i-1)}$ is the derivative of NRMSED. This similarity increases monotonically with the number of iterations. A threshold is set to terminate the algorithm when the evolution curve of R is close to zero: $R \ge \text{threshold}$. We have chosen -0.15% as threshold, since this value, obtained as average over numerous test data sets, gave robust and rapid results for all tested data sets and imaging modalities.

In the implementation of our workflow, a relatively large number (~ e.g., 200) of iterations is preset to limit the maximum reconstruction time. The reconstruction of the tomograms using the MLEM algorithm is performed slice by slice, and the ratio R is calculated for each slice after each iteration. The algorithm terminates automatically for a given virtual slice after reaching the threshold and proceeds to reconstruct the next slice. Optionally, the reconstruction could perform a few additional iterations after reaching the R threshold. According to our experience, the algorithm generally reaches the $R \ge \text{threshold}$ condition around 20 iterations, well before 200 iteration steps.

The process of the optimization of the iteration number can be summarized as follows:

(1) For all the slices, set the threshold to -0.15% and the initial number of iterations to a relatively large value (~ 200);

(2) Start the reconstruction for the first slice using the MLEM algorithm while calculating the ratio R for each iteration;

(3) Stop the iteration when the condition $R \ge \text{threshold}$ is reached;

(4) Repeat (2)-(3) for all slices.

This robust automatic reconstruction process provides exploitable results even for 5 angular projections, and it assures that even inexperienced users in synchrotron-based imaging/tomography can reconstruct their tomographic dataset on-site.

**Spatial resolution evaluation.** Nyquist angular sampling condition. During the experiment, a reasonable trade-off must be reached between the acquisition time and spatial resolution. For the analytical reconstruction methods, the number of projections has to satisfy the Nyquist angular sampling condition⁵⁹ in order to preserve the spatial resolution along the horizontal direction in the reconstructed tomogram:

$$n_{proj} = \frac{\pi}{2} N \quad (3)$$

where $n_{proj}$ is the number of projections over 180°, and N is the number of pixels along the scanning direction. The corresponding maximum resolvable spatial frequency $f_{res}$ in the Fourier domain can be written as:

$$f_{res} = \frac{1}{\pi} \frac{n_{proj}}{PN} \quad (4)$$

where P is the pixel size.

The resolution limit R in direct space is then:

$$R = \frac{\pi}{2} \frac{PN}{n_{proj}} \quad (5)$$

This resolution limit decreases with increasing number of projections.

Fourier ring correlation for spatial resolution estimation. The most common method to estimate spatial resolution is the knife-edge technique. However, in the case of sparse tomography of highly heterogeneous samples, it can be challenging to find a sharp edge for the knife-edge method. To overcome this difficulty, the Fourier ring correlation (FRC)⁶⁰⁻⁶² method can be used as a general approach.

For the iterative reconstruction methods, the estimation of the achievable spatial resolution in the function of the number of projections can be described by the Fourier ring correlation approach. FRC estimates this by measuring the normalized cross-correlation of two independent datasets of the sample. The spatial frequency elements at different radii can be integrated circularly in the frequency domain:

$$FRC_{12}(r) = \frac{\sum_{r_i \in R} F_1(r_i) \cdot F_2(r_i)^*}{\sqrt{\sum_{r_i \in R} |F_1(r_i)|^2 \cdot \sum_{r_i \in R} |F_2(r_i)|^2}} \quad (6)$$

where $r_i$ is the ith frequency element at radius r, $F_1$ and $F_2$ denote the Fourier transform of the two reconstructed tomograms.

At a specified cut-off threshold, the FRC curve drops below the threshold, indicating an indistinguishable signal-to-noise ratio. This defines the spatial resolution.
Theoretically, the two data sets used in FRC calculation should be formed by independent measurement⁶³,
but this can be impractical, especially in the high-resolution scanning experiments, where the acquisition time
is generally several hours at 3rd generation synchrotrons. To overcome this problem, the two subsets used
for FRC calculation were extracted from the high-resolution tomography sinogram, from which independent
tomographic images of the same slice were reconstructed. This resolution estimation method is mostly applied
to 2D high-resolution datasets.

Ethics approval

*Arabidopsis thaliana* seeds: The research on the wild-type Col0 seed *Arabidopsis Thaliana* seed complies with rel-
evant institutional, national, and international guidelines and legislation. The wild-type *Arabidopsis thaliana* seeds
had been provided by C. Curie, and is the same as described by Carrió-Seguí, À., Romero, P., Curie, C. et al. "Cop-
per transporter COPT5 participates in the crosstalk between vacuolar copper and iron pools mobilization." In *Sci
Rep* 9, 4648 (2019) https://rdcu.be/cROAu. Renal papilla samples: All studies were performed in accordance with
the European Union, NIH guidelines (Comité d'Ethique en Experimentation Charles Darwin C2EA-05) and
all of our methods are reported as recommended by ARRIVE guidelines as previously described. The project
was authorized by the Health Ministry and local Ethics Committee (authorization 11420 2017092015335292).

Data availability

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary
Materials. Additional data related to this paper will be available from the authors upon reasonable request and in
accordance with the Journal policy described in the Instructions for Authors (Availability of materials and data).

Received: 5 July 2022; Accepted: 26 September 2022

Published online: 08 October 2022

References

1. Medjoubi, K. et al. Development of fast, simultaneous and multi-technique scanning hard X-ray microscopy at Synchrotron Soleil. *J. Synchrotron. Rad.* **20**, 293–299 (2013).

2. Deng, J. et al. Correlative 3D x-ray fluorescence and ptychographic tomography of frozen-hydrated green algae. *Sci. Adv.* **4**, eaau4548 (2018).

3. Kopittke, P. M. et al. Synchrotron-based X-ray fluorescence microscopy as a technique for imaging of elements in plants. *Plant Physiol.* **178**, 507–523 (2018).

4. Wang, Y. et al. A high-throughput x-ray microtomography system at the Advanced Photon Source. *Rev. Sci. Instrum.* **72**, 2062–2068 (2001).

5. Figon, F. et al. Catabolism of lysosome-related organelles in color-changing spiders supports intracellular turnover of pigments. *Proc. Natl. Acad. Sci. USA* **118**, e2103020118 (2021).

6. Wolf, M., May, B. M. & Cabana, J. Visualization of electrochemical reactions in battery materials with X-ray microscopy and map-
ping. *Chem. Mater.* **29**, 3347–3362 (2017).

7. Dean, C., Le Cabec, A., Spiers, K., Zhang, Y. & Garrevoet, J. Incremental distribution of strontium and zinc in great ape and fossil
hominin cementum using synchrotron X-ray fluorescence mapping. *J. R. Soc. Interface* **15**, 20170626 (2018).

8. Xie, R. et al. Seasonal zinc storage and a strategy for its use in buds of fruit trees. *Plant Physiol.* **183**, 1200–1212 (2020).

9. Yan, H. et al. Multimodality hard-x-ray imaging of a chromosome with nanoscale spatial resolution. *Sci. Rep.* **6**, 20112 (2016).

10. Kosior, E. et al. Combined use of hard X-ray phase contrast imaging and X-ray fluorescence microscopy for sub-cellular metal quantification. *J. Struct. Biol.* **177**, 239–247 (2012).

11. Sforna, M. C. et al. Evidence for arsenic metabolism and cycling by microorganisms 2.7 billion years ago. *Nature Geosci.* **7**, 811–815 (2014).

12. Chaurand, P. et al. Multi-scale X-ray computed tomography to detect and localize metal-based nanomaterials in lung tissues of
in vivo exposed mice. *Sci. Rep.* **8**, 4408 (2018).

13. Serpell, C. J. et al. Carbon nanotubes allow capture of krypton, barium and lead for multichannel biological X-ray fluorescence imaging. *Nat. Commun.* **7**, 13118 (2016).

14. Helliwell, J. R. et al. Applications of X-ray computed tomography for examining biophysical interactions and structural develop-
ment in soil systems: A review: X-ray computed tomography for soil physical properties. *Eur. J. Soil Sci.* **64**, 279–297 (2013).

15. Gao, X., Yang, Y., Yang, S., Ma, Y. & Chen, M. Microstructure evolution of chalcopyrite agglomerates during leaching—A synchro-
tron-based X-ray CT approach combined with a data-constrained modelling (DCM). *Hydrometallurgy* **201**, 105586 (2021).

16. Punshon, T. et al. The Role of CAX1 and CAX3 in elemental distribution and abundance in Arabidopsis seed. *Plant Physiol.* **158**, 352–362 (2012).

17. Parsons, A. D. et al. Automatic processing of multimodal tomography datasets. *J. Synchrotron. Rad.* **24**, 248–256 (2017).

18. Suuronen, J.-P. & Sayab, M. 3D nanopetrography and chemical imaging of datable zircons by synchrotron multimodal X-ray tomography. *Sci. Rep.* **8**, 4747 (2018).

19. Hong, Y. P., Chen, S. & Jacobsen, C. A new workflow for X-ray fluorescence tomography: MAPStoTomoPy. In: (ed. Lai, B.) 95920W (2015). https://doi.org/10.1117/12.2194162

20. Gürsoy, D., Biçer, T., Lanzirotti, A., Newville, M. G. & De Carlo, F. Hyperspectral image reconstruction for X-ray fluorescence tomography. *Opt. Express* **23**, 9014 (2015).

21. Somogyi, A. et al. Optical design and multi-length-scale scanning spectro-microscopy possibilities at the Nanoscopium beamline of synchrotron soleil. *J. Synchrotron. Radiat.* **22**, 1118–1129 (2015).

22. Bourassa, D., Gleber, S.-C., Vogt, S., Shin, C. H. & Fahrni, C. J. MicroXRF tomographic visualization of zinc and iron in the zebrafish embryo at the onset of the hatching period. *Metallomics* **8**, 1122–1130 (2016).

23. Villarraga-Gómez, H. & Smith, S. T. Effect of the number of projections on dimensional measurements with X-ray computed tomography. *Precis. Eng.* **66**, 445–456 (2020).

24. Li, P. et al. 4th generation synchrotron source boosts crystalline imaging at the nanoscale. *Light Sci. Appl.* **11**, 73 (2022).

25. Schwartz, J. et al. Real-time 3D analysis during electron tomography using tomviz. *Nat. Commun.* **13**, 4458 (2022).

26. Atwood, R. C., Bodey, A. J., Price, S. W. T., Basham, M. & Drakopoulos, M. A high-throughput system for high-quality tomographic reconstruction of large datasets at Diamond Light Source. *Phil. Trans. R. Soc. A.* **373**, 20140398 (2015).

27. Liu, Z. et al. TomoGAN: Low-dose synchrotron x-ray tomography with generative adversarial networks: Discussion. *J. Opt. Soc. Am. A* **37**, 422 (2020).
28. Somogyi, A. & Mocuta, C. Possibilities and challenges of scanning hard X-ray spectro-microscopy techniques in material sciences. *AIMS Mater. Sci.* **2**, 122–162 (2015).

29. Chu, Y. S. et al. Multimodal, multidimensional, and multiscale X-ray imaging at the national synchrotron light source II. *Synchrotron Radiat. News* **33**, 29–36 (2020).

30. Odstrcil, M. et al. Ab initio nonrigid X-ray nanotomography. *Nat. Commun.* **10**, 2600 (2019).

31. Bazin, D., Daudon, M., Combes, C. & Rey, C. Characterization and some physicochemical aspects of pathological microcalcifications. *Chem. Rev.* **112**, 5092–5120 (2012).

32. de Jonge, M. D. & Vogt, S. Hard X-ray fluorescence tomography—An emerging tool for structural visualization. *Curr. Opin. Struct. Biol.* **20**, 606–614 (2010).

33. Johnson, R. A. & Wichern, D. W. *Applied Multivariate Statistical Analysis* (Pearson Prentice Hall, 2007).

34. Kim, S. A. et al. Localization of iron in *Arabidopsis* seed requires the vacuolar membrane transporter VIT1. *Science* **314**, 1295–1298 (2006).

35. Eroglu, S. et al. Metal tolerance protein 8 mediates manganese homeostasis and iron reallocation during seed development and germination. *Plant Physiol.* **174**, 1633–1647 (2017).

36. Schofield, R. et al. Image reconstruction: Part 1—Understanding filtered back projection, noise and image acquisition. *J. Cardiovasc. Comput. Tomogr.* **14**, 219–225 (2020).

37. Sorzano, C. O. S. et al. A review of resolution measures and related aspects in 3D electron microscopy. *Prog. Biophys. Mol. Biol.* **124**, 1–30 (2017).

38. Trinchieri, A. Epidemiology of urolithiasis: An update. *Clin. Cases Miner. Bone Metab.* **5**, 101–106 (2008).

39. Randall, A. An hypothesis for the origin of renal calculus. *N. Engl. J. Med.* **214**, 234–242 (1936).

40. Van de Perre, E. et al. Randall's plaque as the origin of idiopathic calcium oxalate stone formation: an update. *C. R. Chim.* **24**, 1–19 (2022).

41. Letvernier, E. et al. ABCC6 deficiency promotes development of Randall plaque. *JASN* **29**, 2337–2347 (2018).

42. Kak, A. C. & Slaney, M. *Principles of Computerized Tomographic Imaging* (Society for Industrial and Applied Mathematics, 2001).

43. Carpentier, X. et al. High Zn content of Randall's plaque: A μ-X-ray fluorescence investigation. *J. Trace Elem. Med Biol.* **25**, 160–165 (2011).

44. Dong, X. et al. A review of hyperspectral imaging for nanoscale materials research. *Appl. Spectrosc. Rev.* **54**, 285–305 (2019).

45. Duncan, K. E., Czymmek, K. J., Jiang, N., Thies, A. C. & Topp, C. N. X-ray microscopy enables multiscale high-resolution 3D imaging of plant cells, tissues, and organs. *Plant Physiol.* **188**, 831–845 (2022).

46. Marone, F., Studer, A., Billich, H., Sala, L. & Stampanoni, M. Towards on-the-fly data post-processing for real-time tomographic imaging at TOMCAT. *Adv. Struct. Chem. Imag.* **3**, 1 (2017).

47. Walsh, C. L. et al. Imaging intact human organs with local resolution of cellular structures using hierarchical phase-contrast tomography. *Nat. Methods* **18**, 1532–1541 (2021).

48. Goldstein, J. *Scanning Electron Microscopy and X-ray Microanalysis* (Springer, 2017).

49. Solé, V. A., Papillon, E., Cotte, M., Walter, Ph. & Susini, J. A multiplatform code for the analysis of energy-dispersive X-ray fluorescence spectra. *Spectrochim. Acta Part B* **62**, 63–68 (2007).

50. Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. *Nat. Methods* **9**, 676–682 (2012).

51. Pettersen, E. F. et al. UCSF Chimera: A visualization system for exploratory research and analysis. *J. Comput. Chem.* **25**, 1605–1612 (2004).

52. Zeng, G. L. *Medical Image Reconstruction: A Conceptual Tutorial* (Springer, 2010).

53. Hsieh, J. *Computed Tomography* 2nd edn. (SPIE, 2009). https://doi.org/10.1117/3.817303.

54. Dempster, A. P., Laird, N. M. & Rubin, D. B. Maximum likelihood from incomplete data via the EM algorithm. *J. R. Stat. Soc.: Ser. B (Methodol.)* **39**, 1–22 (1977).

55. Kazantsev, D., Guo, E., Phillion, A. B., Withers, P. J. & Lee, P. D. Model-based iterative reconstruction using higher-order regularization of dynamic synchrotron data. *Meas. Sci. Technol.* **28**, 094004 (2017).

56. Bergamaschi, A., Medjoubi, K., Messaoudi, C., Marco, S. & Somogyi, A. *MMX-I: Data-processing software for multimodal X-ray imaging and tomography*. *J. Synchrotron. Rad* **23**, 783–794 (2016).

57. Hunter, D. R. & Lange, K. A tutorial on MM algorithms. *Am. Stat.* **58**, 30–37 (2004).

58. Pafilis, C. et al. A methodology for the estimation of the optimal iteration in MLEM-based image reconstruction in PET. In *2011 10th International Workshop on Biomedical Engineering* 1–4 (IEEE, 2011). https://doi.org/10.1109/IWBE.2011.6079046.

59. Epstein, C. L. *Introduction to the Mathematics of Medical Imaging* (Pearson Education/Prentice Hall, 2003).

60. Koho, S. et al. Fourier ring correlation simplifies image restoration in fluorescence microscopy. *Nat. Commun.* **10**, 3103 (2019).

61. van Heel, M. & Schatz, M. Fourier shell correlation threshold criteria. *J. Struct. Biol.* **151**, 250–262 (2005).

62. Banterle, N., Bui, K. H., Lemke, E. A. & Beck, M. Fourier ring correlation as a resolution criterion for super-resolution microscopy. *J. Struct. Biol.* **183**, 363–367 (2013).

63. de Jonge, M. D. et al. Spiral scanning X-ray fluorescence computed tomography. *Opt. Express* **25**, 23424 (2017).

## Acknowledgements

The authors wish to thank E. Chouzenoux (Centre de Vision Numerique, CentraleSupélec, 91190 Gif-sur-Yvette, France) for her help in the validation process of the reconstruction algorithm and G. Correc (Synchrotron Soleil) for his help during the X-ray experiments. The measurements have been performed at the Nanoscopium beamline of Synchrotron Soleil during consecutive in-house research projects and user experiments (20201481 and 20191993).

## Author contributions

Supervision of the project: K.M. Sample preparation: E.B., E.L., M.P.I. Synchrotron measurements, design, preparation, and participation: R.G., D.B., E.B., E.L., M.P.I., K.M., A.S. Development and integration of tomography reconstruction algorithms: R.G., K.M. Provision of the wild-type *Arabidopsis thaliana* seeds: C.C. Validation of the reconstruction algorithms: R.G., K.M. Writing—original draft: R.G., K.M., A.S. with discussion and comments from the other authors.

## Funding

Agence Nationale de la Recherche (Project 'DEFIMAN', ANR-19-CE20-0009-02). French Society of Nephrology (SFNDT) financial support for the research dedicated to kidney calcifications.

## Competing interests

The authors declare no competing interests.
**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-21368-0.

**Correspondence** and requests for materials should be addressed to A.S. or K.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022