Snapshot Multi-Wavelength Birefringence Imaging

A snapshot multi-wavelength birefringence imaging measurement method was proposed in this study. The RGB-LEDs at wavelengths 463 nm, 533 nm, and 629 nm were illuminated with circularly polarized light after passing through a circular polarizer. The transmitted light through the birefringent sample was captured by a color polarization camera. A single imaging process captured light intensity in four polarization directions (0°, 45°, 90°, and 135°) for each of the three RGB spectral wavelength channels, and subsequently measured the first three elements of Stokes vectors (S0, S1, and S2) after the sample. The birefringence retardance and fast-axis azimuthal angle were determined simultaneously. An experimental setup was constructed, and polarization response matrices were calibrated for each spectral wavelength channel to ensure the accurate detection of Stokes vectors. A polymer true zero-order quarter-wave plate was employed to validate measurement accuracy and repeatability. Additionally, stress-induced birefringence in a PMMA arch-shaped workpiece was measured both before and after the application of force. Experimental results revealed that the repeatability of birefringence retardance and fast-axis azimuthal angle was better than 0.67 nm and 0.08°, respectively. This approach enables multispectral wavelength, high-speed, high-precision, and high-repeatability birefringence imaging measurements through a single imaging session.


Introduction
Birefringence commonly manifests in anisotropic materials, encompassing crystals, glass, plastics, thin films, and biological samples [1,2].Additionally, external influences such as force, heat, electric, and magnetic fields can induce anisotropic properties [3][4][5][6].The primary optical effect of birefringence entails the splitting of incident light into two orthogonal vibration directions along the principal axes.This separation leads to disparate propagation velocities for the ordinary and extraordinary light components, consequently causing a retardation effect.The retardance value and the direction of the fast axis comprehensively characterize birefringence.Birefringence measurements not only elucidate the anisotropic parameters of optical materials but also serve as the primary basis for evaluating external physical parameters.Crystal birefringence across a broad spectral dynamic range was employed for the precise calibration of polarization optics, such as waveplates and compensators [7].Moreover, birefringence served as the foundation for essential parameters including stress defects, mechanical strength, thermal stability, imaging quality, and beam transmission quality [8].In the realm of biological samples, birefringence played a pivotal role in biosensing and pathology assessment [9,10].Furthermore, the multi-wavelength birefringence exhibited by novel nanostructured materials such as meta-materials and meta-surfaces offers crucial insights for the development of innovative multifunctional devices [11][12][13][14].These applications increasingly necessitate highly accurate, real-time, and distributed measurements of birefringence across a broad spectral dynamic range.
Sensors 2024, 24, 5174 2 of 13 Numerous methods for measuring birefringence have been employed, with photoelasticity techniques utilizing polariscopes, Senarmont compensation, or Tardy phase shifting being particularly prevalent [15][16][17].These methods offer the advantage of a straightforward instrument structure and a sizable measurement aperture and have reached a state of mature commercialization.However, many photoelasticity techniques necessitate the mechanical rotation of either the polarizing element or the sample for concurrent measurement of birefringence retardance and direction.Such mechanical rotation not only diminishes measurement accuracy due to beam mechanical drift but also reduces measurement speed.More recently, liquid crystal variable retarders (LCVRs) have emerged as alternatives for birefringence measurements, heralding the development of novel photoelasticity techniques and demonstrating promise for multi-wavelength and wide dynamic range birefringence assessments [18,19].LCVRs themselves are birefringence dispersive.Their calibration for multi-wavelength applications proves intricate.Moreover, the liquid crystal molecule orientation within LCVRs is markedly influenced by ambient temperature variations, posing challenges to the long-term stability and reproducibility of this measurement approach.Overcoming these hurdles presents formidable tasks.
Modulators such as Faraday rotator, electro-optical modulator, and photoelastic modulator were utilized to facilitate the detection of optical signal modulation, thereby achieving birefringence measurements with a precision surpassing 1 nm and a time resolution in the millisecond range [20,21].However, polarization modulation methods are predominantly employed for single-point measurements.To obtain birefringence mappings, scanning of samples is necessary, with the dynamic range typically limited to half of the wavelength of the detection light source.Advancements in micro/nano processing technology have led to the utilization of micrometer-level wire-grid polarizers for processing imaging detector pixels.Sony Corporation, for instance, has developed an array polarization CMOS image sensor.This high-speed polarization image sensor has been effectively employed for birefringence distribution measurements, as well as rapid birefringence process assessments for fluid samples [22,23].Presently, the polarization imaging sensor is available in monochromatic as well as color configurations.However, this image sensor is only capable of linearly polarized light detection in four polarization directions (0 • , 45 • , 90 • , and 135 • ), and it does not detect circular polarization, thus constraining the birefringence retardance measurement range within the half-wavelength detection limit.Consequently, birefringence distinguishing between fast and slow axes are not achievable.
To address these limitations, this study employed an integrated approach utilizing multi-wavelength LED light arrays, in conjunction with circular polarizers to constitute a detection light source.Additionally, a color polarization camera was utilized for imaging measurements.This integrated setup was anticipated to enable snapshot multi-wavelength birefringence imaging measurements while simultaneously achieving high precision and fast measurements of birefringence retardance and fast-axis azimuthal angle across a wide dynamic range.Such advancements provide sophisticated means for the measurement of birefringence dispersion and properties of novel materials.

Methods
To address the issue of existing polarization image sensors' insensitivity to circular polarization, this study modified the detection light source to emit circularly polarized light.As depicted in Figure 1a, an integrated red, green, and blue multi-wavelength LED light source (RGB-LED) was transformed into circularly polarized light through a circular polarizer.This circularly polarized light was then directed onto the sample for snapshot imaging measurements using a polarization camera.As illustrated in Figure 1b, the camera utilized in this study was a color polarization camera model DYK 33GX250 from The Imaging Source Asia Co., Ltd.(Shenzhen), featuring a Sony IMX250MYR-C CMOS (color) image sensor.The camera boasted a resolution of 2448 × 2048 pixels (5 MP) with a pixel size of 3.45 µm × 3.45 µm and operated at a frame rate of 75 fps.The camera employed a telecentric lens (model FT03X110, Hikvision Digital Technology Co., Ltd.Hangzhou), providing a field of view for the object of 36.67 mm and a working distance of 110 mm.The camera's imaging system comprised a polarized superimage element formed by four wire-grid polarizer pixels oriented horizontally (0°), vertically (90°), and diagonally (45° and 135°).Each polarized superimage element was equipped with a mosaic of Bayer filters in a sequential arrangement, as depicted in the inset of Figure 1a.The RGB-LEDs utilized in the setup were sourced from KOMA Vision Technology Co., Ltd.(Dongguan), illuminating a surface area of 100 × 100 mm.The LEDs emit natural light.A circular polarizer was fabricated by laminating a linear film polarizer onto a high-contrast polymer achromatic quarter-wave retarder sheet, covering the wavelength range of 400-700 nm, sourced from Edmund Optics (Shenzhen) Co., Ltd.The ellipticity uniformity of the circularly polarized light passing through the polarizer exceeds 99%.To shield the light source between the LEDs and the sample, a fused silica window measuring 100 × 100 mm with a thickness of 18 mm was installed.This window protects the circular polarizer from scratches and contamination during sample measurement while providing a smooth, light-transmitting support platform for the sample.The spectral response characteristics of the color camera were taken into account when configuring the RGB three-color light source in this setup.Concurrent illumination with the RGB light sources was confirmed through spectrogram measurements conducted using a spectrometer (model HR4000CG-UV-NIR, Ocean Optics, Orlando, FL), as presented in Figure 2.  As illustrated in Figure 1b, the camera utilized in this study was a color polarization camera model DYK 33GX250 from The Imaging Source Asia Co., Ltd., (Shenzhen, China), featuring a Sony IMX250MYR-C CMOS (color) image sensor.The camera boasted a resolution of 2448 × 2048 pixels (5 MP) with a pixel size of 3.45 µm × 3.45 µm and operated at a frame rate of 75 fps.The camera employed a telecentric lens (model FT03X110, Hikvision Digital Technology Co., Ltd., Hangzhou, China), providing a field of view for the object of 36.67 mm and a working distance of 110 mm.The camera's imaging system comprised a polarized superimage element formed by four wire-grid polarizer pixels oriented horizontally (0 • ), vertically (90 • ), and diagonally (45 • and 135 • ).Each polarized superimage element was equipped with a mosaic of Bayer filters in a sequential arrangement, as depicted in the inset of Figure 1a.The RGB-LEDs utilized in the setup were sourced from KOMA Vision Technology Co., Ltd., (Dongguan, China), illuminating a surface area of 100 × 100 mm.The LEDs emit natural light.A circular polarizer was fabricated by laminating a linear film polarizer onto a high-contrast polymer achromatic quarter-wave retarder sheet, covering the wavelength range of 400-700 nm, sourced from Edmund Optics (Shenzhen, China) Co., Ltd.The ellipticity uniformity of the circularly polarized light passing through the polarizer exceeds 99%.To shield the light source between the LEDs and the sample, a fused silica window measuring 100 × 100 mm with a thickness of 18 mm was installed.This window protects the circular polarizer from scratches and contamination during sample measurement while providing a smooth, light-transmitting support platform for the sample.The spectral response characteristics of the color camera were taken into account when configuring the RGB three-color light source in this setup.Concurrent illumination with the RGB light sources was confirmed through spectrogram measurements conducted using a spectrometer (model HR4000CG-UV-NIR, Ocean Optics, Orlando, FL, USA), as presented in Figure 2. As illustrated in Figure 1b, the camera utilized in this study was a color polarization camera model DYK 33GX250 from The Imaging Source Asia Co., Ltd.(Shenzhen), featuring a Sony IMX250MYR-C CMOS (color) image sensor.The camera boasted a resolution of 2448 × 2048 pixels (5 MP) with a pixel size of 3.45 µm × 3.45 µm and operated at a frame rate of 75 fps.The camera employed a telecentric lens (model FT03X110, Hikvision Digital Technology Co., Ltd.Hangzhou), providing a field of view for the object of 36.67 mm and a working distance of 110 mm.The camera's imaging system comprised a polarized superimage element formed by four wire-grid polarizer pixels oriented horizontally (0°), vertically (90°), and diagonally (45° and 135°).Each polarized superimage element was equipped with a mosaic of Bayer filters in a sequential arrangement, as depicted in the inset of Figure 1a.The RGB-LEDs utilized in the setup were sourced from KOMA Vision Technology Co., Ltd.(Dongguan), illuminating a surface area of 100 × 100 mm.The LEDs emit natural light.A circular polarizer was fabricated by laminating a linear film polarizer onto a high-contrast polymer achromatic quarter-wave retarder sheet, covering the wavelength range of 400-700 nm, sourced from Edmund Optics (Shenzhen) Co., Ltd.The ellipticity uniformity of the circularly polarized light passing through the polarizer exceeds 99%.To shield the light source between the LEDs and the sample, a fused silica window measuring 100 × 100 mm with a thickness of 18 mm was installed.This window protects the circular polarizer from scratches and contamination during sample measurement while providing a smooth, light-transmitting support platform for the sample.The spectral response characteristics of the color camera were taken into account when configuring the RGB three-color light source in this setup.Concurrent illumination with the RGB light sources was confirmed through spectrogram measurements conducted using a spectrometer (model HR4000CG-UV-NIR, Ocean Optics, Orlando, FL), as presented in Figure 2.  The RGB-LED employed in this study emitted light at three distinct wavelengths: 463 nm, 533 nm, and 629 nm.Following the passage through the circular polarizer, the light emitted by the RGB-LED was transformed into a circularly polarized illumination source.Due to the ellipticity uniformity exceeding 99%, the difference between the intensity of the circularly polarized light component and the total light intensity is less than 1%.This transformation can be characterized using the Stokes vector as where I in(λ) represents the intensity of circularly polarized light at the wavelength λ.The polarization transmission characteristics of the sample can be described using the Mueller matrix [24] as follows: In the aforementioned equation, ρ denotes the fast-axis azimuthal angle of birefringence, and X represents the phase retardance corresponding to the birefringence within the sample.This phase retardance can be expressed for various detection light wavelength, as X (λ) = 2πR/λ = 2π∆n (λ) d/λ.R represents the retardance.The dispersion in the birefringence of the sample for different detection light wavelengths is represented by ∆n (λ) , and d indicates the thickness of the sample.Upon transmission through the sample, the circularly polarized light is transformed into elliptically polarized light.This elliptically polarized light is described by the Stokes vector as S out(λ)(i) = M S S in(λ) .The circularly polarized light S in(λ) after passing through the sample is defined in Equation ( 1).The resulting elliptically polarized light transmitted through the sample is then solved as follows: Utilizing a polarization camera for single imaging detection, each polarized superimage element captured the light intensity parameters in four polarization directions: horizontal (0 • ), vertical (90 • ), and two diagonal (45 • and 135 • ), denoted as I 0 , I 45 , I 90 , and I 135 , respectively.The first three components of the Stokes vector, representing orthogonal radiation intensity, the difference between horizontal and vertical intensities, and the disparity between the two diagonal intensities, were concurrently measured.Therefore, the polarization camera detected the elliptically polarized light transmitted through the sample, expressed as S out(λ) = M[I 0 I 45 I 90 I 135 ] T , where M ij = [0.5, 0.5, 0.5, 0.5; 1, 0, −1, 0; 0, 1, 0, −1] is the polarization response matrix.The light intensity acquired from each polarized superimage element satisfies the following relation: From Equation (3), two relative light intensity values were further obtained as in Equation (4).
The numerical simulation depicting the variation in normal intensity in relation to the retardance value R = λX (λ) /2π across the three RGB wavelength channels, 463 nm, 533 nm, and 629 nm, was shown in Figure 3.  vs. retardance R.
In Figure 3, the relationship between normal intensity and retardance R was periodic.The three RGB wavelengths displayed distinct patterns of transmitted light intensity at identical values of retardance, while varying levels of retardance exhibited graded characteristics.By employing Equation ( 4), the phase retardance X and the fastaxis azimuthal angle  of the sample were concurrently determined.
The phase retardation could be expressed further through the normal intensity for the first-and second-order light intensities, as depicted in Figure 3.

