about summary refs log tree commit diff
diff options
context:
space:
mode:
-rw-r--r--NOTES.md33
-rw-r--r--convert_koordinate.py52
-rw-r--r--main2.py4
3 files changed, 87 insertions, 2 deletions
diff --git a/NOTES.md b/NOTES.md
new file mode 100644
index 0000000..e4e002f
--- /dev/null
+++ b/NOTES.md
@@ -0,0 +1,33 @@
+# Satellite Computation Notes
+
+## TODO:
+
+- speed
+- Kugelkoordinatensysteme (räumliche Polarkoordinaten)
+  * helferklasese
+
+## Notes:
+
+# Umrechnung kartesian - polar
+
+koordinatentripel (r, theta, phi)
+
+- r -> alle positiven Reellen Werte
+- theta -> Intervall [0, pi] bzw. [0, 180°]
+- phi -> Intervall [-pi, pi] bzw. [-180°, 180°] oder [0, 2pi] bzw. [0, 360°]
+
+Umrechnung polar ->  kartesian:
+
+- x = r * sin(theta) * cos(phi)
+- y = r * sin(theta) * sin(phi)
+- z = r * sin(theta)
+
+Umrechnung kartesian -> polar:
+
+- r = sqrt(x^2 + y^2 + z^2)
+- theta = arccos({z}over{sqrt(x^2 + y^2 + z^2)})
+- phi =
+  - wenn x > 0, arctan({y}over{x})
+  - wenn x = 0, sgn(y){pi}over{2}
+  - wenn x < 0 And y >= 0, arctan({y}over{x})+pi
+  - wenn x < 0 And y < 0, arctan({y}over{x})-pi
diff --git a/convert_koordinate.py b/convert_koordinate.py
new file mode 100644
index 0000000..b5e52e3
--- /dev/null
+++ b/convert_koordinate.py
@@ -0,0 +1,52 @@
+# import
+import math
+import numpy as np
+import matplotlib.pyplot as plt
+
+# pol[x, y, z]
+x = 1
+y = 2
+z = 3
+pol = [x, y, z]
+
+# define r, theta, phi
+r = 0
+theta = 0
+phi = 0
+kart = [0, 0, 0]
+
+def poltokart(pol):
+    # split up list
+    x = pol[0]
+    y = pol[1]
+    z = pol[2]
+
+    # radius
+    r = np.sqrt(np.power(x, 2) + np.power(y, 2) + np.power(z, 2))
+
+    # theta
+    a = np.sqrt(x^2 + y^2 + z^2)
+    b = ((z) / a)
+    b = b * math.pi / 180
+    theta = np.arccos(b)
+
+    # phi
+    if x > 0:
+        phi = np.arctan(y/x)
+    elif x == 0:
+        phi = np.sign(y)*(math.pi/2)
+    elif x < 0 and y >= 0:
+        phi = np.arctan(y/x) + math.pi
+    elif x < 0 and y < 0:
+        phi = np.arctan(y/x) - math.pi
+
+    # write to cartesian list
+    kart[0] = r
+    kart[1] = theta
+    kart[2] = phi
+
+# run
+poltokart(pol)
+
+print("{:<15}{:<60}".format("polar:", str(pol) ))
+print("{:<15}{:<60}".format("cartesian:", str(kart) ))
diff --git a/main2.py b/main2.py
index 62a4426..7745d86 100644
--- a/main2.py
+++ b/main2.py
@@ -87,7 +87,7 @@ class TLE:
 start = timer()
 
 # controll values
-category = "iridium"
+category = "argos"
 globalScale = 1
 satSize = 0.5
 orbitSubDivs = 256
@@ -95,7 +95,7 @@ resolution = 100        # get position of sat each x frames
 threshold = 0.001
 
 # if internet connection available:
-TLE.download(category)
+# TLE.download(category)
 
 # define
 sce = bpy.context.scene