Reconstruction of images in non-scanned confocal microscope (NSCM) using speckle imaging

The common formation of images in CSLM assumes mechanically scanned object placed in the common short focus of the objective lenses of the microscope, while in the arrangement under study, the scanning of the object is realized by placing a diffuser behind the collimating lens. A model is suggested in the formation of images in Confocal Scanning Laser Microscope (CSLM) using non-scanned object. Since the illumination and detection are coherent, the obtained image is constructed from the simple product of the Resultant Point Spread Function (RPSF) modulated by the diffuser spread over the object transparency. Hence, the product of the object and the image of the diffuser replace the mechanical scanning of the object. Reconstructed images using this novel arrangement of CNSM are presented using mammographic X-ray image. Convolution of the RPSF and the object is realized by the spreading of the diffuser image over the object. A coherent detector captures the whole image affected by a noisy diffused function. It is noted that image processing is necessary to improve noisy images making use of filtration techniques.


Background
The ordinary confocal laser scanning microscope (CSLM) described early by Sheppard et al. [1][2][3][4][5][6][7][8][9][10][11][12] assumes mechanical scanning of the object placed in the common short focus of the objective lenses of the microscope. Coherent illumination is realized using laser beam, and coherent detection is realized by using pinhole detector to construct the image point by point where the mechanical scanning is synchronized with the electronic scanning during the detection. The detected signal is amplified and localized on CRO. Previous studies [13][14][15][16][17][18] showed an improvement in the lateral resolution using different modulated apertures. A linear, quadratic, combination of linear-quadratic [12], Gaussian and other modulated apertures is considered, while optimization of axial resolution in confocal imaging using annular pupils is investigated in [10]. Recently, a modified Hamming aperture is used in the formation of images in confocal microscope and the lateral resolution is computed from the Point Spread Function (PSF) and compared with the corresponding resolution in case of uniform apertures [19]. Confocal microscope can provide images of thick pieces of tissues which are optically sliced, instead of using microtome. This is realized because the specimen is scanned with the help of point source of laser beam. Hence, out-of-focus light is rejected, and a thin section of the tissue is obtained. Applications of confocal microscope are basically in research laboratory; however, its application in clinical settings has been also reported [20][21][22][23][24][25][26].

Open Access
In this study, we are interested to form images using non-scanned confocal microscope making using diffuser placed behind the collimating lens instead of using a grating [18].

Methods
A collimated beam from He-Ne laser is obtained using spatially filtered techniques as shown in Fig. 1. The collimated parallel beam is incident on the confocal arrangement of the microscope where a diffuser is placed behind the collimating objective lens L. In this arrangement, we assume stationary object where the mechanical scanning is replaced by the image of the diffuser. Hence, the object is covered by the image of the diffuser convoluted by the Point Spread Function corresponding to the 1st objective lens. The following steps for the formation of images using non-scanned object are summarized as follows: 1. Consider unit amplitude of coherent radiation incident upon the diffuser placed in contact with the collimating lens L. The diffuser is represented by a randomly distributed function d (x 1 , y 1 ) and obstructed by the 1st aperture P (x 1 , y 1 ). Then the transmitted complex amplitude is represented as follows: ρ: is the radial coordinate in the aperture plane of the collimating lens L.
(1) A x 1 , y 1 = d x 1 , y 1 · P x 1 , y 1 where P x 1 , y 1 = 1; for ρ ρ 0 ≤ 1 = zero; otherwise 2. In the back focal plane of the collimating lens L, we get by applying the F.T. upon Eq. (1), the following: 3. For simplicity, assume the F.T. corresponding to the collimating lens is replaced by a Dirac-delta function where the illumination is considered coherent. Hence, Eq. (2) is approximately given as: 4. In the focal plane of the 1st objective lens limited by the aperture P 1 (u, v) where the transparency of the object is placed in (x, y) plane, we apply the F.T −1 . upon Eq. (3) to get the following: 5. The 2nd objective is shown conjugate to the 1st objective where the pupil aperture P 2 (u, v) is shown in front of the 2nd objective. The Point Spread Function is equal to that shown for the 1st objective is formed in the common short focus where the object is located. Consequently, the object transparency g (x, y) is multiplied by both the PSF's but convoluted with the diffuser covering the non-scanned object. Hence, we get the following multiplication in the object plane of transmitted complex amplitude C (x, y): And d x 1 , y 1 = rand x 1 , y 1 .
(2) 6. The detected intensity is represented as the modulus square of the complex amplitude C x, y , represented as follows: 7. Since the convolution product of the diffuser and the PSF corresponding to the 1st objective lens of aperture P 1 (x, y) gives truncated image of the diffuser convoluted by the 1st objective PSF formed in the object plane (x, y), then we can write the complex amplitude corresponding to the modulated diffuser image as follows: Substituting from (7) in (6), the detected intensity is rewritten as follows: Or using equation (6), we get the following convolution: Consequently, the detected image intensity is simply affected by the PSF corresponding to the 2nd objective, while the PSF corresponding to the 1st objective is convoluted with the diffuser image giving the modulated diffuser. Hence, the whole image is formed without the object scanning since it is integrated by the modulated diffuser. The Resultant Point Spread Function (RPSF) is computed from the product of the PSF corresponding to each objective and written as follows: 8. Finally, the coherent extended detector is required to capture the whole image covered by the randomly distributed function. 9. Image processing is necessary to extract an improved image.

