The Analysis and Experiment of Pixel-Matching Method for Space-Dimensional Dual-Coded Spectropolarimeter

: In order to meet the high accuracy pixel-matching requirements of space-dimensional dual-coded spectropolarimeter, a dual-coded image pixel-matching method based on dispersion modulation is proposed. The mathematics of the dispersion power and the pixel matching is modeled. The relationship between different pixel-matching coefﬁcients and the peak signal-to-noise ratio (PSNR) and structure similarity index measure (SSIM) of reconstructed images is analyzed. An imaging system experiment consisting of a digital micromirror device (DMD), a micro-polarizer array detector (MPA), and a prism–grating–prism (PGP) is built to reconstruct a spectral linear polarization data cube with 50 spectral channels and linear polarization parameters. The contrast ratio of the reconstructed spectropolarimeter image was raised 68 times against the ground truth. It can be seen from the reconstruction evaluation analysis that the spectral data and polarization data can be matched effectively by optimizing the dispersion coefﬁcient of the PGP. The system can effectively reconstruct when the noise SNR is greater than 15 dB. The PSNR and SSIM of the reconstruction images can be improved by increasing the pixel-matching spacing. The optimal choice of the dual-coded pixel-matching spacing is one super-polarized pixel. The spectral resolution and quality of the spectropolarimeter are improved using the pixel-matching method.


Introduction
Spectral characteristics are inherent characteristics of matter. The intensity and spectral 2D data of an object can be obtained by spectral imaging, and it is sensitive to the type, material, and composition of the target. It can effectively distinguish object material and has applications in remote sensing and medicine [1][2][3][4][5]. Spectral polarization imaging (SPI) can acquire space, spectrum, and polarization data and has been a research hotspot worldwide.
Conventional spectral polarization imaging modes use a combination of polarization devices and spectrometers. The traditional time-sharing imaging method uses a scanning form, which makes it unable to image moving targets [6][7][8]. Compressive sensing is used to obtain spatial and spectral data in a scene, which is called compressive imaging [9,10]. This acquisition method reduces the cost of sensing, storing, and transmitting data cubes. It can be applied to detecting moving objects and effectively captures 3D information using 2D projection measurement. Typical compressive spectral polarization imaging methods include the division of amplitude [11,12], Fourier [13], the division of focal plane [14][15][16][17], and a computer tomography imaging spectrometer (CTIS) [18].
With the development of technology, more and more scholars are using compressive imaging methods to achieve spectral polarization imaging. Tsung-Han Tsai proposed a Photonics 2023, 10 coded aperture snapshot spectral polarization imager (CASSPI) [19]. CASSPI utilizes two additional birefringent crystals behind the encoded aperture to achieve polarization coding, enabling the acquisition of the first two Stokes parameters. Chen et al. proposed a compressive spectrum and polarization imager for linear Stokes spectral image measurements. A compressive-based coded aperture snapshot linear Stokes imaging spectral polarimeter (CASSIS) was proposed, which could recover linear Stokes spectral images from sparse samples [20]. Xu proposed a pixelated polarizers array-based coded aperture snapshot spectral polarization imaging method. The 2D measurement was used to recreate the polarization images of 16 spectral bands [21]. Ning used a combination of a double Amici prism (DAP) and Wollaston prisms for a compressive circular polarization snapshot spectral imager (CCPSSI) [22]. The system used Wollaston prisms to obtain the target dual Stocks polarization information containing 25 spectral channels with circular polarization. However, snapshot spectral polarization imaging technology that can provide multi-dimensional imaging of moving objects still faces many challenges, including stable imaging mechanisms, accurate optimization methods, and efficient reconstruction algorithms. Previously, we proposed a compressive space-dimension dual-coded hyperspectral polarimeter (CSDHP). Hyperspectral polarization data with a single shot were realized [23]. In the compressive spectral polarization imaging method, the coded aperture is used to encode spatial and spectral data, and the micro-polarizer array (MPA) detector is used to encode the polarization data. The pixels of digital micromirror devices (DMDs) and MPA detectors are determined using a pixel-matching method to achieve data unmixing and reconstruction. The spectropolarimeters are referred to as dual-coded methods.
In spatial imaging, pixel matching can be achieved by modeling spectral data and spatial data based on the dispersion power of dispersion components. During the polarization encoding process of MPA detectors, the matching relationship between polarization data and spatial data is known. However, the spatial coordinate axes of the two are different, which makes it impossible to accurately know the matching relationship between spatial data, spectral data, and polarization data. As a result, high-precision pixel matching cannot be achieved. In order to improve the accuracy of pixel matching and enhance the imaging ability of space-dimensional dual-coded systems, we propose a dual-coded pixel-matching method based on the dispersion modulation of the prism-grating-prism (PGP) and model the relationship between the angular dispersion of the PGP and pixel matching. We analyze and evaluate the pixel-matching results of the images through the peak signal-to-noise ratio (PSNR) and structure similarity index measure (SSIM) of the reconstructed results. Finally, experiments are established to verify the pixel-matching accuracy using 50 spectral channels and linear polarization parameters obtained.
In Section 2, we induce a method of pixel matching and unmixing for the spacedimensional dual-coded spectropolarimeter (SDSP). A design process for optimizing the adjacent spectra spacing is proposed. The relationship between the angular dispersion of the PGP and the focal length of the imaging lens with different adjacent spectra spacing is analyzed. In Section 3, we set up tabletop experiments. The spectral polarization images with different spectral spacings are acquired. It shows the advantages of spectral polarization imaging over spectral imaging through comparative experiments. We compare the reconstruction effects of different adjacent spectra spacing. The pixel-matching effect of different adjacent spectral spacing is evaluated through SSIM and PSNR to determine the optimal solution. In Sections 4 and 5, we present a discussion and conclusion.

