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 rainbow_file = open (filenamestr+".tpdat","a")
70 #r = calib_default_total_power(data)
71 inter = globals()["calib_decln"]
72 rainbow_file.write(str(ephem.hours(sidtime))+" "+flt+" "+str(inter)+"\n")
76 def calib_numogate_ridge_observatory_fft(data,l):
81 # PyEphem wants lat/long as strings, rather than floats--took me quite
82 # a long time to figure that out. If they don't arrive as strings,
83 # the calculations for sidereal time are complete garbage
85 me.long = str(-76.043)
89 sidtime = me.sidereal_time()
91 foo = time.localtime()
93 if not "calib_prefix" in globals():
96 pfx = globals()["calib_prefix"]
97 filenamestr = "%s/%04d%02d%02d%02d" % (pfx, foo.tm_year,
98 foo.tm_mon, foo.tm_mday, foo.tm_hour)
102 if not "calib_then" in globals():
103 globals()["calib_then"] = now
107 if (now - globals()["calib_then"]) >= delta:
109 globals()["calib_then"] = now
110 rainbow_file = open (filenamestr+".sdat","a")
112 r = calib_default_fft(data,l)
113 inter = globals()["calib_decln"]
114 rainbow_file.write("data:"+str(ephem.hours(sidtime))+" "+str(inter)+" "+str(r)+"\n")
120 def calib_default_fft(db,l):
124 # We capture various parameters from the receive chain here, because
125 # they can affect the calibration equations.
128 def calib_set_gain(gain):
129 globals()["calib_gain_setting"] = gain
131 def calib_set_integ(integ):
132 globals()["calib_integ_setting"] = integ
134 def calib_set_bw(bw):
135 globals()["calib_bw_setting"] = bw
137 def calib_set_freq(freq):
138 globals()["calib_freq_setting"] = freq
140 def calib_set_avg_alpha(alpha):
141 globals()["calib_avg_alpha"] = alpha
143 def calib_set_interesting(inter):
144 globals()["calib_is_interesting"] = inter
146 def calib_set_decln(dec):
147 globals()["calib_decln"] = dec
149 def calib_set_prefix(pfx):
150 globals()["calib_prefix"] = pfx