parse_gml_patch()

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

static

Parse GML PolygonPatch (3.1.1)

Definition at line 1288 of file lwgeom_in_gml.c.

References gml_lwpgerror(), gml_reproject_pa(), gmlGetProp(), is_gml_namespace(), lwalloc(), lwpoly_construct(), 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_psurface(), and parse_gml_surface().

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

        /* PolygonPatch */
        if (strcmp((char *) xnode->name, "PolygonPatch"))
                gml_lwpgerror("invalid GML representation", 48);

        /* GML SF is resticted to planar interpolation  */
        interpolation = gmlGetProp(xnode, (xmlChar *) "interpolation");
        if (interpolation != NULL)
        {
                if (strcmp((char *) interpolation, "planar"))
                        gml_lwpgerror("invalid GML representation", 48);
                xmlFree(interpolation);
        }

        parse_gml_srs(xnode, &srs);

        /* PolygonPatch/exterior */
        for (xa = xnode->children ; xa != NULL ; xa = xa->next)
        {
                if (!is_gml_namespace(xa, false)) continue;
                if (strcmp((char *) xa->name, "exterior")) continue;

                /* PolygonPatch/exterior/LinearRing */
                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", 48);

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

        /* Interior but no Exterior ! */
        if ( ! ppa )
                gml_lwpgerror("invalid GML representation", 48);

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

                /* PolygonPatch/interior/LinearRing */
                for (xb = xa->children ; xb != NULL ; xb = xb->next)
                {
                        if (xb->type != XML_ELEMENT_NODE) 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", 49);

                        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", 50);

        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: