The coherence-controlled holographic microscope (CCHM)1,2 is an innovative system particularly designed for quantitative phase imaging (QPI) and measurement of live cell dynamics.3 An achromatic off-axis interferometer based on the diffraction grating is used in CCHM. The spatial and temporal coherence of illumination can be widely varied, and in this way, imaging properties of the microscope are substantially modifiable. The coherence gating effect induced by low coherence makes possible formation of optical sections of the sample4 (in reflection mode) or imaging through turbid media5,6 (in transmission mode). Low coherence also improves lateral resolution and the imaging in general.7 For exploiting these effects, CCHM works with a broad polychromatic light source, which makes all the difference to other off-axis holographic microscopes usually equipped with a laser light source.8,9 This is because the high-coherence light source leads to the formation of unwanted artifacts in QPI as a consequence of coherence noise, random interferences, and diffraction of light. Low-coherence illumination, however, requires precise alignment and high stability of the system to be maintained during long-lasting time-lapse QPI studies of activity and reaction of living cells. And just these advantages stand for the contribution of CCHM to cell biology research. For alignment of highly sensitive interferometric systems, a secondary light source is often used. The interferometer state detection is carried out by an auxiliary detector10 or by a system of detectors.1112.–13 These methods are unsuitable for CCHM, because they would cause photo-toxic effect on living cells,10,11,13 and because the alignment of the many additional optical elements would be very complicated.12 For this reason, an original method based on the measurement of the modulus of the reconstructed holographic signal was elaborated for assessing the instrument state and guiding the optimization. The basic realignment three-dimensional algorithm (BReTA) method is suitable for CCHM because it does not require additional optical components and light sources. Moreover, it is amenable to full automation because the value of the reconstructed holographic signal is available from the online image processing and the operation can be robotized. Finally, the verification of the BReTA method applicability is presented.
Methods and Materials
Optical Setup and Image Processing
The optical assembly of CCHM is shown in Fig. 1. A halogen lamp is used as a light source; its light is guided by an optical fiber to the plane light source. This plane is imaged by the relay lens RL into the rear focal plane of condensers and so that the sample SP and the reference object RO are Köhler illuminated. The light beam is split by beam splitter and mirrors , , and to the optical paths of object and reference arm. Changeable aperture stop AS is placed into the beam for setting spatial coherence of the light and bandpass filter F for setting its temporal coherence. The specimen and the reference object are imaged from the object plane to the output plane OP by the objective lens and , tube lens and , and output lens and , respectively. The hologram is captured by the CCD camera D in the output plane OP. Diffraction grating DG and mirror are placed in the intermediate image planes behind the tube lenses and to form an interference pattern with the same spatial frequency of fringes (by the same carrier frequency ) for all wavelengths of light and for all points of the extended source LS in OP.2 In principle, the system provides off-axis holography by the interference of the zeroth and first diffraction orders. Achromaticity and spatial invariance of the interference pattern follow from the fact that the diffraction grating is imaged onto output plane OP by .14
The microscope is capable of imaging at various depths in a sample. Two approaches are available depending on the coherence of illumination.15,16 Full numerical refocusing is possible with highly coherent illumination, while with partially coherent illumination, the axial range of the refocusing gets reduced.15,16 On the other hand, low spatial coherence leads to the formation of an optical section in a broader sense,15 the thickness of which corresponds to the depth of field of the objective.17 Then, refocusing outside the section is only possible optically. This imaging mode is suitable for observation of adhering cells on a coverslip.
The hologram captured by the CCD is processed by the computer in the following way, see Fig. 2. First, the Fourier transform of the hologram is computed. Then, the spatial-frequency spectrum of the object image is obtained by the windowing operation on the sideband of the spatial-frequency spectrum of the hologram (around the carrier frequency ). The zeroth spatial frequency is shifted to the center of the window and the spectrum is multiplied by the Hanning weight function. Finally, the image complex amplitude is computed by the inverse Fourier transform, which provides the image phase and the image amplitude . A quantitative phase image is obtained after applying an unwrapping procedure and background compensation.18,19
Spatially incoherent illumination requires correct lateral alignment of the images in both arms of the microscope. Temporally incoherent light demands adjusting the arms to the same length. In order to prove that the modulus of the reconstructed holographic signal is of significant value for the microscope alignment, we had to describe its dependence on the optical path length difference, , of the arms and on the mutual displacement of the images in the output plane OP that is described by the corresponding displacement vector in the object plane. The theoretical reconstructed holographic signal can be described by the following formula [see Eq. (3.20) in Ref. 15 for zero defocus]:15], is a two-dimensional (2-D) coherent transfer function of CCHM [see Eq. (3.21) in Ref. 15], is the transverse part of scattering vector, and is the Cartesian coordinate vector in the output plane divided by the microscope magnification. For sample-free object space, the relation can be applied. Assuming a broad monochromatic light source and waves that are propagated at small angles to the optical axis, the modified equation can be obtained from Eq. (1) [see Eq. (6.1) in Ref. 15]. In our calculations, we retain the negative second power of the wave number for its subsequent extension to the broadband source, so that 20 with the support of the radius 1, is the 2-D effective pupil function of the illumination, where is the lowest of the numerical apertures of the condenser lenses and the objective lens in the reference arm, the function describes the distribution of the light intensity in the plane of the light source LS for the Köhler illumination,15 is the transverse part of , , and , is the wavelength and is the refraction index in the object space of the objective lenses and . Because , the function can be removed from the integrand in Eq. (2). If the image fields of both arms are mutually shifted by the nonzero displacement vector , the function must be added to Eq. (2) [see Eq. (2.25) in Ref. 15]. Then, 15)
Equation (4) describes the reconstructed holographic signal for a spatially incoherent source characterized by the intensity distribution . To also take into account the limited temporal coherence, suppose that the spectral properties of the source are described by the function , and the difference of the optical length of the arms in the image space is . The phase shift of the beam inclined by is then for the refractive index . Assuming that the sine condition holds for the microscope objectives, we get for the beam of the maximum inclination in the image space, where is the magnification between the object plane and the output plane OP. The maximum ratio for the objectives used in CCHM has been found for NIKON CFI S Fluor , where . Supposing and maximum , the difference of the phase shifts for the central beam and the maximum-inclination beam is . Hence, we can approximate the phase shift by for any inclination of the beam. After completing Eq. (4) by the spectral function of the source and by the complex exponential depending on and by the integration over , we get
For verification of the BReTA methodology, the modulus of the theoretically reconstructed holographic signal described by Eq. (5) was compared with experimental data. No specimen was inserted on the object plane. The measurement was performed with a spatially broad and spectrally narrowband light source. The source spectral function was given by the manufacturer data (Thorlabs) of the interference filter FB650-10 (, ); its radial intensity distribution was approximated by a Gaussian distribution21 with the proportional reciprocal standard deviation that was found by fitting the theoretical curve to the experimental data. To eliminate the noise, the modulus of the measured reconstructed holographic signal was averaged over the whole image field as follows:
Figure 3(a) shows the comparison of dependences of reconstructed holographic signal moduli and on . Measured curves slightly differ mutually and from the theoretical curve. It is due to the optical aberrations of the output lens , which is located behind the DG in the reference arm. Light diffracted to the first order passes through the border of aperture (in comparison to the axial propagation in the object arm). Moreover, shifting the objective lens laterally leads to transversal displacement of the diffracted beam in the aperture . In contrast, the measured values of modulus of the reconstructed holographic signal on Fig. 3(b) perfectly fit the theoretical values. Measurement was performed for the same light source and for by changing .
The theoretical reconstructed holographic signal modulus independence between and is shown in Fig. 3(c). Because of its symmetry, it is displayed only in the first quadrant. Its global maximum is apparent in the origin. Local maximum along the axis is the consequence of the form of the source spectral function . Behavior in the direction of the axis is smooth with no side lobes as a result of the Gaussian form of the intensity . The modulus has the maximum values on the axis for constant and on the axis for constant , see also experimental data in Fig. 4.
To demonstrate that the measured reconstructed holographic signal modulus, , has a strong maximum at and even with broadband light source, a measurement was carried out similar to the previous case, but with the narrowband filter removed. Figure 5 compares measurements containing the maximum value for the case of spectrally narrowband and broadband light source. Broadband illumination leads to a broader peak around [Fig. 5(b)] in comparison with spectrally narrowband illumination [Fig. 5(a)]. Extension of the peak is caused by the superposition of peaks related to different wavelengths, which do not overlap ideally due to chromatic aberration of the optical imaging system. Therefore, the measured reconstructed holographic signal modulus is a significant enough value for BReTA method in the case of spectrally narrowband light, as well as broadband light.
Basic Realignment Three-Dimensional Algorithm Method
Setup adjustment to basic realignment three-dimensional algorithm method
Several essential elements of the microscope setup had to be motorized (see Fig. 1) for experimental implementation of the BReTA method. The XY linear stage of the objective was motorized to vary ; the holder of the mirror was mounted on a linear stage for longitudinal movement . A sighting pattern was placed in the plane conjugated with the image plane, and computer-controlled shutters and were placed in the reference and object arms. The algorithm of the BReTA method described in this part is based on the measurement of the reconstructed holographic signal modulus (hereafter the signal ), see Fig. 6.
This procedure is used for correction of large misalignment typical for the microscope startup. First, the proper shifts of the mirror and objective lens are found with the aim to approach the area of the global maximum of . The final value of the signal must be higher than the threshold value , where is the value of the first side maximum and is the safety factor. In the beginning of the process, the sighting pattern P is inserted into the field plane and imaged to the object plane (the plane of the specimen SP and of the reference object RO) and finally to the output plane OP. The centers of images are mutually shifted by
After this process, the value of the signal is tested. If is greater than , the process ends. If not, is changed in the direction of increasing signal value by moving the mirror , until the condition is met.
Advanced alignment procedure
In this procedure, the elements and are moved sequentially to search for the maximum value of signal . Because the area of the global maximum is smooth, it is possible to choose any sequence of independent processes for search of the value , which are described in this section.
The first process deals with searching the greatest value of signal by changing the microscope objective position while maintaining constant. It can be performed in many ways. For example, by a 2-D scanning of objective lens around its current position; the resulting scans can be seen in Fig. 7(b). This is a very easy and robust way. Another faster possibility is to use a heuristic algorithm with appropriate termination condition. Algorithm of this process is illustrated in Fig. 7(a). It determines the value of signal in the initial position of . Its position is then changed by a small defined step in any direction, and the obtained value of signal is compared with . If is true, algorithm executes the next step right toward the original step direction; otherwise, the step is performed in the original direction. If the exit condition is not true, the value of signal is obtained and again compared with the previous value of the signal . The termination condition can be selected according to the situation, whichever occurs first. Optionally, the combination of the conditions may be required. The termination can be based on multiple detection of the greatest signal in the same position, when the algorithm is repeatedly going through the same coordinates, or on the calculation of the variational coefficient from the last several values of signal . The variation coefficient must be close to one.
The second process finds the value of the signal by moving the mirror in the direction of the increased values of the signal [see Fig. 7(b)]. It seeks such a position of mirror that corresponds to .
Both processes together lead to the alignment defined by and .
Testing of both procedures
For testing the procedures, Nikon objective lenses were used. Broadband light from the source was filtered by the interference filter (, ). Holograms were captured by camera XIMEA MR4021MC-BH. The choice of the minimum alignment steps in the lateral () and axial () directions depends on the targeted signal level, according to the graph in Fig. 3(c). As it will be shown in Sec. 2.4.5, the sufficient signal level resulting in a good-quality QPI is . The corresponding acceptable misalignment is in and in [see Fig. 3(c)]. Hence, for the initial alignment, we chose the step of , while is set approximately to zero by overlapping the sighting pattern images formed in reference and object arm. The robust scanning method runs laterally within positions with the steps of and , while the axial step of could be fined down to . As the movement of the mirror is linear and unidirectional, this correction did not noticeably prolong the axial scanning, in contrast to the situation in lateral directions.
The testing method consists of repeated misalignment of the microscope and subsequent activation of both alignment procedures. The microscope was randomly misaligned by the shift of the objective lens and by moving the mirror by a distance . Twenty independent measurements were performed. The average time required to approach the aligned state of the microscope was 136 s. Figure 8 shows the results of automated alignment in comparison with the manually obtained (i.e., the most precise) result. The relative values of signal , where is the value of the signal obtained by manual alignment, are on the axis. The axis is the count axis.
It is obvious that tested procedures of the BReTA method always resulted without fail with a high success rate.
Alignment procedure for a long-term experiment
During the long-term experiments, the value of signal is influenced by temperature changes or vibrations. Therefore, we extended the BReTA method for long-term maintenance of the maximum value of the signal . The procedure periodically activates both the transverse stages for the movement of the microscope lens and the axial stage to move the mirror . They are moved symmetrically around the current positions of active elements and by given steps. The procedure compares the original value of signal and values in adjacent positions of an active element. The active element is always moved to the position with the maximum value of signal . In Fig. 9(a), the time dependence of the signal values during the run of the procedure is shown. The procedure changes the positions of the microscope objective with the step 100 nm of and and the optical path difference of arms (the position of the mirror ) with the step 200 nm. As the procedure is running during observations, the minimum repeatable steps of the linear stages have been chosen to avoid visible changes of the image recorded. Figure 9(b) shows the time dependence of the signal without the alignment procedure running. The signal is inconsistent and too low most of the time. The decrease of the signal makes the phase noise higher, thereby impairing the QPI resolution. We define the QPI resolution as , where is the standard deviation of the background phase values measured in the window drawn in Fig. 9(c). The decrease of the signal to 90% has no measurable effect on the resolution; for its lower values, the resolution is elevated measurably [see Fig. 9(b)]. The increase of the noise in the phase image of real object is demonstrated in Figs. 9(d) and 9(e).
The microscope was subjected to temperature changes of environment with fluctuations of tenths of centigrade. It is obvious that the procedure successfully maintained the aligned state of the microscope.
We have developed a unique set of procedures constituting the original BReTA. The BReTA method allows for automated alignment of CCHM based on maximizing the value of modulus of the measured holographic signal. For exerting control over this parameter, some alignment elements of the original CCHM setup had to be motorized. The method consists in the initial alignment of the microscope in order to find the required minimal interference signal. The holographic signal is then optimized by searching the best alignment corresponding to the maximum signal. This maximum can be subsequently maintained by small changes in the alignment during a long-term experiment. All procedures were programmed in LabView and C++, and they are being used in the multimodal holographic microscopes produced by TESCAN ORSAY HOLDING a.s. The automatic BReTA method offers easy alignment of the microscope for ordinary users. However, it is essential for managing long-term QPI observations/measurements of live cells activity, which is the primary assignment of holographic microscopy. The BReTA method described in this article is patent pending.22
The authors thank their colleagues from the Experimental Biophotonics Group (CEITEC) for helpful discussions, especially Jana Collakova, Pavel Vesely, and Vera Kollarova. This work was supported by CEITEC-Central European Institute of Technology (CZ.1.05/1.1.00/02.0068) from the European Regional Development Fund.
T. Slaby et al., “Off-axis setup taking full advantage of incoherent illumination in coherence-controlled holographic microscope,” Opt. Express 21, 14747–14762 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.014747Google Scholar
M. Mir et al., “Progress in optics,” in Quantitative Phase Imaging, and E. Wolf, Ed., pp. 133–217, Elsevier, Amsterdam (2012).Google Scholar
L. Lovicar, J. Komrska and R. Chmelík, “Quantitative-phase-contrast imaging of a two-level surface described as a 2D linear filtering process,” Opt. Express 18, 20585–20594 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.020585Google Scholar
N. T. Shaked et al., “Quantitative analysis of biological cells using digital holographic microscopy,”in Holography, Research and Technologies, and J. Rosen, Ed., pp. 219–236, InTech, Rijeka (2011).Google Scholar
H. Iwai et al., “Quantitative phase imaging using actively stabilized phase-shifting low-coherence interferometry,” Opt. Lett. 29, 2399–2401 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.002399Google Scholar
R. Schwiesow, “Stabilization mechanism for optical interferometer,” US Patent 4, 444, 501 (1984).Google Scholar
V. V. Krishnamachari et al., “An active interferometer-stabilization scheme with linear phase control,” Opt. Express 14, 5210–5215 (2006).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.14.005210Google Scholar
R. Chmelik et al., “Progress in optics,” in The Role of Coherence in Image Formation in Holographic Microscopy, and E. Wolf, Ed., pp. 133–217, Elsevier, Amsterdam (2014).Google Scholar
M. Tyc et al., “Numerical refocusing in digital holographic microscopy with extended-sources illumination,” Opt. Express 21, 28258–28271 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.028258Google Scholar
J. Collakova et al., “Coherence-controlled holographic microscopy enabled recognition of necrosis as the mechanism of cancer cells death after exposure to cytopathic turbid emulsion,” J. Biomed. Opt. 20(11), 111213 (2015).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.20.11.111213Google Scholar
R. M. Goldstein, H. A. Zebken and C. L. Werner, “Satellite radar interferometry: two-dimensional phase unwrapping,” Radio Sci. 23(4), 713–720 (1988).RASCAD0048-6604http://dx.doi.org/10.1029/RS023i004p00713Google Scholar
D. C. Ghiglia and M. D. Pritt, Two-Dimensional Phase Unwrapping: Theory, Algorithms and Software, Wiley-Interscience, New York (1998).Google Scholar
J. W. Goodman, Introduction to Fourier Optics, 3rd ed., Roberts & Company, Greenwood Village (2005).Google Scholar
R. Chmelik et al., “Interferometric system with variable optics for incoherent source of radiation,” Czech Republic Patent Application No. PV 2014-538 (2014).Google Scholar
Zbynek Dostal is a PhD student at the Brno University of Technology. He received his BS degree in mechanical engineering in 2007 and MS degree in optics and precise mechanics in 2009 from the Brno University of Technology. He is a coauthor of the coherence-controlled holographic microscope patent for which he received the Werner von Siemens Excellence Award 2013. His current research interests include holographic microscopy, automation in microscopy, and optical and mechanical design.
Tomas Slaby is head of light microscopy R&D in TESCAN Brno company. He received his MS and PhD degrees in holographic microscopy from the Brno University of Technology in 2008 and 2015. He is involved in the development of a new generation of coherence-controlled holographic microscope.
Lukas Kvasnica is a PhD student at the Brno University of Technology. He received his MS degree in optics and precise mechanics in 2008 from the Brno University of Technology. Since 2012, he has been working as a C++ programmer in TESCAN Brno company, where he is involved in the development of the coherence-controlled holographic microscope.
Martin Lostak is a physicist in TESCAN Brno company. He received his MS and PhD degrees in holographic microscopy from the Brno University of Technology in 2008 and 2015. He is involved in developing the new generation of coherence-controlled holographic microscope.
Aneta Krizova is a PhD student at the Brno University of Technology. She received her BS and MS degrees in physical engineering from the Brno University of Technology in 2010 and 2012. Her current research interests include light microscopy, especially holographic microscopy and its applications.
Radim Chmelik is a professor of applied physics at the Brno University of Technology. He received his MS degree in solid-state physics from the Masaryk University in Brno in 1989 and his PhD degree in physical and materials engineering from the Brno University of Technology in 1997. He is the author of more than 40 journal papers. His current research interests include wave optics, imaging theory, advanced and three-dimensional light microscopy, and holographic microscopy.