1. Introduction

3D imaging in the life sciences is of great importance for studying fundamental biology and performing (pre-) clinical studies. For these studies, label free optical imaging methods play an important role. There are various label-free contrast mechanisms such as scattering, absorption, or refractive index (RI). However, in some cases these contrast mechanisms are not sufficiently sensitive to observe the relevant information, hence, there is a need for imaging with alternative types of intrinsic contrast.

Optical diffraction tomography (ODT) has shown to be an effective tool for 3D imaging of RI contrast on the scale of cells [1] or small organisms [2]. More recently, phase contrast ODT was applied on a millimeter scale, where different structural features of a zebrafish larva and a cryo-injured heart could be distinguished in 3D using RI contrast [3]. However, some types of tissue are not visible in conventional phase contrast ODT.

An alternative form of contrast is given by the polarization change of the optical wavefield caused by tissue birefringence. Birefringent samples are not described by a single scalar RI value per voxel that contributes to the optical path length, but the RI value experienced by the wavefield depends on its polarization state. Polarization contrast has been widely applied in microscopy [4,5], digital holography [6], optical coherence tomography [7], and optical projection tomography [8]. Birefringence provides a high-constrast label-free mechanism for imaging fibrous structures such as muscle (collagen) or brain (myelin) tissue. Muscle tissue has been imaged in 3D using polarization sensitive optical projection tomography (OPT), as an extension of brightfield OPT using a white light source [8]. However, with OPT phase information is lost and refractive index contrast cannot be determined.

In this work we show that in addition to phase contrast also polarization contrast is compatible with large scale ODT and offers a significantly higher signal to noise ratio (SNR) compared to conventional phase contrast ODT. We determine under what conditions a birefringent sample can be properly reconstructed using conventional filtered backprojection (FBP). Furthermore, we show that off-axis sample placement, which has been used in conventional ODT [9] for noise reduction, also for polarization ODT offers significant noise reduction and that the same steps of numerical refocusing to correct for defocus can be applied. Finally, we demonstrate 3D multi-contrast imaging of a zebrafish larva using two orthogonal components of the transmitted wavefield, from which a conventional phase contrast and polarization contrast ODT image are reconstructed.

2. Polarization contrast imaging

In conventional ODT, refractive index differences in the sample cause a change in optical path length of the transmitted light wave. Assuming an isotropic medium, each voxel in the sample gives a fixed contribution to the optical path length of a ray traveling through it regardless of its’ polarization. However, when a sample is birefringent this contribution generally depends on the orientation of the polarization of the wave with respect to the medium.

Here we use Jones calculus to calculate light interactions. We assume that the birefringent tissue locally can be described as uniaxial, where the optical axis corresponds to the predominant fiber direction. The birefringent tissue is modeled as a wave retarder that introduces a relative phase shift $\delta$ along the fast axis with respect to the slow axis, and introduces a common phase shift $\epsilon$ (i.e. the average phase of the two components) for both polarization components. The relative phase shift $\delta$ between the two components is then defined as [10]

 

Fig. 1. Schematic of the ODT sample arm geometry. (a) The orientation of the uniaxial sample is defined by the inclination angle $\alpha$ and the direction angle $\varphi $. The sample is rotated around the $x$-axis for tomographic measurement. The input polarization state is along the $x$-axis, after which the parallel polarized $x$-component (a) or the cross-polarized $y$-component (b) of the complex wave is measured for each projection angle. The tomographic angle $\beta$ is defined with respect to the fiber orientation in the $y$-$z$ plane. The angle of the polarizers $\rho$ is defined with respect to the $x$ axis.

Download Full Size | PDF

We assume an incoming beam polarized along the $x$-axis that travels through the sample in the $z$-direction. Both the $x$ and $y$ components are extracted by placing an analyzer in the sample arm that can be rotated to align with the parallel $x$ or cross-polarized $y$-axis. The complex wavefield of an incoming wave polarized along the $x$-axis after transmission through the birefringent medium is

2.1 Parallel-polarization output

The first component in Eq. (3) is the $x$-component of the transmitted field with a polarization parallel to that of the input field. It can be extracted by placing a polarizer aligned along the $x$-axis after the sample. The $x$-component in Eq. (3) contains phase contributions of both the conventional phase contrast $\epsilon$ and the birefringence contrast $\delta$. The phase of this component is defined as the inverse tangent of the imaginary part divided by the real part

2.2 Cross-polarization output

