 Software
 Open Access
 Published:
ODTbrain: a Python library for fullview, dense diffraction tomography
BMC Bioinformaticsvolume 16, Article number: 367 (2015)
Abstract
Background
Analyzing the threedimensional (3D) refractive index distribution of a single cell makes it possible to describe and characterize its inner structure in a markerfree manner. A dense, fullview tomographic data set is a set of images of a cell acquired for multiple rotational positions, densely distributed from 0 to 360 degrees. The reconstruction is commonly realized by projection tomography, which is based on the inversion of the Radon transform. The reconstruction quality of projection tomography is greatly improved when first order scattering, which becomes relevant when the imaging wavelength is comparable to the characteristic object size, is taken into account. This advanced reconstruction technique is called diffraction tomography. While many implementations of projection tomography are available today, there is no publicly available implementation of diffraction tomography so far.
Results
We present a Python library that implements the backpropagation algorithm for diffraction tomography in 3D. By establishing benchmarks based on finitedifference timedomain (FDTD) simulations, we showcase the superiority of the backpropagation algorithm over the backprojection algorithm. Furthermore, we discuss how measurment parameters influence the reconstructed refractive index distribution and we also give insights into the applicability of diffraction tomography to biological cells.
Conclusion
The present software library contains a robust implementation of the backpropagation algorithm. The algorithm is ideally suited for the application to biological cells. Furthermore, the implementation is a dropin replacement for the classical backprojection algorithm and is made available to the large user community of the Python programming language.
Background
The measurement of the refractive index of a biological cell is always connected to the observable phase change of light as it passes through the cell. For example, phase contrast microscopy in combination with refractive index matching is used to obtain an average value for the refractive index of a cell population [1]. Accuracy and flexibility of this approach are greatly improved with quantitative phase imaging techniques, such as digital holographic microscopy (DHM) [2]. Nevertheless, the reconstruction of the refractive index with subcellular resolution and in three dimensions requires a tomographic approach.
Optical projection tomography
Optical projection tomography (OPT) is a wellstudied technique that is used to reconstruct and quantify volumetric data of biological specimens [3, 4]. What distinguishes OPT from conventional 3D imaging techniques in biology, such as e.g. confocal microscopy or selective plane illumination microscopy (SPIM), is the type of measurement data that is acquired and the way it is processed to reconstruct the 3D measurement volume. Tomographic data sets (sinograms) consist of multiple angular projections of the 3D specimen, whereas the techniques mentioned above directly acquire sclices of the reconstruction volume. In OPTbased refractive index measurements, the projections from multiple angles are quantitative phase images. From these phase projections, the 3D refractive index can be reconstructed by means of the filtered backprojection algorithm [5, 6], an efficient algorithm that computes the inverse Radon transform. In the context of this paper, we refer to the filtered backprojection algorithm as the Radon approximation, because it assumes that light travels along straight lines. The Radon approximation is only valid for small wavelengths (e.g. xray radiation). Thus, when light with wavelengths in the visible spectrum (e.g. 400–600 nm) propagates through a biological cell, with structures on similar length scales, diffraction takes place and significant interference effects emerge.
Optical diffraction tomography
Optical diffraction tomography (ODT) takes this wave nature of light into account. The diffraction of waves at objects is described by the Helmholtz equation. To simplify the complex description of waves by the Helmholtz equation, two known approximations are commonly applied; the Born and Rytov approximations. While the Born approximation assumes that only a small fraction of the wave is diffracted by the cell, the Rytov approximation assumes that the local refractive index variation inside the cell is small. The Rytov approximation is known to be superior to the Born approximation in ODT [7], but its implementation requires a phase unwrapping step that is not present in the implementation of the Born approximation. Note that OPT with the Radon approximation requires a similar phase unwrapping step. Furthermore, OPT operates on the measured phase data only, whereas ODT additionally incorporates the amplitude data in the reconstruction process. The principles of ODT were first introduced in [8] and a backpropagation algorithm was then described in [9]. The theoretical basis of ODTbrain, as well as a derivation of the full 3D backpropagation algorithm in the Rytov approximation are given in [6].
Our contribution
We present what to our knowledge is the first publicly available software implementation of the 3D backpropagation algorithm. The algorithm can be used to reconstruct 3D refractive index maps from projections of biological or artificial phase objects. We showcase our implementation by reconstructing simulated cells and comparing the reconstruction qualities of OPT (Radon) and ODT (Born, Rytov). Furthermore, we investigate the contribution of several physical parameters to the overall reconstruction quality, complement previously conducted twodimensional (2D) studies [7, 10, 11] with 3D data, and draw conclusions concerning the validity of the backpropagation algorithm in three dimensions.
Implementation
The ODTbrain library is implemented in the Python programming language. To make the reconstruction process more transparent and to allow userdefined modifications of the reconstruction, we split the reconstruction process into three steps: a) apply a filter to the complex wave sinogram that corresponds to the required approximation (Born, Radon, or Rytov), b) reconstruct the object data from the sinogram (backproject or backpropagate), c) compute the refractive index distribution from the obtained object data. A full documentation, including a method reference and multiple examples, is available at the project home page.
The filtering step (a) is necessary for the application of the Radon and Rytov approximations. For 3D data sets, we rely on the 2D phaseunwrapping algorithm described in [12]. Step (b) is the main part of the reconstruction algorithm. To speed up the computation of Fourier transforms for the 3D reconstruction, ODTbrain employs the FFTW library [13].
In practice, for a successfull reconstruction a preprocessing step is required: the presented implementation requires a sinogram whose projections are distributed equally from 0 to 360° and a rotational axis that is located at the center of each projection. Phase retrieval and image alignment are not part of the reconstruction process and therefore not discussed in this paper.
Results and discussion
Finitedifference timedomain simulations
In order to test the ODTbrain library, we created tomographic data sets from artificial cell phantoms. We performed 2D and 3D finitedifference timedomain (FDTD) simulations with the software MEEP [14] to obtain in silico sinograms containing phase and amplitude of the scattered wave field. We decided to use the FDTD technique, because it describes the vectorial propagation of light according to the Maxwell equations, allowing us to test the approximations against the most realistic in silico data sets for optical tomography. The MEEP implementation of FDTD is also very robust, easy to use, and thus more practical than simulations based on the exact generalized LorenzMie theory [15, 16]. In order to obtain an image as it would be seen through a microscope, we numerically autofocus the measured field from one wavelength (1 λ) behind the phantom to the center of the phantom by minimizing the gradient of the field amplitude [17–19]. Note that the resulting field at the detector is effectively measured with a numerical aperture of N A=1. Thus, the FDTD simulation technique is ideal for the investigation of arbitrary cell phantom geometries.
An exemplary 2D FDTD simulation is shown in Fig. 1. The figure depicts the 2D cell phantom (Fig. 1 a) and the phase of the sinogram which was computed using FDTD simulations (Fig. 1 b). The values of the refractive index for cytoplasm, nucleus, and nucleolus were taken from [20]. A detailed description of the 2D cell phantom is given in Additional file 1. The 2D reconstructions shown in Fig. 1 c, d, e and the corresponding line plots in Fig. 1 f illustrate the superiority of the Rytov approximation over the Born and Radon approximations. The Born approximation entirely fails to reproduce absolute refractive index values [7, 10, 11] and the Radon approximation produces smearedout refractive index distributions, which becomes visible when comparing the shape of the nucleolus (red) in Fig. 1 d and e with a. The reconstruction with the Rytov and Born approximations were performed with ODTbrain, while the reonstruction with the Radon approximation was performed with the backprojection algorithm as implemented in [21]. The Rytov approximation yields the best reproduction of the initial phantom.
The 3D analog to Fig. 1 is shown in Fig. 2. The cell phantom is rotated about the yaxis and the resulting 3D sinogram consists of 2D projections for the different rotational positions ϕ of the cell phantom. The line plots through the nucleolus of the reconstructed volume in Fig. 2 g and h show the same trend as the 2D line plots in Fig. 1 f. A visualization of the 3D cell phantom is given in the Additional file 2. The reconstructed slice shown in Fig. 2 f exhibits a directional blur that is visible at the upper and lower perimeter of the cell. This directional blur is the result of an incomplete data coverage in the object spectrum commonly referred to as missing applecore artifacts [22]. Furthermore, it is known that the values at the center of the reconstruction volume along the axis of rotation are errorprone due to coherent noise in the measured sinogram [17]. We documented the artifacts resulting from the incomplete object spectrum and from the noise introduced during simulation and reconstruction in Additional file 3. As in the 2D case, the Rytov approximation yields the best reconstruction in the 3D case.
Reconstruction quality measures
In order to obtain a measure for the reconstruction quality, we introduce two metrics; the normalized rootmeansquare (RMS) error and the normalized totalvariation (TV) error. We compute the normalized RMS error E_{RMS} according to
where n _{rec} is the realvalued refractive index distribution of the reconstruction and n _{ph} that of the phantom. The summation \(\sum _{\text {vol}}\) runs over all points that are within the circular (2D) or spherical (3D) reconstruction volume. Points that are outside of this reconstruction volume do not have contributions from all projections and are thus not considered. The RMS error quantifies the error in the refractive index at each point of the reconstruction.
As briefly mentioned in the previous section, the reconstruction with the Radon approximation leads to a blurred reconstruction. Furthermore, the Born approximation cannot reconstruct the correct magnitude of the refractive index values, but produces less blurry images than the Radon approximation. We quantify this blurriness by introducing the TV norm of the differences between phantom n _{ph} and reconstruction n _{rec}. We compute the normalized TV error E_{TV} according to
with an averaged TV norm computed for N = 2 or N = 3 dimensions according to the dimensionality of the simulation
The normalization \(\sum _{\text {vol}} \left (n_{\text {ph}}1\right)^{2}\) implies a normalized RMS error of 100 %, when the reconstruction n _{rec} deviates from the phantom n _{ph} at magnitudes that are comparable to the difference between the phantom n _{ph} and a refractive index value of 1. In the same manner, the normalized TV error becomes 100 % if the averaged TV norm is equal to the squared difference between phantom n _{ph} and 1. The chosen RMS and TV metrics lead to error values that allow a direct comparison between the 2D and 3D reconstruction algorithms (see below).
Quality dependence on total number of projections
In order to determine the number of projections that are necessary for fullview and dense diffraction tomography, we performed 2D and 3D simulations with varying total numbers of projections, equally distributed between 0 and 360°. The cell phantoms used for the 2D and 3D simulations are those described above. The RMS and TV error in dependence of the total number of projections are shown in Fig. 3 a and b. At about 160 projections, both RMS and TV errors reach a plateau, which suggests that above this point, the reconstruction cannot be improved any further. For this reason, we set the number of projections to 200 for all subsequent simulations.
Besides the determination of the required number of projections, one can observe two other important facts in Fig. 3. First, the Rytov approximation results in the best reconstruction in 2D as well as in 3D. Second, the RMS and TV errors for the 3D reconstruction follow the RMS and TV errors of the 2D reconstruction. Apparently, the choice of the two error metrics defined in Eqs. 1 (RMS) and 2 (TV) results in similar error estimates for 2D and 3D reconstructions, independent of the number of projections. This similarity of the 2D and 3D errors suggests that the error of a 3D reconstruction may be inferred from the error of a 2D reconstruction. Even though no proof exists for such a connection, this observation potentially saves time or, in the case of limited computational resources, allows to make error estimates in 3D possible at all (see below).
Quality dependence on refractive index variation
In all the cases discussed above, we used fixed refractive index values for the cell phantom. Here, we show how the reconstruction quality varies with different refractive index values for cytoplasm, nucleus and nucleolus. Figure 4 a shows the different combinations of intracellular refractive index values (see figure caption). In the first simulation, the refractive index values of cytoplasm, nucleus, and nucleolus were set to a fixed value of just above water. The resulting RMS and TV errors are very small (Fig. 4 b and c, first simulation) which results in a good reconstruction of the homogeneous refractive index distribution (Fig. 4 d). The RMS and TV errors both increase as the refractive index values of the cell phantom increase. At refractive index values that are above the values for biological cells (1.36 – 1.39 [20]), the reconstruction of the cell phantom with the Rytov approximation is still successful (Fig. 4 e). At very large refractive index values (cytoplasm 1.455, nucleus 1.435, nucleolus 1.543), the reconstruction starts to show artifacts, as shown in Fig. 4 f. Quantitative line plots through the nucleolus are shown in Additional file 4. In comparison to Fig. 3, we again observe that the Rytov approximation yields the lowest reconstruction error and that the 3D error follows the 2D error. Furthermore, the Rytov approximation is valid across a large range of refractive index values, even outside of the regime of biological cells.
Quality dependence on object size
The size of the cell phantom influences the reconstruction quality. For fixed resolution and detector size, small cells are more reliably reconstructed than large cells, because less pixel or voxel values have to be reconstructed. To illustrate this fact, we conducted simulations for different diameters of the cell phantom. In order to be able to compare the different reconstructions in terms of RMS and TV errors, the size of the simulation volume was kept constant. Due to limited computational power, a 3D simulation for cell diameters above 20 wavelengths was not feasible. Figure 5 shows the dependence of the reconstruction quality on the size of the cell. As noted in the previous sections, the RMS and TV errors of the 3D reconstruction follow that of the 2D reconstruction for the number of projections (Fig. 3) and for different refractive index distributions (Fig. 4). We make the same observation for different cell sizes (Fig. 5 a, c) when the simulation volume is small (30 λ). Therefore, we assume that the 2D errors depicted in Fig. 5 b and d may serve as an approximation for the 3D errors that we would observe for very large cells or small cell clusters. The graphs show that up until a cell diameter of 60 wavelengths the reconstruction with the Rytov approximation is very accurate (see also Additional file 5). Above this critical size, phase unwrapping errors lead to reconstruction artifacts. In summary, the correct reconstruction below object sizes of 60 wavelengths justifies the application of the Rytov approximation to 3D arrangements like biological cell clusters.
Conclusions
The presented algorithm is an extension to optical projection tomography that takes into account diffraction of light due to the refractive index of the sample. We verified previous reports of the superiority of the Rytov approximation over the Born and Radon approximations in 2D for the 3D case. As predicted by the theory, the Born approximation breaks down quickly for large cells and for large RI changes [6]. We have shown that fullview, dense data sets must contain at least 160 projections to achieve the best possible reconstruction of a biological cell with a diameter of 17 wavelengths or less. In contrast, fewview and sparse data sets usually require some form of regularization (see [23, 24]). To improve reconstruction quality, the presented implementation in ODTbrain can be extended by regularization methods such as missingangle artifact minimization, artifact removal by total variation minimization, and iterative reconstruction algorithms [24, 25].
The refractive index range covered successfully by the Rytov approximation is large and includes the range of values found in most biological cells (up to 1.40). Furthermore, the Rytov approximation yields good results for large objects (diameter ≈50 λ). This versatile validity makes the Rytov approximation interesting for the investigation of complex biological samples like cell clusters or even small embryos but also nonbiological structures such as optical fibers.
Availability and requirements
Project name: ODTbrainProject home page: http://odtbrain.craban.de Operating system: Platform independentProgramming language: PythonOther requirements: Python 2.7 or Python 3.4License: BSD (3 clause)Any restrictions to use by nonacademics: None
Abbreviations
 2D:

Twodimensional
 3D:

Threedimensional
 DHM:

Digital holographic microscopy (quantitative phase imaging technique)
 FDTD:

Finitedifference timedomain (simulation technique used for forward process)
 ODT:

Optical diffraction tomography (e.g. backpropagation)
 OPT:

Optical projection tomography (e.g. backprojection)
 RMS:

Rootmeansquare (used for error estimation)
 TV:

Total variation (used for error estimation)
References
 1
Barer R. Refractometry and interferometry of living cells. JOSA. 1957; 47(6):545–56. doi:10.1364/josa.47.000545.
 2
Schürmann M, Scholze J, Müller P, Chan CJ, Ekpenyong AE, Chalut KJ, et al. Refractive index measurements of single, spherical cells using digital holographic microscopy. In: Paluch EK, editor. Biophysical Methods in Cell Biology. Methods in Cell Biology, vol. 125. 1st ed. USA: Academic: 2015. p. 143–59, doi:10.1016/bs.mcb.2014.10.016.
 3
Sharpe J. Optical projection tomography. Annu Rev Biomed Eng. 2004; 6(1):209–28. doi:10.1146/annurev.bioeng.6.040803.140210.
 4
Charrière F, Pavillon N, Colomb T, Depeursinge C, Heger TJ, Mitchell EAD, et al. Living specimen tomography by digital holographic microscopy: morphometry of testate amoeba. Opt Express. 2006; 14(16):7005–13. doi:10.1364/oe.14.007005.
 5
Kak AC, Slaney M. Principles of Computerized Tomographic Imaging. Philadepphia, USA: SIAM; 2001, p. 327.
 6
Müller P, Schürmann M, Guck J. The theory of diffraction tomography. 2015. ArXiv eprints. http://arxiv.org/abs/1507.00466
 7
Slaney M, Kak AC, Larsen LE. Limitations of imaging with firstorder diffraction tomography. IEEE Trans Microw Theory Tech. 1984; 32(8):860–74. doi:10.1109/TMTT.1984.1132783.
 8
Wolf E. Threedimensional structure determination of semitransparent objects from holographic data. Opt Commun. 1969; 1(4):153–6. doi:10.1016/00304018(69)900522.
 9
Devaney AJ. A filtered backpropagation algorithm for diffraction tomography. Ultrason Imaging. 1982; 4(4):336–50. doi:10.1016/01617346(82)900177.
 10
Wedberg TC, Stamnes JJ, Singer W. Comparison of the filtered backpropagation and the filtered backprojection algorithms for quantitative tomography. Appl Opt. 1995; 34(28):6575–81. doi:http://dx.doi.org/10.1364/ao.34.006575.
 11
Chen B, Stamnes JJ. Validity of diffraction tomography based on the first born and the first Rytov approximations. Appl Opt. 1998; 37(14):2996–3006. doi:10.1364/AO.37.002996.
 12
