analyzeFullSkyAzElMap_saveFigures

 analyzeFullSkyAzElMap_saveFigures(fname,dirname)

 This function plots the maps saved by analyzeFullSkyAzElMap. Why have a
 separate function? Because some matlab installations (*cough* haggis
 *cough*) save awful looking maps for some reason.

 ogk 18 June 2011

0001 function analyzeFullSkyAzElMap_saveFigures(varargin)
0002 % analyzeFullSkyAzElMap_saveFigures(fname,dirname)
0003 %
0004 % This function plots the maps saved by analyzeFullSkyAzElMap. Why have a
0005 % separate function? Because some matlab installations (*cough* haggis
0006 % *cough*) save awful looking maps for some reason.
0007 %
0008 % ogk 18 June 2011
0009 %
0010 
0011 fname = varargin{1};
0012 dirname = varargin{2};
0013 
0014 eval(sprintf('load %s',fname))
0015 
0016 if nargin == 3
0017     climI = varargin{3};
0018 end
0019 
0020 fprefix(fprefix=='_')='-';
0021 
0022 % If prefix is AllSkyMap-day-month-year-, change it to
0023 % AllSkyMap-year-month-day-
0024 A = strsplit('-',fprefix,'include');
0025 if str2double(A{3}) > 31
0026     yr = A{3};
0027     dy = A{7};
0028     mn = A{5};
0029 else
0030     yr = A{7};
0031     dy = A{3};
0032     mn = A{5};
0033 end
0034 fprefix = ['AllSkyMap-' yr '-' mn '-' dy '-'];
0035 
0036 figure('Position',[1 1000 360*3 90*3])
0037 imagesc(azC(:),elC(:),(CI1-repmat(nanmedian(CI1,2),1,size(CI1,2))),climI)
0038 set(gca,'YDir','normal')
0039 colorbar
0040 xlabel('Azimuth')
0041 ylabel('Elevation')
0042 title([fprefix 'I1'])
0043 set(gcf,'PaperPositionMode','auto')
0044 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'I1.eps']))
0045 
0046 figure('Position',[1 1000 360*3 90*3])
0047 imagesc(azC(:),elC(:),(CI2-repmat(nanmedian(CI2,2),1,size(CI2,2))),climI)
0048 set(gca,'YDir','normal')
0049 colorbar
0050 xlabel('Azimuth')
0051 ylabel('Elevation')
0052 title([fprefix 'I2'])
0053 set(gcf,'PaperPositionMode','auto')
0054 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'I2.eps']))
0055 
0056 figure('Position',[1 1000 360*3 90*3])
0057 imagesc(azC(:),elC(:),(CI-repmat(nanmedian(CI,2),1,size(CI,2))),climI)
0058 set(gca,'YDir','normal')
0059 colorbar
0060 xlabel('Azimuth')
0061 ylabel('Elevation')
0062 title([fprefix 'I=(I1+I2)/2'])
0063 set(gcf,'PaperPositionMode','auto')
0064 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'I.eps']))
0065 
0066 figure('Position',[1 1000 360*3 90*3])
0067 imagesc(azC(:),elC(:),(CV-repmat(nanmedian(CV,2),1,size(CV,2))),climV)
0068 set(gca,'YDir','normal')
0069 colorbar
0070 xlabel('Azimuth')
0071 ylabel('Elevation')
0072 title([fprefix 'V=(I1-I2/2)'])
0073 set(gcf,'PaperPositionMode','auto')
0074 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'V.eps']))
0075 
0076 figure('Position',[1 1000 360*3 90*3])
0077 imagesc(azC(:),elC(:),(CQ-repmat(nanmedian(CQ,2),1,size(CQ,2))),climP)
0078 set(gca,'YDir','normal')
0079 colorbar
0080 xlabel('Azimuth')
0081 ylabel('Elevation')
0082 title([fprefix 'Q'])
0083 set(gcf,'PaperPositionMode','auto')
0084 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'Q.eps']))
0085 
0086 figure('Position',[1 1000 360*3 90*3])
0087 imagesc(azC(:),elC(:),(CU-repmat(nanmedian(CU,2),1,size(CU,2))),climP)
0088 set(gca,'YDir','normal')
0089 colorbar
0090 xlabel('Azimuth')
0091 ylabel('Elevation')
0092 title([fprefix 'U'])
0093 set(gcf,'PaperPositionMode','auto')
0094 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'U.eps']))
0095 
0096 figure('Position',[1 1000 360*3 90*3])
0097 imagesc(azC(:),elC(:),(CP-repmat(nanmedian(CP,2),1,size(CP,2))),[0 max(climP)])
0098 set(gca,'YDir','normal')
0099 colorbar
0100 xlabel('Azimuth')
0101 ylabel('Elevation')
0102 title([fprefix 'P = sqrt(Q^2+U^2)'])
0103 set(gcf,'PaperPositionMode','auto')
0104 eval(sprintf('print -depsc2 -r200 %s',[dirname fprefix 'P.eps']))
0105 
0106 fid = fopen([dirname fprefix 'Plots.tex'],'w');
0107 fprintf(fid,'\\documentclass[10pt]{article}\n');
0108 fprintf(fid,'\\usepackage{graphicx,subfigure}\n');
0109 fprintf(fid,'\\usepackage[letterpaper,inner=1cm,outer=1cm,top=1cm,bottom=1cm,pdftex]{geometry}\n');
0110 fprintf(fid,'\\begin{document}\n');
0111 fprintf(fid,'\\begin{figure}[t!] \\centering\n');
0112 fprintf(fid,'\\subfigure[]{\n');
0113 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'I.pdf']);
0114 fprintf(fid,'\\subfigure[]{\n');
0115 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'Q.pdf']);
0116 fprintf(fid,'\\subfigure[]{\n');
0117 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'U.pdf']);
0118 fprintf(fid,'\\subfigure[]{\n');
0119 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'P.pdf']);
0120 fprintf(fid,'\\end{figure}\n');
0121 fprintf(fid,'\\thispagestyle{empty}\n');
0122 fprintf(fid,'\\clearpage \\newpage\n');
0123 fprintf(fid,'\\begin{figure}[t!] \\centering\n');
0124 fprintf(fid,'\\subfigure[]{\n');
0125 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'V.pdf']);
0126 fprintf(fid,'\\subfigure[]{\n');
0127 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'I1.pdf']);
0128 fprintf(fid,'\\subfigure[]{\n');
0129 fprintf(fid,'\\includegraphics[width=18.5cm]{%s}}\n',[dirname fprefix 'I2.pdf']);
0130 fprintf(fid,'\\end{figure}\n');
0131 fprintf(fid,'\\thispagestyle{empty}\\clearpage\n');
0132 fprintf(fid,'\\end{document}\n');
0133 fclose(fid);
0134 
0135 end