Comments on: Simplifying a map layer using PostGIS topology

By: Wayne

Wayne — Wed, 10 Apr 2013 17:21:09 +0000

I ran into an issue with an exception being thrown using the vanilla code in the blog post. I found that it was resolved by changing the order of the lines that calculate tol and the if condition that throws an error if tol is 0. This forces the code to try tolerance 0 first, then thrown an error if that fails.

SQL/MM Spatial exception is the exception that gets thrown.

I changed the code to the following:
IF tol = 0 THEN
RAISE EXCEPTION ‘%’, SQLERRM;
END IF;
tol := round( (tol/2.0) * 1e8 ) / 1e8; — round to get to zero quicker

Needless to say simplifying with tolerance 0 is not supported.

By: strk

strk — Tue, 19 Mar 2013 14:16:24 +0000

Please use the bug tracker for this: http://trac.osgeo.org/postgis

By: Fabian Zeindl

Fabian Zeindl — Tue, 19 Mar 2013 14:05:12 +0000

when trying to use my own data, I get the following error:

UPDATE rayons SET topogeom = toTopoGeom(geom, ‘france_dept_topo’, 1);
ERROR: lwpoly_from_lwlines: shell must be closed
CONTEXT: SQL statement “WITH ids as ( select row_number() over () as seq, edge from unnest(‘{-160,-104,45,-46,41,-44,43,-53,50,-49,47,-67,66,64,-57,55,-61,60,-84,70,-73,72,79,-76,74,-102,101,-87,85,-96,95,-112,98,-100,93,-92,90,-121,120,-125,124,-119,117,-128,113,-109,107,115,-132,-131,-115,116,-108,114,-127,129,-39,-158,159,-153,155,-137,-138,147,-141,144,-145,152,-148,149,-133,134}’::int[] ) u(edge) ), edges AS ( select CASE WHEN i.edge < 0 THEN ST_Reverse(e.geom) ELSE e.geom END as g FROM ids i left join france_dept_topo.edge_data e ON(e.edge_id = abs(i.edge)) ORDER BY seq) SELECT ST_MakePolygon(ST_MakeLine(g.g)) FROM edges g;"
PL/pgSQL function _st_addfacesplit(character varying,integer,integer,boolean) line 34 at EXECUTE statement
SQL statement "SELECT topology._ST_AddFaceSplit(atopology, -newedge.edge_id, newedge.left_face, false)"
PL/pgSQL function st_addedgemodface(character varying,integer,integer,geometry) line 513 at SQL statement
PL/pgSQL function topogeo_addlinestring(character varying,geometry,double precision) 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(character varying,geometry,double precision) line 24 at assignment
PL/pgSQL function totopogeom(geometry,character varying,integer,double precision) line 100 at FOR over SELECT rows

any ideas?

By: On the fly simplification of topologically defined geometries « Strk's Blog

On the fly simplification of topologically defined geometries « Strk's Blog — Fri, 08 Mar 2013 09:05:48 +0000

[...] saw how to get a topologically consistent simplified version of a full layer, but that method isn’t fast enough for on-demand usage. Also, we saw how to perform a fast [...]

By: strk

strk — Wed, 27 Feb 2013 15:59:51 +0000

I suggest subscribing to the postgis users mailing list for this kind of help

By: Francis Seong

Francis Seong — Wed, 27 Feb 2013 11:48:39 +0000

I run :
SELECT st_astext(topogeom::geometry) FROM idns_link_0 WHERE id(topogeom) = #table.id

result is :
“MULTILINESTRING((14128213.2930175 4505878.66922969,14128214.0420175 4505877.92022969))”

I don’t understant this result.

By: strk

strk — Tue, 26 Feb 2013 09:38:00 +0000

The message is telling you that by dropping that edge it will be impossible for one TopoGeometry object to retain the same shape.

You can see the shape of the TopoGeometry using:

SELECT topogeom::geometry FROM #table WHERE id(topogeom) = #number

By: Francis Seong

Francis Seong — Tue, 26 Feb 2013 04:24:39 +0000

I want to reduce edges and nodes. I run “select st_remedgemodface(‘topo’, edge_id)”, I got follow error:

ERROR: TopoGeom #number in layer 1 (#table.topogeom) cannot be represented dropping edge #edge_id

By: strk

strk — Tue, 23 Oct 2012 07:21:38 +0000

Use: SHOW search_path; to get the current value and SET search_path TO x,y,z for setting the value.

By: peas

peas — Tue, 23 Oct 2012 03:30:15 +0000

quick update: I was able to create a topology by calling topology.CreateTopology()