Polarization Radiometric Calibration in Laboratory for a Channeled Spectropolarimeter

: The process of radiometric calibration would be coupled with the polarization properties of an optical system for spectropolarimetry, which would have signiﬁcant inﬂuences on reconstructed Stokes parameters. In this paper, we propose a novel polarization radiometric calibration model that decouples the radiometric calibration coe ﬃ cient and polarization properties of an optical system. The alignment errors of the polarization module and the variation of the retardations at di ﬀ erent ﬁelds of view are considered and calibrated independently. According to these calibration results, the input Stokes parameters at di ﬀ erent ﬁelds of view can be reconstructed accurately through the proposed model. Simulations are performed for the presented calibration and reconstruction methods, which indicate that the measurement accuracy of polarization information is improved compared with the traditional undecoupled calibration method.


Introduction
Imaging spectrometer and imaging polarimeter are of particular importance in many fields, including atmospheric aerosol characterization [1,2], remote sensing [3][4][5][6], and material characterization [7,8]. These fields require an imaging spectropolarimeter to simultaneously measure the spatial, spectral, and polarization information of a scene. In 1999, Oka and Iannarilli first described the channeled spectropolarimetry technique, an attractive approach for polarimetry [9,10] by which we can eliminate the disadvantages of conventional modulation principles. With a simple optical system and no movable polarization components, the entire wavelength-dependent state of polarization (SOP) and spectral information of a scene can be acquired simultaneously.
In recent publications, researchers have investigated incorporating the channeled spectropolarimetry technique into different types of imaging spectrometers, such as the dispersive imaging spectrometer [3,4,11], Fourier transform imaging spectrometer [6,12], and computed tomography imaging spectrometer [13,14]. Among these incorporations, the channeled spectropolarimetry technique with the dispersive imaging spectrometer, which is called the channeled dispersive imaging spectropolarimeter (CDISP), has a simple structure and is easier to implement in terms of technology. Noteworthy is that the CDISP should go through a series of accurate calibration before quantitative applications. For imaging instruments, the actual spectral radiance of the detected object should be obtained through radiometric calibration. For the imaging polarimeter, the target light from the detected object will be affected by the polarization properties of an optical system when transmitted in the optical system. This will change the polarimetric state of target light and affect the measurement accuracy of polarization information. Therefore, in the process of radiometric calibration for the imaging spectropolarimeter, the polarization properties of an optical system should be considered.
However, to the best of our knowledge, radiometric calibration for the imaging spectropolarimeter have not be well solved with the previous literature. Chunmin Zhang et al. implemented radiometric calibration and polarimetric calibration for the channeled interference imaging spectropolarimeter [15], while the two types of calibration are carried out separately. The radiometric calibration is the same as that of traditional imaging spectrometer, without taking into account the influence of the polarization properties of optical system. Michael W. Kudenov et al. provided the polarimetric calibration of the channeled polarimeter with a pushbroom hyperspectral imaging sensor only [4], without considering the radiometric calibration of the instrument. Bin Yang et al. presented a method of polarimetric calibration and reconstruction for the fieldable CDISP [16]. This method took into account the polarization properties of the optical system, but only if it can accurately measure the spectral radiance of the detected object.
The interaction between the polarized beam and the optical system could be called polarization radiation transmission. Modeling and analyzing the polarization radiation transmission properties of optical system is the core to obtain the polarization properties of optical system and the theoretical basis for calibration of polarization detection system. In this paper, we use a fieldable CDISP designed for remote sensing as an example to analyze the influence of polarization properties of the optical system on radiometric calibration of the imaging spectropolarimeter. Considering simultaneously alignment errors of the polarimetric spectral intensity modulation (PSIM) module and variation of the retardations at different fields of view, we establish a polarization radiometric calibration model for the fieldable CDISP. Formerly, we proposed a model of self-correction with a 3:1 ratio of retarders to calibrate and compensate alignment errors and retardations simultaneously by measuring the SOP of target light in Ref. [17]. In order to maintain the advantages of this calibration method, the polarization radiometric calibration model proposed in this paper still adopts the scheme with a 3:1 ratio of retarders. Based on this model, we present a reconstruction method of Stokes parameters and a series of parameters calibration methods, classical radiometric calibration, self-correction technique, polarization radiometric decoupling calibration, and so on. By using the presented methods, the radiometric calibration coefficient and polarization parameters are determined accurately, and the reconstruction accuracy of the polarization information could be improved obviously, which is significant for the practical quantitative applications of the fieldable CDISP. This paper is organized as follows. Section 2 provides the principle of fieldable channeled dispersive imaging spectropolarimetry, and then, derive a theoretical model for the polarization radiometric calibration and the reconstruction method of Stokes parameters. Section 3 proposes a series of calibration methods of parameters in the polarization radiometric calibration model. In Section 4, we analyze and summary our simulation experiment results. The conclusion is presented in Section 5.

Establishment of Polarization Radiometric Calibration Model of the Fieldable Channeled Dispersive Imaging Spectropolarimetry
In this section, we first introduce the principle of fieldable channeled dispersive imaging spectropolarimetry and its four sub-modules. Then, according to modularity theory and the principle of classical radiometric calibration, we derive a theoretical model for the polarization radiometric calibration.

Principle of the Fieldable Channeled Dispersive Imaging Spectropolarimetry
The optical schematic of the fieldable CDISP is depicted in Figure 1. In order to facilitate the establishment of polarization radiometric transmission model for optical system, we advanced a method modeling to divide into four sub-modules: fore-optics, PSIM module, imaging optics, and dispersive imaging spectrometer. Polarization information from a detected object is collected and collimated by the fore-optics, and then modulated by PSIM module. After passing though the imaging optics, the modulation spectrum is finally received by a dispersive imaging spectrometer. The fore-optics is used to meet the requirement that all incident angles on R 1 and R 2 are limited to less than about 5 • in the PSIM module [3], and to reduce the size of polarization elements and the processing difficulty.

Principle of the Fieldable Channeled Dispersive Imaging Spectropolarimetry
The optical schematic of the fieldable CDISP is depicted in Figure 1. In order to facilitate the establishment of polarization radiometric transmission model for optical system, we advanced a method modeling to divide into four sub-modules: fore-optics, PSIM module, imaging optics, and dispersive imaging spectrometer. Polarization information from a detected object is collected and collimated by the fore-optics, and then modulated by PSIM module. After passing though the imaging optics, the modulation spectrum is finally received by a dispersive imaging spectrometer. The fore-optics is used to meet the requirement that all incident angles on R1 and R2 are limited to less than about 5° in the PSIM module [3], and to reduce the size of polarization elements and the processing difficulty. The PSIM module consists of two high-order retarders, R1 and R2, with thicknesses d1 and d2, respectively, and a linear analyzer, A. The fast axis of R1 is aligned with the transmission axis of A, and R2 is oriented with its fast axis at 45° to the transmission axis of A in an ideal condition. According to Mueller matrix theory and radiation transmission theory, if the incident Stokes vector is The polarization properties of optical system include diattenuation, retardance and depolarization. Thereinto, depolarization is intrinsically associated with scattering and a loss of coherence in the polarization state. A small amount of depolarization is associated with the scattered light from all optical components [18]. Therefore, only diattenuation and retardance were studied as the polarization properties of CDISP.
The fore-optics and imaging optics can be regarded as refractors. According to the Mueller matrix of each lens interface of optics and its series properties [18], the equivalent Mueller matrices of the fore-optics and imaging optics at the field of view θ can be expressed as The PSIM module consists of two high-order retarders, R 1 and R 2 , with thicknesses d 1 and d 2 , respectively, and a linear analyzer, A. The fast axis of R 1 is aligned with the transmission axis of A, and R 2 is oriented with its fast axis at 45 • to the transmission axis of A in an ideal condition. According to Mueller matrix theory and radiation transmission theory, if the incident Stokes vector is S in = [S 0 , S 1 , S 2 , S 3 ] T (where T denotes the transpose operation), the Stokes vector of the target light passing through the PSIM module can be given by where M R 2 , and M R 1 denote the Mueller matrices of A, R 2 , and R 1 , respectively [18], and ϕ 2 (σ) and ϕ 1 (σ) are the ideal phase retardations of R 2 and R 1 , respectively. However, for a fieldable CDISP, ϕ 2 (σ) and ϕ 1 (σ) are related to the incident angle of the ray to the high-order retarders and will change with the increase of the field angle. In addition, alignment errors, ε 1 and ε 2 of R 1 and R 2 , are inevitable for the manufacturing and alignment tolerances, which have significant influences on the reconstructed Stokes parameters. The Mueller matrix of the PSIM module becomes The polarization properties of optical system include diattenuation, retardance and depolarization. Thereinto, depolarization is intrinsically associated with scattering and a loss of coherence in the polarization state. A small amount of depolarization is associated with the scattered light from all optical components [18]. Therefore, only diattenuation and retardance were studied as the polarization properties of CDISP. The fore-optics and imaging optics can be regarded as refractors. According to the Mueller matrix of each lens interface of optics and its series properties [18], the equivalent Mueller matrices of the fore-optics and imaging optics at the field of view θ can be expressed as where A f ore (σ, θ) and A imag (σ, θ) are the transmittance of the fore-optics and imaging optics, respectively. D f ore (σ, θ) and D imag (σ, θ) are the diattenuation of the fore-optics and imaging optics, respectively. δ f ore (σ, θ) and δ imag (σ, θ) are the retardance of the fore-optics and imaging optics, respectively. Similarly, the dispersive imaging spectrometer can be regarded as a reflector, which is composed of many reflecting interfaces. The equivalent Mueller matrices at the field of view θ is given by where A spec (σ, θ), D spec (σ, θ) and δ spec (σ, θ) are the transmittance, diattenuation and retardance of the dispersive imaging spectrometer, respectively. For the fieldable CDISP, these polarization properties increase with the increase of field angle, so they will affect the measurement results and can not be ignored. In this system, the slit is placed along the X-axis, and only the field of view in X-direction exists. The azimuth of incident plane of each viewing chief ray is a fixed value of 0 • , not a variable. Therefore, the Mueller matrices of sub-module would not be affected by azimuth angle.

Establishment of Polarization Radiometric Calibration Model
According to the principle of classical radiometric calibration [19,20], the relationship between the digital number (DN) of CCD detector and the spectral radiance at entrance pupil in the actual operation of the instrument can be given by where DN(σ, θ) is the original DN of the CCD detector, A(σ, θ) represents the detector pixel response, and C(σ, θ) is the DN of CCD dark current. For a dispersive spectropolarimeter, (σ, θ) represents the positions of different pixels on the CCD detector. And is a function of the transmittance and diattenuation of sub-modules. Their values are different at different pixel positions of the CCD detector, that is, they change with the wavenumber and field of view. Therefore, we simply need to determine A k (σ, θ) with no need for determining the transmittance and diattenuation of sub-modules independently in the calibration.
T is the spectral radiance at entrance pupil of the instrument. The spectral radiance at the entrance pupil, contains spectral, intensity and polarization information of the detected object, is modulated by the PSIM module into modulation spectrum. The modulation spectrum radiance could be obtained by radiometric calibration of the DN of CCD detector, which can be expressed as where and the superscript * denotes the complex conjugate. Based on Equation (8), the nine phase terms represent nine distinct channels in Fourier domain which contain the data related to the input Stokes parameters S 0 (σ), S 1 (σ), and S 3 (σ). Calculating the autocorrelation of B(σ, θ) using the Fourier transformation, we can distribute the nine channels with nine phase terms where L denotes the variable conjugate to wavenumber σ under the Fourier transform, and the superscript θ is the field angle. L j = ∆n(σ)d j , is the optical path difference (OPD) introduced by R j , ∆n(σ) is the birefringence of the birefringent crystal. By filtering out nine channels and performing inverse Fourier transformations independently, the whole Stokes parameters of the target light can be described as where X i (σ, θ) can be obtained by reference to the self-correction method. The interested reader can find more information in [17]. From Equation (6), we can find that the polarization properties of imaging optics and dispersive imaging spectrometer will not influence the measurement of the CDISP. The process of measuring polarization information is complicated due to diattenuation and retardance of the fore-optics, D f ore (σ, θ) and δ f ore (σ, θ). In addition, this effect will be more serious with the increase of field angle. However, the fore-optics is indispensable because of its important role in the fieldable CDISP [16], so we need to calibrate the polarization properties of fore-optics accurately. Based on the derived theoretical model, the parameters, , and δ f ore (σ, θ), are needed to be calibrated accurately to reconstruct the Stokes parameters accurately.

