DROP   TYPE IF EXISTS otm_al_grid;
CREATE TYPE otm_al_grid AS 
 (x                      DOUBLE PRECISION, 
  y                      DOUBLE PRECISION,
  is_in_area             BOOLEAN,
  is_near_axis           BOOLEAN,
  dist_to_border         DOUBLE PRECISION,
  dist_to_middle         DOUBLE PRECISION,
  dist_to_start          DOUBLE PRECISION,
  way_to_middle          INTEGER,
  way_to_start           INTEGER
);


--
-- Thanks to Gabor Farkas
-- https://gis.stackexchange.com/questions/56835/how-to-perform-sia-or-bezier-line-smoothing-in-postgis
--
CREATE OR REPLACE FUNCTION OTM_CreateCurve(geom geometry, percent int DEFAULT 40)
    RETURNS geometry AS
$$
DECLARE
    result text;
    p0 geometry;
    p1 geometry;
    p2 geometry;
    intp geometry;
    tempp geometry;
    geomtype text := ST_GeometryType(geom);
    factor double precision := percent::double precision / 200;
    i integer;
BEGIN
    IF percent < 0 OR percent > 100 THEN
        RAISE EXCEPTION 'Smoothing factor must be between 0 and 100';
    END IF;
    IF geomtype != 'ST_LineString' OR factor = 0 THEN
        RETURN geom;
    END IF;
    result := 'COMPOUNDCURVE((';
    p0 := ST_PointN(geom, 1);
    IF ST_NPoints(geom) = 2 THEN
        p1:= ST_PointN(geom, 2);
        result := result || ST_X(p0) || ' ' || ST_Y(p0) || ',' || ST_X(p1) || ' ' || ST_Y(p1) || '))';
    ELSE
        FOR i IN 2..(ST_NPoints(geom) - 1) LOOP
            p1 := ST_PointN(geom, i);
            p2 := ST_PointN(geom, i + 1);
            result := result || ST_X(p0) || ' ' || ST_Y(p0) || ',';
            tempp := ST_Line_Interpolate_Point(ST_MakeLine(p1, p0), factor);
            p0 := ST_Line_Interpolate_Point(ST_MakeLine(p1, p2), factor);
            intp := ST_Line_Interpolate_Point(
                ST_MakeLine(
                    ST_Line_Interpolate_Point(ST_MakeLine(p0, p1), 0.5),
                    ST_Line_Interpolate_Point(ST_MakeLine(tempp, p1), 0.5)
                ), 0.5);
            result := result || ST_X(tempp) || ' ' || ST_Y(tempp) || '),CIRCULARSTRING(' || ST_X(tempp) || ' ' || ST_Y(tempp) || ',' || ST_X(intp) || ' ' ||
            ST_Y(intp) || ',' || ST_X(p0) || ' ' || ST_Y(p0) || '),(';
        END LOOP;
        result := result || ST_X(p0) || ' ' || ST_Y(p0) || ',' || ST_X(p2) || ' ' || ST_Y(p2) || '))';
    END IF;
    RETURN ST_SetSRID(result::geometry, ST_SRID(geom));
END;
$$ LANGUAGE plpgsql;




