15be23d5aee8c0f865a7ebb3115fce8857ff7acb
[debian/gnuradio] / gr-wxgui / src / python / scopesink_gl.py
1 #
2 # Copyright 2008,2010 Free Software Foundation, Inc.
3 #
4 # This file is part of GNU Radio
5 #
6 # GNU Radio is free software; you can redistribute it and/or modify
7 # it under the terms of the GNU General Public License as published by
8 # the Free Software Foundation; either version 3, or (at your option)
9 # any later version.
10 #
11 # GNU Radio is distributed in the hope that it will be useful,
12 # but WITHOUT ANY WARRANTY; without even the implied warranty of
13 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 # GNU General Public License for more details.
15 #
16 # You should have received a copy of the GNU General Public License
17 # along with GNU Radio; see the file COPYING.  If not, write to
18 # the Free Software Foundation, Inc., 51 Franklin Street,
19 # Boston, MA 02110-1301, USA.
20 #
21
22 ##################################################
23 # Imports
24 ##################################################
25 import scope_window
26 import common
27 from gnuradio import gr
28 from pubsub import pubsub
29 from constants import *
30 import math
31
32 class ac_couple_block(gr.hier_block2):
33         """
34         AC couple the incoming stream by subtracting out the low pass signal.
35         Mute the low pass filter to disable ac coupling.
36         """
37
38         def __init__(self, controller, ac_couple_key, sample_rate_key):
39                 gr.hier_block2.__init__(
40                         self,
41                         "ac_couple",
42                         gr.io_signature(1, 1, gr.sizeof_float),
43                         gr.io_signature(1, 1, gr.sizeof_float),
44                 )
45                 #blocks
46                 lpf = gr.single_pole_iir_filter_ff(0.0)
47                 sub = gr.sub_ff()
48                 mute = gr.mute_ff()
49                 #connect
50                 self.connect(self, sub, self)
51                 self.connect(self, lpf, mute, (sub, 1))
52                 #subscribe
53                 controller.subscribe(ac_couple_key, lambda x: mute.set_mute(not x))
54                 controller.subscribe(sample_rate_key, lambda x: lpf.set_taps(0.05))
55                 #initialize
56                 controller[ac_couple_key] = controller[ac_couple_key]
57                 controller[sample_rate_key] = controller[sample_rate_key]
58
59 ##################################################
60 # Scope sink block (wrapper for old wxgui)
61 ##################################################
62 class _scope_sink_base(gr.hier_block2, common.wxgui_hb):
63         """
64         A scope block with a gui window.
65         """
66
67         def __init__(
68                 self,
69                 parent,
70                 title='',
71                 sample_rate=1,
72                 size=scope_window.DEFAULT_WIN_SIZE,
73                 v_scale=0,
74                 t_scale=0,
75                 v_offset=0,
76                 xy_mode=False,
77                 ac_couple=False,
78                 num_inputs=1,
79                 trig_mode=scope_window.DEFAULT_TRIG_MODE,
80                 frame_rate=scope_window.DEFAULT_FRAME_RATE,
81                 use_persistence=False,
82                 persist_alpha=None,
83                 **kwargs #do not end with a comma
84         ):
85                 #ensure analog alpha
86                 if persist_alpha is None: 
87                   actual_frame_rate=float(frame_rate)
88                   analog_cutoff_freq=0.5 # Hertz
89                   #calculate alpha from wanted cutoff freq
90                   persist_alpha = 1.0 - math.exp(-2.0*math.pi*analog_cutoff_freq/actual_frame_rate)
91
92                 if not t_scale: t_scale = 10.0/sample_rate
93                 #init
94                 gr.hier_block2.__init__(
95                         self,
96                         "scope_sink",
97                         gr.io_signature(num_inputs, num_inputs, self._item_size),
98                         gr.io_signature(0, 0, 0),
99                 )
100                 #scope
101                 msgq = gr.msg_queue(2)
102                 scope = gr.oscope_sink_f(sample_rate, msgq)
103                 #controller
104                 self.controller = pubsub()
105                 self.controller.subscribe(SAMPLE_RATE_KEY, scope.set_sample_rate)
106                 self.controller.publish(SAMPLE_RATE_KEY, scope.sample_rate)
107                 self.controller.subscribe(DECIMATION_KEY, scope.set_decimation_count)
108                 self.controller.publish(DECIMATION_KEY, scope.get_decimation_count)
109                 self.controller.subscribe(TRIGGER_LEVEL_KEY, scope.set_trigger_level)
110                 self.controller.publish(TRIGGER_LEVEL_KEY, scope.get_trigger_level)
111                 self.controller.subscribe(TRIGGER_MODE_KEY, scope.set_trigger_mode)
112                 self.controller.publish(TRIGGER_MODE_KEY, scope.get_trigger_mode)
113                 self.controller.subscribe(TRIGGER_SLOPE_KEY, scope.set_trigger_slope)
114                 self.controller.publish(TRIGGER_SLOPE_KEY, scope.get_trigger_slope)
115                 self.controller.subscribe(TRIGGER_CHANNEL_KEY, scope.set_trigger_channel)
116                 self.controller.publish(TRIGGER_CHANNEL_KEY, scope.get_trigger_channel)
117                 actual_num_inputs = self._real and num_inputs or num_inputs*2
118                 #init ac couple
119                 for i in range(actual_num_inputs):
120                         self.controller[common.index_key(AC_COUPLE_KEY, i)] = ac_couple
121                 #start input watcher
122                 common.input_watcher(msgq, self.controller, MSG_KEY)
123                 #create window
124                 self.win = scope_window.scope_window(
125                         parent=parent,
126                         controller=self.controller,
127                         size=size,
128                         title=title,
129                         frame_rate=frame_rate,
130                         num_inputs=actual_num_inputs,
131                         sample_rate_key=SAMPLE_RATE_KEY,
132                         t_scale=t_scale,
133                         v_scale=v_scale,
134                         v_offset=v_offset,
135                         xy_mode=xy_mode,
136                         trig_mode=trig_mode,
137                         ac_couple_key=AC_COUPLE_KEY,
138                         trigger_level_key=TRIGGER_LEVEL_KEY,
139                         trigger_mode_key=TRIGGER_MODE_KEY,
140                         trigger_slope_key=TRIGGER_SLOPE_KEY,
141                         trigger_channel_key=TRIGGER_CHANNEL_KEY,
142                         decimation_key=DECIMATION_KEY,
143                         msg_key=MSG_KEY,
144                         use_persistence=use_persistence,
145                         persist_alpha=persist_alpha,
146                 )
147                 common.register_access_methods(self, self.win)
148                 #connect
149                 if self._real:
150                         for i in range(num_inputs):
151                                 self.wxgui_connect(
152                                         (self, i),
153                                         ac_couple_block(self.controller, common.index_key(AC_COUPLE_KEY, i), SAMPLE_RATE_KEY),
154                                         (scope, i),
155                                 )
156                 else:
157                         for i in range(num_inputs):
158                                 c2f = gr.complex_to_float() 
159                                 self.wxgui_connect((self, i), c2f)
160                                 for j in range(2):
161                                         self.connect(
162                                                 (c2f, j), 
163                                                 ac_couple_block(self.controller, common.index_key(AC_COUPLE_KEY, 2*i+j), SAMPLE_RATE_KEY),
164                                                 (scope, 2*i+j),
165                                         )
166
167 class scope_sink_f(_scope_sink_base):
168         _item_size = gr.sizeof_float
169         _real = True
170
171 class scope_sink_c(_scope_sink_base):
172         _item_size = gr.sizeof_gr_complex
173         _real = False
174
175 # ----------------------------------------------------------------
176 # Stand-alone test application
177 # ----------------------------------------------------------------
178
179 import wx
180 from gnuradio.wxgui import stdgui2
181
182 class test_top_block (stdgui2.std_top_block):
183     def __init__(self, frame, panel, vbox, argv):
184         stdgui2.std_top_block.__init__ (self, frame, panel, vbox, argv)
185
186         default_input_rate = 1e6
187         if len(argv) > 1:
188             input_rate = int(argv[1]) 
189         else:
190             input_rate = default_input_rate
191
192         if len(argv) > 2:
193             v_scale = float(argv[2])  # start up at this v_scale value
194         else:
195             v_scale = None  # start up in autorange mode, default
196
197         if len(argv) > 3:
198             t_scale = float(argv[3])  # start up at this t_scale value
199         else:
200             t_scale = .00003*default_input_rate/input_rate # old behavior
201
202         print "input rate %s  v_scale %s  t_scale %s" % (input_rate,v_scale,t_scale)
203             
204
205         # Generate a complex sinusoid
206         ampl=1.0e3
207         self.src0 = gr.sig_source_c (input_rate, gr.GR_SIN_WAVE, 25.1e3*input_rate/default_input_rate, ampl)
208         self.noise =gr.sig_source_c (input_rate, gr.GR_SIN_WAVE, 11.1*25.1e3*input_rate/default_input_rate, ampl/10) 
209         #self.noise =gr.noise_source_c(gr.GR_GAUSSIAN, ampl/10)
210         self.combine=gr.add_cc()
211
212         # We add this throttle block so that this demo doesn't suck down
213         # all the CPU available.  You normally wouldn't use it...
214         self.thr = gr.throttle(gr.sizeof_gr_complex, input_rate)
215
216         scope = scope_sink_c (panel,"Secret Data",sample_rate=input_rate,
217                               v_scale=v_scale, t_scale=t_scale)
218         vbox.Add (scope.win, 1, wx.EXPAND)
219
220         # Ultimately this will be
221         # self.connect("src0 throttle scope")
222         self.connect(self.src0,(self.combine,0))
223         self.connect(self.noise,(self.combine,1))
224         self.connect(self.combine, self.thr, scope) 
225
226 def main ():
227     app = stdgui2.stdapp (test_top_block, "O'Scope Test App")
228     app.MainLoop ()
229
230 if __name__ == '__main__':
231     main ()