[geos-devel] [GEOS] #323: A point interpolated from a line does not always intersect the same line

GEOS geos-trac at osgeo.org
Thu Jan 28 17:49:17 EST 2010


#323: A point interpolated from a line does not always intersect the same line
--------------------------+-------------------------------------------------
 Reporter:  starsareblue  |       Owner:  geos-devel at lists.osgeo.org
     Type:  defect        |      Status:  new                       
 Priority:  major         |   Milestone:                            
Component:  Default       |     Version:  3.2.0                     
 Severity:  Unassigned    |    Keywords:                            
--------------------------+-------------------------------------------------
 Using the shapely Python bindings to GEOS, the following assertion fails.

 from shapely.wkt import loads

 m = loads("""MULTILINESTRING ((336700.9845364579814486
 1738886.2085458301007748, 343584.4261802079854533
 1738555.8341833299491554), (343584.4261802079854533
 1738555.8341833299491554, 350539.0000000000000000
 1736135.0000000000000000), (350539.0000000000000000
 1736135.0000000000000000, 353532.0000000000000000
 1733674.0000000000000000))""")

 p = loads("""POINT (344705.7704083333374001 1737219.7509562498889863)""")

 assert m.intersects(m.interpolate(m.project(p))) == True

-- 
Ticket URL: <http://trac.osgeo.org/geos/ticket/323>
GEOS <http://geos.refractions.net/>
GEOS (Geometry Engine - Open Source) is a C++ port of the Java Topology Suite (JTS).


More information about the geos-devel mailing list