PostGIS  2.5.0dev-r@@SVN_REVISION@@

◆ parse_gml_triangle()

static LWGEOM* parse_gml_triangle ( xmlNodePtr  xnode,
bool *  hasz,
int *  root_srid 
)
static

Parse GML Triangle (3.1.1)

Definition at line 1221 of file lwgeom_in_gml.c.

References get_xlink_node(), gml_lwpgerror(), gml_reproject_pa(), gmlGetProp(), is_gml_namespace(), is_xlink(), lwalloc(), lwtriangle_as_lwgeom(), lwtriangle_construct(), lwtriangle_construct_empty(), POINTARRAY::npoints, parse_gml_data(), parse_gml_srs(), ptarray_flip_coordinates(), ptarray_is_closed_2d(), ptarray_is_closed_3d(), struct_gmlSrs::reverse_axis, struct_gmlSrs::srid, and SRID_UNKNOWN.

Referenced by parse_gml(), and parse_gml_tin().

1222 {
1223  gmlSrs srs;
1224  LWGEOM *geom;
1225  xmlNodePtr xa, xb;
1226  POINTARRAY *pa = NULL;
1227  xmlChar *interpolation=NULL;
1228 
1229  if (is_xlink(xnode)) xnode = get_xlink_node(xnode);
1230 
1231  if (xnode->children == NULL)
1232  return lwtriangle_as_lwgeom(lwtriangle_construct_empty(*root_srid, 0, 0));
1233 
1234  /* GML SF is restricted to planar interpolation
1235  NOTA: I know Triangle is not part of SF, but
1236  we have to be consistent with other surfaces */
1237  interpolation = gmlGetProp(xnode, (xmlChar *) "interpolation");
1238  if (interpolation != NULL)
1239  {
1240  if (strcmp((char *) interpolation, "planar"))
1241  gml_lwpgerror("invalid GML representation", 45);
1242  xmlFree(interpolation);
1243  }
1244 
1245  parse_gml_srs(xnode, &srs);
1246 
1247  for (xa = xnode->children ; xa != NULL ; xa = xa->next)
1248  {
1249  /* Triangle/exterior */
1250  if (xa->type != XML_ELEMENT_NODE) continue;
1251  if (!is_gml_namespace(xa, false)) continue;
1252  if (strcmp((char *) xa->name, "exterior")) continue;
1253 
1254  for (xb = xa->children ; xb != NULL ; xb = xb->next)
1255  {
1256  /* Triangle/exterior/LinearRing */
1257  if (xb->type != XML_ELEMENT_NODE) continue;
1258  if (!is_gml_namespace(xb, false)) continue;
1259  if (strcmp((char *) xb->name, "LinearRing")) continue;
1260 
1261  pa = (POINTARRAY*) lwalloc(sizeof(POINTARRAY));
1262  pa = parse_gml_data(xb->children, hasz, root_srid);
1263 
1264  if (pa->npoints != 4
1265  || (!*hasz && !ptarray_is_closed_2d(pa))
1266  || (*hasz && !ptarray_is_closed_3d(pa)))
1267  gml_lwpgerror("invalid GML representation", 46);
1268 
1269  if (srs.reverse_axis) pa = ptarray_flip_coordinates(pa);
1270  }
1271  }
1272 
1273  /* Exterior Ring is mandatory */
1274  if (pa == NULL) gml_lwpgerror("invalid GML representation", 47);
1275 
1276  if (srs.srid != *root_srid && *root_srid != SRID_UNKNOWN)
1277  gml_reproject_pa(pa, srs.srid, *root_srid);
1278 
1279  geom = (LWGEOM *) lwtriangle_construct(*root_srid, NULL, pa);
1280 
1281  return geom;
1282 }
static void gml_lwpgerror(char *msg, __attribute__((__unused__)) int error_code)
Definition: lwgeom_in_gml.c:82
LWTRIANGLE * lwtriangle_construct_empty(int srid, char hasz, char hasm)
Definition: lwtriangle.c:58
static void parse_gml_srs(xmlNodePtr xnode, gmlSrs *srs)
Parse gml srsName attribute.
int ptarray_is_closed_3d(const POINTARRAY *pa)
Definition: ptarray.c:701
static xmlChar * gmlGetProp(xmlNodePtr xnode, xmlChar *prop)
Retrieve a GML property from a node or NULL otherwise Respect namespaces if presents in the node elem...
int ptarray_is_closed_2d(const POINTARRAY *pa)
Definition: ptarray.c:688
static POINTARRAY * gml_reproject_pa(POINTARRAY *pa, int srid_in, int srid_out)
Use Proj4 to reproject a given POINTARRAY.
POINTARRAY * ptarray_flip_coordinates(POINTARRAY *pa)
Reverse X and Y axis on a given POINTARRAY.
Definition: ptarray.c:368
static xmlNodePtr get_xlink_node(xmlNodePtr xnode)
Return a xmlNodePtr on a node referenced by a XLink or NULL otherwise.
#define SRID_UNKNOWN
Unknown SRID value.
Definition: liblwgeom.h:187
LWGEOM * lwtriangle_as_lwgeom(const LWTRIANGLE *obj)
Definition: lwgeom.c:325
static bool is_xlink(xmlNodePtr node)
Return true if current node contains a simple XLink Return false otherwise.
LWTRIANGLE * lwtriangle_construct(int srid, GBOX *bbox, POINTARRAY *points)
Definition: lwtriangle.c:40
static POINTARRAY * parse_gml_data(xmlNodePtr xnode, bool *hasz, int *root_srid)
Parse data coordinates.
static bool is_gml_namespace(xmlNodePtr xnode, bool is_strict)
Return false if current element namespace is not a GML one Return true otherwise. ...
void * lwalloc(size_t size)
Definition: lwutil.c:229
uint32_t npoints
Definition: liblwgeom.h:373
Here is the call graph for this function:
Here is the caller graph for this function: