89 lines
2.9 KiB
Python
89 lines
2.9 KiB
Python
import itertools
|
|
import logging
|
|
from datetime import datetime
|
|
from math import pi, sqrt, sin, cos, atan2
|
|
|
|
import pytz
|
|
import requests
|
|
from timezonefinder import TimezoneFinder
|
|
|
|
from ..util import Context, Prediction
|
|
|
|
tf = TimezoneFinder(in_memory=True)
|
|
|
|
|
|
def night_on_iss(context: Context) -> Prediction:
|
|
"""
|
|
It is night if it is night on the ISS and it is currently orbiting above us.
|
|
"""
|
|
p = Prediction()
|
|
|
|
if not context.flat_earth:
|
|
iss_position = requests.get("http://api.open-notify.org/iss-now.json").json()["iss_position"]
|
|
the_iss = "The ISS"
|
|
iss_position_description = "on board the ISS"
|
|
else:
|
|
p.reasons.append("The ISS is (obviously) located in Hollywood")
|
|
the_iss = "Hollywood"
|
|
iss_position = {'latitude': 34.092808, 'longitude': -118.328659} # Hollywood
|
|
iss_position_description = "in the Hollywood studio"
|
|
|
|
phone_position = context.position
|
|
|
|
# Calculate ratio: a number between 0 and 1 saying how close we are to the ISS
|
|
distance = haversine(iss_position, phone_position)
|
|
max_distance = 40075 / 2 # the furthest you can be from any position is half of the earth's circumference
|
|
ratio = distance / max_distance
|
|
|
|
# We're in the same "timezone" as the ISS if we're on the same half of the earth
|
|
on_iss_time = ratio < 0.5
|
|
|
|
side = "same" if on_iss_time else "other"
|
|
p.reasons.append(f"{the_iss} is {int(distance)} km away, so we are on the {side} side of the earth.")
|
|
for i in itertools.count(1):
|
|
iss_tz = tf.closest_timezone_at(lng=float(iss_position["longitude"]),
|
|
lat=float(iss_position["latitude"]),
|
|
delta_degree=i)
|
|
if iss_tz is not None:
|
|
break
|
|
iss_time = datetime.now(pytz.timezone(iss_tz))
|
|
|
|
iss_night = iss_time.hour < 6 or iss_time.hour >= 22
|
|
|
|
# iss_night on_iss_time night
|
|
# 0 0 1
|
|
# 0 1 0
|
|
# 1 0 0
|
|
# 1 1 1
|
|
night = iss_night == on_iss_time
|
|
|
|
iss_time_description = "nighttime" if iss_night else "daytime"
|
|
time_description = "nighttime" if night else "daytime"
|
|
p.probability = float(night)
|
|
p.reasons.append(f"It is {iss_time_description} {iss_position_description}.")
|
|
p.reasons.append(f"Therefore, it must be {time_description} where we are.")
|
|
|
|
return p
|
|
|
|
|
|
def haversine(pos1, pos2):
|
|
"""
|
|
Distance between two GPS coordinates.
|
|
https://stackoverflow.com/a/18144531
|
|
"""
|
|
lat1 = float(pos1["latitude"])
|
|
long1 = float(pos1["longitude"])
|
|
lat2 = float(pos2["latitude"])
|
|
long2 = float(pos2["longitude"])
|
|
|
|
degree_to_rad = float(pi / 180.0)
|
|
|
|
d_lat = (lat2 - lat1) * degree_to_rad
|
|
d_long = (long2 - long1) * degree_to_rad
|
|
|
|
a = pow(sin(d_lat / 2), 2) + cos(lat1 * degree_to_rad) * cos(lat2 * degree_to_rad) * pow(sin(d_long / 2), 2)
|
|
c = 2 * atan2(sqrt(a), sqrt(1 - a))
|
|
km = 6367 * c
|
|
|
|
return km
|