Next Article in Journal
A Robust Brain MRI Segmentation and Bias Field Correction Method Integrating Local Contextual Information into a Clustering Model
Previous Article in Journal
A CNN Model for Human Parsing Based on Capacity Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electro-Optical Imaging Technology Based on Microlens Array and Fiber Interferometer

1
University of Chinese Academy of Sciences, Beijing 100049, China
2
Institute of Optics and Electronics, Chinese Academy of Science, Chengdu 610209, Sichuan, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(7), 1331; https://doi.org/10.3390/app9071331
Submission received: 2 March 2019 / Revised: 24 March 2019 / Accepted: 26 March 2019 / Published: 29 March 2019
(This article belongs to the Section Optics and Lasers)

Abstract

:
To reduce the volume and weight of traditional optical telescopes effectively, this article proposes an electro-optical imaging technology based on a microlens array and fiber interferometer. Pairs of microlenses in the microlens array collect light and couple it into a fiber interferometer to form interference fringes. Then the amplitude and phase of a large number of interferometer baselines are analyzed to generate images. In this work, the principle of electro-optical imaging technology has been analyzed according to the partially coherent light theory. The microlens-array arrangement method and baseline pairing method have been optimized for arbitrary targets. From the simulation results, it was found that the imaging resolution depends on the maximum baseline length, and the imaging quality could be effectively improved by adjusting the Nyquist sampling density and baseline pairing method. This technology can provide an important reference for the miniaturization and complanation of imaging systems.

1. Introduction

Without considering the influence of astronomical observation, the imaging resolution of traditional optical telescopes based on the principle of refraction can only be improved by increasing the aperture. However, the manufacture of large-aperture optical telescopes is an enormously expensive undertaking and also requires significant technical support. To overcome the limitations of the image resolution by the size of a single aperture, many novel optical telescopes based on interference imaging have been designed in recent years. For example, the Cambridge Optical Aperture Synthesis Telescope (COAST) has achieved imaging [1], and J. A. Benson et al. completed the Naval Prototype Optical Interferometer (NPOI) optical synthetic aperture telescope array and successfully performed observation and image reconstruction [2,3]. In addition, several cutting-edge electro-optical imaging systems have been designed to greatly reduce the volume and weight of traditional optical telescopes, like the geosynchronous earth orbit (GEO) telescopes, which are an interferometric imaging system with arrays of telescopes. Other approaches include that of combining observations from a separated-element interferometer with interferometric data obtained by optical masking of a single-dish telescope [4]. The most representative cutting-edge electro-optical imaging system is the Segmented Planar Imaging Detector for Electro-Optical Reconnaissance (SPIDER), which has been developed by Lockheed Martin and can realize a 10× to 100× reduction in size and weight by replacing the traditional lens stacks with silicon chip photonic integrated circuits (PICs) [5,6,7].
Most of the above telescopes are based on astronomical light interference and optical synthetic aperture imaging [8,9,10]. The optical synthetic aperture imaging technique is to improve the resolution of single aperture imaging by measuring the interference result between the observation plane aperture pairs (ie, the spectral value of the observed target). In this technique, the interference fringe of a pair of apertures corresponds to a spectral value of the target, so the complete spectrum of the target is often obtained by changing the relative position between the pairs of apertures or increasing the number of apertures. As we all know, the optical fiber optical path has the characteristics of softness, long transmission distance and strong anti-interference. The use of optical fiber to realize light interference is an important application of interference phenomenon, and will greatly reduce the working distance of the interference system [11,12].
To reduce the volume and weight of traditional optical telescopes effectively, we propose in this paper a cutting-edge electro-optical imaging technique based on a microlens array and fiber interference. In this technique, the light from a scene was first collected by the microlens array and then coupled into the fiber interferometer. Note that each microlens pair was coupled by fiber to form an interference baseline, which had a significant effect on the imaging resolution. Thus, through the optimal microlens arrangement and pairing method, we could complete spectrum sampling of the scene and generate high-resolution images. To establish the optical simulation model, the imaging principle based on the partial coherence theory was first studied first. Then, a complex rectangular sub-aperture arrangement and a symmetrical baseline pairing method were proposed to improve the imaging resolution and quality further. Moreover, we investigated the impact of spatial frequency distribution (u, v) on imaging quality, including the resolution, field of view, and depth of focus of the imaging system. The simulation results provide a valid reference for the design of the actual system.

