<div dir="ltr">Hello,<div><br></div><div>I am running into an interesting version of the bad topology error. </div><div>During a raster/vector overlay I receive the following error: "GEOSIntersects: TopologyException"</div><div><br></div><div>I figure I have some bad topologies in my vector dataset. I run the following code and receive no records. I don't have any invalid topologies.</div><div><div><br></div><div>with data as</div><div>(</div><div>SELECT <a href="http://sgl.id">sgl.id</a> as sample_geog_level_id, <a href="http://gi.id">gi.id</a> as geog_instance_id, gi.label as geog_instance_label, gi.code as geog_instance_code, ST_IsValid(bound.geog::geometry) as valid</div><div>FROM sample_geog_levels sgl</div><div>inner join geog_instances gi on <a href="http://sgl.id">sgl.id</a> = gi.sample_geog_level_id</div><div>inner join boundaries bound on bound.geog_instance_id = <a href="http://gi.id">gi.id</a></div><div>WHERE <a href="http://sgl.id">sgl.id</a> = 1043</div><div>)</div><div>SELECT sample_geog_level_id</div><div>FROM data</div><div>WHERE valid = FALSE</div></div><div><br></div><div>However whenever I run this I receive the invalid topologies error</div><div><br></div><div><div>WITH geographic_boundaries AS</div><div>(</div><div>SELECT <a href="http://sgl.id">sgl.id</a> as sample_geog_level_id, <a href="http://gi.id">gi.id</a> as geog_instance_id, gi.label as geog_instance_label, gi.code as geog_instance_code, bound.geog::geometry as geom</div><div>FROM sample_geog_levels sgl</div><div>inner join geog_instances gi on <a href="http://sgl.id">sgl.id</a> = gi.sample_geog_level_id</div><div>inner join boundaries bound on bound.geog_instance_id = <a href="http://gi.id">gi.id</a></div><div>WHERE <a href="http://sgl.id">sgl.id</a> = 1043</div><div>)</div><div>SELECT p.geog_instance_id as geog_id, p.geog_instance_label as place_name, p.geog_instance_code as place_code, ST_Clip(r.rast, 1, p.geom, True) as rast</div><div>FROM geographic_boundaries p inner join landcover.glc2000_tp  r on ST_Intersects(r.rast, 1, p.geom)</div></div><div><br></div><div>I also went through the process of running ST_MakeValid(geom) on my data into a new table and still encounter the error. What am I missing? This particular error also spits out an x,y location that changes</div><div>-113.02232164290001 43.120535816300048</div><div>-126.41517874290001 47.584821516300046</div><div>-122.12053589950004 48.281250085500034<br></div></div>