References

During the past decade, optical imaging through scattering medium has proved to be a powerful technique for many applications. It is especially effective in medical diagnostic, since it is safe, noninvasive and low-cost compared with the conventional radiation techniques. Based on a similar principle of the fly's visual system, we show a novel method of optical imaging through scattering medium. An image of bones hidden between two biological tissues (chicken breast) is recovered from many noisy speckle pictures obtained on the output of a multi-channeled optical imaging system. The operation of multiple imaging is achieved using a microlens array. Each lens from the array projects a different speckled image on a digital camera. The set of speckled images from the entire array are first shifted to a common center and then accumulated to a single average picture in which the concealed object is exposed.

The compound eye is the most common eye in the nature [1], but species of the fly have special kind of compound eyes equipped with visual system called neural superposition.Neural superposition means that several (seven according to Ref. [2]) identical images from several facets are superposed together to a single common image by neural connections of the photoreceptors.Land and Nilsson [2] explain that the advantage of the neural superposition is in the gain of the observed image which enables the flies to see better than their competitors under faint illumination conditions.The question whether it is the only advantage is still open, but the problem we consider here is whether any engineering visual system following the principle of neural superposition might have any useful purpose.Obviously, the nowadays digital cameras have more efficient ways to gain a single recorded image without capturing multiple images of the same scene.This study shows that by using the fly eye principle, i.e., superposing multiple images from many imaging channels, enables one to see general objects hidden behind scattering layers.In recent years much effort has been devoted to research in the optical imaging of objects embedded in a scattering medium, and different optical imaging techniques have been proposed [3].This topic has many potential applications in medical diagnostics since it is safe, noninvasive, and relatively inexpensive compared with other often-used tomography techniques.
In this work we propose a new scheme of seeing through scattering medium which is similar to the fly's visual system.This system averages over an ensemble of noisy speckle images and therefore we term the system as Noninvasive Optical Imaging by Speckle Ensemble (NOISE).In order to collect different blurred images of the same object, the object should be observed from unrelated spatial parts of the same scattering medium.This goal is achieved by use of microlens array (MLA), where each lens from the array images the same object through different part of the scattering layer.In the experimental system, shown in Fig. 1, the object hidden between two scattering tissues is coherently illuminated from its back by a plane wave.Each lens of the array, together with the spherical lens L, operate as two successive imaging systems.Without the scattering layers, the illuminated coherently system is characterized by a relatively narrow point spread function (PSF) [4] h o (x,y).This PSF is calculated conventionally as an inverse Fourier transform of the aperture of a single microlens.In the present case, the micro-lens imposes the system bandwidth because its numerical aperture is smaller than that of the lens L.
Next, we consider the effect of the back scattering layer S 1 .This layer diffuses the light such that each micro-lens gets almost uniformly part of the illumination.In addition, because of the randomness of the medium S 1 and its uniformity, the object is multiplied by a random phase function with almost constant magnitude.The entire system is modeled as an array of several identical imaging systems, all with the same PSF given by h o (r), where r=(x,y) is the When the front scattering layer S 2 is introduced into the system, the output image is distorted such that the object cannot be recognized.Since each micro-lens observes the object through a different transverse cross-section of the scattering layer, each k-th micro-lens together with the lens L create a linear system characterized by a different random PSF h k (r).Therefore, the output intensity pattern in each coherently illuminated k-th channel is given by 2 It is assumed that although each PSF h k (r) is a random function, wider than h o (r), the ensemble average PSF over the entire K channels satisfies the relation: This assumption is valid for scattering mediums satisfying statistics of Rytov model with weak phase modulation and Born model [5].Having the set of K speckled images { } ( ) we first center each one of them and then sum them to a single average image given by, ( ) ( ) ( ) ( ) Based on Eq. ( 2), it can be shown that the average image is approximately, where σ 2 is the variance of the random set {h k (r)}, defined as The first term of Eq. ( 4) is approximately the desired diffraction-limited image given by Eq. ( 1).The second term is a convolution between the object and the variance functions.Clearly σ 2 (r) is wider than h o (r) because the scattering layer broaden the diffraction-limited image of a point.Therefore, we conclude that the second convolution in Eq. ( 4) blurs the diffraction-limited image of the object.The value of this blurring term is determined by the average value of the variance σ 2 (r).The contrast and the sharpness of the reconstructed object are inversely dependent on the variance.
To demonstrate the proposed technique, an opaque object made of chicken bones in a shape of a cross-junction with the size of 9×9mm was embedded between two layers of Fig. 2. The average picture of the entire array when the cross-junction of chicken bones is positioned in front of layer S 1 and the layer S 2 is removed.
chicken breast separated from each other by a distance of 12mm.Fig. 2 shows the average picture of the exposed cross-junction of bones over the entire array when only layer S 1 exists in the set up, and S 2 is removed.The thickness of the rear tissue S 1 was about 3mm, whereas the thickness of the front tissue S 2 , was about 8mm.The reduced scattering coefficient of the tissues of µ s '=4.5±0.3cm - was measured by the method proposed in Ref. [6].Assuming the anisotropy factor [7] is g=0.965,The scattering coefficient becomes . The rear tissue S 1 was illuminated by a collimated plane wave of 35mW He-Ne laser with λ=632.8nmwavelength.The MLA, placed a distance of b=160mm from the object, was composed of 115×100 hexagonal refractive lenses.A similar MLA has been used in integral imaging systems for purpose of three-dimensional imaging [8].Only the central 132=12×11 lenses were used in the present experiments.The radius of each microlens was r l =250µm and its focal length was 3.3mm.Under these conditions the optical system without the tissues can resolve a minimum size of λb/r l ≅0.4mm.The resolution can be improved using lenses with larger apertures.However, this change may increase the total view angle of the system and thus the various channels may image different perspectives with different shapes of the same object.The image plane of the MLA was projected onto the CCD plane by a single spherical lens L, with a 300mm focal length.The distance Z 1 and Z 2 shown in Fig. 1 were 520mm and 710mm, respectively.In the experiment we used a CCD camera with 1280(H)×1024(V) pixels, within 8.6×6.9mmsquare active area.
Figure 3 shows the array of all speckled images recorded by the CCD, whereas the embedded object was the cross-junction of bones.The white lines indicate the image area improvements in the image happen during the summation of the first 10 pictures, and adding more pictures beyond 70 do not improve much the final image.These observations are quantitatively demonstrated by the plot of the average mean square error (MSE) between the final image, averaged over 132 images, and every partial summation from 1 to 132 (the blue line of Fig. 5).Indeed the improvement after the first 65 pictures is relatively small.The red line of Fig. 5 shows the MSE between the non-scattered object of Fig. 2 and the images obtained from N additions.This plot indicates that after 85 additions the MSE reduction is negligible.These results indicate that we could reduce the number of channels by about half and gain more resolution, or alternatively we could hide the object behind a thicker scattering layer and still distinguish the object.
To verify that it is necessary to use coherent light, we imaged through the same scattering layers with incoherent light.The picture obtained in each channel is a wide blurry spot.All the pictures of the set seem to be more or less the same wide smoothly blurry spots.Averaging over all 132 pictures yields again a wide smoothly blurry spot.Incoherent light means that the object is illuminated by a large number of plane waves with many different angles.The effect on the output image in each channel is an accumulation of many blurred images of the object shifted randomly from the true object center.Therefore the result in each imaging channel is a smoothly blurred unrecognizable image of the object.Accumulating these images along all the channels does not allow us to see through the scattering medium.
In conclusion, by the NOISE technique we have been able to reconstruct the shape of objects embedded between two scattering layers.The weakness of the present setup is the relatively low spatial bandwidth product of the diffraction limited system.The use of a small aperture lens at each imaging channel reduces both the field of view and the system's bandwidth.However, this drawback seems as a reasonable penalty to pay for the ability to see through scattering medium in a simple and robust way.The advantages of the method are relative simplicity, low cost, fast operation and the need of low power CW laser illumination.Because of all these advantages NOISE might be useful for many imaging applications, especially in the medical diagnostic.

Fig. 3 .
Fig. 3. 12×11 blurred images recorded by the CCD when the cross-junction sign is embedded between the two scattering layers.

Fig. 5 .
Fig. 5. Average Mean Square Error versus the number of additions, where N goes from 1 to 132.Blue line -MSE between the final image and the images obtained from N additions.Red line -MSE between the images of Fig. 2 and the image obtained from N additions.