zhetadelta/BoSLOO
1
0
Fork 0
Code Issues Pull requests Actions Packages Projects Releases Wiki Activity

started reporting lat/long on screen during render

Browse source
This commit is contained in:
Sakimori 2022-01-20 15:19:02 -05:00
parent 8d53f693fd
commit 454475ddec
2 changed files with 24 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
OrbitSim.py
View file

@ -52,6 +52,17 @@ class Planet:
if self.rotationPercentage >= 100.0: if self.rotationPercentage >= 100.0:
self.rotationPercentage -= 100.0 self.rotationPercentage -= 100.0
def sphericalToLatLong(self, theta, phi):
"""Converts theta and phi spherical coordinates to latitude and longitude. -> lat, long"""
rotRadian = self.rotationPercentage * 2 * math.pi
lat = -math.degrees(phi - (math.pi/2)) #positive lat is north, negative is south
long = theta - rotRadian #positive long is east, negative is west
if long < -math.pi:
long += math.pi*2
elif long > math.pi:
long -= math.pi*2
return (lat, math.degrees(long))
class DisplayPoint: class DisplayPoint:
"""A single point of any color""" """A single point of any color"""
def __init__(self, location, color): def __init__(self, location, color):

16
renderer.py
View file

@ -11,13 +11,19 @@ class Point:
self.vector = numpy.array([x, y, z]) self.vector = numpy.array([x, y, z])
def polar(self): def polar(self):
"""Converts the vector rectangular coordinates to polar coordinates."""
if self.vector[0] == 0: if self.vector[0] == 0:
self.vector[0] = 0.1 self.vector[0] = 0.1
if self.vector[2] == 0: if self.vector[2] == 0:
self.vector[2] = 0.1 self.vector[2] = 0.1
rho = math.sqrt(int(self.vector[0]) ** 2 + int(self.vector[1]) ** 2 + int(self.vector[2]) ** 2) rho = math.sqrt(int(self.vector[0]) ** 2 + int(self.vector[1]) ** 2 + int(self.vector[2]) ** 2)
theta = math.atan(self.vector[1]/self.vector[0]) theta = math.atan(self.vector[1]/self.vector[0]) #this has a range of -pi/2 to pi/2 but we need 0 to 2pi so more work needed
phi = math.acos(self.vector[2]/rho) phi = math.acos(self.vector[2]/rho)
if self.vector[0] < 0:
if self.vector[1] >= 0: #if x is positive, atan is fine. need to check if x is negative, first.
theta += math.pi
else:
theta -= math.pi
return [rho, theta, phi] return [rho, theta, phi]
def magnitude(self): def magnitude(self):
@ -197,14 +203,16 @@ class Camera:
#generate text #generate text
rho, theta, phi = sat.location.polar() rho, theta, phi = sat.location.polar()
theta = math.degrees(theta)
phi = math.degrees(phi)
if rho < self.target.radius: if rho < self.target.radius:
0 == 0 0 == 0
#textSurface, rect = font.render(f"Speed: {round(sat.velocity.magnitude())} m/s \nAltitude: {round(rho - target.radius)} m", False, (255,255,255)) rawLat, rawLong = self.target.sphericalToLatLong(theta, phi)
latString = f"Latitude: {round(rawLat,4)}⁰ N" if rawLat >= 0 else f"Latitude: {-round(rawLat,4)}⁰ S"
longString = f"Longitude: {round(rawLong,4)}⁰ E" if rawLong >= 0 else f"Longitude: {-round(rawLong,4)}⁰ W"
font.render_to(backSurface, (0,0), f"Speed: {round(sat.velocity.magnitude()/1000,3)} km/s", (255,255,255)) font.render_to(backSurface, (0,0), f"Speed: {round(sat.velocity.magnitude()/1000,3)} km/s", (255,255,255))
font.render_to(backSurface, (0,20), f"Altitude: {round((rho - self.target.radius)/1000)} km", (255,255,255)) font.render_to(backSurface, (0,20), f"Altitude: {round((rho - self.target.radius)/1000)} km", (255,255,255))
font.render_to(backSurface, (0,50), latString, (255,255,255))
font.render_to(backSurface, (0,70), longString, (255,255,255))
self.spriteGroup.update() self.spriteGroup.update()
self.spriteGroup.draw(backSurface) self.spriteGroup.draw(backSurface)