Imaging Configuration
Previously, we proposed a compressive space-dimensional dual-coded hyperspectral polarization imaging method [23]. The structure of SDSP is shown in Figure 1. First, the data cube is imaged at the first image plane through the objective lens, and the 3D data cube is spatially modulated by the DMD. The DMD consists of a series of micromirrors, each of which has two possible states: "ON" (rays reaching this micromirror can be transmitted through the relay lens) and "OFF" (light reaching this micromirror will be reflected out of the optical system). Applying a coding mode on the DMD through computer control. The micromirror tilt angle is ±12 • . The coded data cube reaches the PGP through the relay lens and reaches the detector through the PGP and imaging lens. Then, we acquire 2D spatial, spectral, and polarization aliasing data on the detector. The MPA detector is an array of micro-polarizers integrated directly into the focal plane of the camera sensor, and 2 × 2 adjacent pixels correspond to a polarization super-pixel.

Imaging Configuration
Previously, we proposed a compressive space-dimensional dual-coded hyperspectral polarization imaging method [23]. The structure of SDSP is shown in Figure 1. First, the data cube is imaged at the first image plane through the objective lens, and the 3D data cube is spatially modulated by the DMD. The DMD consists of a series of micromirrors, each of which has two possible states: "ON" (rays reaching this micromirror can be transmitted through the relay lens) and "OFF" (light reaching this micromirror will be reflected out of the optical system). Applying a coding mode on the DMD through computer control. The micromirror tilt angle is ±12°. The coded data cube reaches the PGP through the relay lens and reaches the detector through the PGP and imaging lens. Then, we acquire 2D spatial, spectral, and polarization aliasing data on the detector. The MPA detector is an array of micro-polarizers integrated directly into the focal plane of the camera sensor, and 2 × 2 adjacent pixels correspond to a polarization super-pixel. In the CDSP, a DMD is used as a spatial light modulator, and the data cube is encoded at the primary image plane using the DMD to achieve data compression encoding. The micro-polarizer array detector completes the encoding and information acquisition of the polarization channel simultaneously. DMD is combined with an MPA to achieve controlled mixing of spectrally polarized data cubes, so the data cubes undergo two encodings. The first encoding is the spatial data coded by the coded aperture at the first image plane. Because of the special structure of MPA, its polarization encoding corresponds to the spatial position one by one. So, the second encoding is the polarization data (spatial) versus spectrum data by the MPA. The micromirrors of the DMD and the pixels of the MPA achieve precise encoding of spectral polarization data by means of pixel matching. The interval between adjacent spectral lines is determined by the dispersion power of the PGP and the focal length of the imaging lens. The simulation model of the PGP and the imaging lens is shown in the figure. The material of the two prisms in the PGP is BK7 with top angles of 7.4°. The material of the grating is B270, the number of lines is 360 lins/um, and the dispersion level is +1. Each pixel is coded according to the spatial coordinate position. The four adjacent pixels within a square region form a 2 × 2 super-polarized pixel. The polarized focal plane array divides into several 2 × 2 units to calculate the Stokes vector based on the super-polarized pixels.
The Stokes vector method is represented as In the CDSP, a DMD is used as a spatial light modulator, and the data cube is encoded at the primary image plane using the DMD to achieve data compression encoding. The micro-polarizer array detector completes the encoding and information acquisition of the polarization channel simultaneously. DMD is combined with an MPA to achieve controlled mixing of spectrally polarized data cubes, so the data cubes undergo two encodings. The first encoding is the spatial data coded by the coded aperture at the first image plane. Because of the special structure of MPA, its polarization encoding corresponds to the spatial position one by one. So, the second encoding is the polarization data (spatial) versus spectrum data by the MPA. The micromirrors of the DMD and the pixels of the MPA achieve precise encoding of spectral polarization data by means of pixel matching. The interval between adjacent spectral lines is determined by the dispersion power of the PGP and the focal length of the imaging lens. The simulation model of the PGP and the imaging lens is shown in the figure. The material of the two prisms in the PGP is BK7 with top angles of 7.4 • . The material of the grating is B270, the number of lines is 360 lins/um, and the dispersion level is +1. Each pixel is coded according to the spatial coordinate position. The four adjacent pixels within a square region form a 2 × 2 super-polarized pixel. The polarized focal plane array divides into several 2 × 2 units to calculate the Stokes vector based on the super-polarized pixels.
The Stokes vector method is represented as where I is the total irradiance of the beam, M is the horizontal (I 0 ) polarized flux component minus the vertical (I 90 )  The degree of linear polarization (DOLP), which represents the percentage of the total intensity of polarized rays, is defined as The aliasing data received by the MPA detector, coming from each of the four pixels (0 • , 45 • , 90 • , and 135 • ) composed of one super-pixel, are relatively independent and in different directions. In compressive imaging, the TwIST algorithm is selected to compare structural similarity (SSIM) and peak signal-to-noise ratio (PSNR). The TwIST algorithm used in this study is represented as [25] where τ is a weighting coefficient between fidelity and sparsity and Γ tw ( f ) is the total regularization term. If we define the discrete version of the original continuous 3D data cube f 0 (x, y; λ) as the f 0 (i, j; w), the Total Variation (TV) regularization term becomes The spectral images of the four polarization angles are obtained separately after reconstruction. The pixels are then reconstructed into a spectral polarization image with an instantaneous field-of-view error based on the corresponding positions of the polarization code. Subsequently, using the bilinear interpolation method, the polarization super-pixels are complemented to resolve the instantaneous field-of-view error. Finally, the spectral data cubes of the Stocks are obtained separately. The 4D data cube is deconstructed with a single shot.

Matching Model
The schematic diagram of pixel matching, unmixing, and reconstruction is shown in Figure 2. The target data cube is coded by the DMD placed at the primary image plane of the objective lens. The rays pass through the relay lens to the PGP to generate dispersion. Eventually, the dispersive rays converge to the MPA detector by the imaging lens. The dispersion data cube is selected according to the polarization direction of the MPA at a specific position. A matched dispersive 2D data cube is formed. Subsequently, according to the coded polarization pixels, the super-pixels are separated. The groups of pixels with the same polarization angle are grouped to form four groups of images with different polarization directions. Finally, the matched dispersive 2D data cube is reconstructed. criterion. Therefore, according to the Rayleigh criterion, the linear dispersion power is expressed as where dl d λ denotes the line linear dispersion power, MPA p denotes the size of the MPA detector pixels, and the coefficient γ is introduced to represent a multiple of detector pixels and can only be a positive even number. Convert linear dispersion into character dispersion: According to the design process, the center of the Airy point of the λ n spectrum coincides with the primary dark ring of the Airy point of the λ n+1 spectrum at the pixel. To ensure the aliasing state, the adjacent spectral interval is changed from 2γ(γ ∈ N * ) to γ pixels. At this time, the sensing matrices are still known and conform to the original RIP criterion. Therefore, according to the Rayleigh criterion, the linear dispersion power is expressed as where dl dλ denotes the line linear dispersion power, P MPA denotes the size of the MPA detector pixels, and the coefficient γ is introduced to represent a multiple of detector pixels and can only be a positive even number.
Convert linear dispersion into character dispersion: where f is the focal length of the imaging lens.
Because the role dispersion of the PGP is the sum of the dispersion rates of each device, there is Combine the Equations (6) and (7), From Equation (8), it can be seen that controlling the focal length f can achieve changes γ without changing the overall dispersion of the PGP.
To realize spectral polarization data unmixing and dual-coded pixel matching for the MPA, DMD, and PGP, according to the polarization encoding form of the MPA detecto γ, determines the resolution and imaging quality of the system. So, we adopt a way to optimize the process of adjacent spectral spacing in terms of spectral resolution, spatial resolution, and imaging quality. The optimization process is shown in Figure 3.  Firstly, the DMD, MPA, and PGP introduce spatial coordinate aliasing, which is solved by pixel matching and super-pixel unmixing. Secondly, the pixel-matching theory is proposed based on to achieve controllable mixing of spatial, spectral, and  Firstly, the DMD, MPA, and PGP introduce spatial coordinate aliasing, which is solved by pixel matching and super-pixel unmixing. Secondly, the pixel-matching theory is proposed based on dl dλ = 2γ × P MPA to achieve controllable mixing of spatial, spectral, and polarization data. Thirdly, pixel matching is used to change the γ to 2, 4, 6, 8, . . . . Changing the value to 2, 4, 6, 8, . . . , the aliased data is unmixed and reconstructed. Finally, We analyze the PSNR and SSIM of the reconstructed image under different spectral resolutions and spatial resolutions and then judge the value of γ.
The pixel size of the MPA corresponds to the pixel size of the DMD through scaling of the focal length design of the relay and imaging lens. According to Equations (5)-(8), the distance of line dispersion is shown in Figure 4 (only γ = 2 and γ = 4 are drawn).
ing the value to 2, 4, 6, 8, …, the aliased data is unmixed and reconstructed. Finally, We analyze the PSNR and SSIM of the reconstructed image under different spectral resolutions and spatial resolutions and then judge the value of γ .
The pixel size of the MPA corresponds to the pixel size of the DMD through scaling of the focal length design of the relay and imaging lens. According to Equations (5)-(8), the distance of line dispersion is shown in Figure 4 (only γ = 2 and γ = 4 are drawn). Therefore, we model the discrete energy transport of the system with the spectral density of each polarization Stocks parameter before entering the system, expressed as   Therefore, we model the discrete energy transport of the system with the spectral density of each polarization Stocks parameter before entering the system, expressed as where S 0 (x 0 , y 0 ; λ), S 1 (x 0 , y 0 ; λ), and S 2 (x 0 , y 0 ; λ) represent the corresponding spatial and spectral data cubes of the three covariates of Stocks before entering the system. The spectral density after the coded aperture is where T(x DMD , y DMD ) denotes the spectral polarization density corresponding to the DMD spatial coordinates. After optical design alignment, the spatial coordinates of space, the DMD, and the detector target surface correspond to each other, and so they are where t(i, j) represents a binary value at the (i, j) th micromirror on the DMD and rect x DMD p DMD − i, y DMD p DMD − j denotes the spatial range of the (i, j) th micromirror. After PGP dispersion along the y-axis, the spectral density currently is where the δ function represents the dispersion effect of PGP and λ C is the central wavelength.
The spectral density is dispersed on the y-axis of the detector plane. The encoded data cube is projected onto the sensor of the MPA detector. The final continuous image on the detector plane can be described as where ω(λ) is the spectral response coefficient corresponding to the wavelength. Finally, each pixel of the detector measures the integrated intensity of the spectral density at a specific polarization angle. Detector surfaces are pixelated in space at the size of image elements p MPA such that that the spatial domain at the detector plane g (x 3 , y 3 ) is sampled as (14), the correspondence of the micromirrors and pixels can be derived. Each pixel contains discrete spectral data and three parameters of the polarization data. In the polarization array, each polarization parameter is encoded and filtered to form the spectral polarization aliasing data. Without changing the system structure, the pixelmatching spacing is determined by the angular dispersion of the PGP and the focal length of the imaging lens; the normalized results are shown in Figure 5. From the figure, It can be seen that the angular dispersion is proportional to γ. In addition, the correspondence can be used to provide a basis for the design of PGP and the imaging lens.  Figure 6 depicts the experimental prototype of SDSP. The coded aperture selected DMD (Texas Instruments DLP6500) included 1920 × 1080 elements of a random binary pattern with a 7.65 × 7.65µm mirror pitch. The measured images were taken using an MPA detector (Flir Blackfly BFS-U3-123S6C-C) with a pixel size of 3.45 µm and a resolution of 4120 × 3000. The working wavelength of the prism was 425-675 nm (Edmund #35-788). The central wavelength of the grating was 532 nm (WP-360/532-25.4).

Contrast of Reconstruction Results
According to Equations (5)- (8) and Figure 4 in Section 2, we changed the focal length of the imaging lens to make γ = 2 and 4 to unmix and reconstruct the data cube. As shown in Figure 7a, the abbreviated letters of this institution were chosen for the corresponding target: "C" for red aluminum, "U" for black plastic, "S" for light yellow cardboard, and "T" for blue rubber. The ground truth is shown in Figure 7b. Because the letter "U" was black, it is not obvious in a dark room.  Figure 6 depicts the experimental prototype of SDSP. The coded aperture selected DMD (Texas Instruments DLP6500) included 1920 × 1080 elements of a random binary pattern with a 7.65 × 7.65µm mirror pitch. The measured images were taken using an MPA detector (Flir Blackfly BFS-U3-123S6C-C) with a pixel size of 3.45 µm and a resolution of 4120 × 3000. The working wavelength of the prism was 425-675 nm (Edmund #35-788). The central wavelength of the grating was 532 nm (WP-360/532-25.4).

Contrast of Reconstruction Results
According to Equations (5)- (8) and Figure 4 in Section 2, we changed the focal length of the imaging lens to make γ = 2 and 4 to unmix and reconstruct the data cube. As shown in Figure 7a, the abbreviated letters of this institution were chosen for the corresponding target: "C" for red aluminum, "U" for black plastic, "S" for light yellow cardboard, and "T" for blue rubber. The ground truth is shown in Figure 7b. Because the letter "U" was black, it is not obvious in a dark room.  The reconstructed spectral polarization image is shown in Figure 8a-e. Figure 8f,g show the magnified view at 539 nm, respectively. The variable γ = 4, spatial resolution is 512 × 496, and the 50 channels are from 450 to 650 nm.
A comparison of Figure 8a-d shows that each letter had high spectral resolutions in the spectral images of all four polarization directions. The letters "C" and "T" had distinct peak bands that could be easily distinguished. The letter "S" had a complete spectrum. Since the letter "U" was black and had the same spectral characteristics as the background, it was difficult to distinguish. However, in the polarization images of Figure 8e,g, the letter "U" was clearly distinguishable because the polarization characteristics were different from the background. Comparing the contrast ratio (CR) of the letter "U" in Figure 7b with that in Figure 8g, the contrast ratio improved from 0.57 to 39.55. It can be concluded that polarization can improve the imaging contrast and the system spectral polarization imaging effect is effective in the case of interference due to unfavorable environmental spectral imaging. This proves that spectrally polarized images can be used to distinguish the object material.
Next, we changed the structure of the PGP and the imaging lens to make γ = 2 and performed the unmixing and reconstruction experiment again.
It is known from Figure 9 that each letter had the same high spectral resolution in all four polarization directions of the spectral image. The letters "C", "U", "S", and "T" did not differ much from those in Figure 8, but the signal-to-noise ratio (SNR) was lower. The reconstructed spectral polarization image is shown in Figure 8a-e. Figure 8f,g show the magnified view at 539 nm, respectively. The variable 4 γ = , spatial resolution is 512 × 496, and the 50 channels are from 450 to 650 nm.

Evaluation
To compare the final effects of the two aliasing methods, the reconstruction performance was measured using SSIM and PSNR, thus evaluating the effect of dual-coded aliasing of spectral polarization data cubes.
To analyze the robustness of noise, Gaussian noise was introduced into the standard image to generate contaminated images at 5 dB intervals within an SNR of 5 to 45 dB. The PSNR and SSIM for each polarization angle with γ = 2 and 4 are shown in Figures 10-13. Each graph has three horizontal and vertical auxiliary lines, and the trend is shown above and to the right, where the wavelengths selected were 488, 539, and 632 nm, and SNR was selected as 15, 25, and 35 dB. that polarization can improve the imaging contrast and the system spectral polarization imaging effect is effective in the case of interference due to unfavorable environmental spectral imaging. This proves that spectrally polarized images can be used to distinguish the object material.
Next, we changed the structure of the PGP and the imaging lens to make 2 γ = and performed the unmixing and reconstruction experiment again. It is known from Figure 9 that each letter had the same high spectral resolution in all four polarization directions of the spectral image. The letters "C", "U", "S", and "T" did not differ much from those in Figure 8, but the signal-to-noise ratio (SNR) was lower. DOLP, (f) 135° @ 539 nm, and (g) DOLP @ 539 nm.

Evaluation
To compare the final effects of the two aliasing methods, the reconstruction performance was measured using SSIM and PSNR, thus evaluating the effect of dual-coded aliasing of spectral polarization data cubes.
To analyze the robustness of noise, Gaussian noise was introduced into the standard image to generate contaminated images at 5 dB intervals within an SNR of 5 to 45 dB. The PSNR and SSIM for each polarization angle with γ = 2 and 4 are shown in Figures 10-13.
Each graph has three horizontal and vertical auxiliary lines, and the trend is shown above and to the right, where the wavelengths selected were 488, 539, and 632 nm, and SNR was selected as 15, 25, and 35 dB.      As can be seen from the overall perspective, when the SNR is greater than 25 dB, the PSNR is greater than 20 dB whether =2 γ or =4 γ . The reconstruction performance of the shorter spectral band is lower than that of the longer spectral band. The best reconstruction performance can be achieved with an SNR greater than 30 dB. When the SNR of the contaminated image is greater than 35 dB, the SSIM value is greater than 0.9 dB and the SSIM value is greater than 0.8 dB. When the image SNR of the four reconstructed polarization reaches 45 dB, their SSIMs are all close to 1 dB. When the SNR is greater than 15 dB, all the images can be effectively reconstructed. From the PSNR corresponding to the three wavelengths, there is a slight decrease in the PSNR of =2 γ over =4 γ , but the slope is almost unchanged. It shows that the change in γ only affects the degraded image and does not affect the reconstruction algorithm accuracy. From the PSNR corresponding to the three SNRs, there is a large improvement from 15 to 25 dB, but the improvement becomes slow from 25 to 35 dB, which indicates that the image reconstruction quality is nonlinearly related to the image SNR after the SNR reaches a certain level. From the SSIM corresponding to the three wavelengths, curve slopes vary with different spacing and wavelengths, which indicates that the adjacent wavelengths have completely different aberrations at γ = 2 or γ = 4.

Discussion
For spectral polarization imaging, we propose a pixel-matching imaging structure based on DMD, PGP, and MPA detectors, building on our previous research. Based on the corresponding hardware principles of obtaining spectral data and polarization data, we introduce a new classification mode named dual-coded imaging mode. The SDSP retains the advantages of compressive imaging and uses unmixing for polarization data, As can be seen from the overall perspective, when the SNR is greater than 25 dB, the PSNR is greater than 20 dB whether γ = 2 or γ= 4. The reconstruction performance of the shorter spectral band is lower than that of the longer spectral band. The best reconstruction performance can be achieved with an SNR greater than 30 dB. When the SNR of the contaminated image is greater than 35 dB, the SSIM value is greater than 0.9 dB and the SSIM value is greater than 0.8 dB. When the image SNR of the four reconstructed polarization reaches 45 dB, their SSIMs are all close to 1 dB. When the SNR is greater than 15 dB, all the images can be effectively reconstructed.
From the PSNR corresponding to the three wavelengths, there is a slight decrease in the PSNR of γ = 2 over γ= 4, but the slope is almost unchanged. It shows that the change in γ only affects the degraded image and does not affect the reconstruction algorithm accuracy. From the PSNR corresponding to the three SNRs, there is a large improvement from 15 to 25 dB, but the improvement becomes slow from 25 to 35 dB, which indicates that the image reconstruction quality is nonlinearly related to the image SNR after the SNR reaches a certain level. From the SSIM corresponding to the three wavelengths, curve slopes vary with different spacing and wavelengths, which indicates that the adjacent wavelengths have completely different aberrations at γ = 2 or γ = 4.

Discussion
For spectral polarization imaging, we propose a pixel-matching imaging structure based on DMD, PGP, and MPA detectors, building on our previous research. Based on the corresponding hardware principles of obtaining spectral data and polarization data, we introduce a new classification mode named dual-coded imaging mode. The SDSP retains the advantages of compressive imaging and uses unmixing for polarization data, eliminating the process of reconstruction using algorithms such as Fourier transform and obtaining high spatial resolution and spectral resolution.
In the experiment, data cubes with two groups of 50 spectral channels and three Stocks vectors with different linear dispersion powers were reconstructed. When the spectral curve of the object is similar to the background, the contrast of the spectral polarization reconstruction image is increased from 0.57 to 39.55, which is 68 times higher than the ground truth. The data provided at this time were selected through multiple experiments.
The overall variation of PSNR and SSIM shows that the image reconstruction quality at this time at γ = 2 is not as good as at γ = 4. With the rise in γ, the system aberration decreases, the degraded image is closer to the ideal image, the adjacent spectra spacing channels become larger, and the focal length of the imaging lens becomes longer; thus, the reconstruction quality increases. This also indicates that the spectral smile is the main problem of system reconstruction. Therefore, a larger adjacent spectral spacing should be selected under the premise of ensuring spectral resolution and spatial resolution. Without changing the spatial resolution or structure of the imaging system, we finally conducted a comparative theoretical analysis of the results of different γ values, and the results are shown in Table 1. When the spatial resolution remains unchanged at 512 × 496, theoretical spectral channels correspond to different γ, PSNR, and SSIM, as shown in Table 1. It can be seen that after exceeding γ= 4, no significant increase in PSNR or SSIM occurs. This indicates that the system is sufficient in correcting linear dispersion power at this time.
At present, the resolution of DMD is the main limiting factor for the spatial resolution of our method, which is also a problem that everyone is facing. However, our method can increase the spatial resolution to the DMD limit resolution. In terms of spectral resolution, we can increase the number of channels to 1089 through simulation, but it takes a long time to calculate, so we will start to improve the reconstruction rate based on deep learning in the future.

Conclusions
In summary, we proposed an image pixel-matching approach based on DMD, PGP, and MPA detectors. A coefficient γ is proposed to model the dispersive power and pixel matching. The relationships among γ, the focal length of the imaging lens, and the dispersion power are given. The results indicate that the spectral imaging effect is not ideal when the target and background spectra are close, but the SDSP is still effective. The analysis shows that the system has good robustness to noise under different pixel-matching rules. The data can be effectively reconstructed when the SNR is greater than 15 dB under different matching rules.
The analysis of the imaging results of the matching rules shows that the larger the adjacent spectral spacing, the better the imaging results, and a loss of resolution occurs. From the perspective of PSNR, the changing of matching rules does not affect the accuracy of the reconstruction algorithm, and the image reconstruction quality is nonlinearly related to the image SNR after it reaches a certain level. From the perspective of SSIM, adjacent wavelengths at γ = 2 or γ= 4 have completely different aberrations, showing the spectral smile is the main consideration in the system reconstruction. Finally, the comparative analysis of the results in different γ values theoretically shows that γ= 4 is a relatively balanced result.