Mercurial > hg > orthanc-stone
view Resources/Computations/IntersectSegmentAndVerticalLine.py @ 482:f58fe38c8c04 bgo-commands-codegen
Ongoing work on codegen: ts and cpp enum and struct writing seem to be OK. No file write yet
author | bgo-osimis |
---|---|
date | Thu, 14 Feb 2019 20:58:42 +0100 |
parents | 35c2b85836ce |
children |
line wrap: on
line source
#!/usr/bin/env python from sympy import * # Intersection between the 2D line segment (prevX,prevY)-(curX,curY) and the # vertical line "x = x0" using homogeneous coordinates prevX, prevY, curX, curY, x0 = symbols('prevX prevY curX curY x0') p1 = Matrix([prevX, prevY, 1]) p2 = Matrix([curX, curY, 1]) l1 = p1.cross(p2) h1 = Matrix([x0, 0, 1]) h2 = Matrix([x0, 1, 1]) l2 = h1.cross(h2) a = l1.cross(l2) pprint(a / a[2])