5 August 2015 Quantitative phase imaging through scattering media by means of coherence-controlled holographic microscope
Author Affiliations +
J. of Biomedical Optics, 20(11), 111206 (2015). doi:10.1117/1.JBO.20.11.111206
A coherence-controlled holographic microscope (CCHM) enables quantitative phase imaging with coherent as well as incoherent illumination. The low spatially coherent light induces a coherence gating effect, which makes observation of samples possible also through scattering media. The paper describes theoretically and simulates numerically imaging of a two-dimensional object through a static scattering layer by means of CCHM, with the main focus on the quantitative phase imaging quality. The authors have investigated both strongly and weakly scattering media characterized by different amounts of ballistic and diffuse light. It is demonstrated that the phase information can be revealed also for the case of the static, strongly scattering layer. The dependence of the quality of imaging process on the spatial light coherence is demonstrated. The theoretical calculations and numerical simulations are supported by experimental data gained with a model phase object, as well as living carcinoma cells treated in an optically turbid emulsion.
Kollarova, Collakova, Dostal, Vesely, and Chmelik: Quantitative phase imaging through scattering media by means of coherence-controlled holographic microscope



Imaging through scattering media is an important task in biomedical as well technical research and applications. It is solved by various methods based on, e.g., wavefront shaping,1,2 speckle correlations,3,4 inverse scattering,5 time gating,6,7 or optical coherence tomography,8 which can reveal information about the object.

In recent decades, digital holographic microscopy based on transmission mode has made quantitative phase imaging (QPI) of living cells and other phase objects910.11.12.13 possible, and its potential and applications are further investigated. At the same time, it has been demonstrated that digital holography makes possible imaging of objects also through scattering media. The techniques are often based on a coherence gating effect and utilization of ballistic light. The aim is to sort out nonscattered (ballistic) or slightly scattered (snake-like) light to construct a high-quality image, e.g., using the first arriving light method for pulse lasers14 or spatial filtering of c-w lasers with coherent light.15 Recently, the possibility of imaging by means of a coherent light through a flowing turbid medium based on a Doppler frequency shift16 and through a quasistationary turbid medium by multiple holographic acquisitions was studied.17 A combination of digital holography with two-point intensity correlation18 or phase-shifting digital holography using temporally low-coherence light19 to recover the phase information through a scattering medium was also presented.

Imaging through scattering media is possible also by means of spatially incoherent light as reported by Leith.14 The coherence-controlled holographic microscope20 (CCHM) makes possible imaging through a flowing turbid as well as static scattering medium D by means of the spatially incoherent light, utilizing ballistic and diffuse light. The microscope is a two-arm, off-axis holographic system with a halogen lamp as a source [see Fig. 1(a)]. The temporal and spatial coherence of light can be modified by interference filter IF and by a changeable size of an aperture A. Thus, the quality of the image information can also be changed. The system is particularly used in the spatially incoherent regime. Besides the benefits of better resolution and no speckles in the image,20 observation through a scattering layer or medium is then possible due to the coherence gating effect.21 A point object imaging through a diffuse layer was theoretically described recently.21,22 It was shown that the image can be taken only locally invariant.22 The experiments were carried out with an amplitude object, and it was demonstrated that the imaging of a two-dimensional (2-D) amplitude object is possible not only by means of ballistic light, but also by means of diffuse (scattered) light.21 However, neither theoretical nor experimental study of the influence of the scattering layer on the phase information was performed and discussed.

Fig. 1

(a) Scheme of the microscope: LS, light source; A, aperture of changeable size; IF, interference filter; L, collector lens; C, condenser lens; R, reference plane; S, sample (object) plane; D, diffuser; O, objective lens; TL, tube lens; DG, diffraction grating; OL, output lens; OP, output plane; BS, beam splitter; M, mirror; zD, distance between the diffuser and the object plane; rf, transverse shift of the objective lens. (b) Simplified model of the imaging process through the object arm. A point object is imaged through the thin diffuse layer D by means of an objective lens O to the output plane. The position of the point object in the object plane is described by the coordinate vector rSt=(xS,yS); coordinates in the plane of the diffuser are rDt=(xD,yD); and coordinates in the output (image) plane reduced by objective magnification M are rt=(x,y). The reference arm is optically identical, however, without the diffuser and object. (c) Function s(rt) acts like a filter of pD(rt); the transverse shift rf of s(rt) makes it possible to change the profile of wD(rt) (see the text).


This paper proves the possibility of QPI of 2-D phase objects through a static phase scattering layer by means of CCHM. The imaging of an object with defined phase is numerically as well as experimentally evaluated for weak and strong diffuse layers. The paper shows that the phase values can be quantitatively reconstructed, however, in some cases in a restricted region only. The authors demonstrate here the benefits of the spatially incoherent light for the imaging, especially for the case of strong diffusers. Finally, a time-lapse experiment with a 4-mm-thick layer of turbid medium demonstrates the possibility to observe living cells through scattering media.


Two-Dimensional Object Imaging Through a Scattering Layer—Mathematical Description

The three-dimensional imaging properties of the CCHM [Fig. 1(a)] have been already presented in detail.22 The optical model was simplified in the calculations; only the condenser and objective lenses were assumed. The calculations were based on the Fourier optics concept, assuming elastic and linear scattering. Moreover, a point object imaging through a thin 2-D diffuse layer was discussed there in the parabolic approximation.

This study investigates observation of a 2-D object placed in front of a thin diffuse layer, in the parabolic approximation. The light is supposed to be quasi-monochromatic with a central wavelength λ, however spatially incoherent. According to the microscope arrangement, interference of the object and reference beams is recorded at the output plane. The information about the object amplitude and phase is carried in the mixed term of the interference law.22

According to Chmelik et al.,22 imaging of a single point placed in the position rSt=(xS,yS) at the object plane through a single, thin scattering layer at a distance zD behind the object plane [Fig. 1(b)] can be expressed by the function



Vector rt=(x,y) represents coordinates in the output plane, and rf=(xf,yf) is the mutual transverse shift of the images formed by the reference and object beam at the output plane, both reduced by objective magnification M. The symbol * denotes complex conjugation. Function s(rt) depends on the optical properties of the reference arm and the coherence properties of the source.22 It represents the mutual intensity function and can be approximated by the Bessel function of the first kind J1(μ) for the case of waves propagating at small angles, when the pupil functions of lenses can be approximated by the circ function.22 Then,


with μ=2π(NAS/λ)|rt|. NAS is the minimum of the numerical apertures of the source, condenser lens, and objective lens in the reference arm, and affects the degree of the spatial coherence. The width of the peak of s(rt) decreases with increasing NAS, which is equivalent to decrease of the spatial coherence. Function pD(rt) in Eq. (1) is related to imaging of the sample through the object arm and has the following form:22



According to Chmelik et al.,22 the function p0(rt) is for the assumption of small angles propagation described by p0(rt)=[2J1(ν)]/ν, where ν=2π(NA/λ)|rt|, NA is the numerical aperture of the objective lens in the object arm, and K=1/λ. Function TD represents the angular spectrum of the complex transmission function tD(rDt) of the 2-D diffuse layer defined by the Fourier transform



Here, QDt is a scattering vector of the diffuse layer and rDt=(xD,yD) is the vector of coordinates in the diffuser plane. According to Eqs. (1) and (3) and the diffuser properties defined by tD(rDt), the image wD of a point object depends on the point object position rSt.22 The imaging is not spatially invariant; therefore, the function cannot be understood as an impulse response function. However, imaging can be considered to be locally isoplanatic (linear and spatially invariant). The stronger the scattering properties of the diffuse layer, the more significant the shape change of wD and the smaller the region of its spatial invariance.

Imaging of a 2-D object defined by a transmission function tS(rSt) through a scattering layer can be described by the integral



The mutual intensity s(rt) [Eq. (2)] involved in wD in Eq. (5) reflects the degree of spatial coherence. In the case of spatially incoherent light, s(rt) exhibits a significant narrow peak and acts as a filter of function pD(rt) [see Eq. (1) and Fig. 1(c)]. The light from a certain object point therefore interferes only with the light from the corresponding point of the reference object plane (the so-called coherence gating effect). Due to it, the image can be created by ballistic (for rf=0) as well as by diffuse light (rf0) in the presence of diffuse medium. In the case of a high spatial coherence, when NAS0, the function s(rt) is almost constant and equal to 1 and no filter effect occurs. The light from all points of the reference and image plane then interferes without restriction and speckles are created. The image of a single point object [Eq. (1)] is then determined by pD(rt) [Eq. (3)] completely.


Two-Dimensional Phase Object Imaging Through a Scattering Layer—Numerical Simulation

Numerical simulation of a 2-D phase object imaging through a scattering layer was based on Eq. (5). The function wD in Eq. (1) was numerically calculated in MATLAB® for every point of the object plane and for the defined zD and rf (local invariance was not factored in). It represented a four-dimensional transmission matrix between the input and output planes. If necessary, the calculated phase of the image u(rt) was unwrapped using the MATLAB® unwrap function. In case it failed, Goldstein unwrapping algorithm was used.23

The complex transmission function tD(rDt) of the diffuser was designed with regard to the diffuser properties used in the experimental part. It was modeled by a phase function defined by a matrix of N×N squared facets with the side length d=λ/4 and random phase shifts φDjk:



The symbol * denotes a convolution here, J=(N1)/2, and the rectangle function “rect” is defined as rect(x)=1 for |x|<d/2; rect(x)=0.5 for |x|=d/2; otherwise, it is zero. The phases φDjk generated randomly by the computer were characterized by a Gaussian statistical distribution with a standard deviation σD. It means that various diffuser realizations generated with the same σD were different. This study presents two diffusers simulated with N=800, λ=0.65μm and defined by standard deviation 4π/5 and 9π/5rad. They differ by the amount of ballistic and diffuse light. The first diffuser (denoted as a weak one) transmits a significant amount of the ballistic (unscattered) light with regard to the diffuse (scattered) light. The second diffuser (denoted as a strong one) scatters into all directions without preference of the ballistic light direction.

A numerical simulation of imaging of a 2-D object with the phase φS=arg[tS(rSt)], Fig. 2(a), through a scattering layer placed zD=3mm behind the object by means of CCHM is shown in Fig. 2. The authors considered numerical aperture of the objective lenses and of the source NA=NAS=0.25, λ=0.65μm. In Fig. 2(b), a weak diffuser with σD=4π/5rad is used and reconstructed phase φ=arg[u(rt)] is shown. In the case of weak diffusers, the ballistic (unscattered) light is significant in comparison to the diffuse light. Therefore, the best image is formed always by the ballistic light. The fine details of the image can be lost because of the signal-to-noise ratio decrease; however, the quality of the image is comparable to imaging without a diffuser and the image can be obtained without difficulties.

Fig. 2

Numerical simulation of reconstructed phase φ=arg(u) obtained by coherence-controlled holographic microscope, when a phase object φs=arg(ts) (a) is imaged through a scattering layer placed 3 mm behind the layer. Phase image φ obtained (b) through a weak diffuser with σD=4π/5rad, (c) through a strong diffuser with σD=9π/5rad and zero transverse shift rf=(0,0)μm, (d) through a strong diffuser with σD=9π/5rad and rf=(26.5,7.0)μm. The transverse shift rf can significantly change the quality of the signal in the case of strong diffusors—compare (c) and (d).


In contrast, the quality of imaging with ballistic and diffuse light is similar in the case of strong diffusers, and mostly so noisy that the phase cannot be unwrapped. However, a proper mutual shift rf0 of the images formed by the object and reference beam can induce significant improvement of the signal, as demonstrated in Fig. 2(d) in comparison to Fig. 2(c). Figures 2(c) and 2(d) represent the reconstructed phase image φ=arg[u(rt)] through a strong diffuser with σD=9π/5rad for the zero shift and rf=(26.5,7.0)μm, respectively.

The improvement of the image quality for rf0 happens if the function wD [Eq. (1)] becomes more narrow and smooth, exhibiting a significant principal maximum with low side lobes; thus, leading to a good-quality function wD. In the case of spatially incoherent light, the function s [Eq. (2)] has a significant principal maximum. Hence, the quality of wD is determined by the shape of function pD [Eq. (3)] and by the mutual shift rf of the functions s and p. If there is at least a slightly significant direction of scattering QDtf in TD [Eq. (4)], pD [Eq. (3)] exhibits a small local maximum in some position, see Fig. 1(c). A proper transverse shift rf eliminates QDtfzD/K in Eq. (3) and leads to overlapping of s and p maxima and a good quality of the function wD.

The good-quality function wD can appear only in a part of the field of view, as will be demonstrated in Sec. 4. The phase can then be unwrapped at least in the relevant part of the image with the good-quality wD. The improvement in the phase reconstruction is also reflected by the increase of the amplitude |u(rt)|. Finding the highest reconstructed amplitude image leads to the best possible phase reconstruction in the same region. For the nonzero shift rf0, the phase background is tilted according to Eqs. (1) and (3) and has to be compensated.

The effect of the degree of spatial coherence on the image quality due to the coherence gating effect is illustrated in Fig. 3. It demonstrates the decrease of image quality with increasing spatial coherence of the light in the presence of a diffuser for zD=3mm. Figures 3(a)3(d) represent imaging with a weak diffuser with σD=4π/5rad; Figs. 3(e)3(h) represent imaging with a strong diffuser with σD=9π/5rad and the shift rf=(26.5,7.0)μm, leading to the best possible image quality for this diffuser realization (both diffusers are expressed by the same phase transmission function as in Fig. 2). The effect of the increasing spatial coherence of light on the amplitude |u| is demonstrated in Figs. 3(a) and 3(c) for the weak diffuser and in Figs. 3(e) and 3(g) for the strong diffuser. The same effect on phase φ=arg(u) is simulated in Figs. 3(b) and 3(d), and Figs. 3(f) and 3(h) for the weak and strong diffusers, respectively. The low spatial coherence is determined by NAS=0.25 for Figs. 3(a), 3(b), 3(e), and 3(f), and high spatial coherence by NAS=0.01 for Figs. 3(c), 3(d), 3(g), and 3(h). In the case of high spatial coherence and a weak diffuser, Fig. 3(d) is affected by the coherence noise, the finer details are lost, and the phase unwrapping is a more complicated task. In case the MATLAB® unwrapping function failed, the Goldstein algorithm was successful. For the strong diffuser and high spatial coherence, the phase cannot be reconstructed, as clearly seen in Fig. 3(h). A similar situation also occurs for the weak diffuser if it is placed closer to the object. However, the low spatially coherent light makes QPI possible in all cases [see Figs. 3(b) and 3(f)]. The influence of noise on the phase reconstruction is discussed in more detail in Sec. 4.

Fig. 3

Numerical simulation of the spatial coherence impact on the image amplitude |u| and the image phase φ=arg(u) for NA=0.25, zD=3mm. The object phase is represented by Fig. 2(a). (a) to (d) Images obtained through a weak diffuser with σD=4π/5rad: (a) reconstructed amplitude and (b) phase for spatially incoherent light characterized by NAS=0.25; (c) reconstructed amplitude and (d) phase for almost spatially coherent light characterized by NAS=0.01. (e) to (h) Images obtained through a strong diffuser with σD=9π/5rad and rf=(26.5,7.0)μm leading to the best reconstructed image quality, see Fig. 2(d). (e) Amplitude and (f) phase reconstruction for NAS=0.25; (g) amplitude and (h) phase reconstruction for NAS=0.01. Due to higher coherence noise, phase unwrapping is not successful in the case of spatially coherent light (h).



Quantitative Phase Imaging Through a Scattering Layer—Experiments

The following experiments verify the possibility of quantitative phase imaging (QPI) through static scattering media by means of CCHM20 [Fig. 1(a)]. A halogen lamp was used as a source, interference filter 650 nm with FWHM 10 nm (corresponding to the coherence length 42μm), objective lenses 10×/0.3, and spatially incoherent light characterized by NAS=0.3 [aperture A in Fig. 1(a) was fully open]. The optical magnification M between the object and image plane was 40×. The weak and strong diffusers in the form of a microroughened surface were created on one side of coverslips by different times of sandblasting and placed 0.8 mm behind the object. The phase object was formed in a glass coverslip by squares with defined depths leading to the phase shifts of 0.1, 0.35, 4, 6, and 8 rad for λ=650nm. In the first two columns of Table 1, the mean phase φm and the standard deviation σ of the phase values within the squares measured by means of the CCHM without a diffuser are summed. The surface of the sample was not smooth; therefore, σ is quite high.

Table 1

Comparison of the mean phase φm and the standard deviation σ of the phase values φ measured in the largest squares without and with a weak diffuser by means of a coherence-controlled holographic microscope.

Without diffuserWith a weak diffuser
φm (rad)σ (rad)φm (rad)σ (rad)

