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).
384 {
385
386 int ngeoms;
387 int i;
388 xmlNodePtr datafield;
389 xmlNodePtr subfield;
392 uint8_t geometry_type;
393 uint8_t result_type;
394 char *code;
395 char *literal;
396
397 POSTGIS_DEBUGF(2, "parse_marc21 called: root '<%s>'", xnode->name);
398
405 lwpgerror("invalid MARC21/XML document. Root element <record> expected but <%s> found.",xnode->name);
406
407 result_type = 0;
408 ngeoms = 0;
409
410 for (datafield = xnode->children; datafield != NULL; datafield = datafield->next) {
411
412 char *lw = NULL;
413 char *le = NULL;
414 char *ln = NULL;
415 char *ls = NULL;
416
417 if (datafield->type != XML_ELEMENT_NODE) continue;
418
419 if (!
is_xml_element(datafield,
"datafield") || xmlStrcmp(xmlGetProp(datafield, (xmlChar*)
"tag"),(xmlChar*)
"034") != 0)
continue;
420
421 POSTGIS_DEBUG(3, " datafield found");
422
423 for (subfield = datafield->children; subfield != NULL; subfield = subfield->next) {
424
425 if (subfield->type != XML_ELEMENT_NODE) continue;
427 continue;
428
429 code = (char*) xmlGetProp(subfield, (xmlChar*) "code");
430
431 if ((strcmp(code, "d") != 0 &&
432 strcmp(code, "e") != 0 &&
433 strcmp(code, "f") != 0 &&
434 strcmp(code, "g")) != 0)
435 continue;
436
437 literal = (char*) xmlNodeGetContent(subfield);
438
439 POSTGIS_DEBUGF(3, " subfield code '%s': %s", code, literal);
440
442
443 if (strcmp(code, "d") == 0) lw = literal;
444 else if (strcmp(code, "e") == 0) le = literal;
445 else if (strcmp(code, "f") == 0) ln = literal;
446 else if (strcmp(code, "g") == 0) ls = literal;
447
448 } else {
449
450 lwpgerror("parse error - invalid literal at 034$%s: \"%s\"", code, literal);
451
452 }
453
454 }
455
456 xmlFreeNode(subfield);
457
458 if (lw && le && ln && ls) {
459
464 geometry_type = 0;
465
466 if (ngeoms > 0) lwgeoms = (
LWGEOM**)
468
469 if (fabs(
w - e) < 0.0000001f && fabs(n -
s) < 0.0000001f) {
470
479
480 } else {
481
484
485 }
486
487 if (ngeoms && result_type != geometry_type) {
489 } else {
490 result_type = geometry_type;
491 }
492
493 ngeoms++;
494
495 } else {
496
497 if (lw || le || ln || ls) {
498
499 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.");
500 }
501
502 }
503
504 }
505
506 POSTGIS_DEBUG(5, " xmlFreeNode(datafield)");
507 xmlFreeNode(datafield);
508
509 if (ngeoms == 1) {
510
513 return lwgeoms[0];
514
515 } else if (ngeoms > 1) {
516
518
519 for (i = 0; i < ngeoms; i++) {
520
524
525 }
526
529
530 }
531
535 POSTGIS_DEBUG(2, "=> parse_marc21 returns NULL");
536 return NULL;
537
538}
char result[OUT_DOUBLE_BUFFER_SIZE]
const char * lwtype_name(uint8_t type)
Return the type name string associated with a type number (e.g.
void * lwrealloc(void *mem, size_t size)
LWPOLY * lwpoly_construct_envelope(int32_t srid, double x1, double y1, double x2, double y2)
void * lwalloc(size_t size)
LWPOINT * lwpoint_make2d(int32_t srid, double x, double y)
void lwgeom_force_clockwise(LWGEOM *lwgeom)
LWCOLLECTION * lwcollection_construct_empty(uint8_t type, int32_t srid, char hasz, char hasm)
#define LW_TRUE
Return types for functions with status returns.
LWCOLLECTION * lwcollection_add_lwgeom(LWCOLLECTION *col, const LWGEOM *geom)
Appends geom to the collection managed by col.
#define SRID_UNKNOWN
Unknown SRID value.
static bool is_xml_element(xmlNodePtr xn, const char *xml_name)
static double parse_geo_literal(char *literal)
static int is_literal_valid(const char *literal)
static uint32_t lwgeom_get_type(const LWGEOM *geom)
Return LWTYPE number.