Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Freeform reflector light source used for space traceable spectral radiance calibration on the solar reflected band

Open Access Open Access

Abstract

For purpose of improving the accuracy of in-orbit radiometric calibration of Chinese Space-based Radiometric Benchmark (CSRB) reference payload on the reflected solar band and reduce resource consumption, this paper proposed a freeform reflector radiometric calibration light source system based on the beam shaping characteristics of the freeform surface. The design method of initial structure discretization based on Chebyshev points was used to design and solve the freeform surface, and the feasibility of the design method was verified by optical simulation. The designed freeform surface is machined and tested, and the test results show that the surface roughness root mean square (RMS) of the freeform reflector is 0.61 μm, which indicates that the continuity of the machined surface is good. The optical characteristics of the calibration light source system are measured, the results show that the irradiance uniformity and radiance uniformity are better than 98% in the effective illumination area of 100 mm × 100 mm on the target plane. The constructed freeform reflector calibration light source system can meet the requirements of large area, high uniformity and light weight for onboard calibration of the payload of the radiometric benchmark, improving the measurement accuracy of spectral radiance on the reflected solar band.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Radiation calibration is an important step in the development of space optical remote sensing load and an important prerequisite for obtaining high-quality remote sensing data [1]. The calibration light source is an important part of the onboard radiometric calibration system, which determines the consumption of space, load, energy, and other resources by the whole calibration system to a large extent. At the same time, the characteristics of the target light source also directly limit the accuracy of radiation calibration. At present, there are still some problems in the radiometric calibration of space remote sensors, such as the inconsistency of radiometric calibration values between different remote sensors, the incomparability of data and the low accuracy of the radiometric calibration [2]. To this end, China put forward the concept of a space radiation measurement reference satellite in 2006 and launched the Chinese Space-based Radiometric Benchmark (CSRB) project [3,4]. The purpose of the project is to launch a radiation reference satellite with self-calibration capability, use a new in-orbit calibration system to cross-calibrate the in-orbit satellite, establish a high-precision space radiation reference, and unify the onboard radiation scale, which coincides with the Traceable Radiometry Underpinning Terrestrial-and Helio-Studies (TRUTHS) and the Climate Absolute Radiance and Refractivity Observatory (CLARREO) projects proposed by the National Physical Laboratory (NPL) and National Aeronautics and Space Administration (NASA) respectively [5,6].

Under the support of CSRB, Changchun Institute of Optics, Fine Mechanics and Physics (CIOMP), Chinese Academy of Sciences proposed an onboard radiometric calibration scheme [7], which is traceable to the Space Cryogenic Absolute Radiometer (SCAR) [8], the purpose is to establish an onboard radiance benchmark with an uncertainty of better than 1% on the reflected solar band. The calibration system transfers the reference of SCAR to a Transmitted Radiometer (TR) that can measure both radiant power and radiance [9]. And then the TR and the Earth-Moon Imaging Spectrometer (EMIS) simultaneously observe the calibration light source to complete the full-field and full-aperture radiation calibration. The primary role of the TR is to convert from the spectral power responsivity traceable to the SCAR to the spectral radiance responsivity and transfer it to the EMIS (the detailed traceability chain will be described in section 2.1). In the calibration process, the target light source needs to fill the field of view and aperture of the EMIS and TR, and therefore its size needs to be large enough. At the same time, due to the mismatch between the observation area of the EMIS and TR, to achieve higher calibration accuracy, a very strict requirement is put forward for the uniformity of the whole spectral calibration light source. The existing onboard calibration light source system on reflected solar band usually takes the sunlight or halogen tungsten lamp as the light source and uses the integrating sphere or diffuser to homogenize the light to form a uniform target light source [1013]. When the diffuser is used to homogenize the light, the intensity and stability of the light will change with the change of the incident light field angle, which is easy to introduce the response error of the detector caused by the uneven incident angle and position of the incident light. Using the integrating sphere to homogenize light can achieve high uniformity, but to cover the aperture of EMIS, the volume and weight of the integrating sphere will be very large, which will be a huge consumption of onboard resources. In addition, the off-axis parabolic mirror can obtain a large area extended light source, but the optical path deflection and collimation will require multiple mirrors, which will undoubtedly increase the complexity of the calibration system structure and cannot meet the uniformity requirements. Therefore, due to the current situation that the onboard lighting system cannot take into account both large area, high uniformity and simple structure, we urgently need a new onboard calibration light source system to meet the requirement of radiation calibration accuracy of the reference satellite.

With the development of freeform surfaces, the lighting output enters the era of “custom”. The system precisely controls the light direction to reach the predetermined area, and the design freedom is high [1416]. This not only makes the spatial layout of the optical system more flexible but also improves the energy utilization efficiency and simplifies the lighting system structure. In addition, the size and uniformity of the outgoing spot of the freeform can be accurately controlled by a reasonable surface design. The beam shaping ability of freeform surface brings a new idea for the design of the onboard whole spectral calibration light source system, which enables us to obtain large-size and high-uniformity target light source by using a single freeform surface.

In this paper, a direct discretization method of a freeform surface is proposed based on existing design methods, and the design, simulation and analysis of the freeform reflector of the SCAR-based calibration system are presented. Combining the homogenization ability of the integrating sphere and diffuser and the beam shaping ability of the freeform reflector, a new onboard calibration light source system is constructed. The system not only takes the energy, load and space utilization onboard into account but also provides a large size calibration light source with ultra-high uniformity. The rest of this paper is organized as follows. In Section 2, we introduce the design concept of the calibration light source system and the freeform reflector and present the design requirements of the light source system in combination with the SCAR-based calibration scheme. In Section 3, we propose the design and solving method of the freeform surface. In Section 4, we design, simulate, manufactured and test the freeform calibration light source system in detail. In Sections 5 and 6, we give the discussion and conclusion respectively.

2. Design concept

2.1 Freeform calibration light source system design concept

The light path of the whole spectral calibration light source (composed of halogen tungsten lamp integrating sphere, freeform reflector and diffuser) of the SCAR-based radiation calibration system is shown by the bold arrow in Fig. 1. The freeform calibration light source system adopts an off-axis reflection lighting scheme, which avoids the influence of color difference introduced by the transmission scheme on calibration. We effectively combined the homogenizing light ability of the integrating sphere and the diffuser with the beam shaping ability of the freeform surface, a small diameter integrating sphere with a halogen tungsten lamp is used to vertically illuminate the diffuser by the freeform reflector as the whole spectral calibration light source. The single-mirror lighting method can shorten the optical path effectively and simplify the optical path structure to the greatest extent. The small-diameter integrating sphere can output stronger spectral radiance with a low-power halogen tungsten lamp, which greatly improves the energy utilization rate of the light source system. The off-axis reflection lighting path avoids all the components in the calibration system, making the entire on-board layout very compact and reducing the reflection of light in the system. Therefore, this scheme not only reduces the impact of stray light on the calibration accuracy, but also improves the space utilization of the SCAR-based radiation calibration system.

 figure: Fig. 1.

Fig. 1. Light path of the whole spectral calibration light source system of SCAR-based radiation calibration.

Download Full Size | PDF

