1.IntroductionPhase distributions induced by objects cannot directly measure a camera because an image sensor can measure only intensity distributions. Phase-contrast1 and differential interference contrast2 techniques are proposed for visualizing the phase distributions qualitatively. To quantitatively measure the phase distribution, various quantitative phase imaging (QPI) techniques are proposed3 because quantitative phase distribution has informative physical parameters such as optical path difference,4 dry mass,5 aberrations,6 and shear stress.7 Numerous QPI techniques have been proposed, and they can be categorized as interferometric or non-interferometric methods. Although interferometric techniques accurately measure the phase distribution, special optical instruments are required to obtain interferograms such as a laser light source, phase-shifter, and interferometry.8 Contrarily, non-interferometric methods such as transport-of-intensity equation (TIE)9 and iterative phase retrieval algorithms10,11 require no special device because these methods obtain phase distributions using defocused intensity distributions. Particularly, TIE-based QPI can deterministically obtain phase distributions induced by objects. Due to the simplicity of the setup and the algorithm, TIE-based QPI has been applied to various research fields such as low-photon conditions,12,13 fluorescence imaging,14 and optical memory.15,16 To improve phase reconstruction accuracy, there are many phase reconstruction algorithms are proposed.17–19 Although TIE-based QPI needs multiple defocused intensity distributions to obtain high-quality phase distributions by shifting an image sensor along an optical axis, various configurations have been proposed to mitigate this problem.20–24 As an advantage of TIE-based QPI, the partially coherent light source can be used for phase imaging.25,26 Although using a partially coherent light source can reduce the speckle noise compared with using a laser, the ambiguity of the phase for wavelength increases when the temporal coherence is low.27 TIE-based QPI has been applied to microscopy so far.28,29 In the research field of biology, a confocal laser scanning microscope is usually used for visualizing fluorescence signals from biosamples in real time. To simultaneously measure fluorescence signals and phase distributions, multi-modal microscopy is required. Although the holographic approach measures phase and fluorescence distributions simultaneously,30 a specially designed optical setup is required.31 While a commercially available wide-field microscope would also be useful for measuring phase and fluorescence simultaneously, the axial resolution of the fluorescence image is degraded. There is an attempt to apply TIE to confocal microscopy;32 however, a special varifocal lens is needed to get multiple defocused intensity distributions. In this paper, we propose TIE-based QPI using commercially available confocal microscopy. Using the bright-field module of the commercially available confocal microscope, we have successfully demonstrated that phase retrieval can be achieved without additional optical elements and specially designed equipment. Because this setup is laser scanning microscopy, the light source is temporally coherent and spatially partially coherent,33,34 which means that detected intensity distributions by the proposed method are speckle-free and the ambiguity of the phase for wavelength is low. We demonstrate phase imaging using living cells of the moss Physcomitrium patens and HeLa cells. We also show experimental results of multi-modal fluorescence and phase imaging of living cells under confocal microscopy. 2.MethodsThe schematic of TIE-based QPI is shown in Fig. 1. The TIE is derived from the Helmholtz equation under the condition of paraxial approximation and is described as follows:9 where and denote the gradient operator in the lateral dimension and the wavenumber, respectively. and indicate the phase and intensity distributions of an object at . To solve Eq. (1) for , the calculation of an axial intensity derivative is required.In general, the axial intensity derivative is approximated with multiple defocused intensity distributions and is described as follows:35 where . is the coefficients of . shows a defocus interval. When a target is a pure phase object (phase-only object), the phase distribution can be deterministically obtained using the Fourier transform as follows:36 where and are Fourier and inverse Fourier transform operators. and are the 2D coordinates of spatial frequencies and a regularization parameter to prevent divergence. is an infocus intensity distribution.3.ExperimentsThe schematic of the optical setup of confocal microscopy for the proposed method is shown in Fig. 2. This is an epi-fluorescence confocal microscopy. The light from a laser is scanned at an object plane using Galvano mirrors, and the light from the object goes back to the same optical path and is detected by single-pixel detectors. Due to the pinhole in front of the detector, high axial resolution is achieved in fluorescence imaging. On the other hand, the light through the object is detected by a single-pixel detector without a pinhole. Due to the absence of the pinhole, the detected intensity distributions are similar to bright field images. Because laser scanning microscopy can be considered as a spatially partially coherent technique,33,34 speckle-free images are obtained. By changing the axial position of the object, multiple defocused intensity distributions are obtained. Then, these intensity distributions are used for the axial intensity derivative described in Eq. (2). The object phase distribution can be obtained using the intensity derivative through Eq. (3). Notably, in this study, the samples are mechanically shifted according to the optical axis to get defocus distributions. Although the proposed method can also be applied to real-time imaging, the measurable dynamics are limited by the shift speed. In this study, is set to 1 to reduce the TIE calculation time. The defocus distance in this study is empirically decided to suppress the noise effect. The regularization parameter is also empirically chosen. First, we evaluated the applicability of QPI under the off-the-shelf confocal microscopy using a micro-lens array (Thorlabs, Inc.: MLA300-14AR-M, Newton, United States). The commercially available confocal microscopy (Nikon, Co. Ltd.: A1 Rsi, Tokyo, Japan) was used. The wavelength of the laser was 488 nm. The scanning pitch was set to . The magnification of the objective lens was 10. The number of pixels was . The lenslet pitch and the radius curvature of the micro-lens array were and 14.2 mm, respectively. The defocused distance for the TIE was set to . The measured defocused intensity distributions are shown in Figs. 3(a)–3(c). The phase distribution obtained by the TIE is shown in Fig. 3(d). The micro-lens array was also measured using self-built off-axis digital holographic microscopy with a Mach-Zehnder interferometer. The result of digital holography is shown in Fig. 3(e). The phase distribution obtained by digital holography was unwrapped by the TIE-based phase unwrapping method37 and is shown in Fig. 3(f). The sectional profiles at red and blue lines in Figs. 3(d) and 3(e) are shown in Fig. 3(g). To compare the results of Figs. 3(d) and 3(e), bias phases were manually adjusted. The results of the proposed method and digital holography are almost consistent. The experimental results indicate the proposed method can accurately measure the phase distribution of the object. Second, we applied the proposed method to visualize the phase distribution of fluorescent protein-labeled living plant cells. We used the filamentous tissue called protonema of the transgenic moss P. patens. In the transgenic line, the yellow fluorescent protein gene Citrine was inserted into the histone H3.3 locus (Pp3c18_14481),38,39 allowing for detection of fluorescence in the nuclei.30 We prepared the samples for confocal microscopy as described previously.40 Confocal microscopy (Olympus Co., Ltd.: FV1000, Tokyo, Japan) visualized the plant cells. The wavelength of the laser for phase imaging was 635 nm. The wavelengths of the lasers for exciting the yellow fluorescent protein Citrine localized to the nuclei and chloroplast autofluorescence were 473 and 635 nm, respectively. The scanning pitch was set to . The images were captured at 0.15 frames per second (fps) and the 3D stack of 41 images took 19 min. The used objective lens was UPLSAPO20X (Olympus). The number of pixels was . The defocused distance for the TIE was set to . Fluorescent and phase images were simultaneously observed. Experimental results are shown in Fig. 4. The nuclei visualized with Citrine fluorescence in Fig. 4(a) are unidentifiable in the phase distribution of Fig. 4(d), suggesting no clear phase difference distribution between nuclei and cytoplasm. By contrast, chloroplast autofluorescence in Fig. 4(b) highly corresponded to the phase distribution of Fig. 4(d), indicating the large phase difference between chloroplasts and cytoplasm. Compared with chloroplasts in the bright-field image of Fig. 4(c), chloroplast visualized with the phase distribution in Fig. 4(d) shows in-focus chloroplast as Fig. 4(b). This suggests the shallow depth of focus in our method. In Fig. 4(c), the positions of the cross-wall (septa) but not the side-wall are indecisive due to the crowded chloroplasts. However, the phase distribution in Fig. 4(d) improved the visualization of the positions. Our method might facilitate label-free imaging of specific cell walls of plants. As shown in Fig. 4(c), this sample has absorption, which may affect the phase reconstruction accuracy due to the assumption made in the derivation of Eq. (3). The accuracy of phase measurement in TIE with this moss sample should be compared with that of digital holography in the future. Using a TIE algorithm that takes absorption into account may improve the accuracy of phase measurement.41,42 Third, we simultaneously observe fluorescence and phase distributions of living mammalian cultured cells, HeLa cells (RCB0007, Riken CellBank, Tsukuba, Japan). Cellular nuclei, mitochondria, and microtubules were stained by NucBlue Live ReadyProbes Reagent (Invitrogen), MitoTracker Red CMXRos, and Tubulin Tracker Green, respectively. Confocal microscopy was performed under the A1 Rsi (Nikon, Co. Ltd.) equipped with a water-immersion objective lens (PlanApo-VC 60xA WI) at room temperature. The transmission light images were obtained using a 488 nm laser and used for TIE-based QPI. The fluorescence images of the nucleus, mitochondria, and microtubules were simultaneously acquired using the excitation laser of 405, 488, and 561 nm. The fluorescent signals of for nuclei, for microtubules, and for mitochondria were detected. The scanning pitch was set to . The images were captured at 1 fps and the 3D stack of 81 images took 4.7 min. The number of pixels was . The defocused distance for the TIE was set to . The fluorescence images of the nucleus, microtubules, and mitochondria are shown in Figs. 5(a)–5(c), and phase distribution is shown in Fig. 5(d). In Fig. 5(d), a curved line-like phase difference was observed around the center of the cell. These phase differences coincide with the edge of the nucleus, nuclear membrane, in Fig. 5(a). Furthermore, circular phase differences in the nuclear region [Fig. 5(d)] correspond to nucleoli visualized as lacking fluorescence in the nucleus [Fig. 5(a)]. This result suggests that the nuclear membrane and nucleoli can be visualized without labeling by the proposed method. Although individual mitochondria clearly visualized in fluorescence image were indistinguishable from the phase distribution, the region with a higher density of mitochondria had relatively higher phase difference from other regions [Figs. 5(a) and 5(d)]. By contrast, microtubules visualized with fluorescent dyes [Fig. 5(b)] were not identifiable in Fig. 5(d). These suggest that no clear phase differences were detected between cytosol and microtubules by the proposed method. Note that the processing time for the TIE calculation is about 0.8 s with an Intel Core i7-11700 CPU at 2.50 GHz. 4.ConclusionIn this paper, we proposed TIE-based QPI using commercially available confocal microscopy. We evaluated the feasibility of the proposed method using the microlens array. The experimental result indicates that QPI can be realized in commercially available confocal microscopy using the TIE technique. In addition, we applied this technique to practical cases using plant and mammalian cultured cells. This method will be useful for measuring dry mass,5 viscosity, and temperature of cells and for correcting phase fluctuation to cancel aberration and scattering caused by an object.43 To obtain information such as dry mass or intracellular temperature distribution, the accuracy of this method, including its sensitivity and reliability, must be ensured. In future work, it will be necessary to compare and evaluate the phase detection sensitivity of the proposed method with other established quantitative phase measurement techniques. The proposed method could be further improved through the fusion of diffraction tomography. One limitation of the proposed method is that the general TIE-based analysis assumes samples with small scattering and absorption effects. Therefore, when the effects of scattering and absorption are dominant, the phase distribution obtained using this method is less accurate than the results obtained by general interferometric methods. This effect could be mitigated by applying a diffraction tomography algorithm29 or by combining a reconstruction method that takes absorption effects into account.41,42 An example of an application of confocal microscopy is the surface inspection of substrates. In this paper, we have measured biological samples, but it could be applied to the inspection of objects such as diffractive optical elements. Code and Data AvailabilityAll relevant data, materials, and software code used in this research are available upon request from the corresponding author. AcknowledgmentsThe imaging with a confocal microscope NikonA1 was performed at the Optics and Imaging Facility, Trans-Scale Biology Center, NIBB, and with an FV1000 system on an IX-81 inverted optical microscope was carried out at Center for Instrumental Analysis, Utsunomiya University. A part of this research is supported by JSPS KAKENHI (Grant Nos. 20H05886, 21H04663, 20H05891, and 20H05669) and AMED Brain/MINDS (Grant No. JP19dm0207078). ReferencesF. Zernike,
“How is discovered phase contrast,”
Science, 121
(3141), 345
–349 https://doi.org/10.1126/science.121.3141.345 SCIEAS 0036-8075
(1955).
Google Scholar
R. D. Allen, G. B. David and G. Nomarski,
“The Zeiss-Nomarski differential interference equipment for transmitted-light microscopy,”
Zeitschrift fur wissenschaftliche Mikroskopie und mikroskopische Technik, 69 193
–221
(1969).
Google Scholar
Y. Park, C. Depeursinge and G. Popescu,
“Quantitative phase imaging in biomedicine,”
Nat. Photonics, 12
(10), 578
–589 https://doi.org/10.1038/s41566-018-0253-x NPAHBY 1749-4885
(2018).
Google Scholar
V. H. F. Muñoz et al.,
“Measurement of mean thickness of transparent samples using simultaneous phase shifting interferometry with four interferograms,”
Appl. Opt., 55
(15), 4047
–4051 https://doi.org/10.1364/AO.55.004047 APOPAI 0003-6935
(2016).
Google Scholar
S. Aknoun et al.,
“Living cell dry mass measurement using quantitative phase imaging with quadriwave lateral shearing interferometry: an accuracy and sensitivity discussion,”
J. Biomed. Opt., 20
(12), 126009 https://doi.org/10.1117/1.JBO.20.12.126009 JBOPFO 1083-3668
(2015).
Google Scholar
I. Choi, K. Lee and Y. Park,
“Compensation of aberration in quantitative phase imaging using lateral shifting and spiral phase integration,”
Opt. Express, 25 30771
–30779 https://doi.org/10.1364/OE.25.030771 OPEXFF 1094-4087
(2017).
Google Scholar
W. J. Eldridge et al.,
“Imaging deformation of adherent cells due to shear stress using quantitative phase imaging,”
Opt. Lett., 41 352
–355 https://doi.org/10.1364/OL.41.000352 OPLEDP 0146-9592
(2016).
Google Scholar
I. Yamaguchi and T. Zhang,
“Phase-shifting digital holography,”
Opt. Lett., 22 1268
–1270 https://doi.org/10.1364/OL.22.001268 OPLEDP 0146-9592
(1997).
Google Scholar
M. R. Teague,
“Deterministic phase retrieval: a Green’s function solution,”
J. Opt. Soc. Amer., 73 1434
–1441 https://doi.org/10.1364/JOSA.73.001434
(1983).
Google Scholar
F. Zhang, G. Pedrini and W. Osten,
“Phase retrieval of arbitrary complex-valued fields through aperture-plane modulation,”
Phys. Rev. A, 75 043805 https://doi.org/10.1103/PhysRevA.75.043805
(2007).
Google Scholar
A. Grjasnow, A. Wutting and R. Riesenberg,
“Phase resolving microscopy by multi-plane diffraction detection,”
J. Microsc., 231
(1), 115
–123 https://doi.org/10.1111/j.1365-2818.2008.02022.x JMICAR 0022-2720
(2008).
Google Scholar
K. Komuro, Y. Yamazaki and T. Nomura,
“Transport-of-intensity computational ghost imaging,”
Appl. Opt., 57 4451
–4456 https://doi.org/10.1364/AO.57.004451 APOPAI 0003-6935
(2018).
Google Scholar
A. K. Gupta, N. K. Nishchal and P. P. Banerjee,
“Transport of intensity equation based photon-counting phase imaging,”
OSA Contin., 3 236
–245 https://doi.org/10.1364/OSAC.383527
(2020).
Google Scholar
S. K. Rajput et al.,
“Three-dimensional fluorescence imaging using the transport of intensity equation,”
J. Biomed. Opt., 25
(3), 032004 https://doi.org/10.1117/1.JBO.25.3.032004 JBOPFO 1083-3668
(2019).
Google Scholar
N. Yoneda et al.,
“Transport-of-intensity holographic data storage based on a computer-generated hologram,”
Appl. Opt., 57
(30), 8836
–8840 https://doi.org/10.1364/AO.57.008836 APOPAI 0003-6935
(2018).
Google Scholar
M. Bunsen and S. Tateyama,
“Detection method for the complex amplitude of a signal beam with intensity and phase modulation using the transport of intensity equation for holographic data storage,”
Opt. Express, 27 24029
–24042 https://doi.org/10.1364/OE.27.024029 OPEXFF 1094-4087
(2019).
Google Scholar
C. Zuo et al.,
“Transport-of-intensity phase imaging using Savitzky-Golay differentiation filter - theory and applications,”
Opt. Express, 21
(5), 5346
–5362 https://doi.org/10.1364/OE.21.005346 OPEXFF 1094-4087
(2013).
Google Scholar
Z. Jingshan et al.,
“Transport of intensity phase imaging by intensity spectrum fitting of exponentially spaced defocus planes,”
Opt. Express, 22
(9), 10661
–10674 https://doi.org/10.1364/OE.22.010661 OPEXFF 1094-4087
(2014).
Google Scholar
N. Yoneda et al.,
“Single-shot higher-order transport-of-intensity quantitative phase imaging using deep learning,”
Appl. Opt., 60 8802
–8808 https://doi.org/10.1364/AO.435538 APOPAI 0003-6935
(2021).
Google Scholar
C. Zuo et al.,
“High-speed transport-of-intensity phase microscopy with an electrically tunable lens,”
Opt. Express, 21
(20), 24060
–24075 https://doi.org/10.1364/OE.21.024060 OPEXFF 1094-4087
(2013).
Google Scholar
L. Waller et al.,
“Transport of intensity phase imaging in a volume holographic microscope,”
Opt. Lett., 35 2961
–2963 https://doi.org/10.1364/OL.35.002961 OPLEDP 0146-9592
(2010).
Google Scholar
L. Waller et al.,
“Phase from chromatic aberrations,”
Opt. Express, 18
(22), 22817
–22825 https://doi.org/10.1364/OE.18.022817 OPEXFF 1094-4087
(2010).
Google Scholar
C. Zuo et al.,
“Noninterferometric single-shot quantitative phase microscopy,”
Opt. Lett., 38
(18), 3538
–3541 https://doi.org/10.1364/OL.38.003538 OPLEDP 0146-9592
(2013).
Google Scholar
N. Yoneda et al.,
“Single-shot higher-order transport-of-intensity quantitative phase imaging based on computer-generated holography,”
Opt. Express, 29 4783
–4801 https://doi.org/10.1364/OE.415598 OPEXFF 1094-4087
(2021).
Google Scholar
D. Paganin and K. A. Nugent,
“Noninterferometric phase imaging with partially coherent light,”
Phys. Rev. Lett., 80 2586
–2589 https://doi.org/10.1103/PhysRevLett.80.2586 PRLTAO 0031-9007
(1998).
Google Scholar
C. Zuo et al.,
“High-resolution transport-of-intensity quantitative phase microscopy with annular illumination,”
Sci. Rep., 7 7654 https://doi.org/10.1038/s41598-017-06837-1 SRCEC3 2045-2322
(2017).
Google Scholar
A. M. Zysk et al.,
“Transport of intensity and spectrum for partially coherent fields,”
Opt. Lett., 35 2239
–2241 https://doi.org/10.1364/OL.35.002239 OPLEDP 0146-9592
(2010).
Google Scholar
Y. Li et al.,
“Quantitative phase microscopy for cellular dynamics based on transport of intensity equation,”
Opt. Express, 26
(1), 586
–593 https://doi.org/10.1364/OE.26.000586 OPEXFF 1094-4087
(2018).
Google Scholar
J. Li et al.,
“Transport of intensity diffraction tomography with non-interferometric synthetic aperture for three-dimensional label-free microscopy,”
Light: Sci. Appl., 11 154 https://doi.org/10.1038/s41377-022-00815-7
(2022).
Google Scholar
M. Kumar et al.,
“Common-path multimodal three-dimensional fluorescence and phase imaging system,”
J. Biomed. Opt., 25
(3), 032010 https://doi.org/10.1117/1.JBO.25.3.032010 JBOPFO 1083-3668
(2020).
Google Scholar
M. Schnell et al.,
“Quantitative confocal phase imaging by synthetic optical holography,”
Opt. Express, 22 15267
–15276 https://doi.org/10.1364/OE.22.015267 OPEXFF 1094-4087
(2014).
Google Scholar
J. Zheng et al.,
“Dual-mode phase and fluorescence imaging with a confocal laser scanning microscope,”
Opt. Lett., 43 5689
–5692 https://doi.org/10.1364/OL.43.005689 OPLEDP 0146-9592
(2018).
Google Scholar
K. Itoh and Y. Ohtsuka,
“Coherence control by laser scanning,”
Appl. Opt., 19 3184
–3188 https://doi.org/10.1364/AO.19.003184 APOPAI 0003-6935
(1980).
Google Scholar
Y. Ichioka et al.,
“Digital scanning laser microscope,”
Appl. Opt., 24 691
–696 https://doi.org/10.1364/AO.24.000691 APOPAI 0003-6935
(1985).
Google Scholar
L. Waller, L. Tian and G. Barbastathis,
“Transport of intensity phase-amplitude imaging with higher order intensity derivatives,”
Opt. Express, 18
(12), 12552
–12561 https://doi.org/10.1364/OE.18.012552 OPEXFF 1094-4087
(2010).
Google Scholar
L. Allen and M. Oxley,
“Phase retrieval from series of images obtained by defocus variation,”
Opt. Commun., 199
(1), 65
–75 https://doi.org/10.1016/S0030-4018(01)01556-5 OPCOB8 0030-4018
(2001).
Google Scholar
J. Martinez-Carranza, K. Falaggis and T. Kozacki,
“Fast and accurate phase-unwrapping algorithm based on the transport of intensity equation,”
Appl. Opt., 56
(25), 7079
–7088 https://doi.org/10.1364/AO.56.007079 APOPAI 0003-6935
(2017).
Google Scholar
A. A. Heikal et al.,
“Molecular spectroscopy and dynamics of intrinsically fluorescent proteins: coral red (dsRed) and yellow (Citrine),”
Proc. Natl. Acad. Sci., 97
(22), 11996
–12001 https://doi.org/10.1073/pnas.97.22.11996
(2000).
Google Scholar
D. Lang et al.,
“The physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution,”
Plant J., 93
(3), 515
–533 https://doi.org/10.1111/tpj.13801
(2018).
Google Scholar
T. Tomoi, Y. Coudert and T. Fujita, Tracking Intercellular Movement of Fluorescent Proteins in Bryophytes, 321
–332 Springer US, New York, NY
(2022). Google Scholar
M. Watabe et al.,
“Optical dispersions through intracellular inhomogeneities,”
Phys. Rev. Res., 5 L022043 https://doi.org/10.1103/PhysRevResearch.5.L022043 PRSTCR 1554-9178
(2023)).
Google Scholar
J. Zhang et al.,
“On a universal solution to the transport-of-intensity equation,”
Opt. Lett., 45 3649
–3652 https://doi.org/10.1364/OL.391823 OPLEDP 0146-9592
(2020).
Google Scholar
Y. Ashida et al.,
“Imaging performance of microscopy adaptive-optics system using scene-based wavefront sensing,”
J. Biomed. Opt., 25
(12), 123707 https://doi.org/10.1117/1.JBO.25.12.123707 JBOPFO 1083-3668
(2020).
Google Scholar
|
Phase distribution
Confocal microscopy
Phase imaging
Fluorescence
Microscopes
Phase measurement
Visualization