Herráez MA, Burton DR, Lalor MJ, Gdeisat MA. Fast twodimensional phaseunwrapping algorithm based on sorting by reliability following a noncontinuous path. Appl Opt. 2002; 41(35):7437–44. doi:10.1364/AO.41.007437.
 13
Frigo M, Johnson SG. FFTW: an adaptive software architecture for the FFT. In: Acoustics, Speech and Signal Processing, 1998. Proceedings of the 1998 IEEE International Conference On, vol. 3. Piscataway, New Jersey, USA: IEEE: 1998. p. 1381–43, doi:10.1109/ICASSP.1998.681704.
 14
Oskooi AF, Roundy D, Ibanescu M, Bermel P, Joannopoulos JD, Johnson SG. MEEP: A flexible freesoftware package for electromagnetic simulations by the FDTD method. Comput Phys Commun. 2010; 181:687–702. doi:10.1016/j.cpc.2009.11.008.
 15
Boyde L, Chalut KJ, Guck J. Near and farfield scattering from arbitrary threedimensional aggregates of coated spheres using parallel computing. Phys Rev E. 2011; 83(2):26701. doi:10.1103/PhysRevE.83.026701.
 16
Boyde L, Ekpenyong A, Whyte G, Guck J. Comparison of stresses on homogeneous spheroids in the optical stretcher computed with geometrical optics and generalized LorenzMie theory. Appl Opt. 2012; 51(33):7934–44. doi:http://dx.doi.org/10.1364/ao.51.007934.
 17
