Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 1 | #!/usr/bin/env python |
| 2 | # -*- coding: utf-8 -*- |
| 3 | # @file |
Piotr Krysik | a6268a5 | 2017-08-23 16:02:19 +0200 | [diff] [blame] | 4 | # @author (C) 2015 by Roman Khassraf <rkhassraf@gmail.com> |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 5 | # @section LICENSE |
| 6 | # |
| 7 | # Gr-gsm 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 3, or (at your option) |
| 10 | # any later version. |
| 11 | # |
| 12 | # Gr-gsm 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. |
| 16 | # |
| 17 | # You should have received a copy of the GNU General Public License |
| 18 | # along with gr-gsm; see the file COPYING. If not, write to |
| 19 | # the Free Software Foundation, Inc., 51 Franklin Street, |
| 20 | # Boston, MA 02110-1301, USA. |
| 21 | # |
| 22 | # |
| 23 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 24 | from gnuradio import blocks |
| 25 | from gnuradio import eng_notation |
| 26 | from gnuradio import gr |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 27 | from gnuradio.eng_option import eng_option |
| 28 | from gnuradio.filter import firdes |
| 29 | from math import pi |
| 30 | from optparse import OptionParser |
| 31 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 32 | import grgsm |
| 33 | import osmosdr |
| 34 | import pmt |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 35 | import signal |
| 36 | import sys |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 37 | |
| 38 | |
Piotr Krysik | 3e192c4 | 2016-02-13 07:31:51 +0100 | [diff] [blame] | 39 | class grgsm_capture(gr.top_block): |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 40 | |
Piotr Krysik | df0adf3 | 2017-09-11 09:35:56 +0200 | [diff] [blame] | 41 | def __init__(self, fc, gain, samp_rate, ppm, arfcn, cfile=None, burst_file=None, verbose=False, rec_length=None, args=""): |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 42 | |
Piotr Krysik | 3e192c4 | 2016-02-13 07:31:51 +0100 | [diff] [blame] | 43 | gr.top_block.__init__(self, "Gr-gsm Capture") |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 44 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 45 | ################################################## |
| 46 | # Parameters |
| 47 | ################################################## |
| 48 | self.fc = fc |
| 49 | self.gain = gain |
| 50 | self.samp_rate = samp_rate |
| 51 | self.ppm = ppm |
| 52 | self.arfcn = arfcn |
| 53 | self.cfile = cfile |
| 54 | self.burst_file = burst_file |
| 55 | self.verbose = verbose |
| 56 | self.shiftoff = shiftoff = 400e3 |
Piotr Krysik | 351d217 | 2015-09-20 13:53:06 +0200 | [diff] [blame] | 57 | self.rec_length = rec_length |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 58 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 59 | ################################################## |
| 60 | # Processing Blocks |
| 61 | ################################################## |
| 62 | |
Piotr Krysik | 0ff0d09 | 2016-03-18 11:39:09 +0100 | [diff] [blame] | 63 | self.rtlsdr_source = osmosdr.source( args="numchan=" + str(1) + " " + args ) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 64 | self.rtlsdr_source.set_sample_rate(samp_rate) |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 65 | self.rtlsdr_source.set_center_freq(fc - shiftoff, 0) |
| 66 | self.rtlsdr_source.set_freq_corr(ppm, 0) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 67 | self.rtlsdr_source.set_dc_offset_mode(2, 0) |
| 68 | self.rtlsdr_source.set_iq_balance_mode(2, 0) |
| 69 | self.rtlsdr_source.set_gain_mode(True, 0) |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 70 | self.rtlsdr_source.set_gain(gain, 0) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 71 | self.rtlsdr_source.set_if_gain(20, 0) |
| 72 | self.rtlsdr_source.set_bb_gain(20, 0) |
| 73 | self.rtlsdr_source.set_antenna("", 0) |
| 74 | self.rtlsdr_source.set_bandwidth(250e3+abs(shiftoff), 0) |
| 75 | self.blocks_rotator = blocks.rotator_cc(-2*pi*shiftoff/samp_rate) |
Piotr Krysik | 351d217 | 2015-09-20 13:53:06 +0200 | [diff] [blame] | 76 | |
| 77 | if self.rec_length is not None: |
| 78 | self.blocks_head_0 = blocks.head(gr.sizeof_gr_complex, int(samp_rate*rec_length)) |
| 79 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 80 | if self.verbose or self.burst_file: |
| 81 | self.gsm_receiver = grgsm.receiver(4, ([self.arfcn]), ([])) |
| 82 | self.gsm_input = grgsm.gsm_input( |
| 83 | ppm=0, |
| 84 | osr=4, |
| 85 | fc=fc, |
| 86 | samp_rate_in=samp_rate, |
| 87 | ) |
Roman Khassraf | e87e955 | 2016-10-01 13:00:43 +0200 | [diff] [blame] | 88 | self.gsm_clock_offset_control = grgsm.clock_offset_control(fc-shiftoff, samp_rate, osr=4) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 89 | |
| 90 | if self.burst_file: |
| 91 | self.gsm_burst_file_sink = grgsm.burst_file_sink(self.burst_file) |
| 92 | |
| 93 | if self.cfile: |
| 94 | self.blocks_file_sink = blocks.file_sink(gr.sizeof_gr_complex*1, self.cfile, False) |
| 95 | self.blocks_file_sink.set_unbuffered(False) |
| 96 | |
| 97 | if self.verbose: |
| 98 | self.gsm_bursts_printer_0 = grgsm.bursts_printer(pmt.intern(""), |
| 99 | False, False, False, False) |
| 100 | |
| 101 | ################################################## |
| 102 | # Connections |
| 103 | ################################################## |
| 104 | |
Piotr Krysik | 351d217 | 2015-09-20 13:53:06 +0200 | [diff] [blame] | 105 | if self.rec_length is not None: #if recording length is defined connect head block after the source |
| 106 | self.connect((self.rtlsdr_source, 0), (self.blocks_head_0, 0)) |
| 107 | self.connect((self.blocks_head_0, 0), (self.blocks_rotator, 0)) |
| 108 | else: |
| 109 | self.connect((self.rtlsdr_source, 0), (self.blocks_rotator, 0)) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 110 | |
| 111 | if self.cfile: |
| 112 | self.connect((self.blocks_rotator, 0), (self.blocks_file_sink, 0)) |
| 113 | |
| 114 | if self.verbose or self.burst_file: |
| 115 | self.connect((self.gsm_input, 0), (self.gsm_receiver, 0)) |
| 116 | self.connect((self.blocks_rotator, 0), (self.gsm_input, 0)) |
Roman Khassraf | 1245b22 | 2016-10-01 13:13:05 +0200 | [diff] [blame] | 117 | self.msg_connect(self.gsm_clock_offset_control, "ctrl", self.gsm_input, "ctrl_in") |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 118 | self.msg_connect(self.gsm_receiver, "measurements", self.gsm_clock_offset_control, "measurements") |
| 119 | |
| 120 | if self.burst_file: |
| 121 | self.msg_connect(self.gsm_receiver, "C0", self.gsm_burst_file_sink, "in") |
| 122 | if self.verbose: |
| 123 | self.msg_connect(self.gsm_receiver, "C0", self.gsm_bursts_printer_0, "bursts") |
| 124 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 125 | def get_fc(self): |
| 126 | return self.fc |
| 127 | |
| 128 | def set_fc(self, fc): |
| 129 | self.fc = fc |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 130 | if self.verbose or self.burst_file: |
| 131 | self.gsm_input.set_fc(self.fc) |
Roman Khassraf | 5a80400 | 2015-08-24 20:11:10 +0200 | [diff] [blame] | 132 | |
| 133 | def get_arfcn(self): |
| 134 | return self.arfcn |
| 135 | |
| 136 | def set_arfcn(self, arfcn): |
| 137 | self.arfcn = arfcn |
| 138 | if self.verbose or self.burst_file: |
| 139 | self.gsm_receiver.set_cell_allocation([self.arfcn]) |
leommxj | caa7741 | 2017-09-10 14:49:03 +0800 | [diff] [blame] | 140 | new_freq = grgsm.arfcn.arfcn2downlink(self.arfcn) |
Roman Khassraf | 5a80400 | 2015-08-24 20:11:10 +0200 | [diff] [blame] | 141 | self.set_fc(new_freq) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 142 | |
| 143 | def get_gain(self): |
| 144 | return self.gain |
| 145 | |
| 146 | def set_gain(self, gain): |
| 147 | self.gain = gain |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 148 | |
| 149 | def get_samp_rate(self): |
| 150 | return self.samp_rate |
| 151 | |
| 152 | def set_samp_rate(self, samp_rate): |
| 153 | self.samp_rate = samp_rate |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 154 | self.rtlsdr_source.set_sample_rate(self.samp_rate) |
| 155 | if self.verbose or self.burst_file: |
| 156 | self.gsm_input.set_samp_rate_in(self.samp_rate) |
| 157 | |
| 158 | def get_ppm(self): |
| 159 | return self.ppm |
| 160 | |
| 161 | def set_ppm(self, ppm): |
| 162 | self.ppm = ppm |
| 163 | self.set_ppm_slider(self.ppm) |
| 164 | |
Piotr Krysik | 351d217 | 2015-09-20 13:53:06 +0200 | [diff] [blame] | 165 | def get_rec_length(self): |
| 166 | return self.rec_length |
| 167 | |
| 168 | def set_rec_length(self, rec_length): |
| 169 | self.rec_length = rec_length |
| 170 | self.blocks_head_0.set_length(int(self.samp_rate*self.rec_length)) |
| 171 | |
| 172 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 173 | if __name__ == '__main__': |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 174 | |
| 175 | parser = OptionParser(option_class=eng_option, usage="%prog [options]", |
| 176 | description="RTL-SDR capturing app of gr-gsm.") |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 177 | |
| 178 | parser.add_option("-f", "--fc", dest="fc", type="eng_float", |
| 179 | help="Set frequency [default=%default]") |
| 180 | |
| 181 | parser.add_option("-a", "--arfcn", dest="arfcn", type="intx", |
Piotr Krysik | df0adf3 | 2017-09-11 09:35:56 +0200 | [diff] [blame] | 182 | help="Set ARFCN instead of frequency (for PCS1900 add 0x8000 (2**15) to the ARFCN number)") |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 183 | |
| 184 | parser.add_option("-g", "--gain", dest="gain", type="eng_float", |
| 185 | default=eng_notation.num_to_str(30), |
| 186 | help="Set gain [default=%default]") |
| 187 | |
| 188 | parser.add_option("-s", "--samp-rate", dest="samp_rate", type="eng_float", |
Piotr Krysik | 9648540 | 2016-07-21 16:06:03 +0200 | [diff] [blame] | 189 | default=eng_notation.num_to_str(1000000), |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 190 | help="Set samp_rate [default=%default]") |
| 191 | |
| 192 | parser.add_option("-p", "--ppm", dest="ppm", type="intx", default=0, |
| 193 | help="Set ppm [default=%default]") |
| 194 | |
| 195 | parser.add_option("-b", "--burst-file", dest="burst_file", |
| 196 | help="File where the captured bursts are saved") |
| 197 | |
| 198 | parser.add_option("-c", "--cfile", dest="cfile", |
| 199 | help="File where the captured data are saved") |
| 200 | |
Piotr Krysik | 7185b66 | 2016-02-14 20:24:54 +0100 | [diff] [blame] | 201 | parser.add_option("", "--args", dest="args", type="string", default="", |
| 202 | help="Set device arguments [default=%default]") |
| 203 | |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 204 | parser.add_option("-v", "--verbose", action="store_true", |
| 205 | help="If set, the captured bursts are printed to stdout") |
Piotr Krysik | 351d217 | 2015-09-20 13:53:06 +0200 | [diff] [blame] | 206 | |
| 207 | parser.add_option("-T", "--rec-length", dest="rec_length", type="eng_float", |
| 208 | help="Set length of recording in seconds [default=%default]") |
| 209 | |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 210 | (options, args) = parser.parse_args() |
| 211 | |
| 212 | if options.cfile is None and options.burst_file is None: |
| 213 | parser.error("Please provide a cfile or a burst file (or both) to save the captured data\n") |
| 214 | |
| 215 | if (options.fc is None and options.arfcn is None) or (options.fc is not None and options.arfcn is not None): |
| 216 | parser.error("You have to provide either a frequency or an ARFCN (but not both).\n") |
| 217 | |
| 218 | arfcn = 0 |
| 219 | fc = 939.4e6 |
| 220 | if options.arfcn: |
Piotr Krysik | a80567f | 2017-09-11 10:01:12 +0200 | [diff] [blame] | 221 | if not grgsm.arfcn.is_valid_arfcn(options.arfcn): |
| 222 | parser.error("ARFCN is not valid\n") |
| 223 | else: |
| 224 | arfcn = options.arfcn |
| 225 | fc = grgsm.arfcn.arfcn2downlink(arfcn) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 226 | elif options.fc: |
| 227 | fc = options.fc |
Piotr Krysik | df0adf3 | 2017-09-11 09:35:56 +0200 | [diff] [blame] | 228 | arfcn = grgsm.arfcn.downlink2arfcn(options.fc) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 229 | |
Piotr Krysik | 3e192c4 | 2016-02-13 07:31:51 +0100 | [diff] [blame] | 230 | tb = grgsm_capture(fc=fc, gain=options.gain, samp_rate=options.samp_rate, |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 231 | ppm=options.ppm, arfcn=arfcn, cfile=options.cfile, |
Piotr Krysik | a80567f | 2017-09-11 10:01:12 +0200 | [diff] [blame] | 232 | burst_file=options.burst_file, verbose=options.verbose, |
Piotr Krysik | 0ff0d09 | 2016-03-18 11:39:09 +0100 | [diff] [blame] | 233 | rec_length=options.rec_length, args=options.args) |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 234 | |
| 235 | def signal_handler(signal, frame): |
| 236 | tb.stop() |
| 237 | tb.wait() |
| 238 | sys.exit(0) |
| 239 | |
| 240 | signal.signal(signal.SIGINT, signal_handler) |
Roman Khassraf | 529e215 | 2015-08-22 19:52:41 +0200 | [diff] [blame] | 241 | |
| 242 | tb.start() |
Roman Khassraf | 7eda480 | 2015-08-24 20:00:48 +0200 | [diff] [blame] | 243 | tb.wait() |