parse_gml_curve()

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

static

Parse GML Curve (3.1.1)

Definition at line 994 of file lwgeom_in_gml.c.

References get_xlink_node(), getPoint_internal(), gml_lwpgerror(), gml_reproject_pa(), gmlGetProp(), is_gml_namespace(), is_xlink(), lwalloc(), lwfree(), lwline_construct(), lwrealloc(), POINTARRAY::npoints, parse_gml_data(), parse_gml_srs(), ptarray_construct(), ptarray_flip_coordinates(), ptarray_point_size(), struct_gmlSrs::reverse_axis, struct_gmlSrs::srid, and SRID_UNKNOWN.

Referenced by parse_gml().

{
        xmlNodePtr xa;
        int lss, last, i;
        bool found=false;
        gmlSrs srs;
        LWGEOM *geom=NULL;
        POINTARRAY *pa=NULL;
        POINTARRAY **ppa=NULL;
        uint32 npoints=0;
        xmlChar *interpolation=NULL;

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

        /* Looking for gml:segments */
        for (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, "segments"))
                {
                        found = true;
                        break;
                }
        }
        if (!found) gml_lwpgerror("invalid GML representation", 37);

        ppa = (POINTARRAY**) lwalloc(sizeof(POINTARRAY*));

        /* Processing each gml:LineStringSegment */
        for (xa = xa->children, lss=0; xa != NULL ; xa = xa->next)
        {
                if (xa->type != XML_ELEMENT_NODE) continue;
                if (!is_gml_namespace(xa, false)) continue;
                if (strcmp((char *) xa->name, "LineStringSegment")) continue;

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

                if (lss > 0) ppa = (POINTARRAY**) lwrealloc((POINTARRAY *) ppa,
                                           sizeof(POINTARRAY*) * (lss + 1));

                ppa[lss] = parse_gml_data(xa->children, hasz, root_srid);
                npoints += ppa[lss]->npoints;
                if (ppa[lss]->npoints < 2)
                        gml_lwpgerror("invalid GML representation", 39);
                lss++;
        }
        if (lss == 0) gml_lwpgerror("invalid GML representation", 40);

        /* Most common case, a single segment */
        if (lss == 1) pa = ppa[0];

        /*
         * "The curve segments are connected to one another, with the end point
         *  of each segment except the last being the start point of the next
         *  segment"  from  ISO 19107:2003 -> 6.3.16.1 (p43)
         *
         * So we must aggregate all the segments into a single one and avoid
         * to copy the redundants points
         */
        if (lss > 1)
        {
                pa = ptarray_construct(1, 0, npoints - (lss - 1));
                for (last = npoints = i = 0; i < lss ; i++)
                {
                        if (i + 1 == lss) last = 1;
                        /* Check if segments are not disjoints */
                        if (i > 0 && memcmp( getPoint_internal(pa, npoints),
                                             getPoint_internal(ppa[i], 0),
                                             *hasz ? sizeof(POINT3D) : sizeof(POINT2D)))
                                gml_lwpgerror("invalid GML representation", 41);

                        /* Aggregate stuff */
                        memcpy( getPoint_internal(pa, npoints),
                                getPoint_internal(ppa[i], 0),
                                ptarray_point_size(ppa[i]) * (ppa[i]->npoints + last));

                        npoints += ppa[i]->npoints - 1;
                        lwfree(ppa[i]);
                }
                lwfree(ppa);
        }

        parse_gml_srs(xnode, &srs);
        if (srs.reverse_axis) pa = ptarray_flip_coordinates(pa);
        if (srs.srid != *root_srid && *root_srid != SRID_UNKNOWN)
                gml_reproject_pa(pa, srs.srid, *root_srid);
        geom = (LWGEOM *) lwline_construct(*root_srid, NULL, pa);

        return geom;
}

Here is the call graph for this function:

Here is the caller graph for this function: