Influence of multiple scattering on the resolution of an imaging system: a Cramér-Rao analysis

We revisit the notion of resolution of an imaging system in th e light of a probabilistic concept, the Cram ér-Rao bound (CRB). We show that the CRB provides a simple quantitative estimation of th e accuracy one can expect in measuring an unknown parameter from a scatt ering experiment. We then investigate the influence of multiple sc attering on the CRB for the estimation of the interdistance between two obje cts in a typical two-sphere scattering experiments. We show that, contrari ly to a common belief, the occurence of strong multiple scattering does no t automatically lead to a resolution enhancement. © 2006 Optical Society of America OCIS codes:290.4210 Multiple Scattering; 290.4020 Mie Theory; 110.01 8 Microscopy References and links 1. F. Simonetti, ”Multiple scattering: The key to unravel the subwavelength world from the far-field pattern of a scattered wave,” Phys. Rev. E 73, (2006). 2. F.C. Chen and W.C. Chew, ”Experimental verification of supe r resolution in nonlinear inverse scattering,” Appl. Phys. Lett.72, 3080–3082 (1998). 3. Belkebir K, Sentenac A, Chaumet PC, ”Influence of multiple sc attering on three-dimensional imaging with optical diffraction tomography,” J. Opt. Soc. Am. A 23, 586–595 (2006). 4. P. Blomgren and G. Papanicolaou, ”Super-resolution in time -rev rsal acoustics,” J. Acoust. Soc. Am. 111, 230– 248 (2002). 5. C. Prada and J.L. Thomas, ”Experimental subwavelength loca lization of scatterers by decomposition of timereversal operator,” J. Acoust. Soc. Am. 114, 235–243 (2003). 6. M. Shahram and P. Milanfar, ”Imaging below the diffraction limit: A statistical analysis,” IEEE Trans. Image Proc.13, 677–689 (2004). 7. S. Van Aert, D. Van Dyck, and A.J. den Dekker, ”Resolution o f coherent and incoherent imaging systems reconsidered classical criteria and a statistical alternative ,” Opt. Express14, 3830–3839 (2006). 8. P. Refregier, “Noise theory and application to physics fr om fluctuation to information,” chapter statistical estimation, 167–205, springer 2004. 9. Y.L. Xu, ”Electromagnetic scattering by an aggregate of sp heres,” Appl. Opt.34, 4573–4588 (1995). 10. P.C. Chaumet and M. Nieto-Vesperinas, ”Optical binding o f particles with or without the presence of a flat dielectric surface,” Phys. Rev. B 64, 035422–035429 (2001).


Introduction
Many quantitative imaging systems, in the optical, micro-wave or acoustical domains, resort to numerical inversion procedures to reconstruct the image of the sample from its diffracted far-field. Several authors have pointed out that, accounting for multiple scattering in the inversion model, could ameliorate the image resolution beyond that expected with the Rayleigh criterion λ /(2NA), where λ is the incident radiation and NA is the numerical aperture of the system [1,2,3]. Similarly, it has been observed, in the context of time reversal focusing that, the spot width of the time reversed wave is smaller when the wave propagates in a random medium, underlying multiple scattering, than when it propagates in a homogeneous medium [4,5]. Hence, the question whether multiple scattering within the sample may lead to a better resolution than that expected with a single-scattering analysis stirs considerable interest. The usual argument in favor of this thesis is that multiple scattering permits the conversion of evanescent waves into propagative ones so that the diffracted far-field conveys information on the sample's high spatial frequencies which, hopefully, can be recovered with an adequate inversion procedure [1,3]. Yet, to our knowledge, sole qualitative studies on this subject have been conducted and a quantitative analysis of the influence of multiple scattering on the resolution has never been done. To address this problem, the first step consists in clarifying the notion of resolution. Indeed, in quantitative imaging, the resolution depends strongly on the chosen physical model, the a priori information therein, and the signal to noise ratio (SNR). In this case, the Rayleigh criterion, which depends solely on the configuration of the set-up, is not sufficient and a statistical criterion accounting for the SNR and the fitting physical model appears more adequate [6,7].
In this paper, we quantify the influence of multiple scattering on the resolution of an imaging system by means of the Cramér-Rao bound (CRB). This signal processing notion permits one to evaluate the lower bound of the variance of the parameters retrieved by an inversion procedure from noisy experimental data. The smaller the variance, the more accurate the estimation of the parameter. More precisely, our experiment consists in illuminating two small spheres by a plane wave and collecting their diffracted intensity in far-field with a certain number of detectors within a certain angular cone. To simplify the analysis to its core, we limit the imaging problem to the sole estimation of the interdistance α between the spheres centers. The variance of the estimateα obtained from the data of the diffracted intensity is calculated for several values of α, radius and permittivity of the spheres, the latter being chosen to generate more or less interaction between the spheres.

The Cramér-Rao bound
The Cramér-Rao inequality provides a universal lower bound for the variance of an arbitrary unbiased estimator (e.g. Ref. [8]). This bound is optimal in the sense that the equality can be reached in extreme cases. We will recall the main result of the Cramér-Rao inequality in the context of an optical imaging system. We call I m = [I m 1 , ..I m N ] the vector of the N measured far-field intensities at different angles in one realization of the experiment. We denote α the unknown parameter of the scattering object to be estimated. In our experiment α is the centerto-center distance of two homogeneous identical dielectric spheres of given diameter d and optical index n. The Cramér-Rao analysis also applies to the multi-parameter case, but for sake of simplicity we will restrict ourselves to a single parameter. Now, the measurement of the scattered intensities I m has limited accuracy since it is affected by various sources of noise. We denote I t = [I t 1 , ..I t N ] the "true" values of the intensity, that is the values that would be measured in a noiseless experimental system. An estimatorα for the unknown parameter can only be derived on the basis of the noisy data. Hence, the variance of this estimator depends drastically on the level of noise in the measurement. The Cramér-Rao inequality sets an absolute lower bound for this variance in the case of an unbiased estimator: where E is the mathematical expectation over the realizations of noise. We have introduced the log-likelihood function L(α), which is the logarithm of the density of probability of measuring the values I m for a given value of the parameter α,and its derivative ∂ α L with respect to the parameter α. We see from Eq. (1) that the CRB is a measure of the sensitivity of the experiment to a parameter. Indeed, a sharp variation of the detected intensity, as the parameter α varies, results in large derivative of the log-likelihood function and a small CRB. More explicit calculation requires a model of noise. We will adopt some reasonable assumptions that render the analysis tractable.
Gaussian additive noise The source-independent background noise is often modelled as a Gaussian additive noise. We therefore write the measures I m j as: where N j are Gaussian, centered, and perfectly uncorrelated identical variables: where I 0 is a constant intensity and δ i j the Kronecker symbol. With these hypotheses the CRB can be easily derived: Multiplicative noise Another cause of noise in an optical or micro-wave experiment is the fluctuation of the source. This is in general modelled by a Gamma distributed multiplicative noise: The variables N j will again be assumed perfectly uncorrelated. They will be taken to follow a Gamma law of order L and mean µ: Direct calculations lead to a similar expression as in the additive case: Note that in the case of a multiplicative noise, the SNR is kept constant in all configurations. The CRB is the basis of our definition of the resolution of an imaging system. Precisely, we define the resolution on parameter α as the square root of the relative minimal variance predicted by the Cramér-Rao bound: This quantity gives an estimation of the relative accuracy one can expect from a parameter estimation based on a scattering experiment. With this quantitative tool we will now study the influence of multiple scattering on the resolution of an optical imaging experiment.

A scattering experiment
The experimental set-up is as follows. Two homogeneous dielectric spheres of diameter d and optical index n are placed in vacuum. They are separated by a distance α and aligned along the x-or z-axis. An incoming plane wave with wave vector K 0 = 2π/λ z illuminates the spheres. The plane of observation is ( x, z) and the scattering angle θ j defines the direction of observation K j = cos θ j z − sin θ j x. This incident wave is polarized along x or y. The scattered intensity is recorded in a 30 degrees aperture cone around the forward or around the backward direction.
The different configurations with their nomenclatures are depicted on Fig. 1. These configurations have been chosen to enhance or disminish the role of single scattering in the estimation of α. Fig. 1. The set-up of the scattering experiment. In the first configuration, the spheres are aligned in the direction of the incident wave vector K 0 . The scattered waves (K) are recorded in a 30 degree aperture cone in the scattering plane ( z, x) around the forward (A) or backward (B) direction. The incident polarization E 0 is perpendicular to the scattering plane. In the second configuration, the sphere alignment is perpendicular to K 0 and only the backward direction is investigated. The polarization E 0 is perpendicular (C1) or parallel (C2) to the scattering plane.

