2. Brief review of the APSCI method
Let us assume a scene with a homogeneous circular object surrounded by a homogeneous background. Then, the scene can be modeled into two mutually exclusive regions 𝒪 and ℬ having polarimetric properties characterized by their Mueller matrices MO and MB, respectively for the object and the background. As the scene is considered to be a priori unknown, we need an initial estimation of Mueller matrices of the object M̃O and of the background M̃B by CMI before implementing APSCI method. During the Mueller imaging process, we consider that each pixel of the detector indexed by (u, v) receives an intensity I(u, v) perturbed by a Poisson distribution in order to take into account the shot noise. The Mueller matrix M̃(u, v) at each pixel is then calculated from the noisy detected intensities Ĩ(u, v).
Let us assume a totally polarized Stokes vector
S⃗ is used to illuminate the scene after CMI. The estimations of the Stokes vectors of the field scattered by the object
and background
can be expressed as :
We define the measure of separation of
and
in the Poincaré sphere by the Euclidean distance D between their last three parameters. Then, we determine numerically using a simplex search algorithm the specific incident Stokes vector S̃in that maximizes this Euclidean distance. It is worthy to emphasize that the Stokes vectors
and
are not normalized and can exhibit different rate of depolarization. As a consequence, the maximization of the Euclidean distance mentioned above takes into account two physical entities : the polarization state and the intensity of the polarized part of the scattered field.
Finally, we utilize a Two Channel Imaging (TCI) system that projects the scattered field resulting from the selective excitation
S⃗in, into 2 states of polarization
S⃗out1 and
S⃗out2, that are defined to maximize respectively (
ĨO –
ĨB) and (
ĨB –
ĨO), where
ĨO and
ĨB are the evaluations of the mean intensity detected respectively from the object and background scattering. From simple calculation it can be shown that:
where
with
The APSCI parameter is then defined for each pixel of the detector indexed by the coordinates (
u,
v) as :
where
I1(
u,
v) and
I2(
u,
v) are the detected intensity after projection respectively on the 2 states of polarization
S⃗out1 and
S⃗out2.
In this study, we use the Bhattacharyya distance as a contrast parameter for each physical quantity under investigation that can be, for comparison purposes, either the APSCI parameter as defined above, either the more pertinent parameters extracted from the polar decomposition of the Mueller matrices of the object and background. A more detailed discussion of the APSCI method is proposed in [
1M. Richert, X. Orlik, and A. De Martino, “Adapted Polarization state contrast imaging,” Opt. Express 17, 14199–14210 (2009). [CrossRef] [PubMed]
].
3. Characteristics of the Speckle noise
We have chosen to study an unfavourable situation of imaging regarding both the speckle grain size and its contrast. Thus, we assume a speckle grain with a size similar to that of the pixel of the detector. On the experimental point of view, this situation corresponds to a contrast that is not decreased by the integration of several grains into a single pixel.
Moreover, we choose to study the effect of a completely polarized and developed circular Gaussian speckle because it exhibits a strong contrast and so is susceptible to decrease the performance of the APSCI method. As will be pointed out later in this article, biological applications of the APSCI method seem very promising. So, in order to take into account some possible movement of the object, we consider a dynamical speckle: each intensity acquisition is then submitted to a different speckle pattern.
The effect of a partially polarized speckle is also of interest regarding the APSCI method because it combines 2 antagonist effects : a decrease of the speckle contrast that increases the Bhattacharyya distance of the APSCI parameter and a lower amount of polarized light usable by the APSCI method for the optimization that, on the contrary, is expected to decrease signal to noise ratio and hence this distance. In our simulations, the speckle is taken into account by a modulation of intensity at the image plane that is considered independent of the state of polarization scattered by the object and background. This modulation of intensity is performed according to the probability density function of intensity
pI(
I) of a completely developed circular Gaussian speckle that depends on the degree of polarization
P [
5J. W. Goodman, Speckle Phenomena in optics (Roberts & Company Pub., Colorado, 2007) pp. 48–50.
]:
where
Ī is the average intensity.
4. Results and analysis
We have chosen to study in
Fig. 1 the effect of a completely developed circular Gaussian speckle on three different situations where the object and background are defined to have a difference of 10% in one polarimetric property : the cases (a) and (b) exhibit this difference in the scalar birefringence, the cases (c) and (d) in the scalar dichroism and the cases (e) and (f) in the degree of linear polarization. The situations (a), (c) and (e) consider only the shot noise whereas (b), (d) and (f) take into account an additional speckle noise. For each of these situations, we calculate between the object and background region, the Bhattacharyya distance of the APSCI parameter and of the other pertinent parameters extracted from the polar decomposition. For comparison purposes, the Bhattacharyya distances are plotted versus the signal to noise ratio (SNR) for a same number of intensity acquisition. We would like to point out that, due to their different dichroism, the energy scattered by the object and background and focalized by imaging elements towards the detector can be different, and hence their corresponding classical SNR’s. Thus, we choose to define here a global SNR by considering the shot noise generated by the amount of energy received by the detector without the use of any polarizer and after the back-scattering on a virtual perfectly lambertian and non absorbing object, whose size and position are similar to that of the scene under investigation.
Fig. 1 Bhattacharyya distances obtained from different SNR levels without (case a, c, e) and with (case b, d, f) speckle noise. The object and background have a difference of 10% only in scalar birefringence, scalar dichroism and linear degree of polarization respectively for the cases a & b, c & d and e & f. In case a & b, (R, λR, ɛR) represent respectively the scalar birefringence, azimuth and ellipticity of the birefringence vector R⃗. Similarly for case c & d, (D, λD, ɛD) represent respectively the scalar dichroism, azimuth and ellipticity of the dichroism vector D⃗. For the case e & f, eigen axes of depolarization of the object and background are assumed aligned and ΔDOPL and ΔDOPC represent respectively the difference between the object and background, in their ability to depolarize linear and circular polarized light.
In
Fig. 1 B(
M) is defined to be the selected element between
M̃O and
M̃B that provides the best Bhattacharyya distance over the 16 possible elements. In a similar way,
B(
MR),
B(
MD) and
B(
MΔ) represent respectively the best Bhattacharyya distance obtained from the selected element of the birefringence, the dichroism and the depolarization matrices extracted from
M̃O and
M̃B using the forward polar decomposition.
The Bhattacharyya distances corresponding to scenes exhibiting a difference of scalar birefringence, scalar dichroism and in their ability to depolarize linear polarized light are plotted respectively as B(R), B(D) and B(DOPL). Finally, BAPSCI represents the Bhattacharyya distance of the APSCI parameter as defined in section 2.
Considering the situations (a), (c) and (e) that take into account only the shot noise, we observe that from a SNR threshold, the polar decomposition that isolates the property of interest (red and black curves) brings always better Bhattacharyya distances than
B(
M) (green curve) selected from the raw data of
M̃O and
M̃B. Below this threshold, the noise introduced by this decomposition worsen the situation (as can be seen in case (a)) because
M̃O and
M̃B are insufficiently determined. Secondly, we observe that the parameter
BAPSCI (blue curve) exhibits the highest Bhattacharyya distances for all the SNR studied in (a) (c) and (e). However, for the case (c), we notice that it exhibits also higher uncertainty bars associated to lower mean values of Bhattacharyya distances compared to cases (a) and (e). This lower performance of the AP-SCI method in the case of dichroism is coming from 2 phenomena : the absorption of energy due to the dichroism effect and the cartesian distance between the matrices of the object and background defined here as the square root of the sum of the square of the element-wise differences. Indeed, as previously discussed in [
1M. Richert, X. Orlik, and A. De Martino, “Adapted Polarization state contrast imaging,” Opt. Express 17, 14199–14210 (2009). [CrossRef] [PubMed]
], a 10% difference in one polarimetric property between the object and background gives rise to various cartesian distances in function of the scene studied. For cases (a), (c) and (e), the cartesian distances are respectively : 0.44, 0.09 and 0.14. The lowest value corresponds to the dichroism case and explains the lower performance of the APSCI method in that case.
When adding the speckle noise, we observe in (b) compared to (a), a strong degradation of all the Bhattacharyya distances under study. APSCI still remains the more pertinent parameter to distinguish the object from the background even if its standard deviation noticeably increases due to the presence of speckle. The effect of the same speckle noise on situation (c) is plotted on (d). We observe that B(MD) and B(M) fall to very low values even for high SNR and as a consequence become unusable for imaging. As a result, from the raw data of the Mueller matrices M̃O and M̃B associated to the polar decomposition, only B(D) reaches an order of magnitude similar to BAPSCI. Moreover, we notice that the standard deviation of BAPSCI has considerably increased due to the speckle noise. After a deeper analysis, we have observed that M̃O and M̃B are particularly poorly estimated for the case of dichroism (for the 2 reasons mentioned above) and that the addition of speckle noise worsen noticeably this situation. As a consequence, selective states of excitation Sin spread near all over the Poincaré sphere, showing only a weak increase of density of probability in the theoretical optimum region.
We consider now on
Fig.1 (e) and (f) the effect of a partially polarized speckle noise with degrees of polarization being respectively
Pobj =0.78 and
Pback =0.71 for the object and background that exhibit a difference of 10% in their ability to depolarize linear polarized light. We observe only a weak decrease of all the Bhattacharyya distances due to the fact that the speckle, only partially polarized, exhibits a lower contrast (
Cobj = 0.90 and
Cback = 0.87) than in previous situations. Moreover, the APSCI parameter gives rise to Bhattacharyya distances much higher than using the CMI alone or associated with the polar decomposition.
In all the previous situations, we have studied scenes that exhibit a difference between the object and background in only one polarimetric property. However, in such pure cases, due to the numerical propagation of errors, the interest of using Mueller Imaging can be inappropriate compared to simpler methods such as ellipsometry [
6P. Drude, “Über Oberflächenschichten,” Ann. der Physik 36, 86597 (1889).
][
7I. Jung, M. Vaupel, M. Pelton, R. Piner, D. A. Dikin, S. Stankovich, J. An, and R. S. Ruoff, “Characterization of Thermally Reduced Graphene Oxide by Imaging Ellipsometry,” J. Phys. Chem. C , 112(23), 8499–8506 (2008). [CrossRef]
] or polarization difference imaging methods [
8M. P. Rowe, E. N. Pugh Jr., J. S. Tyo, and N. Engheta, “Polarization-difference imaging: a biologically inspired technique for observation through scattering media,” Opt. Lett. 20, 608–610 (1995). [CrossRef] [PubMed]
][
9J. S. Tyo, M. P. Rowe, E. N. Pugh Jr., and N. Engheta, “Target detection in optically scattering media by polarization-difference imaging,” Appl. Opt. 35, 1855–1870 (1996). [CrossRef] [PubMed]
]. However, Mueller Imaging can be of great interest in the case of scenes exhibiting several polarimetric properties at the same time.
Thus, in order to examine the performance of APSCI method in such case, we consider a more complex scene where the object and background have 10% difference in scalar birefringence, scalar dichroism and in the linear degree of polarization simultaneously. In
Fig. 2, we show the Bhattacharyya distances of the APSCI parameter compared to the best Bhattacharyya distances obtained from the CMI associated to the polar decomposition that is, in this new situation, the Bhattacharyya distance corresponding to scalar dichroism. We observe that both parameters show only a weak decrease of performance (around 5%) due to the speckle noise because it is only partially polarized and exhibits a contrast significantly inferior to 1. Secondly, we see clearly that the Bhattacharyya distances of the APSCI parameter exhibits much higher values than the ones of the scalar dichroism. A visual comparison for a SNR of 3.2 is proposed in
Fig. 2 where the object can clearly be seen only using the APSCI parameter because of having 3.8 times higher Bhattacharyya distances compared to the one of the scalar dichroism.
Fig. 2 Comparison of Bhattacharyya distances vs. SNR curves for the best performing parameter of CMI (in this case the scalar dichroism) vs APSCI parameter with (ws) and without speckle noise (wos). The scene is composed of an object and a background exhibiting 10% difference in scalar birefringence R, scalar dichroism D and in the linear degree of polarisation DOPL. At the same SNR level of 3.2, the embedded images (i) and (ii) are obtained respectively using the APSCI parameter and the best parameter of the CMI.
We would like to point out that the APSCI parameter of this complex scene exhibits better Bhattacharyya distances than the pure cases of dichroism and depolarization studied separately in
Fig. 1.
However, inspite of the same amount of birefringence in the situations of
Fig. 2 and
Fig. 1(a) and that there are additional properties (dichroism and depolarization) in the mixed case which could help us to differentiate the object from the background, the Bhattacharyya distances obtained in the pure case of birefringence are higher than in the mixed case for a given SNR, even if we correct the SNR value by taking into account absorption that occurs in the latter case. This can be explained by the cartesian distances between the object and background matrices which are respectively 0.44 and 0.13 for the pure birefringent and mixed case. It can appear surprising that adding some polarimetric differences between the object and background can reduce the cartesian distance between their Mueller matrices and so our ability to distinguish them.
However, we have to keep in mind that most of the elements of a Mueller matrix describe several properties simultaneously and that from a qualitative point of view, the effects of birefringence and dichroism can produce counter-effects that will decrease the maximum distance achievable between the Stokes vector scattered by the object and background. Such observation is true for all kind of polarimetric measurements, however, as the APSCI method doesn’t consider each polarimetric property separately, and rather takes into account the whole Mueller matrices, it is expected to give always the best polarimetric contrast achievable for a 2 channel imaging system.