RGB=imread('C:\Users\David\Desktop\Shmily\宋家成\左90右0.bmp');
gray=rgb2gray(RGB);
i90=imcrop(gray,[4,210,493,685]);
figure
imshow(i90);
i0_1=imcrop(gray,[756,220,493,685]);
i0 = i0_1(:, end:-1:1);
figure
imshow(i0);
RGB0=imread('C:\Users\David\Desktop\Shmily\宋家成\左135右45.bmp');
gray0=rgb2gray(RGB0);
i135=imcrop(gray0,[4,210,493,685]);
figure
imshow(i135);
i45_1=imcrop(gray0,[756,220,493,685]);
i45 = i45_1(:, end:-1:1);
figure
imshow(i45);
I=i0+i90;
Q=i0-i90;
U=i45-i135;
V=i0+i90-2*i45;
P1=Q.^2+U.^2+V.^2;
P2=double(P1).^0.5;
P=P2./double(I);
A2=U./Q;
A1=atan2(double(U),double(Q));
A=A1./2
E0=Q.^2+U.^2;
E1=double(E0).^0.5;
E2=double(V)./E1;
E3=tan(E2);
E=1./(2.*E3);
figure
subplot(2,2,1);imshow(I);title('强度图像');
subplot(2,2,2);imshow(3*P);title('偏振度图像');
subplot(2,2,3);imshow(A);title('偏振角图像');
subplot(2,2,4);imshow(E);title('椭圆偏振率图像');
imwrite(I, './强度图像.bmp');
imwrite(3*P, './偏振度图像.bmp');
imwrite(A, './偏振角图像.bmp');
imwrite(E, './椭圆偏振率图像.bmp');
imwrite(i0, './i0.bmp');
imwrite(i45, './i45.bmp');
imwrite(i90, './i90.bmp');
imwrite(i135, './i135.bmp');
评论12