mirror of https://github.com/miaowware/qrm2.git
143 lines
6.0 KiB
Python
143 lines
6.0 KiB
Python
"""
|
|
Grid extension for qrm
|
|
---
|
|
Copyright (C) 2019 Abigail Gold, 0x5c
|
|
|
|
This file is part of discord-qrm2 and is released under the terms of the GNU
|
|
General Public License, version 2.
|
|
"""
|
|
|
|
import math
|
|
|
|
import discord.ext.commands as commands
|
|
|
|
import common as cmn
|
|
|
|
|
|
class GridCog(commands.Cog):
|
|
def __init__(self, bot: commands.Bot):
|
|
self.bot = bot
|
|
|
|
@commands.command(name="grid", category=cmn.cat.maps)
|
|
async def _grid_sq_lookup(self, ctx: commands.Context, lat: str, lon: str):
|
|
'''Calculates the grid square for latitude and longitude coordinates,
|
|
with negative being latitude South and longitude West.'''
|
|
with ctx.typing():
|
|
grid = "**"
|
|
try:
|
|
latf = float(lat) + 90
|
|
lonf = float(lon) + 180
|
|
if 0 <= latf <= 180 and 0 <= lonf <= 360:
|
|
grid += chr(ord('A') + int(lonf / 20))
|
|
grid += chr(ord('A') + int(latf / 10))
|
|
grid += chr(ord('0') + int((lonf % 20)/2))
|
|
grid += chr(ord('0') + int((latf % 10)/1))
|
|
grid += chr(ord('a') + int((lonf - (int(lonf/2)*2)) / (5/60)))
|
|
grid += chr(ord('a') + int((latf - (int(latf/1)*1)) / (2.5/60)))
|
|
grid += "**"
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Maidenhead Grid Locator for {float(lat):.6f}, {float(lon):.6f}'
|
|
embed.description = grid
|
|
embed.colour = cmn.colours.good
|
|
else:
|
|
raise ValueError('Out of range.')
|
|
except ValueError as err:
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Error generating grid square for {lat}, {lon}.'
|
|
embed.description = str(err)
|
|
embed.colour = cmn.colours.bad
|
|
await ctx.send(embed=embed)
|
|
|
|
@commands.command(name="ungrid", aliases=['loc'], category=cmn.cat.maps)
|
|
async def _location_lookup(self, ctx: commands.Context, grid: str, grid2: str = None):
|
|
'''Calculates the latitude and longitude for the center of a grid square.
|
|
If two grid squares are given, the distance and azimuth between them is calculated.'''
|
|
with ctx.typing():
|
|
if grid2 is None or grid2 == '':
|
|
try:
|
|
grid = grid.upper()
|
|
loc = get_coords(grid)
|
|
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Latitude and Longitude for {grid}'
|
|
embed.colour = cmn.colours.good
|
|
|
|
if len(grid) >= 6:
|
|
embed.description = f'**{loc[0]:.5f}, {loc[1]:.5f}**'
|
|
embed.url = f'https://www.openstreetmap.org/#map=13/{loc[0]:.5f}/{loc[1]:.5f}'
|
|
else:
|
|
embed.description = f'**{loc[0]:.1f}, {loc[1]:.1f}**'
|
|
embed.url = f'https://www.openstreetmap.org/#map=10/{loc[0]:.1f}/{loc[1]:.1f}'
|
|
except Exception as e:
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Error generating latitude and longitude for grid {grid}.'
|
|
embed.description = str(e)
|
|
embed.colour = cmn.colours.bad
|
|
else:
|
|
radius = 6371
|
|
try:
|
|
grid = grid.upper()
|
|
grid2 = grid2.upper()
|
|
loc = get_coords(grid)
|
|
loc2 = get_coords(grid2)
|
|
# Haversine formula
|
|
d_lat = math.radians(loc2[0] - loc[0])
|
|
d_lon = math.radians(loc2[1] - loc[1])
|
|
a = math.sin(d_lat/2) ** 2 +\
|
|
math.cos(math.radians(loc[0])) * math.cos(math.radians(loc2[0])) *\
|
|
math.sin(d_lon/2) ** 2
|
|
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
|
|
d = radius * c
|
|
d_mi = 0.6213712 * d
|
|
|
|
# Bearing
|
|
y_dist = math.sin(math.radians(loc2[1]-loc[1])) * math.cos(math.radians(loc2[0]))
|
|
x_dist = math.cos(math.radians(loc[0])) * math.sin(math.radians(loc2[0])) -\
|
|
math.sin(math.radians(loc[0])) * math.cos(math.radians(loc2[0])) *\
|
|
math.cos(math.radians(loc2[1] - loc[1]))
|
|
bearing = (math.degrees(math.atan2(y_dist, x_dist)) + 360) % 360
|
|
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Great Circle Distance and Bearing from {grid} to {grid2}'
|
|
embed.description = f'**Distance:** {d:.1f} km ({d_mi:.1f} mi)\n**Bearing:** {bearing:.1f}°'
|
|
embed.colour = cmn.colours.good
|
|
except Exception as e:
|
|
embed = cmn.embed_factory(ctx)
|
|
embed.title = f'Error generating great circle distance and bearing from {grid} and {grid2}.'
|
|
embed.description = str(e)
|
|
embed.colour = cmn.colours.bad
|
|
await ctx.send(embed=embed)
|
|
|
|
|
|
def get_coords(grid: str):
|
|
if len(grid) < 3:
|
|
raise ValueError('The grid locator must be at least 4 characters long.')
|
|
|
|
if not grid[0:2].isalpha() or not grid[2:4].isdigit():
|
|
if len(grid) <= 4:
|
|
raise ValueError('The grid locator must be of the form AA##.')
|
|
if len(grid) >= 6 and not grid[5:7].isalpha():
|
|
raise ValueError('The grid locator must be of the form AA##AA.')
|
|
|
|
lon = ((ord(grid[0]) - ord('A')) * 20) - 180
|
|
lat = ((ord(grid[1]) - ord('A')) * 10) - 90
|
|
lon += ((ord(grid[2]) - ord('0')) * 2)
|
|
lat += ((ord(grid[3]) - ord('0')) * 1)
|
|
|
|
if len(grid) >= 6:
|
|
# have subsquares
|
|
lon += ((ord(grid[4])) - ord('A')) * (5/60)
|
|
lat += ((ord(grid[5])) - ord('A')) * (2.5/60)
|
|
# move to center of subsquare
|
|
lon += (2.5/60)
|
|
lat += (1.25/60)
|
|
return (lat, lon)
|
|
# move to center of square
|
|
lon += 1
|
|
lat += 0.5
|
|
return (lat, lon)
|
|
|
|
|
|
def setup(bot: commands.Bot):
|
|
bot.add_cog(GridCog(bot))
|