Laser Doppler holographic microscopy in transmission: application to fish embryo imaging

Transcription

1 Laser Doppler holographic microscopy in transmission: application to fish embryo imaging Nicolas Verrier, Daniel Alexandre, Michel Gross To cite this version: Nicolas Verrier, Daniel Alexandre, Michel Gross. Laser Doppler holographic microscopy in transmission: application to fish embryo imaging. Optics Express, Optical Society of America, 2014, 22, pp < /OE >. <hal > HAL Id: hal Submitted on 30 Apr 2014 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.

2 Laser Doppler holographic microscopy in transmission: application to fish embryo imaging Nicolas Verrier 1, Daniel Alexandre 1,2 and Michel Gross 1, 1 Laboratoire Charles Coulomb - UMR 5221 CNRS-UM2 CC 026 Université Montpellier II Place Eugène Bataillon Montpellier cedex, France 2 Montpellier Rio Imaging, Centre de Recherche en Biochimie Macromoléculaire - UMR 5237 CNRS-UM1-UM route de Mende Montpellier Cedex 5, France michel.gross@univ-montp2.fr Abstract: We have extended Laser Doppler holographic microscopy to transmission geometry. The technique is validated with living fish embryos imaged by a modified upright bio-microcope. By varying the frequency of the holographic reference beam, and the combination of frames used to calculate the hologram, multimodal imaging has been performed. Doppler images of the blood vessels for different Doppler shifts, images where the flow direction is coded in RGB colors or movies showing blood cells individual motion have been obtained as well. The ability to select the Fourier space zone that is used to calculate the signal, makes the method quantitative Optical Society of America OCIS codes: ( ) Digital holography; ( ) Holographic interferometry; ( ) Laser Doppler velocimetry; ( ) Interference microscopy; ( ) Scattering, particles; ( ) Spectroscopy, heterodyne References and links 1. E. Friedman, S. Krupsky, A. Lane, S. Oak, E. Friedman, K. Egan, and E. Gragoudas, Ocular blood flow velocity in age-related macular degeneration, Ophthalmology 102, (1995). 2. M. P. Pase, N. A. Grima, C. K. Stough, A. Scholey, and A. Pipingas, Cardiovascular disease risk and cerebral blood flow velocity, Stroke 43, (2012). 3. J. Kur, E. A. Newman, and T. Chan-Ling, Cellular and physiological mechanisms underlying blood flow regulation in the retina and choroid in health and disease, Prog. Retin. Eye Res. 31, (2012). 4. O. Sakurada, C. Kennedy, J. Jehle, J. Brown, G. L. Carbin, and L. Sokoloff, Measurement of local cerebral blood flow with iodo [14c] antipyrine, Am. J. Physiol.-Heart C. 234, H59 H66 (1978). 5. I. Kanno, H. Iida, S. Miura, M. Murakami, K. Takahashi, H. Sasaki, A. Inugami, F. Shishido, and K. Uemura, A system for cerebral blood flow measurement using an h215o autoradiographic method and positron emission tomography, J Cerebr. Blood F. Met. 7, (1987). 6. Y. Yeh and H. Cummins, Localized fluid flow measurements with an he-ne laser spectrometer, Appl. Phys. Lett. 4, (1964). 7. J. D. Briers, Laser doppler, speckle and related techniques for blood perfusion mapping and imaging, Physiol. Meas. 22, R35 R66 (2001). 8. J. D. Briers and S. Webster, Laser speckle contrast analysis (lasca): a nonscanning, full-field technique for monitoring capillary blood flow, J. Biomed. Opt. 1, (1996). 9. A. K. Dunn, Laser speckle contrast imaging of cerebral blood flow, Ann. Biomed. Eng. 40, (2012). 10. D. Briers, D. D. Duncan, E. Hirst, S. J. Kirkpatrick, M. Larsson, W. Steenbergen, T. Stromberg, and O. B. Thompson, Laser speckle contrast imaging: theoretical and practical limitations, J. Biomed. Opt. 18, (2013). 11. S. Yuan, A. Devor, D. A. Boas, and A. K. Dunn, Determination of optimal exposure time for imaging of blood flow changes with laser speckle contrast imaging, Appl. Opt. 44, (2005).

3 12. Y. Zeng, M. Wang, G. Feng, X. Liang, and G. Yang, Laser speckle imaging based on intensity fluctuation modulation, Opt. Lett. 38, (2013). 13. D. Gabor, A new microscopic principle, Nature 161, (1948). 14. D. Gabor, Microscopy by reconstructed wave-fronts, Proc. R. Soc. A 197, (1949). 15. E. N. Leith and J. Upatnieks, Reconstructed wavefronts and communication theory, J. Opt. Soc. Am. A 52, (1962). 16. U. Schnars and W. Jüptner, Direct recording of holograms by a ccd target and numerical reconstruction, Appl. Opt. 33, (1994). 17. A. Lozano, J. Kostas, and J. Soria, Use of holography in particle image velocimetry measurements of a swirling flow, Exp. Fluids 27, (1999). 18. Y. Pu and H. Meng, Four-dimensional dynamic flow measurement by holographic particle image velocimetry, Appl. Opt. 44, (2005). 19. J.-M. Desse, P. Picart, and P. Tankam, Digital three-color holographic interferometry for flow analysis, Opt. Express 16, (2008). 20. N. Verrier, S. Coëtmellec, M. Brunel, and D. Lebrun, Digital in-line holography in thick optical systems: application to visualization in pipes, Appl. Opt. 47, (2008). 21. N. Verrier, C. Remacha, M. Brunel, D. Lebrun, and S. Coëtmellec, Micropipe flow visualization using digital in-line holographic microscopy, Opt. Express 18, (2010). 22. F. Charrière, N. Pavillon, T. Colomb, C. Depeursinge, T. J. Heger, E. A. Mitchell, P. Marquet, and B. Rappaz, Living specimen tomography by digital holographic microscopy: morphometry of testate amoeba, Opt. Express 14, (2006). 23. W. Xu, M. Jericho, I. Meinertzhagen, and H. Kreuzer, Digital in-line holography for biological applications, P. Natl. Acad. Sci. USA 98, (2001). 24. B. Kemper and G. von Bally, Digital holographic microscopy for live cell applications and technical inspection, Appl. Opt. 47, A52 A61 (2008). 25. M. K. Kim, Principles and techniques of digital holographic microscopy, SPIE Reviews 1, (2010). 26. K. Lee, K. Kim, J. Jung, J. Heo, S. Cho, S. Lee, G. Chang, Y. Jo, H. Park, and Y. Park, Quantitative phase imaging techniques for the study of cell pathophysiology: from principles to applications, Sensors 13, (2013). 27. P. Picart, J. Leval, D. Mounier, and S. Gougeon, Some opportunities for vibration analysis with time averaging in digital fresnel holography, Appl. Opt. 44, (2005). 28. J. Leval, P. Picart, J. P. Boileau, and J. C. Pascal, Full-field vibrometry with digital fresnel holography, Appl. Opt. 44, (2005). 29. A. Asundi, V. R. Singh et al., Time-averaged in-line digital holographic interferometry for vibration analysis, Appl. Opt. 45, (2006). 30. F. Le Clerc, L. Collot, and M. Gross, Numerical heterodyne holography with two-dimensional photodetector arrays, Opt. Lett. 25, (2000). 31. M. Atlan, M. Gross, and E. Absil, Accurate phase-shifting digital interferometry, Opt. Lett. 32, (2007). 32. F. Verpillat, F. Joud, M. Atlan, and M. Gross, Digital holography at shot noise level, J. Display Technol. 6, (2010). 33. M. Lesaffre, N. Verrier, and M. Gross, Noise and signal scaling factors in digital holography in weak illumination: relationship with shot noise, Appl. Opt. 52, A81 A91 (2013). 34. M. Atlan and M. Gross, Laser doppler imaging, revisited, Rev. of Sci. Instrum. 77, (2006). 35. M. Atlan, M. Gross, and J. Leng, Laser doppler imaging of microflow, J. Europ. Opt. Soc. Rapid pub. 1, (2006). 36. F. Verpillat, F. Joud, M. Atlan, and M. Gross, Imaging velocities of a vibrating object by stroboscopic sideband holography, Opt. Express 20, (2012). 37. N. Verrier and M. Atlan, Absolute measurement of small-amplitude vibrations by time-averaged heterodyne holography with a dual local oscillator, Opt. Lett. 38, (2013). 38. N. Verrier, M. Gross, and M. Atlan, Phase-resolved heterodyne holographic vibrometry with a strobe local oscillator, Opt. Lett. 38, (2013). 39. M. Atlan, M. Gross, B. C. Forget, T. Vitalis, A. Rancillac, and A. K. Dunn, Frequency-domain wide-field laser doppler in vivo imaging, Opt. Lett. 31, (2006). 40. M. Atlan, B. C. Forget, A. C. Boccara, T. Vitalis, A. Rancillac, A. K. Dunn, and M. Gross, Cortical blood flow assessment with frequency-domain laser doppler microscopy, J. Biomed. Opt. 12, (2007). 41. M. Simonutti, M. Paques, J.-A. Sahel, M. Gross, B. Samson, C. Magnain, and M. Atlan, Holographic laser doppler ophthalmoscopy, Opt. Lett. 35, (2010). 42. C. Fang-Yen, S. Oh, Y. Park, W. Choi, S. Song, H. S. Seung, R. R. Dasari, and M. S. Feld, Imaging voltagedependent cell motions with heterodyne mach-zehnder phase microscopy, Opt. Lett. 32, (2007). 43. H. Iwai, T. Yamauchi, M. Miwa, and Y. Yamashita, Doppler-spectrally encoded imaging of translational ob-

4 jects, Opt. Commmun. 319, (2014). 44. S. Joseph, J.-M. Gineste, M. Whelan, and D. Newport, A heterodyne mach-zehnder interferometer employing static and dynamic phase demodulation techniques for live-cell imaging, Proc. SPIE 7554, 75540P 1 9 (2010). 45. M. Westerfield, The Zebrafish Book: a Guide for the Laboratory use of Zebrafish (Danio rerio) (Institute of Neuroscience. University of Oregon, 1995). 46. S. Isogai, M. Horiguchi, and B. M. Weinstein, The vascular anatomy of the developing zebrafish: an atlas of embryonic and early larval development, Dev. Biol. 230, (2001). 47. M. Atlan, M. Gross, B. C. Forget, T. Vitalis, A. Rancillac, and A. K. Dunn, Frequency-domain wide-field laser doppler in vivo imaging, Opt. Lett. 31, (2006). 48. M. Atlan, B. C. Forget, A. C. Boccara, T. Vitalis, A. Rancillac, A. K. Dunn, and M. Gross, Cortical blood flow assessment with frequency-domain laser doppler microscopy, J. Biomed. Opt. 12, (2007). 49. M. Atlan, M. Gross, T. Vitalis, A. Rancillac, J. Rossier, and A. Boccara, High-speed wave-mixing laser doppler imaging in vivo, Opt. Lett. 33, (2008). 50. N. Warnasooriya, F. Joud, P. Bun, G. Tessier, M. Coppey-Moisan, P. Desbiolles, M. Atlan, M. Abboud, and M. Gross, Imaging gold nanoparticles in living cell environments using heterodyne digital holographic microscopy. Opt. Express 18, (2010). 51. F. Verpillat, F. Joud, P. Desbiolles, and M. Gross, Dark-field digital holographic microscopy for 3d-tracking of gold nanoparticles, Opt. Express 19, (2011). 52. E. Cuche, P. Marquet, and C. Depeursinge, Spatial filtering for zero-order and twin-image elimination in digital off-axis holography, Appl. Opt. 39, (2000). 53. J. Gao, J. A. Lyon, D. P. Szeto, and J. Chen, In vivo imaging and quantitative analysis of zebrafish embryos by digital holographic microscopy, Biomed. Opt. Express 3, (2012). 1. Introduction Blood flow imaging techniques are widely used in biomedical studies, since they can assess physiological processes or can be used for early detection of disease [1 3]. However, many blood flow studies require, for imaging purposes, the use of a contrast agent, making the blood flow characterization invasive [4, 5]. Scanning Doppler imaging techniques can be considered to alleviate this issue, but due to the scanning step, acquisition of an image is a time consuming process [6]. An overview of the main techniques has been proposed by J. D. Briers [7]. Among these, one can outline two ways of monitoring the blood flow: performing measurements either in spatial domain (speckle analysis) or in temporal domain (Doppler analysis). Blood flow monitoring has been demonstrated by Laser Speckle Contrast Analysis/Imaging (LASCA/LSCI) [8]. Here, spatial statistics of the dynamic speckle are used to obtain blood flow images [9, 10]. Improvement of the acquired constrast image have been achieved through exposure time optimization [11], or intensity fluctuation analysis [12] resulting in high quality perfusion images. However, these techniques are only limited to perfusion monitoring and do not allow to assess other quantities such as amplitude or phase contrasts. Originally proposed by Gabor [13] as an optimization of electron microscopy, classical holography aims at recording, on a photo-sensitive material, the interference between a reference field and a diffracted object field. However, within the proposed common-path interferometric configuration, holographic imaging suffer from the so-called twin-image noise [14]. In 1962, Leith and Upatnieks propose to introduce an off-axis reference field, which lead to a separation of the object and its twin image contribution in spatial frequencies domain [15]. Development and democratization of high resolution CCD/CMOS sensors paid a major role in holographic imaging spreading out, making possible both digital recording and processing of holographic images [16]. Therefore, due to its intrinsic properties, digital holographic imaging is used in a wide variety of studies in fluid mechanics [17 21], biomedical imaging and microscopy [22 26], or mechanical inspection [27 29]. Further development led to the introduction of heterodyne digital holography [30]. Here, the reference field is dynamically phase shifted with respect to the object field. Therefore, the recorded hologram is time modulated, thus enabling phase-shifted interferometric measurements [31]. Moreover, this data acquisition scheme as been demonstrated to be shot-noise

5 limited [32, 33]. Temporal modulation feature of digital heterodyne holography can also be used to investigate dynamic phenomena, and being considered as a laser Doppler imaging technique [34]. The ability of heterodyne digital holography to perform Doppler imaging has been demonstrated in various domains such as microfluidics [35], vibration motion characterization [36 38], in vivo vasculature assessment, without contrast agent [39 41] and motion of biological objects [42 44]. In this paper, we combine laser Doppler holography [34] and transmission microscopy to analyze blood flow in fish embryos. We have adapted a laser Doppler holographic setup to a standard bio-microscope by carrying the two beams of the holographic interferometer (illumination of the object and reference), whose frequency offset is controlled, by optical fibers. Multimodal acquisition and analysis of the data is made by acting on the frequency offset of the two beams, and on the location of the Fourier space filtered zone. With the same set of data, we have imaged, with amplitude contrast, the whole fish embryo, or only the moving blood vessels. In that last case, we have obtained images where the flow direction is coded in RGB color. For large vessels, the dependance of the holographic signal with the frequency offset yields the Doppler frequency shift, and thus the flow velocity. For small vessels, individual Red Blood Cells (RBC) can be imaged, and movies showing the RBC motion are obtained. 2. Materials 2.1. Heterodyne holographic microscopy setup Fig. 1. Heterodyne digital holographic microscopy experimental arrangement (a) Injection part of the interferometer. HWP: half wave plate; PBS: polarizing beam splitter; AOM1, AOM2: acousto optic modulators (Bragg cells). (b) Classical upright microscope used for the off-axis recombinaison of reference and object beams. BS: cube beam splitter; CCD: CCD camera. Our off-axis digital holography set-up is a Mach-Zehnder interferometer in which the recombining cube beam splitter is angularly tilted. To fit our set-up into regular upright microscope, the interferometer is split into an injection (a) and a recombination (b) part, which are connected together with optical fibers (Fig. 1). The injection part (Fig. 1(a)) is mounted on a separate breadboard, and is used to build both object illumination (field E I ) and reference (or local oscillator (LO) E LO ) beams. Laser light is provided by a 785 nm, 80 mw Sanyo DL S continuous mono-mode laser diode (angular frequency ω I ). Laser beam is split by a polarizing beam splitter (PBS), which in combination with the first half wave plate (HWP) allows to adjust the optical power in both reference and object arms of the interferometer. To enable phase shifting interferometry and

6 frequency scanning of the holographic detection frequency, we control the frequency ω LO of the reference beam by using the heterodyne holography method [30]. The reference arm is thus frequency shifted by using two acousto-optic modulators AOM1 and AOM2 (80 MHz AA electronics ), operating at ω 1 and ω 2 respectively, so that ω LO = ω I + ω with ω = ω 2 ω 1. Reference and object fields are finally injected into two mono-mode optical fibers (Thorlabs P1-780A-FC-2 fiber patch). Both reference (E LO ) and object illumination (E I ) fibers are brought to an upright microscope (Olympus CX41) working in transmission configuration (Fig. 1(b)). The object is imaged by an Olympus DPlan microscope objective (MO: NA= 0.25, G=10, corrected at 160 mm). The illumination field (E I ) is scattered by the studied living object (the zebrafish embryo) yielding the object signal field (E). Signal (E) and reference field (E LO ) are combined together by a cube beam splitter (BS) that is angularly tilted thus introducing an off-axis angle in the interferometric arrangement. Interferences (i.e. E+ e j ωt E LO ) between both fields are recorded on a pixel (6.45 µm square pitch) 12-bits CCD camera (Jenoptik ProgRes MF ) operating at a framerate of ω S /(2π) 10 Hz. Recorded data are cropped to and FFT (Fast Fourier Transform) are made on a calculation grid. The CCD to MO distance is such that the effective MO enlargement factor is G = 10.4 (instead of G=10). In order to perfectly control the phase and frequency of the E LO reference beam, the CCD camera is triggered at frequency ω CCD, and the ω 1, ω 2 and ω CCD signals are generated by digital synthesizers driven by a common 10 MHz clock. Zebrafish (Danio rerio, wild type AB line) were maintained according to standard protocols [45]. Larvae were mounted on a standard glass slide under a 22 mm 1.5 round coverslip with a 0.5 mm thick caoutchouc ring as a spacer and sealant. They were embeded in a drop of cooling 1.5% low melting point agarose at 37 C and oriented before gelation. The chamber was filled with 100 mg/l tricaine in filtered tank water and imaging was performed at room temperature. Blood vessels nomenclature is according to [46]. 3. Laser Doppler holography mode Fig. 2. Detection of zebrafish embryo blood flow in transmission geometry. (a) Illustration of the spectral broadening of light encountering moving scattering objects (here, red blood cells). (b) Principles of the spectral measurement by heterodyne digital holography. Due to motions within the sample (mainly the RBC motion within the embryo vasculature), the incident field E I (of wave vector k I ) that is scattered to wave vector k S (yielding E S ) undergoes a Doppler shift ω S ω I = q.v (1) where v is the velocity of the scatterer and q=k S k I the scattering vector (see Fig. 2(a)) and q.v the inner product of vectors q and v. As detection is made for all wave vectors k S within

7 the MO collection cone, the frequency spectrum of the moving scatterers signal is broadened. One then gets two components: a sharp narrow peak from the photons that are not scattered (ballistic photons) or scattered by objects that do not move, and a broad pedestal from the photons that are scattered by moving objects and that are collected by MO. The width of the pedestal is proportional to the MO numerical aperture NA since q NA k I. With NA=0.25, the Doppler broadening we observe here is much smaller (10 to 100 Hz) than the one (100 to 1000 Hz) observed in previous in vivo Doppler holographic experiments made in reflection configuration [41, 47 49], where q 2 k I. Since our camera operates at ω S /(2π)=10 Hz, the bandwidth of our holographic detection (±5 Hz around ω LO /(2π)) in not enough to acquire the moving scatterer Doppler signal. As in previous works, the frequency ω LO = ω I + ω (with ω = ω 1 ω 2 ) is swept in order to acquire the whole Doppler signal and to analyze its spectrum. 4. Experimental results Fig. 3. Spatial filtering and reconstruction principle. Holograms H(x,y) (a), H 1 (k x,k y ) (b), H 2 (k x,k y ) (c) and H 3 (x,y,z = 0) (d): see Media 1. The display is made in arbitrary Log scale for the average intensity H X 2. This average is calculated made by making the reconstruction with I 0..I 3, with I 1..I 4... and with I 28..I 31, and by averaging the resulting H X 2. Images (a) and (d) show a side view of a 6 days old zebrafish embryo. Ventral side in on the right and anterior is up. The holographic information is acquired by sweeping the frequency shift ω/(2π) from 0 up to 100 Hz by step of ω CCD /2π (i.e. ω = mω CCD with m integer) except in section 4.1 where ω = ω CCD /4. On the other hand, the movies Media 1 and Media 2 have been recorded with ω = 0. To optimize sensitivity and selectivity, the illumination beam versus local oscillator ratio is adjusted in order to have E LO 2 > E 2, with E 2 as large as possible. Exposure time and frame rate are T exp = 50 ms and ω CCD /(2π)=1/T CCD = 10 Hz. For each frequency shift ω/(2π), sequences of N= 32 (or N= 96) CCD camera frames (i.e. I 0, I 1...I N 1 ) are recorded.

8 4.1. Spatial filtering and reconstruction. The reconstruction procedure is similar to the one used in previous works [50, 51]. To illustrate reconstruction, let us consider 4 phase detection of the ballistic peak signal (no Doppler shift). We have thus adjusted ω 1 and ω 2 to have ω = ω CCD /4, and calculated the demodulated hologram by: H(x,y) = [I 0 (x,y) I 2 (x,y)] + j[i 1 (x,y) I 3 (x,y)] (2) Figure 3(a) shows the holographic signal H(x,y) 2 displayed in log scale. To select the +1 holographic grating order by the Cuche et al. method [52], we have performed the holographic reconstruction of the MO pupil by the Schnars et al. method [16] yielding H 1 (k x,k y ) with: [ ] H 1 (k x,k y ) = FFT H(x,y) e j k (x2 +y 2 )/2d (3) where FFT is the Fast Fourier transform operator, and e j k (x2 +y 2 )/2d the phase factor that makes the pupil image sharp (see Fig. 3(b)). For a plane wave reference, the parameter d is equal to the MO pupil to CCD camera distance. Optimal spatial filtering is then made by cropping the MO pupil disk on its sharp circular edge, and by moving the cropped zone in the center of the Fourier space (see Fig. 3(c)) yielding H 2 (k x,k y ). The reconstructed hologram H 3 (x,y,z) is then calculated by the angular spectrum method that involves 2 FFT from FFT 1 [H 2 (k x,k y )]. We have thus: H 3 (x,y,z) = FFT 1[ H 2 (k x,k y ) e j(k2 x+k 2 y)/2z ] (4) where e j(k2 x+k 2 y)/2z is the phase factor that describe the propagation of the field from the CCD plane z = 0 to the location z of the zebrafish image (via MO). In most cases, the embryo is in focus on the camera, and z=0. Here, the selection of the +1 grating order is made by the crop, the compensation of the MO phase curvature by the e j k (x2 +y 2 )/2d phase factor, and the compensation of the off axis tilt by the translation of the cropped zone. We have calculated the so called instantaneous intensity signal H 3 2 by using data of successive sequences of 4 frames, i.e. by calculating H 3 2 with I 0..I 3, then with I 1..I 4... and to the end with I 28..I 31 ). The movie of the instantaneous intensity signal H 3 2 (see Media 1) shows that the time varying components of H 3 (t) 2 are extremely low. We have then averaged H 3 (t) 2 over the sequence of n frames, getting H 3 2 that is displayed on Fig. 3(d). The images show a side view of the zebrafish embryo. Due to the use of coherent light illumination, the image exhibits some speckle noise but anatomical details such as the notochord, somite boundaries and pigment cells are distinguished. In the corresponding movie Media 1, motion of blood can be barely seen in the caudal artery and vein but none is observe in the smaller capillaries Holographic Doppler images of the moving scatterers. To select the moving scatterer signal that corresponds to the Doppler pedestal, we have swept the detection frequency shift ω by step ω CCD, and calculated 2 phase holograms from the recorded data. We have thus: ω = ω LO ω I = mω CCD (5) H(x,y) = I 0 (x,y) I 1 (x,y) where m is integer. Note that H(x,y) is defined by Eq. 2 for 4 phase holograms, and by Eq. 5 for 2 phase holograms. We have kept the same notation (i.e. H(x, y)), because the reconstruction equations Eq. 3 and Eq. 4 are the same in both cases.

9 Fig. 4. Two phase detection efficiency η 2 as a function of the LO versus signal frequency offset: x. Calculation is made by Eq.6 with T CCD = 100 ms and T exp = 50 ms. We have calculated the 2 phase holographic detection efficiency (i.e. η for the field and η 2 for the intensity) as a function of the frequency ω S of the scattered field. This calculation is similar to the one done in previous works [31, 32]. We get: 1 1 ktccd η(x) = 2T exp ( 1) k +T exp /2 e j2πxt dt (6) k=0 t=kt CCD T exp /2 = 1 2 sinc(πxt exp) [ 1 e j2πxt ] CCD where T CCD = 2π/ω CCD is the frame time interval, and x=(ω S ω LO )/2π the LO versus signal frequency offset. The efficiency η 2 is plotted on Fig.4. As seen, the 2 phase detection selects scattered field components with x close to zero (i.e. with ω S ω LO ). Moreover, since η 2 is null for x = 20, 10,0,10,20... Hz, the holographic detection efficiency η 2 remains null at the illumination frequency ω I, whatever the frequency shift ω = mω CCD is. The ballistic peak (whose frequency is ω I ), which is much bigger than the Doppler pedestal, is thus not seen during the sweep. Figure 5 shows the Doppler average reconstructed holograms H 3 (x,y) 2 obtained with frequency shifts ω/(2π) varying from 0 Hz (a) to 80 Hz (i) by step of 10 Hz. The recorded hologram H(x,y), and the reconstructed intensity hologram H 3 (x,y) 2 are calculated from successive sequence of 2 frames (I n and I n+1 ) within the sequence of N= 32 frames (with N = 0 to 30), the displayed signal H 3 (x,y) 2 corresponding to the average of H 3 (x,y) 2 over n. Since the ballistic peak is not seen, the contrast of the Fig. 5 images is reversed with respect to Fig. 3 (d), and vessels where blood flows are therefore seen in white on a black background. In image (a) caudal vein (A) and caudal artery (B,C) are seen as well as intersegmental vessels (D, E) and dorsal longitudinal anastomotic vessels. An abnormal shunt between caudal artery and vein in seen near arrow (C). The reconstructed images (a) to (i) are displayed with the same Log scale for the average intensity H 3 2. One thus clearly see the decrease of the holographic signal with the frequency shift (ω LO ω I )/(2π) that is expected. Figure 6 analyzes this decreases for the zones of the reconstructed image marked by arrows (A) to (F) in Fig. 5(a). All curves are calculated by applying on the holographic signal H 3 2 a gaussian low pass filter (of 10 pixels width for curves A,B,C and F curves, and 4 pixels width for curves D and E), and by summing the filtered holographic signal over 8 8 pixels. On curve (A), which corresponds to the flow in caudal vein, the holographic signal decreases regularly with the frequency offset ω. The same

10 Fig. 5. Reconstructed hologram H 3 (x,y,z=0) made for (ω LO ω I )/(2π)=0 Hz (a): see also Media 2, 10 Hz (b)... to 80 Hz (i). Images (a) to (i) are displayed with the same Log scale for the average intensity H 3 2. The sample is the same as the one in Fig. 3. Fig. 6. Dependance of the Doppler holographic signal H 3 (x,y) 2 with the frequency offset (ω LO ω I )/(2π) for different location A to F of the reconstructed image of Fig. 3(a). Curves are drawn with linear (a) and logarithmic (b) scales.

11 occurs in (B) and (C), but the decrease is slower showing faster velocities for the blood cells in the caudal artery. In intersegmental vessels (curves D and E), the signal H 3 2 is much lower than the one of larger artery and vein. The corresponding frequency spectrum (Fig. 6D and E) exhibit a decrease of the signal at ω/(2π)= 40 and 80 Hz offset followed by an increase at 60 and 100 Hz. This might be related to the instrumental response of our holographic detection [31]. A fine analysis shows that the signal fluctuates a lot within each sequence of N frames, and from one sequence to the next, since the flowing blood cells can be visualized individually on the instantaneous intensity signal H 3 (t) 2 reconstructed with successive sequence of 2 frames. Movie Media 2 shows H 3 (t) 2 for the sequence of N= 32 frames recorded at 0 Hz (i.e. without frequency offset) that has been used to display H 3 2 on Fig. 5(a). Movie Media 3 shows H 3 (t) 2 for a sequence of N = 96 frames that image the embryo shown in Fig. 8(h) and (i). The RBCs individual motion in intersegmental vessels is clearly seen on Media 2 and Media 3. These two movies are very similar to those seen in reference [53]. They can be used to measure the RBC velocity by Particule Image Velocimetry (PIV) RGB holographic Doppler images yielding direction of motion. Fig. 7. Fourier space reconstructed hologram H 2 (k x,k y ) made without (a) and with (b,c) selection of the scattered wave vector k S. In (b) the selected zone is oriented toward u ϕ=0. In (c), three zones oriented along u ϕ with ϕ = 0 (blue), 2π/3 (green) and 4π/3 (red) are selected. Display is made in arbitrary log scale for H 2 2. Frequency shift is(ω LO ω I )= 0. The holographic experiment made in transmission configuration gives here a Doppler signal corresponding to single scattering events like in Dynamic Light Scattering (DLS) experiments. The Doppler shift is thus simply q.v where q = k S k I. The illumination (or incident) wave vector k I is fixed by the geometry of the experimental setup. In transmission geometry, k I is parallel to u z (where u x, u y and u z are the unit vectors in the x,y and z directions). On the other hand, k S =(k x,k y,k z ) with k S = 2π/λ is measured by our holographic experiment. In the reconstruction procedure, the cropped and translated zone seen on Fig. 3(c) and on Fig 7(a) is a map of the holographic signal as a function of k x and k y. It is then possible to select k S by spatial filtering in the Fourier space, and to reconstruct images corresponding to given intervals of k S. To illustrate this idea, we have selected three k S zones covering half of the Fourier space. These zones are oriented toward u ϕ=0, u ϕ=2π/3 and u ϕ=4π/3 with u ϕ = cosϕ u x + sinϕ u y. With the holographic signal H 2 (k x,k y ) of each zone, we have reconstructed a Red, Green or Blue colored image H 3 (x,y,z) 2, and we have combined these three images to get RGB images sensitive to the direction of k S (and thus to the direction of q and v). Figure 7 (b) shows the H 2 (k x,k y ) 2 signal, filtered by the ϕ = 0 spatial filter, which corresponds to the blue image. Figure 7 (c) shows the Fourier space RGB image of the H 2 (k x,k y ) 2 signal obtained by combining the three colored images corresponding to the

12 ϕ = 0, ϕ = 2π/3 and ϕ = 4π/3 filters. One must notice here, that the wave vector of ballistic light remains parallel to u z yielding a bright peak near the center of the Fourier space (k x,k y 0). Although most of the ballistic light is removed in 2 phase holograms recorded without frequency offset, this peak still remains visible in that case as seen on Fig. 7(a). To fully eliminate the ballistic light, the center of the Fourier space (i.e. k x,k y 0) has thus been removed from the selected zones as shown on Figure 7(b) and (c). Fig. 8. Colored reconstructed hologram H 3 (x,y,z=0) made for (ω LO ω I )/(2π)=0Hz (a), 10 Hz (b)... to 60 Hz (g), and for -20 Hz (h) and + 20 Hz (i). Images (a) to (g) are made with same zebrafish sample and same viewpoint as Fig. 3 and Fig. 5. Images (h) to (i) correspond to another fish embryo. Images (a) to (g) (and (h) and (i)) are displayed with the same color RGB Log scale for the average intensity H 3 2. Media 3 corresponds to images (h) and (i). Figure 8(a) to (g) shows the RGB colored reconstructed image of the zebrafish sample obtained with the holographic data of Fig. 5(a) to (g). Without frequency offset (Fig. 8(a) with ω = 0), the RGB colored reconstructed image is very similar to the black and white one (Fig. 5(a)), but with lower resolution, since half of the Fourier space information is used in the reconstruction. The images of Fig. 8 are obtained from a sequence of 32 images (i.e. with an acquisition time of 3.2 s). The reconstruction is made on GPU in less than one second for the 32 images. Since the detection shift ω is null, the detection efficiency does not depend on the sign of the Doppler shift of the signal ω S ω I. Detection is not sensitive to the sign of q.v, nor to the v direction (even if q direction is selected in the Fourier space). The colors in the reconstructed image of Fig. 8(a) are thus the same whatever the motion direction is. With increasing

13 frequency shifts ω, signal decreases but colors change, depending on the flow direction. On figure 8(d), with ω = 30 Hz, caudal artery appears in cyan, while caudal vein, where blood flows in opposite direction, is seen in orange. The flow direction corresponds here to the Fig. 8 color wheel. This wheel corresponds to the colored image of the Fourier space filtered signal H 2 (k x,k y ) 2 of Fig. 7(d) rotated by 180. To verify that the RGB colored images give a realistic information on the direction of the velocity v, we have recorded, with another embryo, holograms for both sign of the frequency offset ω. Figure 8(h) and (i) shows the reconstructed RGB images obtained for ω/(2π) = ±20 Hz. Here again, the colors of the reconstructed image depend on the flow direction, but the color coding is reversed for negative frequency offsets: ω < 0. The correspondence of the motion direction with the color wheel is only valid for positive shifts ω > Conclusion In this paper, we have shown that an upright commercial bio-microscope can be transformed into a powerful holographic setup, by coupling the microscope to an heterodyne holographic [30] breadboard with optical fibers. By controlling the frequency offset ω of the reference beam, and by using different combinations of frames to calculate the hologram H, we have imaged a living zebrafish embryo under different modalities. With ω = 0 and H =(I 0 I 2 )+ j(i 1 I 3 ), we have calculated the amplitude and the phase of the transmitted light. Although the phase can be qualitatively explored, we have not displayed phase images because the phase variation is too fast (the thickness of the embryo is hundreds of microns). With ω = mω CCD (where m 0 is integer) and H =(I 0 I 1 ), we have selected the signal whose Doppler shift ω S ω I is close to ω according to Eq. 4 and Fig. 4) and imaged the blood flow. With ω = 0 and H =(I 0 I 1 ), we have visualized individual moving bloods cells in intersegmental vessels, getting movies that show their instantaneous motion. Since the scattering is moderate like in DLS (Dynamic Light Scattering) experiment, the Doppler shift is directly related to the scatterers (i.e. blood cells) velocity v by ω S ω I = q.v with q=k S k I. Since k I u z is known, and since k S can be selected in the Fourier space, the scatterers velocity v can be measured quantitatively by a proper choice of both the frequency offset ω and the Fourier space selected zone. To illustrate one of the numerous possibilities of Doppler holography in transmission, we have imaged the blood cell motion direction in RGB colors. In the present paper, we have not used the possibilities for numerical refocusing to image different z layers of our sample. This could be useful, since all blood vessels are not in the same z plane. For instance, in Fig. 5 and Fig. 8(a)-8(g), the caudal artery is roughly on focus, while the caudal vein is out of focus by about 100 µm. Refocusing could thus be used to reconstruct the embryo vascular system in 3D. We intend to further improve the method by using a faster camera (to better analyze the motion of individual blood cells) and a higher numerical aperture objective (for better optical sectioning). Acknowledgments We acknowledge OSEO-ISI Datadiag grant and ANR Blanc Simi 10 (n 11 BS ) grant for funding, and Rémy Vialla for technical support.