Table 4 and 6 of this paper https://contactproject.org/wp-content/uploads/2025/06/Earths-Evolving-Radio-Footprint_-A-Megawatt-Analysis-of-Anthropogenic-Emissions-into-Outer-Space-1900-2025.pdf turned into a graph

space.csv

space.py

"""
Earth's detectable radio footprint: TV power escaping into space (1920–2025).
Data: Table 3 cleaned (space.csv). Contact Project PDF:
Earth's Evolving Radio Footprint: A Megawatt Analysis (1900-2025).
"""
import pandas as pd
import numpy as np
import matplotlib
matplotlib.use("Agg")
import matplotlib.pyplot as plt

# space.csv: decade, seti_detectable_mw, total_anthropogenic_rf_mw, ...
df = pd.read_csv("space.csv")
df["seti_detectable_mw"] = pd.to_numeric(df["seti_detectable_mw"], errors="coerce")
df["total_anthropogenic_rf_mw"] = pd.to_numeric(df["total_anthropogenic_rf_mw"], errors="coerce")
decades = df["decade"].astype(str).tolist()
seti_mw = df["seti_detectable_mw"].values
total_rf_mw = df["total_anthropogenic_rf_mw"].values

x_pos = np.arange(len(decades))

fig, ax = plt.subplots(figsize=(10, 6), facecolor="#FFFFFF")
ax.set_facecolor("#FFFFFF")

ax.set_xlabel("Year / Decade", fontsize=11)
ax.set_ylabel("Estimated Average RF Power Escaping (MW)", fontsize=11)
ax.set_xticks(x_pos)
ax.set_xticklabels(decades, rotation=0)
ax.set_xlim(-0.5, len(decades) - 0.5)
ax.set_ylim(0, 750)  # total_anthropogenic_rf_mw max ~700
ax.plot(x_pos, seti_mw, color="#1E8449", marker="o", markersize=8, linewidth=2.5, label="Detectable with classical SETI (TV-like narrowband leakage)")
ax.plot(x_pos, total_rf_mw, color="#C0392B", marker="o", markersize=8, linewidth=2.5, label="Total anthropogenic RF emissions (all technologies)")

ax.grid(True, linestyle="-", alpha=0.3)
ax.set_axisbelow(True)

# Title and subtitle (figure coords) close to plot
fig.text(0.5, 0.94, "Earth’s Radio Emissions: Total vs SETI-Detectable Signals", ha="center", fontsize=16, fontweight="bold", color="#1a1a2e")
fig.text(0.5, 0.90, "The last value in the Drake Equation is how long a civilisation releases signals", ha="center", fontsize=10, color="#444", style="italic")
ax.legend(loc="upper left", framealpha=0.95, fontsize=10)

# Credit line (figure coordinates)
plt.figtext(
    0.01, 0.02,
    "Data contactproject.org  •  Graph by @iamreddave",
    ha="left",
    fontsize=9,
    style="italic",
    color="gray",
)

plt.tight_layout(rect=[0, 0.04, 1, 0.90])
out_path = "space_chart.png"
fig.savefig(out_path, dpi=200, bbox_inches="tight", facecolor="#FFFFFF")
print(f"Saved {out_path}")