Okay... I'm not so good at maths (didn't even finish our version of high school) and is really tired right now, but I need to orbit four objects elliptically around the centroid and have been stuck for a while, but it's not like I haven't tried or come anywhere by myself. This is what I've come up with so far:
from math import cos,sin,pi,sqrt,atan2
def orbit(p, deg, pivot=(0.32563325, 0.123498),ellipse=(0.5, 0.743992)):
# p = current (x,y)-position of current object,
# Pivot = (x,y)-position of point to orbit around,
# retrieved by centroid((x1,y1),(x2,y2),(x3,y3),(x4,y4))
# Ellipse = (width,height) of ellipse to rotate around
# retrieved by ellipse((x1,y1),(x2,y2),(x3,y3),(x4,y4))
px,py = pivot
if ellipse == (0.0,0.0):
o = polar(p[0]-px,p[1]-py)[0]
xr,yr = 1.0,1.0
else:
ew,eh = ellipse
if ew < eh:
o = eh/2 # Distance to the point most far away from the middle of the ellipse (Outer radius)
xr = ew/eh # Horizontal ratio of ellipse so we can move it back properly into ellipse after performing circular orbit
yr = 1.0 # Verical movement will not be affected because it's on the large axis
else:
o = ew/2
xr = 1.0
yr = eh/ew
x,y = p[0]-px,p[1]-py # Subtract pivot's position (that's the the point I want to orbit around)
d,a = polar(x,y) # Get angle
x,y = rect(o,a+deg) # Move point as far away from middle as it will be as most and make circular orbit around pivot by deg degrees
x,y = x*xr,y*yr # Move points back to elliptic shape by multiplying positions with according ratio <--- I guess it's here something goes wrong
x,y = x+px,y+py # Move point back to original position by adding pivot's positions
return x,y
# Other functions
def centroid(*points):
x,y = izip(*points)
return (sum(x) / len(points), sum(y) / len(points))
def ellipse(*points):
x,y = izip(*points)
xd,yd = max(x)-min(x),max(y)-min(y)
return (xd,yd)
def polar(x,y):
d = sqrt(x**2 + y**2)
a = atan2(y,x) * (180.0/pi)
return d, a
def rect(d,a):
x = d * cos(a*pi/180.0)
y = d * sin(a*pi/180.0)
return x, y
If i use ellipse=(0.0,0.0) and orbit everything in an "ordinary" circle instead of an elliptic shape it works just fine, so I guess it's when I try to multiply the x/y-position with the ratio ew/eh or eh/ew I do something wrong but I can't figure out what at the moment.
I am a bit tired and will try and get some sleep now and see if I can solve it tomorrow but some helpt would really be appriciated here.
This encapsulates the math for each Ellipse, and shows sample usage; you will have to add code for whatever output you want.