cleanpolys

cleanpolys.sql

-- Function: gkdb.repairgeom(geometry)

-- DROP FUNCTION gkdb.repairgeom(geometry);

CREATE OR REPLACE FUNCTION gkdb.repairgeom(polygon geometry)
  RETURNS geometry AS
$BODY$

DECLARE
	linework geometry;
	cleanpolygon geometry;
BEGIN
	cleanpolygon := polygon;
	LOOP
		linework := ST_SnapToGrid(
			ST_Collect(
				array(select ST_ExteriorRing((ST_DumpRings((ST_Dump(cleanpolygon)).geom)).geom)
			)
		), 0.001);

		IF linework IS NULL THEN RETURN  NULL; END IF;
	
		cleanpolygon := ST_SnapToGrid(ST_BuildArea(ST_Node(lines)), 0.001);

		IF cleanpolygon IS NULL THEN RETURN  NULL; END IF;
		
		EXIT WHEN ST_IsValid(cleanpolygon) IS TRUE;
	END LOOP;

	RETURN cleanpolygon;
	
--EXCEPTION
--	WHEN OTHERS THEN
--		RAISE NOTICE 'Error poly: %', ST_AsText(polygon);
--		RETURN null;

END;

$BODY$
  LANGUAGE plpgsql VOLATILE
  COST 100;
ALTER FUNCTION gkdb.repairgeom(geometry)
  OWNER TO postgres;