Method of Polarization Radiometric Calibration
In the presented methods, we proposed a series of calibration methods of parameters in the polarization radiometric calibration model. We first briefly expound the classical radiometric calibration method to obtain the dark current response C(σ, θ) and modulated radiometric calibration coefficient A k (σ, θ). Then the polarization parameters independent of field of view are calibrated with self-correction technique. Different from the radiometric calibration of traditional imaging spectrometer, the radiometric calibration coefficient and polarization parameters would be coupled in the calibration process for the fieldable CDISP. Therefore, we come up with a polarization radiometric decoupling calibration method to achieve the accurate calibration of the radiometric calibration coefficient A k (σ, θ) and diattenuation of fore-optics D f ore (σ, θ). Finally, exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] and δ f ore (σ, θ) are acquired by using the 0 • and 45 • linearly polarized reference beam.

Classical Radiometric Calibration
According to the classical radiometric calibration method [19], the dark current response C(σ, θ) and modulated radiometric calibration coefficient A k (σ, θ) can be calibrated by using the lights generated by integrating sphere. The integrating sphere provides unpolarized polychromatic lights that S 0 (σ) 0,S 1 (σ) = 0,S 2 (σ) = 0, and S 3 (σ) = 0. Substituting these equations into Equation (6), we know that, due to the existence of polarization properties, the calibrated A k (σ, θ) is not A k (σ, θ) as we needed, but the result of its coupling with the polarization parameters. A k (σ, θ) is given by Based on Equation (16), the radiometric calibration coefficient is coupled with the polarization properties in the process of radiometric calibration for fieldable CDISP. It is necessary to decouple them to achieve accurate polarization radiometric calibration. Before decoupling them, the basic polarization parameters need to be calibrated.

Self-Correction of Alignment Errors ε 1 ,ε 2 and Retardation exp[iϕ 2 (σ)]
The alignment errors of two high-order retarders do not change with the field of view. In addition, the variation of the phase retardation of R 2 with the field of view is too small, so that it can be ignored [16], so we can use the value of ϕ 2 (σ) at the central fields of view (FOV) as the values of the entire FOV. At central FOV, the beam light is incident vertically to the interface of each lens of fore-optics, so its polarization properties can be negligible, D f ore (σ, 0) ≈ 0 and δ f ore (σ, 0) ≈ 0.
Substituting these into the Equation (16), we obtain that A k (σ, 0) = A k (σ, 0). Therefore, the modulated radiometric calibration coefficient of central FOV can be used to acquire the modulation spectrum radiance after passing the PSIM module. Based on the calibrated modulation spectrum radiance, the above polarization parameters, ε 1 ,ε 2 , and exp[iϕ 2 (σ)], can be obtained by using the self-correction method [17]. A k (σ, θ) and Diattenuation of Fore-Optics D f ore (σ, θ)

Decoupling of Radiation Calibration Coefficient
In the process of radiometric calibration, calibration coefficient A k (σ, θ) and polarization property D f ore (σ, θ) are modulated to different frequencies by two high-order retarders. The Fourier transformation of A k (σ, θ) gives the autocorrelation function, which can be expressed as By filtering out the desired channels and performing inverse Fourier transformations independently, the radiometric calibration coefficient A k (σ, θ) and diattenuation of fore-optics D f ore (σ, θ) can be described as where a i , b i , and c i can be computed through ε 1 and ε 2 calibrated in Section 3.2, and exp[iϕ 2 (σ)] as computed in Section 3.2.
For the traditional imaging spectropolarimeter, the influence of polarization properties of optical system on radiometric calibration is not considered. A k (σ, θ) is directly taken as the radiometric calibration coefficient, so that there is a deviation between the spectral radiance obtained by calibration and the spectral radiance at the entrance pupil.

Calibration of Retardation exp
and Retardance of Fore-Optics δ f ore (σ, θ) The variation of the phase retardation of R 1 with the field of view will significantly influence the measurement accuracy of the fieldable CDISP [16], so ϕ 1 (σ, θ) cannot be regarded as a constant across the entire FOV. In addition, because the phase retardations of R 1 is always in the form of exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] in the calibration and reconstruction of Stokes parameters, we just need to To calibrate exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] with the change of field angle, a 0 • linearly polarized beam is entered into the system and filled the entire entrance pupil. For a 0 • linearly polarized beam, S 0 (σ) = S 1 (σ), S 2 (σ) = 0, and S 3 (σ) = 0. 0 and the 4 channels in Equation (11) then become where C θ 0,0 • and C θ 4,0 • denote the channels obtained when the 0 • linearly polarized beam is entered into entire entrance pupil of the instrument. By filtering out the desired channels and performing inverse Fourier transformations independently, exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] can be described as For the fieldable CDISP, the polarization properties caused by the fore-optics is more serious with the increase of the field angle. The retardance of fore-optics, is different with the change of field angle. To calibrate δ f ore (σ, θ), a 45 • linearly polarized beam is entered into the system and filled the entire entrance pupil. For a 45 • linearly polarized beam, S 0 (σ) = S 2 (σ), S 1 (σ) = 0, and S 3 (σ) = 0, Substituting these equations into Equation (11), we can obtain The channels C θ 0,45 • and C θ 4,45 • are extracted by the frequency filtering technique and performed Fourier transformation independently. Then, the expression of δ f ore (σ, θ) can be obtained where Im means the operator to take the imaginary part, D f ore (σ, θ) and exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] was obtained in the previous calibration.
, and δ f ore (σ, θ) will change with field of view, we should calibrate them at different field of view independently. To improve the efficiency of the calibration, we can calibrate them at part of all field of view, and get the calibration results across the entire FOV through curve fitting method. At this point, the accurate calibration of the parameters in the polarization radiometric calibration model is completed.

