The depolarization optical engine is mainly composed of LCOS, depolarization spectroscopic prism, and high-uniform illumination optical system, which can increase the resolution of the image plane and reduce the stray light in the optical engine.
2.2.2. Stray Light Analysis and Suppression of Dim Target Simulation System
- (1)
Causes and effects of stray light
The stray light in the simulation system of dim targets mainly comes from two aspects. First, the dark state light leakage phenomenon is caused by the difference between the response ability of LCOS to incident light of different wavelengths. Second, the illumination optical system does not match the spectroscopic prism. The spectroscopic prism is sensitive to the incident light angle [
12], and its outgoing light cannot maintain a single-line polarization state, so more imaging polarized light is converted into polarized stray light. After the polarized stray light is emitted by PBS, the image surface of the detector will be covered with a layer of uneven stray light, which will reduce the contrast of the image surface, cause the centroid of the star point to deviate from the geometric center, and reduce the position accuracy of the star point.
Figure 5 shows the grayscale extraction of the star point under different stray light distributions.
Figure 5 shows the star point grayscale extraction diagram under different stray light distribution. If there is no stray light affecting the simulated star point, the theoretical position (x, y) of the simulated star point should be (5.5, 5.5) and the theoretical star point grayscale on the camera image plane should be 189.0. When there is a large amount of stray light in the simulation system, the overall background grayscale of the camera image plane will be increased. If the stray light is more uniform, the signal-to-noise ratio of the simulated star point will be reduced, resulting in difficulty in extracting the centroid position of the star point. If the stray light is not uniform, some maximum values of the stray light will easily affect the centroid position of the nearby star point, the grayscale of the star point will increase to 199.0, and the star point position will be mistakenly received by the camera as (6, 4).
- (2)
Stray light suppression method
- (a)
Dark state light leakage suppression
When the LCOS is in the off state, the incoming polarized light is not regulated by the liquid crystal molecules, so it cannot exit the simulation system to achieve the off state. However, when the liquid crystal molecules appear under the dark state light leakage phenomenon, the light on the surface of the LCOS with the emission closed will have different polarization states, resulting in a partially polarized stray light emission simulation system. The quarter wave plate can be placed on the surface of LCOS to phase compensate the LCOS outgoing light, reduce the stray light of the outgoing simulation system, and achieve the effect of suppressing the dark state light leakage. For the quarter wave plate compensation method [
13], the dark state light leakage
R1 of the LCOS outgoing light is shown in Equation (2):
where
α is the angle between the transmission axis of the polarizer and the orientation of the liquid crystal molecule at the incident point,
is the Jones vector of the incoming polarized light,
M is the Jones matrix of the LCOS liquid crystal molecule, and
W is the Jones matrix of the wave plate [
14], as shown in Equations (3) and (4):
where
is the distortion angle of the liquid crystal molecule,
is the angle between the optical axis of the wave plate and the direction of the liquid crystal molecule,
is the central wavelength of the incident light,
is the difference of the refractive index of the material to the light o and e,
d is the thickness of the liquid crystal molecule, and
X is the intermediate variable, as shown in Equation (5).
After one-quarter slide compensation, the outgoing light wavelength and dark state light leakage of the system are shown in
Figure 6, and the maximum value of LCOS dark state light leakage is 0.0035%.
- (b)
Calculation of critical incident angle matching of spectroscopic film
The incident beam angle of the illumination optical system affects the effective thickness of the spectroscopic film in PBS [
15]. If the incident angle of the illumination beam does not match the critical incident angle of the spectroscopic film, S light and P light in the transmitted and reflected light will be doped with each other, which greatly reduces the stray light suppression ability of the optical engine. Therefore, it is necessary to determine the critical incidence angle of the spectroscopic film to match the beam angle of the illumination system.
Starting from the application scenario of the target simulator [
16], the stray light is equivalent to the detector noise. When the detector receives the signal, the number of electrons
produced is represented by Equation (6):
where
is the signal equivalent magnitude number, the optical system entry pupil diameter is
, the optical system transmittance is
, the detector exposure time is
, the signal wavelength is
, the spectral efficiency of the detector is
,
is Planck’s constant, and
is the propagation speed of light in vacuum.
The system noise equivalent number of electrons
is represented by Equation (7):
where
is the incident illuminance,
is the point source transmittance of the optical engine, both of which are functions of the incident azimuth angle
and incident space angle
, and
is the dark current noise detector.
Therefore, the signal-to-noise ratio (SNR) is expressed by Equation (8):
Taking +11 magnitude star as an example (+11 magnitude star illuminance is about 1.05 × 10
−10 lx), when
SNR = 5, the star sensor can identify the target [
17], and the corresponding optical engine PST of the target simulator should be
. For an optical engine, its
can be represented by the transmittance
of P light, the transmittance
of S light, the reflectance
of P light, and the reflectance
of S light, as shown in Equation (9):
Based on the analysis of the working principle of the illumination system and the optical engine, combined with the propagation path of the light in the optical engine, it can be seen that the beam forming two LCOS image planes can be regarded as the spectroscopic film passing through PBS twice, so the process is simplified by using the characteristic matrix of the double-layer film system [
18], as shown in Equation (10):
where
is the characteristic matrix of a single-layer spectroscopic film determined by film system and substrate parameters, P component film admittance is
, S component film admittance is
,
is the refractive index of the substrate film,
is the beam incident azimuth, film phase thickness is
, and
is the film layer spacing.
Then, the reflectance R and transmittance T of the film system can be shown as Equation (11):
where
is the admittance of the incident medium.
Considering that the polarized stray light in the optical engine propagates in three-dimensional space, and for the spectroscopic film system, the angle affecting its reflectivity and transmittance is the angle between the incident light and the normal of the film layer; the light azimuth angle
of the second incident spectroscopic film can be represented by the azimuth angle
of the incident optical engine and the space angle
, as shown in Equation (12):
Assuming that the light source emits a uniform beam that meets the critical incidence angle of the depolarization spectroscopic prism, the corresponding PST size can be obtained when the incident azimuth varies from 0° to 20° at different incident space angles, as shown in
Figure 7.
When PST = 6.88 × 10−9, the critical incidence angle of the depolarizing spectroscopic prism is about 9.2°. When the incident angle is greater than the critical incident angle, the PST of the optical engine will continue to increase, which will cause the stray light suppression ability of the optical engine to decrease, and PST can not meet the requirement of signal recognition. Therefore, the divergence angle of the illumination optical system should match the critical incidence angle of the depolarization spectroscopic prism, so that the depolarization optical engine can effectively suppress stray light.
2.2.3. Analysis of Illumination Optical System Design Objectives
According to the optimization results of the depolarizing spectroscopic prism, the illumination optical system of the dim target simulator needs to meet the constraint condition that the prism incidence angle is less than ±9.2°. At the same time, the uniformity of irradiation surface should be taken into account. When it is better than 92%, the uniformity has little influence on the accuracy of magnitude simulation. The LCOS side length is 16.24 mm, so the diameter of the illumination area should be greater than the bevel side length after splicing 22.96 mm. Considering the allowance for installation, the illumination area is
. The design objectives of the illumination optical system are shown in
Table 3.
- (1)
Composition and structure of illumination optical system
The illumination optical system mainly consists of a light source, CPC, homogenizing rod, and illumination lens, and its system structure is shown in
Figure 8. The light beam emitted by the light source is distributed through the CPC twice, the angle of the output beam is reduced and then homogenized by the light rod, and finally, the illumination lens is used to achieve small-angle and uniform illumination.
- (2)
Selection of light source and design of each component unit
- (a)
Light source selection
We selected LED as the target simulator light source and determined the LED luminous flux. The energy transfer function E of an optical system is defined as the ratio of the outgoing light energy to the incident light energy, and the outgoing luminous flux φ
′ can be expressed by Equation (13):
where
is the incident luminous flux,
Ecpc is the CPC energy transfer function,
E1 is the homogenizing rod energy transfer function,
E2 is the energy transfer function of the collimation system,
E3 is the energy transfer function of the depolarization spectroscopic prism,
E4 is the LCOS energy transfer function, and
E5 is the energy transfer function of the collimation optical system. According to the experimental and empirical data,
KCPC = 0.115,
E1 = 0.69,
E2 = 0.89,
E3 = 0.25,
E4 = 0.95, and
E5 = 0.851. In order to meet the target illuminance of 10
−10 lx and diameter of 25 mm, the minimum luminous flux
of the incident light is shown in Equation (14).
The result is .
- (b)
CPC design
The CPC is designed to collect the large-angle outgoing light of the LED, and after angle modulation, the outgoing light is superimposed to realize the control of the luminous intensity distribution in the irradiation center and the edge area and improve the uniformity and increase the utilization rate of light energy [
19]. The working principle is shown in
Figure 9.
The composite parabolic mirror is mainly formed by the interceptor segments
and
of parabola
and parabola
, rotating around the symmetry axis, which can be obtained by analyzing the edge ray principle in the non-imaging optics theory. Parabola
is shown in Equation (15):
where
,
,
is the focal length,
,
is the radius of the LED,
,
is the radius of the light outlet hole
, and
is the maximum beam absorption angle of the composite parabolic mirror.
The luminous intensity
at the exit of the composite paraboloid mirror includes two parts: the luminous intensity
directly emitted by the LED and the luminous intensity
reflected by the composite paraboloid mirror. If the luminous intensity in the normal direction of the LED is
, the luminous intensity at the angle
is
and the luminous intensity at the exit of the composite paraboloid mirror can be shown by Equation (16):
where
is the function of
related to
and
,
is the outgoing light angle of the LED corresponding to
,
,
is the ordinate of the intersection point of the light with the incident angle of
emitted by the LED and the CPC,
is the reflectance of the reflector surface of the composite parabolic mirror,
is the length of the composite parabolic mirror, and
.
When is larger, the range of is larger, and the reflected light is reflected and compensated on both sides of the irradiation surface. The radius of the LED is a = 1.88 mm and the incident light angle is 60°. In order to ensure that the CPC length is appropriate, and are used to calculate the CPC focal length f = 2.042 mm. CPC length L = 22.844 mm; CPC light output diameter 2b = 10.81 mm.
- (c)
The design of the homogenizing rod
We designed a homogenizing rod for high uniform lighting. The homogenizing rod uses the principle of total reflection to reflect the light incident inside the rod many times and then emit it. The aperture ratio of the homogenizing rod determines the total reflection times, and the incident irradiation is divided according to the incident angle, so as to improve the irradiation uniformity of the system [
20]. The principle of the homogenizing rod is shown in
Figure 10, in which ① to ④ respectively represent the equivalent position of the edge light of
to
light from different incident angles.
The homogenizing rod shall satisfy Equation (17):
where
m is the maximum number of total reflections in the homogenizing rod,
n is the refractive index of the homogenizing rod,
is the length of the homogenizing rod,
is the cross-section width of the homogenizing rod,
is the division angle of the homogenizing rod, and
is the incidence angle of the homogenizing rod.
Since the CPC incident angle in front of the homogenizing rod is 17° and the outlet diameter is 10.81 mm, . Considering that the effective irradiation surface size is 25 mm, the maximum number of total reflections should be greater than or equal to 25/10.81 = 2.3. The maximum number of total reflections is 3, and since JGS3 heat-resistant quartz glass n = 1.487 is selected as the homogenizing rod, can be calculated.
- (d)
Illumination lens design
The system spot diagrams and axial aberration curve are shown in
Figure 11. According to the spot diagrams, the color difference correction of each field of view is good, the dispersion spot of the optimized image surface is not large, and the energy concentration is high. According to the axial aberration map, the color difference is corrected at the maximum aperture, although there is a certain secondary spectrum; but because this system is an illumination system, there is no need to correct the secondary spectrum, and the system aberration meets the requirements of use.
- (3)
Analysis of simulation results of illumination optical system
Tracepro 7.3 software was used to simulate the illumination system. The light tracing diagram and the outgoing light irradiance diagram of the illumination system are shown in
Figure 12 and
Figure 13.
The uniformity of luminous intensity
over a range of
can be expressed by Equation (18):
where
is the maximum luminous intensity and
is the minimum luminous intensity. The uniformity was calculated to be 95%.
Therefore, the design parameters of the illumination optical system are better than ±9.2°, the uniformity is better than 92%, and the irradiation surface is , which meets the requirements of use.