aplot_CJC

load('AllSkyMap_25-Aug-2013:21:05:47_26-Aug-2013:02:05:00_data.mat')

0001 function aplot(filename,cLim,cLimP)
0002 
0003 %load('AllSkyMap_25-Aug-2013:21:05:47_26-Aug-2013:02:05:00_data.mat')
0004 load(filename)
0005 
0006 
0007 
0008 fprefix = filename;
0009 %fprefix = sprintf('AllSkyMap_25-Aug-2013:21:30:00_26-Aug-2013:17:00:00_2_');
0010 
0011 
0012 
0013 %
0014 % figure('Position',[1 1000 360*3 90*3])
0015 %      imagesc(azC(:),elC(:),t2,[4.5 4.8])
0016 %      %imagesc(azC(:),elC(:),(CI-repmat(nanmedian(CI,2),1,size(CI,2))))
0017 %      set(gca,'YDir','normal')
0018 %      colorbar
0019 %      xlabel('Azimuth')
0020 %      ylabel('Elevation')
0021 %      title('I=(I1+I2)/2')
0022 %      set(gcf,'PaperPositionMode','auto')
0023 %      eval(sprintf('print -depsc2 -r200 %s',[filename 'I.eps']))
0024 %
0025 % else
0026     figure('Position',[1 1000 360*3 90*3])
0027      %imagesc(azC(:),elC(:),log10(t1),[4.5 4.8])
0028      imagesc(azC(:),elC(:),(CI1-repmat(nanmedian(CI1,2),1,size(CI1,2))),cLim)
0029      set(gca,'YDir','normal')
0030      colorbar
0031      xlabel('Azimuth')
0032      ylabel('Elevation')
0033      title('I=(I1+I2)/2')
0034      set(gcf,'PaperPositionMode','auto')
0035      eval(sprintf('print -depsc2 -r200 %s',[filename 'I1.eps']))
0036 
0037       figure('Position',[1 1000 360*3 90*3])
0038      %imagesc(azC(:),elC(:),log10(t1),[4.5 4.8])
0039      imagesc(azC(:),elC(:),(CI2-repmat(nanmedian(CI2,2),1,size(CI2,2))),cLim)
0040      set(gca,'YDir','normal')
0041      colorbar
0042      xlabel('Azimuth')
0043      ylabel('Elevation')
0044      title('I=(I1+I2)/2')
0045      set(gcf,'PaperPositionMode','auto')
0046      eval(sprintf('print -depsc2 -r200 %s',[filename 'I2.eps']))
0047      
0048 % end
0049 
0050     figure('Position',[1 1000 360*3 90*3])
0051     imagesc(azC(:),elC(:),(CQ1-repmat(nanmedian(CQ1,2),1,size(CQ1,2))),cLimP)
0052      set(gca,'YDir','normal')
0053      colorbar
0054      xlabel('Azimuth')
0055      ylabel('Elevation')
0056      title('Q')
0057      set(gcf,'PaperPositionMode','auto')
0058      eval(sprintf('print -depsc2 -r200 %s',[filename 'Q1.eps']))
0059  
0060      figure('Position',[1 1000 360*3 90*3])
0061      imagesc(azC(:),elC(:),(CU1-repmat(nanmedian(CU1,2),1,size(CU1,2))),cLimP)
0062      set(gca,'YDir','normal')
0063      colorbar
0064      xlabel('Azimuth')
0065      ylabel('Elevation')
0066      title('U')
0067      set(gcf,'PaperPositionMode','auto')
0068      eval(sprintf('print -depsc2 -r200 %s',[filename 'U1.eps']))
0069 
0070 figure('Position',[1 1000 360*3 90*3])
0071      imagesc(azC(:),elC(:),(CP1-repmat(nanmedian(CP1,2),1,size(CP1,2))),cLimP)
0072      set(gca,'YDir','normal')
0073      colorbar
0074      xlabel('Azimuth')
0075      ylabel('Elevation')
0076      title('P = sqrt(Q^2+U^2)')
0077      set(gcf,'PaperPositionMode','auto')
0078      eval(sprintf('print -depsc2 -r200 %s',[filename 'P1.eps']))