Results
The input mammographic image used in the processing is shown in Fig. 2.
We construct random diffuser d (u, v), of dimensions 512 × 512 pixels as shown in Fig. 3. It is placed in front of the collimating lens L. This diffuser incident upon a circular aperture of diameter 128 pixels is given in Fig. 4. The PSF corresponding to the aperture of the collimating h r x, y = h 1 x, y h 2 x, y lens convoluted with the Fourier spectrum of the diffuser is shown in Fig. 5. It is called modulated speckle pattern since the F.T. of the diffuser is named ordinary speckle assuming high numerical aperture (NA).
The modulated speckle pattern shown in Fig. 5 is truncated by the 1st microscope objective is shown in Fig. 6.
The FT of the modulated speckle multiplied by the transmittance from 1st objective of the microscope will give the image of the diffuser convoluted by the PSF of the 1st objective located in the object plane as shown in Fig. 7. The modulus square of the multiplication of the above convolution with the object and the PSF corresponding to the 2nd objective forming the detected noisy image is shown in Fig. 8a. The detected images in absence of the diffuser are shown as in Fig. 8b for aperture radius = 128 pixels, while Fig. 8c shows the image for aperture radius = 16 pixels. The figures from 2 up to 8 represent images of two-dimensional matrix of 512 × 512 pixels starting from (0,0) in the upper left corner and ending with (512, 512) in the lower right corner.

Discussion
The proposal of NSCM using a diffuser gives reconstructed images in the detection plane affected by the diffuser noise. The reconstructed image may be improved by filtration technique. This new arrangement of NSCM is compared with the ordinary confocal microscope provided with the mechanical scanning of the object placed in the confocal plane of the objectives. It is shown equal resolution in both cases with and without diffuser. The different confocal microscope arrangements gave equal values of resolution since the Resultant Point Spread Function (RPSF) is only dependent on the objectives and the laser beam of wavelength λ. It is computed from the PSF corresponding to each objective lens represented by Eq. (10). The PSF corresponding to uniform circular aperture has the known Airy disc [2]. It is shown a noisy image affected by the diffuser as in Fig. 8a, of resolution dependent on the PSF corresponding The modulated speckle pattern which is the FT of the multiplication of the diffuser and the circular aperture. It is located before the 1st microscope objective Fig. 6 The modulated speckle pattern shown in Fig. 5 is truncated by the 1st microscope objective. The aperture radius = 128 pixels Fig. 7 The FT corresponding to the multiplication of the modulated speckle shown and the 1st objective of the microscope shown in Fig. 6. The formed pattern which is the diffuser image convoluted with the PSF corresponding to the 1st objective located in the object plane Page 5 of 7 Hamed and Al-Saeed Beni-Suef Univ J Basic Appl Sci (2021) 10:67 to the objectives. The detected intensity for greater radius (128 pixels) shown in Fig. 8b has better resolution than that shown in Fig. 8c for smaller radius = 16 pixels. This is attributed to the inverse relation between resolution and the aperture radius for certain focal length which is determined from the spatial radial cutoff in the focal plane. The PSF using uniform circular aperture versus radial distance in the Fourier plane is plotted as in Fig. 9. The cutoff values computed to represent resolution are as follows: The cutoff radial distance (r c ) = 0.5039 µm for aperture radius = 16 pixels, while (r c ) = 0.06 µm for aperture radius = 128 pixels. Another values of cutoff are given as follows: (r c ) = 0.2581 µm for aperture radius = 32 pixels, and (r c ) = 0.1352 µm for aperture radius = 64 pixels. It is shown inverse relation between the cut-off value r c and the aperture radius as expected from the resolution limit, assuming monochromatic light for the microscope illumination of wavelength λ. The theoretical resolution limit is given by: resolution limit = NA .
Consequently, the highly resolved images are obtained for sharper PSF hence lower cutoff value as shown in Fig. 8b.

Conclusions
A confocal microscope based on substituting the mechanical scanning of the object by the image of the diffuser formed in the object plane. The diffuser is placed before the collimating lens of the spatial filter. Hence, speckle pattern is formed in the back focal plane of the collimating lens which is the Fourier transform of the transmitted diffuser. Then, operating the inverse Fourier transform upon the 1st objective lens limited by the aperture P 1 to get the convolution of the diffuser image and the PSF corresponding to the 1st objective. This convolution product is formed in the object plane. Consequently, the detected intensity is computed from the multiplication of the object with the diffuser convoluted with the PSF corresponding to both microscope objectives. The detected image is affected by a noise originated from the  Fig. 8 a The detected intensity originated from the multiplication of the modulated diffuser with the object transparency affected by the PSF corresponding to both objectives. b The detected intensity obtained in absence of the diffuser using mechanical scanning of the object synchronized with the electronic scanning in the detection plane. The microscope used is called confocal laser scanning microscope (CLSM). The aperture radius = 128 pixels. c The detected intensity obtained in absence of the diffuser. The microscope used is called confocal laser scanning microscope (CLSM). The aperture radius = 16 pixels radial coordinate in microns The PSF using uniform circular aperture of radius = 16 pixels. The cutoff radial distance (r c ) = 0.5039 µm as shown in the curve