The obtained holograms were processed using the CCHM software, in which the Goldstein algorithm23 was chosen to unwrap the reconstructed phase and the phase background was then compensated.24 Although the phase shift of 8 rad is greater than 2π, the algorithm was mostly successful in the phase unwrapping (it depended also on the quality of focusing). There was obviously some evidence about the character of the phase changes at the square edge; therefore, the unwrapping was correct. In the case of an abrupt change, there is no certainty about the correct phase. However, it is possible to verify it by another measurement, e.g., with two different central wavelengths and to use the results for the multiwavelength digital holography.25 CCHM is an achromatic system allowing this kind of measurement.

Imaging of the phase object through a weak diffuser is presented in Fig. 4. Figure 4(a) illustrates the wrapped image phase; Fig. 4(c) depicts the unwrapped and compensated image phase of the object with phase shifts of 0.1, 0.35, and 4 rad. The phase distribution along the black line is compared with that obtained by CCHM without the diffuser [Fig. 4(e)]. Similarly, Fig. 4(b) illustrates the wrapped image phase, while Fig. 4(d) shows the unwrapped and compensated image phase of squares with phase shifts of 4, 6, and 8 rad through the weak diffuser. Fig. 4(f) compares the measured phase along the black line in Fig. 4(d) with and without a diffuser. The third and fourth columns of Table 1 sum up the mean phases φm and the standard deviations σ of the phases measured in the interior of the biggest squares in Figs. 4(c) and 4(d). It is clear from Fig. 4 and Table 1 that the phase shifts of 0.35 rad and greater are imaged well, however, with a higher standard deviation of the phase. Finer details are lost, as obvious for the phase of 0.098 rad. The QPI is possible in the whole field of view without difficulties and the best image is always formed by ballistic light.

Fig. 4

Experimental verification of quantitative phase imaging (QPI) of square objects through a weak diffuser. (a) The wrapped and (c) compensated unwrapped phase image of squares with phase shifts of 0.1, 0.35, and 4 rad. (e) Comparison of the measured values along the black line in (c) with QPI values without the diffuser. (b) The wrapped and (d) compensated unwrapped phase image of the squares with phase shifts of 4, 6, and 8 rad. (f) Comparison of the values along the black line in (d) with the reference QPI values obtained without the diffuser.


Imaging through a strong diffuser is presented in Fig. 5. Figures 5(a)5(c) are the reconstructed amplitudes, Figs. 5(d)5(f) are the wrapped phases, Figs. 5(g)5(i) are compensated unwrapped phases, and Figs. 5(j)5(l) are phase distributions along the white line in Figs. 5(g)5(i) compared with the values obtained by CCHM without the diffuser. The first column corresponds to the zero mutual shift rf (imaging by means of the ballistic light); the second and third columns show imaging by means of diffuse light with a shift of rf=(3,12)μm and rf=(10,3)μm, respectively. In the experiments, rf was realized by the shift of the reference objective lens.

Fig. 5

Experimental verification of QPI of square objects introducing phase shifts 4, 6, and 8 rad through a strong diffuser. The reconstructed amplitude image |u| for (a) rf=(0,0)μm, (b) rf=(3,12)μm, and (c) rf=(10,3)μm. The wrapped phase image for (d) rf=(0,0)μm, (e) rf=(3,12)μm, and (f) rf=(10,3)μm. The compensated unwrapped phase image for (g) rf=(0,0)μm, (h) rf=(3,12)μm, and (i) rf=(10,3)μm. (j) to (l) Comparison of the phase values obtained without the diffuser with the values along the white line in (g) to (i).


In the case of imaging through strongly scattering media, the reconstructed amplitude is weak and more noisy. The areas with a weak reconstructed amplitude correspond to the areas with high phase noise.26 If the phase differences in the adjacent points due to noise are greater than π, the 2π shifts of the wrapped phase corresponding to the image are badly localizable. The unwrapping then fails, the compensated phase is destroyed, and information about the image phase is lost [see Figs. 5(a), 5(d), and 5(g)]. This is equivalent to the situation when the function wD does not exhibit a significant principal maximum.

However, as discussed in Sec. 3, CCHM makes it possible to change the quality of the function wD (and therefore also of the image) by the mutual shift rf of the reference and object beam at the output plane. Then, in the region with the good-quality function wD, the local improvements of the modulus of amplitude, Figs. 5(b)5(c), and decrease of the amount of the noise in the wrapped phase, Figs. 5(e)5(f), appear. Then, the phase can be unwrapped and compensated more successfully [see Figs. 5(e), 5(f), 5(h), and 5(i)]. The measurement is influenced by a greater error than in the case of weak diffuser; however, the mean values are relevant [see the regions x=(0:100)μm in Fig. 5(k) and x=(100:200)μm in Fig. 5(l)]. This region of the QPI quality is surrounded by a part of field in which the phase is unwrapped and visualized, however is not relevant [see the regions x=(125:260)μm in Fig. 5(k) and x=(0:100)μm and x=(200:280)μm um in Fig. 5(l)]. Some parts of Figs. 5(h) and 5(i) cannot be unwrapped because the quality of the function wD is not sufficient here.

Figure 6 illustrates the amplitude and phase imaging through a weak diffuser placed 0.8 mm behind the sample for the case of spatially incoherent light, and almost fully coherent light. Figures 6(a)6(c) present reconstructed amplitude, wrapped phase, and compensated unwrapped phase, respectively, observed by means of CCHM without a diffuser and with spatially incoherent light. In the other cases, a weak diffuser was used. Figures 6(d)6(f) show the reconstructed amplitude [Fig. 6(d)], wrapped phase [Fig. 6(e)], and compensated wrapped phase [Fig. 6(f)] with spatially incoherent light. In both cases, the changeable aperture was fully open and the illuminating beam filled the condenser pupil, i.e., NA=NAS=0.3; the exposure time was 100 ms. The measured phase [Fig. 6(f)] was comparable to the results [Fig. 6(c)] for all three square objects [see Table 1]. Figures 6(g)6(i) illustrate the amplitude, wrapped, and compensated unwrapped image phase with increased spatial coherence of light. The area of the aperture was reduced to 0.1%, which corresponds to NAS=0.02; the exposure time was 8 s. The speckles visible in Fig. 6(g) made unwrapping of the image phase [Fig. 6(h)] impossible. Obviously, the spatially incoherent light can significantly improve the quality of the obtained CCHM image due to the coherence gating effect (in contrast to the spatially coherent illumination) and thus make QPI possible also through scattering media.

Fig. 6

Experimental verification of amplitude imaging and quantitative phase imaging through a weak diffuser in dependence on the spatial coherence of light. Observation of a square object defined by phase shifts of 4, 6, and 8 rad and the square side length of 50μm is presented. (a) to (c) Amplitude, wrapped phase, and compensated unwrapped phase observed without the diffuser. (d) to (f) Amplitude, wrapped phase, and compensated unwrapped phase observed with the diffuser in spatially incoherent light, NA=NAS=0.3. (g) to (i) Amplitude, wrapped phase, and compensated unwrapped phase observed with a diffuser and spatially coherent light, NAS=0.02.



Experiment with Living Cells

CCHM is primarily used for experiments with biological objects. The coherence gating effect makes possible observation of these objects when they are treated in turbid or stationary scattering media. The time-lapse QPI of colorectal carcinoma cells by CCHM and corresponding unwrapped phase images are shown in Fig. 7. Living cells under the 4-mm-thick layer of 0.15% bioactive phospholipid emulsion (BAP) in the cultivation medium are barely visible in the Zernike phase contrast, while in CCHM QPI, one can recognize intracellular components. In the series, it is clearly seen that BAP affects not only the cells at the colony edge, but also the cells inside the colony are hit almost at the same time. This indicates that BAP strikes cancer cells indiscriminately of the protection provided by full contact with neighbor cells. A halogen lamp with interference filter 650 nm with FWHM 10 nm, objective lenses 20×/0.5, and spatially incoherent light characterized by NAS=0.29 were used in the experiment.

Fig. 7

Colorectal carcinoma cells DLD1 treated with bioactive phospholipids (0.15% bioactive phospholipid emulsion) with anticancer properties emulsified in growth medium. The cells were observed under a 4-mm-thick layer of turbid medium: (a) just after the application and (b) after 6.5 h of exposure. (c) Final situation showing dying cells after 19.5 h. In the calibration scale, 1 rad is equivalent to 0.57pg/μm2. (d), (e), and (f) demonstrate corresponding wrapped phase reconstructed from the holograms for the cases (a), (b), and (c).




This study presents a mathematical description of 2-D object imaging through a thin phase scattering layer by means of CCHM. The results of quantitative phase imaging were numerically simulated and experimentally verified for a weak and strong diffuser characterized by a different ratio of amounts of ballistic and diffuse light. The paper has shown that the phase image can be formed both by ballistic and diffuse light, depending on the diffuser properties. In the case of a weak diffuser, the best signal is always determined by the ballistic light, with resolution comparable or slightly worse than without the diffuser. The quality of the image behind a strong diffuser depends on the mutual shift rf of the reference and object beam images at the output plane. In some cases, the object can be imaged by both ballistic (rf=0) and diffuse (rf0) light, and in other cases by the diffuse light only with surprisingly good image quality. The signal is, however, more noisy, the resolution worse, and in some cases, the imaging can be rather a phase visualization technique. Even the best choice of rf can usually lead to a good-quality QPI in a limited area of the image field. However, this area can be extended recombining images for several rf.

The authors have also demonstrated the influence of the degree of spatial coherence on the quality of the image created through a thin scattering layer, for weak and strong diffusers. The numerical simulations and experimental results confirm the benefits of the spatially incoherent light for imaging through scattering media by means of CCHM. The results can be useful in biomedical applications and research, for example, for assessment of reactions of carcinoma cells to various drugs, which exhibit scattering properties.


This work was supported by the project EXCELLENT TEAMS (CZ.1.07/2.3.00/30.0005) from European Social Fund, “CEITEC—Central European Institute of Technology” (CZ.1.05/1.1.00/02.0068) from European Regional Development Fund, and by the project 15-14612S from the Czech Science Foundation.



I. Freund, “Looking through walls and around corners,” Phys. A 168(1), 49–65 (1990).http://dx.doi.org/10.1016/0378-4371(90)90357-XGoogle Scholar


O. Katz et al., “Focusing and compression of ultrashort pulses through scattering media,” Nat. Photon. 5, 372–377 (2011).NPAHBY1749-4885http://dx.doi.org/10.1038/nphoton.2011.72Google Scholar


J. Bertolotti et al., “Non-invasive imaging through opaque scattering layers,” Nature 491(7423), 232–234 (2012).http://dx.doi.org/10.1038/nature11578Google Scholar


M. Vellekoop and C. M. Aegerter, “Scattered light fluorescence microscopy: imaging through turbid layers,” Opt. Lett. 35(8), 1245–1247 (2010).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.35.001245Google Scholar


H. P. Baltes, Inverse Scattering Problems in Optics, Springer-Verlag, Berlin (1980).Google Scholar


C. Dunsby and P. M. W. French, “Techniques for depth-resolved imaging through turbid media including coherence-gated imaging,” J. Phys. D: Appl. Phys. 36, 207–227 (2003).JPAPBE0022-3727http://dx.doi.org/10.1088/0022-3727/36/14/201Google Scholar


V. T. Willwacher et al., “Recovering three-dimensional shape around a corner using ultrafast time-of-flight imaging,” Nat. Commun. 3, 745 (2012).NCAOBW2041-1723http://dx.doi.org/10.1038/ncomms1747Google Scholar


M. Sticker et al., “Quantitative differential phase measurement and imaging in transparent and turbid media by optical coherence tomography,” Opt. Lett. 26(8), 518–520 (2001).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.26.000518Google Scholar


P. Marquet et al., “Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy,” Opt. Lett. 30(5), 468–470 (2005).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.30.000468Google Scholar


G. Popescu et al., “Erythrocyte structure and dynamics quantified by Hilbert phase microscopy,” J. Biomed. Opt. 10(6), 060503 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2149847Google Scholar


W. C. WargerII et al., “Phase-subtraction cell-counting method for live mouse embryos for live mouse embryos beyond the eight-cell stage,” J. Biomed. Opt. 13(3), 034005 (2008).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2937468Google Scholar


H. Janeckova, P. Vesely and R. Chmelik, “Proving tumour cells by acute nutritional/energy deprivation as a survival threat: a task for microscopy,” Anticancer Res. 29(6), 2339–2345 (2009).ANTRD40250-7005Google Scholar


J. Balvan et al., “Multimodal holographic microscopy: distinction between apoptosis and oncosis,” PLoS One 10(3), e0121674 (2015).POLNCL1932-6203http://dx.doi.org/10.1371/journal.pone.0121674Google Scholar


