Mercurial > hg > orthanc-stone

view Resources/Computations/IntersectSegmentAndVerticalLine.py @ 1013:53cc787bd7bc toa2019092301

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
- Added an optimized ProjectPoint2 to CoordinateSystem3D. It has *not* replaced the ProjectPoint method because more tests need to be written. - ProjectPointOntoPlane2 is a faster version of GeometryToolbox::ProjectPointOntoPlane. Same remark as above. - DicomStructureSet.cpp now uses this optimized call.
author Benjamin Golinvaux <bgo@osimis.io>
date Mon, 23 Sep 2019 15:18:33 +0200
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])