checkImages

0001 function d = checkImages(filename)
0002 
0003 d = [];
0004 
0005 txt= sprintf('[s w] = unix(''ls --color=never -lh /mnt/data/cbass/fits/%s*'');', filename);
0006 
0007 eval(txt);
0008 
0009 
0010 entVals = strfind(w, '.fit');
0011 
0012 for m=1:length(entVals)
0013   entry{m} = w(entVals(m)-15:entVals(m)+3);
0014 end
0015 
0016 
0017 for m=1:length(entry)
0018   eval(sprintf('thisFits = fitsread2(''/mnt/data/cbass/fits/%s'');', entry{m}));
0019 imagesc(thisFits);
0020 
0021 eval(sprintf('title(''%s'');', entry{m}(10:15)));
0022 pause(0.1);
0023 
0024 end
0025 
0026 
0027 % now we read in the data
0028 %d = read_arc('11-feb-2010:04:43:31', '11-feb-2010:11:34:00');
0029 %badTimes = [date2mjd(2010, 2, 11, 5, 55) date2mjd(2010, 2, 11, 6, 05);
0030 %  date2mjd(2010, 2, 11, 6, 10) date2mjd(2010, 2, 11, 6, 17);
0031 %  date2mjd(2010, 2, 11, 6, 18) date2mjd(2010, 2, 11, 6, 20);
0032 %  date2mjd(2010, 2, 11, 6, 59) date2mjd(2010, 2, 11, 7, 01);
0033 %  date2mjd(2010, 2, 11, 7, 36) date2mjd(2010, 2, 11, 7, 37);
0034 %  date2mjd(2010, 2, 11, 8, 04) date2mjd(2010, 2, 11, 8, 08);
0035 %  date2mjd(2010, 2, 11, 8, 36) date2mjd(2010, 2, 11, 8, 41);
0036 %  date2mjd(2010, 2, 11, 9, 02) date2mjd(2010, 2, 11, 9, 10);
0037 %  date2mjd(2010, 2, 11, 9, 30) date2mjd(2010, 2, 11, 9, 36, 30);
0038 %  date2mjd(2010, 2, 11, 9, 54) date2mjd(2010, 2, 11, 9, 57);
0039 %  date2mjd(2010, 2, 11,10, 15) date2mjd(2010, 2, 11, 10,19);
0040 %  date2mjd(2010, 2, 11,10, 34) date2mjd(2010, 2, 11, 10, 38);
0041 %  date2mjd(2010, 2, 11,10, 51) date2mjd(2010, 2, 11, 11, 01);
0042 %  date2mjd(2010, 2, 11,11, 06) date2mjd(2010, 2, 11, 11, 10);
0043 %  date2mjd(2010, 2, 11,11, 16) date2mjd(2010, 2, 11, 11, 22);
0044 %  date2mjd(2010, 2, 11,11, 24) date2mjd(2010, 2, 11, 11, 27);
0045 %  date2mjd(2010, 2, 11,11, 30) date2mjd(2010, 2, 11, 11, 33);
0046 %  ];
0047 %ind = ones(size(d.array.frame.utc));
0048 %for m=1:size(badTimes,1)
0049 %  f = find(d.array.frame.utc>badTimes(m,1) & ...
0050 %      d.array.frame.utc<badTimes(m,2));
0051 %  ind(f) = 0;
0052 %end
0053 %d = framecut(d, logical(ind));
0054 
0055     
0056   
0057   
0058 
0059 
0060 return
0061 
0062 % data Apr 26, 2011
0063 d = read_arc('26-Apr-2011:04:42:54',  '26-Apr-2011:12:51:00');
0064 % data to remove
0065 badTimes = [date2mjd(2011, 04, 26, 10, 25) date2mjd(2011, 04, 26, 10, 29);
0066   date2mjd(2011, 04, 26, 10, 36) date2mjd(2011, 04, 26, 10, 40);
0067   date2mjd(2011, 04, 26, 11, 52) date2mjd(2011, 04, 26, 12, 31);];
0068 ind = ones(size(d.array.frame.utc));
0069 for m=1:size(badTimes,1)
0070   f = find(d.array.frame.utc>badTimes(m,1) & ...
0071       d.array.frame.utc<badTimes(m,2));
0072   ind(f) = 0;
0073 end
0074 d = framecut(d, logical(ind));
0075 dc = framecut(d, d.array.frame.features>0);
0076 save opt_point_apr26.mat dc
0077 
0078 
0079 
0080 
0081 % data from October 14, 2010
0082 d = read_arc('14-oct-2010:05:05:54', '14-oct-2010:13:38:00');
0083 
0084 % data to remove
0085 badTimes = [date2mjd(2010, 10, 14, 3, 33) date2mjd(2010, 10, 14, 3, 37);
0086   date2mjd(2010, 10, 14, 4, 20) date2mjd(2010, 10, 14, 4, 23);
0087   date2mjd(2010, 10, 14, 5, 35) date2mjd(2010, 10, 14, 5, 38);
0088   date2mjd(2010, 10, 14, 5, 42) date2mjd(2010, 10, 14, 5, 45);
0089   date2mjd(2010, 10, 14, 6, 40) date2mjd(2010, 10, 14, 6, 42);
0090   date2mjd(2010, 10, 14, 7, 31) date2mjd(2010, 10, 14, 7, 35);
0091   date2mjd(2010, 10, 14, 9, 30) date2mjd(2010, 10, 14, 9, 35);
0092   date2mjd(2010, 10, 14,10, 11) date2mjd(2010, 10, 14,10, 16);
0093   date2mjd(2010, 10, 14,11, 03) date2mjd(2010, 10, 14,11, 18);
0094   date2mjd(2010, 10, 14,11, 59) date2mjd(2010, 10, 14,12, 14);
0095   date2mjd(2010, 10, 14,12, 45) date2mjd(2010, 10, 14, 12,54);
0096   date2mjd(2010, 10, 14,13, 14) date2mjd(2010, 10, 14, 13,16);
0097   date2mjd(2010, 10, 14,13, 21) date2mjd(2010, 10, 14, 13, 24);
0098   ];
0099 ind = ones(size(d.array.frame.utc));
0100 for m=1:size(badTimes,1)
0101   f = find(d.array.frame.utc>badTimes(m,1) & ...
0102       d.array.frame.utc<badTimes(m,2));
0103   ind(f) = 0;
0104 end
0105 d = framecut(d, logical(ind));
0106 dc = framecut(d, d.array.frame.features>0);
0107 save point_oct14.mat dc;
0108 
0109 %------------------------------------------------------------
0110 
0111 d = read_arc('15-Oct-2010:02:54:02', '15-oct-2010:10:00');
0112 
0113 badTimes = [date2mjd(2010, 10, 15,10, 50) date2mjd(2010, 10, 15,10, 53);
0114   date2mjd(2010, 10, 15,10, 57) date2mjd(2010, 10, 15,11, 00)];
0115 ind = ones(size(d.array.frame.utc));
0116 for m=1:size(badTimes,1)
0117   f = find(d.array.frame.utc>badTimes(m,1) & ...
0118       d.array.frame.utc<badTimes(m,2));
0119   ind(f) = 0;
0120 end
0121 d = framecut(d, logical(ind));
0122 dc = framecut(d, d.array.frame.features>0);
0123 
0124 d = read_arc('15-Oct-2010:10:00:00', '15-oct-2010:13:39');
0125 badTimes = [date2mjd(2010, 10, 15,10, 50) date2mjd(2010, 10, 15,10, 53);
0126   date2mjd(2010, 10, 15,10, 57) date2mjd(2010, 10, 15,11, 00)];
0127 ind = ones(size(d.array.frame.utc));
0128 for m=1:size(badTimes,1)
0129   f = find(d.array.frame.utc>badTimes(m,1) & ...
0130       d.array.frame.utc<badTimes(m,2));
0131   ind(f) = 0;
0132 end
0133 d = framecut(d, logical(ind));
0134 dc2 = framecut(d, d.array.frame.features>0);
0135 
0136 dc = catstruct(1, [dc dc2]);
0137 
0138 ind = diff(dc.array.frame.utc)*24*60*60>2;
0139 ind = [1; ind];
0140 
0141 dc = framecut(dc, ind);
0142 
0143 save point_oct15.mat dc;
0144 
0145 
0146 %------------------------------------------------------------
0147 d1 = read_arc('25-Nov-2010:01:17:51', '25-nov-2010:04:27:47');
0148 d2 = read_arc('25-Nov-2010:05:55:00', '25-nov-2010:09:12:02');
0149 d = catstruct(1, [d1 d2]);
0150 clear d1; 
0151 clear d2;
0152 badTimes = [date2mjd(2010, 11, 25,03, 30) date2mjd(2010, 11, 25,03, 35);
0153   date2mjd(2010, 11, 25, 7, 45) date2mjd(2010, 11, 25, 7, 52)];
0154 ind = ones(size(d.array.frame.utc));
0155 for m=1:size(badTimes,1)
0156   f = find(d.array.frame.utc>badTimes(m,1) & ...
0157       d.array.frame.utc<badTimes(m,2));
0158   ind(f) = 0;
0159 end
0160 d = framecut(d, logical(ind));
0161 dc = framecut(d, d.array.frame.features>0);
0162 
0163 ind = diff(dc.array.frame.utc)*24*60*60>2;
0164 ind = [1; ind];
0165 
0166 dc = framecut(dc, logical(ind));
0167 
0168 save point_nov25.mat dc;
0169 
0170 
0171 %------------------------------------------------------------
0172 
0173 d = read_arc('30-May-2011:06:01:10',  '30-May-2011:12:02:00');
0174 badTimes = [...
0175   date2mjd(2011, 05, 30, 7, 31) date2mjd(2011, 05, 30, 7, 40);
0176   date2mjd(2011, 05, 30, 7, 43) date2mjd(2011, 05, 30, 7, 48);
0177   date2mjd(2011, 05, 30, 7, 52) date2mjd(2011, 05, 30, 8, 06);
0178   date2mjd(2011, 05, 30, 8, 21) date2mjd(2011, 05, 30, 8, 45);
0179   date2mjd(2011, 05, 30, 8, 50) date2mjd(2011, 05, 30, 8, 59);
0180   date2mjd(2011, 05, 30, 9, 08) date2mjd(2011, 05, 30, 9, 15);
0181   date2mjd(2011, 05, 30, 9, 25) date2mjd(2011, 05, 30, 9, 34);
0182   date2mjd(2011, 05, 30, 9, 39) date2mjd(2011, 05, 30, 10,03);
0183   date2mjd(2011, 05, 30,10, 08) date2mjd(2011, 05, 30, 10,13);
0184   date2mjd(2011, 05, 30,10, 32) date2mjd(2011, 05, 30, 10,32) ; 
0185   date2mjd(2011, 05, 30,10, 50) date2mjd(2011, 05, 30, 10,57);
0186   date2mjd(2011, 05, 30,11, 03) date2mjd(2011, 05, 30, 11, 26);
0187   date2mjd(2011, 05, 30,11, 38) date2mjd(2011, 05, 30, 11, 48)];
0188 
0189 ind = ones(size(d.array.frame.utc));
0190 for m=1:size(badTimes,1)
0191   f = find(d.array.frame.utc>badTimes(m,1) & ...
0192       d.array.frame.utc<badTimes(m,2));
0193   ind(f) = 0;
0194 end
0195 d = framecut(d, logical(ind));
0196 dc = framecut(d, d.array.frame.features>0);
0197 
0198 ind = diff(dc.array.frame.utc)*24*60*60>2;
0199 ind = [1; ind];
0200 
0201 dc = framecut(dc, logical(ind));
0202 
0203 ind = dc.antenna0.tracker.horiz_off(:,2)<-0.05;
0204 dc = framecut(dc, ~ind);
0205 
0206 save point_may30.mat dc;
0207 
0208 %------------------------------------------------------------
0209 
0210 d = read_arc('28-May-2011:05:11:17' , '28-May-2011:11:23:58');
0211 dc = framecut(d, d.array.frame.features>0);
0212 
0213 ind = diff(dc.array.frame.utc)*24*60*60>2;
0214 ind = [1; ind];
0215 
0216 dc = framecut(dc, logical(ind));
0217 
0218 save point_may28_lst.mat dc;
0219 
0220 %------------------------------------------------------------
0221 
0222 d = read_arc('29-May-2011:05:01:30',  '29-May-2011:12:05:44');
0223 
0224 badTimes = [...
0225   date2mjd(2011, 05, 29, 5, 28) date2mjd(2011, 05, 29, 5, 31);
0226   date2mjd(2011, 05, 29, 5, 34) date2mjd(2011, 05, 29, 5, 40);
0227   date2mjd(2011, 05, 29, 6, 17) date2mjd(2011, 05, 29, 6, 19);
0228   date2mjd(2011, 05, 29, 6, 29) date2mjd(2011, 05, 29, 6, 38);
0229   date2mjd(2011, 05, 29, 6, 43) date2mjd(2011, 05, 29, 6, 58);
0230   date2mjd(2011, 05, 29, 7, 09) date2mjd(2011, 05, 29, 7, 11);
0231   date2mjd(2011, 05, 29, 7, 16) date2mjd(2011, 05, 29, 7, 32);
0232   date2mjd(2011, 05, 29, 7, 37) date2mjd(2011, 05, 29, 7, 50);
0233   date2mjd(2011, 05, 29, 7, 56) date2mjd(2011, 05, 29,11, 50);]
0234 
0235 ind = ones(size(d.array.frame.utc));
0236 for m=1:size(badTimes,1)
0237   f = find(d.array.frame.utc>badTimes(m,1) & ...
0238       d.array.frame.utc<badTimes(m,2));
0239   ind(f) = 0;
0240 end
0241 d = framecut(d, logical(ind));
0242 dc = framecut(d, d.array.frame.features>0);
0243 ind = diff(dc.array.frame.utc)*24*60*60>2;
0244 ind = [1; ind];
0245 
0246 dc = framecut(dc, logical(ind));
0247 
0248 save point_may29.mat dc;
0249 
0250 %------------------------------------------------------------
0251 regs={'array.frame.received'...
0252     'array.frame.utc double',...
0253     'array.frame.features',...
0254     'array.weather.utc double',...
0255     'array.weather.pressure double',...
0256     'array.weather.airTemperature double',...
0257     'array.weather.windSpeed[0] double',...
0258     'array.weather.windDirection[0] double',...
0259     'array.weather.status double',...
0260     'array.weather.relativeHumidity double',...
0261     'antenna0.thermal.ccHeaterCurrent double',...
0262     'antenna0.thermal.lsTemperatureSensors double',...
0263     'antenna0.thermal.ccTemperatureLoad double',...
0264     'antenna0.thermal.utc double',...
0265     'antenna0.thermal.dlpTemperatureSensors double',...
0266     'antenna0.tracker.lst double',...
0267     'antenna0.tracker.lacking double',...
0268     'antenna0.tracker.equat_geoc double',...
0269     'antenna0.tracker.horiz_topo double',...
0270     'antenna0.tracker.horiz_mount double',...
0271     'antenna0.tracker.flexure double',...
0272     'antenna0.tracker.horiz_off double',...
0273     'antenna0.tracker.tilts double',...
0274     'antenna0.tracker.fixedCollimation double',...
0275     'antenna0.tracker.encoder_off double',...
0276     'antenna0.tracker.sky_xy_off double',...
0277     'antenna0.tracker.source string',...
0278     'antenna0.tracker.scan_off double',...
0279     'antenna0.tracker.refraction double',...
0280     'antenna0.tracker.ut1utc double',...
0281     'antenna0.tracker.eqneqx double',...
0282     'antenna0.tracker.utc double',...
0283     'antenna0.tracker.time_diff double',...
0284     'antenna0.tracker.offSource double',...
0285     'antenna0.tracker.siteActual double',...
0286     'antenna0.tracker.siteFiducial double',...
0287     'antenna0.servo.utc double',...
0288     'antenna0.servo.fast_az_pos double',...
0289     'antenna0.servo.fast_el_pos double',...
0290     'antenna0.servo.fast_az_err double',...
0291     'antenna0.servo.fast_el_err double',...
0292     'antenna0.servo.actual_current_az1[0] double',...
0293     'antenna0.servo.actual_current_az2[0] double',...
0294     'antenna0.servo.actual_current_el1[0] double',...
0295     'antenna0.frame.utc double',...
0296     'antenna0.frame.received double',...
0297     'antenna0.receiver.flags double',...
0298     'antenna0.receiver.utc double',...
0299     'antenna0.receiver.data double',...
0300     'antenna0.receiver.diagnostics double',...
0301     'antenna0.receiver.drainCurrent double',...
0302     'antenna0.receiver.drainVoltage double',...
0303     'antenna0.receiver.gateVoltage double',...
0304     };
0305                 
0306 d = read_arc('28-oct-2011:02:15:23',  '28-Oct-2011:11:51:00', regs);
0307 
0308 badTimes = [...
0309   date2mjd(2011, 10, 28, 4, 05) date2mjd(2011, 10, 28, 4, 08);
0310   date2mjd(2011, 10, 28, 7, 03) date2mjd(2011, 10, 28, 7, 07);
0311   date2mjd(2011, 10, 28, 7, 10) date2mjd(2011, 10, 28, 7, 14);
0312   date2mjd(2011, 10, 28, 8, 40) date2mjd(2011, 10, 28, 8, 49, 30);
0313   date2mjd(2011, 10, 28, 8, 55) date2mjd(2011, 10, 28, 9, 32);
0314   date2mjd(2011, 10, 28, 9, 46) date2mjd(2011, 10, 28,10, 36);
0315   date2mjd(2011, 10, 28,11, 08) date2mjd(2011, 10, 28,11, 11);
0316   date2mjd(2011, 10, 28,11, 31) date2mjd(2011, 10, 28,11, 36);
0317   date2mjd(2011, 10, 28,11, 39) date2mjd(2011, 10, 28,11, 48);]
0318 
0319 ind = ones(size(d.array.frame.utc));
0320 for m=1:size(badTimes,1)
0321   f = find(d.array.frame.utc>badTimes(m,1) & ...
0322       d.array.frame.utc<badTimes(m,2));
0323   ind(f) = 0;
0324 end
0325 d = framecut(d, logical(ind));
0326 dc = framecut(d, d.array.frame.features>0);
0327 ind = diff(dc.array.frame.utc)*24*60*60>2;
0328 ind = [1; ind];
0329 
0330 dc = framecut(dc, logical(ind));
0331 
0332 dc1 = framecut(dc, dc.array.frame.features>0);
0333 
0334 save point_20111028.mat dc;
0335 
0336 
0337 d = read_arc('09-Feb-2012:02:51:45', '09-Feb-2012:12:32:00', regs);
0338 
0339 badTimes = [...
0340   date2mjd(2012, 02, 09, 5, 43) date2mjd(2012, 02, 09, 5, 47);
0341   date2mjd(2012, 02, 09, 6, 11) date2mjd(2012, 02, 09, 6, 23);
0342   date2mjd(2012, 02, 09, 6, 31) date2mjd(2012, 02, 09, 6, 39);
0343   date2mjd(2012, 02, 09, 6, 42) date2mjd(2012, 02, 09, 6, 45, 30);
0344   date2mjd(2012, 02, 09, 6, 50) date2mjd(2012, 02, 09, 7, 03);
0345   date2mjd(2012, 02, 09, 9, 15) date2mjd(2012, 02, 09, 9, 19);
0346   date2mjd(2012, 02, 09, 9, 25) date2mjd(2012, 02, 09, 9, 29);
0347   date2mjd(2012, 02, 09, 9, 32) date2mjd(2012, 02, 09, 9, 41);
0348   date2mjd(2012, 02, 09,10, 32) date2mjd(2012, 02, 09,10, 37);
0349   date2mjd(2012, 02, 09,12, 24) date2mjd(2012, 02, 09,12, 29);]
0350 
0351 ind = ones(size(d.array.frame.utc));
0352 for m=1:size(badTimes,1)
0353   f = find(d.array.frame.utc>badTimes(m,1) & ...
0354       d.array.frame.utc<badTimes(m,2));
0355   ind(f) = 0;
0356 end
0357 d = framecut(d, logical(ind));
0358 dc = framecut(d, d.array.frame.features>0);
0359 ind = diff(dc.array.frame.utc)*24*60*60>2;
0360 ind = [1; ind];
0361 
0362 dc = framecut(dc, logical(ind));
0363 
0364 dc1 = framecut(dc, dc.array.frame.features>0);
0365 
0366 save point_20120209.mat dc1;
0367