#! /usr/bin/env python

"""
Usage:
python mecartor-projection <latitude> <longitude>
Example Bundestags Kooridnaten: 52°31'7.22" N 13°22'13.69" E
./mecartor-projection.py 52.5186729836 13.3704687765
x = 5.55928
y = 18.86398
"""
from math import log, tan, radians, pi
import sys
from typing import Final

# Sternwarte Peterberg WGS84 (49° 34′ N, 7° 0′ E)
MAP_CENTER_LONGITUDE: Final[float] = radians(7)

RADIUS: Final[float] = 50  # cm

latitude = float(sys.argv[1])  # Degree
longitude = float(sys.argv[2])  # Degree

# Hilfe-Variablen
lat = radians(latitude)
lon = radians(longitude)
# Kartesische Koordinaten
x: float = RADIUS * (lon - MAP_CENTER_LONGITUDE)
y: float = RADIUS * log(tan(pi/4) + lat/2)

print(f"x = {x:.5f}")
print(f"y = {y:.5f}")