To assume compatibility with my survey software, I would like to lighten geometries of my linestrings by cutting them every 1000 points without creating or interpolating points.
I transformed a function found on the net (resample) but it cuts only the first 1000 points...
Does someone knows how to loop it correctly?
CREATE OR REPLACE FUNCTION simplify_npoints(inGeom geometry, maxPoints integer) RETURNS geometry AS $$
DECLARE nPoints integer; outGeom geometry; fraction float; points geometry;
nPoints := ST_NPoints(inGeom); outGeom := inGeom;
IF nPoints > maxPoints THEN FOR i IN 1..maxPoints LOOP points := array_append(points, ST_PointN(inGeom,i)); END LOOP; outGeom := ST_MakeLine(points);
END; $$ LANGUAGE plpgsql IMMUTABLE;