parse_marc21()

static LWGEOM * parse_marc21 ( xmlNodePtr  xnode )

static

MARC21/XML documents must have <record> as top level element. https://www.loc.gov/standards/marcxml/xml/spy/spy.html

If the coordinates are given in terms of a center point rather than outside limits, the longitude and latitude which form the central axis are recorded twice (in subfields $d and $e and in $f and $g, respectively).

MARC21/XML Document has a <record> but no <datafield:034>

Definition at line 364 of file lwgeom_in_marc21.c.

                               {
 
        int ngeoms;
        int i;
        xmlNodePtr datafield;
        xmlNodePtr subfield;
        LWGEOM *result;
        LWGEOM **lwgeoms = (LWGEOM**) lwalloc(sizeof(LWGEOM*));
        uint8_t geometry_type;
        uint8_t result_type;
        char *code;
        char *literal;
 
        POSTGIS_DEBUGF(2, "parse_marc21 called: root '<%s>'", xnode->name);
 
        if (xmlStrcmp(xnode->name, (xmlChar*) "record")) lwpgerror("invalid MARC21/XML document. Root element <record> expected but <%s> found.",xnode->name);
 
        result_type = 0;
        ngeoms = 0;
 
        for (datafield = xnode->children; datafield != NULL; datafield = datafield->next) {
 
                char *lw = NULL;
                char *le = NULL;
                char *ln = NULL;
                char *ls = NULL;
 
                if (datafield->type != XML_ELEMENT_NODE) continue;
 
                if (xmlStrcmp(datafield->name, (xmlChar*) "datafield") != 0     || xmlStrcmp(xmlGetProp(datafield, (xmlChar*) "tag"),(xmlChar*) "034") != 0) continue;
 
                POSTGIS_DEBUG(3, "  datafield found");
 
                for (subfield = datafield->children; subfield != NULL; subfield = subfield->next) {
 
                        if (subfield->type != XML_ELEMENT_NODE) continue;
                        if (xmlStrcmp(subfield->name, (xmlChar*) "subfield") != 0) continue;
 
                        code = (char*) xmlGetProp(subfield, (xmlChar*) "code");
 
                        if ((strcmp(code, "d") != 0 && strcmp(code, "e") != 0 && strcmp(code, "f") != 0 && strcmp(code, "g")) != 0)     continue;
 
                        literal = (char*) xmlNodeGetContent(subfield);
 
                        POSTGIS_DEBUGF(3, "    subfield code '%s': %s", code, literal);
 
                        if (is_literal_valid(literal) == LW_TRUE) {
 
                                if (strcmp(code, "d") == 0) lw = literal;
                                else if (strcmp(code, "e") == 0) le = literal;
                                else if (strcmp(code, "f") == 0) ln = literal;
                                else if (strcmp(code, "g") == 0) ls = literal;
 
                        } else {
 
                                lwpgerror("parse error - invalid literal at 034$%s: \"%s\"", code, literal);
 
                        }
 
                }
 
                xmlFreeNode(subfield);
 
                if (lw && le && ln && ls) {
 
                        double w = parse_geo_literal(lw);
                        double e = parse_geo_literal(le);
                        double n = parse_geo_literal(ln);
                        double s = parse_geo_literal(ls);
                        geometry_type = 0;
 
                        if (ngeoms > 0) lwgeoms = (LWGEOM**) lwrealloc(lwgeoms, sizeof(LWGEOM*) * (ngeoms + 1));
 
                        if (fabs(w - e) < 0.0000001f && fabs(n - s) < 0.0000001f) {
 
                                lwgeoms[ngeoms] = (LWGEOM*) lwpoint_make2d(SRID_UNKNOWN, w, s);
                                geometry_type = MULTIPOINTTYPE;
 
                        } else {
 
                                lwgeoms[ngeoms] = (LWGEOM*) lwpoly_construct_envelope(SRID_UNKNOWN, w, n, e, s);
                                geometry_type = MULTIPOLYGONTYPE;
 
                        }
 
                        if (ngeoms && result_type != geometry_type) {
                                result_type = COLLECTIONTYPE;
                        } else {
                                result_type = geometry_type;
                        }
 
                        ngeoms++;
 
                } else {
 
                        if (lw || le || ln || ls) {
 
                                lwpgerror("parse error - the Coded Cartographic Mathematical Data (datafield:034) in the given MARC21/XML is incomplete. Coordinates for subfields \"$d\",\"$e\",\"$f\" and \"$g\" are expected.");
                        }
 
                }
 
        }
 
        POSTGIS_DEBUG(5, "  xmlFreeNode(datafield)");
        xmlFreeNode(datafield);
 
        if (ngeoms == 1) {
 
                POSTGIS_DEBUGF(2, "=> parse_marc21 returns single geometry: %s",lwtype_name(lwgeom_get_type(lwgeoms[0])));
                lwgeom_force_clockwise(lwgeoms[0]);
                return lwgeoms[0];
 
        } else if (ngeoms > 1) {
 
                result = (LWGEOM*) lwcollection_construct_empty(result_type,SRID_UNKNOWN, 0, 0);
 
                for (i = 0; i < ngeoms; i++) {
 
                        POSTGIS_DEBUGF(3, "  adding geometry to result set: %s",lwtype_name(lwgeom_get_type(lwgeoms[i])));
                        lwgeom_force_clockwise(lwgeoms[i]);
                        result = (LWGEOM*) lwcollection_add_lwgeom((LWCOLLECTION*) result,lwgeoms[i]);
 
                }
 
                POSTGIS_DEBUGF(2, "=> parse_marc21 returns a collection: %s", lwtype_name(lwgeom_get_type(result)));
                return result;
 
        }
 
        POSTGIS_DEBUG(2, "=> parse_marc21 returns NULL");
        return NULL;
 
}

References COLLECTIONTYPE, is_literal_valid(), LW_TRUE, lwalloc(), lwcollection_add_lwgeom(), lwcollection_construct_empty(), lwgeom_force_clockwise(), lwgeom_get_type(), lwpoint_make2d(), lwpoly_construct_envelope(), lwrealloc(), lwtype_name(), MULTIPOINTTYPE, MULTIPOLYGONTYPE, parse_geo_literal(), result, s, SRID_UNKNOWN, and w.

Referenced by ST_GeomFromMARC21().

Here is the call graph for this function:

Here is the caller graph for this function: