%?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