<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#333399">
Interestingly, if I put a "LIMIT 1" right before the semicolon, the
query works as expected. But let it do more than one record and it
fails with that error! This smells like a bug to me! I've used the
IsValid() and all geometries in all records check out ("t"), and all
geometries are polygons. Please help! I'm at a standstill.<br>
Mark<br>
<br>
Mark Fenbers wrote:
<blockquote cite="mid4288F3F7.80802@noaa.gov" type="cite">
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
I'm consistently getting an error that I cannot make sense of. Here's
the query I run:<br>
<br>
<tt>SELECT 'ILN' AS "base radar",radid, <br>
npoints(<br>
Difference(<br>
base_geom,<br>
x_geom<br>
)<br>
)<br>
FROM MJF WHERE radid != 'ILN';<br>
<br>
</tt>Here's the error I get::<br>
<br>
<tt>NOTICE: TopologyException: no outgoing dirEdge found
(34380.5,521104)<br>
ERROR: GEOS difference() threw an error!<br>
</tt><br>
If I change my WHERE clause such that only one row is returned instead
(e.g., unnegate the equality), then I do NOT get the error. I get the
error whether I use Difference(), Intersection(), GeomUnion(), or the
like. The 'base_geom' contains a circular geometry returned from the
Buffer() function (to create a disk around point). The x_geom is a
rectangle that I know intersects the circle (base_geom). Any ideas how
to interpret this error?<br>
<br>
Mark<br>
<br>
<br>
<pre wrap="">_______________________________________________
postgis-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>
<a class="moz-txt-link-freetext" href="http://postgis.refractions.net/mailman/listinfo/postgis-users">http://postgis.refractions.net/mailman/listinfo/postgis-users</a>
</pre>
</blockquote>
</body>
</html>