apm_files = glob(opts.apm_files) assert(len(apm_files)>0) print "found",len(apm_files),"apm files" #load the apm data positiontimes,positions,angletimes,angles,cmdtimes,CMDnums = read_apm_logs(apm_files) #load the receiver data rx_files = glob(opts.rx_files) assert(len(rx_files)>0) if rx_files[0].find('satpower')!=-1:#we're in orbcomm town pol = opts.pol.split('_')[1] rxtimes,freqs,rxspectrum = read_orbcomm_spectrum(rx_files,opts.rxant,pol) else: #read the default echo rx spectrum rxtimes,freqs,rxspectrum = read_echo_spectrum(rx_files) #get the power in our current channel rx_power = channel_select(freqs,rxspectrum,opts.freq) freqs = np.array([opts.freq]) #interpolate the rx data down to match the gps times rx_interp = interp_rx(positiontimes,rxtimes,rx_power) #flag based on yaws print "flagging based on yaw" yawmask,badyawtimes = flag_angles(angletimes,angles,2) posmask = apply_flagtimes(positiontimes,badyawtimes,1.0) rx_interp = np.ma.masked_where(posmask,rx_interp) print "total flags after yaw flagging:",np.sum(rx_interp.mask)
#! /usr/bin/env python from ECHO.read_utils import read_apm_log,apm_version,read_echo_spectrum import sys from pytz import timezone TZ = 'US/Arizona' tz = timezone(TZ) for filename in sys.argv[1:]: print filename apm_times,freqs,data =read_echo_spectrum([filename]) print "GPS start (UTC):",apm_times[0].iso print "GPS end (UTC):",apm_times[-1].iso print "GPS start ({z}):".format(z=TZ),apm_times[0].to_datetime(timezone=tz) print "GPS end ({z}):".format(z=TZ),apm_times[-1].to_datetime(timezone=tz) print "freqs (min,max) [MHz]:",freqs.min(),freqs.max()