[postgis-users] problem with intersection

[dblasby at openplans.org]

Martin Davis,

Here's another test case for intersection() for your JTS testcase pile. 
This one's from the PostGIS list.

com.vividsolutions.jts.geom.TopologyException: no outgoing dirEdge found
[ (588928.6000000002, 2432123.7, NaN) ]
	at
com.vividsolutions.jts.geomgraph.DirectedEdgeStar.linkResultDirectedEdges(DirectedEdgeStar.java:241)
..



public class test {

	public static void main(String[] args)
	{
		try
		{
			WKTReader r = new WKTReader();
			Geometry g1 = r.read("POLYGON((588912.7 2432104.5,588927.7 2432124,
588928.9 2432123.6, 588923.7 2432132.3,588912.7 2432104.5))");
			Geometry g2 = r.read("POLYGON((588916.4 2432124.4, 588920.7
2432129.9,588928.6 2432123.7,588930.1 2432122.6,588916.4 2432124.4))");
			System.out.println("isvalid(g1) = "+ g1.isValid());
			System.out.println("isvalid(g2) = "+ g2.isValid());
			Geometry g3 = g1.intersection(g2);
			System.out.println("g3="+g3.toText());
		}
		catch(Exception e)
		{
			e.printStackTrace();
		}
	}
}

----------------------------------------------------------
This mail sent through IMP: https://webmail.limegroup.com/