makePlots

0001 [s w] = unix('ls --color=never -lh /mnt/data/sjcm/jup_sec');
0002 
0003 entVals = strfind(w, '.fit');
0004 
0005 for m=1:length(entVals)
0006   entry{m} = w(entVals(m)-15:entVals(m)+3);
0007 end
0008 
0009 %eloff = 0.75:0.25:4.76;
0010 %azoff = 0.15:0.25:4.16;
0011 eloff = -2:0.15:2;
0012 azoff = -2:0.15:2;
0013 
0014 keyboard;
0015 index = 1;
0016 fullMap = [];
0017 for m=1:length(eloff)
0018   thisRow = [];
0019   for n=1:length(azoff)
0020     eval(sprintf('thisFits = fitsread2(''/mnt/data/sjcm/jup_sec/%s'');', entry{index}));
0021 %thisFits(thisFits<0) = nan;
0022 %thisFits(thisFits>50) = nan;
0023     tf(:,:,n) = thisFits;
0024     index  = index+1;
0025   end
0026     if m==16
0027     display('at right row');
0028 keyboard;
0029 end
0030 
0031     
0032   badPix = mean(tf,3)>100 | std(tf,0,3)==0;
0033 flat = ones(size(thisFits));
0034 flat(badPix) = nan;
0035 flat = repmat(flat, [1 1 n]);
0036 tf = tf./flat;
0037 thisRow = [];
0038 for n=1:length(azoff)
0039   thisRow = [thisRow tf(:,:,n)];
0040 end
0041 imagesc(thisRow);
0042 eval(sprintf('title(''Elevation Row %d: %2.2f'');', m, eloff(m)));
0043 meanRowAz(m) = mean(tf(:));
0044 colorbar;
0045 
0046 pause(1);
0047 end
0048 
0049 keyboard;