Compare commits

..

3 commits

Author SHA1 Message Date
e6dc6b06cb
Change report format slightly and add options 2024-04-25 21:33:38 +02:00
11dd849b5c
Add shebang 2024-04-25 21:32:57 +02:00
0baab314ed
Improve plotting script 2024-04-25 21:31:59 +02:00
4 changed files with 70 additions and 40 deletions

1
.gitignore vendored
View file

@ -3,3 +3,4 @@ legacy/
__pycache__/ __pycache__/
*.pyc *.pyc
*.pyo *.pyo
*.png

2
numberdealers/users.py Normal file → Executable file
View file

@ -1,3 +1,5 @@
#!/usr/bin/env python3
from os.path import expanduser from os.path import expanduser
import json import json

View file

@ -1,7 +1,7 @@
#!/usr/bin/env python3 #!/usr/bin/env python3
from datetime import timezone, datetime from datetime import timezone, datetime
import parse_numberdealers from numberdealers import parse_numberdealers
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import numpy as np import numpy as np
@ -10,8 +10,8 @@ def xy(messages):
x = [] x = []
y = [] y = []
for msg in messages: for msg in messages:
x.append(datetime.fromtimestamp(msg.create_at / 1000, timezone.utc)) # x.append(datetime.fromtimestamp(msg.create_at / 1000, timezone.utc))
# x.append(msg.create_at / 1000) x.append(msg.create_at / 1000)
y.append(msg.recognized_number) y.append(msg.recognized_number)
return x, y return x, y
@ -22,44 +22,71 @@ def find(xs, condition):
return i return i
def abline(slope, intercept): def abline(slope, intercept, point1_x, point2_x):
"""Plot a line from slope and intercept""" """Plot a line from slope and intercept"""
axes = plt.gca() axes = plt.gca()
x_vals = np.array(axes.get_xlim()) x_vals = np.array([point1_x, point2_x])
y_vals = slope * x_vals + intercept y_vals = slope * x_vals + intercept
plt.plot(x_vals, y_vals, '--') plt.plot(x_vals, y_vals, linewidth=1, color='grey')
def main(): def main():
import sys import sys
if len(sys.argv) <= 2:
print("Usage: plot_numberdealers.py <numberdealers.json> <numberdealers-ng.json>")
return 1
with open(sys.argv[1], encoding="utf-8") as fh: with open(sys.argv[1], encoding="utf-8") as fh:
numbers_og, _errors = parse_numberdealers.parse(fh) numbers_og, _errors = parse_numberdealers.parse(fh)
with open(sys.argv[2], encoding="utf-8") as fh: with open(sys.argv[2], encoding="utf-8") as fh:
numbers_ng, _errors = parse_numberdealers.parse(fh) numbers_ng, _errors = parse_numberdealers.parse(fh)
# start_of_current_slope = find(numbers_og, lambda msg: msg.recognized_number <= 10464) start_of_current_slope = find(numbers_og, lambda msg: msg.create_at >= numbers_ng[0].create_at)
# slope_og = ( # a = (y2-y1) / (x2-x1)
# numbers_og[-1].recognized_number - numbers_og[start_of_current_slope].recognized_number slope_og = (
# ) / ( numbers_og[-1].recognized_number - numbers_og[start_of_current_slope].recognized_number
# numbers_og[-1].create_at/1000 - numbers_og[start_of_current_slope].create_at/1000 ) / (
# ) numbers_og[-1].create_at/1000 - numbers_og[start_of_current_slope].create_at/1000
# # b = y - a x )
# intercept_og = numbers_og[-1].recognized_number - slope_og * numbers_og[-1].create_at # b = y - a x
intercept_og = numbers_og[-1].recognized_number - slope_og * numbers_og[-1].create_at/1000
# slope_ng = ( slope_ng = (
# numbers_ng[-1].recognized_number - numbers_ng[0].recognized_number numbers_ng[-1].recognized_number - numbers_ng[0].recognized_number
# ) / ( ) / (
# numbers_ng[-1].create_at/1000 - numbers_ng[0].create_at/1000 numbers_ng[-1].create_at/1000 - numbers_ng[0].create_at/1000
# ) )
# intercept_ng = numbers_ng[-1].recognized_number - slope_ng * numbers_ng[-1].create_at intercept_ng = numbers_ng[-1].recognized_number - slope_ng * numbers_ng[-1].create_at/1000
fig, ax = plt.subplots() # a1 x + b1 = a2 x + b2
# (a1 - a2) x = b2 - b1
# x = (b2 - b1) / (a1 - a2)
intercept_x = (intercept_ng - intercept_og) / (slope_og - slope_ng)
intercept_y = slope_ng * intercept_x + intercept_ng
fig, ax = plt.subplots(figsize=(15, 6))
ax.plot(*xy(numbers_og)) ax.plot(*xy(numbers_og))
ax.plot(*xy(numbers_ng)) ax.plot(*xy(numbers_ng))
# abline(slope_og, -120000) abline(slope_og, intercept_og, numbers_og[-1].create_at/1000, intercept_x)
# abline(slope_ng, -500000) abline(slope_ng, intercept_ng, numbers_ng[-1].create_at/1000, intercept_x)
plt.show()
plt.margins(x=0)
locs = np.append(
np.linspace(numbers_og[0].create_at/1000, numbers_ng[0].create_at/1000, 6),
[numbers_ng[-1].create_at/1000, intercept_x]
)
# func = lambda x: x < intercept_x
# locs = np.append(locs[func(locs)][:-1], [intercept_x])
labels = ["{:%Y-%m-%d}".format(datetime.fromtimestamp(x, timezone.utc)) for x in locs]
plt.xticks(locs, labels)
plt.grid()
# plt.show()
plt.savefig("plot.png")
return 0
if __name__ == "__main__": if __name__ == "__main__":
main() import sys
sys.exit(main())

View file

@ -5,6 +5,7 @@ from numberdealers import parse_numberdealers, times, report_errors, numbers_per
def main(): def main():
channel = sys.argv[1]
numbers, errors = parse_numberdealers.parse(sys.stdin) numbers, errors = parse_numberdealers.parse(sys.stdin)
if numbers == [] and errors == []: if numbers == [] and errors == []:
@ -14,21 +15,20 @@ def main():
print("No valid number messages!") print("No valid number messages!")
return return
print(f"Checked from {numbers[0].recognized_number} up to {numbers[-1].recognized_number}") print(f"##### Checked ~{channel} from {numbers[0].recognized_number} up to {numbers[-1].recognized_number}")
print() if "--no-errors" not in sys.argv[2:]:
report_errors.report_errors(errors) print()
print() report_errors.report_errors(errors)
print("---") if "--no-times" not in sys.argv[2:]:
print() print()
print("```") print("```")
times.report_times(numbers) times.report_times(numbers)
print("```") print("```")
print() if "--no-users" not in sys.argv[2:]:
print("---") print()
print() print("```")
print("```") numbers_per_user.report_users(numbers)
numbers_per_user.report_users(numbers) print("```")
print("```")
if __name__ == "__main__": if __name__ == "__main__":