meteor collection

quicklook.py

import numpy as np
from matplotlib.pylab import plt
plt.ion()

def to_cplx( r ):
    c = ( r[0::2]-127.5 ) + ( r[1::2]-127.5 )*1.0J
    return c

d = np.fromfile( 'file2019-11-17-21', dtype=np.uint8, count = 2**32 )

dc = to_cplx( d[0:2**16] )
fdc = np.fft.fft(dc)
plt.figure(); plt.plot( np.fft.fftshift( 10.0*np.log10(np.abs(fdc**2) ) ) )

nfft = 1024
pdc,freqs = plt.psd( dc,
                     NFFT = nfft,
                     Fs = 2.4e6,
                     Fc = 101.2e6,
                     window = np.hamming( nfft ),
                     noverlap = nfft/2 )

run.py

import schedule
import time
import datetime

import os

def job():
    print( datetime.datetime.now() )
    os.system( '/home/apn/proj/meteor/run.sh' )

#schedule.every().day.at('12:00').do( job )
#schedule.every().day.at('13:00').do( job )
#schedule.every().day.at('14:00').do( job )
    
#schedule.every().day.at('15:00').do( job )
#schedule.every().day.at('16:00').do( job )
#schedule.every().day.at('17:00').do( job )
#schedule.every().day.at('18:00').do( job )
#schedule.every().day.at('19:00').do( job )

#schedule.every().day.at('20:00').do( job )
schedule.every().day.at('21:00').do( job )
schedule.every().day.at('22:00').do( job )
schedule.every().day.at('23:00').do( job )
schedule.every().day.at('00:00').do( job )

schedule.every().day.at('01:00').do( job )
schedule.every().day.at('02:00').do( job )
schedule.every().day.at('03:00').do( job )
schedule.every().day.at('04:00').do( job )
schedule.every().day.at('05:00').do( job )

schedule.every().day.at('06:00').do( job )
#schedule.every().day.at('07:30').do( job )
#schedule.every().minute.do( job )
#schedule.every().hour.at(":00:").do(job)

while True:
    schedule.run_pending()
    time.sleep(1)

run.sh

#!/bin/sh

FN=`date +%F-%H`

rtl_sdr -f 101200000 -s 2400000 -g 496 -n 1440000000 "/home/apn/proj/meteor/file$FN"