结构光正弦条纹 格雷码图案matlab生成程序

clc;close all;clear all; % *************************************载入原始图像*************************************% row=1140;%行 col=912;%列 n=2; f=n/228;%频率，条纹数为n T=1/f/2; for i=1:row for j=1:col % *************************************格雷码1*************************************% if(j<=T*8) Img1(i,j)=0; else Img1(i,j)=1; end % *************************************格雷码2*************************************% if(j<=T*4) Img2(i,j)=0; else if(j<=T*8) Img2(i,j)=1; else if(j<=T*12) Img2(i,j)=1; else Img2(i,j)=0; end end end % *************************************格雷码3*************************************% if(j<=T*2) Img3(i,j)=0; else if(j<=T*4) Img3(i,j)=1; else if(j<=T*6) Img3(i,j)=1; else if(j<=T*8) Img3(i,j)=0; else if(j<=T*10) Img3(i,j)=0; else if(j<=T*12) Img3(i,j)=1; else if(j<=T*14) Img3(i,j)=1; else Img3(i,j)=0; end end end end end end end % *************************************格雷码4*************************************% if(j<=T) Img4(i,j)=0; else if(j<T*2) Img4(i,j)=1; else if(j<=T*3) Img4(i,j)=1; else if(j<=T*4) Img4(i,j)=0; else if(j<=T*5) Img4(i,j)=0; else if(j<=T*6) Img4(i,j)=1; else if(j<=T*7) Img4(i,j)=1; else if(j<=T*8) Img4(i,j)=0; else if(j<=T*9) Img4(i,j)=0; else if(j<=T*10) Img4(i,j)=1; else if(j<=T*11) Img4(i,j)=1; else if(j<=T*12) Img4(i,j)=0; else if(j<=T*13) Img4(i,j)=0; else if(j<=T*14) Img4(i,j)=1; else if(j<=T*15) Img4(i,j)=1; else Img4(i,j)=0; end end end end end end end end end end end end end end end end end for i=1:row for j=1:col % *************************************反码1*************************************% if (Img1(i,j)==0) Img5(i,j) =1; else Img5(i,j) =0; end % *************************************反码2*************************************% if (Img2(i,j)==0) Img6(i,j) =1; else Img6(i,j) =0; end % *************************************反码3*************************************% if (Img3(i,j)==0) Img7(i,j) =1; else Img7(i,j) =0; end % *************************************反码4*************************************% if (Img4(i,j)==0) Img8(i,j) =1; else Img8(i,j) =0; end end end for i=1:3 MLM1(:,:,i)=Img1; MLM2(:,:,i)=Img2; MLM3(:,:,i)=Img3; MLM4(:,:,i)=Img4; MLM5(:,:,i)=Img5; MLM6(:,:,i)=Img6; MLM7(:,:,i)=Img7; MLM8(:,:,i)=Img8; end imwrite( MLM1,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\GLM1.bmp'); imwrite( MLM2,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\GLM2.bmp'); imwrite( MLM3,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\GLM3.bmp'); imwrite( MLM4,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\GLM4.bmp'); imwrite( MLM5,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\FM1.bmp'); imwrite( MLM6,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\FM2.bmp'); imwrite( MLM7,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\FM3.bmp'); imwrite( MLM8,'C:\\Users\\白景湘\\Desktop\\工作总结材料\\仿真\\FM4.bmp');