Comment #6 on issue 1121 by [email protected]: Add Azimut
(North-Orientation) to sympy.geometry.Polygon
http://code.google.com/p/sympy/issues/detail?id=1121
A property of polygon object.
@property
def azimuts(self):
"""By Christian Karrie"""
def tarea(a, b, c):
return (b[0] - a[0])*(c[1] - a[1]) - (c[0] - a[0])*(b[1] - a[1])
def isright(a, b, c):
return bool(tarea(a, b, c) <= 0)
# Determine orientation of points
cw = isright(self.vertices[-1], self.vertices[0], self.vertices[1])
ret = {}
for i in xrange(0, len(self.vertices)):
a,b,c = self.vertices[i-2], self.vertices[i-1], self.vertices[i]
ang = Line.angle_between(Line(b, Point(1000000000,b[1])),
Line(b, c))
if cw ^ isright(a, b, c):
ret[b] = 2*S.Pi - ang
else:
ret[b] = ang
return ret
--
You received this message because this project is configured to send all
issue notifications to this address.
You may adjust your notification preferences at:
https://code.google.com/hosting/settings
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sympy-issues.
For more options, visit https://groups.google.com/groups/opt_out.
