_lwt_EdgeMotionArea()

static GEOSGeometry* _lwt_EdgeMotionArea ( LWLINE *  geom, int  isclosed  )

static

Definition at line 3242 of file lwgeom_topo.c.

References getPoint4d_p(), LW_FAILURE, LW_TRUE, lwalloc(), lwerror(), lwfree(), LWGEOM2GEOS(), lwgeom_free(), lwgeom_geos_errmsg, lwgeom_geos_error(), lwgeom_make_valid(), lwnotice(), lwpoly_as_lwgeom(), lwpoly_construct(), lwpoly_free(), LWLINE::points, ptarray_append_point(), ptarray_clone_deep(), and ptarray_free().

Referenced by lwt_ChangeEdgeGeom().

{
  GEOSGeometry *gg;
  POINT4D p4d;
  POINTARRAY *pa;
  POINTARRAY **pas;
  LWPOLY *poly;
  LWGEOM *g;

  pas = lwalloc(sizeof(POINTARRAY*));

  initGEOS(lwnotice, lwgeom_geos_error);

  if ( isclosed )
  {
    pas[0] = ptarray_clone_deep( geom->points );
    poly = lwpoly_construct(0, 0, 1, pas);
    gg = LWGEOM2GEOS( lwpoly_as_lwgeom(poly), 0 );
    lwpoly_free(poly); /* should also delete the pointarrays */
  }
  else
  {
    pa = geom->points;
    getPoint4d_p(pa, 0, &p4d);
    pas[0] = ptarray_clone_deep( pa );
    /* don't bother dup check */
    if ( LW_FAILURE == ptarray_append_point(pas[0], &p4d, LW_TRUE) )
    {
      ptarray_free(pas[0]);
      lwfree(pas);
      lwerror("Could not append point to pointarray");
      return NULL;
    }
    poly = lwpoly_construct(0, NULL, 1, pas);
    /* make valid, in case the edge self-intersects on its first-last
     * vertex segment */
    g = lwgeom_make_valid(lwpoly_as_lwgeom(poly));
    lwpoly_free(poly); /* should also delete the pointarrays */
    if ( ! g )
    {
      lwerror("Could not make edge motion area valid");
      return NULL;
    }
    gg = LWGEOM2GEOS(g, 0);
    lwgeom_free(g);
  }
  if ( ! gg )
  {
    lwerror("Could not convert old edge area geometry to GEOS: %s",
            lwgeom_geos_errmsg);
    return NULL;
  }
  return gg;
}

Here is the call graph for this function:

Here is the caller graph for this function: