# Determination of the Size of Irregular Particles Using Interferometric Out-of-Focus Imaging

Authors: Marc Brunel; Huanhuan Shen; Sébastien Coëtmellec; Gérard Gréhan; Tony Delobel
Journal: International Journal of Optics
 (2014)

---

## Abstract

We present a mathematical formalism to predict speckle-like interferometric out-of-focus patterns created by irregular scattering objects. We describe the objects by an ensemble of Dirac emitters. We show that it is not necessary to describe rigorously the scattering properties of an elliptical irregular object to predict some physical properties of the interferometric out-of-focus pattern. The fit of the central peak of the 2D autocorrelation of the pattern allows the prediction of the size of the scattering element. The method can be applied to particles in a size range from a tenth of micrometers to the millimeter.


---

## Body

## 1. Introduction

Interferometric laser imaging for droplet or bubble sizing is a very convenient technique. It is relatively simple to develop and gives accurate evaluations of the size of these scattering particles in a wide range of configurations [1–7]. In addition, real-time algorithms can be developed to investigate a whole image directly using Fourier transforms instead of investigating the interferogram droplet per droplet or bubble per bubble [8]. Unfortunately, this technique is limited to the analysis of simple scattering objects as spherical droplets or bubbles. The deformation of these objects can affect significantly the interferometric patterns (orientation and frequency of the fringes) and the analysis of complex irregular objects does not seem to be possible. In this paper, we show however that interferometric out-of-focus imaging can allow the determination of some simple morphological information of irregular particles (size and global form of the envelope).It is well established that interferometric out-of-focus imaging of droplets is based on the observation of interferences produced by the two “glare points” produced by the illumination beam on the droplet [9, 10]. Under a geometrical optics approximation, these emitters correspond to the ray reflected on the particle and the ray refracted within the droplet [11]. These emitters can be further defined using Lorenz-Mie theory expressed in terms of Debye series [10]. In the case of droplets in air, the scattering angle of 66° is currently chosen. It is indeed the angle leading to the amplitudes of the reflected beam and the refracted beam having the same magnitude. The interferences that are recorded present then a contrast reaching unity.For irregular particles, there is unfortunately no general solution to the field scattered by the particles when they are illuminated by a coherent plane wave, Gaussian beam, or laser sheet. The interferometric out-of-focus images are then more or less a speckle-like pattern. It is not possible to represent the scattering objects by a precise number of well-defined glare points as can be done with droplets or bubbles. To illustrate this, Figure1 shows the interferometric out-of-focus image recorded on a CCD sensor. The experimental setup is as shown in Figure 2. A frequency-doubled Nd:YAG laser sheet illuminates NaCl salt crystals (the laser wavelength is 532 nm). When observed with a microscope, these crystals have a global ellipsoidal form. The field scattered by these particles around the angle Θ=66° is imaged with a lens (focus length f=10 cm). The distance between the crystals and the aperture is z1=8.5 cm, the distance from the aperture to the lens is z2=6.5 cm, and the distance from the lens to the CCD sensor is z3=20 cm. The radius of the aperture is R0=4.8 mm. We can observe in Figure 1 a speckle-like pattern. Although it is not possible to describe rigorously the field scattered by the irregular scattering NaCl salt crystals, it does not mean that the speckle pattern does not contain rigorous information concerning the scattering particle.Figure 1
Experimental interferometric out-of-focus image of a NaCl salt crystal. All axes are in meters.Figure 2
Experimental setup (schematic representation; the setup is not to scale).In this paper, we will first develop in Section2 a mathematical formalism to predict speckle-like interferometric out-of-focus patterns created by irregular objects. The objects that we will describe are simplified to an ensemble of Ngp Dirac emitters located on a half sphere or more generally on a half ellipsoid. We will show that we can predict speckle-like patterns that are observed experimentally. We will show that the size of the bright speck of the speckle patterns does not depend on the relative location of the glare points but depends on the dimension of the global half sphere or of the half ellipsoid. Finally, using this formalism, we will show in Section 3 that we can deduce the size of irregular NaCl salt crystal from the out-of-focus pattern of Figure 1.


## 2. Simulation of Speckle-Like Interferometric Out-of-Focus Patterns

From a theoretical point of view, the field emitted by a point Dirac emitteri, whose transverse position in the plane of the scattering particle is (ai,bi), can be evaluated in the plane of the CCD sensor by generalized Huygens Fresnel integrals [12–15]. Their coefficients depend on the coefficients of optical transfer matrices describing the optical elements encountered between the particle and the sensor. The procedure has been detailed in references [7, 16] in the case of droplets represented by a pair of Dirac emitters and has been validated in all cases that could be tested experimentally [7, 17, 18]. For a unique Dirac emitter, we obtain the following expression of the electric field in the plane of the CCD sensor [7]:
(1)Gi(x,y,z)=(×((iλ)2B1,ixB1,iyB2xB2x)-1(αieiφiei(2π/λ)(n1z1,i+n2(z2+z3))ei(π/λ)((D2xx2/B2x)+(D2yy2/B2y)))×((iλ)2B1,ixB1,iyB2xB2y)-1)×∑k=1NPkπeβi(k)γx,i(k)γy,i(k),
where z is the optical longitudinal axis of the imaging setup, x and y are the transverse coordinates in the plane of the CCD sensor, and
(2)βi(k)=iπA1,ixai2λB1,ix+iπA1,iybi2λB1,iy-ϕx,i24γx,i(k)-ϕy,i24γy,i(k),γx,i(k)=QkR02-iπλ(D1,ixB1,ix+A2xB2x),γy,i(k)=QkR02-iπλ(D1,iyB1,iy+A2yB2y),ϕx,i(k)=2πλ(aiB1,ix+xB2x),ϕy,i(k)=2πλ(biB1,iy+yB2y).
The different A1,ip, B1,ip, C1,ip, and D1,ip (with p=x or p=y) are the coefficients of the optical transfer matrices between the emitter i and the plane where the aperture is located. n1z1,i is the optical path between the emitter i and the plane of the aperture. We can define two different transfer matrices for both transverse axes x and y, which allows us to describe cylindrical geometries. The transmission coefficient of the aperture equals 1 within a disk of radius R0 and equals 0 outside. It is expressed as a superposition of 10 Gaussian functions, as was proposed by Wen and Breazeale [19]. The corresponding coefficients of this decomposition Pk and Qk can be found in [19]. The A2p, B2p, C2p, and D2p (with p=x or p=y) are the coefficients of the optical transfer matrices between the plane of the aperture and the plane of the CCD sensor. n2(z2+z3) is the optical path between the plane of the aperture and the plane of the CCD sensor.In ILIDS, it has been shown that a droplet can be well represented by two emitters [7, 16]. Let us now consider the case of a multiemitter composed of Ngp Dirac emitters randomly located on a half sphere, as shown in Figure 3(a). The pattern observed on the CCD sensor can then be theoretically calculated. According to previous relations, the expression of the electric field in the CCD sensor plane can be written as follows:
(3)Gtot(x,y,z)=∑i=1NgpGi(x,y,z).
In the next simulations, we will assume that the amplitude of the electric field is the same for all emitters (αi=1 in normalized units). In addition, we are interested in the characterization of particles in the size domain from a tenth of micrometers to hundreds of micrometers. For spherical droplets, the phase shift introduced between the two glare points (which represent the reflected beam and the beam refracted by the droplet without any internal reflection) is then much higher than 2π. In the case of irregular particles, we will thus assume that the phase φi of each emitter is a random variable in the range [0,2π].Half sphere composed of 40 randomly located Dirac emitters (a) and interferometric out-of-focus image created by this object (b). The radius of the half sphere is 20 micrometers. All axes are in meters.
(a)
(b)We present in Figure3 the geometrical repartition of Ngp=40 Dirac emitters randomly located on a half sphere of radius of 20 micrometers (Figure 3(a)) and the interferometric out-of-focus image created by this object (Figure 3(b)) using previous relations (1)–(3). The imaging system as shown in Figure 2 is a lens whose focus length and radius are f=10 cm and R0=10 mm, respectively. The distance between the center of the half sphere and the aperture is z1=13 cm, the distance from the aperture to the lens is z2=0 cm (the lens is the aperture in this case), and the distance from the lens to the CCD sensor is z3=39.3 cm. We obtain a speckle-like pattern as in the case of the experimental image of Figure 1, although the dimension of the speck is different. For comparison, Figure 4 shows the geometrical repartition of Ngp=40 Dirac emitters randomly located on a half sphere of radius of 40 micrometers (Figure 4(a)) and the interferometric out-of-focus image created by this object (Figure 4(b)). Finally, Figure 5 shows the geometrical repartition of Ngp=40 Dirac emitters randomly located on a half sphere of radius 80 micrometers (Figure 5(a)) and the interferometric out-of-focus image created by this object (Figure 5(b)). The imaging system is the same in all three cases of Figures 3, 4, and 5. To characterize the size of the speck of a speckle pattern, it is interesting to plot the 2D-autocorrelation of this speckle pattern. Figures 6(a), 6(b), and 6(c) show the 2D autocorrelations of the three speckle patterns presented in Figures 3(b), 4(b), and 5(b), respectively. We observe a thin central peak which characterizes the size of the bright speck. We can clearly observe that the higher the dimension of the sphere, the smaller the size of the speck of the speckle pattern.Half sphere composed of 40 randomly located Dirac emitters (a) and interferometric out-of-focus image created by this object (b). The radius of the half sphere is 40 micrometers. All axes are in meters.
(a)
(b)Half sphere composed of 40 randomly located Dirac emitters (a) and interferometric out-of-focus image created by this object (b). The radius of the half sphere is 80 micrometers. All axes are in meters.
(a)
(b)2D autocorrelations of the interferometric out-of-focus patterns of Figures3(b), 4(b), and 5(b), respectively. All axes are expressed in meters.
(a)
(b)
(c)It is now important to know whether this characteristic of the speckle-like patterns, that is, the size of the bright speck, depends on the relative location of the emitters to each other. We present thus in Figures7 and 8 a comparison of the results obtained when the Dirac emitters are located randomly or uniformly on a half sphere of radius of 40 micrometers. The setup parameters are those of previous Figures 3, 4, 5, and 6. First, Figures 7(a) and 7(b) present the location of the Ngp=40 glare points on the half sphere when the emitters are randomly located (Figure 7(a)) or uniformly located (Figure 7(b)). We have then calculated the out-of-focus patterns created by these objects using relations (1)–(3). Their 2D autocorrelations are presented in Figures 7(c) and 7(d), respectively. Although both 2D autocorrelations are not identical, we observe a central bright peak whose dimensions seem to be the same in both cases. To quantify this observation exactly, Figure 8(a) shows the x-cut of the central peak of the 2D-autocorrelation in the case of a random location of the emitters (black curve) and in the case of a uniform location of the emitters (red curve). In a similar way, Figure 8(b) shows the y-cut of the central peak of the 2D-autocorrelation in the case of a random location of the emitters (black curve) and in the case of a uniform location of the emitters (red curve). Although the pedestals of the curves are not identical, the width of the central peak is the same for a random or a uniform location of the emitters. This is true for both x-cuts and y-cuts of the peaks. In conclusion, the relative location of the emitters does not significantly affect the dimension of the central peak of the 2D-autocorrelation.Half spheres composed of 40 Dirac emitters randomly located (a) or uniformly located over the half sphere (b). Autocorrelations of the speckle-like interferometric out-of-focus patterns that they create: (c) and (d), respectively. The radius of the half spheres is 40 micrometers. All axes are expressed in meters.
(a)
(b)
(c)
(d)x-cut (a) and y-cut (b) of the central peak of the 2D-autocorrelation of the simulated speckle-like interferometric out-of-focus patterns for a random location of the emitters (black dashed lines) or a regular location of the emitters (red dotted lines). The x-axes are expressed in meters, while the y-axes represent normalized amplitudes (no unit).
(a)
(b)In a similar way, we could further illustrate that the number of emitters to be considered (20, 30, 40,…) does not affect significantly the dimension of the central peak of the 2D-autocorrelation. In conclusion, it is not necessary to describe rigorously the scattering properties of an irregular object to predict some physical properties of the scattering pattern observed in an out-of-focus plane. In particular, the observation of the central peak of the 2D-autocorrelation of the interferometric out-of-focus pattern allows us to predict the dimension of the scattering element by comparison with simulated patterns using our relations (1)–(3).


## 3. Comparison with Experimental Data

Finally, let us consider the experimental out-of-focus pattern presented in Figure1. We have calculated the 2D-autocorrelation of this pattern (not reported here). As previously observed with our simulations, we observe a central peak on a large pedestal. We report in Figure 9(a) the x-cut of the central peak of this 2D-autocorrelation (red curve), while Figure 9(b) shows the y-cut of the central peak of this 2D-autocorrelation (red curve). We can see that the central peak of the 2D-autocorrelation is not circular. A zoom of this central peak shows that it is actually elliptical (with principal axes x and y). It indicates that the scattering object (a NaCl salt crystal) cannot be viewed as a spherical object but as a rough ellipsoidal object as observed with a microscope. In order to simulate these results, we have used our relations (1)–(3) with the parameters that are necessary to simulate the experimental setup of Figure 2. Let us recall these parameters. The illumination wavelength is λ=532 nm. The imaging lens (SL in Figure 2) has a focus length f=10 cm. The aperture D has a radius R0=4.8 mm. The distance between the NaCl salt crystal and the aperture is z1=8.5 cm, the distance from the aperture to the lens is z2=6.5 cm, and the distance from the lens to the CCD sensor is z3=20 cm. In the simulations, the matrices to be considered are matrices M1,i (coefficients A1,i, B1,i, C1,i, and D1,i) describing the first section of the imaging system, which takes into account propagation from the different emitters of the object to the aperture (see Figure 2), and a matrix M2 (coefficients A2, B2, C2, and D2) which takes into account the propagation from the aperture to the lens, the transmission through the lens, and the propagation from the lens to the CCD sensor. The NaCl salt crystal is modeled by Ngp=60 Dirac emitters uniformly located on a half ellipsoid. It is presented in Figure 10(a). The best fit of the experimental data is obtained for dimensions of the half ellipsoid: 800 micrometers along the x-axis and 380 micrometers along the y- and z-axis. Note that the dimension along the z-axis (which is the optical axis of the experiment) has no significant influence on the planar interferometric out-of-focus pattern which is sensitive to the x- and y-dimensions. Figure 10(b) shows a zoom of the central peak of the 2D-autocorrelation of the interferometric pattern that is then calculated in the CCD sensor plane. The central peak is clearly elliptical, due to the ellipticity of the object around the x- and y-axes. Finally, we report in black curves in Figure 9 the x-cut of the central peak of this 2D-autocorrelation (Figure 9(a)) and the y-cut of the central peak of this 2D-autocorrelation (Figure 9(b)). We can see the very good concordance between experimental and simulated results concerning the dimensions of the central peak.x-cuts (a) and y-cuts (b) of the central peak of the 2D autocorrelations of the experimental (dotted red lines) and simulated (solid black lines) speckle-like interferometric out-of-focus patterns. The x-axes are expressed in meters, while the y-axes represent normalized amplitudes (no unit).
(a)
(b)Half ellipsoid composed of 60 Dirac emitters uniformly located (a) and zoom of the central peak of the 2D-autocorrelation of the interferometric out-of-focus image created by this object (b). All axes are in meters.
(a)
(b)In conclusion, we confirm experimentally that the observation of the central peak of the 2D-autocorrelation of the interferometric out-of-focus scattering pattern allows us to predict the dimensions of the scattering element by comparison with simulated patterns using our relations (1)–(3). This can be done despite not being able to describe rigorously the scattering properties of the rough irregular object. These first experiments show the potentiality of this technique with relatively simple irregular objects: ellipsoidal rough salt crystals. Future studies should show whether it is possible to characterize the morphology of more complex irregular objects.It is very important to note that the defocus parameter needs to be sufficiently large, in such a way that the defocused images of all emitters overlap in the global out-of-focus image of the particle observed in the CCD sensor plane. This point is similar to classical ILIDS, where the defocus parameter needs to be sufficiently large, in such a way that the defocused images of both reflected and refracted glare points (in the case of a droplet) overlap to deliver interference fringes. This condition is respected in the results presented here.


## 4. Conclusion

In this paper, we have developed a mathematical formalism to predict speckle-like interferometric out-of-focus patterns created by ellipsoidal rough scattering objects. The objects that we describe are simplified to an ensemble ofNgp Dirac emitters located on global forms as half spheres or half ellipsoids. We could confirm experimentally that it is not necessary to describe rigorously the scattering properties of an elliptical irregular object to predict some physical properties of the scattering pattern observed in an out-of-focus plane. In particular, the observation of the central peak of the 2D-autocorrelation of the interferometric out-of-focus pattern allows us to predict the dimensions of the scattering element by comparison with simulated patterns using our formalism. These first experiments show the possibility to characterize relatively simple irregular objects: ellipsoidal rough salt crystals. Future studies should show whether it is possible to characterize the morphology of more complex irregular objects. The formalism developed allows us to describe any imaging system as it is expressed versus the coefficients of the optical transfer matrices that describe the system. This method could have a wide number of applications, for example, the characterization of ashes in the atmosphere, or sediments in water.



---

*Source: 143904-2014-01-09.xml*