In-Orbit Calibration of Phase Retardance for Channeled Spectropolarimeter

The phase retardance of the optical system (PROS) is a crucial factor limiting the accuracy of the Stokes vector reconstruction for the channeled spectropolarimeter. The dependence on reference light with a specific angle of polarization (AOP) and the sensitivity to environmental disturbance brings challenges to the in-orbit calibration of PROS. In this work, we propose an instant calibration scheme with a simple program. A function with a monitoring role is constructed to precisely acquire a reference beam with a specific AOP. Combined with numerical analysis, high-precision calibration without the onboard calibrator is realized. The simulation and experiments prove the effectiveness and anti-interference characteristics of the scheme. Our research under the framework of fieldable channeled spectropolarimeter shows that the reconstruction accuracy of S2 and S3 in the whole wavenumber domain are 7.2 × 10−3 and 3.3 × 10−3, respectively. The highlight of the scheme is to simplify the calibration program and ensure that the PROS high-precision calibration is not disturbed by the orbital environment.


Introduction
Polarimetric spectral imaging is a powerful tool with various applications [1], such as atmospheric aerosol characterization [2,3] and material analysis [4][5][6][7]. Channeled spectropolarimetry, proposed by Oka et al. [8,9], enables full Stokes vector acquisition with simple optical structure and fixed polarization components. All these advantages promote the development of channeled spectropolarimeters based on different structures [10,11]. Most of the reported schemes follow the assumption that the polarization aberrations of the optical systems have negligible impact on the measurement process and reconstruction by default. However, this idealized simplification is a concrete source of error in the context of large fields of view and the complexity of optical film design [12]. The phase retardance of the optical system (PROS) is a crucial factor. Specifically, the PROS characterizes the maximum phase difference in the eigenpolarization state. Due to the difference, the element produces different phase changes for beams with different polarization states. This effect is superimposed along the beam propagation path, limiting the accuracy of the reconstructed Stokes vectors. The PROS is quite sensitive to environmental disturbance [13]. Considering that the orbital environment will be affected by uncertainties such as temperature, launch vibration, and satellite attitude, in-orbit instant calibration of PROS is essential and meaningful.
So far, the quantitative analysis of PROS for the channeled spectropolarimeter has focused on the laboratory stage, requiring the assistance of a reference beam with a known polarization state. Stephen et al. [14] characterized the non-ideal factors of the system using linear operators to record the model response of the system by inputting reference beams with different angles of polarization (AOP). The corresponding data acquisition and processing are redundant. Yang, Xing et al. [15,16] simplified the process and considered the issue of polarization effects with radiation coupling. Linearly polarized beams at 22.5 • and 45 • were used for the quantitative calibration of PROS in the laboratory. However, it is challenging to achieve quantitative calibration of PROS in orbit. One factor is the dependence of the high-precision calibration on the reference beam. In existing in-orbit polarization calibration methods, the onboard calibrator with a built-in standard illuminant is the most effective way to acquire a specific reference beam. The supporting mechanical structure and control system adds to the bulk of the device and gets in the way of a compact implementation. The complexity of the components is also detrimental to long-term in-orbit maintenance. Another factor is that the onboard calibrator is independent of the observing system and can only achieve periodic calibrations. Environmental disturbances are out-of-synchrony for calibrations of different parameters.
In this work, we propose a novel scheme for the in-orbit calibration of the PROS for the channeled spectropolarimeter. A dynamic function with a monitoring role is constructed to precisely acquire a reference beam with a specific AOP at high frequencies to achieve the high-precision calibration of PROS without additional calibrators.

