当前位置：首页 > news >正文

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

news 2026/5/11 21:35:38

根据四个偏振角度的偏振图像计算偏振斯托克斯矢量（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);  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');