Hi all, I've got a raster table and a polygon table, and I am trying to
get the intersection between the two. A bit of the SQL is
SELECT studyid
, st_intersection(a.the_geom_2926, r.rast) AS geomval
FROM (
SELECT studyid
, the_geom_2926
FROM buffers
-- test on one location
WHERE studyid = 'pat1020'
) AS a
, rslope_usgs_gt8pct as r
WHERE st_intersects(a.the_geom_2926, r.rast)
I get an error message:
ERROR: GEOSIntersects: TopologyException: side location conflict at
1251947.8569092923 201701.02349212748
CONTEXT: PL/pgSQL function _st_intersects(geometry,raster,integer) line
22 at RETURN
PL/pgSQL function st_intersection(geometry,raster,integer) line 5 at
assignment
********** Error **********
ERROR: GEOSIntersects: TopologyException: side location conflict at
1251947.8569092923 201701.02349212748
SQL state: XX000
Context: PL/pgSQL function _st_intersects(geometry,raster,integer) line
22 at RETURN
PL/pgSQL function st_intersection(geometry,raster,integer) line 5 at
assignment
I have checked the validity of the buffer polygons and PostGIS reports
no errors.
Any idea why the intersection is not working?
PostgreSQL version:
PostgreSQL 9.4.5 on x86_64-unknown-linux-gnu, compiled by gcc (GCC)
4.4.7 20120313 (Red Hat 4.4.7-11), 64-bit
PostGIS version:
POSTGIS="2.2.0 r14208" GEOS="3.5.0-CAPI-1.9.0 r4084" PROJ="Rel. 4.7.1,
23 September 2009" GDAL="GDAL 1.11.2, released 2015/02/10"
LIBXML="2.7.6" RASTER
Thanks for any assistance!
--
-P.
**************************************************************
Philip M. Hurvitz, PhD | Research Assistant Professor | UW-CBE
Urban Form Lab | 1107 NE 45th Street, Suite 535 | Box 354802
University of Washington, Seattle, Washington 98195-4802, USA
[email protected] | http://gis.washington.edu/phurvitz
"What is essential is invisible to the eye." -de Saint-Exupéry
**************************************************************
_______________________________________________
postgis-users mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/postgis-users