Polarization Radiometric Calibration Model of Channeled Spectropolarimeter
The optical schematic of the fieldable channeled spectropolarimeter designed for airborne remote sensing is depicted in Figure 1 [16]. A polychromatic beam is modulated by the polarimetric-spectral intensity modulation (P-SIM) [8] after the collection and collimation of the fore-optics. The P-SIM consists of high-order retarders R 1 , R 2 and a polarizer A. A dispersive imaging spectrometer receives the modulation spectrum. The phase retardance of the fore-optics δ f ore is analyzed as the only PROS source for this structure. polarization state. Stephen et al. [14] characterized the non-ideal factors of the system using linear operators to record the model response of the system by inputting reference beams with different angles of polarization (AOP). The corresponding data acquisition and processing are redundant. Yang, Xing et al. [15,16] simplified the process and considered the issue of polarization effects with radiation coupling. Linearly polarized beams at 22.5° and 45° were used for the quantitative calibration of PROS in the laboratory. However, it is challenging to achieve quantitative calibration of PROS in orbit. One factor is the dependence of the high-precision calibration on the reference beam. In existing inorbit polarization calibration methods, the onboard calibrator with a built-in standard illuminant is the most effective way to acquire a specific reference beam. The supporting mechanical structure and control system adds to the bulk of the device and gets in the way of a compact implementation. The complexity of the components is also detrimenta to long-term in-orbit maintenance. Another factor is that the onboard calibrator is independent of the observing system and can only achieve periodic calibrations. Environmental disturbances are out-of-synchrony for calibrations of different parameters.
In this work, we propose a novel scheme for the in-orbit calibration of the PROS for the channeled spectropolarimeter. A dynamic function with a monitoring role is constructed to precisely acquire a reference beam with a specific AOP at high frequencies to achieve the high-precision calibration of PROS without additional calibrators.

Polarization Radiometric Calibration Model of Channeled Spectropolarimeter
The optical schematic of the fieldable channeled spectropolarimeter designed for airborne remote sensing is depicted in Figure 1 [16]. A polychromatic beam is modulated by the polarimetric-spectral intensity modulation (P-SIM) [8] after the collection and collimation of the fore-optics. The P-SIM consists of high-order retarders  is the in-orbit acquisition of the reference beam with a specific AOP. We chose to "capture" the reference beam rather than "create" it. In the framework of the polarization radiometric transmission model, the corresponding relationship between the grayscale value of the spectrometer and the incident pupil radiance is expressed as: is the modulation spectrum radiance, characterized by [12,13,17]: The crux of the calibration of δ f ore is the in-orbit acquisition of the reference beam with a specific AOP. We chose to "capture" the reference beam rather than "create" it. In the framework of the polarization radiometric transmission model, the corresponding relationship between the grayscale value of the spectrometer and the incident pupil radiance is expressed as: where DN orbit is the original digital number DN; A k and C are the radiometric calibration coefficient and the DN of dark current, both of which can be accurately calibrated prior to orbit by mature absolute radiometric calibration. Furthermore, B orbit is the modulation spectrum radiance, characterized by [12,13,17]: where nA, linearly related to A k , is the constant of the system effective gain coefficient; M f ore , M PSI M , M imag , and M spec , correspond to the cascaded Muller matrix of four subsystems: fore-optics, P-SIM, imaging optics, and spectrometer, respectively. S in−orbit = [S 0 , S 1 , S 2 , S 3 ] T is the Stokes vector of the incident pupil beam.

