pianzhen.zip_偏振椭圆_偏振角_偏振角图像_四 偏振图像_椭圆偏振

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');