The vertical $y$-component is the second component of the field in Eq. (3) and is perpendicular to the input polarization. The amplitude of this component is given by

However, from Eq. (9) it follows that in case the maximum value of $\delta$ in the projection does not exceed $\pi$, the signal can be directly retrieved by taking the inverse sine and no further processing is necessary. Even more, if $\delta$ is small, the amplitude of the $y$-component of Eq. (3) can be approximated as a linear function of $\delta$, since for small values of $\delta$ it holds that

 

Fig. 2. Phase shift $\delta$ between the two orthogonal polarizations in the case of a uniaxial birefringent cylinder with maximum projected phase shift $\delta =18$ radians as quantified with polarization contrast imaging.

Download Full Size | PDF

The red line shows two times the inverse sine of the measured $|$sin($\delta /2$)$|$ term. The green line is obtained by flipping the inverse sine function in the appropriate domains and adding $\pi$ according to Eq. (9). The function $\delta$ can then be retrieved with standard phase unwrapping and is plotted in magenta and corresponds with the original birefringence distribution. Thus, in theory the projection function $\delta$ can be retrieved. However, in practice this may not be possible, for example, when the data is noisy or the jumps in the sinusoidal signal of the transmitted field $U_{y}$ are not properly sampled due to large increase of $\delta$.

2.3 Polarization tomography

In 3D polarization sensitive tomographic imaging, the sample is rotated and the $x$ (parallel) and $y$ (cross) components of the wave are recorded for each angle for phase and polarization contrast respectively. Due to the small contribution of the birefringence contrast in the $x$-component phase it can be used for conventional ODT. However, it should be noted that in order to preserve the linear relationship between the projection and the $y$-component, it can be seen from Eq. (10) that not the intensity (amplitude squared) of the wavefield should be taken as the projection, but the square root of the intensity (amplitude).

However, in general $\delta$ itself depends on the tomographic rotation angle $\beta$ through $\alpha$ in Eq. (1) and Eq. (2). Furthermore, the angle $\varphi$ in Eq. (10) depends on $\beta$ as well through Eq. (2). Using these dependencies we find that for small $\delta$ the $y$-component of the field is

The question arises whether the phase of the crossed-polarizer component can be used to do the conventional phase reconstruction, so that capturing of $U_{x}$ is not necessary. In cross polarization, the phase $\epsilon$ of the transmitted $y$-component is defined for any path through the birefringent sample where the field amplitude is not zero. Hence, this component cannot be used to reconstruct the conventional RI contrast $\epsilon$ across the whole sample. However, the phase of the $y$-component can be used in order to propagate the wavefield. This can be used to numerically refocus the wavefield if necessary, for example in the case of off-axis placement of the sample for noise suppression [3,9], or to extend the depth of field of the imaging system [2].

3. Materials and methods

3.1 Acquisition of projections

In ODT, the scattered field is recorded from multiple angles using digital holography. The digital holography setup is shown in Fig. 3 and consists of a Mach-Zehnder interferometer operated in transmission. The light source is a HeNe laser with a wavelength of 633 nm and an output power of 3 mW. Two lenses (Thorlabs, LD2568 and LA1979) are used to expand and collimate the illuminating laser beam to a full width at half maximum (FWHM) of approximately 15 mm.

 

Fig. 3. Experimental setup for acquiring the digital holograms. HeNe: Helium Neon laser, BE: Beam expander, BS: Beam splitter, IML: Index matching liquid, S: Sample rotated around the z-axis, MO: Microscope objective, M: Mirror, TL: Tube lens, PZT: Mirror mounted on piezo stage, C: Camera, P: Polarizer, WP: Half-wave plate.

Download Full Size | PDF

