aboutsummaryrefslogtreecommitdiffstats
path: root/apps
diff options
context:
space:
mode:
authorRoman Khassraf <roman@khassraf.at>2016-09-23 07:26:13 +0200
committerRoman Khassraf <roman@khassraf.at>2016-09-23 07:26:13 +0200
commit98af2c97e1a3015d83cafd4467b5756c6da08424 (patch)
treef3bd6b104597982e2212a6c48659ebe39d0dcc89 /apps
parentfb772f8228023ed47ad2870ada8405c725752961 (diff)
Reflecting changes of arfcn module, re-ordering of output code to prevent messages about reattached kernel driver
Diffstat (limited to 'apps')
-rwxr-xr-xapps/grgsm_scanner240
1 files changed, 120 insertions, 120 deletions
diff --git a/apps/grgsm_scanner b/apps/grgsm_scanner
index 169a9e5..acd48d9 100755
--- a/apps/grgsm_scanner
+++ b/apps/grgsm_scanner
@@ -38,14 +38,13 @@ import pmt
import time
-#from wideband_receiver import *
+# from wideband_receiver import *
class receiver_with_decoder(grgsm.hier_block):
-
def __init__(self, OSR=4, chan_num=0, fc=939.4e6, ppm=0, samp_rate=0.2e6):
grgsm.hier_block.__init__(
self, "Receiver With Decoder",
- gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
+ gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(0, 0, 0),
)
self.message_port_register_hier_out("bursts")
@@ -63,7 +62,7 @@ class receiver_with_decoder(grgsm.hier_block):
##################################################
# Variables
##################################################
- self.samp_rate_out = samp_rate_out = 1625000.0/6.0*OSR
+ self.samp_rate_out = samp_rate_out = 1625000.0 / 6.0 * OSR
##################################################
# Blocks
@@ -88,16 +87,15 @@ class receiver_with_decoder(grgsm.hier_block):
self.msg_connect(self.gsm_control_channels_decoder_0, 'msgs', self, 'msgs')
self.msg_connect(self.gsm_receiver_0, 'C0', self.gsm_bcch_ccch_demapper_0, 'bursts')
self.msg_connect(self.gsm_receiver_0, 'measurements', self.gsm_clock_offset_control_0, 'measurements')
- self.connect((self.gsm_input_0, 0), (self.gsm_receiver_0, 0))
- self.connect((self, 0), (self.gsm_input_0, 0))
-
+ self.connect((self.gsm_input_0, 0), (self.gsm_receiver_0, 0))
+ self.connect((self, 0), (self.gsm_input_0, 0))
def get_OSR(self):
return self.OSR
def set_OSR(self, OSR):
self.OSR = OSR
- self.set_samp_rate_out(1625000.0/6.0*self.OSR)
+ self.set_samp_rate_out(1625000.0 / 6.0 * self.OSR)
self.gsm_input_0.set_osr(self.OSR)
def get_chan_num(self):
@@ -135,17 +133,16 @@ class receiver_with_decoder(grgsm.hier_block):
class wideband_receiver(grgsm.hier_block):
-
def __init__(self, OSR=4, fc=939.4e6, samp_rate=0.4e6):
grgsm.hier_block.__init__(
self, "Wideband receiver",
- gr.io_signature(1, 1, gr.sizeof_gr_complex*1),
+ gr.io_signature(1, 1, gr.sizeof_gr_complex * 1),
gr.io_signature(0, 0, 0),
)
self.message_port_register_hier_out("bursts")
self.message_port_register_hier_out("msgs")
self.__init(OSR, fc, samp_rate)
-
+
def __init(self, OSR=4, fc=939.4e6, samp_rate=0.4e6):
##################################################
# Parameters
@@ -153,9 +150,9 @@ class wideband_receiver(grgsm.hier_block):
self.OSR = OSR
self.fc = fc
self.samp_rate = samp_rate
- self.channels_num = int(samp_rate/0.2e6)
+ self.channels_num = int(samp_rate / 0.2e6)
self.OSR_PFB = 2
-
+
##################################################
# Blocks
##################################################
@@ -171,15 +168,16 @@ class wideband_receiver(grgsm.hier_block):
# Connections
##################################################
self.connect((self, 0), (self.pfb_channelizer_ccf_0, 0))
- for chan in xrange(0,self.channels_num):
+ for chan in xrange(0, self.channels_num):
self.connect((self.pfb_channelizer_ccf_0, chan), (self.receivers_with_decoders[chan], 0))
self.msg_connect(self.receivers_with_decoders[chan], 'bursts', self, 'bursts')
self.msg_connect(self.receivers_with_decoders[chan], 'msgs', self, 'msgs')
def create_receivers(self):
self.receivers_with_decoders = {}
- for chan in xrange(0,self.channels_num):
- self.receivers_with_decoders[chan] = receiver_with_decoder(fc=self.fc, OSR=self.OSR, chan_num=chan, samp_rate=self.OSR_PFB*0.2e6)
+ for chan in xrange(0, self.channels_num):
+ self.receivers_with_decoders[chan] = receiver_with_decoder(fc=self.fc, OSR=self.OSR, chan_num=chan,
+ samp_rate=self.OSR_PFB * 0.2e6)
def get_OSR(self):
return self.OSR
@@ -200,59 +198,54 @@ class wideband_receiver(grgsm.hier_block):
class wideband_scanner(gr.top_block):
-
def __init__(self, rec_len=3, sample_rate=2e6, carrier_frequency=939e6, ppm=0, args=""):
-
gr.top_block.__init__(self, "Wideband Scanner")
-
+
self.rec_len = rec_len
self.sample_rate = sample_rate
self.carrier_frequency = carrier_frequency
self.ppm = ppm
-
+
# if no file name is given process data from rtl_sdr source
- print "Args=",args
- self.rtlsdr_source = osmosdr.source( args="numchan=" + str(1) + " " + args )
-
+ print "Args=", args
+ self.rtlsdr_source = osmosdr.source(args="numchan=" + str(1) + " " + args)
+
self.rtlsdr_source.set_sample_rate(sample_rate)
-
+
# capture half of GSM channel lower than channel center (-0.1MHz)
# this is needed when even number of channels is captured in order to process full captured bandwidth
-
+
self.rtlsdr_source.set_center_freq(carrier_frequency - 0.1e6, 0)
-
+
# correction of central frequency
- # if the receiver has large frequency offset
+ # if the receiver has large frequency offset
# the value of this variable should be set close to that offset in ppm
self.rtlsdr_source.set_freq_corr(options.ppm, 0)
-
+
self.rtlsdr_source.set_dc_offset_mode(2, 0)
self.rtlsdr_source.set_iq_balance_mode(0, 0)
self.rtlsdr_source.set_gain_mode(True, 0)
self.rtlsdr_source.set_bandwidth(sample_rate, 0)
-
+
self.head = blocks.head(gr.sizeof_gr_complex * 1, int(rec_len * sample_rate))
-
+
# shift again by -0.1MHz in order to align channel center in 0Hz
- self.blocks_rotator_cc = blocks.rotator_cc(-2 * pi * 0.1e6 / options.samp_rate)
-
+ self.blocks_rotator_cc = blocks.rotator_cc(-2 * pi * 0.1e6 / options.samp_rate)
+
self.wideband_receiver = wideband_receiver(OSR=4, fc=carrier_frequency, samp_rate=sample_rate)
self.gsm_extract_system_info = grgsm.extract_system_info()
-
-
+
self.connect((self.rtlsdr_source, 0), (self.head, 0))
self.connect((self.head, 0), (self.blocks_rotator_cc, 0))
- self.connect((self.blocks_rotator_cc, 0), (self.wideband_receiver,0))
+ self.connect((self.blocks_rotator_cc, 0), (self.wideband_receiver, 0))
self.msg_connect(self.wideband_receiver, 'msgs', self.gsm_extract_system_info, 'msgs')
-
+
def set_carrier_frequency(self, carrier_frequency):
self.carrier_frequency = carrier_frequency
self.rtlsdr_source.set_center_freq(carrier_frequency - 0.1e6, 0)
-
-
+
class channel_info(object):
-
def __init__(self, arfcn, freq, cid, lac, mcc, mnc, ccch_conf, power, neighbours, cell_arfcns):
self.arfcn = arfcn
self.freq = freq
@@ -264,13 +257,13 @@ class channel_info(object):
self.power = power
self.neighbours = neighbours
self.cell_arfcns = cell_arfcns
-
+
def get_verbose_info(self):
i = " |---- Configuration: %s\n" % self.get_ccch_conf()
i += " |---- Cell ARFCNs: " + ", ".join(map(str, self.cell_arfcns)) + "\n"
i += " |---- Neighbour Cells: " + ", ".join(map(str, self.neighbours)) + "\n"
return i
-
+
def get_ccch_conf(self):
if self.ccch_conf == 0:
return "1 CCCH, not combined"
@@ -284,17 +277,18 @@ class channel_info(object):
return "4 CCCH, not combined"
else:
return "Unknown"
-
+
def getKey(self):
return self.arfcn
-
+
def __cmp__(self, other):
if hasattr(other, 'getKey'):
return self.getKey().__cmp__(other.getKey())
-
+
def __repr__(self):
- return "ARFCN: %4u, Freq: %6.1fM, CID: %5u, LAC: %5u, MCC: %3u, MNC: %3u, Pwr: %3i" % (self.arfcn, self.freq/1e6, self.cid, self.lac, self.mcc, self.mnc, self.power)
-
+ return "ARFCN: %4u, Freq: %6.1fM, CID: %5u, LAC: %5u, MCC: %3u, MNC: %3u, Pwr: %3i" % (
+ self.arfcn, self.freq / 1e6, self.cid, self.lac, self.mcc, self.mnc, self.power)
+
if __name__ == '__main__':
parser = OptionParser(option_class=eng_option, usage="%prog: [options]")
@@ -302,98 +296,104 @@ if __name__ == '__main__':
parser.add_option("-b", "--band", dest="band", default="P-GSM",
help="Specify the GSM band for the frequency.\nAvailable bands are: " + bands_list)
parser.add_option("-s", "--samp-rate", dest="samp_rate", type="float", default=2e6,
- help="Set sample rate [default=%default] - allowed values even_number*0.2e6")
+ help="Set sample rate [default=%default] - allowed values even_number*0.2e6")
parser.add_option("-p", "--ppm", dest="ppm", type="intx", default=0,
- help="Set frequency correction in ppm [default=%default]")
+ help="Set frequency correction in ppm [default=%default]")
parser.add_option("-g", "--gain", dest="gain", type="eng_float", default=24.0,
- help="Set gain [default=%default]")
+ help="Set gain [default=%default]")
parser.add_option("", "--args", dest="args", type="string", default="",
- help="Set device arguments [default=%default]")
+ help="Set device arguments [default=%default]")
parser.add_option("--speed", dest="speed", type="intx", default=4,
- help="Scan speed [default=%default]. Value range 0-5.")
- parser.add_option("-v", "--verbose", action="store_true",
+ help="Scan speed [default=%default]. Value range 0-5.")
+ parser.add_option("-v", "--verbose", action="store_true",
help="If set, verbose information output is printed: ccch configuration, cell ARFCN's, neighbour ARFCN's")
-
+
"""
Dont forget: sudo sysctl kernel.shmmni=32000
"""
-
+
(options, args) = parser.parse_args()
-
+
if options.band not in grgsm.arfcn.get_bands():
parser.error("Invalid GSM band\n")
-
+
if options.speed < 0 or options.speed > 5:
parser.error("Invalid scan speed.\n")
-
+
if (options.samp_rate / 0.2e6) % 2 != 0:
parser.error("Invalid sample rate. Sample rate must be an even numer * 0.2e6")
-
- channels_num = int(options.samp_rate/0.2e6)
-
- first_arfcn = grgsm.arfcn.get_first_arfcn(options.band)
- last_arfcn = grgsm.arfcn.get_last_arfcn(options.band)
- last_center_arfcn = last_arfcn - int((channels_num / 2) - 1)
-
- current_freq = grgsm.arfcn.arfcn2downlink(first_arfcn + int(channels_num / 2) - 1, options.band)
- last_freq = grgsm.arfcn.arfcn2downlink(last_center_arfcn, options.band)
- stop_freq = last_freq + 0.2e6 * channels_num
-
- while current_freq < stop_freq:
-
- # silence rtl_sdr output:
- # open 2 fds
- null_fds = [os.open(os.devnull, os.O_RDWR) for x in xrange(2)]
- # save the current file descriptors to a tuple
- save = os.dup(1), os.dup(2)
- # put /dev/null fds on 1 and 2
- os.dup2(null_fds[0], 1)
- os.dup2(null_fds[1], 2)
-
- # instantiate scanner and processor
- scanner = wideband_scanner(rec_len=6-options.speed,
- sample_rate=options.samp_rate,
- carrier_frequency=current_freq,
- ppm=options.ppm, args=options.args)
-
- # start recording
- scanner.start()
- scanner.wait()
- scanner.stop()
-
- # restore file descriptors so we can print the results
- os.dup2(save[0], 1)
- os.dup2(save[1], 2)
- # close the temporary fds
- os.close(null_fds[0])
- os.close(null_fds[1])
-
- freq_offsets = numpy.fft.ifftshift(numpy.array(range(int(-numpy.floor(channels_num/2)),int(numpy.floor((channels_num+1)/2))))*2e5)
- detected_c0_channels = scanner.gsm_extract_system_info.get_chans()
-
- if detected_c0_channels:
- chans = numpy.array(scanner.gsm_extract_system_info.get_chans())
- found_freqs = current_freq + freq_offsets[(chans)]
-
- cell_ids = numpy.array(scanner.gsm_extract_system_info.get_cell_id())
- lacs = numpy.array(scanner.gsm_extract_system_info.get_lac())
- mccs = numpy.array(scanner.gsm_extract_system_info.get_mcc())
- mncs = numpy.array(scanner.gsm_extract_system_info.get_mnc())
- ccch_confs = numpy.array(scanner.gsm_extract_system_info.get_ccch_conf())
- powers = numpy.array(scanner.gsm_extract_system_info.get_pwrs())
-
+
+ channels_num = int(options.samp_rate / 0.2e6)
+
+ for arfcn_range in grgsm.arfcn.get_arfcn_ranges(options.band):
+ first_arfcn = arfcn_range[0]
+ last_arfcn = arfcn_range[1]
+ last_center_arfcn = last_arfcn - int((channels_num / 2) - 1)
+
+ current_freq = grgsm.arfcn.arfcn2downlink(first_arfcn + int(channels_num / 2) - 1, options.band)
+ last_freq = grgsm.arfcn.arfcn2downlink(last_center_arfcn, options.band)
+ stop_freq = last_freq + 0.2e6 * channels_num
+
+ while current_freq < stop_freq:
+
+ # silence rtl_sdr output:
+ # open 2 fds
+ null_fds = [os.open(os.devnull, os.O_RDWR) for x in xrange(2)]
+ # save the current file descriptors to a tuple
+ save = os.dup(1), os.dup(2)
+ # put /dev/null fds on 1 and 2
+ os.dup2(null_fds[0], 1)
+ os.dup2(null_fds[1], 2)
+
+ # instantiate scanner and processor
+ scanner = wideband_scanner(rec_len=6 - options.speed,
+ sample_rate=options.samp_rate,
+ carrier_frequency=current_freq,
+ ppm=options.ppm, args=options.args)
+
+ # start recording
+ scanner.start()
+ scanner.wait()
+ scanner.stop()
+
+ freq_offsets = numpy.fft.ifftshift(
+ numpy.array(range(int(-numpy.floor(channels_num / 2)), int(numpy.floor((channels_num + 1) / 2)))) * 2e5)
+ detected_c0_channels = scanner.gsm_extract_system_info.get_chans()
+
found_list = []
- for i in range(0, len(chans)):
- cell_arfcn_list = scanner.gsm_extract_system_info.get_cell_arfcns(chans[i])
- neighbour_list = scanner.gsm_extract_system_info.get_neighbours(chans[i])
- info = channel_info(grgsm.arfcn.downlink2arfcn(found_freqs[i], options.band), found_freqs[i], cell_ids[i], lacs[i], mccs[i], mncs[i], ccch_confs[i], powers[i], neighbour_list, cell_arfcn_list)
- found_list.append(info)
-
+ if detected_c0_channels:
+ chans = numpy.array(scanner.gsm_extract_system_info.get_chans())
+ found_freqs = current_freq + freq_offsets[(chans)]
+
+ cell_ids = numpy.array(scanner.gsm_extract_system_info.get_cell_id())
+ lacs = numpy.array(scanner.gsm_extract_system_info.get_lac())
+ mccs = numpy.array(scanner.gsm_extract_system_info.get_mcc())
+ mncs = numpy.array(scanner.gsm_extract_system_info.get_mnc())
+ ccch_confs = numpy.array(scanner.gsm_extract_system_info.get_ccch_conf())
+ powers = numpy.array(scanner.gsm_extract_system_info.get_pwrs())
+
+ for i in range(0, len(chans)):
+ cell_arfcn_list = scanner.gsm_extract_system_info.get_cell_arfcns(chans[i])
+ neighbour_list = scanner.gsm_extract_system_info.get_neighbours(chans[i])
+
+ info = channel_info(grgsm.arfcn.downlink2arfcn(found_freqs[i], options.band), found_freqs[i],
+ cell_ids[i], lacs[i], mccs[i], mncs[i], ccch_confs[i], powers[i],
+ neighbour_list, cell_arfcn_list)
+ found_list.append(info)
+
+ scanner = None
+
+ # restore file descriptors so we can print the results
+ os.dup2(save[0], 1)
+ os.dup2(save[1], 2)
+ # close the temporary fds
+ os.close(null_fds[0])
+ os.close(null_fds[1])
+
for info in sorted(found_list):
print info
if options.verbose:
print info.get_verbose_info()
-
- scanner = None
- current_freq += channels_num * 0.2e6
+
+ current_freq += channels_num * 0.2e6