parse_geo_literal()

static double parse_geo_literal ( char *  literal )

static

Coordinate formats supported (from the MARC21/XML documentation):

-> hdddmmss (hemisphere-degrees-minutes-seconds) -> hddd.dddddd (hemisphere-degrees.decimal degrees) -> +-ddd.dddddd (hemisphere[+/-]-degrees.decimal degrees) (“+” for N and E, “-“ for S and W; the plus sign is optional) -> hdddmm.mmmm (hemisphere-degrees-minutes.decimal minutes): -> hdddmmss.sss (hemisphere-degrees-minutes-seconds.decimal seconds)

degrees/minutes/seconds: hdddmmss (hemisphere-degrees-minutes-seconds) Ex. 0793235 +0793235 E0793235 || | |-> seconds (2 digits; left padded with 0) || |-> minutes (2 digits; left padded with 0) ||-> degrees (3 digits; left padded with 0) |-> start_literal (1 digit; optional)

decimal degrees: hddd.dddddd (hemisphere-degrees.decimal degrees) Ex. E079.533265 +079.533265 || |-> indicates decimal degrees ||-> degrees (3 digits; left padded with 0) |-> start_literal (1 digit; optional)

decimal minutes: hdddmm.mmmm (cardinal direction|degrees|minutes.decimal minutes) Ex. E07932.5332 || | |-> indicates decimal minutes || |-> minutes (2 digits; left padded with 0) ||-> degrees (3 digits; left padded with 0) |-> start_literal (1 digit; optional)

decimal seconds: hdddmmss.sss (hemisphere-degrees-minutes-seconds.decimal seconds) Ex. E0793235.575 || | | |-> indicates decimal seconds || | |-> seconds (2 digits; left padded with 0) || |-> minutes (2 digits; left padded with 0) ||-> degrees (3 digits; left padded with 0) |-> start_literal (1 digit; optional)

“+” for N and E (the plus sign is optional) “-“ for S and W

Definition at line 171 of file lwgeom_in_marc21.c.

                                               {
 
        char *dgr;
        char *min;
        char *sec;
        size_t literal_length;
 
        char start_character = literal[0];
        int start_literal = 0;
        double result = 0.0;
 
        const size_t numdigits_degrees = 3;
        const size_t numdigits_minutes = 2;
        const size_t numdigits_seconds = 2;
 
        POSTGIS_DEBUGF(2, "parse_geo_literal called (%s)", literal);
        POSTGIS_DEBUGF(2, "  start character: %c", start_character);
 
        literal_length = strlen(literal);
 
        if (!isdigit(start_character)) start_literal = 1;
 
        POSTGIS_DEBUGF(2, "    start_literal=%d", start_literal);
 
        dgr = palloc(sizeof(char)*numdigits_degrees+1);
        snprintf(dgr, numdigits_degrees+1, "%s", &literal[start_literal]);
 
        if (strchr(literal, '.') == NULL && strchr(literal, ',') == NULL) {
 
                POSTGIS_DEBUG(2, "  lat/lon integer coordinates detected");
                POSTGIS_DEBUGF(2, "    parsed degrees (lon/lat): %s", dgr);
 
                /* literal contain at least degrees.
                 * minutes and seconds are optional */
                result = atof(dgr);
 
                /* checks if the literal contains minutes */
                if (literal_length > (start_literal + numdigits_degrees)) {
 
                        min = palloc(sizeof(char)*numdigits_minutes+1);
                        snprintf(min, numdigits_minutes+1, "%s", &literal[start_literal+numdigits_degrees]);
                        POSTGIS_DEBUGF(2, "    parsed minutes (lon/lat): %s", min);
                        result = result + atof(min) / 60;
                        pfree(min);
 
                        /* checks if the literal contains seconds */
                        if (literal_length >= (start_literal + numdigits_degrees + numdigits_minutes)) {
 
                                sec = palloc(sizeof(char)*numdigits_seconds+1);
                                snprintf(sec, numdigits_seconds+1, "%s", &literal[start_literal+numdigits_degrees+numdigits_minutes]);
                                POSTGIS_DEBUGF(2, "    parsed seconds (lon/lat): %s", sec);
 
                                result = result + atof(sec) / 3600;
                                pfree(sec);
 
                        }
 
 
                }
 
 
        } else {
 
                POSTGIS_DEBUG(2, "  decimal coordinates detected");
 
                if (strchr(literal, ',')) {
 
                        /* changes the literal decimal sign from comma to period to avoid problems with atof.
                         * from the docs "In MARC21/XML coordinates, the decimal sign may be either a period or a comma." */
 
                        literal[literal_length-strlen(strchr(literal, ','))]='.';
                        POSTGIS_DEBUGF(2, "  decimal separator changed to '.': %s",literal);
 
                }
 
                /* checks if the literal is encoded in decimal degrees */
                if (literal[start_literal + numdigits_degrees] == '.') {
 
                        char *dec = palloc(sizeof(char)*literal_length+1);
                        snprintf(dec, literal_length+1, "%s", &literal[start_literal]);
                        result = atof(dec);
 
                        POSTGIS_DEBUGF(2, "    parsed decimal degrees: %s", dec);
                        pfree(dec);
 
                /* checks if the literal is encoded in decimal minutes */
                } else if (literal[start_literal + numdigits_degrees + numdigits_minutes] == '.') {
 
                        size_t len_decimal_minutes = literal_length - (start_literal + numdigits_degrees);
 
                        min = palloc(sizeof(char)*len_decimal_minutes+1);
                        snprintf(min, len_decimal_minutes+1, "%s", &literal[start_literal + numdigits_degrees]);
 
                        POSTGIS_DEBUGF(2, "    parsed degrees: %s", dgr);
                        POSTGIS_DEBUGF(2, "    parsed decimal minutes: %s", min);
 
                        result = atof(dgr) + (atof(min) / 60);
 
                        pfree(min);
 
                /* checks if the literal is encoded in decimal seconds */
                } else if (literal[start_literal + numdigits_degrees + numdigits_minutes + numdigits_seconds] == '.') {
 
                        size_t len_decimal_seconds = literal_length - (start_literal + numdigits_degrees + numdigits_minutes);
 
                        min = palloc(sizeof(char)*numdigits_minutes+1);
                        snprintf(min, numdigits_minutes+1, "%s", &literal[start_literal + numdigits_degrees]);
 
                        sec = palloc(sizeof(char)*len_decimal_seconds+1);
                        snprintf(sec, len_decimal_seconds+1, "%s", &literal[start_literal + numdigits_degrees + numdigits_minutes]);
 
                        result = atof(dgr) + (atof(min) / 60) + (atof(sec) / 3600);
 
                        POSTGIS_DEBUGF(2, "    parsed degrees: %s", dgr);
                        POSTGIS_DEBUGF(2, "    parsed minutes: %s", min);
                        POSTGIS_DEBUGF(2, "    parsed decimal seconds: %s", sec);
                        pfree(min);
                        pfree(sec);
 
                }
 
        }
 
        pfree(dgr);
 
        if (start_character == 'S' || start_character == 'W' || start_character == '-') {
 
                POSTGIS_DEBUGF(2, "  switching sign due to start character: '%c'", start_character);
                result = -result;
 
        }
 
        POSTGIS_DEBUGF(2, "=> parse_geo_literal returns: %.*f (in decimal degrees)", literal_length-(3+start_literal), result);
        return result;
}

References result.

Referenced by parse_marc21().

Here is the caller graph for this function: