PostGIS 3.7.0dev-r@@SVN_REVISION@@
Loading...
Searching...
No Matches

◆ parse_gml_polygon()

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

Parse GML Polygon (2.1.2, 3.1.1)

Definition at line 1235 of file lwgeom_in_gml.c.

1236{
1237 gmlSrs srs;
1238 int i, ring;
1239 LWGEOM *geom;
1240 xmlNodePtr xa, xb;
1241 POINTARRAY **ppa = NULL;
1242
1243 if (is_xlink(xnode)) xnode = get_xlink_node(xnode);
1244 if (xnode == NULL)
1245 gml_lwpgerror("invalid GML representation", 30);
1246
1247 if (xnode->children == NULL)
1248 return lwpoly_as_lwgeom(lwpoly_construct_empty(*root_srid, 0, 0));
1249
1250 parse_gml_srs(xnode, &srs);
1251
1252 for (xa = xnode->children ; xa != NULL ; xa = xa->next)
1253 {
1254 /* Polygon/outerBoundaryIs -> GML 2.1.2 */
1255 /* Polygon/exterior -> GML 3.1.1 */
1256 if (xa->type != XML_ELEMENT_NODE) continue;
1257 if (!is_gml_namespace(xa, false)) continue;
1258 if (!(is_gml_element(xa, "outerBoundaryIs") ||
1259 is_gml_element(xa, "exterior")))
1260 continue;
1261
1262 for (xb = xa->children ; xb != NULL ; xb = xb->next)
1263 {
1264 if (xb->type != XML_ELEMENT_NODE) continue;
1265 if (!is_gml_namespace(xb, false)) continue;
1266 if (!is_gml_element(xb, "LinearRing")) continue;
1267
1268 ppa = (POINTARRAY**) lwalloc(sizeof(POINTARRAY*));
1269 ppa[0] = parse_gml_data(xb->children, hasz, root_srid);
1270
1271 if (ppa[0]->npoints < 4
1272 || (!*hasz && !ptarray_is_closed_2d(ppa[0]))
1273 || (*hasz && !ptarray_is_closed_3d(ppa[0])))
1274 gml_lwpgerror("invalid GML representation", 43);
1275
1276 if (srs.reverse_axis) ppa[0] = ptarray_flip_coordinates(ppa[0]);
1277 }
1278 }
1279
1280 /* Found an <exterior> or <outerBoundaryIs> but no rings?!? We're outa here! */
1281 if ( ! ppa )
1282 gml_lwpgerror("invalid GML representation", 43);
1283
1284 for (ring=1, xa = xnode->children ; xa != NULL ; xa = xa->next)
1285 {
1286 /* Polygon/innerBoundaryIs -> GML 2.1.2 */
1287 /* Polygon/interior -> GML 3.1.1 */
1288 if (xa->type != XML_ELEMENT_NODE) continue;
1289 if (!is_gml_namespace(xa, false)) continue;
1290 if (!(is_gml_element(xa, "innerBoundaryIs") ||
1291 is_gml_element(xa, "interior")))
1292 continue;
1293
1294 for (xb = xa->children ; xb != NULL ; xb = xb->next)
1295 {
1296 if (xb->type != XML_ELEMENT_NODE) continue;
1297 if (!is_gml_namespace(xb, false)) continue;
1298 if (!is_gml_element(xb, "LinearRing")) continue;
1299
1300 ppa = (POINTARRAY**) lwrealloc((POINTARRAY *) ppa,
1301 sizeof(POINTARRAY*) * (ring + 1));
1302 ppa[ring] = parse_gml_data(xb->children, hasz, root_srid);
1303
1304 if (ppa[ring]->npoints < 4
1305 || (!*hasz && !ptarray_is_closed_2d(ppa[ring]))
1306 || (*hasz && !ptarray_is_closed_3d(ppa[ring])))
1307 gml_lwpgerror("invalid GML representation", 43);
1308
1309 if (srs.reverse_axis) ppa[ring] = ptarray_flip_coordinates(ppa[ring]);
1310 ring++;
1311 }
1312 }
1313
1314 /* Exterior Ring is mandatory */
1315 if (ppa == NULL || ppa[0] == NULL) gml_lwpgerror("invalid GML representation", 44);
1316
1317 if (srs.srid != *root_srid && *root_srid != SRID_UNKNOWN)
1318 {
1319 for (i=0 ; i < ring ; i++)
1320 gml_reproject_pa(ppa[i], srs.srid, *root_srid);
1321 }
1322 geom = (LWGEOM *) lwpoly_construct(*root_srid, NULL, ring, ppa);
1323
1324 return geom;
1325}
void * lwrealloc(void *mem, size_t size)
Definition lwutil.c:242
int ptarray_is_closed_3d(const POINTARRAY *pa)
Definition ptarray.c:723
void * lwalloc(size_t size)
Definition lwutil.c:227
POINTARRAY * ptarray_flip_coordinates(POINTARRAY *pa)
Reverse X and Y axis on a given POINTARRAY.
Definition ptarray.c:368
LWPOLY * lwpoly_construct(int32_t srid, GBOX *bbox, uint32_t nrings, POINTARRAY **points)
Definition lwpoly.c:43
int ptarray_is_closed_2d(const POINTARRAY *pa)
Definition ptarray.c:710
LWPOLY * lwpoly_construct_empty(int32_t srid, char hasz, char hasm)
Definition lwpoly.c:161
#define SRID_UNKNOWN
Unknown SRID value.
Definition liblwgeom.h:215
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
Definition lwgeom.c:357
static xmlNodePtr get_xlink_node(xmlNodePtr xnode)
Return a xmlNodePtr on a node referenced by a XLink or NULL otherwise.
static bool is_xlink(xmlNodePtr node)
Return true if current node contains a simple XLink Return false otherwise.
static POINTARRAY * gml_reproject_pa(POINTARRAY *pa, int32_t epsg_in, int32_t epsg_out)
Use Proj to reproject a given POINTARRAY.
static bool is_gml_namespace(xmlNodePtr xnode, bool is_strict)
Return false if current element namespace is not a GML one Return true otherwise.
static void parse_gml_srs(xmlNodePtr xnode, gmlSrs *srs)
Parse gml srsName attribute.
static void gml_lwpgerror(char *msg, __attribute__((__unused__)) int error_code)
static POINTARRAY * parse_gml_data(xmlNodePtr xnode, bool *hasz, int *root_srid)
Parse data coordinates.
static bool is_gml_element(xmlNodePtr xn, const char *gml_name)

References get_xlink_node(), gml_lwpgerror(), gml_reproject_pa(), is_gml_element(), is_gml_namespace(), is_xlink(), lwalloc(), lwpoly_as_lwgeom(), lwpoly_construct(), lwpoly_construct_empty(), lwrealloc(), 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().

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