CREATE OR REPLACE FUNCTION arealabel(myosm_id IN BIGINT,myway IN GEOMETRY,myname IN TEXT) RETURNS GEOMETRY AS $$

 DECLARE
  bbox            GEOMETRY;
  tmpway          GEOMETRY;
  retway          GEOMETRY;
  smoothway       GEOMETRY;
  tmpbox          GEOMETRY;
  tmppoint        GEOMETRY;
  centroid        GEOMETRY;
  tmplinestring   VARCHAR;
  linestring      VARCHAR;
  t               VARCHAR;
  rname           TEXT;
  areashape       VARCHAR:='unknown';
  x               DOUBLE PRECISION;
  y               DOUBLE PRECISION;
  x1              DOUBLE PRECISION;
  y1              DOUBLE PRECISION;
  xmin            DOUBLE PRECISION;
  xmax            DOUBLE PRECISION;
  ymin            DOUBLE PRECISION;
  ymax            DOUBLE PRECISION;
  i               INTEGER;
  j               INTEGER;
  q               INTEGER;
  r               INTEGER;
  i1              INTEGER;
  j1              INTEGER;
  imin            INTEGER;
  jmin            INTEGER;
  imax            INTEGER;
  jmax            INTEGER;
  p               DOUBLE PRECISION;
  k               DOUBLE PRECISION;
  f               DOUBLE PRECISION;
  g               DOUBLE PRECISION;
  m               DOUBLE PRECISION;
  n               DOUBLE PRECISION;
  changes         INTEGER;
  done            INTEGER;
  gridsize        CONSTANT INTEGER:=25;
  grid            otm_al_grid[625];
  gridpoint       otm_al_grid;
  targetpoint     otm_al_grid;
  gridwidth       DOUBLE PRECISION;
  gridheight      DOUBLE PRECISION;
  griddiag        DOUBLE PRECISION;
  maybehorizontal INTEGER:=0;
  lastindex       INTEGER;
  middleindex     INTEGER;
  startindex      INTEGER;
  endindex        INTEGER;
 BEGIN

--
-- get bounding box, adjust grid to same grid width for x and y
--
  done=0;
  bbox:=ST_SetSRID(ST_Envelope(myway),3857);
  xmin:=ST_XMin(bbox);xmax:=ST_XMax(bbox);
  ymin:=ST_YMin(bbox);ymax:=ST_YMax(bbox);
  x=xmax-xmin;y:=ymax-ymin;
  rname:=myname;
  gridwidth:=x/(gridsize-1);
  gridheight:=y/(gridsize-1);
  griddiag:=SQRT(gridwidth*gridwidth+gridheight*gridheight);
  retway:=NULL;
--
-- check for circles, sqares and rectangles in W-E-direction. 
-- (= any shape, wider than high or max 1.4x higher than wide and filling ists bbox>50%)
-- 
  IF ((((y/x)<1.4)) AND (ST_Area(myway)/(x*y)>0.5)) THEN
   maybehorizontal:=1;
   centroid=st_centroid(myway);
--
-- test if a simple horizontal line at the centroid is long enough
--
   IF (ST_Within(centroid,myway)) THEN
    x1:=ST_X(centroid);
    y:=ST_Y(centroid);
    WHILE  (ST_Within(ST_SetSRID(ST_MakePoint(x1+gridwidth/2,y),3857),ST_SetSRID(myway,3857))) LOOP
     x1:=x1+gridwidth/2;
    END LOOP;
    x:=ST_X(centroid);
    WHILE  (ST_Within(ST_SetSRID(ST_MakePoint(x-gridwidth/2,y),3857),ST_SetSRID(myway,3857))) LOOP
     x:=x-gridwidth/2;
    END LOOP;
    IF ((x1-x)/(xmax-xmin)>0.66)  THEN
     areashape:='simple';
     retway=ST_SetSRID(ST_MakeLine(ST_MakePoint(x,y),ST_MakePoint(x1,y)),3857);
    END IF;
   END IF;
  END IF;


--
-- initialize grid point with coordinates, check which grid point is in the area
--
  IF ( retway IS NULL ) THEN
   x:=xmin;i:=0;j:=0;
   WHILE (x<=xmax) LOOP
    y:=ymin;j:=0;
    WHILE (y<=ymax) LOOP
     gridpoint.x=x;
     gridpoint.y=y;
     tmpbox:=ST_SetSRID(ST_Envelope(ST_MakeLine(ST_MakePoint(x-gridwidth/2,y-gridheight/2), ST_MakePoint(x+gridwidth/2,y+gridheight/2))),3857);
     gridpoint.is_in_area=ST_Within(tmpbox,ST_SetSRID(myway,3857));
     gridpoint.is_near_axis:=false;
     gridpoint.dist_to_border:=0;
     gridpoint.dist_to_middle:=0;
     gridpoint.way_to_middle:=0;
     gridpoint.dist_to_start:=0;
     gridpoint.way_to_start:=0;
     grid[i*gridsize+j]:=gridpoint;
     y:=y+gridheight;
     j:=j+1;
    END LOOP; 
    x:=x+gridwidth; 
    i:=i+1;
   END LOOP;  
--
-- calc distance to borders
--
   FOR i IN 1..(gridsize-2) LOOP
    FOR j IN 1..(gridsize-2) LOOP
     gridpoint=grid[i*gridsize+j];
     IF (gridpoint.is_in_area) THEN
      gridpoint.dist_to_border=LEAST(grid[(i-1)*gridsize+j].dist_to_border+gridwidth,grid[(i)*gridsize+(j-1)].dist_to_border+gridheight);
     END IF;
     grid[i*gridsize+j]:=gridpoint;
    END LOOP;
   END LOOP;
   FOR i IN REVERSE (gridsize-2)..1 LOOP
    FOR j IN REVERSE (gridsize-2)..1 LOOP
     gridpoint=grid[i*gridsize+j];
     IF (gridpoint.is_in_area) THEN
      gridpoint.dist_to_border=LEAST(gridpoint.dist_to_border,LEAST(grid[(i+1)*gridsize+j].dist_to_border+gridwidth,grid[(i)*gridsize+(j+1)].dist_to_border+gridheight));
     END IF;
     grid[i*gridsize+j]:=gridpoint;
    END LOOP;
   END LOOP;
--
-- for horizontal: find the longest row in that grid, must be in the medium 1/3 of the area and at least 50% of the grid
--  
   i1:=0;j1:=0;m:=gridsize;g:=0;
   IF (maybehorizontal=1) THEN
    FOR j IN gridsize/3..(gridsize-2)*2/3 LOOP
     i:=1;
     WHILE NOT ((grid[i*gridsize+j].is_in_area) OR (i>gridsize-2)) LOOP i:=i+1; END LOOP;
     q=i;
     WHILE ((grid[i*gridsize+j].is_in_area) AND (i<=gridsize-1)) LOOP i:=i+1; END LOOP;
     r:=i-1;
     IF ((r-q>1) AND ((r-q>g) OR ((r-q=g) AND (ABS((j+1-(gridsize::FLOAT/2)))<m)))) THEN 
      g:=r-q;
      m:=ABS(j+1-(gridsize::FLOAT/2));
      startindex:=q*gridsize+j;
      endindex:=r*gridsize+j;
     END IF;
    END LOOP;
   END IF;
   IF (g>0.5*gridsize) THEN
    linestring:='LINESTRING(';
    gridpoint=grid[startindex];
    linestring:=linestring || gridpoint.x || ' ' || gridpoint.y || ',';
    gridpoint=grid[endindex];
    linestring:=linestring || gridpoint.x || ' ' || gridpoint.y || ')';
    rname:='h_' || rname;
    retway:=ST_GeomFromText(linestring,3857);
    areashape:='horizontal';
   END IF;
--
-- Horizontal labels are done now
--
   IF ( retway IS NULL ) THEN
--
-- mark points at a "slope" in all 4 directions
--
    middleindex:=0;p:=0;imin=gridsize;jmin=gridsize;imax:=-1;jmax:=-1;
    FOR i IN 1..(gridsize-2) LOOP
     FOR j IN 1..(gridsize-2) LOOP
      gridpoint=grid[i*gridsize+j];
      IF (gridpoint.is_in_area) THEN
       m:=1;n:=1;
       IF ((grid[(i-1)*gridsize+j].dist_to_border<=gridpoint.dist_to_border)AND(grid[(i+1)*gridsize+j].dist_to_border>gridpoint.dist_to_border))  THEN m:=0; END IF;
       IF ((grid[(i-1)*gridsize+j].dist_to_border>gridpoint.dist_to_border) AND(grid[(i+1)*gridsize+j].dist_to_border<=gridpoint.dist_to_border)) THEN m:=0; END IF; 
       IF ((grid[(i)*gridsize+j-1].dist_to_border<=gridpoint.dist_to_border)AND(grid[(i)*gridsize+j+1].dist_to_border>gridpoint.dist_to_border))  THEN n:=0; END IF;
       IF ((grid[(i)*gridsize+j-1].dist_to_border>gridpoint.dist_to_border) AND(grid[(i)*gridsize+j+1].dist_to_border<=gridpoint.dist_to_border)) THEN n:=0; END IF; 
       IF ((m=1)OR(n=1)) THEN
        gridpoint.is_near_axis:=true;
        if(imin>i) THEN imin:=i; END IF;
        if(jmin>j) THEN jmin:=j; END IF;
        if(imax<i) THEN imax:=i; END IF;
        if(jmax<j) THEN jmax:=j; END IF;
        grid[i*gridsize+j]:=gridpoint;
        IF (gridpoint.dist_to_border>p) THEN p:=gridpoint.dist_to_border;middleindex:=(i)*gridsize+j; END IF;
       END IF;
      END IF;
     END LOOP;
    END LOOP;
--
-- calc distance to the center point =(one of) the point(s) with the max dist to the border 
-- mark (one of the) node(s) with max. distance
-- the "distance" is not the lenght through the grid, but the straight line distance
--
    gridpoint=grid[middleindex];
    gridpoint.dist_to_middle:=1;
    grid[middleindex]:=gridpoint;
    startindex=middleindex;
    x:=gridpoint.x;y:=gridpoint.y;
    changes:=1;g:=0.0;m:=0.0;
    WHILE (changes>0) LOOP
      changes:=0;m:=m+1;
     FOR i IN imin..imax LOOP
      FOR j IN jmin..jmax LOOP
       gridpoint=grid[i*gridsize+j];
       IF ((gridpoint.is_near_axis) AND (gridpoint.dist_to_middle>0.0)) THEN
        FOR i1 IN (i-1)..(i+1) LOOP
         FOR j1 IN (j-1)..(j+1) LOOP
          IF ((i1!=i)OR(j1!=j)) THEN
           targetpoint:=grid[i1*gridsize+j1];
           IF (targetpoint.is_near_axis) THEN
            IF (j=j1) THEN f:=gridwidth; END IF;
            IF (i=i1) THEN f:=gridheight; END IF;
            IF ((i!=i1) AND (j!=j1)) THEN f:=griddiag ; END IF;
            IF ((targetpoint.dist_to_middle=0.0) OR (targetpoint.dist_to_middle>gridpoint.dist_to_middle+f)) THEN
             targetpoint.dist_to_middle:=gridpoint.dist_to_middle+f;
             targetpoint.way_to_middle:=i*gridsize+j;
             grid[i1*gridsize+j1]:=targetpoint;
             x1:=targetpoint.x; 
             y1:=targetpoint.y;
             p:=(x1-x)*(x1-x)+(y1-y)*(y1-y);
             IF (p>g) THEN g:=p;startindex:=i1*gridsize+j1; END IF;
             changes:=changes+1;
            END IF; 
           END IF;
          END IF;
         END LOOP;
        END LOOP; 
       END IF;
      END LOOP;
     END LOOP;
     FOR i IN REVERSE imax..imin LOOP
      FOR j IN REVERSE jmax..jmin LOOP
       gridpoint=grid[i*gridsize+j];
       IF ((gridpoint.is_near_axis) AND (gridpoint.dist_to_middle>0.0)) THEN
        FOR i1 IN (i-1)..(i+1) LOOP
         FOR j1 IN (j-1)..(j+1) LOOP
          IF ((i1!=i)OR(j1!=j)) THEN
           targetpoint:=grid[i1*gridsize+j1];
           IF (targetpoint.is_near_axis) THEN
            IF (j=j1) THEN f:=gridwidth; END IF;
            IF (i=i1) THEN f:=gridheight; END IF;
            IF ((i!=i1) AND (j!=j1)) THEN f:=griddiag ; END IF;
            IF ((targetpoint.dist_to_middle=0.0) OR (targetpoint.dist_to_middle>gridpoint.dist_to_middle+f)) THEN
             targetpoint.dist_to_middle:=gridpoint.dist_to_middle+f;
             targetpoint.way_to_middle:=i*gridsize+j;
             grid[i1*gridsize+j1]:=targetpoint;
             x1:=targetpoint.x; 
             y1:=targetpoint.y; 
             p:=(x1-x)*(x1-x)+(y1-y)*(y1-y);
             IF (p>g) THEN g:=p;startindex:=i1*gridsize+j1; END IF;
             changes:=changes+1;
            END IF; 
           END IF;
          END IF;
         END LOOP;
        END LOOP; 
       END IF;
      END LOOP;
     END LOOP;
    END LOOP;
--
-- test, if a straight line from middle to this point is possible
--

    tmplinestring:='LINESTRING(';
    gridpoint=grid[startindex];
    tmplinestring:=tmplinestring || gridpoint.x || ' ' || gridpoint.y || ',';
    x1:=gridpoint.x;y1:=gridpoint.y;
    gridpoint=grid[middleindex];
    tmplinestring:=tmplinestring || gridpoint.x || ' ' || gridpoint.y || ')';
    x:=gridpoint.x;y:=gridpoint.y;
    x1:=(x1-x)/(gridsize);
    y1:=(y1-y)/(gridsize);
    tmpway:=ST_GeomFromText(tmplinestring,3857);
    tmppoint:=ST_SetSRID(ST_MakePoint(x,y),3857);
    IF (ST_Within(tmpway,ST_SetSRID(myway,3857))) THEN 
     tmppoint:=ST_SetSRID(ST_MakePoint(x-x1,y-y1),3857);
     WHILE (ST_Within(tmppoint,ST_SetSRID(myway,3857))) LOOP
      x:=x-x1;y:=y-y1;
      tmpway:=ST_AddPoint(tmpway,tmppoint);
      tmppoint:=ST_SetSRID(ST_MakePoint(x-x1,y-y1),3857);
     END LOOP;
     x:=ST_X(ST_EndPoint(tmpway))-ST_X(ST_StartPoint(tmpway));
     y:=ST_Y(ST_EndPoint(tmpway))-ST_Y(ST_StartPoint(tmpway));
--
-- test, if this line covers at least 3/4 of the bounding box (in any direction)
--
     p:=ABS(x/(xmax-xmin));
     m:=ABS(y/(ymax-ymin));
     IF ((p>0.7) OR (m>0.7)) THEN
      rname:='d_' || rname;
      linestring:=ST_astext(tmpway);
      endindex:=middleindex;
      areashape:='diagonal';
      retway:=ST_GeomFromText(linestring,3857);
     END IF; 
    END IF;

    IF (retway IS NULL) THEN
--
-- calc distance to the start point, mark (one of the) node(s) with max. distance
-- the "distance" is not the lenght through the grid, but the straight line distance
--
     gridpoint=grid[startindex];
     gridpoint.dist_to_start:=1;
     grid[startindex]:=gridpoint;
     changes:=1;g:=0.0;m:=0.0;p:=0;x:=gridpoint.x;y:=gridpoint.y;
     WHILE (changes>0) LOOP
      changes:=0;m:=m+1;
      FOR i IN imin..imax LOOP
       FOR j IN jmin..jmax LOOP
        gridpoint=grid[i*gridsize+j];
        IF ((gridpoint.is_near_axis) AND (gridpoint.dist_to_start>0.0)) THEN
         FOR i1 IN (i-1)..(i+1) LOOP
          FOR j1 IN (j-1)..(j+1) LOOP
           IF ((i1!=i)OR(j1!=j)) THEN
            targetpoint:=grid[i1*gridsize+j1];
            IF (targetpoint.is_near_axis) THEN
             IF (j=j1) THEN f:=gridwidth; END IF;
             IF (i=i1) THEN f:=gridheight; END IF;
             IF ((i!=i1) AND (j!=j1)) THEN f:=griddiag ; END IF;
             IF ((targetpoint.dist_to_start=0.0) OR (targetpoint.dist_to_start>gridpoint.dist_to_start+f)) THEN
              targetpoint.dist_to_start:=gridpoint.dist_to_start+f;
              targetpoint.way_to_middle:=i*gridsize+j;
              grid[i1*gridsize+j1]:=targetpoint;
              x1:=targetpoint.x;
              y1:=targetpoint.y;
              p:=(x1-x)*(x1-x)+(y1-y)*(y1-y);
              if(p>g) THEN g:=p;endindex:=i1*gridsize+j1; END IF;
              changes:=changes+1;
             END IF; 
            END IF;
           END IF;
          END LOOP;
         END LOOP; 
        END IF;
       END LOOP;
      END LOOP;
      FOR i IN REVERSE imax..imin LOOP
       FOR j IN REVERSE jmax..jmin LOOP
        gridpoint=grid[i*gridsize+j];
        IF ((gridpoint.is_near_axis) AND (gridpoint.dist_to_start>0.0)) THEN
         FOR i1 IN (i-1)..(i+1) LOOP
          FOR j1 IN (j-1)..(j+1) LOOP
           IF ((i1!=i)OR(j1!=j)) THEN
            targetpoint:=grid[i1*gridsize+j1];
            IF (targetpoint.is_near_axis) THEN
             IF (j=j1) THEN f:=gridwidth; END IF;
             IF (i=i1) THEN f:=gridheight; END IF;
             IF ((i!=i1) AND (j!=j1)) THEN f:=griddiag ; END IF;
             IF ((targetpoint.dist_to_start=0.0) OR (targetpoint.dist_to_start>gridpoint.dist_to_start+f)) THEN
              targetpoint.dist_to_start:=gridpoint.dist_to_start+f;
              targetpoint.way_to_middle:=i*gridsize+j;
              grid[i1*gridsize+j1]:=targetpoint;
              x1:=targetpoint.x;
              y1:=targetpoint.y;
              p:=(x1-x)*(x1-x)+(y1-y)*(y1-y);
              if(p>g) THEN g:=p;endindex:=i1*gridsize+j1; END IF;
              changes:=changes+1;
             END IF; 
            END IF;
           END IF;
          END LOOP;
         END LOOP; 
        END IF;
       END LOOP;
      END LOOP;
     END LOOP;
--
-- Build linestring endindex to startindex
--
     i:=endindex;
     linestring:='LINESTRING(';
     WHILE (i!=startindex) LOOP
      gridpoint=grid[i];
      linestring:=linestring || gridpoint.x || ' ' || gridpoint.y || ',';
      i:=gridpoint.way_to_middle;
     END LOOP;
     gridpoint=grid[i];
     linestring:=linestring || gridpoint.x || ' ' || gridpoint.y || ')';
--
-- smooth line
-- 
     tmpway:=ST_GeomFromText(linestring,3857);
     IF (ST_NPoints(tmpway)>3) THEN
      tmppoint:=ST_PointN(tmpway,1);
      FOR i in 2..(ST_NPoints(tmpway)-1) LOOP
       x:=(ST_X(ST_PointN(tmpway,i-1))+ST_X(ST_PointN(tmpway,i))+ST_X(ST_PointN(tmpway,i+1)))/3;
       y:=(ST_Y(ST_PointN(tmpway,i-1))+ST_Y(ST_PointN(tmpway,i))+ST_Y(ST_PointN(tmpway,i+1)))/3;
       IF (i=2) THEN 
        smoothway:=ST_MakeLine(tmppoint,ST_SetSRID(ST_Point(x,y),3857));
       ELSE 
        smoothway:=ST_AddPoint(smoothway,ST_SetSRID(ST_Point(x,y),3857));
       END IF;
      END LOOP; 
      x:=ST_X(ST_PointN(tmpway,ST_NPoints(tmpway)));
      y:=ST_Y(ST_PointN(tmpway,ST_NPoints(tmpway)));
      smoothway:=ST_AddPoint(smoothway,ST_SetSRID(ST_Point(x,y),3857));
      smoothway:=OTM_CreateCurve(smoothway,90);
      retway:=smoothway;
      areashape:='curve';
      rname:='c_'||rname;
      done:=1;
     END IF;
    END IF;
   END IF;
  END IF;
--
-- DEBUG: insert nodes and lines
--
  FOR i IN 0..(gridsize-1) LOOP
   FOR j IN 0..(gridsize-1) LOOP
    IF (grid[(i)*gridsize+j].is_in_area) THEN
     gridpoint=grid[i*gridsize+j];
     x:=gridpoint.x;
     y:=gridpoint.y;     
     tmpway:=ST_SetSRID(ST_Point(x,y),3857);
     tmpbox:=ST_SetSRID(ST_Envelope(ST_MakeLine(ST_MakePoint(x-gridwidth/2,y-gridheight/2), ST_MakePoint(x+gridwidth/2,y+gridheight/2))),3857);
     t:=round(grid[(i)*gridsize+j].dist_to_border)::TEXT;
     IF (i*gridsize+j=middleindex) THEN 
      INSERT into planet_osm_point (osm_id,way,man_made,amenity,name) values (myosm_id,tmpway,'lakepoint','nearaxis',t);
      INSERT into planet_osm_polygon (osm_id,way,man_made,amenity,name) values (myosm_id,tmpbox,'lakepoint','middle',t);
     ELSE 
      IF (gridpoint.is_near_axis) THEN 
       INSERT into planet_osm_point (osm_id,way,man_made,amenity,name) values (myosm_id,tmpway,'lakepoint','nearaxis',t);
       INSERT into planet_osm_polygon (osm_id,way,man_made,amenity,name) values (myosm_id,tmpbox,'lakepoint','nearaxis',t);
      ELSE
       INSERT into planet_osm_point (osm_id,way,man_made,amenity,name) values (myosm_id,tmpway,'lakepoint','notaxis',t);
       INSERT into planet_osm_polygon (osm_id,way,man_made,amenity,name) values (myosm_id,tmpbox,'lakepoint','notaxis',t);      
      END IF;
     END IF;
    END IF;
   END LOOP;
  END LOOP;
  IF ( retway IS NOT NULL ) THEN
   x:=ST_XMax(ST_Envelope(retway))-ST_XMin(ST_Envelope(retway));
   y:=ST_YMax(ST_Envelope(retway))-ST_YMin(ST_Envelope(retway));
   x1:=ST_XMax(ST_Envelope(myway))-ST_XMin(ST_Envelope(myway));
   y1:=ST_YMax(ST_Envelope(myway))-ST_YMin(ST_Envelope(myway));
   IF ( GREATEST(x/x1,y/y1)<0.333 ) THEN
    RAISE NOTICE ' % is % with q=%',myosm_id,areashape,GREATEST(x/x1,y/y1);
   END IF;
   INSERT into planet_osm_line (osm_id,way,name,man_made,amenity) values (myosm_id,retway,rname,'lakeaxis','1');
  ELSE
   RAISE NOTICE 'unable to find way in area %',myosm_id;
  END IF;

  
  RETURN NULL;
 END;
$$ LANGUAGE plpgsql;


delete from planet_osm_point where man_made='lakepoint';
delete from planet_osm_line where man_made='lakeaxis';
delete from planet_osm_polygon where man_made='lakepoint';

-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where (osm_id=8102889)  and "natural"='water';
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where (osm_id=-32246)  and "natural"='water';
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where (osm_id=-32246 or osm_id=-3597692 or  osm_id=-1156846 or osm_id=-2066758 or osm_id=-154667 or osm_id=-2068020 or osm_id=-168897 or osm_id=-168892)  and "natural"='water';
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where (way && st_expand((select way from planet_osm_polygon where osm_id=-32246 and "natural"='water'),20000) OR (osm_id=-32246 or osm_id=-3597692 or  osm_id=-1156846 or osm_id=-2066758 or osm_id=-154667 or osm_id=-2068020 or osm_id=-168897 or osm_id=-168892)) and name is not null and "natural"='water' or ((osm_id=-32246 or osm_id=-3597692 or  osm_id=-1156846 or osm_id=-2066758 or osm_id=-154667 or osm_id=-2068020 or osm_id=-168897 or osm_id=-168892)  and "natural"='water');
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where ((osm_id=-32246 or osm_id=-3597692 or  osm_id=-1156846 or osm_id=-2066758 or osm_id=-154667 or osm_id=-2068020 or osm_id=-168897 or osm_id=-168892) or (way && st_expand((select way from planet_osm_polygon where "natural"='water' and osm_id=-32246),20000)) or (way && st_expand((select way from planet_osm_polygon where "natural"='water' and osm_id=24764378),20000))) and name is not null and "natural"='water';
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where "natural"='water' and osm_id=8116779;
-- SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where "natural"='water' ;
 SELECT osm_id,name,arealabel(osm_id,way,name) from planet_osm_polygon where "natural"='water'  and name  is not null;