In the object arm a 10X objective lens (NA=0.3) is used in combination with a 200 mm focal length tube lens (Thorlabs) to image the sample in close proximity to the detector of a CMOS camera, (Basler beA4000-62kc) with $4096 \times 3072$ pixels and a pixel pitch of 5.5 μm. A rotation mount (Thorlabs CR1) rotates the sample stepwise over 360$^{\circ }$. One polarizer is placed in front of the sample (P1), and a second one is placed behind the sample (P2). For acquisition of the regular phase contrast projections, the optical axes of the polarizers are made parallel and a acquisition of 720 projections over 360$^{\circ }$ is performed. For the polarization contrast projections, the relative angle between both polarizers is kept constant at 90$^{\circ }$. The complete tomographic measurement is then carried out as before. The polarization contrast measurement is then repeated after simultaneous rotation of both the polarizers by 30$^{\circ }$ and 60$^{\circ }$, respectively. In the reference arm, a polarizer (P3) is placed in order to maximize the fringe contrast at the detector; this polarizer is rotated simultaneously with the polarizers in the object arm. A half-wave plate is placed behind the beam expander in order to maximize the signal at the detector. In the reference arm, a 10X Olympus microscope objective partly compensates for the object wave curvature to avoid the presence of too high spatial frequencies on the camera. The mirror in the reference arm is mounted onto a piezoelectric transducer (Thorlabs, KPZ 101) controlled by a computer for phase-shifting the digital hologram. We capture four holograms with reference arm phase shift increments of $\pi /2$ between each subsequent hologram. From a linear combination of these holograms a complex hologram is formed where the zeroth and out of focus conjugate orders are removed [11]. In this way we maximize the lateral resolution in the reconstructed image. This is specifically important for large scale ODT where magnification is low but an as high as possible NA is desired.

3.2 Phase and polarization projections

Autofocus correction is applied on the digital hologram in order to obtain the wavefield in the object region. The object position is determined by calculating a focus metric (grayscale variance) as a function of the reconstruction distance. For transparent objects the gray scale variance has a minimum value when the reconstruction distance is located at the object. For polarization contrast projections, the gray scale variance has a maximum value when reconstructed in focus. For both cases separately the minimum/maximum is determined for ten samples of a full rotation acquisition (i.e. 0$^{\circ }$, 36$^{\circ }$, 72$^{\circ }$, etc.). A sinusoidal function is then fitted to the minimum/maximum as a function of the projection angle to determine the object distance as a function of projection angle. For every angle the hologram is reconstructed for both the phase and polarization contrast data, with the object in focus by propagating the field to the object plane using the angular spectrum method for diffraction calculation, which is exact and valid for small propagation distances. In case of the phase projections, the phase is then calculated by taking the argument of the reconstructed wavefield. The phase projections are unwrapped using a least squares phase unwrapping algorithm [12].

For the polarization contrast projections, the amplitude of the cross-polarized component is calculated. This amplitude then gives a direct, but scaled measure for the birefringence: scaled $n_e-n_o$. For the different (cross) polarizations, the projections are misaligned horizontally by a few pixels. This is corrected by determining the center of rotation from the maximum variance of the tomographic reconstruction as a function of the shift for each polarization contrast sinogram individually. The projections are then shifted to the correct location using the circular shift function of MATLAB. The wavefield amplitude of the projections for the three angles are stacked, and the maximum value for each camera coordinate is extracted to form a single maximum birefringence projection sinogram. Tomographic imaging is performed with 720 projections over 360$^{\circ }$ (steps of 0.5$^{\circ }$) with four phase steps per projection. At every projection angle and phase step, four measurements are taken (one for phase, three for polarization contrast) in total. The net acquisition time for a full 3D measurement is approximately 7 minutes with the total acquired data around 160 GB.

3.3 Tomographic image reconstruction and visualization

For reconstruction of the phase contrast, assuming that RI variation in the sample is sufficiently small so that refraction does not occur, a phase projection is a scaled integral over the RI variation with respect to the background medium along the illumination direction. The average refractive index difference $\Delta n_{avg}$ is calculated from the phase by using the system magnification and the pixel pitch [9]. Subsequently, the $\Delta n_{avg}$ object is reconstructed using the FBP algorithm on a slice by slice basis. For polarization contrast, the maximum birefringence projection sinogram $\delta$ is reconstructed using the FBP algorithm as $n_e - n_o$. We used the Drishti software package [13] to visualize and merge the phase and polarization contrast reconstructions with a non-linear transfer function.

3.4 Noise suppression in polarization sensitive ODT

The sample is displaced from the center of rotation by approximately 0.5 mm. Figure 4 shows the noise distribution, in standard deviation $\sigma$, in a tomographic ODT reconstruction of both the polarization contrast (a) and (b) and the phase contrast (c) and (d). The polarization contrast ODT reconstruction suffers from increased noise in the region of the center of rotation, similar to what has been shown to be the case with phase contrast ODT. The noise at the center of rotation is approximately a factor 7 higher than outside of the center. This also shows that the on-axis noise reduction is even more significant in the case of polarization contrast ODT than in phase contrast ODT, where the noise reduction by off-axis placement was found to be in the order of a factor 2 for 720 projections [3].

 

Fig. 4. (a) Logarithm of the standard deviation $\sigma (n_e - n_o)$ of a single polarization contrast reconstructed slice. (b) Cross-section along the dashed line in figure (a) and the average standard deviation over all slices of the stack (red). (c) Logarithm of the standard deviation $\sigma (\Delta n_{avg})$ of the phase contrast reconstructed slice. (d) Cross-section along the dashed line in figure (c) and the average standard deviation over all slices of the stack (red).

Download Full Size | PDF

3.5 Zebrafish sample preparation

The sample is a 3 day old zebrafish embryo (wild type). The eggs are grown on a petridish and subsequently placed in PTU (1-phenyl 2-thiourea) to prevent pigment formation. At 72 hours, the eggs are dechorionated and fixated in 4% paraformaldehyde. Then, the eggs are washed with Phosphate buffered saline three times, after which it is replaced with 100% MeOH in two cycles for dehydration. The embryos are placed in small cylinders (4 mm diameter) and mixed with agarose (2% mass-percentage). After the agarose is dry, the agarose containing the embryo’s is removed from the cylinders and as a whole placed in BABB, a mixture of benzyl alcohol (Sigma B-1042) and benzyl benzoate (Sigma B-6630) in a 1:2 ratio, which makes the sample completely transparent [14]. During this process, the RI of the sample becomes almost that of the BABB clearing solution. We used a clearing time of 3 hours (similar to [3]) that ensures that the sample is transparent enough for optical phase tomography, while at the same time maximizing remaining RI contrast in order to keep a good signal (RI contrast in the reconstruction) to noise (background) ratio in the final reconstruction.

4. Results

The polarization and phase contrast projections of a 3 day old zebrafish tail are shown in Fig. 5(a)-(b) and (d)-(e), respectively. The phase contrast projections are similar to our earlier work on ODT applied to zebrafish larvae [3]. In the polarization contrast projections most of the larva appears dark, due to the absence of birefringent tissue, except in the tail where the developing highly birefringent muscle tissue (myotome) is located. The polarization contrast results are found to be similar in comparison with 2D polarization contrast measurements of Jacoby et al. [15]. The histograms of the 3D polarization and phase contrast reconstructions are shown in Fig. 5(c) and (f) respectively. The polarization contrast histogram of the scaled birefringence shows two components, namely the background and the myotome tissue. In the phase contrast histogram of the polarization averaged refractive index multiple peaks, corresponding to different organs, are visible [3].

 

Fig. 5. Reconstructed amplitude (a) and (b) and phase projections (d) and (e) from two different angles of a 3 day old optically cleared zebrafish larva, illustrating the different contrasts obtained through polarization and phase contrast respectively. In (c) and (f) the histograms of the full 3D data set are plotted for the polarization and phase contrasts, respectively. The background contribution is indicated in both histograms, and the myotome and interstitial tissue for the polarization and phase contrast respectively.

Download Full Size | PDF

A 3D visualization of the phase contrast, the polarization contrast, the merged datasets and transverse cross-sections after tomographic reconstruction using FBP are shown in Fig. 6. It can be clearly seen from the visibility of the developing muscle tissue (myotome) that the phase and polarization contrast offer complementary contrasts, even though they spatially overlap. The anatomical structures are annotated based on reference data from microscopy [15] and OPT [16]. A striking result is the high contrast obtained in the polarization contrast projections compared to the phase projections. We quantify this by calculating the standard deviation of a background region outside of the center (since the level of noise is lower there), and estimate the mean of the signal in the tail at the same location for both the polarization and phase contrast reconstructions. For polarization contrast ODT, this yields an SNR of approximately SNR=$73$, and for the phase contrast ODT we obtain an SNR of approximately SNR=$15$. Polarization contrast ODT thus yields significantly higher SNR than phase contrast ODT for imaging the zebrafish tail.

 

Fig. 6. 3D visualization of the phase (a) and polarization (b) contrast, and combined (c) ODT reconstructions of a 3 day old zebrafish larva tail. In phase contrast, the tail (in red) and the spinal cord (in purple) appear, but not the developing muscle tissue (myotome), which is birefringent. In the polarization contrast reconstruction the structure of the myotome can be clearly discerned. Insets show transverse cross sections in linear intensity scale taken at the dashed line. Scalebar for 3D reconstruction corresponds to 200 µm.

Download Full Size | PDF

5. Discussion and conclusion

We demonstrate 3D polarization contrast ODT, which has previously only been achieved only with OPT. Applying it within the framework of ODT makes it possible to image both phase and polarization contrast and make use of the benefits of ODT such as numerical refocusing and extended depth of field, due to the fact that both phase and amplitude of the polarization contrast field are measured.

5.1 Polarization ODT contrast

Coherent speckle causes increased noise levels close to the center of rotation in polarization contrast ODT similar as in conventional phase contrast ODT and the same strategy of off-axis placement and numerical refocusing can be applied to reduce the noise level up to a factor of 7. The polarization contrast ODT reconstruction yields a significantly higher signal to noise ratio compared to the phase contrast reconstruction. We attribute this to the fact that in phase contrast ODT the refractive index differences decrease during clearing, leading to a reduction of the signal to noise ratio in the reconstructed images. For polarization contrast ODT, the background is zero (no transmission in the absence of birefringence) and consequently leads to a relatively high contrast when birefringent tissue is present. Besides this qualitative argument, also quantitatively, the value of the average refractive index, which is proportional to $n_e + n_o$, and the birefringence $n_e - n_o$ may vary during the clearing process [17] and thus influence the image contrast in both ODT modes.

5.2 Limit on maximum projected $\delta$

Straightforward tomographic reconstruction only yields valid results for polarization contrast ODT in case $\delta$ is small. In phase projections of highly birefringent materials, such as a FEP (fluorinated ethylene propylene) tube, phase wrapping is clearly visible as a dense amplitude modulation. For cleared biological samples we have not observed dense amplitude modulation and, for all practical purposes, the wrapping problem is absent. Even for uncleared samples with 0.5 mm of birefringent tissue, phase wrapping is absent for birefringence lower than $n_e - n_o=6\cdot 10^{-4}$, which is still smaller than the typical birefringence of uncleared tissue [7].

For application outside of biomedicine, the wrapping of $\delta$ places a practical limitation on the amount of birefringence and/or the maximum sample thickness that can be imaged using conventional reconstruction. In principle, the correct projection and reconstruction can be retrieved in case the linearity requirement is violated using a modified unwrapping procedure based on the forward model. However, further research is needed for application of this procedure on experimental data.

5.3 Absolute quantification of birefringence

A limitation of the current method is that the polarization contrast is qualitative. Absolute quantification of the birefringence is challenging as the magnitude of the signal is dependent on the incident field distribution, sample optical absorption, the light to electron conversion, and the fiber orientation. In principle the first three factors can be divided out using a reference measurement, e.g., from the amplitude of the parallel polarization state projection.

A further complication comes from the tomographic angle dependence of $\delta$ that causes a modulation outside of a continuous region of birefringence. In case the macroscopic assumption of uniform birefringence across a region does not apply, but the fiber orientation changes significantly on small length scales this may cause reconstruction artifacts. The case for quantitative birefringence tomography (quantification of optic axis, $n_e$, and $n_o$) is more complicated as it requires more information per projection angle and a non-linear inversion scheme. This is outside of the scope of the current work.

5.4 Applicability of the uniaxial model

The analysis and simulations in this paper are based on the assumption of uniaxial birefringence. The uniaxial model is a simple and widely used model in polarization microscopy, and applicable to fibrous structures such as myelin, elastin, and collagen. The well-defined fiber orientation as is present in the uniaxial model would be of importance to extract from the data. Further research is needed to determine whether fiber orientation can be retrieved in 3D, for example by performing more measurements under different input polarizations and using a full vectorial reconstruction [18].

Although the uniaxial model works for a large class of tissues, some types of tissues exhibit biaxial birefringence [19]. In addition, in some voxels there may be overlapping tissue fibers. Incorporating this in the tomographic reconstruction requires a more elaborate birefringence model.

5.5 Conclusion

We demonstrated 3D polarization contrast ODT. The developing muscle tissue in the tail of the zebrafish larva is known to be birefringent and cannot be discerned in conventional phase contrast ODT reconstruction. By illuminating the sample with a single polarization input state and measuring both the parallel (for the phase) and the orthogonal component (for the polarization contrast) with digital holography a conventional and polarization contrast ODT reconstruction of the same object can be obtained.

A. Appendix

Here we demonstrate the effect of the angular dependency in the amplitude projection on the tomographic reconstruction. The object we consider is a cylinder according to the orientation outlined in the theory section of this paper. The cylinder has radius $R$, and birefringence $n_{e}-n_{o}=\delta n$. For a plane wave traveling along the $z$-axis linearly polarized along the $x$-axis, the polarization contrast is retrieved from the cross-polarized transmitted component after traveling through the sample. This component $U_{y}$ is given by the $y$ component of Eq. (A2). Using the relations $\delta =k \Delta \cos ^2(\alpha (\beta ))$, $\alpha =\gamma \sin (\beta )$ and $\varphi =\gamma \cos (\beta )$, the full expression for $U_{y}$ on $\beta$ becomes

(A1)$$U_{y}(\beta)={-}i\textrm{e}^{i \epsilon} \sin (2\rho - 2 \gamma \cos (\beta )) \sin \left(\frac{1}{2} k \Delta \cos ^2(\gamma \sin (\beta ))\right) \, ,$$

and the amplitude of $U_{y}(\beta )$ is

(A2)$$|U_{y}(\beta)|=\left| \sin (2\rho - 2\gamma \cos (\beta )) \right| \left| \sin \left(\frac{1}{2} k \Delta \cos ^2(\gamma \sin (\beta ))\right)\right| \, .$$

For a cylinder located at the origin with a tilt $\gamma$ with respect to the $x$-axis of tomographic rotation, the cross-section seen by a wave traveling along the $z$-axis is an ellipse $f(y,z)$, with semi-major and semi-minor axes $a=R \sec (\gamma )$ and $b=R$, respectively. The Radon transform $\Re (f)$ for a 2D slice of the ellipse gives the path length experienced by the probing wave per projection angle $\beta$ and is given by [20]

(A3)$$\Re(f)= \begin{cases} \frac{2 R^2 \sec (\gamma ) \sqrt{A-p^2}}{A} & p^2\leq A \\ 0 & \textrm{otherwise} \end{cases}$$

where

(A4)$$A=R^2 \cos ^2(\beta ) \sec^2(\gamma )+R^2 \sin^2(\beta ) \, ,$$

and $p$ is the transverse coordinate along the projection. Replacing $\Delta$ in Eq. (A2) with $\Re (f) \delta n$, the effective amplitude projection function measured at the detector becomes

(A5)$$|U_{y}(p,\beta)|=\begin{cases} \left| \sin (2 \rho -2 \gamma \cos (\beta )) \sin \left(\frac{R^{2}\,\delta n\,k \sec (\gamma ) \cos ^2(\gamma \sin (\beta )) \sqrt{A-p^2}}{A}\right)\right| & p^2\leq A \\ 0 & \textrm{otherwise} \end{cases}$$

The projection functions $|U_{y}(p,\beta )|$ along with the resulting tomographic reconstructions are plotted in Fig. 7 for tilt angles $\gamma =0^{\circ }$ (a-b) and $\gamma =54^{\circ }$ (c-d). The simulation parameters are cross-polarizer angle $\rho =27^{\circ }$, $\delta n=1\cdot 10^{-5}$, $R=1$ mm and $\lambda =633 \cdot 10^{-9}$. For comparison, the case for a non-birefringent cylinder at $\gamma =54^{\circ }$ is shown (e-f). It can be seen that the angular dependency of the amplitude projections results in a modulation along the horizontal projection angle axis in Fig. 7(c). Since this does not cause modulation along the transverse coordinate axis and the amplitude is zero outside of the projection of the birefringent object, the contrast inside the birefringent sample is not modulated (Fig. 7(d)). Instead, it gives a slowly varying angular modulation in the background.

 

Fig. 7. Plot of the projection functions $|U_{y}(p,\beta )|$ along with the resulting tomographic reconstructions for tilt angles $\gamma =0^{\circ }$ (a-b) and $\gamma =54^{\circ }$ (c-d). The simulation parameters are cross-polarizer angle $\rho =27^{\circ }$, $\delta n=1 \cdot 10^{-5}$, $R=1$ mm and $\lambda =633 \cdot 10^{-9}$. For comparison, the case for $\gamma =54^{\circ }$ for a non-birefringent cylinder is shown (e-f).

Download Full Size | PDF

Acknowledgements

We would like to thank Sonja Chocron and Jeroen Bakkers from the Hubrecht Institute for providing zebrafish larvae. We thank Dr. Miriam Menzel for useful discussions.

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. K. Kim, K. S. Kim, H. Park, J. C. Ye, and Y. Park, “Real-time visualization of 3-D dynamic microscopic objects using optical diffraction tomography,” Opt. Express 21(26), 32269 (2013). [CrossRef]  