I( )
I( ) For the majority of optical materials, birefringent dispersion is typically minimal [25], and the variance in birefringence retardance across shorter spectral ranges is sufficiently minute to serve as a constraint for precise retardance determination at each wavelength.Moreover, the range of the fast-axis azimuthal angle can be ascertained concurrently with retardance determination.When (0, / 2) R   , according to light intensity described in Equation (3), whether the light intensity is positive or negative, the range of fast-axis azimuthal angle can be determined as follows: In Figure 3, the relationship between normal intensity γ I(λ) and retardance R was periodic.The three RGB wavelengths displayed distinct patterns of transmitted light intensity at identical values of retardance, while varying levels of retardance exhibited graded characteristics.By employing Equation ( 4), the phase retardance X and the fast-axis azimuthal angle ρ of the sample were concurrently determined.
The phase retardation could be expressed further through the normal intensity γ I(λ) for the first-and second-order light intensities, as depicted in Figure 3.
For the majority of optical materials, birefringent dispersion is typically minimal [25], and the variance in birefringence retardance across shorter spectral ranges is sufficiently minute to serve as a constraint for precise retardance determination at each wavelength.Moreover, the range of the fast-axis azimuthal angle can be ascertained concurrently with retardance determination.When R ∈ (0, λ/2), according to light intensity described in Equation (3), whether the light intensity is positive or negative, the range of fast-axis azimuthal angle can be determined as follows: When R ∈ (λ/2, λ), according to light intensity described in Equation (3), the range of fast-axis azimuthal angle can be determined as follows: Thus, a single imaging measurement enables the simultaneous determination of birefringence retardance and fast-axis azimuthal angle at multi-wavelengths.

