Description of aplotFrequency

0001 function aplotFrequency(filename,cLim,cLimP)
0002 
0003 %function to plot data produced by analyzeFullSkyAzElMapSouthFrequency
0004 % CJC 10 June 2014
0005 
0006 %filename is the data
0007 %cLim is the intensity colour scal elimits e.g. [0 1000]
0008 %cLimP is the polarization colour scale e.g [-10 +10]
0009 load(filename)
0010 
0011 
0012 for map=1:numberMaps
0013 %    a=figure('visible','off')
0014 %    b=figure('visible','off')
0015 %    set(a,'renderer','zbuffer')
0016 %    set(b,'renderer','zbuffer')
0017     Cmap = T{map};
0018     Nmap = Tn{map};
0019     Cmap = Cmap./repmat(Nmap(:,:,1),[1,1,8]);
0020     CILL1 = Cmap(:,:,1); %%LL
0021     CIRR1 = Cmap(:,:,4); %%RR
0022     CV1 = (CILL1-CIRR1)/2; 
0023     CQ1 = Cmap(:,:,2);
0024     CU1 = Cmap(:,:,3);
0025     CP1 = sqrt(CQ1.^2+CU1.^2);
0026     CI1 = (CILL1+CIRR1)/2;
0027 
0028     CILL2 = Cmap(:,:,5);
0029     CIRR2 = Cmap(:,:,8);
0030     CV2 = (CILL2-CIRR2)/2;
0031     CQ2 = Cmap(:,:,6);
0032     CU2 = Cmap(:,:,7);
0033     CP2 = sqrt(CQ2.^2+CU2.^2);
0034     CI2 = (CILL2+CIRR2)/2;
0035 
0036 
0037      fprefix = filename;
0038      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0039      %imagesc(azC(:),elC(:),log10(t1),[4.5 4.8])
0040      imagesc(azC(:),elC(:),(CI1-repmat(nanmean(CI1,2),1,size(CI1,2))),cLim)
0041      set(gca,'YDir','normal')
0042      colorbar
0043      xlabel('Azimuth')
0044      ylabel('Elevation')
0045      title('Roach 1 I=(I1+I2)/2 ')
0046      set(gcf,'PaperPositionMode','auto')
0047 %     eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach1_I_Chan_' num2str(map) '.eps']))
0048      eval(sprintf('print -dpng -r200 %s',[filename 'Roach1_I_Chan_' num2str(map) '.png']))
0049      close(a)
0050      
0051 % end
0052 
0053      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0054     imagesc(azC(:),elC(:),(CQ1-repmat(nanmean(CQ1,2),1,size(CQ1,2))),cLimP)
0055      set(gca,'YDir','normal')
0056      colorbar
0057      xlabel('Azimuth')
0058      ylabel('Elevation')
0059      title('Roach1 Q')
0060      set(gcf,'PaperPositionMode','auto')
0061  %    eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach1_Q1_Chan_' num2str(map) '.eps']))
0062      eval(sprintf('print -dpng -r200 %s',[filename 'Roach1_Q1_Chan_' num2str(map) '.png']))
0063     close(a) 
0064 
0065      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0066      imagesc(azC(:),elC(:),(CU1-repmat(nanmean(CU1,2),1,size(CU1,2))),cLimP)
0067      set(gca,'YDir','normal')
0068      colorbar
0069      xlabel('Azimuth')
0070      ylabel('Elevation')
0071      title('Roach 1 U')
0072      set(gcf,'PaperPositionMode','auto')
0073   %   eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach1_U1_Chan_' num2str(map) '.eps']))
0074      eval(sprintf('print -dpng -r200 %s',[filename 'Roach1_U1_Chan_' num2str(map) '.png']))
0075     close(a) 
0076 
0077      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0078      imagesc(azC(:),elC(:),(CP1-repmat(nanmean(CP1,2),1,size(CP1,2))),cLimP)
0079      set(gca,'YDir','normal')
0080      colorbar
0081      xlabel('Azimuth')
0082      ylabel('Elevation')
0083      title('Roach 1 P = sqrt(Q^2+U^2)')
0084      set(gcf,'PaperPositionMode','auto')
0085    %  eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach1_P1_Chan_' num2str(map) '.eps']))
0086      eval(sprintf('print -dpng -r200 %s',[filename 'Roach1_P1_Chan_' num2str(map) '.png']))
0087     close(a) 
0088     
0089      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0090      %imagesc(azC(:),elC(:),log10(t1),[4.5 4.8])
0091      imagesc(azC(:),elC(:),(CI2-repmat(nanmean(CI2,2),1,size(CI2,2))),cLim)
0092      set(gca,'YDir','normal')
0093      colorbar
0094      xlabel('Azimuth')
0095      ylabel('Elevation')
0096      title('Roach 2 I=(I1+I2)/2')
0097      set(gcf,'PaperPositionMode','auto')
0098     % eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach2_I1_Chan_' num2str(map) '.eps']))
0099      eval(sprintf('print -dpng -r200 %s',[filename 'Roach2_I1_Chan_' num2str(map) '.png']))
0100     close(a) 
0101 
0102      
0103 % end
0104 
0105      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0106     imagesc(azC(:),elC(:),(CQ2-repmat(nanmean(CQ2,2),1,size(CQ2,2))),cLimP)
0107      set(gca,'YDir','normal')
0108      colorbar
0109      xlabel('Azimuth')
0110      ylabel('Elevation')
0111      title('Roach2 Q')
0112      set(gcf,'PaperPositionMode','auto')
0113     % eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach2_Q1_Chan_' num2str(map) '.eps']))
0114      eval(sprintf('print -dpng -r200 %s',[filename 'Roach2_Q1_Chan_' num2str(map) '.png']))
0115     close(a) 
0116  
0117      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0118      imagesc(azC(:),elC(:),(CU2-repmat(nanmean(CU2,2),1,size(CU2,2))),cLimP)
0119      set(gca,'YDir','normal')
0120      colorbar
0121      xlabel('Azimuth')
0122      ylabel('Elevation')
0123      title('Roach2 U')
0124      set(gcf,'PaperPositionMode','auto')
0125     % eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach2_U1_Chan_' num2str(map) '.eps']))
0126      eval(sprintf('print -dpng -r200 %s',[filename 'Roach2_U1_Chan_' num2str(map) '.png']))
0127     close(a) 
0128 
0129      a=figure('Position',[1 1000 360*3 90*3],'visible','off')
0130      imagesc(azC(:),elC(:),(CP2-repmat(nanmean(CP2,2),1,size(CP2,2))),cLimP)
0131      set(gca,'YDir','normal')
0132      colorbar
0133      xlabel('Azimuth')
0134      ylabel('Elevation')
0135      title('Roach2 P = sqrt(Q^2+U^2)')
0136      set(gcf,'PaperPositionMode','auto')
0137     % eval(sprintf('print -depsc2 -r200 %s',[filename 'Roach2_P1_Chan_' num2str(map) '.eps']))
0138      eval(sprintf('print -dpng -r200 %s',[filename 'Roach2_P1_Chan_' num2str(map) '.png']))
0139     close(a) 
0140 
0141 end
aplotFrequency

PURPOSE

SYNOPSIS

DESCRIPTION

CROSS-REFERENCE INFORMATION

SOURCE CODE
