Hi,
I'm building a topology from a big set of features, around 100k with 37
million vertexes. Since loading the whole set at once with toTopoGeom
wasn't feasible,
I partitioned my data country by country first, then I tried to load
polygon by polygon. No matter the method, after around 25k polygons I
get this error:
DBD::Pg::db do failed: ERROR: GEOSContains: TopologyException: side
location conflict at 0.097239999999999993 48.088819999999998
CONTEXT: SQL statement "UPDATE eu_geom_base_topo.edge_data SET
left_face = CASE WHEN left_face = 25444 THEN 25612 ELSE left_face END,
right_face = CASE WHEN right_face = 25444 THEN 25612 ELSE right_face END
WHERE ( left_face = 25444 OR right_face = 25444) AND NOT edge_id = ANY
('{87330,87328,87153,87155,87036,87152,87327,87328,87174,87245,87246,87279,87283,87278,87282,87167,87168,87170,87018,87020,87052,87054,87055,87247,87248,87250,87301,87303,87293,87226,87292,87329}')
AND NOT
ST_Contains('0103000020E610000001000000C30B0000B858518369189E3F336DFFC.........'::geometry,
geom)"
PL/pgSQL function "_st_addfacesplit" line 114 at instruction EXECUTE
SQL statement "SELECT topology._ST_AddFaceSplit(atopology,
newedge.edge_id, newedge.left_face, false)"
PL/pgSQL function "st_addedgemodface" line 506 at instruction SQL
PL/pgSQL function "topogeo_addlinestring" line 124 at assignment
SQL statement "SELECT array_cat(edges, array_agg(x)) FROM ( select
topology.TopoGeo_addLinestring(atopology, rec.geom, tol) as x ) as foo"
PL/pgSQL function "topogeo_addpolygon" line 24 at assignment
PL/pgSQL function "totopogeom" line 100 at loop FOR su righe di SELECT
at ./processWorkDb.pl line 797.
Had to truncate the geometry in ST_contains because of the message size
limit, the hexdump of the geometry is around 200k.
The instruction that fails is something like:
UPDATE eu_geom_base
SET topogeom = toTopoGeom(geom, 'eu_geom_base_topo', 1)
WHERE eu_geom_base.named_place_id= 20077019
It doesn't fail with just a specific geometry, this example is the last
one I found.
If I try to issue it in pgadmin it keeps failing, so it's not something
related to the loader script.
I tried to look for the error, it seems it's failing at this line
result = GEOSContains( g1, g2);
of postgis/lwgeom_geos.c
<http://postgis.net/docs/doxygen/2.1/d4/d06/postgis_2lwgeom__geos_8c_source.html>.
I'd guess it's geos throwing a topologyexception because of the geometry.
My postgis full version is: "POSTGIS="2.0.1 r9979"
GEOS="3.3.6-CAPI-1.7.6" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL
1.9.2, released 2012/10/08" LIBXML="2.7.8" TOPOLOGY RASTER"
Would upgrade to postgis 2.0.2 help for this issue?
I'm loading geometries with a "default" 1 tolerance, would play with
this parameter help to solve the issue? Would it dissolve geometries?
Thanks for any feedback.
Regards,
Paolo
--
Paolo Crosato
_______________________________________________
postgis-users mailing list
[email protected]
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users
