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 1129 of file lwgeom_in_gml.c.

References get_xlink_node(), gml_lwpgerror(), gml_reproject_pa(), 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().

{
        gmlSrs srs;
        int i, ring;
        LWGEOM *geom;
        xmlNodePtr xa, xb;
        POINTARRAY **ppa = NULL;

        if (is_xlink(xnode)) xnode = get_xlink_node(xnode);

        if (xnode->children == NULL)
                return lwpoly_as_lwgeom(lwpoly_construct_empty(*root_srid, 0, 0));

        parse_gml_srs(xnode, &srs);

        for (xa = xnode->children ; xa != NULL ; xa = xa->next)
        {
                /* Polygon/outerBoundaryIs -> GML 2.1.2 */
                /* Polygon/exterior        -> GML 3.1.1 */
                if (xa->type != XML_ELEMENT_NODE) continue;
                if (!is_gml_namespace(xa, false)) continue;
                if  (strcmp((char *) xa->name, "outerBoundaryIs") &&
                        strcmp((char *) xa->name, "exterior")) continue;

                for (xb = xa->children ; xb != NULL ; xb = xb->next)
                {
                        if (xb->type != XML_ELEMENT_NODE) continue;
                        if (!is_gml_namespace(xb, false)) continue;
                        if (strcmp((char *) xb->name, "LinearRing")) continue;

                        ppa = (POINTARRAY**) lwalloc(sizeof(POINTARRAY*));
                        ppa[0] = parse_gml_data(xb->children, hasz, root_srid);

                        if (ppa[0]->npoints < 4
                                || (!*hasz && !ptarray_is_closed_2d(ppa[0]))
                                ||  (*hasz && !ptarray_is_closed_3d(ppa[0])))
                                gml_lwpgerror("invalid GML representation", 43);

                        if (srs.reverse_axis) ppa[0] = ptarray_flip_coordinates(ppa[0]);
                }
        }

        /* Found an <exterior> or <outerBoundaryIs> but no rings?!? We're outa here! */
        if ( ! ppa )
                gml_lwpgerror("invalid GML representation", 43);

        for (ring=1, xa = xnode->children ; xa != NULL ; xa = xa->next)
        {
                /* Polygon/innerBoundaryIs -> GML 2.1.2 */
                /* Polygon/interior        -> GML 3.1.1 */
                if (xa->type != XML_ELEMENT_NODE) continue;
                if (!is_gml_namespace(xa, false)) continue;
                if  (strcmp((char *) xa->name, "innerBoundaryIs") &&
                        strcmp((char *) xa->name, "interior")) continue;

                for (xb = xa->children ; xb != NULL ; xb = xb->next)
                {
                        if (xb->type != XML_ELEMENT_NODE) continue;
                        if (!is_gml_namespace(xb, false)) continue;
                        if (strcmp((char *) xb->name, "LinearRing")) continue;

                        ppa = (POINTARRAY**) lwrealloc((POINTARRAY *) ppa,
                                                       sizeof(POINTARRAY*) * (ring + 1));
                        ppa[ring] = parse_gml_data(xb->children, hasz, root_srid);

                        if (ppa[ring]->npoints < 4
                                || (!*hasz && !ptarray_is_closed_2d(ppa[ring]))
                                ||  (*hasz && !ptarray_is_closed_3d(ppa[ring])))
                                gml_lwpgerror("invalid GML representation", 43);

                        if (srs.reverse_axis) ppa[ring] = ptarray_flip_coordinates(ppa[ring]);
                        ring++;
                }
        }

        /* Exterior Ring is mandatory */
        if (ppa == NULL || ppa[0] == NULL) gml_lwpgerror("invalid GML representation", 44);

        if (srs.srid != *root_srid && *root_srid != SRID_UNKNOWN)
        {
                for (i=0 ; i < ring ; i++)
                        gml_reproject_pa(ppa[i], srs.srid, *root_srid);
        }
        geom = (LWGEOM *) lwpoly_construct(*root_srid, NULL, ring, ppa);

        return geom;
}

Here is the call graph for this function:

Here is the caller graph for this function: