StellarDominion/prototype/lunar-cross-visualization.py at master · lispmeister/StellarDominion · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
"""Scientifically accurate lunar laser cross visualization at exact naked-eye angular scales."""
# pylint: disable=invalid-name
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.patches import Circle

# ================================================
# Scientifically Accurate Lunar Laser Cross Visualization
# Exact naked-eye angular scales (no artistic enlargement)
# ================================================

fig = plt.figure(figsize=(14, 9), facecolor='black')
ax = fig.add_subplot(111)
ax.set_facecolor('black')

# Set angular coordinate system (degrees on sky)
# x = horizontal angular position (degrees)
# y = vertical angular position (degrees above horizon)
ax.set_xlim(-2.0, 2.0)
ax.set_ylim(0.0, 1.5)          # Moon low at ~15° but scaled for visual clarity
ax.set_aspect('equal')         # True angular proportions

# Remove axes spines and ticks for clean sky look
ax.axis('off')

# ==================== STARS ====================
np.random.seed(42)  # reproducible star field
N_STARS = 300
star_x = np.random.uniform(-2.0, 2.0, N_STARS)
star_y = np.random.uniform(0.0, 1.5, N_STARS)
ax.scatter(star_x, star_y, color='white', s=np.random.uniform(0.5, 3, N_STARS),
           alpha=np.random.uniform(0.6, 1.0, N_STARS), edgecolor='none')

# ==================== HORIZON ====================
ax.axhline(y=0.0, color='#1e3a5f', linewidth=6, zorder=1, label='Horizon')
ax.text(-1.95, 0.03, 'Horizon', color='white', fontsize=12, fontfamily='monospace')

# ==================== SCALE BAR (1°) ====================
ax.plot([-1.0, 0.0], [0.08, 0.08], color='white', linewidth=2)
ax.plot([-1.0, -1.0], [0.06, 0.10], color='white', linewidth=2)
ax.plot([0.0, 0.0], [0.06, 0.10], color='white', linewidth=2)
ax.text(-0.5, 0.12, '1°', color='white', fontsize=14, ha='center', fontfamily='monospace')
ax.text(-1.95, 0.15, 'Observers on Earth\n(looking up)', color='white', fontsize=11,
        ha='left', va='top', fontfamily='monospace')

# ==================== FULL MOON (exactly 0.5° diameter) ====================
MOON_CENTER_X = 0.0
MOON_CENTER_Y = 0.35
MOON_RADIUS = 0.25   # = 0.5° diameter

moon = Circle((MOON_CENTER_X, MOON_CENTER_Y), MOON_RADIUS,
              facecolor='#f8f1e9', edgecolor='#e5d9c8', linewidth=2, zorder=10)
ax.add_patch(moon)

# Slight lunar surface shading for realism
ax.add_patch(Circle((MOON_CENTER_X - 0.08, MOON_CENTER_Y + 0.08), 0.12,
                    facecolor='#e5d9c8', alpha=0.4, zorder=11))

# Label for Moon
ax.text(MOON_CENTER_X, MOON_CENTER_Y - 0.42, 'FULL MOON\n(exactly 0.5° apparent diameter)',
        color='white', fontsize=13, ha='center', va='top', fontfamily='monospace')

# ==================== EMERALD-GREEN LASER CROSS ====================
# Center of cross = lunar northern limb (top of Moon)
CROSS_CENTER_X = MOON_CENTER_X
CROSS_CENTER_Y = MOON_CENTER_Y + MOON_RADIUS   # exactly at top limb

# Cross total span = 0.7° (matches continental US angular width from Moon)
# Each arm: 0.7° long × 0.05° thick (matches README solid-angle calculation)
CROSS_HALF_SPAN = 0.35        # half of 0.7°
CROSS_HALF_THICKNESS = 0.025  # half of 0.05°

# Horizontal arm (filled rectangle in data coordinates)
ax.fill([CROSS_CENTER_X - CROSS_HALF_SPAN, CROSS_CENTER_X + CROSS_HALF_SPAN,
         CROSS_CENTER_X + CROSS_HALF_SPAN, CROSS_CENTER_X - CROSS_HALF_SPAN],
        [CROSS_CENTER_Y - CROSS_HALF_THICKNESS, CROSS_CENTER_Y - CROSS_HALF_THICKNESS,
         CROSS_CENTER_Y + CROSS_HALF_THICKNESS, CROSS_CENTER_Y + CROSS_HALF_THICKNESS],
        color='#00ff9d', zorder=20)

# Vertical arm (filled rectangle in data coordinates, symmetric 0.7° span)
ax.fill([CROSS_CENTER_X - CROSS_HALF_THICKNESS, CROSS_CENTER_X + CROSS_HALF_THICKNESS,
         CROSS_CENTER_X + CROSS_HALF_THICKNESS, CROSS_CENTER_X - CROSS_HALF_THICKNESS],
        [CROSS_CENTER_Y - CROSS_HALF_SPAN, CROSS_CENTER_Y - CROSS_HALF_SPAN,
         CROSS_CENTER_Y + CROSS_HALF_SPAN, CROSS_CENTER_Y + CROSS_HALF_SPAN],
        color='#00ff9d', zorder=20)

# Label for the cross
ax.text(CROSS_CENTER_X + 0.05, CROSS_CENTER_Y + CROSS_HALF_SPAN + 0.28,
        'EMERALD-GREEN LASER CROSS\n(~0.7° angular span from lunar north pole)\n'
        'Visible simultaneously across the entire continental US',
        color='#00ff9d', fontsize=14, ha='left', va='top', fontfamily='monospace',
        bbox={"boxstyle": "round,pad=0.5", "facecolor": 'black', "alpha": 0.7})

# Title
ax.text(0.0, 1.42, 'Scientifically Accurate Visualization\n'
                   'Lunar Laser Cross Projection from Northern Pole\n'
                   '(Exact naked-eye angular scales — no artistic enlargement)',
        color='white', fontsize=16, ha='center', fontfamily='monospace', weight='bold')

# ==================== SAVE / SHOW ====================
plt.tight_layout()
plt.savefig('lunar_cross_scientific_visualization_proper_cross.png', dpi=300, bbox_inches='tight')
plt.show()