E. Leith et al., “Imaging through scattering media with holography,” J. Opt. Soc. Am. A 9(7), 1148–1153 (1992).http://dx.doi.org/10.1364/JOSAA.9.001148Google Scholar


S. Li and J. Zhong, “Dynamic imaging through turbid media based on digital holography,” J. Opt. Soc. Am. A 31(3), 480–486 (2014).http://dx.doi.org/10.1364/JOSAA.31.000480Google Scholar


M. Paturzo et al., “Microscopy imaging and quantitative phase contrast mapping in turbid microfluidic channels by digital holography,” Lab Chip 12(17), 3073–3076 (2012).LCAHAM1473-0197http://dx.doi.org/10.1039/C2LC40114BGoogle Scholar


V. Bianco et al., “Clear coherent imaging in turbid microfluidics by multiple holographic acquisitions,” Opt. Lett. 37(20), 4212–4214 (2012).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.37.004212Google Scholar


R. K. Singh, A. M. Sharma and B. Das, “Quantitative phase-contrast imaging through a scattering media,” Opt. Lett. 39(17), 5054–5057 (2014).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.39.005054Google Scholar


S. Tamano, Y. Hayasaki and N. Nishida, “Phase-shifting digital holography with a low-coherence light source for reconstruction of a digital relief object hidden behind a light-scattering medium,” Appl. Opt. 45(5), 953–959 (2006).APOPAI0003-6935http://dx.doi.org/10.1364/AO.45.000953Google Scholar


T. Slabý et al., “Off-axis setup taking full advantage of incoherent illumination in coherence-controlled holographic microscope,” Opt. Express 21(12), 14747–14762 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.014747Google Scholar


M. Lošt’ák et al., “Coherence-controlled holographic microscopy in diffuse media,” Opt. Express 22(4), 4180–4195 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.004180Google Scholar


R. Chmelik et al., “The role of coherence in image formation in holographic microscopy,” Prog. Opt. 59, 267–335 (2014).POPTAN0079-6638Google 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


.T. Zikmund et al., “Sequential processing of quantitative phase images for the study of cell behavior in real-time digital holographic microscopy,” J. Microsc. 256(2), 117–125 (2014).JMICAR0022-2720http://dx.doi.org/11.1111/jmi.12165Google Scholar


J. Gass, A. Dakoff and M. K. Kim, “Phase imaging without 2π ambiguity by multiwavelength digital holography,” Opt. Lett. 28(13), 1141–1143 (2003).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.28.001141Google Scholar


J. W. Goodman, Statistical Optics, Wiley-Interscience, New York, pp. 498–500 (2000).Google Scholar


Vera Kollarova is a postdoctoral fellow at Brno University of Technology in the Czech Republic. She received her BS and MS degrees in optics from Palacký University in Olomouc in 2001 and 2003, respectively, and a PhD in optics from the same university in 2011. Her current research interests include holographic microscopy and imaging through scattering media.

Jana Collakova is a PhD student at the Brno University of Technology. She received her BS degree in optics from the Palacky University in Olomouc in 2009 and her MS degree in optics and precise mechanics from the Brno University of Technology in 2011. Her current research interests include digital holographic microscopy in turbid media and living cells imaging.

Zbynek Dostal is a PhD student at the Brno University of Technology. He received his BS degree in mechanical engineering in 2007 and his MS degree in optics and precise mechanics in 2009 from the Brno University of Technology. He is a co-author of the CCHM patent for which he received the Werner von Siemens Excellence Award 2013. His current research interests include holographic microscopy, automation in microscopy, optical and mechanical design.

Pavel Vesely is a senior research scientist at the CEITEC, Brno University of Technology. He received his MD degree from the Charles University in Prague in 1961 and his PhD in experimental biology from the Czechoslovak Academy of Sciences in 1965. He is the author of more than 100 journal papers. His life long research interests in cancer cell biology currently include domestication of coherence controlled holographic microscopy in cell biology research.

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 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 3D light microscopy, and holographic microscopy.

Vera Kollarova, Jana Collakova, Zbynek Dostal, Pavel Vesely, Radim Chmelik, "Quantitative phase imaging through scattering media by means of coherence-controlled holographic microscope," Journal of Biomedical Optics 20(11), 111206 (5 August 2015). http://dx.doi.org/10.1117/1.JBO.20.11.111206


Scattering media

Image quality


Phase imaging

Light scattering

Spatial coherence

Back to Top