plot_trigpoint_files

 21/2/2014 - ACT
 raster and large az and el tracks on the sun

0001 % 21/2/2014 - ACT
0002 % raster and large az and el tracks on the sun
0003 
0004 %d = read_arcSouth('21-Feb-2014:13:23:03 ','     21-Feb-2014:13:57:26 ')
0005 
0006 d = read_arcSouth('23-Feb-2014:13:20:20','23-Feb-2014:13:54:44')
0007 
0008 % Horn at 45 deg, 0dBm-30dB attenuation
0009 close all;
0010 clear;
0011 %d = read_arcSouth('23-Feb-2014:15:09','23-Feb-2014:15:30')
0012 d = read_arcSouth('23-Feb-2014:15:07:51','23-Feb-2014:15:24:02') 
0013 
0014 %Horn at 0deg, 0dbm-30dB attenuation
0015 close all;clear;
0016 d = read_arcSouth('23-Feb-2014:15:27','23-Feb-2014:15:40')
0017 
0018 %Horn at 0deg, 0dbm-60dB attenuation
0019 close all;clear;
0020 d = read_arcSouth('23-Feb-2014:16:05','23-Feb-2014:16:15')
0021 
0022 
0023 %%
0024 % Put data into Cbass-N style
0025 trigpoint=1
0026 
0027 if (trigpoint==1)
0028     chan=77
0029     d = cbassS2N(d,chan);
0030     
0031 else
0032     d = cbassS2N(d);
0033 end
0034 
0035 % nchan=4;
0036 % LL = [d.antenna0.roach1.LL,fliplr(d.antenna0.roach2.LL)]; % **Check which order the roaches are in **
0037 % RR = [d.antenna0.roach1.RR,fliplr(d.antenna0.roach2.RR)]; % **Check which order the roaches are in **
0038 % QQ = [d.antenna0.roach1.Q,fliplr(d.antenna0.roach2.Q)];
0039 % UU = [d.antenna0.roach1.U,fliplr(d.antenna0.roach2.U)];
0040 % plot_labels = {'LL','RR','QQ','UU'};
0041 %
0042 % % Make a data structure d.antenna0.recevier.data(:,1:4)  LL,RR,QQ,UU
0043 % % averaged over frequency
0044 % % saner piece of code so don't have to do this! Relic of history...
0045 %  d.antenna0.receiver.data(:,1)=  LL(:,78); %LL
0046 %  d.antenna0.receiver.data(:,2)=  RR(:,78); %RR - channel 78 only
0047 %  d.antenna0.receiver.data(:,3)=  QQ(:,78); %QQ
0048 %  d.antenna0.receiver.data(:,4)=  UU(:,78); %UU
0049 %  d.antenna0.receiver.utc = d.antenna0.roach1.utc;
0050  %%
0051  
0052  nchan=4
0053  plot_labels = {'LL','RR','QQ','UU'};
0054  chan_order=[1 6 2 3];
0055 %%
0056 if (trigpoint==1)
0057     
0058     % Calculate the az and El offset form the source
0059     %az_trig = 17-1
0060     %el_trig = 11.5-4.5
0061     az_trig=16
0062     el_trig=6.5
0063     
0064     
0065     d.azOffSave = (d.antenna0.servo.apparent(:,1)-az_trig).*cosd(d.antenna0.servo.apparent(:,2));
0066     d.elOffSave = d.antenna0.servo.apparent(:,2)-el_trig;
0067     d.angle = sqrt(d.elOffSave.^2 + d.azOffSave.^2);
0068     
0069 else
0070     azApp = interp1(d.antenna0.tracker.utc, ...
0071         d.antenna0.tracker.horiz_topo(:,1), d.antenna0.roach2.utc);
0072     azOffSave = azApp - d.antenna0.servo.apparent(:,1);
0073     azOffSave = -azOffSave;
0074     azOffSave = azOffSave.*cos(d.antenna0.servo.el*pi/180);% Scale azimuth by cos(elevation)???CHECK
0075     
0076     elApp = interp1(d.antenna0.tracker.utc, ...
0077         d.antenna0.tracker.horiz_topo(:,2), d.antenna0.roach2.utc);
0078     elOffSave = elApp - d.antenna0.servo.apparent(:,2);
0079     elOffSave = -elOffSave;
0080     d.angle = sqrt(elOffSave.^2 + azOffSave.^2);%
0081     d.elOffSave = elOffSave;
0082     d.azOffSave= azOffSave;
0083     
0084 end
0085 %%
0086  % Separate data in to raster, elscan and azscan
0087 
0088 draster = framecut(d,d.index.beammap.fast);
0089 dazscan = framecut(d,d.index.radio_point_scan.fast);
0090 delscan = framecut(d,d.index.elscan.fast);
0091 
0092 
0093 
0094 %%
0095  
0096 % Now analyse the Raster data
0097 
0098 pixel_size = 0.1; % degree
0099 
0100 azmin = -3;
0101 azmax = +3;
0102 elmin = -3;
0103 elmax = +3;
0104 nanflag = isnan(draster.antenna0.receiver.data(:,1));
0105 flag = ~nanflag;
0106 
0107 % Set colour scale mapping
0108 %clim = [0 3.5;-0.2 0.2;-0.2 0.2;0 0.2]
0109 
0110 %for m= 1:length(s)
0111  %   n = n+1;
0112  %   ind = zeros(size(FLAG));
0113  %   ind(s(m):e(m)) = 1;
0114  %   ind = logical(ind);
0115  %   flagvector3 = (~nanflag & ~FLAG & ind);
0116    
0117     for i=1:4
0118         j=chan_order(i);
0119         [TauAmap{i}, xidx, yidx, TauAmapIdx{i}] = bin_quickly2d(azmin, azmax, elmin, elmax, pixel_size,...
0120             draster.azOffSave(flag), draster.elOffSave(flag),draster.antenna0.receiver.data(flag,j));
0121         figure(3)
0122         subplot(2,2,i)
0123         
0124         %imagesc(fliplr(xidx),yidx,TauAmap{i},clim(i,:))
0125         %clim=[5e3 1.5e4]
0126        imagesc(fliplr(xidx),fliplr(yidx),TauAmap{i})
0127         %imagesc(fliplr(xidx),fliplr(yidx),TauAmap{i},clim)
0128         
0129         %title([source_name,' map:', num2str(i)])
0130         colormap('default')
0131         colorbar
0132         title(plot_labels{i})
0133         axis square
0134         figure(4)
0135         subplot(2,2,i)
0136         %clines = [1:-0.1:0]*max(max(TauAmap{i}));
0137         contour(fliplr(xidx),yidx,TauAmap{i})
0138 %        contour(fliplr(xidx),yidx,TauAmap{i},clines)
0139         title(plot_labels{i})
0140         axis square
0141         
0142     end
0143 
0144 %%
0145 
0146 % Plot the elevation scan.
0147 % Make sure that we only plot data were azOff = 0 or -180
0148 az0 = (delscan.antenna0.servo.az<90);
0149 figure
0150 h(1)=subplot(3,1,1)
0151 plot(delscan.antenna0.servo.az(az0))
0152 title('Az')
0153 h(2)=subplot(3,1,2)
0154 plot(delscan.antenna0.servo.el(az0))
0155 title('El')
0156 h(3)=subplot(3,1,3)
0157 title('ll 30')
0158 plot(delscan.antenna0.receiver.data(az0,1))
0159 grid on
0160 linkaxes(h,'x')
0161