Camera Polarization Response Calibration
The precise calibration of the camera's polarization response is crucial due to several factors: the non-uniform photoelectric response exhibited by the CMOS image sensor pixels and the variations in extinction ratio and transmittance among the wire-grid polarizer pixels.Initially, the circular polarizer was removed from the detection light source.Subsequently, without introducing any sample, a linear polarizer was positioned between the RGB-LED and the polarization camera.The calibration procedure followed the methodology delineated in reference [26].The angle of the transmittance axis (θ) was systematically adjusted from 0 • to 180 • in 2 • increments.Calibration was carried out to derive the response matrices for each polarized superimage element, as illustrated in Figure 4a-c.Based on the findings presented in Figure 4, the response matrices for each polarized superimage element were calibrated for the RGB color polarization camera.Each color polarization yielded 512 × 612 polarization superimages, with each superimage corresponding to four measurement matrix subunits containing 3 × 4 polarization directions.The calibration process effectively mitigated measurement errors in light intensity.To illustrate, the polarized light detection results before and after calibration for the three color channels are documented in Figure 5, with the pixel (256, 306) superimage element serving as an exemplar.Based on the findings presented in Figure 4, the response matrices for each polarized superimage element were calibrated for the RGB color polarization camera.Each color polarization yielded 512 × 612 polarization superimages, with each superimage corresponding to four measurement matrix subunits containing 3 × 4 polarization directions.The calibration process effectively mitigated measurement errors in light intensity.To illustrate, the polarized light detection results before and after calibration for the three color channels are documented in Figure 5, with the pixel (256, 306) superimage element serving as an exemplar.The results presented in Figure 5(a1-c1) indicate that the light intensity measurements were notably more accurate when employing the calibrated response matrices.For instance, considering the incident light at the 0° direction of the polarizer, the degree of linear polarization (DOLP) for the three color measurements-(a2,a3) for B-, (b2,b3) for G-, and (c2,c3) for R-color channel-before and after the calibration of the response matrices were examined.Specifically, the B-color DOLP improved from 0.972 to 1.005, the G-color DOLP improved from 0.992 to 1.003, and the R-color DOLP improved from 0.957 to 1.004.Moreover, the response fluctuation of the entire superimage arrays enhanced from 0.001 to 0.0004, thus ensuring the high accuracy of birefringence measurement applications.

Birefringence Measurements Conducted on a Quarter-Wave Plate and a PMMA Arch-Shaped Workpiece
The accuracy and repeatability of the system were evaluated by employing a polymer true zero-order quarter-wave plate as a standard sample.The wave plate was a QWP25-633A-M quarter-wave plate manufactured by Shenzhen LUBON Technology Co., Ltd., and designed for operation at a working wavelength of 633 nm, featuring an optical aperture diameter of 21.5 mm.According to the manufacturer's measurement and calibration The results presented in Figure 5(a1-c1) indicate that the light intensity measurements were notably more accurate when employing the calibrated response matrices.For instance, considering the incident light at the 0 • direction of the polarizer, the degree of linear polarization (DOLP) for the three color measurements-(a2,a3) for B-, (b2,b3) for G-, and (c2,c3) for R-color channel-before and after the calibration of the response matrices were examined.Specifically, the B-color DOLP improved from 0.972 to 1.005, the G-color DOLP improved from 0.992 to 1.003, and the R-color DOLP improved from 0.957 to 1.004.Moreover, the response fluctuation of the entire superimage arrays enhanced from 0.001 to 0.0004, thus ensuring the high accuracy of birefringence measurement applications.

Birefringence Measurements Conducted on a Quarter-Wave Plate and a PMMA Arch-Shaped Workpiece
The accuracy and repeatability of the system were evaluated by employing a polymer true zero-order quarter-wave plate as a standard sample.The wave plate was a QWP25-633A-M quarter-wave plate manufactured by Shenzhen LUBON Technology Co., Ltd., Shenzhen, China and designed for operation at a working wavelength of 633 nm, featuring an optical aperture diameter of 21.5 mm.According to the manufacturer's measurement and calibration data, the retardance of the wave plate is 188.9 nm at 463 nm (B), 170.0 nm at 533 nm (G), and 157.3 nm at 629 nm (R) across the three wavelength channels [27].
The wave plate was placed in our system, its fast-axis orientation was adjusted to approximately 22.5 • , and then the measurement was started.The measurement results for the three RGB wavelength channels obtained using our system are shown in Figure 6.
Sensors 2024, 24, x FOR PEER REVIEW 9 of 14 data, the retardance of the wave plate is 188.9 nm at 463 nm (B), 170.0 nm at 533 nm (G), and 157.3 nm at 629 nm (R) across the three wavelength channels [27].
The wave plate was placed in our system, its fast-axis orientation was adjusted to approximately 22.5°, and then the measurement was started.The measurement results for the three RGB wavelength channels obtained using our system are shown in Figure 6.In Figure 6, the retardance and fast-axis azimuthal angle measurements are uniformly distributed across the entire aperture of the wave plate.A single imaging measurement sufficed to obtain the birefringence retardance and fast-axis azimuthal angles for all three spectral wavelength channels of RGB.The average and standard deviation of the retardance and fast-axis azimuthal angle within the effective aperture area of the wave plate (horizontal: 86 to 426 pixels; vertical: 145 to 485 pixels) are presented in Table 1 for the three spectral wavelength channels.We also included the calibration data provided by the manufacturer in Table 1.Based on the test results of the true zero-order quarter-wave plate, noticeable differences in birefringence were observed across the three wavelength channels, primarily attributable to spectral dispersion birefringence.The measurement value for the R channel was very close to the retardance value at 633 nm wavelength because their wavelengths were similar.The retardance value decreased as the wavelength increased.The maximum relative deviation between the measured values in the three spectral channels and the calibration values provided by the manufacturer is 1.1% in the R channel, calculated as |159.06− 157.3|/157.3= 1.1%.This indicated that the system achieved high measurement accuracy.The standard deviation of birefringence measurements did not exceed 0.67 nm in any of the spectral wavelength channels, indicating a high level of repeatability.Additionally, the fast-axis azimuthal angle remained consistent across multispectral measurements, with a standard deviation of the angle not exceeding 0.08°.In Figure 6, the retardance and fast-axis azimuthal angle measurements are uniformly distributed across the entire aperture of the wave plate.A single imaging measurement sufficed to obtain the birefringence retardance and fast-axis azimuthal angles for all three spectral wavelength channels of RGB.The average and standard deviation of the retardance and fast-axis azimuthal angle within the effective aperture area of the wave plate (horizontal: 86 to 426 pixels; vertical: 145 to 485 pixels) are presented in Table 1 for the three spectral wavelength channels.We also included the calibration data provided by the manufacturer in Table 1.Based on the test results of the true zero-order quarter-wave plate, noticeable differences in birefringence were observed across the three wavelength channels, primarily attributable to spectral dispersion birefringence.The measurement value for the R channel was very close to the retardance value at 633 nm wavelength because their wavelengths were similar.The retardance value decreased as the wavelength increased.The maximum relative deviation between the measured values in the three spectral channels and the calibration values provided by the manufacturer is 1.1% in the R channel, calculated as |159.06− 157.3|/157.3= 1.1%.This indicated that the system achieved high measurement accuracy.The standard deviation of birefringence measurements did not exceed 0.67 nm in any of the spectral wavelength channels, indicating a high level of repeatability.Additionally, the fast-axis azimuthal angle remained consistent across multispectral measurements, with a standard deviation of the angle not exceeding 0.08 • .
To further validate the feasibility of snapshot multispectral birefringence measurements, stress birefringence in an arch-shaped workpiece fabricated from polymethyl methacrylate (PMMA) was assessed both before and after the application of force, as depicted in Figure 7.
Sensors 2024, 24, x FOR PEER REVIEW 10 of 14 To further validate the feasibility of snapshot multispectral birefringence measurements, stress birefringence in an arch-shaped workpiece fabricated from polymethyl methacrylate (PMMA) was assessed both before and after the application of force, as depicted in Figure 7.The PMMA arch-shaped workpiece, as depicted in Figure 7a, had a thickness of 5 mm.Birefringence distribution was measured for both unclamped and clamped states, with clamping applied at the bottom of the arch-shaped workpiece.The distribution of birefringence within the workpiece was non-uniform prior to the application of force.Figures 7(b1-b3) distinctly illustrate significant and highly heterogeneous edge stresses, primarily attributed to residual stresses resulting from the laser ablation processing of the workpiece.The birefringence retardance within the workpiece did not exceed 120 nm. Figure 7(b4-b6) show the birefringence distribution within the workpiece being divided The PMMA arch-shaped workpiece, as depicted in Figure 7a, had a thickness of 5 mm.Birefringence distribution was measured for both unclamped and clamped states, with clamping applied at the bottom of the arch-shaped workpiece.The distribution of birefringence within the workpiece was non-uniform prior to the application of force.Figure 7(b1-b3) distinctly illustrate significant and highly heterogeneous edge stresses, primarily attributed to residual stresses resulting from the laser ablation processing of the workpiece.The birefringence retardance within the workpiece did not exceed 120 nm.
Figure 7(b4-b6) show the birefringence distribution within the workpiece being divided into two sections from the midpoint of the arch.At the midpoint, the fast-axis azimuthal angle of the birefringence was approximately 90 • .Subsequently, the left half of the workpiece exhibited a gradual transition from 90 • to 0 in accordance with the shape of the arch, while the right half transitioned from 90 • to 180 • following the arch's curvature.
Figure 7(c1-c3) clearly illustrate an increase in birefringence following the application of force.However, there was a smaller region of birefringence distribution along the midpoint of the arch.Consequently, the arch exhibited a layered sub-regional arrangement along its inner and outer edges.Notably, the fast-axis azimuthal angle of birefringence at the midpoint of the outer gable arch approached 0, and the intersection point at 90 • shifted to where the arch intersects the vertical straight line.Meanwhile, the inner 90 • remained centrally positioned.We recorded the change in birefringence value at the outer eave vertex of the arched workpiece, approximately 1 mm from the region of interest designated as point A. The region of interest measuring about 3 × 3 pixels was recorded for further analysis.Birefringence measurements are detailed in Table 2. From the statistical data recorded in Table 2, it was evident that the birefringence value at the region of interest, point A, on the arch-shaped workpiece increased following the application of external force.Additionally, there was a discernible dispersion in the distribution of birefringence retardance across the three spectral wavelength channels.Notably, the direction of the fast-axis distribution shifted from its original vertical orientation by approximately 90 • to a horizontal orientation of about 0 • (180 • ).

Discussion
Through the precise calibration of the response matrix for each polarized superimage element, the polarization directions (0 • , 45 • , 90 • , and 135 • ) of light intensity were accurately detected across the three RGB wavelength channels (463 nm, 533 nm, and 629 nm) in a single imaging session.This calibration facilitated highly accurate and sensitive birefringent snapshot imaging measurements.The measurement process, exemplified by the true zero-order polymer quarter-wave calibration, resulted in standard deviation values for birefringence retardance and fast-axis azimuthal angle of better than 0.67 nm and 0.08 • , respectively, indicating high repeatability and sensitivity.The birefringence distribution of materials such as wave-plates and PMMA could be rapidly assessed using the snapshot imaging techniques.The RGB three-channel birefringence retardance measurements exhibited spectral dispersion characteristics.However, the directional distribution of birefringence fast-axis azimuthal angle across the three-channel measurements remained consistent and did not follow spectral dispersion patterns.Given the camera's frame rate of 75 fps, the measurement time could be as short as 13 ms if adequate light intensity for illumination is provided.In addition, the lens utilized in this investigation had a field of view of 36.67 mm.Enhancing the system's lens with a broader field of view would enable an extension of the aperture for birefringence imaging detection.The test results of this study are compared with those in the existing literature and are listed in Table 3.This approach enabled high-precision and high-repeatability birefringence imaging measurements across three wavelength channels in a single imaging session.In the singlewavelength channel measurement method, the birefringence retardance range is limited to 0-λ/2, whereas multi-wavelength measurements can extend the retardance measurement range to 0-λ.This multi-wavelength measurement technique is applicable to the study of material birefringence dispersion.The snapshot single-imaging birefringence measurement achieved in this work makes the method particularly suitable for measuring high-speed dynamic changes in birefringence distribution.

Conclusions
In summary, a study focusing on birefringence imaging measurements in RGB multiwavelength channels was conducted.The experimental setup included an RGB-LED array incorporating a circular polarizer and a fused silica window, which generated a uniform circularly polarized surface-array detection light source for illuminating the sample.The outgoing light was then captured by a color polarized camera.Imaging measurements of the outgoing light (S 0 , S 1 , and S 2 ) were performed through a single measurement of light intensity across the polarization directions of the three RGB spectral channels (0 • , 45 • , 90 • , and 135 • ).The birefringence retardance and fast-axis azimuthal angle were determined through a combination of measurements from the three wavelength channels simultaneously.We conducted a thorough analysis of the measurement principle and constructed the experimental system devices.In the experiments, we initially completed the precise calibration of the measurement response matrices of the polarization camera across the RGB channels.Subsequently, the scheme's measurement accuracy and repeatability were demonstrated using the polymer true zero-order quarter-wave plate.We further investigated the stress birefringence distribution before and after applying force to the PMMA workpiece.The experimental results indicated standard deviations of birefringence retardance and fastaxis azimuthal angle better than 0.67 nm and 0.08 • , respectively, affirming the high level of repeatability and sensitivity.The method showcased that a single imaging measurement could effectively capture the birefringence distribution across multiple-wavelength channels of RGB.The achieved multispectral, high-speed, high-precision, and high-repeatability birefringence imaging measurements offer an effective method for conducting birefringence assessments of optical materials and components.Additionally, they facilitate analysis and evaluation of dispersion.
The numerical simulation depicting the variation in normal intensity in relation to the retardance value () /2 RX   = across the three RGB wavelength channels, 463 nm, 533 nm, and 629 nm, was shown in Figure3.

Figure 6 .
Figure 6.Results of quarter-wave plate: (a-c) retardance at B-, G-, and R-wavelength channels and (d-f) fast-axis azimuthal angle at B-, G-, and R-wavelength channels.

Figure 6 .
Figure 6.Results of quarter-wave plate: (a-c) retardance at B-, G-, and R-wavelength channels and (d-f) fast-axis azimuthal angle at B-, G-, and R-wavelength channels.

Table 3 .
Comparison with reported birefringence measurement methods.