Methods for In-Orbit Calibration of the PROS
Here, we set some new parameters for the intuitive demonstration of derivation [16]: ( It should be noted that: ε 1 and ε 2 characterize the angular difference between the actual setting position and the preset position of R 1 and R 2 . In current mainstream cognition, researchers tend to treat them as error sources. In practice, however, these two parameters are deeply involved in the decoupling process, which will directly impact calibration accuracy. The following work will prove that they cannot simply be set to small quantities. D f ore−orbit and δ f ore−orbit correspond to the diattenuation and phase retardance of the fore-optics set in the orbital environment.
To intuitively display the superposition state of the modulation spectrum in the optical path difference domain, we express it in the form of phase factor accumulation, as follows [16]: where δ 1,orbit and δ 2,orbit are the ideal retardances of R 1 and R 2 and the nine phase factors correspond to the nine distinct channels of the modulation spectrum. The channels can be obtained through Fourier transform: It should be noted that D f ore−orbit can be calibrated with the beams without a polarizer, as in the laboratory method, which is expressed as follows [16]: where G 0 and G 1 are the corresponding channels with the beams without a polarizer. The water cloud zero degree of linear polarization (DOLP) calibration source [18,19] or the sun (usually regarded as lambertian) [13,20] are ideal calibration sources for obtaining the beams without a polarizer. ε 1 , ε 2 , δ 1,orbit , and δ 2,orbit can be calibrated accurately by the method in Refs. [13,21]. The overall calibration scheme consists of multiple parameter calibrations, which are usually coupled. Based on the calibration sequence, these parameters, as the observed specific values, are used for numerical analysis of δ f ore−orbit . The Stokes vector reconstruction model is expressed as [16]: The overall calibration is based on the numerical operation of channel observations. Further, the construction sequence is unilinear. The channel observations and the calibrated parameters, can be used for subsequent analysis. Fortunately, we find a function to express the characteristic relationship between the AOP of the beam and the known quantity. For any partially linearly polarized beam S in−orbit = S 0 [1, P sin 2θ, P cos 2θ, 0] T with a known degree of polarization (DOP) P, the function is written as: We name Z as Target-beam Monitoring Function (TMF). Figure 2 shows the numerical solution of TMF plotted within the unit value interval of the independent variables ε 1 and θ.
The overall calibration is based on the numerical operation of channel observations. Further, the construction sequence is unilinear. The channel observations and the calibrated parameters, can be used for subsequent analysis. Fortunately, we find a function to express the characteristic relationship between the AOP of the beam and the known We name Z as Target-beam Monitoring Function (TMF). Figure 2 shows the numerical solution of TMF plotted within the unit value interval of the independent variables 1 ε and θ: Furthermore, the maximum is the specific value 1, which is direct feedback to "capture" the target beam with the specific AOP. Even though there are reasonable errors in the parameter The only maximum condition of TMF is θ = ε 1 , excluding the impossibility of ε 1 = 0 • (perfect installation) and ε 1 = 180 • (reverse installation) in practical projects. Furthermore, the maximum is the specific value 1, which is direct feedback to "capture" the target beam with the specific AOP. Even though there are reasonable errors in the parameter calibrations, the acquisition of the beam with a specific AOP is not affected, which ensures the in-orbit calibration accuracy of the PROS.
To increase the calibration frequency, we choose the marine flares, common in orbit, as the calibration area. The flares, as ideal reflective media, have high radiation utilization. Moreover, according to Fresnel reflection, when the beam's incident angle to the flares is greater than Brewster's angle, the reflected beam is fully linearly polarized. Further, considering natural factors such as wind speed and foam, a more accurate rough sea surface model is required in practice. In the relevant research of many scholars, we select the anisotropic Breon and Henriot model (BHA) [22]. The BHA uses reflectance data over the global oceans taken by the Polarization and Directionality of the Earth's Reflectances (POLDER) and wind data from NASA Scatterometer instruments. It has the highest correlation coefficient with the analytic DOP in the operating band. Assisted by the modified vector radiative transfer model and uncertainty analysis, the analytic DOP of the received polarized beam is fed back. The accuracy of the analytic DOP is controlled within 2%, meeting the calibration requirements. The above conditions endow the analytic solution of Z with instant parameters.
Using the "captured" reference beam, we obtain the in-orbit calibration value of δ f ore as follows:

Numerical Simulation
In the numerical simulation, the subsystem is modularized based on the cascade rule of the Muller matrix. The Muller matrix of each subsystem can be obtained by polarization ray tracing. The operating band range, the high-order retarder's thickness, the fore-optic's structure, and the detector selection is referred to in Ref. [16]. The relevant data is used as the input value of the simulation. To strictly depict the integrity of Z, the tuning range of the DOP and the AOP of the simulated illuminant covers all possible conditions. The demonstration of calibration in a continuous band is complex and abstract. Therefore, the center of the operating band is selected as the observation wavelength if necessary. In addition, the direct effect of pupil distribution on polarization aberration cannot be ignored. Therefore, according to the pupil distribution of δ f ore , we choose the location with a high value (normalized field of view [0, 0.5]) as the observation point. ε 1 is set to 1 • , 3 • , 30 • , and 45 • .
The array of Z analytical values varying with the simulated illuminant is shown in Figure 3. In different cases, the analytical trend of Z is highly consistent with the calculation results based on the theoretical model. Furthermore, we reached some interesting conclusions: (a) The DOP is not coupled with other parameters as shown in Figure 3e. The error of the analytic DOP is not a factor affecting the analytical characteristics of Z; (b) Compared with Figure 3c,d, the maximum of Z in Figure 3a,b are too "dense" to be distinguished, which corresponds to ε 1 being considered as only caused by the alignment error. It proves that tiny ε 1 will directly impact the sensitivity and effectiveness of the TMF. It is necessary to consider making the angular difference between the actual setting position and the preset position of high-order retarders easier to distinguish. We compare the simulated in-orbit calibration method of δ f ore with the conventional laboratory calibration method at ε 1 = 30 • . Figure 4 compares calibration errors at 10,500-23,450 nm −1 . Sensors 2023, 23, x FOR PEER REVIEW 7 of 10  Due to the periodicity of the Fourier transform algorithm, the apodization function is used to suppress the ringing effect that may occur in the reconstruction. This processing inevitably leads to edge data exceptions. When this interference is eliminated, the calibration accuracy of fore orbit δ − and fore lab δ − are highly consistent.

Experimental Verification
The instrument is in the pre-research stage. The approach to in-orbit test data is not yet available. Therefore, we built the configuration shown in Figure 5 for the calibration experiments. The functional unit consists of the PSIM and a spectrometer (Field Spec 3, Analytical Spectral Devices). The input data is consistent with the simulation.   Due to the periodicity of the Fourier transform algorithm, the apodization function is used to suppress the ringing effect that may occur in the reconstruction. This processing inevitably leads to edge data exceptions. When this interference is eliminated, the calibration accuracy of fore orbit δ − and fore lab δ − are highly consistent.

Experimental Verification
The instrument is in the pre-research stage. The approach to in-orbit test data is not yet available. Therefore, we built the configuration shown in Figure 5 for the calibration experiments. The functional unit consists of the PSIM and a spectrometer (Field Spec 3, Analytical Spectral Devices). The input data is consistent with the simulation.
(e) Due to the periodicity of the Fourier transform algorithm, the apodization function is used to suppress the ringing effect that may occur in the reconstruction. This processing inevitably leads to edge data exceptions. When this interference is eliminated, the calibration accuracy of δ f ore−orbit and δ f ore−lab are highly consistent.

Experimental Verification
The instrument is in the pre-research stage. The approach to in-orbit test data is not yet available. Therefore, we built the configuration shown in Figure 5 for the calibration experiments. The functional unit consists of the PSIM and a spectrometer (Field Spec 3, Analytical Spectral Devices). The input data is consistent with the simulation.
The in-orbit scene is restored as much as possible. A variable polarizer is added to the illuminant-generating device. The polarization state of the target beam is obtained by Stokes polarimeter. In addition, possible ambient temperature disturbance in orbit is of consideration. The configuration is placed in the temperature control system. The regulating range is 20 • C to 30 • C, which is consistent with the temperature control requirements of the aircraft.
In the validation experiment of TMF, we set ε 1 and θ as adjusting variables. Based on the conclusion of the numerical simulation, Figure 6 shows the measurement results with a fully linearly polarized beam (P = 1 ). When the polarization state of the target light changes continuously, the response distribution of TMF is consistent with the simulation results. In other words, TMF strictly has the only maximum condition "θ = ε 1 , Z= 1". The The in-orbit scene is restored as much as possible. A variable polarizer is added to the illuminant-generating device. The polarization state of the target beam is obtained by Stokes polarimeter. In addition, possible ambient temperature disturbance in orbit is of consideration. The configuration is placed in the temperature control system. The regulating range is 20 °C to 30 °C, which is consistent with the temperature control requirements of the aircraft.
In the validation experiment of TMF, we set 1 ε and θ as adjusting variables. Based on the conclusion of the numerical simulation, Figure 6 shows the measurement results with a fully linearly polarized beam ( =1 P ). When the polarization state of the target light changes continuously, the response distribution of TMF is consistent with the simulation results. In other words, TMF strictly has the only maximum condition " ". The maximum is the specific value 1. The proposed method of "capturing" the target beam with a specific AOP is further verified. In the multi-parameter coupling calibration scheme, it is not easy to highlight the influence of a single parameter. Therefore, we take the reconstruction accuracy of the Stokes vector as the traceability benchmark. The fore δ is coupled to the reconstruction model of S2 and S3. Figure 7 compares the normalized reconstructed Stokes vector with the input value corresponding to different temperature settings. The reconstruction accuracy of S2 and S3 in the whole wavenumber domain by RMSE is 7.2 × 10 −3 and 3.3 × 10 −3 .  The in-orbit scene is restored as much as possible. A variable polarizer is added to the illuminant-generating device. The polarization state of the target beam is obtained by Stokes polarimeter. In addition, possible ambient temperature disturbance in orbit is of consideration. The configuration is placed in the temperature control system. The regulating range is 20 °C to 30 °C, which is consistent with the temperature control requirements of the aircraft.
In the validation experiment of TMF, we set 1 ε and θ as adjusting variables. Based on the conclusion of the numerical simulation, Figure 6 shows the measurement results with a fully linearly polarized beam ( =1 P ). When the polarization state of the target light changes continuously, the response distribution of TMF is consistent with the simulation results. In other words, TMF strictly has the only maximum condition " ". The maximum is the specific value 1. The proposed method of "capturing" the target beam with a specific AOP is further verified. In the multi-parameter coupling calibration scheme, it is not easy to highlight the influence of a single parameter. Therefore, we take the reconstruction accuracy of the Stokes vector as the traceability benchmark. The fore δ is coupled to the reconstruction model of S2 and S3. Figure 7 compares the normalized reconstructed Stokes vector with the input value corresponding to different temperature settings. The reconstruction accuracy of S2 and S3 in the whole wavenumber domain by RMSE is 7.2 × 10 −3 and 3.3 × 10 −3 . In the multi-parameter coupling calibration scheme, it is not easy to highlight the influence of a single parameter. Therefore, we take the reconstruction accuracy of the Stokes vector as the traceability benchmark. The δ f ore is coupled to the reconstruction model of S2 and S3. Figure 7 compares the normalized reconstructed Stokes vector with the input value corresponding to different temperature settings. The reconstruction accuracy of S2 and S3 in the whole wavenumber domain by RMSE is 7.2 × 10 −3 and 3.3 × 10 −3 .
The applicability of the scheme is analyzed as necessary. The scheme is not limited to adapting specific channeled spectropolarimeters or PROS in-orbit calibration. The grayscale value received by the spectrometer is added to the linear calibration operation as the only external input. It is ruled out that the polarization information of the polychromatic beam may cause crosstalk to channel construction. Therefore, for all the channeled spectropolarimeters characterized with the Muller matrix, TMF can be constructed similarly to accurately obtain beams with specific AOP. The beams with different AOPs can be obtained by adjusting the setting angle of high-order retarders as required. It provides a novel idea for constructing the overall in-orbit calibration scheme. It should be noted that a single reference beam is only valid in one calibration cycle. The effective way to improve the calibration accuracy is to increase the calibration frequency. The corresponding verification work will be supported by the in-orbit experimental data. The applicability of the scheme is analyzed as necessary. The scheme is not l to adapting specific channeled spectropolarimeters or PROS in-orbit calibration. The scale value received by the spectrometer is added to the linear calibration operation only external input. It is ruled out that the polarization information of the polychro beam may cause crosstalk to channel construction. Therefore, for all the channeled tropolarimeters characterized with the Muller matrix, TMF can be constructed sim to accurately obtain beams with specific AOP. The beams with different AOPs can tained by adjusting the setting angle of high-order retarders as required. It prov novel idea for constructing the overall in-orbit calibration scheme. It should be note a single reference beam is only valid in one calibration cycle. The effective way to im the calibration accuracy is to increase the calibration frequency. The corresponding fication work will be supported by the in-orbit experimental data.

Conclusions
In sum, a novel scheme of PROS in-orbit calibration for the channeled spectro imeter is presented. To avoid the burden of the volume and design difficulty onboard calibrator, we build the Target-beam Monitoring Function to capture the with a specific AOP. In conjunction with the general derivation of the PROS, hig quency calibration is realized. The numerical simulation and experiments prove the tion's validity and the high calibration precision. The highlight of the scheme is to sim the calibration program and ensure that the PROS high-precision calibration is n turbed by the orbital environment.

Conclusions
In sum, a novel scheme of PROS in-orbit calibration for the channeled spectropolarimeter is presented. To avoid the burden of the volume and design difficulty of the onboard calibrator, we build the Target-beam Monitoring Function to capture the beam with a specific AOP. In conjunction with the general derivation of the PROS, high-frequency calibration is realized. The numerical simulation and experiments prove the function's validity and the high calibration precision. The highlight of the scheme is to simplify the calibration program and ensure that the PROS high-precision calibration is not disturbed by the orbital environment.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.