Kostencka J, Kozacki T, Dudek M, Kujawińska M. Noise suppressed optical diffraction tomography with autofocus correction. Opt Express. 2014; 22(5):5731–45. doi:10.1364/OE.22.005731.
 18
Wu SY, Dugan N, Hennelly BM, Vol. 9131. Investigation of autofocus algorithms for Brightfield microscopy of Unstained cells. Bellingham, Washington USA: SPIE; 2014, pp. 91310–12.
 19
Müller P. nrefocus: Python algorithms for numerical focusing (Version 0.1.4) [Software]. 2013. https://pypi.python.org/pypi/nrefocus/. Accessed date 19 September 2015.
 20
Choi WJ, Park KS, Eom TJ, Oh MK, Lee BH. Tomographic imaging of a suspending single live cell using optical tweezercombined fullfield optical coherence tomography. Opt Lett. 2012; 37(14):2784–6. doi:10.1364/OL.37.002784.
 21
Müller P. radontea: Python algorithms for the inversion of the Radon transform (Version 0.1.7) [Software]. 2013. https://pypi.python.org/pypi/radontea/. Accessed date 19 September 2015.
 22
Vertu S, Delaunay JJ, Yamada I, Haeberlé O. Diffraction microtomography with sample rotation: influence of a missing apple core in the recorded frequency space. Central Eur J Phys. 2009; 7(1):22–31. doi:10.2478/s1153400801546.
 23
LaRoque SJ, Sidky EY, Pan X. Accurate image reconstruction from fewview and limitedangle data in diffraction tomography. J Opt Soc Am A. 2008; 25(7):1772–82. doi:10.1364/JOSAA.25.001772.
 24
Su JW, Hsu WC, Chou CY, Chang CH, Sung KB. Digital holographic microtomography for highresolution refractive index mapping of live cells. J Biophoton. 2013; 6(5):416–24. doi:10.1002/jbio.201200022.
 25
Behrooz A, Zhou HM, Eftekhar AA, Adibi A. Total variation regularization for 3D reconstruction in fluorescence tomography: experimental phantom studies. Appl Opt. 2012; 51(34):8216–27. doi:http://dx.doi.org/10.1364/ao.51.008216.
Acknowledgements
We acknowledge helpful discussions with Ivo Sbalzarini, Pavel Tomancak, Moritz Kreysing, Martin Weigert (Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany), and Kevin Chalut (Cavendish Laboratory, University of Cambridge, Cambridge, UK; Wellcome Trust  Medical Research Council, Cambridge Stem Cell Institute, Cambridge, UK). This project has received funding from the European Union’s Seventh Framework Programme for research, technological development and demonstration under grant agreement no 282060. We thank the Dresden International Graduate School for Biomedicine and Bioengineering (DIGSBB) and the German Research Foundation (DFG) for their support with the publication costs in the context of the Excellence Initiative.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
PM designed and implemented the software, performed the simulations, and carried out the data analysis. MS contributed to key mechanisms of the reconstruction methods and revised the documentation. MS and JG contributed to the design of the study. PM and MS drafted the manuscript. JG supervised the study and helped to draft the manuscript. All authors read and approved the final manuscript.
Additional files
Additional file 1
Description of cell phantom for 2D FDTD simulations. The schematic drawing shows the 2D refractive index phantom for the FDTD simulations. The cytoplasm with a refractive index of 1.365 (orange) is centered in the simulation volume (P _{1}=(0,0)) and has major and minor axis of a _{1}=8.5 λ and b _{1}=7.0 λ. The nucleus with a refractive index of 1.360 is positioned offcenter at P _{2}=(2 λ,1 λ) with a _{2}=4.5 λ and b _{2}=3.5 λ. The nucleus is rotated with respect to the coordinate system at a fixed angle of θ=0.5 rad. The nucleolus with a refractive index of 1.387 is positioned at P _{3}=(2 λ,2 λ) with a radius of r _{3}=1 λ. (PNG 774 kb)
Additional file 2
Description of cell phantom for 3D FDTD simulations. The video is a visualization of the 3D refractive index phantom for the FDTD simulations. The values for the refractive indices, sizes, and positions for the 3D phantom are identical to that of the 2D phantom. The cytoplasm and nucleus of the 3D version of the phantom are prolate ellipsoids. (AVI 5140 kb)
Additional file 3
Central crosssections of 3D reconstruction with the Rytov approximation. This figure supplements Fig. 2 with crosssections and line plots through the center of the 3D volume. Along the axis of rotation (highlighted by two white lines), the reconstruction exhibits strong variations. (PNG 1003 kb)
Additional file 4
Line plots of the 3D reconstruction for different magnitudes of the refractive index variation. The figure is a quantitative representation of the three distributions of refractive index highlighted in (Fig. 4 d, e, f) with line plots through the nucleolus parallel to the minor (a, c, e) and the major axis (b, d, f). The reconstruction with the Born, Radon, and Rytov approximations are plotted. (PNG 1280 kb)
Additional file 5
Line plots of the 2D reconstruction for different cell sizes. The figure is a quantitative representation of the three distributions of refractive index highlighted in (Fig. 5 e, f, g) with line plots through the nucleolus as indicated in (Fig. 1 c, d, e). The reconstruction with the Born, Radon, and Rytov approximations are plotted. (PNG 1065 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Refractive index
 Singlecell analysis
 Diffraction tomography
 Backprojection
 Backpropagation
 Rytov
 Born
 Radon