3 # Copyright 2003,2004,2005 Free Software Foundation, Inc.
5 # This file is part of GNU Radio
7 # GNU Radio is free software; you can redistribute it and/or modify
8 # it under the terms of the GNU General Public License as published by
9 # the Free Software Foundation; either version 2, or (at your option)
12 # GNU Radio is distributed in the hope that it will be useful,
13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 # GNU General Public License for more details.
17 # You should have received a copy of the GNU General Public License
18 # along with GNU Radio; see the file COPYING. If not, write to
19 # the Free Software Foundation, Inc., 51 Franklin Street,
20 # Boston, MA 02110-1301, USA.
29 # Simple class for allowing local definition of a calibration function
30 # for raw samples coming from the RA detector chain. Each observatory
31 # is different, and rather than hacking up the main code in usrp_ra_receiver
32 # we define the appropriate function here.
34 # For example, one could calibrate the output in Janskys, rather than
39 def calib_default_total_power(data):
40 r = 10.0*math.log10(data)
43 def calib_numogate_ridge_observatory_total_power(data):
48 # PyEphem wants lat/long as strings, rather than floats--took me quite
49 # a long time to figure that out. If they don't arrive as strings,
50 # the calculations for sidereal time are complete garbage
52 me.long = str(-76.043)
56 sidtime = me.sidereal_time()
58 foo = time.localtime()
59 if not "calib_prefix" in globals():
62 pfx = globals()["calib_prefix"]
63 filenamestr = "%s/%04d%02d%02d%02d" % (pfx, foo.tm_year,
64 foo.tm_mon, foo.tm_mday, foo.tm_hour)
66 numogate_file = open (filenamestr+".tpdat","a")
70 #r = calib_default_total_power(data)
71 inter = globals()["calib_decln"]
72 integ = globals()["calib_integ_setting"]
73 fc = globals()["calib_freq_setting"]
75 bw = globals()["calib_bw_setting"]
77 ga = globals()["calib_gain_setting"]
81 if not "calib_then_tpdat" in globals():
82 globals()["calib_then_tpdat"] = now
84 if (now - globals()["calib_then_tpdat"]) >= 20:
85 globals()["calib_then_tpdat"] = now
87 numogate_file.write(str(ephem.hours(sidtime))+" "+flt+" Dn="+str(inter)+",")
88 numogate_file.write("Ti="+str(integ)+",Fc="+str(fc)+",Bw="+str(bw))
89 numogate_file.write(",Ga="+str(ga)+"\n")
91 numogate_file.write(str(ephem.hours(sidtime))+" "+flt+"\n")
96 def calib_numogate_ridge_observatory_fft(data,l):
101 # PyEphem wants lat/long as strings, rather than floats--took me quite
102 # a long time to figure that out. If they don't arrive as strings,
103 # the calculations for sidereal time are complete garbage
105 me.long = str(-76.043)
108 me.date = ephem.now()
109 sidtime = me.sidereal_time()
111 foo = time.localtime()
113 if not "calib_prefix" in globals():
116 pfx = globals()["calib_prefix"]
117 filenamestr = "%s/%04d%02d%02d%02d" % (pfx, foo.tm_year,
118 foo.tm_mon, foo.tm_mday, foo.tm_hour)
122 if not "calib_then" in globals():
123 globals()["calib_then"] = now
127 if (now - globals()["calib_then"]) >= delta:
129 globals()["calib_then"] = now
130 numogate_file = open (filenamestr+".sdat","a")
132 r = calib_default_fft(data,l)
133 inter = globals()["calib_decln"]
134 fc = globals()["calib_freq_setting"]
136 bw = globals()["calib_bw_setting"]
138 av = globals()["calib_avg_alpha"]
139 numogate_file.write("data:"+str(ephem.hours(sidtime))+" Dn="+str(inter)+",Fc="+str(fc)+",Bw="+str(bw)+",Av="+str(av))
140 numogate_file.write(" "+str(r)+"\n")
141 numogate_file.close()
146 def calib_default_fft(db,l):
150 # We capture various parameters from the receive chain here, because
151 # they can affect the calibration equations.
154 def calib_set_gain(gain):
155 globals()["calib_gain_setting"] = gain
156 globals()["calib_then_tpdat"] = time.time() - 50
158 def calib_set_integ(integ):
159 globals()["calib_integ_setting"] = integ
160 globals()["calib_then_tpdat"] = time.time() - 50
162 def calib_set_bw(bw):
163 globals()["calib_bw_setting"] = bw
164 globals()["calib_then_tpdat"] = time.time() - 50
166 def calib_set_freq(freq):
167 globals()["calib_freq_setting"] = freq
168 globals()["calib_then_tpdat"] = time.time() - 50
170 def calib_set_avg_alpha(alpha):
171 globals()["calib_avg_alpha"] = alpha
173 def calib_set_interesting(inter):
174 globals()["calib_is_interesting"] = inter
176 def calib_set_decln(dec):
177 globals()["calib_decln"] = dec
178 globals()["calib_then_tpdat"] = time.time() - 50
180 def calib_set_prefix(pfx):
181 globals()["calib_prefix"] = pfx