差不多符合老板对画图的要求了。
mark一下:
h_x=1:331;
h_x_l=0.065*h_x;
r=165;
cx1=1997;
cy1=1295;
figure;
subplot(2, 1,1);plot(h_x_l, syn_height(cx1,cy1-r:cy1+r),'r','linewidth',3);xlabel(['X-axis Position/','\mu','m'],'Fontsize',20);ylabel(['Height/','\mu','m'],'fontsize',20);
hold on, plot(h_x_l, Oheight(cx1,cy1-r:cy1+r),'b','linewidth',3);
axis([0 20 -1 11]);set(gca,'FontSize',15);set(gca,'XTick',[0:2:20]);set(gca,'YTick',[-2:2:12]);
legend('Aperture Synthesized Profile','Profile Under Normal Illumination','Location','NorthWest' );
subplot(2, 1,2);plot(h_x_l, syn_height(cx1-r:cx1+r,cy1),'r','linewidth',3);xlabel(['Y-axis Position/','\mu','m'],'Fontsize',20);ylabel(['Height/','\mu','m'],'fontsize',20);
hold on, plot(h_x_l, Oheight(cx1-r:cx1+r,cy1),'b','linewidth',3);
axis([0 20 -1 11]);set(gca,'FontSize',15);set(gca,'XTick',[0:2:20]);set(gca,'YTick',[-2:2:12]);
legend('Aperture Synthesized Profile','Profile Under Normal Illumination','Location','NorthWest' );
s1=Oheight(cx1-r:cx1+r,cy1-r:cy1+r);
sc1=syn_height(cx1-r:cx1+r,cy1-r:cy1+r);
A1=std2(s1(1:60,1:60));
As=std2(sc1(1:60,1:60));
B1=std2(s1(end-60:end,end-60:end));
Bs=std2(sc1(end-60:end,end-60:end));
pb=syn_height;
pb(cx1-r:cx1+r,cy1-r:cy1+r)=255;
figure;imagesc(s1);colormap jet,colorbar('fontsize',30);title('Oheight'),axis equal,axis off,caxis([-2,11]);
figure;imagesc(sc1);colormap jet,colorbar('fontsize',30);title('height'),axis equal,axis off,caxis([-2,11]);
figure;imagesc(pb);colormap jet,colorbar('fontsize',30);title('height'),axis equal,axis off,caxis([-2,11]);
另外polystyrene beads的折射率可以在网站查到:refractiveIndex.info