[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Commit-gnuradio] r4447 - gnuradio/branches/developers/nldudok1/general-
From: |
nldudok1 |
Subject: |
[Commit-gnuradio] r4447 - gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp |
Date: |
Sun, 11 Feb 2007 16:30:27 -0700 (MST) |
Author: nldudok1
Date: 2007-02-11 16:30:27 -0700 (Sun, 11 Feb 2007)
New Revision: 4447
Added:
gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp/usrp_fmtv_rcv.py
Log:
added fmtv sample (no sync color or audio yet). Works on ASTRA satellite, not
tested with amateur TV yet.
Added:
gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp/usrp_fmtv_rcv.py
===================================================================
---
gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp/usrp_fmtv_rcv.py
(rev 0)
+++
gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp/usrp_fmtv_rcv.py
2007-02-11 23:30:27 UTC (rev 4447)
@@ -0,0 +1,438 @@
+#!/usr/bin/env python
+"""
+Realtime capture and display of analog Tv stations.
+Can also use a file as source or sink
+When you use an output file you can show the results frame-by-frame using
ImageMagick
+When you want to use the realtime sdl display window you must first install
gr-video-sdl (is in gnuradio cvs).
+When you use a file source, in stead of the usrp, make sure you capture
interleaved shorts.
+(Use usrp_rx_file.py, or use usrp_rx_cfile.py --output-shorts if you have a
recent enough usrp_rx_cfile.py)
+There is no synchronisation yet. The sync blocks are in development but not
yet in cvs.
+
+"""
+from gnuradio import gr, gru, eng_notation, optfir
+try:
+ from gnuradio import video_sdl
+except:
+ print "FYI: gr-video-sdl is not installed"
+ print "realtime SDL video output window will not be available"
+from gnuradio import usrp
+from gnuradio import blks
+from gnuradio.eng_option import eng_option
+from gnuradio.wxgui import slider, powermate
+from gnuradio.wxgui import stdgui, fftsink, form
+from optparse import OptionParser
+import usrp_dbid
+import sys
+import math
+import wx
+import math
+
+# To debug, insert this in your test code...
+#import os
+#print 'Blocked waiting for GDB attach (pid = %d)' % (os.getpid(),)
+#raw_input ('Press Enter to continue: ')
+# remainder of your test code follows...
+
+def pick_subdevice(u):
+ """
+ The user didn't specify a subdevice on the command line.
+ Try for one of these, in order: TV_RX, BASIC_RX, whatever is on side A.
+
+ @return a subdev_spec
+ """
+ return usrp.pick_subdev(u, (usrp_dbid.TV_RX,
+ usrp_dbid.TV_RX_REV_2,
+ usrp_dbid.BASIC_RX))
+
+
+class tv_rx_graph (stdgui.gui_flow_graph):
+ def __init__(self,frame,panel,vbox,argv):
+ stdgui.gui_flow_graph.__init__ (self,frame,panel,vbox,argv)
+
+ usage="%prog: [options] [input_filename]. \n If you don't specify an
input filename the usrp will be used as source\n " \
+ "Make sure your input capture file containes interleaved shorts
not complex floats"
+ parser=OptionParser(option_class=eng_option)
+ parser.add_option("-R", "--rx-subdev-spec", type="subdev",
default=None,
+ help="select USRP Rx side A or B (default=A)")
+ parser.add_option("-d", "--decim", type="int", default=64,
+ help="set fgpa decimation rate to DECIM
[default=%default]")
+ parser.add_option("-f", "--freq", type="eng_float", default=519.25e6,
+ help="set frequency to FREQ", metavar="FREQ")
+ parser.add_option("-g", "--gain", type="eng_float", default=None,
+ help="set gain in dB (default is midpoint)")
+ parser.add_option("-c", "--contrast", type="eng_float", default=1.0,
+ help="set contrast (default is 1.0)")
+ parser.add_option("-b", "--brightness", type="eng_float", default=0.0,
+ help="set brightness (default is 0)")
+ parser.add_option("-8", "--width-8", action="store_true",
default=False,
+ help="Enable 8-bit samples across USB")
+ parser.add_option("-p", "--pal", action="store_true", default=False,
+ help="PAL video format (this is the default)")
+ parser.add_option("-n", "--ntsc", action="store_true", default=False,
+ help="NTSC video format")
+ parser.add_option("-F", "--nframes", type="int", default=2,
+ help="number of frames to use for sync detection and
framerate")
+ parser.add_option("-S", "--search_window", type="int", default=4,
+ help="number of samples to use for correlation
search window")
+ parser.add_option("-V", "--video_alpha", type="eng_float", default=1.0,
+ help="set video alpha (multi-frame averaging)")
+ parser.add_option("-C", "--corr_alpha", type="eng_float", default=1.0,
+ help="set correlation alpha")
+ parser.add_option("-o", "--out-filename", type="string", default="sdl",
+ help="For example out_raw_uchar.gray. If you don't
specify an output filename you will get a video_sink_sdl realtime output
window. You then need to have gr-video-sdl installed)")
+ parser.add_option("-P", "--pga", type="eng_float", default=None,
+ help="set pga")
+ parser.add_option("-t", "--type", type="string",
default="do_fm_no_sync",
+ help="do_no_sync, do_fm_no_sync, do_fm_sync_corr ")
+ parser.add_option("-r", "--repeat", action="store_false", default=True,
+ help="do not repeat file in a loop")
+ parser.add_option( "--debug", action="store_true", default=False,
+ help="output debugging info")
+ parser.add_option("-N", "--no-hb", action="store_true", default=False,
+ help="don't use halfband filter in usrp")
+
+ (options, args) = parser.parse_args()
+ if not ((len(args) == 1) or (len(args) == 0)):
+ parser.print_help()
+ sys.exit(1)
+
+ if len(args) == 1:
+ filename = args[0]
+ else:
+ filename = None
+
+ self.frame = frame
+ self.panel = panel
+
+ self.contrast = options.contrast
+ self.brightness = options.brightness
+ self.state = "FREQ"
+ self.freq = 0
+ self.pga=options.pga
+
+ # build graph
+
+ self.u=None
+
+ usrp_decim = options.decim # 32
+
+ if not (options.out_filename=="sdl"):
+ options.repeat=False
+
+ if not ((filename is None) or (filename=="usrp")):
+ self.filesource =
gr.file_source(gr.sizeof_short,filename,options.repeat) # file is data source
+ self.istoc = gr.interleaved_short_to_complex()
+ self.connect(self.filesource,self.istoc)
+ self.src=self.istoc
+ options.gain=0.0
+ self.gain=0.0
+ else:
+ if options.no_hb or (options.decim<8):
+ self.fpga_filename="std_4rx_0tx.rbf" #contains 4 Rx paths without
halfbands and 0 tx paths
+ else:
+ self.fpga_filename="std_2rxhb_2tx.rbf" # contains 2 Rx paths with
halfband filters and 2 tx paths (the default)
+ self.u = usrp.source_c(0,fpga_filename=self.fpga_filename)
# usrp is data source
+ if options.width_8:
+ sample_width = 8
+ sample_shift = 8
+ format = self.u.make_format(sample_width, sample_shift)
+ r = self.u.set_format(format)
+ self.u.set_decim_rate(usrp_decim)
+ if options.rx_subdev_spec is None:
+ options.rx_subdev_spec = pick_subdevice(self.u)
+ self.u.set_mux(usrp.determine_rx_mux_value(self.u,
options.rx_subdev_spec))
+ self.subdev = usrp.selected_subdev(self.u, options.rx_subdev_spec)
+ try:
+ adc_rate = self.u.adc_rate() # 64 MS/s
+ except:
+ adc_rate=64.0e6
+ pass
+ try:
+ self.subdev.set_bw(0.8*adc_rate/float(options.decim)) #30.0e6)
#17.0e6)
+ except:
+ pass
+ print "Using RX d'board %s" % (self.subdev.side_and_name(),)
+ if options.gain is None:
+ # if no gain was specified, use the mid-point in dB
+ g = self.subdev.gain_range()
+ options.gain = float(g[0]+g[1])/2
+ self.src=self.u
+
+ usrp_rate = adc_rate / usrp_decim # 320 kS/s
+
+ f2uc=gr.float_to_uchar()
+ # sdl window as final sink
+ if not (options.pal or options.ntsc):
+ options.pal=True #set default to PAL
+ if options.pal:
+ lines_per_frame=625.0
+ frames_per_sec=25.0
+ show_width=768
+ elif options.ntsc:
+ lines_per_frame=525.0
+ frames_per_sec=29.97002997
+ show_width=640
+ width=int(usrp_rate/(lines_per_frame*frames_per_sec))
+ height=int(lines_per_frame)
+
+ if (options.out_filename=="sdl"):
+ #Here comes the tv screen, you have to build and install
gr-video-sdl for this (subproject of gnuradio, only in cvs for now)
+ try:
+ video_sink = video_sdl.sink_uc ( frames_per_sec, width,
height,0,show_width,height)
+ except:
+ print "gr-video-sdl is not installed"
+ print "realtime \"sdl\" video output window is not available"
+ raise SystemExit, 1
+ self.dst=video_sink
+ else:
+ print "You can use the imagemagick display tool to show the
resulting imagesequence"
+ print "use the following line to show the demodulated TV-signal:"
+ print "display -depth 8 -size " +str(width)+ "x" + str(height) + "
gray:" + options.out_filename
+ print "(Use the spacebar to advance to next frames)"
+ options.repeat=False
+ file_sink=gr.file_sink(gr.sizeof_char, options.out_filename)
+ self.dst =file_sink
+
+ self.agc=gr.agc_cc(1e-7,1.0,1.0) #1e-7
+ self.am_demod = gr.complex_to_mag ()
+ self.set_blacklevel=gr.add_const_ff(0.0)
+ self.invert_and_scale = gr.multiply_const_ff (0.0) #-self.contrast
*128.0*255.0/(200.0)
+
+ # now wire it all together
+ #sample_rate=options.width*options.height*options.framerate
+
+ nframes=options.nframes
+ search_window=options.search_window*nframes
+ debug=options.debug
+ video_alpha=options.video_alpha
+ corr_alpha=options.corr_alpha
+ #process_type='do_no_sync'
+ #process_type='do_fm_no_sync'
+ #process_type='do_fm_sync_corr'
+ process_type=options.type
+ if process_type=='do_no_sync':
+ self.connect (self.src,
self.agc,self.am_demod,self.invert_and_scale, self.set_blacklevel,f2uc,self.dst)
+ elif process_type=='do_fm_no_sync':
+ max_dev = 4.0e6
+ fm_demod_gain = 1.0 * 255.0 * usrp_rate/(2*math.pi*max_dev) #have to
invert because is inverted again later on in the contrat block
+ self.fm_demod = gr.quadrature_demod_cf (fm_demod_gain)
+ self.connect (self.src, self.fm_demod, self.invert_and_scale,
self.set_blacklevel,f2uc,self.dst) #self.invert_and_scale
+ elif process_type=='do_fm_sync_corr':
+ frame_size=width*height #int(usrp_rate/25.0)
+ tv_corr=gr.tv_correlator_ff(frame_size,nframes, search_window,
video_alpha, corr_alpha,debug) #Note: this block is not yet in cvs
+ max_dev = 4.0e6
+ fm_demod_gain = 1.0*usrp_rate/(2*math.pi*max_dev) #have to invert
because is inverted again later on in the contrat block
+ self.fm_demod = gr.quadrature_demod_cf (fm_demod_gain)
+ self.connect (self.src, self.fm_demod,tv_corr,self.invert_and_scale,
self.set_blacklevel,f2uc,self.dst)
+ elif process_type=='do_tv_sync_adv':
+ #defaults: gr.tv_sync_adv (double sampling_freq, unsigned int
tv_format,bool output_active_video_only=false, bool do_invert=false, double
wanted_black_level=0.0, double wanted_white_level=255.0, double avg_alpha=0.1,
double initial_gain=1.0, double initial_offset=0.0,bool debug=false)
+
self.tv_sync_adv=gr.tv_sync_adv(usrp_rate,0,False,False,0.0,255.0,0.01,1.0,0.0,False)
#note, this block is not yet in cvs
+ self.connect (self.src,
self.am_demod,self.invert_and_scale,self.tv_sync_adv,s2f,f2uc,self.dst)
+ elif process_type=='do_nullsink': #test usb bandwidth
+ c2r=gr.complex_to_real()
+ nullsink=gr.null_sink(gr.sizeof_float)
+ self.connect (self.src, c2r,nullsink)
+ elif process_type=='do_tv_sync_corr':
+ frame_size=width*height #int(usrp_rate/25.0)
+ tv_corr=gr.tv_correlator_ff(frame_size,nframes, search_window,
video_alpha, corr_alpha,debug) #Note: this block is not yet in cvs
+ shift=gr.add_const_ff(-0.7)
+ self.connect (self.src,
self.agc,self.am_demod,tv_corr,self.invert_and_scale,
self.set_blacklevel,f2uc,self.dst) #self.agc,
+ else: # process_type=='do_test_image':
+ src_vertical_bars = gr.sig_source_f (usrp_rate, gr.GR_SIN_WAVE, 10.0
*usrp_rate/320, 255,128)
+ self.connect(src_vertical_bars,f2uc,self.dst)
+
+ self._build_gui(vbox, usrp_rate, usrp_rate, usrp_rate)
+
+
+ if abs(options.freq) < 1e6:
+ options.freq *= 1e6
+
+ # set initial values
+ self.set_gain(options.gain)
+ self.set_contrast(self.contrast)
+ self.set_brightness(options.brightness)
+ if not(self.set_freq(options.freq)):
+ self._set_status_msg("Failed to set initial frequency")
+
+
+ def _set_status_msg(self, msg, which=0):
+ self.frame.GetStatusBar().SetStatusText(msg, which)
+
+
+ def _build_gui(self, vbox, usrp_rate, demod_rate, audio_rate):
+
+ def _form_set_freq(kv):
+ return self.set_freq(kv['freq'])
+
+
+ if 0:
+ self.src_fft = fftsink.fft_sink_c (self, self.panel, title="Data
from USRP",
+ fft_size=512,
sample_rate=usrp_rate)
+ self.connect (self.src, self.src_fft)
+ vbox.Add (self.src_fft.win, 4, wx.EXPAND)
+
+ if 0:
+ post_demod_fft = fftsink.fft_sink_f (self, self.panel, title="Post
Demod",
+ fft_size=512,
sample_rate=demod_rate,
+ y_per_div=10, ref_level=-40)
+ self.connect (self.am_demod, post_demod_fft)
+ vbox.Add (post_demod_fft.win, 4, wx.EXPAND)
+
+ if 0:
+ post_filt_fft = fftsink.fft_sink_f (self, self.panel, title="Post
Filter",
+ fft_size=512,
sample_rate=audio_rate,
+ y_per_div=10, ref_level=-40)
+ self.connect (self.set_blacklevel, post_filt)
+ vbox.Add (fft_win4, 4, wx.EXPAND)
+
+
+ # control area form at bottom
+ self.myform = myform = form.form()
+
+ if not (self.u is None):
+ hbox = wx.BoxSizer(wx.HORIZONTAL)
+ hbox.Add((5,0), 0)
+ myform['freq'] = form.float_field(
+ parent=self.panel, sizer=hbox, label="Freq", weight=1,
+ callback=myform.check_input_and_call(_form_set_freq,
self._set_status_msg))
+
+ hbox.Add((5,0), 0)
+ myform['freq_slider'] = \
+ form.quantized_slider_field(parent=self.panel, sizer=hbox,
weight=3,
+ range=(800.0e6, 2600.0e6, 1.0e6),
+ callback=self.set_freq)
+ hbox.Add((5,0), 0)
+ vbox.Add(hbox, 0, wx.EXPAND)
+
+ hbox = wx.BoxSizer(wx.HORIZONTAL)
+ hbox.Add((5,0), 0)
+
+ myform['contrast'] = \
+ form.quantized_slider_field(parent=self.panel, sizer=hbox,
label="Contrast",
+ weight=3, range=(-4.0, 4.0, 0.1),
+ callback=self.set_contrast)
+ hbox.Add((5,0), 1)
+
+ myform['brightness'] = \
+ form.quantized_slider_field(parent=self.panel, sizer=hbox,
label="Brightness",
+ weight=3, range=(-255.0, 255.0, 1.0),
+ callback=self.set_brightness)
+ hbox.Add((5,0), 0)
+
+ if not (self.u is None):
+ myform['gain'] = \
+ form.quantized_slider_field(parent=self.panel, sizer=hbox,
label="Gain",
+ weight=3,
range=self.subdev.gain_range(),
+ callback=self.set_gain)
+ hbox.Add((5,0), 0)
+ vbox.Add(hbox, 0, wx.EXPAND)
+
+ try:
+ self.knob = powermate.powermate(self.frame)
+ self.rot = 0
+ powermate.EVT_POWERMATE_ROTATE (self.frame, self.on_rotate)
+ powermate.EVT_POWERMATE_BUTTON (self.frame, self.on_button)
+ except:
+ print "FYI: No Powermate or Contour Knob found"
+
+
+ def on_rotate (self, event):
+ self.rot += event.delta
+ if (self.state == "FREQ"):
+ if self.rot >= 3:
+ self.set_freq(self.freq + .1e6)
+ self.rot -= 3
+ elif self.rot <=-3:
+ self.set_freq(self.freq - .1e6)
+ self.rot += 3
+ elif (self.state == "CONTRAST"):
+ step = 0.1
+ if self.rot >= 3:
+ self.set_contrast(self.contrast + step)
+ self.rot -= 3
+ elif self.rot <=-3:
+ self.set_contrast(self.contrast - step)
+ self.rot += 3
+ else:
+ step = 1
+ if self.rot >= 3:
+ self.set_brightness(self.brightness + step)
+ self.rot -= 3
+ elif self.rot <=-3:
+ self.set_brightness(self.brightness - step)
+ self.rot += 3
+
+ def on_button (self, event):
+ if event.value == 0: # button up
+ return
+ self.rot = 0
+ if self.state == "FREQ":
+ self.state = "CONTRAST"
+ elif self.state == "CONTRAST":
+ self.state = "BRIGHTNESS"
+ else:
+ self.state = "FREQ"
+ self.update_status_bar ()
+
+
+ def set_contrast (self, contrast):
+ self.contrast = contrast
+ self.invert_and_scale.set_k(self.contrast)
+ self.myform['contrast'].set_value(self.contrast)
+ self.update_status_bar ()
+
+ def set_brightness (self, brightness):
+ self.brightness = brightness
+ self.set_blacklevel.set_k(self.brightness +127.0) #255.0
+ self.myform['brightness'].set_value(self.brightness)
+ self.update_status_bar ()
+
+ def set_freq(self, target_freq):
+ """
+ Set the center frequency we're interested in.
+
+ @param target_freq: frequency in Hz
+ @rypte: bool
+
+ Tuning is a two step process. First we ask the front-end to
+ tune as close to the desired frequency as it can. Then we use
+ the result of that operation and our target_frequency to
+ determine the value for the digital down converter.
+ """
+ if not (self.u is None):
+ r = usrp.tune(self.u, 0, self.subdev, target_freq)
+ print target_freq
+ print r
+ print self.freq
+ if r:
+ self.freq = target_freq
+ self.myform['freq'].set_value(target_freq) # update
displayed value
+ self.myform['freq_slider'].set_value(target_freq) # update
displayed value
+ self.update_status_bar()
+ self._set_status_msg("OK", 0)
+ return True
+
+ self._set_status_msg("Failed", 0)
+ return False
+
+ def set_gain(self, gain):
+ if not (self.u is None):
+ self.gain=gain
+ self.myform['gain'].set_value(gain) # update displayed value
+ self.subdev.set_gain(gain)
+ if not (self.pga is None):
+ print 'old pga', self.u.pga(0)
+ self.u.set_pga(0,self.pga) #self.u.pga_max()
+ print 'overridden pga', self.u.pga(0)
+ self.update_status_bar()
+
+ def update_status_bar (self):
+ msg = "Setting:%s Freq:%d Mhz Contrast:%r Brightness:%r Gain: %r" %
(self.state,int(self.freq/1e6), self.contrast,self.brightness,self.gain)
+ self._set_status_msg(msg, 1)
+ #self.src_fft.set_baseband_freq(self.freq)
+
+
+
+if __name__ == '__main__':
+ app = stdgui.stdapp (tv_rx_graph, "USRP TV RX black-and-white")
+ app.MainLoop ()
Property changes on:
gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp/usrp_fmtv_rcv.py
___________________________________________________________________
Name: svn:executable
+ *
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Commit-gnuradio] r4447 - gnuradio/branches/developers/nldudok1/general-wip/gnuradio-examples/python/usrp,
nldudok1 <=