FreeDATA/test/test_tx.py

176 lines
6.2 KiB
Python
Raw Normal View History

2021-12-15 21:20:24 +00:00
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import ctypes
from ctypes import *
import pathlib
import sounddevice as sd
import time
import argparse
import sys
sys.path.insert(0,'..')
from tnc import codec2
import numpy as np
2021-12-12 11:09:13 +00:00
# GET PARAMETER INPUTS
parser = argparse.ArgumentParser(description='Simons TEST TNC')
parser.add_argument('--bursts', dest="N_BURSTS", default=1, type=int)
parser.add_argument('--framesperburst', dest="N_FRAMES_PER_BURST", default=1, type=int)
parser.add_argument('--delay', dest="DELAY_BETWEEN_BURSTS", default=500, type=int,
help="delay between bursts in ms")
2022-04-15 09:29:40 +00:00
parser.add_argument('--mode', dest="FREEDV_MODE", type=str, choices=['datac0', 'datac1', 'datac3'])
2021-12-14 22:56:55 +00:00
parser.add_argument('--audiodev', dest="AUDIO_OUTPUT_DEVICE", default=-1, type=int,
help="audio output device number to use, use -2 to automatically select a loopback device")
parser.add_argument('--list', dest="LIST", action="store_true", help="list audio devices by number and exit")
parser.add_argument('--testframes', dest="TESTFRAMES", action="store_true", default=False, help="list audio devices by number and exit")
args = parser.parse_args()
if args.LIST:
devices = sd.query_devices(device=None, kind=None)
index = 0
for device in devices:
print(f"{index} {device['name']}")
index += 1
sd._terminate()
quit()
N_BURSTS = args.N_BURSTS
N_FRAMES_PER_BURST = args.N_FRAMES_PER_BURST
DELAY_BETWEEN_BURSTS = args.DELAY_BETWEEN_BURSTS/1000
AUDIO_OUTPUT_DEVICE = args.AUDIO_OUTPUT_DEVICE
MODE = codec2.FREEDV_MODE[args.FREEDV_MODE].value
2021-12-12 11:09:13 +00:00
# AUDIO PARAMETERS
AUDIO_FRAMES_PER_BUFFER = 2400
MODEM_SAMPLE_RATE = codec2.api.FREEDV_FS_8000
AUDIO_SAMPLE_RATE_TX = 48000
assert (AUDIO_SAMPLE_RATE_TX % MODEM_SAMPLE_RATE) == 0
2021-12-12 11:09:13 +00:00
# check if we want to use an audio device then do an pyaudio init
if AUDIO_OUTPUT_DEVICE != -1:
2021-12-14 22:07:01 +00:00
# auto search for loopback devices
if AUDIO_OUTPUT_DEVICE == -2:
loopback_list = []
devices = sd.query_devices(device=None, kind=None)
index = 0
for device in devices:
if 'Loopback: PCM' in device['name']:
print(index)
loopback_list.append(index)
index += 1
if len(loopback_list) >= 1:
AUDIO_OUTPUT_DEVICE = loopback_list[len(loopback_list)-1] #0 = RX 1 = TX
2021-12-14 22:07:01 +00:00
print(f"loopback_list tx: {loopback_list}", file=sys.stderr)
else:
print("not enough audio loopback devices ready...")
print("you should wait about 30 seconds...")
sd._terminate()
2021-12-14 22:07:01 +00:00
quit()
print(f"AUDIO OUTPUT DEVICE: {AUDIO_OUTPUT_DEVICE}", file=sys.stderr)
# audio stream init
stream_tx = sd.RawStream(channels=1, dtype='int16', device=(0, AUDIO_OUTPUT_DEVICE), samplerate = AUDIO_SAMPLE_RATE_TX, blocksize=4800)
resampler = codec2.resampler()
2021-12-12 11:09:13 +00:00
# data binary string
if args.TESTFRAMES:
data_out = bytearray(14)
data_out[:1] = bytes([255])
data_out[1:2] = bytes([1])
data_out[2:] = b'HELLO WORLD'
else:
data_out = b'HELLO WORLD!'
2022-04-15 09:29:40 +00:00
2021-12-12 11:09:13 +00:00
# ----------------------------------------------------------------
2022-04-15 09:29:40 +00:00
# open codec2 instance
freedv = cast(codec2.api.freedv_open(MODE), c_void_p)
2021-12-12 11:09:13 +00:00
# get number of bytes per frame for mode
2021-12-12 20:52:03 +00:00
bytes_per_frame = int(codec2.api.freedv_get_bits_per_modem_frame(freedv)/8)
payload_bytes_per_frame = bytes_per_frame -2
2022-04-15 09:29:40 +00:00
2021-12-12 11:09:13 +00:00
# init buffer for data
2021-12-12 20:52:03 +00:00
n_tx_modem_samples = codec2.api.freedv_get_n_tx_modem_samples(freedv)
2021-12-12 11:09:13 +00:00
mod_out = create_string_buffer(n_tx_modem_samples * 2)
# init buffer for preample
2021-12-12 20:52:03 +00:00
n_tx_preamble_modem_samples = codec2.api.freedv_get_n_tx_preamble_modem_samples(freedv)
2021-12-12 11:09:13 +00:00
mod_out_preamble = create_string_buffer(n_tx_preamble_modem_samples * 2)
2021-12-12 11:09:13 +00:00
# init buffer for postamble
2021-12-12 20:52:03 +00:00
n_tx_postamble_modem_samples = codec2.api.freedv_get_n_tx_postamble_modem_samples(freedv)
2021-12-12 11:09:13 +00:00
mod_out_postamble = create_string_buffer(n_tx_postamble_modem_samples * 2)
2022-04-15 09:29:40 +00:00
2021-12-12 11:09:13 +00:00
# create buffer for data
buffer = bytearray(payload_bytes_per_frame) # use this if CRC16 checksum is required ( DATA1-3)
buffer[:len(data_out)] = data_out # set buffersize to length of data which will be send
2021-12-12 11:09:13 +00:00
# create crc for data frame - we are using the crc function shipped with codec2 to avoid
# crc algorithm incompatibilities
2021-12-12 20:52:03 +00:00
crc = ctypes.c_ushort(codec2.api.freedv_gen_crc16(bytes(buffer), payload_bytes_per_frame)) # generate CRC16
crc = crc.value.to_bytes(2, byteorder='big') # convert crc to 2 byte hex string
buffer += crc # append crc16 to buffer
2021-12-12 11:09:13 +00:00
print(f"TOTAL BURSTS: {N_BURSTS} TOTAL FRAMES_PER_BURST: {N_FRAMES_PER_BURST}", file=sys.stderr)
2021-12-12 11:09:13 +00:00
for i in range(1,N_BURSTS+1):
2021-12-12 11:09:13 +00:00
# write preamble to txbuffer
2021-12-12 20:52:03 +00:00
codec2.api.freedv_rawdatapreambletx(freedv, mod_out_preamble)
2022-04-15 09:29:40 +00:00
txbuffer = bytes(mod_out_preamble)
2021-12-12 11:09:13 +00:00
# create modulaton for N = FRAMESPERBURST and append it to txbuffer
for n in range(1,N_FRAMES_PER_BURST+1):
data = (ctypes.c_ubyte * bytes_per_frame).from_buffer_copy(buffer)
2021-12-12 20:52:03 +00:00
codec2.api.freedv_rawdatatx(freedv,mod_out,data) # modulate DATA and save it into mod_out pointer
txbuffer += bytes(mod_out)
2022-04-15 09:29:40 +00:00
2021-12-14 22:07:01 +00:00
print(f"TX BURST: {i}/{N_BURSTS} FRAME: {n}/{N_FRAMES_PER_BURST}", file=sys.stderr)
2021-12-12 11:09:13 +00:00
# append postamble to txbuffer
2021-12-12 20:52:03 +00:00
codec2.api.freedv_rawdatapostambletx(freedv, mod_out_postamble)
2022-04-15 09:29:40 +00:00
txbuffer += bytes(mod_out_postamble)
2021-12-12 11:09:13 +00:00
# append a delay between bursts as audio silence
samples_delay = int(MODEM_SAMPLE_RATE*DELAY_BETWEEN_BURSTS)
mod_out_silence = create_string_buffer(samples_delay*2)
txbuffer += bytes(mod_out_silence)
#print(f"samples_delay: {samples_delay} DELAY_BETWEEN_BURSTS: {DELAY_BETWEEN_BURSTS}", file=sys.stderr)
# resample up to 48k (resampler works on np.int16)
x = np.frombuffer(txbuffer, dtype=np.int16)
txbuffer_48k = resampler.resample8_to_48(x)
2021-12-12 11:09:13 +00:00
# check if we want to use an audio device or stdout
if AUDIO_OUTPUT_DEVICE != -1:
stream_tx.start()
stream_tx.write(txbuffer_48k)
else:
2021-12-12 11:09:13 +00:00
# print data to terminal for piping the output to other programs
sys.stdout.buffer.write(txbuffer_48k)
sys.stdout.flush()
2021-12-12 11:09:13 +00:00
# and at last check if we had an opened audio instance and close it
if AUDIO_OUTPUT_DEVICE != -1:
sd._terminate()
2022-04-15 09:29:40 +00:00