2020-05-18 18:58:17 +00:00
|
|
|
# NanoVNASaver
|
2020-06-25 17:52:30 +00:00
|
|
|
#
|
2020-05-18 18:58:17 +00:00
|
|
|
# A python program to view and export Touchstone data from a NanoVNA
|
2020-06-25 17:52:30 +00:00
|
|
|
# Copyright (C) 2019, 2020 Rune B. Broberg
|
2021-06-30 05:21:14 +00:00
|
|
|
# Copyright (C) 2020,2021 NanoVNA-Saver Authors
|
2019-09-01 21:13:21 +00:00
|
|
|
#
|
|
|
|
# This program is free software: you can redistribute it and/or modify
|
|
|
|
# it under the terms of the GNU General Public License as published by
|
|
|
|
# the Free Software Foundation, either version 3 of the License, or
|
|
|
|
# (at your option) any later version.
|
|
|
|
#
|
|
|
|
# This program is distributed in the hope that it will be useful,
|
|
|
|
# but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
|
# GNU General Public License for more details.
|
|
|
|
#
|
|
|
|
# You should have received a copy of the GNU General Public License
|
|
|
|
# along with this program. If not, see <https://www.gnu.org/licenses/>.
|
2019-09-05 12:56:40 +00:00
|
|
|
import math
|
2020-05-18 18:58:17 +00:00
|
|
|
import logging
|
|
|
|
from typing import List
|
2019-11-08 12:58:14 +00:00
|
|
|
|
2020-05-18 18:58:17 +00:00
|
|
|
import numpy as np
|
2019-11-08 12:58:14 +00:00
|
|
|
|
2021-07-01 19:36:17 +00:00
|
|
|
from PyQt5 import QtGui
|
2019-11-08 12:58:14 +00:00
|
|
|
|
2021-07-06 15:01:20 +00:00
|
|
|
from NanoVNASaver.Charts.Chart import Chart
|
2020-05-18 18:58:17 +00:00
|
|
|
from NanoVNASaver.RFTools import Datapoint
|
|
|
|
from .Frequency import FrequencyChart
|
|
|
|
logger = logging.getLogger(__name__)
|
2019-11-08 12:58:14 +00:00
|
|
|
|
|
|
|
|
|
|
|
class GroupDelayChart(FrequencyChart):
|
2019-11-11 22:41:23 +00:00
|
|
|
def __init__(self, name="", reflective=True):
|
2019-11-08 12:58:14 +00:00
|
|
|
super().__init__(name)
|
2021-06-27 13:08:01 +00:00
|
|
|
|
|
|
|
self.name_unit = "ns"
|
|
|
|
|
2019-11-08 12:58:14 +00:00
|
|
|
self.leftMargin = 40
|
2021-06-26 22:55:43 +00:00
|
|
|
self.dim.width = 250
|
|
|
|
self.dim.height = 250
|
2019-11-08 12:58:14 +00:00
|
|
|
self.fstart = 0
|
|
|
|
self.fstop = 0
|
|
|
|
self.minDelay = 0
|
|
|
|
self.maxDelay = 0
|
|
|
|
self.span = 0
|
|
|
|
|
2019-11-11 22:41:23 +00:00
|
|
|
self.reflective = reflective
|
|
|
|
|
2019-11-08 12:58:14 +00:00
|
|
|
self.groupDelay = []
|
|
|
|
self.groupDelayReference = []
|
|
|
|
|
|
|
|
self.minDisplayValue = -180
|
|
|
|
self.maxDisplayValue = 180
|
|
|
|
|
|
|
|
def copy(self):
|
|
|
|
new_chart: GroupDelayChart = super().copy()
|
2019-11-11 22:41:23 +00:00
|
|
|
new_chart.reflective = self.reflective
|
2019-12-04 10:08:07 +00:00
|
|
|
new_chart.groupDelay = self.groupDelay.copy()
|
|
|
|
new_chart.groupDelayReference = self.groupDelay.copy()
|
2019-11-08 12:58:14 +00:00
|
|
|
return new_chart
|
|
|
|
|
2019-11-10 15:02:42 +00:00
|
|
|
def setReference(self, data):
|
|
|
|
self.reference = data
|
|
|
|
self.calculateGroupDelay()
|
|
|
|
|
2019-11-08 12:58:14 +00:00
|
|
|
def setData(self, data):
|
|
|
|
self.data = data
|
2019-11-10 15:02:42 +00:00
|
|
|
self.calculateGroupDelay()
|
|
|
|
|
|
|
|
def calculateGroupDelay(self):
|
2022-05-14 09:00:34 +00:00
|
|
|
self.groupDelay = self.calc_data(self.data)
|
|
|
|
self.groupDelayReference = self.calc_data(self.reference)
|
2019-11-08 12:58:14 +00:00
|
|
|
self.update()
|
|
|
|
|
2022-05-14 09:00:34 +00:00
|
|
|
def calc_data(self, data: List[Datapoint]):
|
|
|
|
data_len = len(data)
|
|
|
|
if data_len <= 1:
|
|
|
|
return []
|
|
|
|
unwrapped = np.degrees(np.unwrap([d.phase for d in data]))
|
|
|
|
delay_data = []
|
|
|
|
for i, d in enumerate(data):
|
|
|
|
# TODO: Replace with call to RFTools.groupDelay
|
|
|
|
if i == 0:
|
|
|
|
phase_change = unwrapped[1] - unwrapped[i]
|
|
|
|
freq_change = data[1].freq - d.freq
|
|
|
|
elif i == data_len - 1:
|
|
|
|
phase_change = unwrapped[-1] - unwrapped[-2]
|
|
|
|
freq_change = d.freq - data[-2].freq
|
|
|
|
else:
|
|
|
|
phase_change = unwrapped[i+1] - unwrapped[i-1]
|
|
|
|
freq_change = data[i+1].freq - data[i-1].freq
|
|
|
|
delay = (-phase_change / (freq_change * 360)) * 10e8
|
|
|
|
if not self.reflective:
|
|
|
|
delay /= 2
|
|
|
|
delay_data.append(delay)
|
|
|
|
return delay_data
|
|
|
|
|
2019-11-08 12:58:14 +00:00
|
|
|
def drawValues(self, qp: QtGui.QPainter):
|
|
|
|
if len(self.data) == 0 and len(self.reference) == 0:
|
|
|
|
return
|
2021-07-06 15:01:20 +00:00
|
|
|
pen = QtGui.QPen(Chart.color.sweep)
|
2021-06-26 22:55:43 +00:00
|
|
|
pen.setWidth(self.dim.point)
|
2021-07-06 15:01:20 +00:00
|
|
|
line_pen = QtGui.QPen(Chart.color.sweep)
|
2021-06-26 22:55:43 +00:00
|
|
|
line_pen.setWidth(self.dim.line)
|
2019-11-08 12:58:14 +00:00
|
|
|
|
|
|
|
if self.fixedValues:
|
|
|
|
min_delay = self.minDisplayValue
|
|
|
|
max_delay = self.maxDisplayValue
|
|
|
|
elif self.data:
|
|
|
|
min_delay = math.floor(np.min(self.groupDelay))
|
|
|
|
max_delay = math.ceil(np.max(self.groupDelay))
|
|
|
|
elif self.reference:
|
|
|
|
min_delay = math.floor(np.min(self.groupDelayReference))
|
|
|
|
max_delay = math.ceil(np.max(self.groupDelayReference))
|
|
|
|
|
|
|
|
span = max_delay - min_delay
|
|
|
|
if span == 0:
|
|
|
|
span = 0.01
|
|
|
|
self.minDelay = min_delay
|
|
|
|
self.maxDelay = max_delay
|
|
|
|
self.span = span
|
|
|
|
|
2021-06-26 22:55:43 +00:00
|
|
|
tickcount = math.floor(self.dim.height / 60)
|
2019-11-08 12:58:14 +00:00
|
|
|
for i in range(tickcount):
|
|
|
|
delay = min_delay + span * i / tickcount
|
2021-06-26 22:55:43 +00:00
|
|
|
y = self.topMargin + round((self.maxDelay - delay) / self.span * self.dim.height)
|
2022-05-14 09:00:34 +00:00
|
|
|
if delay not in {min_delay, max_delay}:
|
2021-07-06 15:01:20 +00:00
|
|
|
qp.setPen(QtGui.QPen(Chart.color.text))
|
2022-05-14 09:00:34 +00:00
|
|
|
# TODO use format class
|
|
|
|
digits = 0 if delay == 0 else max(
|
|
|
|
0, min(2, math.floor(3 - math.log10(abs(delay)))))
|
|
|
|
delaystr = str(round(delay, digits if digits != 0 else None))
|
2019-11-08 12:58:14 +00:00
|
|
|
qp.drawText(3, y + 3, delaystr)
|
2021-07-06 15:01:20 +00:00
|
|
|
qp.setPen(QtGui.QPen(Chart.color.foreground))
|
2021-06-26 22:55:43 +00:00
|
|
|
qp.drawLine(self.leftMargin - 5, y, self.leftMargin + self.dim.width, y)
|
2022-05-14 09:00:34 +00:00
|
|
|
|
2020-05-18 18:58:17 +00:00
|
|
|
qp.drawLine(self.leftMargin - 5,
|
|
|
|
self.topMargin,
|
2021-06-26 22:55:43 +00:00
|
|
|
self.leftMargin + self.dim.width,
|
2020-05-18 18:58:17 +00:00
|
|
|
self.topMargin)
|
2021-07-06 15:01:20 +00:00
|
|
|
qp.setPen(Chart.color.text)
|
2019-11-08 12:58:14 +00:00
|
|
|
qp.drawText(3, self.topMargin + 5, str(max_delay))
|
2021-06-26 22:55:43 +00:00
|
|
|
qp.drawText(3, self.dim.height + self.topMargin, str(min_delay))
|
2019-11-08 12:58:14 +00:00
|
|
|
|
2021-06-27 09:45:22 +00:00
|
|
|
self._set_start_stop()
|
2019-11-08 12:58:14 +00:00
|
|
|
|
|
|
|
# Draw bands if required
|
|
|
|
if self.bands.enabled:
|
2021-06-27 09:45:22 +00:00
|
|
|
self.drawBands(qp, self.fstart, self.fstop)
|
2019-11-08 12:58:14 +00:00
|
|
|
|
|
|
|
self.drawFrequencyTicks(qp)
|
|
|
|
|
2022-05-14 09:00:34 +00:00
|
|
|
self.draw_data(qp, Chart.color.sweep,
|
|
|
|
self.data, self.groupDelay)
|
|
|
|
self.draw_data(qp, Chart.color.reference,
|
|
|
|
self.reference, self.groupDelayReference)
|
2019-12-09 15:22:25 +00:00
|
|
|
|
2022-05-14 09:00:34 +00:00
|
|
|
self.drawMarkers(qp)
|
|
|
|
|
|
|
|
def draw_data(self, qp: QtGui.QPainter, color: QtGui.QColor,
|
|
|
|
data: List[Datapoint], delay: List[Datapoint]):
|
2019-12-09 15:22:25 +00:00
|
|
|
pen = QtGui.QPen(color)
|
2021-06-26 22:55:43 +00:00
|
|
|
pen.setWidth(self.dim.point)
|
2019-12-09 15:22:25 +00:00
|
|
|
line_pen = QtGui.QPen(color)
|
2021-06-26 22:55:43 +00:00
|
|
|
line_pen.setWidth(self.dim.line)
|
2019-12-09 15:22:25 +00:00
|
|
|
qp.setPen(pen)
|
2022-05-14 09:00:34 +00:00
|
|
|
for i, d in enumerate(data):
|
|
|
|
x = self.getXPosition(d)
|
|
|
|
y = self.getYPositionFromDelay(delay[i])
|
2019-12-09 15:22:25 +00:00
|
|
|
if self.isPlotable(x, y):
|
|
|
|
qp.drawPoint(int(x), int(y))
|
2021-06-27 08:59:07 +00:00
|
|
|
if self.flag.draw_lines and i > 0:
|
2022-05-14 09:00:34 +00:00
|
|
|
prevx = self.getXPosition(data[i - 1])
|
|
|
|
prevy = self.getYPositionFromDelay(delay[i - 1])
|
2019-12-09 15:22:25 +00:00
|
|
|
qp.setPen(line_pen)
|
|
|
|
if self.isPlotable(x, y) and self.isPlotable(prevx, prevy):
|
|
|
|
qp.drawLine(x, y, prevx, prevy)
|
|
|
|
elif self.isPlotable(x, y) and not self.isPlotable(prevx, prevy):
|
|
|
|
new_x, new_y = self.getPlotable(x, y, prevx, prevy)
|
|
|
|
qp.drawLine(x, y, new_x, new_y)
|
|
|
|
elif not self.isPlotable(x, y) and self.isPlotable(prevx, prevy):
|
|
|
|
new_x, new_y = self.getPlotable(prevx, prevy, x, y)
|
|
|
|
qp.drawLine(prevx, prevy, new_x, new_y)
|
|
|
|
qp.setPen(pen)
|
|
|
|
|
2019-11-08 12:58:14 +00:00
|
|
|
def getYPosition(self, d: Datapoint) -> int:
|
2022-05-14 09:00:34 +00:00
|
|
|
# TODO: Find a faster way than these expensive "d in data" lookups
|
|
|
|
try:
|
2019-11-08 12:58:14 +00:00
|
|
|
delay = self.groupDelay[self.data.index(d)]
|
2022-05-14 09:00:34 +00:00
|
|
|
except ValueError:
|
|
|
|
try:
|
|
|
|
delay = self.groupDelayReference[self.reference.index(d)]
|
|
|
|
except ValueError:
|
|
|
|
delay = 0
|
2019-12-09 15:22:25 +00:00
|
|
|
return self.getYPositionFromDelay(delay)
|
|
|
|
|
|
|
|
def getYPositionFromDelay(self, delay: float):
|
2021-06-26 22:55:43 +00:00
|
|
|
return self.topMargin + round((self.maxDelay - delay) / self.span * self.dim.height)
|
2019-11-08 12:58:14 +00:00
|
|
|
|
|
|
|
def valueAtPosition(self, y) -> List[float]:
|
|
|
|
absy = y - self.topMargin
|
2021-06-26 22:55:43 +00:00
|
|
|
val = -1 * ((absy / self.dim.height * self.span) - self.maxDelay)
|
2019-11-08 12:58:14 +00:00
|
|
|
return [val]
|