2. W. Choi, C. Fang-Yen, K. Badizadegan, R. R. Dasari, and M. S. Feld, “Extended depth of focus in tomographic phase microscopy using a propagation algorithm,” Opt. Lett. 33(2), 171 (2008). [CrossRef]  

3. J. van Rooij and J. Kalkman, “Large-scale high-sensitivity optical diffraction tomography of zebrafish,” Biomed. Opt. Express 10(4), 1782 (2019). [CrossRef]  

4. R. Oldenbourg, “A new view on polarization microscopy,” Nature 381(6585), 811–812 (1996). [CrossRef]  

5. F. Massoumia, R. Juškaitis, M. A. A. Neil, and T. Wilson, “Quantitative polarized light microscopy,” J. Microsc. 209(1), 13–22 (2003). [CrossRef]  

6. Z. Wang, L. J. Millet, M. U. Gillette, and G. Popescu, “Jones phase microscopy of transparent and anisotropic samples,” Opt. Lett. 33(11), 1270 (2008). [CrossRef]  

7. M. J. Everett, K. Schoenenberger, B. W. Colston, and L. B. D. Silva, “Birefringence characterization of biological tissue by use of optical coherence tomography,” Opt. Lett. 23(3), 228 (1998). [CrossRef]  

8. M. Fang, D. Dong, C. Zeng, X. Liang, X. Yang, A. Arranz, J. Ripoll, H. Hui, and J. Tian, “Polarization-sensitive optical projection tomography for muscle fiber imaging,” Sci. Rep. 6(1), 19241 (2016). [CrossRef]  

9. J. Kostencka, T. Kozacki, M. Dudek, and M. Kujawińska, “Noise suppressed optical diffraction tomography with autofocus correction,” Opt. Express 22(5), 5731–5745 (2014). [CrossRef]  

10. M. Menzel, K. Michielsen, H. D. Raedt, J. Reckfort, K. Amunts, and M. Axer, “A Jones matrix formalism for simulating three-dimensional polarized light imaging of brain tissue,” J. R. Soc., Interface 12(111), 20150734 (2015). [CrossRef]  

11. I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268 (1997). [CrossRef]  

12. D. C. Ghiglia and L. A. Romero, “Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods,” J. Opt. Soc. Am. A 11(1), 107 (1994). [CrossRef]  

13. A. Limaye, “Drishti: a volume exploration and presentation tool,” Proc. SPIE 8506, 85060X (2012). [CrossRef]  

14. A. d’Esposito, D. Nikitichev, A. Desjardins, S. Walker-Samuel, and M. F. Lythgoea, “Quantification of light attenuation in optically cleared mouse brains,” J. Biomed. Opt. 20(8), 080503 (2015). [CrossRef]  

15. A. S. Jacoby, E. Busch-Nentwich, R. J. Bryson-Richardson, T. E. Hall, J. Berger, S. Berger, C. Sonntag, C. Sachs, R. Geisler, D. L. Stemple, and P. D. Currie, “The zebrafish dystrophic mutant softy maintains muscle fibre viability despite basement membrane rupture and muscle detachment,” Development 136(19), 3367–3376 (2009). [CrossRef]  

16. R. J. Bryson-Richardson, S. Berger, T. F. Schilling, T. E. Hall, N. J. Cole, A. J. Gibson, J. Sharpe, and P. D. Currie, “Fishnet: an online database of zebrafish anatomy,” BMC Biol. 5(1), 34 (2007). [CrossRef]  

17. D. Chen, N. Zeng, Q. Xie, H. He, V. V. Tuchin, and H. Ma, “Mueller matrix polarimetry for characterizing microstructural variation of nude mouse skin during tissue optical clearing,” Biomed. Opt. Express 8(8), 3559 (2017). [CrossRef]  

18. V. Lauer, “New approach to optical diffraction tomography yielding a vector equation of diffraction tomography and a novel tomographic microscope,” J. Microsc. 205(2), 165–176 (2002). [CrossRef]  

19. M. Ravanfar and G. Yao, “Measurement of biaxial optical birefringence in articular cartilage,” Appl. Opt. 58(8), 2021 (2019). [CrossRef]  

20. A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging (IEEE Press, 2001).