GeosIntersects Error

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

GeosIntersects Error

haynes
Hello,

I am running into an interesting version of the bad topology error. 
During a raster/vector overlay I receive the following error: "GEOSIntersects: TopologyException"

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.

with data as
(
SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label, gi.code as geog_instance_code, ST_IsValid(bound.geog::geometry) as valid
FROM sample_geog_levels sgl
inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
inner join boundaries bound on bound.geog_instance_id = gi.id
WHERE sgl.id = 1043
)
SELECT sample_geog_level_id
FROM data
WHERE valid = FALSE

However whenever I run this I receive the invalid topologies error

WITH geographic_boundaries AS
(
SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label, gi.code as geog_instance_code, bound.geog::geometry as geom
FROM sample_geog_levels sgl
inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
inner join boundaries bound on bound.geog_instance_id = gi.id
WHERE sgl.id = 1043
)
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
FROM geographic_boundaries p inner join landcover.glc2000_tp  r on ST_Intersects(r.rast, 1, p.geom)

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
-113.02232164290001 43.120535816300048
-126.41517874290001 47.584821516300046
-122.12053589950004 48.281250085500034

_______________________________________________
postgis-users mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/postgis-users