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

Construction method for designing a spectrometer with variable spectral resolution and wide bandwidth using multiple off-axis convex gratings

Open Access Open Access

Abstract

The classic Offner spectrometer has the advantages of having a simple structure, light weight, and high imaging quality. It is easily to achieve a fixed spectral resolution but cannot meet usage requirements. Therefore, we present a practical method for designing a spectrometer with variable spectral resolution. Multiple off-axis convex (OAC) gratings are used to replace the convex grating in the classic Offner spectrometer. We derive the principle through ray tracing and establish an optimization process for the basic parameters of multiple OAC gratings. To demonstrate this method, a corresponding system is designed. The results show that a variable spectral resolution, with a variation ratio close to 4, of 0.45–1.91 nm is achieved over a wide bandwidth of 460–900 nm. Additionally, the smile and keystone of the system are well corrected.

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

1. Introduction

In recent years, spectral imaging technology has been widely used in various fields such as food safety [1,2], environmental monitoring [3,4], archaeological and art preservation [5,6], agriculture [7], and biomedicine [8]. While spectrometers are used to identify and classify a wide range of complex objects, there are higher demands on the speed of spectral data acquisition and processing. In addition, spectrometers with a large field of view and high resolution are desired. However, such spectrometer systems will generate huge image data cubes, which will reduce the efficiency of information processing and transmission.

In practical applications, in the observation area of the full field, different targets require different spectral resolutions to be identified. If the spectrum data cube of all objects in the full field is acquired only with the highest resolution, the unnecessary data volume and transmission pressure are increased. If the spectrum data cube of all objects in the full field is acquired only at the lowest resolution, targets that require hyperspectral cannot be recognized. More specifically, we need a spectrometer that can simultaneously obtain high spectral resolution imaging of the target of interest and low spectral resolution imaging of the full field. When we only need the monitoring information of the full field, the low-resolution output data can meet the requirements. If we want to accurately analyze the spectral characteristics of the interested target, the high-resolution output data can meet the requirements. Therefore, a spectrometer with variable spectral resolution is urgently required to realize the above conditions.

Currently, spectrometers with variable spectral resolution can be broadly classified into three categories: First, spatial light modulation devices such as DMD [9,10] are used to modulate light at the slit to achieve changes in spatial and spectral resolutions. This approach can only obtain the spectral information of one modulation mode per imaging scan. Specifically, spectral information at different spectral resolutions cannot be obtained using a single scan. In the second category, different spectral resolutions are obtained using multiple channel imaging systems, which require multiple detectors [11] or multiple slits, gratings, and reflectors stitched together [12]. This approach is structurally complex and does not meet the compactness requirements. In the third category, a spectrometer with variable spectral resolution is achieved by designing and machining varied-line-spacing gratings [13,14], which involves a complex and costly machining process.

Classic Offner spectrometers have a high imaging quality, simple structure, and smaller volume and weight than those of plane grating imaging spectrometers [15]. In this study, we took advantage of the simplicity and compactness of the traditional Offner spectrometer and made improvements accordingly. Herein, we propose a design concept for a spectrometer with variable spectral resolution based on multiple off-axis convex (OAC) gratings, which are used to replace the convex spherical grating in the classical Offner structure.

The remainder of this paper is organized as follows: Section 2 presents the structural principle of the spectrometer with variable spectral resolution and multiple OAC gratings. According to this principle, each beam of light diffracted from an OAC grating generates a sub-spectral image at the image plane. However, each sub-spectral image overlaps, so a method is required to avoid this. Section 3 presents the automatic optimization process established to prevent overlapping of sub-spectral images, and calculations of the initial parameters of a well example. Section 4 describes further optimization of the design example by considering the distortion and increasing the optimizable parameters of the system. The spectral resolution and distortion are also analyzed. As a result, the spectral resolution variation ratio is found to be approximately 4. To prove the practicability of the system, we analyzed the SNR, as shown in Section 5. Finally, conclusions are presented in Section 6.

This novel structure provides different spectral resolution imaging through a single scan and has the advantage of being both simple and compact. This method does not cause a large number of problems while obtaining hyperspectral data. In addition, compared with the flexibility of digital spectrometers, the proposed spectrometer does not require additional digital components, and has strong environmental adaptability and compact structure, and can be applied to the infrared band. Thus, it can be used for agricultural monitoring, crop identification, and classification.

2. Structural principle

We substituted multiple OAC gratings for each subregion of the convex grating in the classic Offner spectrometer. The parameters of each OAC grating include the curvature, grating density, decenter, diameter, angle of incidence, and diffraction angle. The algorithms used to solve these parameters are described in this section. Moreover, the spectrometer is uniformly defined in the global coordinate system through coordinate transformations to improve the convenience of construction.

2.1 Parameter definitions for the spectrometer

A classic Offner spectrometer in the YZ plane is displayed in Fig. 1. The beam emitted from the slit is reflected by the primary mirror and then diffracted at the convex grating. The tertiary mirror focuses the dispersive beam onto the image plane. The primary and tertiary mirrors are on the same spherical surface. Figure 2 presents a diagram of the ray tracing of the Offner structure in the XYZ three-dimensional space. The origin of the coordinate system O is at the center of the surface where the primary and tertiary mirrors are located. The coordinate system is a Cartesian three-dimensional coordinate system. The line, i.e., the Z-axis, between the center of the grating surface and the origin of the coordinate system is the optical axis of the Offner spectrometer, which is positive from the grating toward the origin. The XOY plane is perpendicular to the Z-axis.

 figure: Fig. 1.

Fig. 1. Classic Offner spectrometer

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. Schematic of light tracing *Coordinate system XYZ is the global coordinate and O is the origin of the global coordinate system. Coordinate system X'Y'Z’ represented by the dotted line at the slit is the local coordinates. The center of the slit is the origin of the local coordinate system, the direction of which is the same as that of the global coordinate system. (0, y0, z0) are the coordinates of the slit center in the global coordinate system.

Download Full Size | PDF

At the slit, the X-axis indicates the slit length and the Y-axis indicates the slit width. As shown in Fig. 2, one of the rays originating from the slit is traced. To facilitate the corresponding calculation, we take the center point at the slit x0 = 0 as the starting point of the ray tracing. (0, y0, z0), (xp, yp, zp), (xg, yg, zg), $({{x_{\textrm{t},{\lambda_n}}},{y_{\textrm{t},{\lambda_n}}},{z_{\textrm{t},{\lambda_n}}}} )$, and $({{x_{\textrm{i},{\lambda_n}}},{y_{\textrm{i},{\lambda_n}}},{z_{\textrm{i},{\lambda_n}}}} )$ denote the coordinates of the ray at the center of the slit, the primary mirror, the grating, the tertiary mirror, and the image plane, respectively. $\overrightarrow {{l_{\textrm{sp}}}} $, $\overrightarrow {{l_{\textrm{pg}}}} $, $\overrightarrow {{l_{\textrm{gt},{\lambda _n}}}} $, and $\overrightarrow {{l_{\textrm{ti},{\lambda _n}}}} $ denote the light vector from the slit to the primary mirror, the primary mirror to convex grating, the convex grating to the tertiary mirror, and the tertiary mirror to the image plane, respectively. In the image plane, the Y-coordinate represents the coordinates of the spectral information at a certain spatial position, and the X-coordinate represents the coordinates of the spectral information at different adjacent spatial positions. , $\overrightarrow {{n_g}} $, and $\overrightarrow {{n_{t,{\lambda _n}}}} $ denote the unit normal vector of the primary mirror in (xp, yp, zp), convex grating in (xg, yg, zg), and tertiary mirror in $({{x_{\textrm{t},{\lambda_n}}},{y_{\textrm{t},{\lambda_n}}},{z_{\textrm{t},{\lambda_n}}}} )$, respectively.

θ is the angle between $\overrightarrow {{l_{\textrm{sp}}}} $ and the Z-axis, and φ is the angle between the projection of $\overrightarrow {{l_{\textrm{sp}}}} $ in the XY plane and the X-axis. Based on the relationship of the space geometry, θ and φ can be expressed by some parameters of the global coordinate system, as follows:

$$\tan \varphi = \frac{{{y_p} - {y_0}}}{{{x_p}}}, $$
$$\tan \theta = \frac{{\sqrt {{{({{y_p} - {y_0}} )}^2} + x_p^2} }}{{{z_p} - {z_0}}}. $$

The primary mirror is spherical with curvature cp. The substrate of the convex grating is a spherical surface with curvature cg and D is the distance from the primary mirror to the convex grating. Thus, zp and zg are given by:

$${z_p} = \frac{{{c_p}(x_p^2 + y_p^2)}}{{1 + \sqrt {1 - c_p^2(x_p^2 + y_p^2)} }}, $$
$${z_g} = \frac{{{c_g}(x_g^2 + y_g^2)}}{{1 + \sqrt {1 - c_g^2(x_g^2 + y_g^2)} }}\textrm{ - }D. $$

Equations (5) and (6) are obtained according to the law of reflection.

$$\left( {\frac{{\overrightarrow {{l_{pg}}} }}{{|{\overrightarrow {{l_{pg}}} } |}} - \frac{{\overrightarrow {{l_{sp}}} }}{{|{\overrightarrow {{l_{sp}}} } |}}} \right) \bullet \overrightarrow {{n_p}} = 0$$
$$\overrightarrow {{l_{pg}}} \cdot (\overrightarrow {{l_{sp}}} \times \overrightarrow {{n_p}} ) = 0$$
$\overrightarrow {{l_{\textrm{pg}}}} $, $\overrightarrow {{l_{\textrm{sp}}}} $, and $\overrightarrow {{n_\textrm{p}}} $ can be expressed using y0, z0, xp, yp, zp, xg, yg, zg, and cp as:
$$\overrightarrow {{n_p}} = \left( {{c_p}{x_p},{c_p}{y_p}, - \sqrt {1 - c_p^2(x_p^2 + y_p^2)} } \right), $$
$$\overrightarrow {{l_{sp}}} = ({{x_p} - {x_0},{y_p} - {y_0},{z_p} - {z_0}} ), $$
$$\overrightarrow {{l_{pg}}} = ({{x_g} - {x_p},{y_g} - {y_p},{z_g} - {z_p}} ). $$

Similarly, $\overrightarrow {{n_\textrm{g}}} $ is given by:

$$\overrightarrow {{n_g}} = \left( {{c_g}{x_g},{c_g}{y_g}, - \sqrt {1 - c_g^2(x_g^2 + y_g^2)} } \right). $$

According to Eqs. (1)–(10), the angle of incidence θi on the convex grating can be obtained from Eq. (11).

$${\theta _i} = \arccos \left( {\frac{{\overrightarrow {{l_{pg}}} \bullet \overrightarrow {{n_g}} }}{{|{\overrightarrow {{l_{pg}}} } |}}} \right)$$

