根据四个偏振角度的偏振图像计算偏振斯托克斯矢量-编程阁

根据四个偏振角度的偏振图像计算偏振斯托克斯矢量（Stokes Vector），遵循图像预处理→强度提取→斯托克斯分量计算的流程。

斯托克斯矢量是描述偏振光状态的核心工具，通过四个实数参数（S0,S1,S2,S3）完整表征光的总强度、线偏振/圆偏振成分及方向。

一、斯托克斯矢量的基本概念

斯托克斯矢量的四个分量定义如下（主流约定）：

S0S_0S0：总光强，即四个偏振角度图像的平均强度（或总和，取决于定义）；
S1S_1S1：水平与垂直线偏振强度差（I0−I90I_0−I_{90}I0−I90），反映线偏振在水平/垂直方向的偏好；
S2S_2S2：±45°线偏振强度差（I45−I135I_{45}−I_{135}I45−I135），反映线偏振在±45°方向的偏好；
S3S_3S3：右旋与左旋圆偏振强度差（IR−ILI_R−I_LIR−IL），反映圆偏振的手性（本题中未涉及圆偏振片，故S3=0S_3=0S3=0）。

其中，I0,I45,I90,I135I_0,I_{45},I_{90},I_{135}I0,I45,I90,I135分别为0°（水平）、45°、90°（垂直）、135°偏振角度下的图像强度。

二、实现步骤与MATLAB代码

1. 图像预处理

读取图像：使用imread读取四个偏振角度的偏振图像（假设为JPG/PNG格式）；
灰度转换：将彩色图像转换为灰度图像（偏振信息主要包含在强度中）；
归一化：将灰度值归一化到[0,1]区间（避免亮度差异影响计算）。

2. 强度提取

直接从预处理后的灰度图像中提取每个像素的强度值（I0,I45,I90,I135I_0,I_{45},I_{90},I_{135}I0,I45,I90,I135）。

3. 斯托克斯分量计算

根据斯托克斯矢量的定义，计算每个像素的四个分量：

S0=I0+I902,S1=I0−I90,S2=I45−I135,S3=0S_0=\frac{I_0+I_{90}}{2},S_1=I_0−I_{90},S_2=I_{45}−I_{135},S_3=0S0=2I0+I90,S1=I0−I90,S2=I45−I135,S3=0

三、MATLAB代码实现

MATLAB代码，包含图像读取、预处理、斯托克斯分量计算及结果可视化：

%% 1. 参数设置偏振角度列表（0°、45°、90°、135°） angles=[0,45,90,135];图像文件名（需根据实际路径修改） image_files={'polar_0deg.jpg','polar_45deg.jpg','polar_90deg.jpg','polar_135deg.jpg'};%% 2. 图像预处理与强度提取初始化存储强度的矩阵 I=zeros([size(imread(image_files{1}),1),size(imread(image_files{1}),2),4]);fori=1:4% 读取图像img=imread(image_files{i});% 转换为灰度图像（如果是彩色图像）ifsize(img,3)==3img_gray=rgb2gray(img);elseimg_gray=img;end% 归一化到[0,1]区间I(:,:,i)=double(img_gray)/255;end% 提取各角度强度I0=I(:,:,1);% 0°偏振强度I45=I(:,:,2);% 45°偏振强度I90=I(:,:,3);% 90°偏振强度I135=I(:,:,4);% 135°偏振强度%% 3. 计算斯托克斯分量S0=(I0+I90)/2;% 总光强S1=I0-I90;% 水平-垂直线偏振差S2=I45-I135;% ±45°线偏振差S3=zeros(size(S0));% 圆偏振差（无圆偏振片时设为0）%% 4. 结果可视化figure('Position',[100,100,1200,800]);% 显示原始偏振图像subplot(2,3,1);imshow(I0);title('0°偏振图像');subplot(2,3,2);imshow(I45);title('45°偏振图像');subplot(2,3,3);imshow(I90);title('90°偏振图像');subplot(2,3,4);imshow(I135);title('135°偏振图像');% 显示斯托克斯分量subplot(2,3,5);imshow(S0,[]);title('斯托克斯分量S0（总光强）');subplot(2,3,6);imshow(S1,[]);title('斯托克斯分量S1（水平-垂直差）');figure('Position',[100,100,800,400]);subplot(1,2,1);imshow(S2,[]);title('斯托克斯分量S2（±45°差）');subplot(1,2,2);imshow(S3,[]);title('斯托克斯分量S3（圆偏振差）');%% 5. 保存结果（可选）% 将斯托克斯分量保存为MAT文件save('stokes_vector.mat','S0','S1','S2','S3');% 将斯托克斯分量保存为图像文件imwrite(S0,'S0.png');imwrite(S1,'S1.png');imwrite(S2,'S2.png');imwrite(S3,'S3.png');