2. Structure Design

The structure of the electro-optical imaging system, which consists of a microlens array, fiber interferometer, and a processing module, is shown in Figure 1. The light information is collected and coupled into the fiber interferometer by the microlens array. Then, the interference fringes are detected and analyzed by the processing module. Finally, a computer generates an image based on the amplitude and phase information extracted from the interference fringes. Each fiber interferometer includes a pair of microlenses for measuring the spectrum of one baseline.
The spectral distribution measurement was determined by the arrangement and pairing of the microlenses. As shown in Figure 1a, the rectangular microlens array uses a centrally symmetric baseline pairing method (e.g., pentagram and pentagram, diamond and diamond, etc.). We also proposed a complex rectangular sub-aperture arrangement method to optimize the microlens-array arrangement (Figure 1b). The complex rectangular sub-aperture arrangement achieved different sampling intervals at high, intermediate, and low frequencies.
The performance of the fiber interferometer is closely related to the quality of interference fringes, whose core device is a 3-dB fiber coupler. Because achieving a fiber coupler coupling efficiency of 100% is difficult, the split ratio cannot be 1:1 [13]. Interference arm length differences of the fiber optic interferometer will also introduce additional phase differences, which will eventually affect the true reflection of the stripe contrast and phase.

3. Theory and Algorithm

