More

analysis.py

@@ -312,4 +312,32 @@ def base_analysis():
     hours_per_model()
     per_month()
 
-base_analysis()
+
+def distribution():
+    x = duckdb.sql(
+        f"""
+INSTALL spatial;
+LOAD spatial;
+SELECT
+    epoch("end" - "start") / 60 / 60 AS "flying_time_k_hours"
+    , ST_Distance(ST_Point(from_lat,from_lon),ST_Point(to_lat,to_lon))
+FROM '{PATH}'
+USING SAMPLE 10%
+    """
+    ).fetchall()
+
+    y = [x[1] for x in x]
+    x = [x[0] for x in x]
+
+    plt.figure(figsize=FIG_SIZE)
+    plt.plot(x, y, '.')
+    plt.xlabel("Flying time (hours)")
+    plt.ylabel("Distance (feet)")
+    plt.grid(linestyle="dotted")
+    plt.tight_layout()
+    plt.savefig("results/dist.png", dpi=300)
+
+
+distribution()
+
+# base_analysis()