As shown in Fig. 1, a multi-channel monochromatic light source was used to achieve the multi-spectral power calibration of the Transfer Radiometer (TR) traced to SCAR. The multi-channel monochromatic light source is connected via fibers to a collimating swing arm, which can introduce the monochromatic light into TR and SCAR in turn (the fibers are not shown in Fig. 1). The TR is used to achieve high-precision conversion and measurement from power to radiance, and establish secondary reference. The uncertainty of radiance measurement is better than 0.55% [9]. The multi-spectral radiance of the whole spectral light source can be calibrated by the TR. Since the spectral curve of the halogen tungsten lamp in the whole spectral light source is smooth without mutation, the secondary reference can be converted into the hyperspectral radiance reference by hyperspectral curve reconstruction. As shown by the thin arrow in Fig. 1, TR and hyperspectral imaging system, the Earth-Moon Imaging Spectrometer (EMIS), simultaneously observe the diffuser to achieve high precision radiometric calibration on the reflected solar band. The placement of the diffuser ensures that EMIS and TR can observe the diffuser in a symmetric direction of ±45°, which minimizes the difference between the reflectance of the diffuser at different angles. However, because the observation area of EMIS and TR on the diffuser are different, the irradiance of the reflected light from the freeform reflector to the diffuser is supposed to very uniform to ensure higher radiation value transfer accuracy. The size of the diffuser is 120 mm × 120 mm. EMIS and TR have a maximum observation area of 90 mm on the diffuser, so 100 mm × 100 mm is sufficient to cover their aperture and field of view. After analysis, the radiance uniformity is supposed to be better than 98% in the effective illumination area of 100 mm × 100 mm. The radiance uniformity of the diffuser is affected by the uniformity of its reflectance and the uniformity of the irradiance of the incident light. The diffuser we used is of very good reflectance uniformity, so we want the irradiance uniformity of the incident light of the diffuser, that is, the reflected light of the freeform reflector to be also better than 98%.

2.2 Freeform reflector design concept

Many ideas have been put forward for the construction and solution of freeform surface. However, the current freeform surface design methods are difficult to form a smooth and continuous freeform reflector that can be used for onboard calibration light source system and generate large area uniform light spots. Schruben [17] constructed an illumination system based on mathematical principles and obtained the mathematical model of the elliptic MA (Monge-Ampère) equation of the asymmetric freeform reflector illumination system. But the partial differential equation has high nonlinearity, making its solution very complicated [18,19]. Oliker et al. [20,21] proposed the supported cone method, which is based on illumination discretization. Simultaneous Multiple Surface (SMS) is used to construct the freeform surface by concatenating multiple surfaces, but this makes the resulting surface shape not smooth and continuous. Wu et al. [2224] solved the elliptic MA equation of freeform reflector by using the Newton descent method and obtained a continuous smooth surface model. However, this method only applies to the design of transmissive freeform surfaces. Luo [16,25] et al. first constructed a point light source model by separating variables, then modified the preset distribution according to the actual irradiance distribution, and finally established a new irradiance distribution by iterative function, thus completing the design of freeform mirror. A rectangular spot distribution with an irradiance uniformity of 81.08% was achieved by this method. Based on the existing freeform surface design ideas, we propose a new discretization method for the initial structure of the freeform surface, based on the definition of Chebyshev points, a two-dimensional extension is carried out to make it a better match between the projection grid of the light emitted from the light source, to improve the continuity of the surface and the uniformity of the reflected light. In the next section, we introduce the method of establishing and solving the freeform surface mathematical model.

3. Freeform reflector design principle

3.1 Establishment of the mathematical model

The essence of the design and solving method proposed in this section is to transfer the energy from the light source to the target plane in a reasonable way by a freeform reflector. As shown in Fig. 2 is the diagram of energy transfer for a freeform reflector. According to the law of conservation of energy, the energy received by the target surface is equal to the energy of the light source incident on the freeform reflector multiplied by the reflectance of the reflector. Then the Jacobi matrix can be used to achieve the energy conversion between each face. As long as the irradiance in each illumination area unit is equal, an irradiance uniformity spot can be obtained on the target plane.

 figure: Fig. 2.

Fig. 2. Diagram of energy transfer for a freeform reflector.

Download Full Size | PDF

As shown in Fig. 3(a) is the 3D diagram of the illumination coordinate system, Fig. 3(b) is the corresponding 2D diagram. According to the design and solving method, first, a spherical coordinate system with a radius of 1 is constructed. All subsequent variables are normalized to facilitate calculation. Assuming that the surface type region of the freeform reflector is ${\Omega _0}$, and the projection region of the luminous region emitted by the light source on the unit sphere is $\mathrm{\rho} (\rm{u},\rm{v})$, where $\mathrm{\rho}$ is the distance between the light source and the freeform reflector, and the region of reflection projection on the target surface is ${\Omega _1}$. The ${\Omega _0}$, ${\Omega _1}$ and point O have been marked in Fig. 1. The topological plane ${\textrm{u}_1} - {\rm{v}_1}$ is titled $\mathrm{\beta}$ from plane $\rm{u} - \rm{v}$. The target surface ${\textrm{t}_{\rm x}} - {{\rm t}_{\rm y}}$ is shifted upward 2 units form plane $\rm{u} - \rm{v}$; ${\textrm{t}_{\rm {x}}} - {\rm{t}_{\rm{y}}}$ is called the target plane of the incident area of the light source. The purpose is to convert the irregular region of incident light on the unit sphere into a plane to facilitate calculation.

 figure: Fig. 3.

Fig. 3. Diagram of the illumination coordinate system: (a) 3D (b) 2D.

Download Full Size | PDF

Assuming that the intensity of the input light source $\textrm{OP}$ is $\textrm{I}$, M is the projection of point P on the plane z = 1 along the line through the south pole (0, 0, -1) of the unit sphere and point P. With coordinates u and v, M can be expressed as $\textrm{I(u,v)}$. We assume that there is no energy loss in the optical system during the design stage. According to the law of conservation of energy, we have:

$$\int\!\!\!\int\limits_{{\Omega _1}} {\textrm{E(}{\rm{t}_{\rm x}},{\rm{t}_{\rm y}})\rm{d}{\rm{t}_{\rm x}}} \rm{d}{\rm{t}_{\rm y}} = \int\!\!\!\int\limits_{{\Omega _{}}} {I(\textrm{u},\textrm{v})d} \Omega .$$
where E is the irradiance of a point at the coordinate $({\rm{t}_{\rm x}},{t_\textrm{y}})$, ${\Omega _1}$ is the illumination area on the target plane which can also be regarded as the diffuser as shown in Fig. 1, and $\Omega $ is the solid angle of the output light. Assuming that plane ${\textrm{u}_1} - {\rm{v}_1}$ is titled along the x-axis, and the obtained point in plane ${\textrm{u}_1} - {\rm{v}_1}$ is S$({u_1},{\rm{v}_1},0)$, as the plane ${\textrm{u}_1} - {\rm{v}_1}$ is titled $\mathrm{\beta}$ from the plane $\rm{u} - \rm{v}$, then point S in the coordinate $\rm{u} - \rm{v}$can be expressed as:
$$\left[ {\begin{array}{c} \rm{u}\\ \rm{v}\\ \rm{w} \end{array}} \right] = \left[ {\begin{array}{ccc} 1&0&0\\ 0&{\cos \mathrm{\beta} }&{\sin \mathrm{\beta} }\\ 0&{ - \sin \mathrm{\beta} }&{\cos \mathrm{\beta} } \end{array}} \right]\left[ {\begin{array}{c} {{\rm{u}_1}}\\ {{\rm{v}_1}}\\ 0 \end{array}} \right].$$

By solving Eq. (2), we can obtain point S in the coordinate $\rm{u} - \rm{v}$ is $({{\textrm{u}_\textrm{1}},{\textrm{v}_\textrm{1}}\cos \mathrm{\beta} , - {\textrm{v}_\textrm{1}}\sin \mathrm{\beta} } )$.

As the plane is shifted upward 1 unit from the plane, the coordinate of point S in the coordinate $\rm{x} - \rm{y}$ is $({{\rm{u}_1},{\rm{v}_1}\cos \mathrm{\beta} ,1 - {\rm{v}_1}\sin \mathrm{\beta} } )$. At this point, the equation of the line $\rm{O}^{\prime}\rm{SP}$ can be expressed as:

$$\frac{\rm x}{{{\textrm{u}_\textrm{1}}}} = \frac{\rm y}{{{\textrm{v}_\textrm{1}}\cos \mathrm{\beta} }} = \frac{{\rm{z} + 1}}{{2 - {\textrm{v}_\textrm{1}}\sin \mathrm{\beta} }} = \rm{k}.$$

Since point P is on the sphere, it needs to satisfy the following equation:

$${{\rm x}^2} + {{\rm y}^2} + {{\rm{z}}^2} = 1.$$

Solving Eqs. (3) and (4),

$${\rm k} = \frac{{2({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{\rm{u}_{1}^{2} + \rm{v}_{1}^{2} - 4{\rm{v}_1}\sin \mathrm{\beta} + 4}}.$$

The vector of the output light ray is OP, considering the points S (u, v, 1) and P (x, y, z). Defining OP to be the vector m, it can then be expressed as

$${\bf m}({{\rm{u}_1},{\rm{v}_1}} )= \frac{1}{{\rm{u}_1^2 + \rm{v}_1^2 - 4{\rm{v}_1}\sin \mathrm{\beta} + 4}}\left( {\begin{array}{c} {2{\rm{u}_1}({2 - {\rm{v}_1}\sin \mathrm{\beta} } ),}\\ {2{\rm{v}_1}\cos \mathrm{\beta} ({2 - {\rm{v}_1}\sin \mathrm{\beta} } ),}\\ {{{({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}^2} - \rm{v}_1^2{{\cos }^2}\mathrm{\beta} - \rm{u}_1^2} \end{array}} \right).$$

As the energy of the output light is contained within the solid angle, we obtain

$${\rm d}\Omega = |{{{\bf m}_{{{\mathbf u}_{\mathbf 1}}}} \times ({{{\bf m}_{{{\mathbf v}_{\mathbf 1}}}}\cos \mathrm{\beta} } )} |{\rm d}{\rm{u}_1}{\rm d}{\rm{v}_1}.$$
where ${{\bf m}_{{\mathbf u}1}}$ and ${{\bf m}_{{\mathbf v}1}}$are the partial derivatives of ${\bf m}$ with respect to u and v, respectively. ${{\bf m}_{{\mathbf u}1}}$ and ${{\bf m}_{{\mathbf v}1}}$ can be expressed as
$$\begin{aligned} &{{\bf m}_{{\mathbf u1}}} = [ - \frac{{4{\rm{u}_1}^2({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}} + \frac{{2({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} }},\\ &- \frac{{4{\rm{u}_1}{\rm{v}_1}\cos \mathrm{\beta} ({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}},\\ &- \frac{{2{\rm{u}_1}}}{{4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} }} - \frac{{2{\rm{u}_1}({ - {\rm{u}_1}^2 - {\rm{v}_1}^2{{\cos }^2}\mathrm{\beta} + {{({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}^2}} )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}}]. \end{aligned}$$
$$\begin{aligned} &{{\bf m}_{{{\mathbf v}_{\mathbf 1}}}} = [ - \frac{{4{\rm{u}_1}({{\rm{v}_1} - 2\sin \mathrm{\beta} } )({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}} - \frac{{2{\rm{u}_1}\sin \mathrm{\beta} }}{{4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} }},\\ &- \frac{{4{\rm{v}_1}\cos \mathrm{\beta} ({{\rm{v}_1} - 2\sin \mathrm{\beta} } )({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}} + \frac{{4\cos \mathrm{\beta} ({1 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} }},\\ &- \frac{{2{\rm{v}_1}{{\cos }^2}\mathrm{\beta} + 2\sin \mathrm{\beta} ({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}}{{4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} }} - \frac{{2({{\rm{v}_1} - 2\sin \mathrm{\beta} } )({ - {\rm{u}_1}^2 - {\rm{v}_1}^2{{\cos }^2}\mathrm{\beta} + {{({2 - {\rm{v}_1}\sin \mathrm{\beta} } )}^2}} )}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}}]. \end{aligned}$$
using Eqs. (8) and (9), Eq. (7) then becomes
$${\rm d}\Omega = \frac{{8{{\cos }^2}\mathrm{\beta} |{ - 2 + {\rm{v}_1}\sin \mathrm{\beta} } |}}{{{{({4 + {\rm{u}_1}^2 + {\rm{v}_1}^2 - 4{\rm{v}_1}\sin \mathrm{\beta} } )}^2}}}{\rm d}{\rm{u}_1}{\rm d}{\rm{v}_1}.$$
when ${\mathrm{\beta}} = 0$,
$${\rm d}\Omega = {(1 + 0.25{\textrm{w}^2})^{ - 2}}{\rm d}\textrm{u}{\rm d}\textrm{v}\textrm{.}$$
where ${\textrm{w}^2} = {{\rm u}^2} + {{\rm v}^2}$. The directions of vectors OF and OP are aligned, and OP is the unit vector. If the distance between the source position O and the point F is $\mathrm{\rho}$, the vector OF can be expressed as:
$${\mathbf A} = {\mathrm{\rho}}{\bf m} = {\mathrm{\rho}}(\frac{\textrm{u}}{{1 + 0.25{\textrm{w}^2}}},\frac{\textrm{v}}{{1 + 0.25{\textrm{w}^2}}},\frac{{1 - 0.25{\textrm{w}^2}}}{{1 + 0.25{\textrm{w}^2}}}).$$

According to the law of reflection, the relationship between vectors A and T is shown in Fig. 4, T represents the vector OT; vector A of the emergent ray and the unit normal vector N at point F can be expressed as:

$$\frac{{\mathbf A}}{{|{\mathbf A} |}} + \frac{{{\mathbf A} - {\mathbf T}}}{{|{{\mathbf T} - {\mathbf A}} |}} = \frac{{2{\mathbf{N}(\mathbf{A}}\cdot {\mathbf N)}}}{{|{\mathbf A} |{{|{\mathbf N} |}^2}}}.$$
and
$$\begin{aligned} {\mathbf N} &= {{\bf A}_{\bf u}} \times {{\bf A}_{\bf v}} = ({{\mathrm{\rho}}_u}{\bf m} + {\mathrm{\rho}}{{\bf m}_{\bf u}}) \times ({{\mathrm{\rho}}_{\rm v}}{\bf m} + {\mathrm{\rho}}{{\bf m}_{\bf v}})\\ &= {{\mathrm{\rho}}^2}{(1 + 0.25{\textrm{w}^2})^{ - 2}}{\bf m} - {{\mathrm{\rho}}_{\rm u}}{\mathrm{\rho}}{{\bf m}_{\bf u}} - {{\mathrm{\rho}}_{\rm v}}{\mathrm{\rho}}{{\bf m}_{\bf v}}. \end{aligned}$$
where ${{\bf A}_{\mathbf u}}$ and ${{\bf A}_{\mathbf v}}$are the partial derivatives of ${\bf A}$ with respect to u and v, respectively, ${{\mathrm{\rho}} _{\mathbf u}}$ and ${{\mathrm{\rho}} _{\mathbf v}}$ are the partial derivatives of $\mathrm{\rho}$ with respect to u and v, respectively, and F and G can be defined as follows:
$$\begin{aligned} &{\bf F}\textrm{ = }|{{\bf T} - {\mathbf A}} |{(1 + 0.25{\textrm{w}^2})^{ - 1}}\\ &{[{{{\mathrm{\rho}}}^2}{(1 + 0.25{\textrm{w}^2})^{ - 2}} + {{{\mathrm{\rho}}}_{\rm u}}^2 + {{{\mathrm{\rho}}}_{\rm v}}^2]^{ - 1}}. \end{aligned}$$
$$\begin{aligned} &{\bf G} = {(1 + 0.25{\textrm{w}^2})^{ - 1}} + {\bf F}[ - {{{\mathrm{\rho}}}^2}{(1 + 0.25{\textrm{w}^2})^{ - 2}}\\ &+ {{{\mathrm{\rho}}}_{\rm u}}^2 + {{{\mathrm{\rho}}}_{\rm v}}^2 - {{\mathrm{\rho}}}({{{\mathrm{\rho}}}_{\rm u}}\textrm{u} + {{{\mathrm{\rho}}}_{\rm v}}\textrm{v}){(1 + 0.25{\textrm{w}^2})^{ - 1}}]. \end{aligned}$$

Finally, the coordinates of point T can be calculated as follows:

$${\rm{t}_{\rm x}} = \textrm{u}{\bf G} + 2{{{\mathrm{\rho}}}_{\rm u}}{{{\mathrm{\rho}}}_{\rm v}}{\bf F,}\,{\rm{t}_{\rm y}} = \textrm{v}{\bf G} + 2{{{\mathrm{\rho}}}_{\rm u}}{{{\mathrm{\rho}}}_{\rm v}}{\bf F}.$$

This equation can be rewritten as:

$${\rm{t}_{\rm x}} = {\rm{t}_{\rm x}}(\textrm{u},\textrm{v},{{\mathrm{\rho}}},{{{\mathrm{\rho}}}_{\rm u}},{{{\mathrm{\rho}}}_{\rm v}}),\,{\rm{t}_{\rm y}} = {\rm{t}_{\rm y}}(\mathrm{u,v,\mathrm{\rho} ,}{{{\mathrm{\rho}}}_\textrm{u}}\textrm{,}{{{\mathrm{\rho}}}_\textrm{v}}).$$

Based on the law of energy conservation, the relationship between the total energy of the outgoing light and the total energy of the target surface irradiation region can be obtained as follows:

$$\int\!\!\!\int\limits_{{\Omega _1}} {\textrm{E}({{\rm{t}_{\rm x}},{\rm{t}_{\rm y}}} )\rm{d}{\rm{t}_{\rm x}}} \rm{d}{\rm{t}_{\rm y}} = \int\!\!\!\int\limits_{{\Omega _{}}} {I({\textrm{u},\textrm{v}} )d} \Omega .$$
where $\Omega $ is the solid angle of the outgoing light. According to the mathematical model established by Schruben,
$$\rm{d}{\rm{t}_{\rm x}}\rm{d}{\rm{t}_{\rm y}} = |{J(T)} |d\textrm{u}d\textrm{v},|{\textrm{J}(T)} |= \left|{\begin{array}{cc} {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{\mathrm{\partial} \textrm{u}}}}&{\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{\mathrm{\partial} \textrm{v}}}}\\ {\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{\mathrm{\partial} \textrm{u}}}}&{\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{\mathrm{\partial} \textrm{v}}}} \end{array}} \right|.$$

Then Eq. (19) can be rewritten:

$$\textrm{E}({{\rm{t}_{\rm x}},{\rm{t}_{\rm y}}} )\left|{\begin{array}{cc} {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{\mathrm{\partial} \textrm{u}}}}&{\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{\mathrm{\partial} \textrm{v}}}}\\ {\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{\mathrm{\partial} \textrm{u}}}}&{\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{\mathrm{\partial} \textrm{v}}}} \end{array}} \right|\rm{d}{\rm{t}_{\rm x}}\rm{d}{\rm{t}_{\rm y}} = I({\textrm{u},\textrm{v}} )({1 + 0.25{\textrm{w}^2}^{ - 2}} )\rm{d}{\rm{t}_{\rm x}}\rm{d}{\rm{t}_{\rm y}}.$$

Based on the optimal transmission theory, the MA equation of the system can be obtained as follows:

$${\rm{a}_1}({{{{\mathrm{\rho}}}_{\rm{uu}}}{{{\mathrm{\rho}}}_{\rm{vv}}} - {{\mathrm{\rho}}}_{\rm{uv}}^2} )+ {\rm{a}_2}{{{\mathrm{\rho}}}_{\rm{uu}}} + {\rm{a}_3}{{{\mathrm{\rho}}}_{\rm{uv}}} + {\rm{a}_4}{{{\mathrm{\rho}}}_{\rm{vv}}} + {\rm{a}_5} = 0.$$
$$\begin{aligned} {{\rm a}_1} &= \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_\textrm{u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}\\ {{\rm a}_2} &= \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{v}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{v}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}} + {{{\mathrm{\rho}}}_{\rm v}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}} \right)\\ {{\rm a}_3} &= \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{v}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{v}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}} + \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{u}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{u}} + {{{\mathrm{\rho}}}_{\rm u}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm u}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}}} \right) + {{{\mathrm{\rho}}}_{\rm v}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}} \right)\\ {{\rm a}_4} &= \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{u}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{u}} + {{{\mathrm{\rho}}}_{\rm u}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{{{\mathrm{\rho}}}_{\rm v}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}}} \right)\\ {{\rm a}_5} &= \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{u}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{v}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{v}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{u}} + {{{\mathrm{\rho}}}_{\rm u}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{v}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{v}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}}} \right) + {{{\mathrm{\rho}}}_{\rm u}}\left( {\frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{\textrm{u}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{{{\mathrm{\rho}}}} - \frac{{\mathrm{\partial} {\rm{t}_{\rm x}}}}{{{\mathrm{\rho}}}}\frac{{\mathrm{\partial} {\rm{t}_{\rm y}}}}{\textrm{u}}} \right). \end{aligned}$$
where ${{\mathrm{\rho}} _{{\mathbf u}{\bf u}}}$ and ${{\mathrm{\rho}} _{{\bf vv}}}$are the second partial derivative of $\mathrm{\rho}$ with respect to u and v, respectively, ${{\mathrm{\rho}} _{{\mathbf u}{\bf v}}}$ is the mixed partial derivative of $\mathrm{\rho}$ with respect to u and v. Thus, the mathematical model is completed.
$$\left\{ {\begin{array}{c} {{{\rm a}_1}({{{\mathrm{\rho}}}_{\rm{uu}}}{{{\mathrm{\rho}}}_{{\rm{vv}}}} - {{\mathrm{\rho}}}_{{\rm{uv}}}^2) + {{\rm a}_2}{{{\mathrm{\rho}}}_{{\rm{uu}}}} + {{\rm a}_3}{{{\mathrm{\rho}}}_{{\rm{uv}}}} + {{\rm a}_4}{{{\mathrm{\rho}}}_{{\rm{vv}}}} + {{\rm a}_5} = 0}\\ {\textrm{BC}:\left\{ {\begin{array}{{c}} {{\rm{t}_{\rm x}} = {\rm{t}_{\rm x}}(\textrm{u},\textrm{v},{{\mathrm{\rho}}},{{{\mathrm{\rho}}}_{\rm u}},{{{\mathrm{\rho}}}_{\rm v}})}\\ {{\rm{t}_{\rm y}} = {\rm{t}_{\rm y}}(\textrm{u},\textrm{v},{{\mathrm{\rho}}},{{{\mathrm{\rho}}}_{\rm u}},{{{\mathrm{\rho}}}_{\rm v}})} \end{array}} \right.\begin{array}{{ccc}} {}&:&{\mathrm{\partial} {\Omega _0} \to \mathrm{\partial} {\Omega _1}} \end{array}.} \end{array}} \right.$$

 figure: Fig. 4.

Fig. 4. .Law of reflection.

Download Full Size | PDF

3.2 Chebyshev point and Chebyshev grid

In a freeform reflector illumination system, when the light source has rotational symmetry, the projection area of the light source is circular, whereas the sampling grid is usually topological squares or rectangles [26,27]; this results in a mismatch between the two. Therefore, in this paper, the method of discretization of the circular region is proposed. The specific implementation methods are as follows: 1) The system is designed by using the mapping from target irradiation to light source emission; 2) The boundary condition is the mapping between the boundary of the irradiation output region and the boundary of the light source emission region; 3) Select the luminous region of the light source, and set its projection boundary to match the mapping boundary. On this basis, the solution of the mathematical model of the optical system is obtained, and the problem of mismatch between the projected area of the light source and the sampling grid is solved. When the area of the mirror is large, 41 × 41 or 81 × 81 points can be used as the discretization condition, as shown in Fig. 5.

 figure: Fig. 5.

Fig. 5. Chebyshev points on freeform reflector: (a) 41 × 41 points (b) 81 × 81 points.

Download Full Size | PDF

Here, only a numerical solution can be used; hence, it is impossible to obtain an accurate answer. First, we must discretize the elliptic Monge-Ampère equation and the nonlinear boundary condition [28]. ${\rm{h}_{\rm i}} = {{\mathrm{\rho}} _{\rm{i} + 1,\rm{j}}} - {{\mathrm{\rho}} _{\rm{i},\rm{j}}}$ and ${\rm{h}_{\rm j}} = {{\mathrm{\rho}} _{\rm{i},{\rm j} + 1}} - {{\mathrm{\rho}} _{\rm{i},\rm{j}}}$ are the spacings in u- and v-axes, respectively, as shown in Fig. 6. The grid points along ${\Omega _0}$ are called the boundary points, and the grid points inside ${\Omega _0}$ are the interior points. Therefore, each interior point should satisfy the elliptic Monge-Ampère equation, and each boundary point should satisfy the boundary condition. Subsequently, we use the difference formula for derivatives to discretize the elliptic Monge-Ampère equation and the nonlinear boundary condition.

 figure: Fig. 6.

Fig. 6. Nine-point finite difference scheme for the interior points and boundary points.

Download Full Size | PDF

For the interior point, we define the scheme shown in Fig. 6 as:

$$\begin{aligned} {{\mathrm{\rho}} _{\rm u}} &= \frac{{{{\mathrm{\rho}} _{\rm{i} + 1,\rm{j}}} - {{\mathrm{\rho}} _{\rm{i} - 1,\rm{j}}}}}{{2{\rm{h}_{\rm i}}}}\textrm{,}{{\mathrm{\rho}} _{\rm v}} = \frac{{{{\mathrm{\rho}} _{\rm{i},{\rm j} + 1}} - {{\mathrm{\rho}} _{\rm{i},{\rm j} - 1}}}}{{2{\rm{h}_{\rm j}}}},{{\mathrm{\rho}} _{{\rm{uu}}}} = \frac{{{{\mathrm{\rho}} _{\rm{i} + 1,\rm{j}}} - \mathrm{\rho} {\rm{z}_{\rm{i},\rm{j}}} + {{\mathrm{\rho}} _{\rm{i} - 1,\rm{j}}}}}{{\rm{h}_{\rm i}^2}}\\ {{\mathrm{\rho}} _{{\rm{vv}}}} &= \frac{{{{\mathrm{\rho}} _{\rm{i},{\rm j} + 1}} - 2{{\mathrm{\rho}} _{\rm{i},\rm{j}}} + {{\mathrm{\rho}} _{\rm{i},{\rm j} - 1}}}}{{h_j^2}}\textrm{,}{{\mathrm{\rho}} _{{\rm{uv}}}} = \frac{{{{\mathrm{\rho}} _{\rm{i} + 1,{\rm j} + 1}} - {{\mathrm{\rho}} _{\rm{i} + 1,{\rm j} - 1}} - {{\mathrm{\rho}} _{\rm{i} - 1,{\rm j} + 1}} + {{\mathrm{\rho}} _{\rm{i} - 1,{\rm j} - 1}}}}{{4{\rm{h}_{\rm i}}{\rm{h}_{\rm j}}}}. \end{aligned}$$

Taking ${\rm u} = {{\rm u}_{\max }}$ as an example of the boundary line, as shown in Fig. 6, by one-sided finite differences, we find that:

$${{\mathrm{\rho}} _{\rm u}} = \frac{{3{{\mathrm{\rho}} _{\rm{m},\rm{j}}} - 4{{\mathrm{\rho}} _{\rm{m} - 1,\rm{j}}} + {{\mathrm{\rho}} _{\rm{m} - 2,\rm{j}}}}}{{2{\rm{h}_{\rm i}}}},{{\mathrm{\rho}} _{\rm v}} = \frac{{{{\mathrm{\rho}} _{\rm{m},{\rm j} + 1}} - {{\mathrm{\rho}} _{\rm{m},{\rm j} - 1}}}}{{2{\rm{h}_{\rm j}}}}.$$

Equation (22) can be rewritten in (m + 1) (n + 1) equation sets and there are (m + 1) (n + 1) Chebyshev points. The intensity of the point source obeys the lambert characteristic, the intensity $\rm{I}(\rm{u,v}){(1 + 0.25{\rm{w}^2})^{ - 2}}$can then be obtained. Solving the (m + 1) (n + 1) equation sets, the coordinates of the points on the freeform reflector can be determined. Newton downhill method is used to solve all the remaining points. Finally, the freeform surface can be constructed, and the Bezier surface pieces proposed by Wu et al. [29] are used to suture each piece to obtain the final freeform surface model.

4. Simulation and testing

4.1 Off-axis reflected illumination model

Based on the freeform surface model establishment and solution method proposed in the previous section, we designed the SCAR-based radiation calibration system freeform reflector in detail. The design process is shown in Fig. 7.

 figure: Fig. 7.

Fig. 7. Design flowchart.

Download Full Size | PDF

The structural parameters of each part of the radiometric calibration light source system are shown in Table 1.

Tables Icon

Table 1. Design parameters of the freeform reflector.

Here θ is the luminous angle of the light source, β is the deflection angle of the freeform reflector, γ is the tilt angle of the diffuser, r is the distance between the light source and the freeform reflector, Z is the distance between the target surface and the freeform reflector, and F# is the F-number of the system, which is represented by the ratio of the distance between the reflector and the light source and the maximum aperture of the reflector. The discrete point data are obtained by solving using MATLAB. At this time, the discrete graph model of the freeform reflector is obtained by using Chebyshev points, as shown in Fig. 8; 81 × 81 points are used. Equation (5) is solved to obtain the freeform reflector model, as shown in Fig. 9.

 figure: Fig. 8.

Fig. 8. Schematic of Chebyshev points.

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. Schematic of the freeform reflector simulation.

Download Full Size | PDF

The results are imported into the LightTools optical software. The optical simulation is carried out according to the optical path of the whole spectral calibration light source, the schematic of the illumination is shown in Fig. 10. It is worth noting that the previous derivations are based on point light source, but the light source in Fig. 10 is a circular outlet with a diameter of 22 mm of a small integrating sphere. If an extended light source is used to replace the point light source, the irradiance of the boundary region on the target surface will show a large trapezoidal decline, as shown in Fig. 11(a). Therefore, the freeform reflector solved based on the point source model cannot be directly applied to the actual calibration task. Here we apply a feedback optimization method which we have used before to optimize the freeform reflector. The feasibility of this method has been demonstrated previously in coaxial systems [30].

 figure: Fig. 10.

Fig. 10. The model of light path.

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. A square illuminated area with a side length of 150 mm: (a) before optimization (b) after optimization.

Download Full Size | PDF

The irradiance simulation result of the target plane after optimized is shown in Fig. 11(b), it can be seen that the illumination area formed by the freeform reflector on the target plane is larger than 120 mm × 120 mm which is sufficient to cover the diffuser and is of good uniformity. The simulation results show that the freeform reflector design method based on Chebyshev points can achieve ideal uniformity and illumination area, and it is feasible to be applied to turn and collimate the optical path in the SCAR-based radiation calibration system. To verify this result, detailed tests on the freeform reflector and calibration light source system were conducted.

4.2 Actual test

According to the above freeform reflector simulation model established based on the off-axis reflection illumination model of the SCAR-based radiation calibration system, the freeform reflector structure model was constructed as shown in Fig. 12(a). The surface of the model is smooth and continuous, which results in good machinability. Moore Nanotech 250UPL precision single-point diamond machining machine was used to process the freeform reflector, and the real object is shown in Fig. 12(b). The whole freeform reflector is made of aluminum. In order not to damage the surface of the reflector, the surface is not coated. And then a series of tests were conducted on the roughness of the freeform reflector and the optical properties of the freeform calibration light source system.

 figure: Fig. 12.

Fig. 12. (a) Freeform reflector model pattern, (b) freeform reflector physical pattern.

Download Full Size | PDF

LUPHO Scan scanner was used to detect the surface shape of the freeform reflector and obtained the roughness data. The test results are shown in Fig. 13, where the RMS is 0.61 μm. The test results show that the processed freeform reflector has no obvious tool marks and the continuity of the processed surface is good, which can meet the practical application.

 figure: Fig. 13.

Fig. 13. Test results of surface roughness of the freeform reflector.

Download Full Size | PDF

As shown in Fig. 14, the freeform reflector was integrated into the SCAR-based radiation calibration system for light path bending and collimation. As expected before, after the bending and collimation of the freeform reflector, a light spot is formed on the diffuser, which is sufficient to cover the 120 mm × 120 mm diffuser. The uniformity of the spot will be evaluated by experimental tests.

 figure: Fig. 14.

Fig. 14. Uniform whole spectral calibration light source.

Download Full Size | PDF

As mentioned before, the size of the light spot created by the freeform reflector on the diffuser needs to be large enough because the diffuser needs to cover the observation area of EMIS and TR, and the irradiance uniformity of the light spot at the diffuser position (Target Plane) is supposed to be better than 98% in the 100 mm × 100 mm region. The location of the target plane is shown in Fig. 15. In this regard, the spot size and its uniformity were tested in detail.

 figure: Fig. 15.

Fig. 15. The location of the target plane.

Download Full Size | PDF

The grid-scanning method is used to detect the light spot formed by the freeform reflector on the target plane, and the scheme is shown in Fig. 16. We fixed the Hamamatsu S2281 Si-photodiode on a electric two-dimensional guide rail, and placed the two-dimensional guide rail on the position of the diffuser to drive S2281 to perform point-by-point scanning in the target plane, this is used to measure the irradiance in each small area received by the diffuser. After connecting the detector to an amplification box and a calibrated 6½ digit digital multimeter, three groups of data were measured at each measuring point and the average value was taken as the irradiance data at that point. The scan area is 150 mm × 150 mm, and the scan step is 10 mm.

 figure: Fig. 16.

Fig. 16. Irradiance uniformity test scheme.

Download Full Size | PDF

Denote the relative position of the test point to the center of the diffuser in two directions by X and Y, respectively. The normalized test result of irradiance is shown in Fig. 17. It can be seen from the data that due to the smoothness and continuity of the freeform reflector, the irradiance distribution of its reflected light on the target plane is also continuous without obvious mutation point. The diameter of the photosensitive surface of S2281 is 11.3 mm, so a 10 × 10 scanning grid can cover the test area of 100 mm × 100 mm. As shown in Table 2, we extracted 10 × 10 data points in the center of the scanned grid for analysis, Eq. (27) is used to calculate the irradiance nonuniformity u, where Emax is the maximum irradiance value and Emin is the minimum irradiance value in grid data. After calculation, the irradiance uniformity in the test area is 98.2%, which meets the requirement of uniformity better than 98%.

$${\rm u} = \frac{{{\rm{E}_{\max }} - {{\rm E}_{\min }}}}{{{{\rm E}_{\max }} + {{\rm E}_{\min }}}}$$

In addition, to verify the radiance uniformity of the freeform calibration light source system, we use SVC HR1024i spectrometer to measure the radiance in the way of grid scanning, the test diagram is shown in Fig. 18. The SVC was placed in the same observation direction as the EMIS, and was driven by a two-dimensional guide rail for point-by-point scanning. The diameter of the SVC observation area on the diffuser is about 45 mm, so a 5 by 5 grid test with a step size of 15 mm is sufficient to cover the 100 mm × 100 mm illuminated area.

 figure: Fig. 17.

Fig. 17. Normalized test data for irradiance uniformity.

Download Full Size | PDF

 figure: Fig. 18.

Fig. 18. Radiance uniformity test scheme.

Download Full Size | PDF

Tables Icon

Table 2. Normalized test data for irradiance.

The test data of two spectral channels were selected of SVC for uniformity analysis. As shown in Fig. 19(a) and Fig. 19(b) is normalized radiance test data at 1036.9 nm and 780.8 nm, respectively. Equation (27) was also used to evaluate the uniformity of grid test data. After calculation, the uniformity of 1036.9 nm and 780.8 nm are both 98.5%, which meets the requirement of uniformity better than 2%.

 figure: Fig. 19.

Fig. 19. Normalized test data for radiance: (a) 1036.9 nm (b) 780.8 nm.

Download Full Size | PDF

5. Discussion

At present, the radiometric calibration of space remote sensing is faced with difficult problems such as poor consistency, poor data comparability of different remote sensors, and low radiometric calibration accuracy, this has greatly affected our monitoring and forecasting of the Earth's climate and environmental changes. The SCAR-based onboard radiation calibration system can greatly improve the credibility of space remote sensing data in both temporal and spatial dimensions. However, the load, space, and energy consumption of the calibration system are the factors that restrict its final performance on board. A freeform reflector can precisely control the light to reach a specified area, this can not only make the space layout of the whole calibration system more flexible and increase the space utilization but also improve the energy utilization and reduce the energy consumption. Most importantly, it can provide a large area of uniform light source for the calibration system, reduce the error introduced by the radiation transfer link, and improve the accuracy of the onboard calibration. Therefore, compared with the current illumination scheme, the onboard illumination mode based on the freeform reflector adopted in this paper more perfectly meets the requirements of onboard calibration in terms of both space, energy utilization and uniform light source. The freeform reflector light source has the potential to be used in more onboard calibration systems.

The Chebyshev spectrum method can be used for differential equations with boundary conditions and computation intervals showing non-periodicity. The initial structure discretization design method of freeform surface based on Chebyshev points proposed in this paper, starting from Chebyshev's basic definition, it is extended and derived in two dimensions, so that it can better match the projection grid of the light source. The advantage of this method lies in that not only the emission region of the light source is effectively discretized, but also the partial differential equation is decomposed into the ordinary differential equation by this method, and the solution of complex (MA) Monge-Ampère equation is transformed into solving ordinary differential equations, this also makes the design of freeform mirror easier and more universal. The freeform reflectors established based on this idea are more smooth and more continuous. As the core device of the calibration light source system, it provides important support for irradiance uniformity of reflected light.

6. Conclusion

This paper proposed a new onboard calibration light source system based on a freeform reflector for the SCAR-based radiation calibration system. A new design idea for the freeform surface is also put forward according to the present freeform surface design method. Based on the characteristics of Chebyshev meshes, the topological meshes of the outgoing light can be discretized inhomogeneous. At the same time, the freeform surface model with a good initial structure is directly generated by the Optimal mass Transport (OMT) theory. Finally, the Bessel curve is used to control and adjust the discrete points, and a continuous smooth freeform surface is constructed. The constructed freeform reflector is simulated and analyzed according to the light path of the SCAR-based radiation calibration system, the size and uniformity of the spot formed in the target plane can meet the application requirements. Then we processed the designed freeform reflector and tested its surface roughness, the result shows that the RMS is 0.61 μm and the continuity of the machining surface is good. We integrate it into the calibration system and use S2281 Si-photodiode and SVC spectrometer respectively to measure the irradiance of the effective illumination area in the target plane and the radiance of the diffuser with grid-scanning method. The experimental results show that in the effective illumination area of 100 mm × 100 mm of the target plane, the irradiance uniformity of the diffuser and the light emitted from the freeform reflector are 98.5% and 98.2%, respectively, which are both better than 98% and meet the practical application requirements. The designed freeform reflector can be perfectly applied to the optical path bending and collimation of the SCAR-based radiation calibration system. The freeform reflector calibration light source system not only provides high uniformity calibration light source for onboard radiometric calibration, but also improves the space utilization of the calibration system, reduces on-board resource consumption, and also provides the possibility to achieve the goal of the spectral radiance measurement uncertainty better than 1% on the reflected solar band of Chinese Space-based Radiometric Benchmark (CSRB).

Funding

Strategic Priority Research Program of the Chinese Academy of Sciences (XDA28050103); National Key Research and Development Program of China (2018YFB0504600, 2018YFB0504603).

Acknowledgments

XinYe thanks Strategic Priority Research Program of the Chinese Academy of Sciences and National Key R&D Program of China under grant for help identifying collaborators for this work.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. J. Xiong, C. Cao, and G. Chander, “An overview of sensor calibration inter-comparison and applications,” Front. Earth Sci. China 4(2), 237–252 (2010). [CrossRef]  

2. X. Q. Hu, Y. X. Zhang, Z. Q. Liu, G. S. Zhang, Y. B. Huang, K. M. Qiu, Y. K. Wang, L. J. Zhang, X. B. Zhu, and Z. G. Rong, “Optical characteristics of China Radiometric Calibration Site for Remote Sensing Satellite Sensors (CRCSRSSS),” Proc. SPIE 4151, 77–86 (2001). [CrossRef]  

3. Y. P. Wang, X. Q. Hu, H. R. Wang, X. Ye, and W. Fang, “Standard transfer chain for radiometric calibration of optical sensing instruments with traceability in solar reflective bands,” Opt. Precis. Eng. 23(7), 1807–1812 (2015). [CrossRef]  

4. P. Zhang, N. Lu, C. Li, L. Ding, and J. Schmetz, “Development of the Chinese space-based radiometric benchmark mission LIBRA,” Remote Sens. 12(14), 2179 (2020). [CrossRef]  

5. N. Fox, A. K. Weiss, W. Schmutz, K. Thome, D. Young, B. A. Wielichi, R. Winkler, and E. Woolliams, “Accurate radiometry from space: An essential tool for climate studies,” Phil. Trans. R. Soc. A 369(1953), 4028–4063 (2011). [CrossRef]  

6. W. Sun and C. Lukashin, “Modeling polarized solar radiation from the ocean - atmosphere system for CLARREO inter-calibration applications,” Atmos. Chem. Phys. 13(20), 10303–10324 (2013). [CrossRef]  

7. X. Ye, X. L. Yi, C. Lin, W. Fang, K. Wang, Z. W. Xia, Z. H. Ji, Y. Q. Zheng, D. Sun, and J. Quan, “Instrument Development: Chinese Radiometric Benchmark of Reflected Solar Band Based on Space Cryogenic Absolute Radiometer,” Remote Sens. 12(17), 2856 (2020). [CrossRef]  

8. X. Ye, X. L. Yi, W. Fang, K. Wang, Y. Luo, Z. W. Xia, and Y. P. Wang, “Design and Investigation of Absolute Radiance Calibration Primary Radiometer,” IET Sci. Meas. Technol. 12(8), 994–1000 (2018). [CrossRef]  

9. K. C. Lei, X. Ye, Z. W. Xia, N. Xu, S. Q. Li, Y. C. Zhang, Y. W. Wang, Z. W. Liu, and Z. G. Li, “Design of a Compact Transfer Radiometer for a Radiometric Benchmark Transfer Chain,” Sensors 22(18), 6795 (2022). [CrossRef]  

10. W. W. Pang, X. B. Zheng, J. J. Li, and X. S. Shi, “High Accuracy Calibration Technology of Remote Sensor Tracing to Cryogenic Radiometer,” J. Atmos. Environ. Opt. 9(2), 138 (2014). [CrossRef]  

11. B. Y. Li, J. W. Ren, Z. Wan, X. S. Li, and J. X. Sun, “Research of large aperture integrating sphere used in the radiative calibration for space remote sensor,” J. Optoelectron. Laser 24(3), 464–469 (2013).

12. L. Hanssen, “Integrating-sphere system and method for absolute measurement of transmittance, reflectance, and absorptance of specular samples,” Appl. Opt. 40(19), 3196–3204 (2001). [CrossRef]  

13. A. Rabl and J. M. Gordon, “Reflector design for illumination with extended sources: the basic solutions,” Appl. Opt. 33(25), 6012–6021 (1994). [CrossRef]  

14. P. Davies, “Edge-Ray Principle Of Non-imaging Optics,” J. Opt. Soc. Am. A 11(4), 1256–1259 (1994). [CrossRef]  

15. T. L. R. Davenport, T. A. Hougha, and W. J. Cassarlya, “Optimi-zation for illumination systems: the next level of design,” Proc. SPIE 5456, 81–90 (2004). [CrossRef]  

16. Y. Luo, Z. Feng, Y. Han, and H. Li, “Design of compact and smooth freeform optical system with uniform illuminance for LED source,” Opt. Express 18(9), 9055–9063 (2010). [CrossRef]  

17. J. Schruben, “Formulation of a reflector-design problem for a lighting fixture,” J. Opt. Soc. Am. 62(12), 1498–1501 (1972). [CrossRef]  

18. C. R. Prins, J. H. M. Ten Thije Boonkkamp, J. Van Roosmalen, L. W. Ijzerman, and W. T. Tukker, “A Monge-Ampère-Solver for freeform Reflector Design,” SIAM J. Sci. Comput. 36(3), B640–B660 (2014). [CrossRef]  

19. T. Yang, Y. Duan, D. Cheng, and Y. Wang, “Freeform imaging optical system design: theories, development and applications,” Acta Opt. Sin. 41(1), 0108001 (2021). [CrossRef]  

20. V. Oliker, “Geometric and variational methods in optical design of reflecting surfaces with prescribed illuminance properties,” Proc. SPIE 5942, 594207 (2005). [CrossRef]  

21. V. Oliker, “Optical design of freeform ywo-mirror beam-shaping system,” J. Opt. Soc. Am. A 24(12), 3741–3752 (2007). [CrossRef]  

22. R. M. Wu, S. Q. Chang, Z. R. Zheng, L. F. Zhao, and X. Liu, “Formulating the design of two freeform lens surfaces for point-like light sources,” Opt. Lett. 43(7), 1619–1622 (2018). [CrossRef]  

23. R. M. Wu and H. Hua, “Direct design of aspherical lenses for extended non-Lambertian sources in three-dimensional rotational geometry,” Opt. Express 24(2), 1017–1030 (2016). [CrossRef]  

24. R. M. Wu, L. Xu, P. Liu, Y. Q. Zhang, Z. R. Zheng, H. F. Li, and X. Liu, “Freeform illumination design: a nonlinear boundary problem for the elliptic Monge–Ampére equation,” Opt. Lett. 38(2), 229–231 (2013). [CrossRef]  

25. L. Wang, K. Y. Qian, and L. Yi, “Discontinuous free-form lens design for prescribed irradiance,” Appl. Opt. 46(18), 3716–3723 (2007). [CrossRef]  

26. X. Mao, H. Li, Y. Han, and Y. Luo, “Polar-grids based source-target mapping construction method for designing freeform illumination system for a lighting target with arbitrary shape,” Opt. Express 23(4), 4313–4328 (2015). [CrossRef]  

27. J. Chaves, Introduction To Nonimaging Optics (CRC press, 2008).

28. M. J. H. Anthonissen, L. B. Romijn, J. H. M. Ten Thije Boonkkamp, and W. L. IJzerman, “Unified mathematical framework for a class of fundamental freeform optical systems,” Opt. Express 29(20), 31650–31664 (2021). [CrossRef]  

29. R. M. Wu, Z. R. Zheng, H. F. Li, and X. Liu, “Constructing optical freeform surfaces using unit tangent vectors of feature data points,” J. Opt. Soc. Am. A 28(9), 1880–1888 (2011). [CrossRef]  

30. J. W. Tian, X. Ye, and W. Fang, “Freeform reflector design for uniform illumination in a space-based remote sensor calibration system,” Appl. Opt. 61(16), 4678–4686 (2022). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (19)

Fig. 1.
Fig. 1. Light path of the whole spectral calibration light source system of SCAR-based radiation calibration.
Fig. 2.
Fig. 2. Diagram of energy transfer for a freeform reflector.
Fig. 3.
Fig. 3. Diagram of the illumination coordinate system: (a) 3D (b) 2D.
Fig. 4.
Fig. 4. .Law of reflection.
Fig. 5.
Fig. 5. Chebyshev points on freeform reflector: (a) 41 × 41 points (b) 81 × 81 points.
Fig. 6.
Fig. 6. Nine-point finite difference scheme for the interior points and boundary points.
Fig. 7.
Fig. 7. Design flowchart.
Fig. 8.
Fig. 8. Schematic of Chebyshev points.
Fig. 9.
Fig. 9. Schematic of the freeform reflector simulation.
Fig. 10.
Fig. 10. The model of light path.
Fig. 11.
Fig. 11. A square illuminated area with a side length of 150 mm: (a) before optimization (b) after optimization.
Fig. 12.
Fig. 12. (a) Freeform reflector model pattern, (b) freeform reflector physical pattern.
Fig. 13.
Fig. 13. Test results of surface roughness of the freeform reflector.
Fig. 14.
Fig. 14. Uniform whole spectral calibration light source.
Fig. 15.
Fig. 15. The location of the target plane.
Fig. 16.
Fig. 16. Irradiance uniformity test scheme.
Fig. 17.
Fig. 17. Normalized test data for irradiance uniformity.
Fig. 18.
Fig. 18. Radiance uniformity test scheme.
Fig. 19.
Fig. 19. Normalized test data for radiance: (a) 1036.9 nm (b) 780.8 nm.

Tables (2)

Tables Icon

Table 1. Design parameters of the freeform reflector.

Tables Icon

Table 2. Normalized test data for irradiance.

Equations (27)

Equations on this page are rendered with MathJax. Learn more.

Ω 1 E( t x , t y ) d t x d t y = Ω I ( u , v ) d Ω .
[ u v w ] = [ 1 0 0 0 cos β sin β 0 sin β cos β ] [ u 1 v 1 0 ] .
x u 1 = y v 1 cos β = z + 1 2 v 1 sin β = k .
x 2 + y 2 + z 2 = 1.
k = 2 ( 2 v 1 sin β ) u 1 2 + v 1 2 4 v 1 sin β + 4 .
m ( u 1 , v 1 ) = 1 u 1 2 + v 1 2 4 v 1 sin β + 4 ( 2 u 1 ( 2 v 1 sin β ) , 2 v 1 cos β ( 2 v 1 sin β ) , ( 2 v 1 sin β ) 2 v 1 2 cos 2 β u 1 2 ) .
d Ω = | m u 1 × ( m v 1 cos β ) | d u 1 d v 1 .
m u 1 = [ 4 u 1 2 ( 2 v 1 sin β ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 + 2 ( 2 v 1 sin β ) 4 + u 1 2 + v 1 2 4 v 1 sin β , 4 u 1 v 1 cos β ( 2 v 1 sin β ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 , 2 u 1 4 + u 1 2 + v 1 2 4 v 1 sin β 2 u 1 ( u 1 2 v 1 2 cos 2 β + ( 2 v 1 sin β ) 2 ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 ] .
m v 1 = [ 4 u 1 ( v 1 2 sin β ) ( 2 v 1 sin β ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 2 u 1 sin β 4 + u 1 2 + v 1 2 4 v 1 sin β , 4 v 1 cos β ( v 1 2 sin β ) ( 2 v 1 sin β ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 + 4 cos β ( 1 v 1 sin β ) 4 + u 1 2 + v 1 2 4 v 1 sin β , 2 v 1 cos 2 β + 2 sin β ( 2 v 1 sin β ) 4 + u 1 2 + v 1 2 4 v 1 sin β 2 ( v 1 2 sin β ) ( u 1 2 v 1 2 cos 2 β + ( 2 v 1 sin β ) 2 ) ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 ] .
d Ω = 8 cos 2 β | 2 + v 1 sin β | ( 4 + u 1 2 + v 1 2 4 v 1 sin β ) 2 d u 1 d v 1 .
d Ω = ( 1 + 0.25 w 2 ) 2 d u d v .
A = ρ m = ρ ( u 1 + 0.25 w 2 , v 1 + 0.25 w 2 , 1 0.25 w 2 1 + 0.25 w 2 ) .
A | A | + A T | T A | = 2 N ( A N ) | A | | N | 2 .
N = A u × A v = ( ρ u m + ρ m u ) × ( ρ v m + ρ m v ) = ρ 2 ( 1 + 0.25 w 2 ) 2 m ρ u ρ m u ρ v ρ m v .
F  =  | T A | ( 1 + 0.25 w 2 ) 1 [ ρ 2 ( 1 + 0.25 w 2 ) 2 + ρ u 2 + ρ v 2 ] 1 .
G = ( 1 + 0.25 w 2 ) 1 + F [ ρ 2 ( 1 + 0.25 w 2 ) 2 + ρ u 2 + ρ v 2 ρ ( ρ u u + ρ v v ) ( 1 + 0.25 w 2 ) 1 ] .
t x = u G + 2 ρ u ρ v F , t y = v G + 2 ρ u ρ v F .
t x = t x ( u , v , ρ , ρ u , ρ v ) , t y = t y ( u , v , ρ , ρ u , ρ v ) .
Ω 1 E ( t x , t y ) d t x d t y = Ω I ( u , v ) d Ω .
d t x d t y = | J ( T ) | d u d v , | J ( T ) | = | t x u t x v t y u t y v | .
E ( t x , t y ) | t x u t x v t y u t y v | d t x d t y = I ( u , v ) ( 1 + 0.25 w 2 2 ) d t x d t y .
a 1 ( ρ u u ρ v v ρ u v 2 ) + a 2 ρ u u + a 3 ρ u v + a 4 ρ v v + a 5 = 0.
a 1 = t x ρ u t y ρ v t x ρ u t y ρ v a 2 = t x ρ u t y v t x v t y ρ u + ρ v ( t x ρ u t y ρ t x ρ t y ρ u ) a 3 = t x ρ v t y v t x v t y ρ v + t x u t y ρ u t x ρ u t y u + ρ u ( t x ρ t y ρ u t x ρ u t y ρ ) + ρ v ( t x ρ v t y ρ t x ρ t y ρ v ) a 4 = t x u t y ρ v t x ρ v t y u + ρ u ( t x ρ t y ρ v t x ρ v t y ρ ) a 5 = t x u t y v t x v t y u + ρ u ( t x ρ t y v t x v t y ρ ) + ρ u ( t x u t y ρ t x ρ t y u ) .
{ a 1 ( ρ u u ρ v v ρ u v 2 ) + a 2 ρ u u + a 3 ρ u v + a 4 ρ v v + a 5 = 0 BC : { t x = t x ( u , v , ρ , ρ u , ρ v ) t y = t y ( u , v , ρ , ρ u , ρ v ) : Ω 0 Ω 1 .
ρ u = ρ i + 1 , j ρ i 1 , j 2 h i , ρ v = ρ i , j + 1 ρ i , j 1 2 h j , ρ u u = ρ i + 1 , j ρ z i , j + ρ i 1 , j h i 2 ρ v v = ρ i , j + 1 2 ρ i , j + ρ i , j 1 h j 2 , ρ u v = ρ i + 1 , j + 1 ρ i + 1 , j 1 ρ i 1 , j + 1 + ρ i 1 , j 1 4 h i h j .
ρ u = 3 ρ m , j 4 ρ m 1 , j + ρ m 2 , j 2 h i , ρ v = ρ m , j + 1 ρ m , j 1 2 h j .
u = E max E min E max + E min
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.