For a convex grating, ${\theta _{e,{\lambda _n}}}$ is the diffraction angle. The relationship between θi and ${\theta _{e,{\lambda _n}}}$ is determined by the grating density T in lines per millimeter (g/mm), given by:

$$\sin {\theta _i}\textrm{ - }\sin {\theta _{e,{\lambda _n}}} = M{\lambda _n}T, $$
where M is the diffraction order, equal to 1(the diffraction order other than the utilized first positive from one grating will not affect the neighboring configuration), and λn is the wavelength of the dispersed beam. Based on the relationship of space geometry, ${\theta _{e,{\lambda _n}}}$ can also be expressed as:
$${\theta _{e,{\lambda _n}}} = \arccos \left( {\frac{{\overrightarrow {{l_{gt}}} \bullet{-} \overrightarrow {{n_g}} }}{{|{\overrightarrow {{l_{gt}}} } |}}} \right). $$

When a single ray is reflected by the grating, the relationship between the incident light, reflected light, and the unit normal vector of the grating at the incident point can be expressed as:

$$\overrightarrow {{l_{gt}}} \cdot (\overrightarrow {{l_{pg}}} \times \overrightarrow {{n_g}} ) = 0. $$

The tertiary mirror is spherical with curvature ct. Then, ${z_{\textrm{t},{\lambda _n}}}$ is given by:

$${z_{t,{\lambda _n}}} = \frac{{{c_t}(x_{t,{\lambda _n}}^2 + y_{t,{\lambda _n}}^2)}}{{1 + \sqrt {1 - c_t^2(x_{t,{\lambda _n}}^2 + y_{t,{\lambda _n}}^2)} }}. $$