Verification by Simulation Experiment
To verify the feasibility and advantages of the presented method, simulations experiment for the calibration of parameters in the polarization radiometric calibration model and reconstruction of Stokes parameters are performed. In the simulations, we use a fieldable CDISP designed for airborne remote sensing in advance, whose schematic layout is shown in Figure 2. It operates over 13,000-18,000 cm −1 , the FOV is 40 • and the F-number is 4. R 1 and R 2 are made of quartz, whose birefringence in the selected wave band can be referred to Ref. [21], and their thicknesses are 6mm and 2 mm, respectively. We analyze the polarization properties of each sub-modules through polarization ray tracing [22,23] and calculate the phase retardations of R 2 and R 1 at different field of view. These results are used as the input values of the simulations. We compare the calibration results of parameters in the polarization radiometric calibration model and reconstructed Stokes parameters with two methods, i.e., the presented method in this paper, and the traditional calibration method of spectropolarimeter. Both methods take into account the alignment errors of polarization module and the variation of the retardations at different field of view in the process of radiometric calibration. The only difference is that in the radiometric calibration, the presented method considers the influence of polarization properties of optical system and decouples it, while the traditional calibration method does not consider its influence.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 17 respectively. We analyze the polarization properties of each sub-modules through polarization ray tracing [22,23] and calculate the phase retardations of R2 and R1 at different field of view. These results are used as the input values of the simulations. We compare the calibration results of parameters in the polarization radiometric calibration model and reconstructed Stokes parameters with two methods, i.e., the presented method in this paper, and the traditional calibration method of spectropolarimeter. Both methods take into account the alignment errors of polarization module and the variation of the retardations at different field of view in the process of radiometric calibration. The only difference is that in the radiometric calibration, the presented method considers the influence of polarization properties of optical system and decouples it, while the traditional calibration method does not consider its influence. We first simulate the process of radiometric calibration of the fieldable CDISP. The detector uses SARNOFF's CAM512 CCD as an example, whose area array is 512 × 512, the pixel size is 18 18 , μ μ × m m the frame frequency is 100 fps, and the quantum efficiency curve is shown in Figure 3.
According to the imaging chain analysis of the radiance from optical system entrance pupil to detector, the DN recorded by detector can be obtained by simulation. Using the radiance of the optical system entrance pupil and the simulated DN of detector, the dark current response ( )  We first simulate the process of radiometric calibration of the fieldable CDISP. The detector uses SARNOFF's CAM512 CCD as an example, whose area array is 512 × 512, the pixel size is 18 µm × 18 µm, the frame frequency is 100 fps, and the quantum efficiency curve is shown in Figure 3. According to the imaging chain analysis of the radiance from optical system entrance pupil to detector, the DN recorded by detector can be obtained by simulation. Using the radiance of the optical system entrance pupil and the simulated DN of detector, the dark current response C(σ, θ) and modulated radiometric calibration coefficient A k (σ, θ) can be obtained on the basis of the classical radiometric calibration principle.
Next, we calibrate the alignment errors ε 1 , ε 2 and retardation ϕ 2 (σ) of the PSIM module by using the self-correction method. The alignment errors of R 1 and R 2 are ε 1 = −0.5 • , ε 2 = 0.5 • , respectively. Because the alignment errors are independent of wavenumber, we can use the averages of the calculated values in different wavenumbers as the final determination results, which are both −0.5029 • and 0.5007 • by the presented method and the traditional calibration method. The input values and calibration results of ϕ 2 (σ) at different field of view is shown in Figure 4. The field of view is normalized and the wavenumber is 16,000 cm −1 . In order to facilitate quantitative analysis of the calibration results, we calculate ϕ 2 (σ) from exp[iϕ 2 (σ)] by using a phase unwrapping algorithm, though it is not needed in the calibration and reconstruction. The maximum residual error of ϕ 2 (σ) is less than 1.0 × 10 −3 rad. These calibration results indicate that the self-correction method is applicative for the fieldable CDISP. Since ε 1 , ε 2 , and ϕ 2 (σ) do not change with the field of view, it is feasible to calibrate them only by using the data at central FOV. The reference beam at central FOV is shined on the system along the optical axis at the center of the aperture, therefore, the effect of polarization properties on calibration results can be ignored. The validity of the self-correction method for alignment errors and retardation calibration has been demonstrated by experiment [17]. We first simulate the process of radiometric calibration of the fieldable CDISP. The detector uses SARNOFF's CAM512 CCD as an example, whose area array is 512 × 512, the pixel size is 18 18 , μ μ × m m the frame frequency is 100 fps, and the quantum efficiency curve is shown in Figure 3.
According to the imaging chain analysis of the radiance from optical system entrance pupil to detector, the DN recorded by detector can be obtained by simulation. Using the radiance of the optical system entrance pupil and the simulated DN of detector, the dark current response ( )   it is feasible to calibrate them only by using the data at central FOV. The reference beam at central FOV is shined on the system along the optical axis at the center of the aperture, therefore, the effect of polarization properties on calibration results can be ignored. The validity of the self-correction method for alignment errors and retardation calibration has been demonstrated by experiment [17]. According to the calibration results of the alignment errors ε 1 , ε 2 and retardation ϕ 2 (σ) of the PSIM module, we could decouple the radiometric calibration coefficient A k (σ, θ) and diattenuation of fore-optics D f ore (σ, θ) from modulated radiometric calibration coefficient A k (σ, θ). Figure 5 gives the magnitude of the autocorrelation function of A k (σ, θ), which is at the edge FOV, the modulation channels of other FOV are consistent with this, only the amplitude is different. The figure only shows the values of channels 0 and 1, and the values of other channels are too small. Fortunately, only channels 0 and 1 are needed in the decoupling of A k (σ, θ) and D f ore (σ, θ). The input values and calibration results of diattenuation of the fore-optics at different field of view are shown in Figure 6. And the maximum relative error of the calibration results of D f ore (σ, θ) is less than 0.07% by using the presented method, while its maximum relative error is 63.57% by using the traditional calibration method. Figures 7 and 8 give the input values and calibration results of the retardation ϕ 1 (σ, θ) and retardance of the fore-optics δ f ore (σ, θ) at different fields of view. In order to facilitate quantitative analysis of the calibration results, we calculate ϕ 1 (σ, θ) from exp −i[ϕ 1 (σ, θ) + ϕ 2 (σ)] by using the calibration results ϕ 2 (σ) and a phase unwrapping algorithm. The calibration results of ϕ 1 (σ, θ) and δ f ore (σ, θ) of two methods are overlapped well with the input values, and the maximum residual errors are less than 1.0 × 10 −4 rad. It can be seen that, compared with the traditional calibration method, the presented method can achieve more accurate calibration when calibrating diattenuation of the fore-optics. For the calibration of alignment errors, phase retardances in the PSIM module, and retardance of the fore-optics, the two methods can both get good calibration results.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 12 of 17 are overlapped well with the input values, and the maximum residual errors are less than 1.0 × 10 −4 rad. It can be seen that, compared with the traditional calibration method, the presented method can achieve more accurate calibration when calibrating diattenuation of the fore-optics. For the calibration of alignment errors, phase retardances in the PSIM module, and retardance of the foreoptics, the two methods can both get good calibration results.   are overlapped well with the input values, and the maximum residual errors are less than 1.0 × 10 −4 rad. It can be seen that, compared with the traditional calibration method, the presented method can achieve more accurate calibration when calibrating diattenuation of the fore-optics. For the calibration of alignment errors, phase retardances in the PSIM module, and retardance of the foreoptics, the two methods can both get good calibration results.    After calibrating the parameters in the polarization radiometric calibration model, we substitute these calibration results into Equations (12)- (15) to reconstruct the Stokes parameters of the target light. In the simulation, we employ a 30° linearly polarized beam as the target light, whose theoretical values of normalized Stokes vectors are S1/S0 = 0.5, S2/S0 = 0.866, and S3/S0 = 0, respectively. Figure 9ac give the normalized Stokes vectors reconstructed by two methods, wherein the solid and broken line respectively represent the presented calibration method and the traditional calibration method, and the different colors represent reconstruction results at different field of view. Through contrastive analysis, the deviations of the Stokes parameters reconstruction results obtained by the traditional calibration method from the input values are large, and they are different with the change of field of view. While the reconstruction results obtained by the presented method are closer to the input values and overlapped well with different field of view. The reason is analyzed mainly because the radiometric calibration of the traditional calibration method does not take into account the influence of the polarization properties of optical system, and the spectral radiance obtained by the radiometric calibration deviates from that at the entrance pupil, which lead to the great inaccuracy of the reconstruction results of Stokes parameters. In order to further quantitatively analyze the measurement accuracy of the polarization information of the fieldable CDISP, we use the degree of  After calibrating the parameters in the polarization radiometric calibration model, we substitute these calibration results into Equations (12)- (15) to reconstruct the Stokes parameters of the target light. In the simulation, we employ a 30° linearly polarized beam as the target light, whose theoretical values of normalized Stokes vectors are S1/S0 = 0.5, S2/S0 = 0.866, and S3/S0 = 0, respectively. Figure 9ac give the normalized Stokes vectors reconstructed by two methods, wherein the solid and broken line respectively represent the presented calibration method and the traditional calibration method, and the different colors represent reconstruction results at different field of view. Through contrastive analysis, the deviations of the Stokes parameters reconstruction results obtained by the traditional calibration method from the input values are large, and they are different with the change of field of view. While the reconstruction results obtained by the presented method are closer to the input values and overlapped well with different field of view. The reason is analyzed mainly because the radiometric calibration of the traditional calibration method does not take into account the influence of the polarization properties of optical system, and the spectral radiance obtained by the radiometric calibration deviates from that at the entrance pupil, which lead to the great inaccuracy of the reconstruction results of Stokes parameters. In order to further quantitatively analyze the measurement accuracy of the polarization information of the fieldable CDISP, we use the degree of After calibrating the parameters in the polarization radiometric calibration model, we substitute these calibration results into Equations (12)- (15) to reconstruct the Stokes parameters of the target light. In the simulation, we employ a 30 • linearly polarized beam as the target light, whose theoretical values of normalized Stokes vectors are S 1 /S 0 = 0.5, S 2 /S 0 = 0.866, and S 3 /S 0 = 0, respectively. Figure 9a-c give the normalized Stokes vectors reconstructed by two methods, wherein the solid and broken line respectively represent the presented calibration method and the traditional calibration method, and the different colors represent reconstruction results at different field of view. Through contrastive analysis, the deviations of the Stokes parameters reconstruction results obtained by the traditional calibration method from the input values are large, and they are different with the change of field of view. While the reconstruction results obtained by the presented method are closer to the input values and overlapped well with different field of view. The reason is analyzed mainly because the radiometric calibration of the traditional calibration method does not take into account the influence of the polarization properties of optical system, and the spectral radiance obtained by the radiometric calibration deviates from that at the entrance pupil, which lead to the great inaccuracy of the reconstruction results of Stokes parameters. In order to further quantitatively analyze the measurement accuracy of the polarization information of the fieldable CDISP, we use the degree of polarization (DOP) [24], given by DOP(σ) = S 1 2 (σ) + S 2 2 (σ) + S 3 2 (σ)/S 0 (σ), to further analyze the reconstruction results, which are shown in Figure 9d.  Figure 9d. With these two calibration methods, the measurement accuracy of normalized Stokes vectors and DOP obtained from central to edge FOV are shown in Figure 10. It can be seen that the measurement accuracy of polarization information obtained by the traditional calibration method is decreased as the FOV increases. It is noteworthy that the measurement accuracy is about the same at different field of view by using the presented calibration method, and they are reduced by approximately an order of magnitude compared with the traditional calibration method. Though the measurement accuracy of S3/S0 by the presented method is less accurate than the traditional method, they are on an order of magnitude and are all very precise. The influence of the polarization properties of an optical system on the radiometric calibration of the fieldable CDISP is effectively compensated. With these two calibration methods, the measurement accuracy of normalized Stokes vectors and DOP obtained from central to edge FOV are shown in Figure 10. It can be seen that the measurement accuracy of polarization information obtained by the traditional calibration method is decreased as the FOV increases. It is noteworthy that the measurement accuracy is about the same at different field of view by using the presented calibration method, and they are reduced by approximately an order of magnitude compared with the traditional calibration method. Though the measurement accuracy of S 3 /S 0 by the presented method is less accurate than the traditional method, they are on an order of magnitude and are all very precise. The influence of the polarization properties of an optical system on the radiometric calibration of the fieldable CDISP is effectively compensated. Since the CDISP studied in this paper is mainly used for airborne remote sensing, the circular polarization component of the target light is usually small enough to be ignored in practical applications, so we use the linearly polarized beam as the target light to verify the validity of the presented calibration method in simulation experiment.