Figure 2 shows the principle and processes of the electro-optical imaging technology; they consist of four parts: interference between microlens pairs, calculation of complex degree of coherence according to fringes, spectral coverage, and image reconstruction.
Using the Van Cittert–Zernike theorem in partial coherence theory [8], this study analyzes the mutual intensity of a light field illuminated by an incoherent source and reaches an important conclusion: when the linearity of the source and observation area are much smaller than the distance between them, the complex coherence coefficient over the observation area is proportional to the normalized Fourier transform of the source intensity. According to the principle shown in Figure 2, the complex coherence coefficient between Q1 and Q2 in aperture planes can be expressed as [14]:
μ ( Q 1 , Q 2 ) = + o ( α , β ) e x p [ j 2 π λ z ( ( x 2 x 1 ) α + ( y 2 y 1 ) β ) ] d α d β + o ( α , β ) d α d β
{ u = Δ x λ z = ( x 2 x 1 ) λ z v = Δ y λ z = ( y 2 y 1 ) λ z
According to Equation (2), the distance between the sub-apertures Q1 and Q2 in a pair is called the interferometer baseline, which is closely related to the spatial frequency ( u , v ) of the detection scene. Therefore, we can optimally arrange the microlens array and use a better baseline-pairing method to realize spatial frequency coverage. Equations (1) and (2) provide a theoretical basis for interference imaging, which is different from the aberration imaging theory and the diffraction imaging principle.
The first step in reconstructing the image is to solve the complex coherence coefficient of a pair of lenses. According to the partial coherence theory, the complex coherence coefficient of two apertures on the observation plane can be determined by measuring the interference fringes. The second step is to change the interference baseline and repeat the first step to achieve spectral coverage. Finally, the solution result is substituted into the Equations (1) and (2), and the inverse intensity of the target scene is obtained by using the inverse Fourier transform. The detailed process and the simplified version of the program are shown in Algorithm 1.
Algorithm 1. Imaging process and reconstruction algorithm. Solving fringes and image reconstruction
1:  Initial(obj,aperture, …);
2:  for i,j = 1: (size(aperture))               //Step2: Change the interference baseline
3:    for s,t = 1: max max(size(obj))          //Step1: Solve one complex coherence coefficient
4:      Δphase = 2π/λ*d(s,t);
5:      I1 = Σ(1 + I0*cos(Δphase)); I2 = Σ(1 + I0*cos(Δphase + π/2));      //Interference fringe expression
6:      I3 = Σ(1 + I0*cos(Δphase + π)); I4 = Σ(1 + I0*cos(Δphase + 3π/2));
7:    end
8:    absIM(i,j) = actan((I4 − I2)/(I1 − I3));              //Amplitude and phase
9:    angIM(i,j) = (I4 − I2)/2/sin(angIM);
10:  end
11:  IM = absIM*exp(angIM);
12:  im = fftshift(ifft2(IM));                  //Step3: Image reconstruction
If the interference fringe coordinate system is omitted, the fiber interference fringe intensity of any interference baseline can be expressed as [15]:
I = I 0 ( 1 + | μ | cos ( φ ) )
In current interferometric measurements, the most commonly used phase calculation method is the four-step phase shift algorithm. This type of algorithm has extremely high measurement accuracy, anti-interference, and high calculation speed [16,17,18]. In this study, the visibility and phase of interference fringes are calculated by a four-step phase shift formula. The result can be expressed as [17,18]:
{ I 1 = I 0 ( 1 + | μ | cos ( φ ) ) I 2 = I 0 ( 1 + | μ | cos ( φ + π / 2 ) ) I 3 = I 0 ( 1 + | μ | cos ( φ + π ) ) I 4 = I 0 ( 1 + | μ | cos ( φ + 3 π / 2 ) )
{ φ = tan 1 ( I 4 I 2 I 1 I 3 ) | μ | = I 4 I 2 2 sin φ
We adopt the parameters in Table 1 to simulate a point source imaging process rectangular microlens array, and the results are shown in Figure 3. For the numerical simulation, we used a square-arranged microlens array and a centrally symmetric arrangement. Figure 3a presents the light intensity distribution of the target source. Figure 3b shows the detected spatial frequency distribution. Figure 3c shows the image obtained by an inverse Fourier transform. The results validate the imaging principle and simulation model.

4. Numerical Simulation of the Imaging System

According to the point-source imaging simulation, as compared with the target source (Figure 3a), the reconstructed image (Figure 3c) is obviously poor. Therefore, we perform further simulations by changing the system parameters in order to evaluate the imaging resolution and quality.
First, we only change Nob, to simulate a resolution panel. With the same sampling interval, Nob determines the maximum interference baseline length (Bmax). According to imaging Equations (1) and (2), the larger the interference baseline length, the richer the high-frequency information contained, and the higher the imaging resolution. The imaging achieved by the system when the maximum baseline (Bmax) is 1.5, 1.7, and 1.9 is shown in Figure 4a–c. It is obvious that the image resolution was up to the length of the maximum baseline. This result also provides a basis for improving the resolution of the imaging system.
The Nyquist sampling density affects the detected spatial frequency distribution. An extremely small sampling density can result in the loss of target information. In contrast, a high sampling density leads to information redundancy. However, the minimum sampling interval is limited to the sub-aperture size (r). Therefore, we have performed several repeated digital simulation experiments on the relationship between sampling density and imaging quality.
In this regard, we introduce a sharpness evaluation function to evaluate the image quality. The commonly used spatial domain functions for sharpness evaluation are listed in the literature, among which the spatial domain functions represented by the variance function and the entropy function perform better [19]. We use two imaging quality evaluation methods based on gray variance and entropy function.
We use the resolution board as the source target and change the sampling interval. The relationship between the value of the gray-scale variance and entropy and the sampling interval is shown in Figure 5. We see that the imaging quality can be improved by increasing the sampling density, although these two evaluation functions are not sensitive to small local changes. In addition, the imaging system has no oversampling when the minimum sampling density is equal to the sub-aperture size (△u = 0.33). Thereafter, we used a photograph of an airplane (Figure 6) as the original target to perform the same simulation. Based on the results (Figure 7), we draw the same conclusion, although we compare the image quality visually in this case.
It is well known that different sources have different frequency information. For some sources, the frequency information was concentrated in the low-frequency part; therefore, the sampling density was increased to obtain a high-quality image. This resulted in a corresponding increase in system complexity. A satellite was used as the detection target; from the simulation results (Figure 8), it can be seen that the spectrum information of the satellite image is concentrated in the low-frequency part. We attempted to improve the image quality without increasing the number of microlenses, i.e., the complexity of the system. Then, the complex rectangular sub-aperture arrangement method was adopted to adjust the sampling density at high and low frequencies separately.
In the results presented in Figure 8, the sampling interval is uniformly set to △u = 0.33. In Figure 9, the sampling spectrum is divided into two parts, a high frequency and a low frequency, and the sampling intervals are respectively set to △u1 = 0.19, △u2 = 0.38. Similarly, the sampling intervals are set to △u1 = 0.33, △u2 = 0.66, △u3 = 1.3 in Figure 10. The above settings ensure that the number of microlenses used in the three sampling modes is nearly equal.
The simulation results are shown in Figure 8, Figure 9 and Figure 10. Visual comparison shows no difference in image quality. This shows that the complex rectangular sub-aperture arrangement method proposed in this paper is feasible. We use variance and entropy to evaluate the subtle differences between the images. The results are presented in Table 2, which proves that higher-quality images could be obtained with the same number of microlenses.

5. Conclusions

This paper proposes a cutting-edge electro-optical imaging technology based on a microlens array and fiber interferometer, a centrally symmetric baseline pairing method, and a complex rectangular sub-aperture arrangement method. We then briefly analyze the theoretical basis of the imaging system and the method of analyzing the fringes. Finally, the imaging system was verified by digital simulation.
The main conclusions are as follows: (1) the resolution of the imaging system can only be improved by increasing the maximum baseline length; (2) increasing the Nyquist sampling density increases the density of the spatial frequency points and improves the imaging quality; (3) by adopting the complex rectangular sub-aperture arrangement, better image quality can be obtained without increasing the complexity of the imaging system.

Author Contributions

F.Y. conceived and designed the experiments; F.L. and F.Y. performed the experiments; F.P. have done a lot of literature research work on the research background; P.T. and W.Y. wrote the paper.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baldwin, J.E.; Beckett, M.G.; Boysen, R.C.; Burns, D.; Buscher, D.F.; Cox, G.C.; Haniff, C.A.; Mackay, C.D.; Nightingale, N.S.; Rogers, J.; et al. The first images from an optical aperture synthesis array: Mapping of Capella with COAST at two epochs. Astron. Astrophys. 1996, 306, L13–L16. [Google Scholar]
  2. Benson, J.A.; Mozurkewich, D.; Jefferies, S.M. Active optical fringe tracking at the NPOI. In Proceedings of the SPIE 3350, Kona, HI, USA, 20–28 March 1998. [Google Scholar] [CrossRef]
  3. Benson, J.A.; Hutter, D.J.; Johnston, K.J.; Zavala, R.T.; White, N.M.; Pauls, T.A.; Gilbreath, G.C.; Armstrong, J.T.; Hindsley, R.B. NPOI: Recent technology and science. In Proceedings of the International Society for Optics and Photonics, Glasgow, UK, 21–25 June 2004; pp. 464–471. [Google Scholar]
  4. Young, J.; Haniff, C.; Buscher, D. Interferometric imaging of geo-synchronous satellites with ground-based telescopes. In Proceedings of the 2013 IEEE Aerospace Conference, Big Sky, MT, USA, 2–9 March 2013; pp. 1–9. [Google Scholar] [CrossRef]
  5. Duncan, A.L.; Kendrick, R.L. Segmented Planar Imaging Detector for Electro-Optic Reconnaissance. U.S. Patent No. 8,913,859, 16 December 2014. [Google Scholar]
  6. Kendrick, R.; Thurman, S.T.; Duncan, A.; Wilm, J.; Ogden, C. Segmented planar imaging detector for eo reconnaissance. In Computational Optical Sensing and Imaging; Optical Society of America: Washington, DC, USA, 2013. [Google Scholar] [CrossRef]
  7. Thurman, S.T.; Kendrick, R.L.; Duncan, A.; Wuchenich, D.; Ogden, C. System Design for a SPIDER Imager. In Frontiers in Optics; Optical Society of America: Washington, DC, USA, 2015. [Google Scholar] [CrossRef]
  8. Bashkansky, M.; Lucke, R.L.; Funk, E.; Rickard, L.J.; Reintjes, J. Two-dimensional synthetic aperture imaging in the optical domain. Opt. Lett. 2002, 27, 1983–1985. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, S.; Rao, C.; Jiang, W. Diffraction-limited performance and aberration tolerance of a sparse optical synthetic aperture imaging system with four sub-apertures. J. Mod. Opt. 2008, 55, 1771–1785. [Google Scholar] [CrossRef]
  10. Lee, D.J.; Weiner, A.M. Optics Express. Optical phase imaging using a synthetic aperture phase retrieval technique. Opt. Express 2014, 22, 9380–9394. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, P.; Brambilla, G.; Ding, M.; Semenova, Y.; Wu, Q.; Farrell, G. High-sensitivity, evanescent field refractometric sensor based on a tapered, multimode fiber interference. Opt. Lett. 2011, 36, 2233–2235. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Q.; Farrell, G.J. All-fiber multimode-interference-based refractometer sensor: Proposal and design. Opt. Lett. 2006, 31, 317–319. [Google Scholar] [CrossRef] [PubMed]
  13. Edwards, C.A.; Presby, H.M.; Dragone, C. Ideal microlenses for laser to fiber coupling. J. Lightw. Technol. 1993, 11, 252–257. [Google Scholar] [CrossRef]
  14. Cristóbal, G.; Javidi, B.; Vallmitjana, S. Information Optics; Springer: Berlin, Germany, 2009. [Google Scholar]
  15. Lauterborn, W.; Kurz, R.N.T.; Wiesenfeldt, D.P.M. Fourier Optics; Springer: Berlin, Germany, 1995. [Google Scholar]
  16. Yi, J.; Huang, S. Modified fourier transform profilometry for the measurement of 3-D steep shapes. Opt. Lasers Eng. 1997, 27, 493–505. [Google Scholar] [CrossRef]
  17. Zhang, H.; Lalor, M.J.; Burton, D.R. Robust, accurate seven-sample phase-shifting algorithm insensitive to nonlinear phase-shift error and second-harmonic distortion: A comparative study. Opt. Eng. 1999, 38, 1524–1533. [Google Scholar] [CrossRef]
  18. Schwider, J.; Burow, R.; Elssner, K.E.; Grzanna, J.; Spolaczyk, R.; Merkel, K. Applied Optics. Digital wave-front measuring interferometry: Some systematic error sources. Appl. Opt. 1983, 22, 3421. [Google Scholar] [CrossRef] [PubMed]
  19. Eskicioglu, A.M.; Fisher, P.S.; Chen, S.Y. Image quality measures and their performance. IEEE Trans. Commun. 1995, 43, 2959–2965. [Google Scholar] [CrossRef]
Figure 1. (a) Rectangular microlens array using a symmetric pairing method; (b) microlens array using a complex rectangular sub-aperture arrangement method; (c) structure of the fiber interferometer.
Figure 1. (a) Rectangular microlens array using a symmetric pairing method; (b) microlens array using a complex rectangular sub-aperture arrangement method; (c) structure of the fiber interferometer.
Applsci 09 01331 g001
Figure 2. Schematic of the cutting-edge electro-optical imaging technology.
Figure 2. Schematic of the cutting-edge electro-optical imaging technology.
Applsci 09 01331 g002
Figure 3. (a) Light intensity distribution of target source; (b) detected spatial frequency distribution; (c) image obtained by inverse Fourier transform.
Figure 3. (a) Light intensity distribution of target source; (b) detected spatial frequency distribution; (c) image obtained by inverse Fourier transform.
Applsci 09 01331 g003
Figure 4. (a) Image obtained when Bmax = 1.5 m; (b) Image obtained when Bmax = 1.7 m; (c) Image obtained when Bmax = 1.9 m.
Figure 4. (a) Image obtained when Bmax = 1.5 m; (b) Image obtained when Bmax = 1.7 m; (c) Image obtained when Bmax = 1.9 m.
Applsci 09 01331 g004
Figure 5. Relationship between the variance and entropy of the imaging result and the sampling interval.
Figure 5. Relationship between the variance and entropy of the imaging result and the sampling interval.
Applsci 09 01331 g005
Figure 6. (a) Light intensity distribution of the source; (b) airplane image spectrum center slice.
Figure 6. (a) Light intensity distribution of the source; (b) airplane image spectrum center slice.
Applsci 09 01331 g006
Figure 7. Image reconstruction result of digital simulation experiment when (a) △u = 0.33, Nap = 200; (b) △u = 0.40, Nap = 166; (c) △u = 0.50, Nap = 134.
Figure 7. Image reconstruction result of digital simulation experiment when (a) △u = 0.33, Nap = 200; (b) △u = 0.40, Nap = 166; (c) △u = 0.50, Nap = 134.
Applsci 09 01331 g007
Figure 8. (a) Light intensity distribution of the source; (b) mutual intensity spectrum detected with the rectangular arrangement (△u = 0.33, N = 200); (c) actual imaging of the system.
Figure 8. (a) Light intensity distribution of the source; (b) mutual intensity spectrum detected with the rectangular arrangement (△u = 0.33, N = 200); (c) actual imaging of the system.
Applsci 09 01331 g008
Figure 9. (a) Mutual intensity spectrum detected with the complex rectangular sub-aperture arrangement method (△u1 = 0.66, △u2 = 0.99, N = 200); (b) partial enlargement of the spectrum; (c) actual imaging of the system.
Figure 9. (a) Mutual intensity spectrum detected with the complex rectangular sub-aperture arrangement method (△u1 = 0.66, △u2 = 0.99, N = 200); (b) partial enlargement of the spectrum; (c) actual imaging of the system.
Applsci 09 01331 g009
Figure 10. (a) Mutual intensity spectrum detected with the complex rectangular sub-aperture arrangement method (△u1 = 0.33, △u2 = 0.66, △u3 = 1.3, N = 200); (b) partial enlargement of the spectrum; (c) actual imaging of the system.
Figure 10. (a) Mutual intensity spectrum detected with the complex rectangular sub-aperture arrangement method (△u1 = 0.33, △u2 = 0.66, △u3 = 1.3, N = 200); (b) partial enlargement of the spectrum; (c) actual imaging of the system.
Applsci 09 01331 g010
Table 1. Simulation parameters of the imaging system.
Table 1. Simulation parameters of the imaging system.
ParametersSymbolValue
Object planeScene side lengthD30 m
Number of samplesNob200 × 200
Aperture planeSampling intervalΔuΔv0.33 m−1
Number of microlensesNap200 × 200
Diameter of microlensr1 mm
Image planeNumber of samplesNim1200 × 1200
SystemWavelengthλ600 nm
Object distancez105 m
Table 2. Imaging process and reconstruction algorithm.
Table 2. Imaging process and reconstruction algorithm.
The Number of Complex ArrangementsParameterVarianceEntropy
1△u = 0.33, N = 2001.5087 × 1083.7735
2△u1 = 0.66, △u2 = 0.99, N = 2001.4983 × 1083.7765
3△u1 = 0.33, △u2 = 0.66, u3 = 1.3, N = 2001.6364 × 1083.8652

Share and Cite

MDPI and ACS Style

Yang, F.; Yan, W.; Tian, P.; Li, F.; Peng, F. Electro-Optical Imaging Technology Based on Microlens Array and Fiber Interferometer. Appl. Sci. 2019, 9, 1331. https://doi.org/10.3390/app9071331

AMA Style

Yang F, Yan W, Tian P, Li F, Peng F. Electro-Optical Imaging Technology Based on Microlens Array and Fiber Interferometer. Applied Sciences. 2019; 9(7):1331. https://doi.org/10.3390/app9071331

Chicago/Turabian Style

Yang, Fan, Wei Yan, Peng Tian, Fanxin Li, and Fupin Peng. 2019. "Electro-Optical Imaging Technology Based on Microlens Array and Fiber Interferometer" Applied Sciences 9, no. 7: 1331. https://doi.org/10.3390/app9071331

APA Style

Yang, F., Yan, W., Tian, P., Li, F., & Peng, F. (2019). Electro-Optical Imaging Technology Based on Microlens Array and Fiber Interferometer. Applied Sciences, 9(7), 1331. https://doi.org/10.3390/app9071331

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop