PostGIS  3.1.6dev-r@@SVN_REVISION@@

◆ lwt_GetFaceGeometry()

LWGEOM* lwt_GetFaceGeometry ( LWT_TOPOLOGY topo,
LWT_ELEMID  face 
)

Return the geometry of a face.

For ST_GetFaceGeometry

Parameters
topothe topology to operate on
faceidentifier of the face
Returns
a polygon geometry representing the face, ownership to caller, to be released with lwgeom_release, or NULL on error (liblwgeom error handler will be invoked with error message)

Definition at line 2824 of file lwgeom_topo.c.

2825 {
2826  uint64_t numfaceedges;
2827  LWT_ISO_EDGE *edges;
2828  LWT_ISO_FACE *face;
2829  LWPOLY *out;
2830  LWGEOM *outg;
2831  uint64_t i, edgeid;
2832  int fields;
2833 
2834  if (faceid == 0)
2835  {
2836  lwerror("SQL/MM Spatial exception - universal face has no geometry");
2837  return NULL;
2838  }
2839 
2840  /* Construct the face geometry */
2841  numfaceedges = 1;
2842  fields = LWT_COL_EDGE_GEOM |
2846  ;
2847  edges = lwt_be_getEdgeByFace( topo, &faceid, &numfaceedges, fields, NULL );
2848  if (numfaceedges == UINT64_MAX)
2849  {
2850  lwerror("Backend error: %s", lwt_be_lastErrorMessage(topo->be_iface));
2851  return NULL;
2852  }
2853 
2854  if ( numfaceedges == 0 )
2855  {
2856  i = 1;
2857  face = lwt_be_getFaceById(topo, &faceid, &i, LWT_COL_FACE_FACE_ID);
2858  if (i == UINT64_MAX)
2859  {
2860  lwerror("Backend error: %s", lwt_be_lastErrorMessage(topo->be_iface));
2861  return NULL;
2862  }
2863  if ( i == 0 ) {
2864  lwerror("SQL/MM Spatial exception - non-existent face.");
2865  return NULL;
2866  }
2867  lwfree( face );
2868  if ( i > 1 ) {
2869  lwerror("Corrupted topology: multiple face records have face_id=%"
2870  LWTFMT_ELEMID, faceid);
2871  return NULL;
2872  }
2873  /* Face has no boundary edges, we'll return EMPTY, see
2874  * https://trac.osgeo.org/postgis/ticket/3221 */
2875  lwnotice("Corrupted topology: face %"
2876  LWTFMT_ELEMID " has no associated edges.", faceid);
2877  out = lwpoly_construct_empty(topo->srid, topo->hasZ, 0);
2878  return lwpoly_as_lwgeom(out);
2879  }
2880  edgeid = edges[0].edge_id;
2881 
2882  outg = _lwt_FaceByEdges( topo, edges, numfaceedges );
2883  _lwt_release_edges(edges, numfaceedges);
2884 
2885  if ( ! outg )
2886  {
2887  /* Face did have edges but no polygon could be constructed
2888  * with that material, sounds like a corrupted topology..
2889  *
2890  * We'll return EMPTY, see
2891  * https://trac.osgeo.org/postgis/ticket/3221 */
2892  lwnotice("Corrupted topology: face %"
2893  LWTFMT_ELEMID " could not be constructed only from edges "
2894  "knowing about it (like edge %" LWTFMT_ELEMID ").",
2895  faceid, edgeid);
2896  out = lwpoly_construct_empty(topo->srid, topo->hasZ, 0);
2897  return lwpoly_as_lwgeom(out);
2898  }
2899 
2900  return outg;
2901 }
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
Definition: lwgeom.c:312
void lwfree(void *mem)
Definition: lwutil.c:242
LWPOLY * lwpoly_construct_empty(int32_t srid, char hasz, char hasm)
Definition: lwpoly.c:161
#define LWT_COL_EDGE_FACE_RIGHT
#define LWT_COL_FACE_FACE_ID
Face fields.
#define LWT_COL_EDGE_FACE_LEFT
#define LWT_COL_EDGE_EDGE_ID
Edge fields.
#define LWT_COL_EDGE_GEOM
void lwerror(const char *fmt,...)
Write a notice out to the error handler.
Definition: lwutil.c:190
void lwnotice(const char *fmt,...)
Write a notice out to the notice handler.
Definition: lwutil.c:177
const char * lwt_be_lastErrorMessage(const LWT_BE_IFACE *be)
Definition: lwgeom_topo.c:119
static LWT_ISO_EDGE * lwt_be_getEdgeByFace(LWT_TOPOLOGY *topo, const LWT_ELEMID *ids, uint64_t *numelems, int fields, const GBOX *box)
Definition: lwgeom_topo.c:238
static LWGEOM * _lwt_FaceByEdges(LWT_TOPOLOGY *topo, LWT_ISO_EDGE *edges, int numfaceedges)
Definition: lwgeom_topo.c:2777
static LWT_ISO_FACE * lwt_be_getFaceById(LWT_TOPOLOGY *topo, const LWT_ELEMID *ids, uint64_t *numelems, int fields)
Definition: lwgeom_topo.c:226
#define LWTFMT_ELEMID
Definition: lwgeom_topo.c:43
static void _lwt_release_edges(LWT_ISO_EDGE *edges, int num_edges)
Definition: lwgeom_topo.c:451
LWT_ELEMID edge_id
const LWT_BE_IFACE * be_iface

References _lwt_FaceByEdges(), _lwt_release_edges(), LWT_TOPOLOGY_T::be_iface, LWT_ISO_EDGE::edge_id, LWT_TOPOLOGY_T::hasZ, lwerror(), lwfree(), lwnotice(), lwpoly_as_lwgeom(), lwpoly_construct_empty(), lwt_be_getEdgeByFace(), lwt_be_getFaceById(), lwt_be_lastErrorMessage(), LWT_COL_EDGE_EDGE_ID, LWT_COL_EDGE_FACE_LEFT, LWT_COL_EDGE_FACE_RIGHT, LWT_COL_EDGE_GEOM, LWT_COL_FACE_FACE_ID, LWTFMT_ELEMID, and LWT_TOPOLOGY_T::srid.

Referenced by lwt_AddPolygon(), and lwt_ChangeEdgeGeom().

Here is the call graph for this function:
Here is the caller graph for this function: