Special Section on Quantitative Phase Imaging in Biomedicine

Automated alignment method for coherence-controlled holographic microscope

[+] Author Affiliations
Zbynek Dostal, Radim Chmelik

Brno University of Technology, Institute of Physical Engineering, Faculty of Mechanical Engineering, Technicka 2896/2, Brno 61600, Czech Republic

Brno University of Technology, CEITEC-Central European Institute of Technology, Technicka 3058/10, Brno 61600, Czech Republic

Tomas Slaby, Martin Lostak

TESCAN ORSAY HOLDING a.s., Libusina trida 21, Brno 623 00, Czech Republic

Lukas Kvasnica, Aneta Krizova

Brno University of Technology, CEITEC-Central European Institute of Technology, Technicka 3058/10, Brno 61600, Czech Republic

TESCAN ORSAY HOLDING a.s., Libusina trida 21, Brno 623 00, Czech Republic

J. Biomed. Opt. 20(11), 111215 (Oct 28, 2015). doi:10.1117/1.JBO.20.11.111215
History: Received March 31, 2015; Accepted October 1, 2015
Text Size: A A A

Open Access Open Access

Abstract.  A coherence-controlled holographic microscope (CCHM) was developed particularly for quantitative phase imaging and measurement of live cell dynamics, which is the proper subject of digital holographic microscopy (DHM). CCHM in low-coherence mode extends DHM in the study of living cells. However, this advantage is compensated by sensitivity of the system to easily become misaligned, which is a serious hindrance to wanted performance. Therefore, it became clear that introduction of a self-correcting system is inevitable. Accordingly, we had to devise a theory of a suitable control and design an automated alignment system for CCHM. The modulus of the reconstructed holographic signal was identified as a significant variable for guiding the alignment procedures. From this, we derived the original basic realignment three-dimensional algorithm, which encompasses a unique set of procedures for automated alignment that contains processes for initial and advanced alignment as well as long-term maintenance of microscope tuning. All of these procedures were applied to a functioning microscope and the tested processes were successfully validated. Finally, in such a way, CCHM is enabled to substantially contribute to study of biology, particularly of cancer cells in vitro.

Figures in this Article

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.1113 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 w¯ 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.

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 C1 and C2 so that the sample SP and the reference object RO are Köhler illuminated. The light beam is split by beam splitter BS1 and mirrors M1, M4, and M5 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 O1 and O2, tube lens TL1 and TL2, and output lens OL1 and OL2, respectively. The hologram is captured by the CCD camera D in the output plane OP. Diffraction grating DG and mirror M2 are placed in the intermediate image planes behind the tube lenses TL1 and TL2 to form an interference pattern with the same spatial frequency of fringes (by the same carrier frequency fc) 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 OL2.14

Graphic Jump Location
Fig. 1
F1 :

Optical setup of the coherence-controlled holographic microscope (CCHM). LS, light source; AS, aperture stop; F, filter; RL, relay lens; M, mirrors; C, condensers; SP, specimen; RO, reference object; O, objectives; DG, diffraction grating; BS, beam splitters; OL, output lenses; OP, output plane; D, detector (CCD camera). The motorized parts of CCHM: mechanism for inserting a sighting pattern P, computer-controlled shutters S, two-dimensional transversal linear stage of objective O2, and linear stage with the axial movement direction of the mirror M2.

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 fc). 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 wD(i,j) is computed by the inverse Fourier transform, which provides the image phase Arg[wD(i,j)] and the image amplitude |wD(i,j)|. A quantitative phase image is obtained after applying an unwrapping procedure and background compensation.18,19

Graphic Jump Location
Fig. 2
F2 :

Image processing: (a) hologram of live cell, (b) the spatial-frequency spectrum with the windowing operation around the carrier frequency fc, (c) amplitude image, (d) phase image, and (e) quantitative phase image.

Theory

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, ΔL, of the arms and on the mutual displacement of the images in the output plane OP that is described by the corresponding displacement vector qf=(xf,yf) 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]: Display Formula

w(qt)=Tt(Qt)Ht(Qt)exp(2πiQt·qt)d2Qt,(1)
where Tt(Qt) is the Fourier transform of the object transmission function [see Eq. (3.19) in Ref. 15], Ht(Qt) is a two-dimensional (2-D) coherent transfer function of CCHM [see Eq. (3.21) in Ref. 15], Qt=(X,Y) is the transverse part of scattering vector, and qt=(x,y) is the Cartesian coordinate vector in the output plane divided by the microscope magnification. For sample-free object space, the relation Tt(Qt)=δ(X)δ(Y) 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 K for its subsequent extension to the broadband source, so that Display Formula
w(qt)=K2Ht(0,0)=K2St*(Kt;K)Pt(Kt;K)d2Kt,(2)
where Pt(Kt;K)=circ[(Ktn)/(KNA)] is the 2-D pupil function of the objective lens O1 with numerical aperture NA, circ(r) is a rotationally symmetrical function20 with the support of the radius 1, St(Kt;K)=it(Kt)circ[(Ktn)/(KNAs)] is the 2-D effective pupil function of the illumination, where NAs is the lowest of the numerical apertures of the condenser lenses and the objective lens O2 in the reference arm, the function it(Kt) describes the distribution of the light intensity in the plane of the light source LS for the Köhler illumination,15Kt is the transverse part of K, Kt=|Kt|, and |K|=K=1/λ, λ is the wavelength and n is the refraction index in the object space of the objective lenses O1 and O2. Because NAsNA, the function Pt(Kt;K) can be removed from the integrand in Eq. (2). If the image fields of both arms are mutually shifted by the nonzero displacement vector qf=(xf,yf), the function exp(2πiKt·qf) must be added to Eq. (2) [see Eq. (2.25) in Ref. 15]. Then, Display Formula
w(qt;qf)=K2it(Kt)circ(KtnKNAs)exp(2πiKt·qf)d2Kt.(3)
If we assume rotational symmetry of the intensity distribution it(Kt) of the light source LS, the Fourier transform [Eq. (3)] can be converted into the Hankel transform (see Ref. 15) Display Formula
w(qt;qf)=2π0NAs/niT(κ)J0(2πKκ·qf)κdκ,(4)
where iT(κ) is a radial intensity distribution, κ=Kt/K, κ=|κ|.

Equation (4) describes the reconstructed holographic signal for a spatially incoherent source characterized by the intensity distribution iT(κ). To also take into account the limited temporal coherence, suppose that the spectral properties of the source are described by the function iK(K), and the difference of the optical length of the arms in the image space is ΔL. The phase shift of the beam inclined by ψ is then ΔΦψ=2πKΔLcosψ for the refractive index n1. Assuming that the sine condition holds for the microscope objectives, we get sinψ=NAs/MNA/M for the beam of the maximum inclination ψ in the image space, where M is the magnification between the object plane and the output plane OP. The maximum ratio NA/M for the objectives used in CCHM has been found for NIKON CFI S Fluor 20×/0.75, where ψ=0.038rad. Supposing λ=650nm and maximum ΔL=100μm, the difference of the phase shifts for the central beam and the maximum-inclination beam is ΔΦ0ΔΦψ=2πKΔL(1cosψ)=0.21πrad. Hence, we can approximate the phase shift by ΔΦ0 for any inclination of the beam. After completing Eq. (4) by the spectral function iK(K) of the source and by the complex exponential depending on ΔΦ0 and by the integration over K, we get Display Formula

w(qt;qf,ΔL)=2π0iK(K)exp(2πiKΔL)0NAsiT(κ)J0(2πKκ·qf)κdκdK.(5)
Obviously, modulus w¯=|w| of the theoretical reconstructed holographic signal reaches its maximum for ΔL=0μm and for |qf|=0μm. For this reason, the modulus w¯ emerged as a significant enough variable for the BReTA method of CCHM alignment.

Experiment

For verification of the BReTA methodology, the modulus w¯ 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 iK(K) was given by the manufacturer data (Thorlabs) of the interference filter FB650-10 (λ=650nm, FWHM=10nm); its radial intensity distribution iT(κ) was approximated by a Gaussian distribution21 with the proportional reciprocal standard deviation α=2.5 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: Display Formula

w¯D=i=0NI1j=0NJ1|wD(i,j)|NINJ,(6)
where NI and NJ are sizes of wD(i,j).

Figure 3(a) shows the comparison of dependences of reconstructed holographic signal moduli w¯ and w¯D on |qf|. Measured curves slightly differ mutually and from the theoretical curve. It is due to the optical aberrations of the output lens OL2, which is located behind the DG in the reference arm. Light diffracted to the first order passes through the border of aperture OL2 (in comparison to the axial propagation in the object arm). Moreover, shifting the objective lens O2 laterally leads to transversal displacement of the diffracted beam in the aperture OL2. 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 |qf|=0μm by changing ΔL.

Graphic Jump Location
Fig. 3
F3 :

Comparison between measured values of moduli of theoretically reconstructed holographic signal w¯ (solid line) and measured reconstructed holographic signal w¯D (symbol lines). (a) Dependence on |qf| for identically long arms of the interferometer ΔL=0μm. Profiles from the qf=(0,0)μm to four directions are shown by the symbol lines. (b) Dependence on ΔL for |qf|=0μm. (c) Contour graph of the dependence of theoretically reconstructed holographic signal modulus w¯ on |qf| and ΔL. The values of w¯ and w¯D are normalized.

The theoretical reconstructed holographic signal modulus w¯ independence between |qf| and ΔL is shown in Fig. 3(c). Because of its symmetry, it is displayed only in the first quadrant. Its global maximum w¯max is apparent in the origin. Local maximum w¯l,i along the axis ΔL is the consequence of the form of the source spectral function iK(K). Behavior in the direction of the axis |qf| is smooth with no side lobes as a result of the Gaussian form of the intensity iT(κ). The modulus w¯ has the maximum values on the axis ΔL for constant |qf| and on the axis |qf| for constant ΔL, see also experimental data in Fig. 4.

Graphic Jump Location
Fig. 4
F4 :

Measured dependence of the modulus w¯D of reconstructed holographic signal on the displacement vector qf=(xf,yf) and on the optical path difference ΔL. Gray levels in the pixels (small squares) of each subimage represent w¯D average on the field of view for various qf and constant ΔL indicated in the subimage. Average values of w¯D are normalized by its maximum over the complete set of measurements.

To demonstrate that the measured reconstructed holographic signal modulus, w¯D, has a strong maximum at |qf|=0μm and ΔL=0μm 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 w¯D,max for the case of spectrally narrowband and broadband light source. Broadband illumination leads to a broader peak around w¯D,max [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 w¯D is a significant enough value for BReTA method in the case of spectrally narrowband light, as well as broadband light.

Graphic Jump Location
Fig. 5
F5 :

Measured dependence of the modulus w¯D of reconstructed holographic signal on the displacement vector qf=(xf,yf) (w¯D measured and represented as in Fig. 4, ΔL=0μm) for (a) spectrally narrowband light source (λ=650nm, FWHM=10nm) and (b) broadband light source. The values of w¯D are normalized; the normalized maximum value w¯D<1 for ΔL0μm.

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 O2 was motorized to vary qf; the holder of the mirror M2 was mounted on a linear stage for longitudinal movement ΔL. A sighting pattern was placed in the plane conjugated with the image plane, and computer-controlled shutters S1 and S2 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 w¯D (hereafter the signal w¯D), see Fig. 6.

Graphic Jump Location
Fig. 6
F6 :

The flow chart of the procedure of initial alignment.

Initial alignment

This procedure is used for correction of large misalignment typical for the microscope startup. First, the proper shifts of the mirror M2 and objective lens O2 are found with the aim to approach the area of the global maximum of w¯D. The final value of the signal must be higher than the threshold value w¯D=kw¯D,l1, where w¯D,l1 is the value of the first side maximum and k 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 Display Formula

qf=(xPrxPo,yPryPo).(7)
The coordinates of the centers of the object and reference arm images are found by closing the shutters S1 and S2, respectively, and measuring the position of the sighting pattern centers Po=[xPo,yPo], Pr=[xPr,yPr], respectively. The required shift qf is then expressed as follows: Display Formula
qf=qfM=(xPrxPo,yPryPo)M,(8)
where M is the magnification between the object plane and the output plane OP. Subsequently, the objective lens O2 is shifted by qf. This sequence is repeated until |qf|0μm. Then, the elements P, S1, and S2 are removed.

After this process, the value of the signal w¯D is tested. If w¯D is greater than w¯D, the process ends. If not, ΔL is changed in the direction of increasing signal value w¯D by moving the mirror M2, until the condition is met.

Advanced alignment procedure

In this procedure, the elements O2 and M2 are moved sequentially to search for the maximum value of signal w¯D,max. Because the area of the global maximum is smooth, it is possible to choose any sequence of independent processes for search of the value w¯D,max, which are described in this section.

The first process deals with searching the greatest value of signal w¯D by changing the microscope objective O2 position while maintaining ΔL constant. It can be performed in many ways. For example, by a 2-D scanning of objective lens O2 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 w¯D,1 in the initial position of O2. Its position qf is then changed by a small defined step in any direction, and the obtained value of signal w¯D,2 is compared with w¯D,1. If w¯D,2<w¯D,1 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 w¯D,3 is obtained and again compared with the previous value of the signal w¯D,2. 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 w¯D 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 w¯D. The variation coefficient must be close to one.

Graphic Jump Location
Fig. 7
F7 :

Illustration of the advanced alignment procedure. The modulus w¯D of holographic signal is measured, normalized, and represented as in Fig. 4. (a) First, the heuristic algorithm searches for the local maximum value of w¯D changing qf and keeping ΔL constant. This is done by lateral movements of the microscope objective O2. (b) Then, the global maximum value w¯D,max is found changing ΔL by the axial movement of the mirror M2 in the direction of increasing values of w¯D. The values of the reconstructed holographic signal modulus (i.e., values of the signal) w¯D are normalized.

The second process finds the value of the signal w¯D,max by moving the mirror M2 in the direction of the increased values of the signal w¯D [see Fig. 7(b)]. It seeks such a position of mirror M2 that corresponds to ΔL=0μm.

Both processes together lead to the alignment defined by ΔL=0μm and |qf|=0μm.

Testing of both procedures

For testing the procedures, Nikon objective lenses 10×/0.30 were used. Broadband light from the source was filtered by the interference filter (λ=650nm, FWHM=10nm). Holograms were captured by camera XIMEA MR4021MC-BH. The choice of the minimum alignment steps in the lateral (qf) and axial (ΔL) 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 90%. The corresponding acceptable misalignment is 0.5μm in qf and 10μm in ΔL [see Fig. 3(c)]. Hence, for the initial alignment, we chose the step 10μm of ΔL, while qf is set approximately to zero by overlapping the sighting pattern images formed in reference and object arm. The robust scanning method runs laterally within 5×5 positions with the steps 0.5μm of xf and yf, while the axial step of ΔL could be fined down to 1μm. As the movement of the mirror M2 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 O2 and by moving the mirror M2 by a distance 100μm. 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 w¯D,s=(w¯D,max/w¯D,max,m)·100%, where w¯D,max,m is the value of the signal obtained by manual alignment, are on the x axis. The y axis is the count axis.

Graphic Jump Location
Fig. 8
F8 :

Histogram demonstrating success rate of the tested procedures.

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 w¯D is influenced by temperature changes or vibrations. Therefore, we extended the BReTA method for long-term maintenance of the maximum value of the signal w¯D. The procedure periodically activates both the transverse stages for the movement of the microscope lens O2 and the axial stage to move the mirror M2. They are moved symmetrically around the current positions of active elements O2 and M2 by given steps. The procedure compares the original value of signal w¯D,max and values w¯D in adjacent positions of an active element. The active element is always moved to the position with the maximum value of signal w¯D. In Fig. 9(a), the time dependence of the signal values w¯D during the run of the procedure is shown. The procedure changes the positions of the microscope objective O2 with the step 100 nm of xf and yf and the optical path difference ΔL of arms (the position of the mirror M2) 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 w¯D 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 r=2σϕ, 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).

Graphic Jump Location
Fig. 9
F9 :

(a) Long-term maintenance of high values of the measured holographic signal w¯D and changes of qf=(xf,yf) and ΔL, which show the corresponding shifts of the objective lens O2 and of the mirror M2, respectively. The values of the signal w¯D are normalized. (b) The time dependence of the signal w¯D without the long-term maintenance procedure. Measured phase resolution r in selected points of the graph is indicated. (c) Quantitative phase image of rat sarcoma cells obtained with 90% value of the holographic signal [r=0.038rad, see the graph (b)]. The measurement of σϕ was performed within the window (22,000 pixels) plotted on the image. (d) Phase distribution along the line depicted in (c) for 100% value of the holographic signal with the long-term maintenance procedure activated. (e) Phase distribution along the same line with the apparent noise increase after the procedure deactivation and the holographic signal drop to 70% value.

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 w¯ 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.

Kolman  P., and Chmelik  R., “Coherence-controlled holographic microscope,” Opt. Express. 18, , 21990 –22003 (2010). 1094-4087 CrossRef
Slaby  T.  et al., “Off-axis setup taking full advantage of incoherent illumination in coherence-controlled holographic microscope,” Opt. Express. 21, , 14747 –14762 (2013). 1094-4087 CrossRef
Mir  M.  et al., “Progress in optics,” in Quantitative Phase Imaging. , and Wolf  E., Ed., pp. 133 –217,  Elsevier ,  Amsterdam  (2012).
Chmelik  R., and Harna  Z., “Parallel-mode confocal microscope,” Opt. Eng.. 38, (10 ), 1635 –1639 (1999).CrossRef
Lostak  M.  et al., “Diffuse light imaging with a coherence controlled holographic microscope,” Proc. SPIE. 7746, , 77461N  (2010). 0277-786X CrossRef
Lostak  M.  et al., “Coherence-controlled holographic microscopy in diffuse media,” Opt. Express. 22, , 4180 –4195 (2014). 1094-4087 CrossRef
Lovicar  L., , Komrska  J., and Chmelík  R., “Quantitative-phase-contrast imaging of a two-level surface described as a 2D linear filtering process,” Opt. Express. 18, , 20585 –20594 (2010). 1094-4087 CrossRef
Shaked  N. T.  et al., “Quantitative analysis of biological cells using digital holographic microscopy,”in Holography, Research and Technologies. , and Rosen  J., Ed., pp. 219 –236,  InTech ,  Rijeka  (2011).
Charrie  F.  et al., “Cell refractive index tomography by digital holographic microscopy,” Opt. Lett.. 31, , 178 –180 (2006). 0146-9592 CrossRef
Iwai  H.  et al., “Quantitative phase imaging using actively stabilized phase-shifting low-coherence interferometry,” Opt. Lett.. 29, , 2399 –2401 (2004). 0146-9592 CrossRef
Kalamatianos  D.  et al., “Active alignment for two-beam interferometers,” Rev. Sci. Instrum.. 77, (1 ), 013103  (2006). 0034-6748 CrossRef
Schwiesow  R., “Stabilization mechanism for optical interferometer,” US Patent 4, 444, 501 (1984).
Krishnamachari  V. V.  et al., “An active interferometer-stabilization scheme with linear phase control,” Opt. Express. 14, , 5210 –5215 (2006). 1094-4087 CrossRef
Leith  E. N., and Upatnikes  J., “Holography with achromatic-fringe systems,” J. Opt. Soc. Am.. 57, , 975 –979 (1967). 0030-3941 CrossRef
Chmelik  R.  et al., “Progress in optics,” in The Role of Coherence in Image Formation in Holographic Microscopy. , and Wolf  E., Ed., pp. 133 –217,  Elsevier ,  Amsterdam  (2014).
Tyc  M.  et al., “Numerical refocusing in digital holographic microscopy with extended-sources illumination,” Opt. Express. 21, , 28258 –28271 (2013). 1094-4087 CrossRef
Collakova  J.  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). 1083-3668 CrossRef
Goldstein  R. M., , Zebken  H. A., and Werner  C. L., “Satellite radar interferometry: two-dimensional phase unwrapping,” Radio Sci.. 23, (4 ), 713 –720 (1988). 0048-6604 CrossRef
Ghiglia  D. C., and Pritt  M. D., Two-Dimensional Phase Unwrapping: Theory, Algorithms and Software. ,  Wiley-Interscience ,  New York  (1998).
Goodman  J. W., Introduction to Fourier Optics. , 3rd ed.,  Roberts & Company ,  Greenwood Village  (2005).
Harris  F., “On the use of windows for harmonic analysis with the discrete Fourier transform,” Proc. IEEE. 66, , 51 –83 (1978). 0018-9219 CrossRef
Chmelik  R.  et al., “Interferometric system with variable optics for incoherent source of radiation,” Czech Republic Patent Application No. PV 2014-538 (2014).

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.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Zbynek Dostal ; Tomas Slaby ; Lukas Kvasnica ; Martin Lostak ; Aneta Krizova, et al.
"Automated alignment method for coherence-controlled holographic microscope", J. Biomed. Opt. 20(11), 111215 (Oct 28, 2015). ; http://dx.doi.org/10.1117/1.JBO.20.11.111215


Figures

Graphic Jump Location
Fig. 5
F5 :

Measured dependence of the modulus w¯D of reconstructed holographic signal on the displacement vector qf=(xf,yf) (w¯D measured and represented as in Fig. 4, ΔL=0μm) for (a) spectrally narrowband light source (λ=650nm, FWHM=10nm) and (b) broadband light source. The values of w¯D are normalized; the normalized maximum value w¯D<1 for ΔL0μm.

Graphic Jump Location
Fig. 4
F4 :

Measured dependence of the modulus w¯D of reconstructed holographic signal on the displacement vector qf=(xf,yf) and on the optical path difference ΔL. Gray levels in the pixels (small squares) of each subimage represent w¯D average on the field of view for various qf and constant ΔL indicated in the subimage. Average values of w¯D are normalized by its maximum over the complete set of measurements.

Graphic Jump Location
Fig. 3
F3 :

Comparison between measured values of moduli of theoretically reconstructed holographic signal w¯ (solid line) and measured reconstructed holographic signal w¯D (symbol lines). (a) Dependence on |qf| for identically long arms of the interferometer ΔL=0μm. Profiles from the qf=(0,0)μm to four directions are shown by the symbol lines. (b) Dependence on ΔL for |qf|=0μm. (c) Contour graph of the dependence of theoretically reconstructed holographic signal modulus w¯ on |qf| and ΔL. The values of w¯ and w¯D are normalized.

Graphic Jump Location
Fig. 2
F2 :

Image processing: (a) hologram of live cell, (b) the spatial-frequency spectrum with the windowing operation around the carrier frequency fc, (c) amplitude image, (d) phase image, and (e) quantitative phase image.

Graphic Jump Location
Fig. 1
F1 :

Optical setup of the coherence-controlled holographic microscope (CCHM). LS, light source; AS, aperture stop; F, filter; RL, relay lens; M, mirrors; C, condensers; SP, specimen; RO, reference object; O, objectives; DG, diffraction grating; BS, beam splitters; OL, output lenses; OP, output plane; D, detector (CCD camera). The motorized parts of CCHM: mechanism for inserting a sighting pattern P, computer-controlled shutters S, two-dimensional transversal linear stage of objective O2, and linear stage with the axial movement direction of the mirror M2.

Graphic Jump Location
Fig. 6
F6 :

The flow chart of the procedure of initial alignment.

Graphic Jump Location
Fig. 7
F7 :

Illustration of the advanced alignment procedure. The modulus w¯D of holographic signal is measured, normalized, and represented as in Fig. 4. (a) First, the heuristic algorithm searches for the local maximum value of w¯D changing qf and keeping ΔL constant. This is done by lateral movements of the microscope objective O2. (b) Then, the global maximum value w¯D,max is found changing ΔL by the axial movement of the mirror M2 in the direction of increasing values of w¯D. The values of the reconstructed holographic signal modulus (i.e., values of the signal) w¯D are normalized.

Graphic Jump Location
Fig. 8
F8 :

Histogram demonstrating success rate of the tested procedures.

Graphic Jump Location
Fig. 9
F9 :

(a) Long-term maintenance of high values of the measured holographic signal w¯D and changes of qf=(xf,yf) and ΔL, which show the corresponding shifts of the objective lens O2 and of the mirror M2, respectively. The values of the signal w¯D are normalized. (b) The time dependence of the signal w¯D without the long-term maintenance procedure. Measured phase resolution r in selected points of the graph is indicated. (c) Quantitative phase image of rat sarcoma cells obtained with 90% value of the holographic signal [r=0.038rad, see the graph (b)]. The measurement of σϕ was performed within the window (22,000 pixels) plotted on the image. (d) Phase distribution along the line depicted in (c) for 100% value of the holographic signal with the long-term maintenance procedure activated. (e) Phase distribution along the same line with the apparent noise increase after the procedure deactivation and the holographic signal drop to 70% value.

Tables

References

Kolman  P., and Chmelik  R., “Coherence-controlled holographic microscope,” Opt. Express. 18, , 21990 –22003 (2010). 1094-4087 CrossRef
Slaby  T.  et al., “Off-axis setup taking full advantage of incoherent illumination in coherence-controlled holographic microscope,” Opt. Express. 21, , 14747 –14762 (2013). 1094-4087 CrossRef
Mir  M.  et al., “Progress in optics,” in Quantitative Phase Imaging. , and Wolf  E., Ed., pp. 133 –217,  Elsevier ,  Amsterdam  (2012).
Chmelik  R., and Harna  Z., “Parallel-mode confocal microscope,” Opt. Eng.. 38, (10 ), 1635 –1639 (1999).CrossRef
Lostak  M.  et al., “Diffuse light imaging with a coherence controlled holographic microscope,” Proc. SPIE. 7746, , 77461N  (2010). 0277-786X CrossRef
Lostak  M.  et al., “Coherence-controlled holographic microscopy in diffuse media,” Opt. Express. 22, , 4180 –4195 (2014). 1094-4087 CrossRef
Lovicar  L., , Komrska  J., and Chmelík  R., “Quantitative-phase-contrast imaging of a two-level surface described as a 2D linear filtering process,” Opt. Express. 18, , 20585 –20594 (2010). 1094-4087 CrossRef
Shaked  N. T.  et al., “Quantitative analysis of biological cells using digital holographic microscopy,”in Holography, Research and Technologies. , and Rosen  J., Ed., pp. 219 –236,  InTech ,  Rijeka  (2011).
Charrie  F.  et al., “Cell refractive index tomography by digital holographic microscopy,” Opt. Lett.. 31, , 178 –180 (2006). 0146-9592 CrossRef
Iwai  H.  et al., “Quantitative phase imaging using actively stabilized phase-shifting low-coherence interferometry,” Opt. Lett.. 29, , 2399 –2401 (2004). 0146-9592 CrossRef
Kalamatianos  D.  et al., “Active alignment for two-beam interferometers,” Rev. Sci. Instrum.. 77, (1 ), 013103  (2006). 0034-6748 CrossRef
Schwiesow  R., “Stabilization mechanism for optical interferometer,” US Patent 4, 444, 501 (1984).
Krishnamachari  V. V.  et al., “An active interferometer-stabilization scheme with linear phase control,” Opt. Express. 14, , 5210 –5215 (2006). 1094-4087 CrossRef
Leith  E. N., and Upatnikes  J., “Holography with achromatic-fringe systems,” J. Opt. Soc. Am.. 57, , 975 –979 (1967). 0030-3941 CrossRef
Chmelik  R.  et al., “Progress in optics,” in The Role of Coherence in Image Formation in Holographic Microscopy. , and Wolf  E., Ed., pp. 133 –217,  Elsevier ,  Amsterdam  (2014).
Tyc  M.  et al., “Numerical refocusing in digital holographic microscopy with extended-sources illumination,” Opt. Express. 21, , 28258 –28271 (2013). 1094-4087 CrossRef
Collakova  J.  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). 1083-3668 CrossRef
Goldstein  R. M., , Zebken  H. A., and Werner  C. L., “Satellite radar interferometry: two-dimensional phase unwrapping,” Radio Sci.. 23, (4 ), 713 –720 (1988). 0048-6604 CrossRef
Ghiglia  D. C., and Pritt  M. D., Two-Dimensional Phase Unwrapping: Theory, Algorithms and Software. ,  Wiley-Interscience ,  New York  (1998).
Goodman  J. W., Introduction to Fourier Optics. , 3rd ed.,  Roberts & Company ,  Greenwood Village  (2005).
Harris  F., “On the use of windows for harmonic analysis with the discrete Fourier transform,” Proc. IEEE. 66, , 51 –83 (1978). 0018-9219 CrossRef
Chmelik  R.  et al., “Interferometric system with variable optics for incoherent source of radiation,” Czech Republic Patent Application No. PV 2014-538 (2014).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.