Equations (16 and (17) were obtained according to the law of reflection.

$$\overrightarrow {{l_{ti}}} \cdot (\overrightarrow {{l_{gt}}} \times \overrightarrow {{n_{t,{\lambda _n}}}} ) = 0$$
$$\left( {\frac{{ - \overrightarrow {{l_{gt}}} }}{{|{ - \overrightarrow {{l_{gt}}} } |}} - \frac{{\overrightarrow {{l_{ti}}} }}{{|{\overrightarrow {{l_{ti}}} } |}}} \right) \bullet \overrightarrow {{n_{t,{\lambda _n}}}} = 0$$
$\overrightarrow {{l_{\textrm{gt}}}} $, $\overrightarrow {{l_{\textrm{ti}}}} $, and $\overrightarrow {{n_{\textrm{t},{\lambda _n}}}} $ can be expressed using xg, yg, zg, ${x_{\textrm{t},{\lambda _n}}}$, ${y_{\textrm{t},{\lambda _n}}}$, ${z_{\textrm{t},{\lambda _n}}}$, ${x_{\textrm{i},{\lambda _n}}}$, ${y_{\textrm{i},{\lambda _n}}}$, ${z_{\textrm{i},{\lambda _n}}}$, and ct as Eqs. (18)–(20).
$$\overrightarrow {{l_{gt}}} \textrm{ = }({{x_{t,{\lambda_n}}} - {x_g},{y_{t,{\lambda_n}}} - {y_g},{z_{t,{\lambda_n}}} - {z_g}} )$$
$$\overrightarrow {{l_{ti}}} = ({{x_{i,{\lambda_n}}} - {x_{t,{\lambda_n}}},{y_{i,{\lambda_n}}} - {y_{t,{\lambda_n}}},{z_{i,{\lambda_n}}} - {z_{t,{\lambda_n}}}} )$$
$$\overrightarrow {{n_{t,{\lambda _n}}}} = \left( {{c_t}{x_{t,{\lambda_n}}},{c_t}{y_{t,{\lambda_n}}}, - \sqrt {1 - c_t^2(x_{t,{\lambda_n}}^2 + y_{t,{\lambda_n}}^2)} } \right)$$

According to Eq. (12), when T is determined, it follows that ${\theta _{\textrm{e},{\lambda _n}}}$. Further, according to Eqs. (13)–(20), $({{x_{\textrm{i},{\lambda_n}}},{y_{\textrm{i},{\lambda_n}}},{z_{\textrm{i},{\lambda_n}}}} )$ can be obtained. More specifically, using the model comprising Eqs. (1)–(20), the initial structure of interest can be finalized once the relevant limitation factor is input into the model.

2.2 Substitution of multiple OAC gratings

There are m OAC gratings: g1, …, gm−1, gm. Figure 3(a) shows a diagram of the OAC gratings in the YZ plane. An OAC grating is simply a side section of a parent convex grating with a decenter d (Y offset in the tangential plane). g1, …, gm−1, gm were fitted to a sphere of curvature cg. The diameter of the OAC gratings is ϕy. The Y-axis coordinates of the edge ray and the center ray through the multiple OAC gratings are yp,max, yp,min, and yp,O. Figure 3(b) presents a diagram of the multiple OAC gratings in the XY plane. Each OAC grating has a different grating density, labeled T1, …, Tm−1, Tm. Therefore, m OAC gratings diffract m sub-spectral images with different dispersion widths.

 figure: Fig. 3.

Fig. 3. Diagram of multiple OAC gratings

Download Full Size | PDF

Figure 4 shows the replacement of the convex grating in the classic Offner structure with the proposed multiple OAC gratings. The beam is diffracted at multiple OAC gratings with different grating densities Tm. According to Eq. (12), the diffraction angle $\theta_{e,{\lambda_n},{g_m}}$ of each OAC grating is influenced by the wavelength $\lambda_{n}$ and Tm. Moreover, according to Eqs. (13)–(20), the coordinates of the beam on the image plane also depend on λn and Tm. Figure 5 shows the image plane in the XY plane. The beam in the λ1λn bandwidth (central wavelength λc) is diffracted at OAC grating gm and focused on the image plane by the tertiary mirror; its dispersion width Δpm is expressed as:

$$\varDelta {p_m} = |{{y_{i,{\lambda_1},{g_m}}} - {y_{i,{\lambda_n},{g_m}}}} |, $$
where ${y_{\textrm{i},{\lambda _n},{g_m}}}$ and ${y_{\textrm{i},{\lambda _1},{g_m}}}$ are the Y-coordinates of the spectral lines with λn and λ1 generated by gm in the image plane. Similarly, the dispersion width Δpm−1 generated by gm−1 is expressed as:
$$\varDelta {p_{m - 1}} = |{{y_{i,{\lambda_1},{g_{m - 1}}}} - {y_{i,{\lambda_n},{g_{m - 1}}}}} |, $$
where ${y_{\textrm{i},{\lambda _n},{g_{m - 1}}}}$ and ${y_{\textrm{i},{\lambda _1},{g_{m - 1}}}}$ are the Y-coordinates of the spectral lines with λn and λ1 generated by gm−1 in the image plane. The offset Δhm,m−1 between the sub-spectral images generated by gm and gm−1 in the image plane is given by:
$$\Delta {h_{m,m - 1}} = |{{y_{i,{\lambda_c},{g_{m - 1}}}} - {y_{i,{\lambda_c},{g_m}}}} |, $$
where ${y_{\textrm{i},{\lambda _c},{g_m}}}$ and ${y_{\textrm{i},{\lambda _c},{g_{m - 1}}}}$ are the Y-coordinates of the central spectral line λc generated by gm and gm−1 in the image plane, respectively. To ensure that no overlaps occur between the sub-spectral images, Δpm, Δpm−1, and Δhm,m−1 should satisfy the following equation:
$$2\Delta {h_{m,m - 1}} \ge \varDelta {p_{m - 1}}\textrm{ + }\varDelta {p_m}. $$

 figure: Fig. 4.

Fig. 4. Schematic showing ray tracing of multiple OAC gratings

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. XY plane of spectral imaging

Download Full Size | PDF

Δpm, Δpm−1, and Δhm,m−1 can be expressed by Y-coordinates on the image plane, which depend on λn and Tm. Therefore, the occurrence of sub-spectral image overlap is mainly determined by the wavelength λ1λn and the grating density T1, …, Tm−1, Tm of the OAC gratings.

2.3 Distortion definition

Owing to the dispersive behavior, two types of distortion require special attention when evaluating the spectrometer performance. These are spectral distortions (smiles), which are bent monochromatic line images along the X-axis, and spatial distortion (keystones), which refer to wavelength-dependent bent slit images along the Y-axis [16,17]. In this study, we used Eqs. (25) and (26) to define these two distortions [18]. Smile distortion is defined in Eq. (25), where ${y_{\textrm{max},{\lambda _n},{g_m}}}$ and ${y_{\textrm{min},{\lambda _n},{g_m}}}$ represent the maximum and minimum Y-coordinates of the bent monochromatic line in the sub-spectral image, respectively. Keystone distortion is defined in Eq. (26), where ${x_{\textrm{max},{\lambda _n},{g_m}}}$ and ${x_{\textrm{min},{\lambda _n},{g_m}}}$ represent the maximum and minimum X-coordinates in the same slit imaging in the sub-spectral image, respectively, both of which are wavelength dependent.

$${D_{smile,{\lambda _n},{g_m}}} = {y_{\max ,{\lambda _n},}}_{{g_m}} - {y_{\min ,{\lambda _n},}}_{{g_m}}$$
$${D_{keystone,{\lambda _n},{g_m}}} = {x_{\max ,{\lambda _n},{g_m}}} - {x_{\min ,{\lambda _n},{g_m}}}$$

2.4 Definition of spectral resolution

The spectral resolution is determined by the detector pixel size aD, RMS spot radius aR,λ, dispersion width Δpm, and bandwidth λ1λn. When aR,λ is smaller than aD, the spectral resolution is limited by aD and is expressed as Eq. (27). The variation ratio βD of the spectral resolution of the spectrometer with multiple OAC gratings is given by Eq. (28). When aD and λ1λn of the system are determined, βD can be obtained by the maximum dispersion widths Δpmax and minimum dispersion widths Δpmin of the image plane. When aR,λ is larger than aD, the spectral resolution is limited by aR,λ, as expressed in Eq. (29), and βD can be realized using Eq. (30). The RMS spot radius varies with wavelength, so the spectral resolution variation ratio is determined by not only Δpmax and Δpmin, but also aR,λ at different wavelengths. In Section 3, for facilitation, Eq. (28) is used as the basis for the initial calculation in the optimization process of the spectrometer with multiple OAC gratings.

$$\delta {\lambda _{D,m}} = \frac{{{a_D}}}{{\Delta {p_m}}}|{{\lambda_1} - {\lambda_n}} |$$
$${\beta _D}\textrm{ = }\frac{{\delta {\lambda _{D,\min }}}}{{\delta {\lambda _{D,\max }}}} = \frac{{\Delta {p_{\min }}}}{{\Delta {p_{\max }}}}$$
$$\delta {\lambda _{R,m}} = \frac{{{a_{R,\lambda }}}}{{\Delta {p_m}}}|{{\lambda_1} - {\lambda_n}} |$$
$${\beta _D}\textrm{ = }\frac{{\delta {\lambda _{R,\min }}}}{{\delta {\lambda _{R,\max }}}}$$

3. Spectrometer construction

3.1 Design specification

The following example details the methodology of a spectrometer with variable spectral resolution. A classic Offner spectrometer is designed as the initial system for the OAC grating spectrometer. A wide bandwidth of 460–900 nm was considered. The spectrum in this range was divided into thirteen sampling wavelengths, labeled as λ1λ13, and the central wavelength is 680 nm, labeled as λ7. The numerical aperture for the divergent beam from the entrance slit was set to 0.145, which matched the optical system with an F value of 3.4. A detector with 2,048 px × 2,048 px was chosen, with the pixel size of 6.5 µm × 6.5 µm. A convex grating with the grating density T of 150 g/mm, and the first positive diffraction order were utilized. The basic parameters of the initial system are listed in Table 1.

Tables Icon

Table 1. Basic parameters of initial system

3.2 Optimization process of the spectrometer with multiple OAC gratings

We replaced the convex grating with multiple OAC gratings based on the initial system. In this example, the convex grating was split into three OAC gratings with grating densities T1, T2, and T3, and decenters ${d_{{g_1}}}$, ${d_{{g_2}}}$, and ${d_{{g_3}}}$. The optimization process is depicted in Fig. 6.

 figure: Fig. 6.

Fig. 6. Flowchart of construction stage of the spectrometer with multiple OAC gratings

Download Full Size | PDF

First, the parameters of the initial system were fed into the optimization process. For easy manufacturing, the grating density should be greater than 40 lp/mm and less than 1,000 lp/mm. The throughput aperture of each OAC grating should be approximately the same to ensure that the energy of each sub-spectral image is balanced. Therefore, we define the decenters of the OAC grating to satisfy ${d_{{g_2}}} = 0$ and $|{{d_{{g_1}}}} |= |{{d_{{g_3}}}} |= \frac{1}{3}{\phi _y}$. The Y-coordinates ${y_{{\textrm{p}_\textrm{c}},{g_1}}}$, ${y_{{\textrm{p}_\textrm{c}},{g_2}}}$, and ${y_{{\textrm{p}_\textrm{c}},{g_3}}}$ of the central rays of the three OAC gratings can be expressed as Eqs. (31)–(33).

$${y_{{p_c},{g_1}}} = \frac{1}{3}{\phi _y} + \frac{1}{2}({y_{p\max }} + {y_{p\min }})$$
$${y_{{p_c},{g_2}}} = 0$$
$${y_{{p_c},{g_3}}} = \textrm{ - }\frac{1}{3}{\phi _y} + \frac{1}{2}({y_{p\max }} + {y_{p\min }})$$

Here, ${\phi _y} = |{{y_{\textrm{pmax}}} - {y_{\textrm{pmin}}}} |$. According to Eqs. (1)–(11), the incidence angles of the central ray in the OAC grating can be obtained. Second, the bandwidth and central wavelength are input into the optimization process. According to Eqs. (21)–(23), the dispersion widths Δp1, Δp2, and Δp3, and offsets Δh12 and Δh23 of the three sub-spectral images can be obtained. Then, the optimization process determines if there are overlaps between the sub-spectral images according to Eq. (24). In addition, according to the effective detector size L, the sum of the three dispersion widths must satisfy Eq. (34). If the optimization result does not satisfy Eqs. (24) and (34), the process returns to reset the grating density. If it is satisfied, the process outputs results for T1, T2, T3, ${d_{{g_1}}}$, ${d_{{g_2}}}$, ${d_{{g_3}}}$, βD, and ${y_{\textrm{i},{\lambda _n},{g_1}}} - {y_{\textrm{i},{\lambda _1},{g_3}}}$. The spectral resolution variation ratio βD can be derived using Eq. (28). ${y_{\textrm{i},{\lambda _n},{g_1}}} - {y_{\textrm{i},{\lambda _1},{g_3}}}$ denotes the Y-axis space occupied by all the dispersive spectral lines in the image plane.

$$\varDelta {p_1} + \varDelta {p_2}\textrm{ + }\varDelta {p_3} \le L$$

Through the above optimization process, the five sets of data shown in Table 2 satisfying Eqs. (24) and (34) are output. In the first and fourth sets of data, βD is higher than 4. However, the ${y_{\textrm{i},{\lambda _n},{g_1}}} - {y_{\textrm{i},{\lambda _1},{g_3}}}$ value of 8.45 mm in the first dataset is much smaller than the effective detector size L of 13 mm, which reduces the detector utilization. Meanwhile, the fifth dataset shows a higher detector utilization but a lower spectral resolution variation ratio. In this study, we chose the fourth set of data for the initial parameters of the multiple OAC gratings.

Tables Icon

Table 2. Five sets of optimization results

4. Performance analysis

Further optimization was performed using the constructed spectrometer. To increase the optimization parameters of the system, the primary and tertiary mirrors of the initial system were separated so that they are not in one mirror and have different curvatures. Meanwhile, the smile and keystone are controlled according to Eqs. (25) and (26). The results are presented in Fig. 7. The colors of the rays in Fig. 7(a) indicate different wavelengths. The colors of the light in Fig. 7(b) indicate the dispersive spectral lines diffracted from each OAC grating. The parameters of the spectrometer with multiple OAC gratings are listed in Table 3, and those of the multiple OAC gratings are provided in Table 4. The grating comprised three OAC gratings with grating densities of 300, 150, and 70 g/mm, respectively. The decenters of the three OAC gratings were adjusted slightly to improve the detector utilization, to 5.5, 0.45, and −4.5 mm, respectively. The grating aperture is 15.4 mm.

 figure: Fig. 7.

Fig. 7. Spectrometer with multiple OAC gratings

Download Full Size | PDF

Tables Icon

Table 3. Parameters of spectrometer with multiple OAC gratings

Tables Icon

Table 4. Parameters of multiple OAC gratings

For our optimization, we reduced the distortion by optimizing the curvature and distances of the optical elements. The spectral smiles and keystones are shown in Fig. 8. The maximum smile, produced by OAC grating 1 at the end of the slit at the wavelength of 900 nm, is less than 3.5 µm, which is approximately half of the pixel size. In addition, the maximum keystone produced by the OAC grating 1 at the end of the slit at 900 nm is less than 3.2 µm, which is also approximately half of the pixel size.

 figure: Fig. 8.

Fig. 8. Distortion of spectrometer with multiple OAC gratings. (a)–(c) Spectral smiles for different wavelengths; (d)–(f) Spectral keystones for 0, 0.2, 0.4, 0.6, 0.8, and 1 normalized fields of view

Download Full Size | PDF

Figure 9 displays the full-field spot diagram in the image plane. The X-axis indicates the imaging of the slit, and the coordinates are represented by the normalized field of view. The Y-axis indicates the dispersion spectral lines of the system, and the colors represent different dispersion wavelengths. Figure 9 shows that the system generates three sub-spectral-images with dispersion widths of Δp1, Δp2, and Δp3 in the image plane.

 figure: Fig. 9.

Fig. 9. Full-field spot diagram of the imaging

Download Full Size | PDF

Figure 10 presents a comparison of the RMS spot radius for the three sub-spectral images. The RMS spot radius diffracted by OAC grating 3 from 460 to 900 nm is smaller than the detector pixel size of 6.5 µm. The RMS spot radius diffracted by OAC gratings 1 and 2 are smaller than the detector pixel size of 6.5 µm from 695 to 900 nm and from 460 to 812 nm, respectively.

 figure: Fig. 10.

Fig. 10. RMS spot radius diffracted by three OAC gratings

Download Full Size | PDF

The spectral resolution of the system, according to the definition provided in Section 2.4, is shown in Fig. 11(a). When the wavelength is 460–636 nm and 900 nm, the spectral resolution diffracted by OAC grating 1 is between 0.46 and 0.49 nm, and is limited by the RMS spot radius. When the wavelength lies between 665 and 871 nm, the spectral resolution, to which the detector pixel size is the main limitation factor, has the value of 0.45 nm. Moreover, when diffracted by OAC grating 2, the spectral resolution varies between 0.90 and 0.93 nm with the wavelength of 841–900 nm. The spectral resolution is 0.89 nm with wavelength of 460–812 nm. For these two spectral resolutions, the former is limited by the RMS spot radius, while the latter is limited by the detector pixel size. When OAC grating 3 is working, the spectral resolution is 1.91 nm at the full wavelength range of 460–900 nm and is limited by the detector pixel size.

 figure: Fig. 11.

Fig. 11. (a) Spectral resolution diffracted by three OAC gratings; (b) Spectral resolution variation ratio

Download Full Size | PDF

According to Eqs. (28) and (30), the spectral resolution variation ratio of the system was calculated, as shown in Fig. 11(b). β is in the range of 3.86–3.94 at the band of 489–607 nm, and greater than the spectral resolution variation ratio of 4.02 for the initial parameters presented in Section 3.2, at the rest of the bands. The highest β value is 4.27.

5 Prediction of the signal-to-noise ratio

The two main factors affecting the SNR of each sub-image with different spectral resolutions are the diffraction efficiency of the OAC grating and the relative aperture of the optical system. According to the processing level of the convex blazed gratings [19,20], the average diffraction efficiencies of the three OAC gratings can reach 60%, 50%, and 40%. Because the convex grating is divided into three OAC gratings, the F value corresponding to each OAC grating is approximately 10. To reduce the etalon effect in the hyperspectral system, we used a detector from HORIBA, model 4.2MP-U-6.5-BI. In the same integration of 0.3 s, the SNR shown in Fig. 12, which was determined while considering the atmospheric effects, is not less than 80 for the three different spectral resolution images. The SNRs of the sub-images generated by the three OAC gratings all meet the detection requirements. Therefore, the proposed system can obtain images with different spectral resolutions at the same integration time. Moreover, the SNR will be higher in the case of overlaying the frame frequency or spatial merging. In addition, we can choose a larger detector and design an optical system with a smaller F number according to the design method in this study. The detailed process of SNR analysis is provided in Supplement 1.

 figure: Fig. 12.

Fig. 12. SNR of each sub-imaging produced by the three OAC gratings

Download Full Size | PDF

6 Conclusions

We introduced a construction method for designing a spectrometer with variable spectral resolution and multiple OAC gratings. We derived the relationship between the image plane coordinates, the basic parameters of the multiple OAC gratings, and the system wavebands via ray tracing. Based on the theory presented in Section 2, we established an optimization process to calculate the optimal initial parameters of the multiple OAC gratings and designed an example system. We also controlled the system smile and keystone. The results demonstrate that the system distortion was well corrected. Finally, we obtained a spectral resolution variation ratio close to 4. Specifically, this system can obtain spectral imaging with the spectral resolution variable from 0.45 to 1.91 nm. Users can choose spectral imaging with different resolutions according to their own needs.

In the example provided in Section 3, the grating comprised three OAC gratings. In future applications, the system can obtain more different resolutions of spectral imaging using more than three OAC gratings. To ensure a high SNR, larger detector sizes will be used when the grating consists of more OAC gratings. Alternatively, if the basic parameters in the example are used, the optimization process can output a higher spectral resolution variation ratio when the detector size is larger. In this paper, we provide an example of a 4-fold spectral resolution variation ratio, but our proposed construction method also enables design of a spectrometer with variable spectral resolution having different specifications according to the requirements. Thus, the proposed system can be widely used in various fields.

Disclosures

The authors declare that there are no conflicts of interest related to this article.

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.

Supplemental document

See Supplement 1 for supporting content.

References

1. Z. Xiong, D. Sun, X. Zeng, and A. Xie, “Recent developments of hyperspectral imaging systems and their applications in detecting quality attributes of red meats: A review,” J. J. Eng. 132(1), 1–13 (2014). [CrossRef]  

2. A. Gowen, C. P. O’Donnell, P. J. Cullen, G. Downey, and J. M. Frias, “Hyperspectral imaging – An emerging process analytical tool for food quality and safety control,” Trends Food Sci. Technol. 18(12), 590–598 (2007). [CrossRef]  

3. E Honkavaara, T Hakala, L Markelin, A Jaakkola, H Saari, H Ojanen, I Pölönen, S Tuominen, R Näsi, T Rosnell, and N Viljanen, “Autonomous hyperspectral UAS photogrammetry for environmental monitoring applications,” J. ISPRS - International Archives of the Photogrammetry. XL-1:155–159 (2014).

4. H. Zhu, H. Liu, Y. Xu, and Y. Guijun, “UAV-based hyperspectral analysis and spectral indices constructing for quantitatively monitoring leaf nitrogen content of winter wheat,” Appl. Opt. 57(27), 7722–7732 (2018). [CrossRef]  

5. P. Xu, H. Xu, C. Diao, and Z. Ye, “Self-training-based spectral image reconstruction for art paintings with multispectral imaging,” Appl. Opt. 56(30), 8461 (2017). [CrossRef]  

6. H. Liang, “Advances in multispectral and hyperspectral imaging for archaeology and art conservation,” Appl. Phys. A. 106(2), 309–323 (2012). [CrossRef]  

7. L. M. Dale, A. Thewis, C. Boudry, I. Rotar, P. Dardenne, V. Baeten, and J. A. F. Pierna, “Hyperspectral imaging applications in agriculture and agro-food product quality and safety control: A review,” Appl. Spectrosc. Rev. 48(2), 142–159 (2013). [CrossRef]  

8. P. Y. Mazzer, C. H. Barbieri, N. Mazzer, and V. P. Fazan, “Morphologic and morphometric evaluation of experimental acute crush injuries of the sciatic nerve of rats,” J. Neurosci. Methods. 173(2), 249–258 (2008). [CrossRef]  

9. A Jung, R Michels, and R Graser, “Hyperspectral camera with spatial and spectral resolution and method,” P. EP20140001857 (2015).

10. X. Dong, X. Xiao, Y. Pan, G. Wang, and Y. Yu, “DMD-based hyperspectral imaging system with tunable spatial and spectral resolution,” Opt. Express. 27(12), 16995–17006 (2019). [CrossRef]  

11. L. Feng, J. Zhou, L. Wei, X. He, Y. Li, J. Jing, and B. Xiangli, “Design of a compact wide-spectrum double-channel prism imaging spectrometer with freeform surface,” Appl. Opt. 57(31), 9512–9522 (2018). [CrossRef]  

12. Q. Xue, F. Lu, M. Duan, Y. Zheng, X. Wang, D. Cao, G. Lin, and J. Tian, “Optical design of double-grating and double wave band spectrometers using a common CCD,” Appl. Opt. 57(23), 6823–6830 (2018). [CrossRef]  

13. W. R. Palmer and W. R. McKinney, “Imaging theory of plane-symmetric varied line-space grating systems,” Opt. Eng. 33(3), 820–829 (1994). [CrossRef]  

14. J. Chen, W. Chen, G. Zhang, H. Lin, and S. Chen, “High-resolution compact spectrometer based on a custom-printed varied-line-spacing concave blazed grating,” Opt. Express. 25(11), 12446–12454 (2017). [CrossRef]  

15. D. R. Lobb, “Theory of concentric designs for grating spectrometers,” Appl. Opt. 33(13), 2648–2658 (1994). [CrossRef]  

16. L. Yuan, J. Xie, Z. He, Y. Wang, and J. Wang, “Optical design and evaluation of airborne prism-grating imaging spectrometer,” Opt. Express. 27(13), 17686–17700 (2019). [CrossRef]  

17. T. Skauli, “Quantifying coregistration errors in Spectral Imaging,” Proc. SPIE. 8158, 81580A (2011). [CrossRef]  

18. J. Hong, Y. Kim, B. Choi, S. Hwang, D. Jeong, J. H. Lee, Y. Kim, and H. Kim, “Efficient method to measure the spectral distortions using periodically distributed slit in hyperspectral imager,” Opt. Express. 25(17), 20340–20351 (2017). [CrossRef]  

19. J Fisher, J A Antoniades, C Rollins, and X LianQin, “Hyperspectral imaging sensor for the coastal environment,” Proc. SPIE 3482, International Optical Design Conference, 1998.

20. F Zamkotsian, I Zhurminsky, P Lanzoni, T Nicolas, and C Vincent, “Convex blazed gratings for high throughput spectrographs in space missions,”. International Conference on Space Optics - ICSO 2018, 2019.

Supplementary Material (1)

NameDescription
Supplement 1       SNR Analysis Detailed Process

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 (12)

Fig. 1.
Fig. 1. Classic Offner spectrometer
Fig. 2.
Fig. 2. Schematic of light tracing *Coordinate system XYZ is the global coordinate and O is the origin of the global coordinate system. Coordinate system X'Y'Z’ represented by the dotted line at the slit is the local coordinates. The center of the slit is the origin of the local coordinate system, the direction of which is the same as that of the global coordinate system. (0, y0, z0) are the coordinates of the slit center in the global coordinate system.
Fig. 3.
Fig. 3. Diagram of multiple OAC gratings
Fig. 4.
Fig. 4. Schematic showing ray tracing of multiple OAC gratings
Fig. 5.
Fig. 5. XY plane of spectral imaging
Fig. 6.
Fig. 6. Flowchart of construction stage of the spectrometer with multiple OAC gratings
Fig. 7.
Fig. 7. Spectrometer with multiple OAC gratings
Fig. 8.
Fig. 8. Distortion of spectrometer with multiple OAC gratings. (a)–(c) Spectral smiles for different wavelengths; (d)–(f) Spectral keystones for 0, 0.2, 0.4, 0.6, 0.8, and 1 normalized fields of view
Fig. 9.
Fig. 9. Full-field spot diagram of the imaging
Fig. 10.
Fig. 10. RMS spot radius diffracted by three OAC gratings
Fig. 11.
Fig. 11. (a) Spectral resolution diffracted by three OAC gratings; (b) Spectral resolution variation ratio
Fig. 12.
Fig. 12. SNR of each sub-imaging produced by the three OAC gratings

Tables (4)

Tables Icon

Table 1. Basic parameters of initial system

Tables Icon

Table 2. Five sets of optimization results

Tables Icon

Table 3. Parameters of spectrometer with multiple OAC gratings

Tables Icon

Table 4. Parameters of multiple OAC gratings

Equations (34)

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

tan φ = y p y 0 x p ,
tan θ = ( y p y 0 ) 2 + x p 2 z p z 0 .
z p = c p ( x p 2 + y p 2 ) 1 + 1 c p 2 ( x p 2 + y p 2 ) ,
z g = c g ( x g 2 + y g 2 ) 1 + 1 c g 2 ( x g 2 + y g 2 )  -  D .
( l p g | l p g | l s p | l s p | ) n p = 0
l p g ( l s p × n p ) = 0
n p = ( c p x p , c p y p , 1 c p 2 ( x p 2 + y p 2 ) ) ,
l s p = ( x p x 0 , y p y 0 , z p z 0 ) ,
l p g = ( x g x p , y g y p , z g z p ) .
n g = ( c g x g , c g y g , 1 c g 2 ( x g 2 + y g 2 ) ) .
θ i = arccos ( l p g n g | l p g | )
sin θ i  -  sin θ e , λ n = M λ n T ,
θ e , λ n = arccos ( l g t n g | l g t | ) .
l g t ( l p g × n g ) = 0.
z t , λ n = c t ( x t , λ n 2 + y t , λ n 2 ) 1 + 1 c t 2 ( x t , λ n 2 + y t , λ n 2 ) .
l t i ( l g t × n t , λ n ) = 0
( l g t | l g t | l t i | l t i | ) n t , λ n = 0
l g t  =  ( x t , λ n x g , y t , λ n y g , z t , λ n z g )
l t i = ( x i , λ n x t , λ n , y i , λ n y t , λ n , z i , λ n z t , λ n )
n t , λ n = ( c t x t , λ n , c t y t , λ n , 1 c t 2 ( x t , λ n 2 + y t , λ n 2 ) )
Δ p m = | y i , λ 1 , g m y i , λ n , g m | ,
Δ p m 1 = | y i , λ 1 , g m 1 y i , λ n , g m 1 | ,
Δ h m , m 1 = | y i , λ c , g m 1 y i , λ c , g m | ,
2 Δ h m , m 1 Δ p m 1  +  Δ p m .
D s m i l e , λ n , g m = y max , λ n , g m y min , λ n , g m
D k e y s t o n e , λ n , g m = x max , λ n , g m x min , λ n , g m
δ λ D , m = a D Δ p m | λ 1 λ n |
β D  =  δ λ D , min δ λ D , max = Δ p min Δ p max
δ λ R , m = a R , λ Δ p m | λ 1 λ n |
β D  =  δ λ R , min δ λ R , max
y p c , g 1 = 1 3 ϕ y + 1 2 ( y p max + y p min )
y p c , g 2 = 0
y p c , g 3 =  -  1 3 ϕ y + 1 2 ( y p max + y p min )
Δ p 1 + Δ p 2  +  Δ p 3 L
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.