lukesTrigPoint

 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=78
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  
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 pixel_size = 0.1; % degree
0098 
0099 azmin = -3;
0100 azmax = +3;
0101 elmin = -3;
0102 elmax = +3;
0103 nanflag = isnan(draster.antenna0.receiver.data(:,1));
0104 flag = ~nanflag;
0105 % Set colour scale mapping
0106 %clim = [0 3.5;-0.2 0.2;-0.2 0.2;0 0.2]
0107 
0108 %for m= 1:length(s)
0109  %   n = n+1;
0110  %   ind = zeros(size(FLAG));
0111  %   ind(s(m):e(m)) = 1;
0112  %   ind = logical(ind);
0113  %   flagvector3 = (~nanflag & ~FLAG & ind);
0114     
0115     for i=1:4
0116         [TauAmap{i}, xidx, yidx, TauAmapIdx{i}] = bin_quickly2d(azmin, azmax, elmin, elmax, pixel_size,...
0117             draster.azOffSave(flag), draster.elOffSave(flag),draster.antenna0.receiver.data(flag,i));
0118         figure(3)
0119         subplot(2,2,i)
0120         
0121         %imagesc(fliplr(xidx),yidx,TauAmap{i},clim(i,:))
0122         clim=[0 1*max(max(TauAmap{i}))]
0123         imagesc(fliplr(xidx),fliplr(yidx),TauAmap{i},clim)
0124         
0125         %title([source_name,' map:', num2str(i)])
0126         colormap('default')
0127         colorbar
0128         title(plot_labels{i})
0129         axis square
0130         figure(4)
0131         subplot(2,2,i)
0132         clines = [1:-0.1:0]*max(max(TauAmap{i}));
0133         contour(fliplr(xidx),yidx,TauAmap{i},clines)
0134         title(plot_labels{i})
0135         axis square
0136         
0137     end
0138 
0139 %%
0140 
0141 % Plot the elevation scan.
0142 % Make sure that we only plot data were azOff = 0 or -180
0143 az0 = (delscan.antenna0.servo.az<90);
0144 figure
0145 h(1)=subplot(3,1,1)
0146 plot(delscan.antenna0.servo.az(az0))
0147 title('Az')
0148 h(2)=subplot(3,1,2)
0149 plot(delscan.antenna0.servo.el(az0))
0150 title('El')
0151 h(3)=subplot(3,1,3)
0152 title('ll 30')
0153 plot(delscan.antenna0.receiver.data(az0,1))
0154 grid on
0155 linkaxes(h,'x')
0156