Mercurial > hg > orthanc-stone

comparison Resources/Computations/IntersectSegmentAndHorizontalLine.py @ 132:35c2b85836ce wasm

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix rtstruct projections
author Sebastien Jodogne <s.jodogne@gmail.com>
date Fri, 17 Nov 2017 18:01:31 +0100
parents
children
comparison
equal deleted inserted replaced
131:3e6163a53b16 132:35c2b85836ce
1 #!/usr/bin/env python
2
3 from sympy import *
4
5 # Intersection between the 2D line segment (prevX,prevY)-(curX,curY) and the
6 # horizontal line "y = y0" using homogeneous coordinates
7
8 prevX, prevY, curX, curY, y0 = symbols('prevX prevY curX curY y0')
9
10 p1 = Matrix([prevX, prevY, 1])
11 p2 = Matrix([curX, curY, 1])
12 l1 = p1.cross(p2)
13
14 h1 = Matrix([0, y0, 1])
15 h2 = Matrix([1, y0, 1])
16 l2 = h1.cross(h2)
17
18 a = l1.cross(l2)
19
20 #pprint(cse(a/a[2], symbols = symbols('a b')))
21 pprint(a / a[2])