Mercurial > hg > nsaunier > traffic-intelligence

view trafficintelligence/tests/moving_shapely.txt @ 1136:30171d4fd3df Tertuis-Ou-draogo/movingpy-issue-22-from-etienne-beauchamp-1584900380975

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
moving.py [Issue #22 from Etienne Beauchamp In annotationTool Correction]
author Tertuis Ouédraogo <tertuis95@gmail.com>
date Sun, 22 Mar 2020 18:06:22 +0000
parents cc5cb04b04b0
children
line wrap: on
line source

>>> from moving import *
>>> from shapely.geometry import Polygon
>>> from shapely.prepared import prep

>>> t1 = Trajectory([[0.5,1.5,2.5],[0.5,3.5,6.5]])
>>> poly = Polygon([[0,0],[4,0],[4,3],[0,3]])
>>> sub1, sub2 = t1.getTrajectoryInPolygon(poly)
>>> sub1
(0.500000,0.500000)
>>> sub1, sub2 = t1.getTrajectoryInPolygon(Polygon([[10,10],[14,10],[14,13],[10,13]]))
>>> sub1.length()
0
>>> sub1, sub2 = t1.getTrajectoryInPolygon(prep(poly))
>>> sub1
(0.500000,0.500000)
>>> t2 = t1.differentiate(True)
>>> sub1, sub2 = t1.getTrajectoryInPolygon(prep(poly), t2)
>>> sub1.length() == sub2.length()
True
>>> sub1
(0.500000,0.500000)
>>> sub2
(1.000000,3.000000)

>>> t1.proportionInPolygon(poly, 0.5)
False
>>> t1.proportionInPolygon(poly, 0.3)
True