偏振图像处理_斯托克斯

可见光下拍摄偏振角度为0度 60度 120度图像，经斯托克斯处理

% clc;
close all;
clear all;
%%%%%%%%%%% 获取图像 %%%%%%%%%
% M11 = imread（‘ayc.bmp‘）;
% M2 = imread（‘a0c.bmp‘）;
% M3 = imread（‘a60c.bmp‘）;
% M4 = imread（‘a120c.bmp‘）;
M11 = imread（‘a11.bmp‘）;
M2 = imread（‘a110.bmp‘）;
M3 = imread（‘a1260.bmp‘）;
M4 = imread（‘a11120.bmp‘）;

M1 = rgb2gray（M11）;
M2 = rgb2gray（M2）;
M3 = rgb2gray（M3）;
M4 = rgb2gray（M4）;

figure;
subplot（221）;imshow（M1）;title（‘原图灰度图像‘）;
subplot（222）;imshow（M2）;title（‘0度偏振灰度图像‘）;
subplot（223）;imshow（M3）;title（‘60度偏振灰度图像‘）;
subplot（224）;imshow（M4）;title（‘120度偏振灰度图像‘）;
%%%%%%%%%%% 处理图像 %%%%%%%%%%
% I = 2/3.*（M2+M3+M4）;
% Q = 2/3.*（2.*M2-M3-M4）;
% U = 2/（sqrt（3））.*（M3-M4）;
% DOLP = （sqrt（double（Q.^2+U.^2）））./（double（I））;
% q = 1/2*atan（double（U./Q））;

i =2/3*（im2double（M2）+im2double（M3）+im2double（M4））;
I = mat2gray（i）;
q = 2/3*（2*im2double（M2）-im2double（M3）-im2double（M4））;
Q = mat2gray（q）;
u = 2/sqrt（3）*（im2double（M3）-im2double（M4））;
U = mat2gray（u）;
dolp = sqrt（q.*q+u.*u）;
Dolp = dolp./i;
DOLP = mat2gray（Dolp）;
aop = （1/2）*atan（u./q）;
AOP = mat2gray（aop）;

figure;imshow（I）;title（‘I‘）;
figure;imshow（Q）;title（‘Q‘）;
figure;imshow（U）;title（‘U‘）;
figure;imshow（DOLP）;title（‘DOLP‘）;
figure;imshow（AOP）;title（‘AOP‘）;
figure;imshow（M11）;title（‘yuan‘）;

% imgMean（AOP）;
% imgMean（DOLP）;
% 
% X1 = imread（‘aaa.bmp‘）;
% X2 = imread（‘bbb.bmp‘）;
% 
% SimilarityFun（DOLPX1）;
% SimilarityFun（AOPX2）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件    15117622  2019-07-20 10:55  a110.bmp
     文件    15117622  2019-07-20 10:55  a11120.bmp
     文件    15117622  2019-07-20 10:55  a1260.bmp
     文件        1456  2019-08-28 20:45  polarization.m
     文件         104  2019-08-28 20:49  说明.txt
     文件    15117622  2019-07-20 10:55  a11.bmp