Single scattering analysis
When the interaction between the spheres is neglected, the scattered intensity in the K direction, I( K), is given by the classical interference formula: where I s ( K) is the scattering intensity of one single sphere, Φ( K) = 2π λ ( z − K) · u and u is the direction of alignment of the spheres. Assuming a large number of angular observations, so that the discrete summation can be replaced by a continous one, we obtain a simple analytical expression for the CRB in the case of a multiplicative noise: where the proportionality constant is independent of the configuration and the bounds of the integral are the extreme values assumed by the phase Φ in the inspected angular interval. The remarkable feature of this expression is that it does not depend on the shape and size of the scattering objects, since the one-particle intensity I s is cancelled out in the calculation. A similar formula can be derived in the additive case, assuming in addition that I s is quasi-constant in the angular region of observation (this is actually the case for small spheres): At this stage we can draw several conclusions from the single scattering analysis of the CRB. In the case of an additive noise, the resolution is trivially ameliorated as the scattering power I s of one single sphere increases. This not surprising since it corresponds to an increase of the SNR. For small spheres and small contrasts it is well known that the scattering power grows like I s ∼ (n 2 − 1)/(n 2 + 2) 2 as the optical index is increased. Hence, doubling the optical constrast n 2 − 1 will result in an almost four times smaller CRB. In the case of a multiplicative noise, the resolution is not related to the scattering power of the spheres.
This simple study shows that a fair comparison of the resolution of different scattering configurations requires a constant SNR, which amounts to consider multiplicative noise only. To mimick the classical experiment of point source imaging, we consider two spheres whose diameter is much smaller than the electromagnetic wavelength (d = 0.06λ ). Figure 2 shows the evolution of the inverse relative CRB as a function of the normalized interdistance α/λ for the three scattering configurations depicted on Fig. 1 (C 1 and C 2 give almost identical results) and the two types of noise. It appears from these plots that the configurations B,C and A listed in this order correspond to an increasing degree of resolution. This is in agreement with the widely accepted rule that finer details of the scatterer can be obtained as the Ewald vector K − K 0 is increased, since this corresponds to the spatial frequency probed by the interrogating wave in the Rayleigh or Born regime. Indeed we recall that: where χ is the Fourier transform of the permittivity function of the object. Note that the peaks observed for the CRB curves in the B configuration are basically an artefact of the multiplicative noise. Indeed, they correspond to interdistances α = λ /4 + pλ /2, where p is an integer, for which the scattered intensity in the backward z direction is null due to destructive interference. Thus, the noise is also null for this direction.
We will now consider the following question: at a given SNR, is a configuration with strong multiple scattering better resolved than a configuration with single scattering ? To answer this question, we consider the same diffraction experiments as previously but we now account for the interaction between the spheres. Since we now cannot calculate the CRB analytically, we use a rigorous method [9] to simulate the intensity scattered by the two spheres and we build numerically the CRB.

Multiple scattering analysis
In Fig. 3 the CRB is plotted for multiplicative noise in the C configuration for the three different values of the sphere optical index. The multiplicative constant µ 2 /L has been set to 1. We consider two incident polarizations, orthogonal (polarization 1, (3a)) or parallel (polarization 2, (3b)) to the plane of observation. Not surprisingly, we observe that as α/λ is increased, the CRB obtained in the rigorous experiments rejoins that obtained when the interaction between the spheres is neglected. On the contrary, when the coupling between the spheres becomes important, the CRB gets lower than that given by the single scattering analysis. We recall that, for a given α, the electromagnetic interaction between the spheres is more important in polarization 1 than in polarization 2 [10] and that it increases with the optical index of the spheres. This analysis clearly shows that, in this configuration, the coupling between the spheres can ameliorate the resolution of the image. A similar conclusion is obtained with the analysis of the CRB in the A configuration as seen in Fig. (4b). On the other hand, the analysis of the CRB in the B configuration leads to a quite different result. Indeed, we observe in Fig. (4a) that, in this case, the interaction between the spheres does not have any effect on the CRB. Hence, we see from these different exemples that the potential amelioration of the resolution due to multiple scattering is only obtained when the experimental configuration is incomplete, i.e. when the Ewald vector scanned in the experiment does not reach the highest possible spatial frequencies.
To increase further the coupling strength between the objects, we have conducted the same experiments as in Fig. 4 with spheres of larger diameter (d = 0.3λ ). The analysis of the CRB in this case (Fig. 5) stresses even more the influence of the configurations. Indeed, in configuration A (Fig. 5b) the amelioration of the resolution due to multiple scattering is patent. Moreover, we notice by comparing Fig. (4b) and (5b) that the estimation of α can be more accurate for larger spheres than for smaller ones. This result confirms the study of Belkebir et al [3] in which it was shown that, with a non-linear inversion algorithm and an experimental set-up close to the A configuration, two big cubes will be distinguished more easily than two small cubes with the same center to center distance. On the other hand, it is seen that, in the B configuration (Fig. (5a)), multiple scattering does not permit an amelioration of the resolution and can even deteriorate it.

Conclusion
We have proposed a statistical analysis of the resolution of a non-linear quantitative imaging system accounting for the noise and the fitting model. This approach has allowed us to point out the role of multiple scattering in the resolution of the system. More precisely, we have studied different angular configurations of an optical diffraction tomography experiment plagued with additive or multiplicative noise. We have studied the accuracy of the estimation, from the scattered intensities, of the center to center distance of two spheres. We have shown that when the Ewald vectors, (or momentum transfer) scanned in the experiment are small, multiple scattering permits a clear improvement of the interdistance estimate. On the other hand, if the Ewald vectors are large, multiple scattering does not modify the accuracy with which the interdistance is evaluated. This analysis is in agreement with previous studies in which an amelioration of the resolution, due to multiple scattering, was observed in incomplete experimental configurations.  Fig. 4 for spheres of diameter 0.3λ . The coupling between the spheres is visible on a wider range as for small spheres. In (a), the peaks occur, again, at the distances that lead to destructive interference between the spheres.