Conclusions
The difficulty of the calibration for the fieldable CDISP lies in the decoupling of radiometric calibration coefficient and diattenuation of fore-optics. In this paper, we propose a polarization radiometric calibration model to decouple them. We first theoretically derive the calibration and compensation model considering the alignment errors of PSIM, the variation of retardations, and the polarization properties at different field of view of optical system. Before decoupling the radiometric calibration coefficient and the polarization properties, we perform classical radiometric calibration and self-correction calibration for the system. Finally, we use 0° and 45° linearly polarized beams to calibrate other parameters. By using the decoupled radiometric calibration coefficient and abovementioned calibration results, we can accurately obtain the spectral radiance and reconstruct the polarization information of detected object. The validity and feasibility of the presented method have been confirmed by simulation experiments, and the reconstruction accuracy of the polarization information is improved compared with the traditional calibration method. The presented calibration method is significant for the quantitative application of the CDISP.
For the integral calibration process for the practical application of CDISP, the spectral calibration is not discussed in this paper. This is because the spectral calibration of CDISP is consistent with the calibration method of traditional dispersive imaging spectrometer, which can be completed in front of the polarization radiometric calibration. Since the CDISP studied in this paper is mainly used for airborne remote sensing, the circular polarization component of the target light is usually small enough to be ignored in practical applications, so we use the linearly polarized beam as the target light to verify the validity of the presented calibration method in simulation experiment.

Conclusions
The difficulty of the calibration for the fieldable CDISP lies in the decoupling of radiometric calibration coefficient and diattenuation of fore-optics. In this paper, we propose a polarization radiometric calibration model to decouple them. We first theoretically derive the calibration and compensation model considering the alignment errors of PSIM, the variation of retardations, and the polarization properties at different field of view of optical system. Before decoupling the radiometric calibration coefficient and the polarization properties, we perform classical radiometric calibration and self-correction calibration for the system. Finally, we use 0 • and 45 • linearly polarized beams to calibrate other parameters. By using the decoupled radiometric calibration coefficient and above-mentioned calibration results, we can accurately obtain the spectral radiance and reconstruct the polarization information of detected object. The validity and feasibility of the presented method have been confirmed by simulation experiments, and the reconstruction accuracy of the polarization information is improved compared with the traditional calibration method. The presented calibration method is significant for the quantitative application of the CDISP.
For the integral calibration process for the practical application of CDISP, the spectral calibration is not discussed in this paper. This is because the spectral calibration of CDISP is consistent with the