Mercurial > hg > orthanc-stone
view Resources/Computations/IntersectSegmentAndHorizontalLine.py @ 1141:7681f3943748
Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to
the DoseGridScaling value if rescale* tags weren't present AND in the case of a
dose. This caused issues in dose files where all three tags (DoseGridScaling,
RescaleSlope and RescaleIntercept) were present, which previously led to the
DoseGridScaling tag to be ignored. Now, the rescale* tags are *ignored* in dose
files and DoseGridScaling is always taken into account.
author | Benjamin Golinvaux <bgo@osimis.io> |
---|---|
date | Fri, 08 Nov 2019 14:25:35 +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 # horizontal line "y = y0" using homogeneous coordinates prevX, prevY, curX, curY, y0 = symbols('prevX prevY curX curY y0') p1 = Matrix([prevX, prevY, 1]) p2 = Matrix([curX, curY, 1]) l1 = p1.cross(p2) h1 = Matrix([0, y0, 1]) h2 = Matrix([1, y0, 1]) l2 = h1.cross(h2) a = l1.cross(l2) #pprint(cse(a/a[2], symbols = symbols('a b'))) pprint(a / a[2])