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

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

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

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

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

$S_0$ ：总光强，即四个偏振角度图像的平均强度（或总和，取决于定义）；
$S_1$ ：水平与垂直线偏振强度差（ $I_0−I_{90}$ ），反映线偏振在水平/垂直方向的偏好；
$S_2$ ：±45°线偏振强度差（ $I_{45}−I_{135}$ ），反映线偏振在±45°方向的偏好；
$S_3$ ：右旋与左旋圆偏振强度差（ $I_R−I_L$ ），反映圆偏振的手性（本题中未涉及圆偏振片，故 $S_3=0$ ）。

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

二、实现步骤与MATLAB代码

1. 图像预处理

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

2. 强度提取

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

3. 斯托克斯分量计算

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

$S_0=\frac{I_0+I_{90}}{2},S_1=I_0−I_{90},S_2=I_{45}−I_{135},S_3=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]);  

for i = 1:4
    % 读取图像
    img = imread(image_files{i});  
    % 转换为灰度图像（如果是彩色图像）
    if size(img, 3) == 3  
        img_gray = rgb2gray(img);  
    else
        img_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');

四、结果分析

1. 斯托克斯分量的物理意义

S0：总光强图像，反映场景的整体亮度分布（与原始偏振图像的平均强度一致）；
S1：水平与垂直线偏振差，正值表示水平偏振占优，负值表示垂直偏振占优；
S2：±45°线偏振差，正值表示45°偏振占优，负值表示135°偏振占优；
S3：圆偏振差，本题中因未使用圆偏振片，故所有像素值为0（若使用圆偏振片，可提取右旋/左旋圆偏振强度差）。

2. 示例结果

假设输入的偏振图像为室内场景（包含水平桌面、垂直墙面），则：

$S_0$ 图像显示桌面的高亮度（因桌面反射的水平偏振光强）；
$S_1$ 图像中，桌面区域为正值（水平偏振占优），墙面区域为负值（垂直偏振占优）；
$S_2$ 图像中，45°倾斜的物体（如书本）会有明显的正负值分布（取决于其偏振方向）。

五、注意事项

图像配准：四个偏振图像需严格配准（即同一场景、同一视角），否则会导致斯托克斯分量计算误差；
噪声处理：若图像存在噪声，可使用高斯滤波（imgaussfilt）或中值滤波（medfilt2）预处理；
归一化：必须归一化到[0,1]区间，否则不同图像的亮度差异会影响斯托克斯分量的准确性；
圆偏振片：若需计算S3，需添加圆偏振片（如1/4波片+线性偏振片），并采集右旋/左旋圆偏振图像。

参考代码根据四个偏振角度的偏振图像计算偏振斯托斯矢量 www.youwenfan.com/contentcnq/60050.html

六、拓展应用

斯托克斯矢量是偏振成像的核心工具，可用于：

材质识别：不同材质的偏振特性（S1,S2）差异显著（如金属 vs 塑料）；
表面检测：通过S1,S2的分布识别表面划痕、凹陷等缺陷；
目标跟踪：利用偏振信息的稳定性（不受光照强度变化影响），提升目标跟踪的鲁棒性。

总结

通过四个偏振角度的图像计算斯托克斯矢量，需遵循预处理→强度提取→分量计算的流程。MATLAB代码实现了这一流程，并可可视化结果。斯托克斯矢量为偏振成像提供了定量的偏振状态描述，是材质识别、表面检测等领域的关键技术。