lwgeom_out_marc21.c

Go to the documentation of this file.

/**********************************************************************
 *
 * PostGIS - Spatial Types for PostgreSQL
 * http://postgis.net
 *
 * PostGIS is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 2 of the License, or
 * (at your option) any later version.
 *
 * PostGIS is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with PostGIS.  If not, see <http://www.gnu.org/licenses/>.
 **********************************************************************/
 
#include "postgres.h"
#include "fmgr.h"
#include "utils/builtins.h"
 
#include "../postgis_config.h"
#include "liblwgeom.h"
#include "lwgeom_transform.h"
 
#include "stringbuffer.h"
#include "liblwgeom_internal.h"
 
#define MARC21_NS ((char *) "http://www.loc.gov/MARC21/slim")
/**********************************************************************
 * Ability to serialise MARC21/XML Records from geometries. Coordinates
 * are returned in decimal degrees.
 *
 * MARC21/XML version supported: 1.1
 * MARC21/XML Cartographic Mathematical Data Definition:
 *    https://www.loc.gov/marc/bibliographic/bd034.html
 *
 * Copyright (C) 2021 University of Münster (WWU), Germany
 * Written by Jim Jones <jim.jones@uni-muenster.de>
 *
 **********************************************************************/
 
 
static int gbox_to_marc21_sb(const GBOX box, const char* format, stringbuffer_t *sb);
static int is_format_valid(const char* format);
lwvarlena_t* lwgeom_to_marc21(const LWGEOM *geom, const char* format);
 
PG_FUNCTION_INFO_V1(ST_AsMARC21);
Datum ST_AsMARC21(PG_FUNCTION_ARGS) {
        lwvarlena_t *marc21;
        int32_t srid;
        LWPROJ *lwproj;
        LWGEOM *lwgeom;
        //uint8_t is_latlong;
        GSERIALIZED *gs = PG_GETARG_GSERIALIZED_P(0);
        text *format_text_input =  PG_GETARG_TEXT_P(1);
        const char *format = text_to_cstring(format_text_input);
 
        srid = gserialized_get_srid(gs);
 
        if (srid == SRID_UNKNOWN) {
 
                PG_FREE_IF_COPY(gs, 0);
                lwpgerror("ST_AsMARC21: Input geometry has unknown (%d) SRID", srid);
                PG_RETURN_NULL();
 
        }
 
        if (lwproj_lookup(srid, srid, &lwproj) == LW_FAILURE) {
 
                PG_FREE_IF_COPY(gs, 0);
                lwpgerror("ST_AsMARC21: Failure reading projections from spatial_ref_sys.");
                PG_RETURN_NULL();
 
        }
 
//      if (GetLWPROJ(srid, srid, &lwproj) == LW_FAILURE) {
//
//              PG_FREE_IF_COPY(gs, 0);
//              lwpgerror("ST_AsMARC21: Failure reading projections from spatial_ref_sys.");
//              PG_RETURN_NULL();
//
//      }
 
//#if POSTGIS_PROJ_VERSION < 61
//      is_latlong = pj_is_latlong(lwproj->pj_from);
//#else
//      is_latlong = lwproj->source_is_latlong;
//#endif
 
        //is_latlong = lwproj_is_latlong(lwproj);
 
        if (!lwproj_is_latlong(lwproj)) {
 
                PG_FREE_IF_COPY(gs, 0);
                lwpgerror("ST_AsMARC21: Unsupported SRID (%d). Only lon/lat coordinate systems are supported in MARC21/XML Documents.", srid);
                PG_RETURN_NULL();
        }
 
 
        lwgeom = lwgeom_from_gserialized(gs);
        marc21 = lwgeom_to_marc21(lwgeom, format);
 
        if (marc21)     PG_RETURN_TEXT_P(marc21);
 
        PG_RETURN_NULL();
}
 
lwvarlena_t*
lwgeom_to_marc21(const LWGEOM *geom, const char* format) {
 
        stringbuffer_t *sb;
        GBOX box;
        lwvarlena_t *v;
 
        POSTGIS_DEBUGF(2, "lwgeom_to_marc21 called: %s with format %s ", lwtype_name(lwgeom_get_type(geom)), format);
 
        if (lwgeom_is_empty(geom))      return NULL;
        if (is_format_valid(format)==LW_FALSE) lwerror("invalid output format: \"%s\"", format);
 
        POSTGIS_DEBUG(3, "creating stringbuffer");
        sb = stringbuffer_create();
 
        POSTGIS_DEBUGF(3, "opening MARC21/XML record: %s", lwtype_name(lwgeom_get_type(geom)));
 
        if (stringbuffer_aprintf(sb, "<record xmlns=\"%s\">", MARC21_NS) < 0) {
 
                stringbuffer_destroy(sb);
                return NULL;
 
        }
 
 
        if (lwgeom_is_collection(geom)) {
 
                LWCOLLECTION *coll = (LWCOLLECTION*) geom;
 
                POSTGIS_DEBUGF(3, "  collection detected: %s", lwtype_name(lwgeom_get_type(geom)));
 
                for (uint32_t i = 0; i < coll->ngeoms; i++) {
 
                        if (lwgeom_calculate_gbox(coll->geoms[i], &box) == LW_FAILURE) {
 
                                stringbuffer_destroy(sb);
                                lwpgerror("failed to calculate bbox for a geometry in the collection: %s", lwtype_name(lwgeom_get_type(coll->geoms[i])));
 
                        }
 
                        if (gbox_to_marc21_sb(box, format, sb) == LW_FAILURE) {
 
                                stringbuffer_destroy(sb);
                                lwpgerror("failed to create MARC21/XML for a geometry in the collection: %s",   lwtype_name(lwgeom_get_type(coll->geoms[i])));
 
                        }
 
                }
 
        } else {
 
                POSTGIS_DEBUGF(3, "  calculating gbox: %s", lwtype_name(lwgeom_get_type(geom)));
 
                if (lwgeom_calculate_gbox(geom, &box) == LW_FAILURE) {
 
                        stringbuffer_destroy(sb);
                        lwpgerror("failed to calculate bbox for %s", lwtype_name(lwgeom_get_type(geom)));
 
                }
 
                POSTGIS_DEBUGF(3, "  creating MARC21/XML datafield: %s", lwtype_name(lwgeom_get_type(geom)));
 
                if (gbox_to_marc21_sb(box, format, sb) == LW_FAILURE) {
 
                        stringbuffer_destroy(sb);
                        lwpgerror("failed to create MARC21/XML for %s", lwtype_name(lwgeom_get_type(geom)));
 
                }
 
        }
 
        POSTGIS_DEBUG(3, "  closing MARC21/XML record");
 
        if (stringbuffer_aprintf(sb, "</record>") < 0) {
 
                stringbuffer_destroy(sb);
                return LW_FAILURE;
 
        }
 
        v = stringbuffer_getvarlenacopy(sb);
        stringbuffer_destroy(sb);
 
        return v;
}
 
 
static int is_format_valid(const char* format){
 
        char *int_part;
        char *dec_part = strchr(format, '.');
        if(!dec_part) dec_part = strchr(format, ',');
 
        if(!dec_part) {
 
                if (strcmp(format, "hdddmmss") && strcmp(format, "dddmmss")) {
 
                        return LW_FALSE;
 
                }
 
        } else {
 
                if(strlen(dec_part)<2)  return LW_FALSE;
 
                int_part = palloc(sizeof(char)*strlen(format));
                memcpy(int_part, &format[0], strlen(format) - strlen(dec_part));
                int_part[strlen(format) - strlen(dec_part)]='\0';
 
                if (strcmp(int_part,"hddd") && strcmp(int_part,"ddd") &&
                        strcmp(int_part,"hdddmm") && strcmp(int_part,"dddmm") &&
                        strcmp(int_part,"hdddmmss") && strcmp(int_part,"dddmmss")) {
 
                        pfree(int_part);
                        return LW_FALSE;
 
                }
 
                for (size_t i = 1; i < strlen(dec_part); i++) {
 
                        if(dec_part[i]!=int_part[strlen(int_part)-1]) {
 
                                pfree(int_part);
                                return LW_FALSE;
                        }
                }
 
                pfree(int_part);
 
        }
 
        return LW_TRUE;
 
}
 
static int corner_to_subfield_sb(stringbuffer_t *sb, double decimal_degrees, const char*format, char subfield) {
 
        char cardinal_direction;
        char decimal_separator;
 
        int degrees = (int) decimal_degrees;
        double minutes = fabs((decimal_degrees-degrees)*60);
        double seconds = fabs((minutes-(int)minutes) *60);
 
        int has_cardinal_direction = 0;
        int num_decimals = 0;
        char* res = palloc(sizeof(char)*strlen(format)+2);
 
        /* size of the buffer for the output snprintf calls.
         * the output strings must have the same length as the format.
         * +1 to make room for the null character '\0' */
        size_t buffer_size = strlen(format)+1;
 
        /* +1 one digit to the buffer size in case of negative
         * numbers to account for the "-" sign */
        if(degrees < 0) buffer_size = buffer_size+1;
 
 
        POSTGIS_DEBUGF(2,"corner_to_subfield_sb called with coordinates: %f and format: %s",decimal_degrees,format);
 
        if((int)(seconds + 0.5)>=60) {
 
                seconds = seconds-60;
                minutes = minutes +1;
 
        }
 
 
        if(strchr(format, '.')) {
                num_decimals = strlen(strchr(format, '.'))-1;
                decimal_separator = '.';
        }
 
        if(strchr(format, ',')) {
                num_decimals = strlen(strchr(format, ','))-1;
                decimal_separator = ',';
        }
 
        if(format[0]=='h'){
 
                has_cardinal_direction = 1;
 
                if(subfield=='d'||subfield=='e'){
 
                        if(decimal_degrees>0){
 
                                cardinal_direction='E';
 
                        } else {
 
                                cardinal_direction='W';
                                degrees=abs(degrees);
                                decimal_degrees= fabs(decimal_degrees);
 
                        }
                }
 
                if(subfield=='f'||subfield=='g'){
 
                        if(decimal_degrees>0){
 
                                cardinal_direction='N';
 
                        } else {
 
                                cardinal_direction='S';
                                degrees=abs(degrees);
                                decimal_degrees= fabs(decimal_degrees);
 
                        }
                }
 
 
        }
 
        if(format[3+has_cardinal_direction]=='.' || format[3+has_cardinal_direction]==',' )     {
 
                int pad_degrees = (int)strlen(format);
 
                if(decimal_degrees <0 && decimal_degrees>-100) pad_degrees=strlen(format)+1;
 
                if(has_cardinal_direction) pad_degrees=pad_degrees-1;
 
                snprintf(res,buffer_size,"%0*.*f",pad_degrees,num_decimals,decimal_degrees);
 
 
        } else if(format[5+has_cardinal_direction]=='.' || format[5+has_cardinal_direction]==',' )      {
 
                int pad_minutes = 0;
 
                if(minutes<10) pad_minutes = (int)strlen(format)-has_cardinal_direction-3;
 
                snprintf(res,buffer_size,"%.3d%0*.*f",degrees,pad_minutes,num_decimals,fabs(minutes));
 
        }
 
        else if(format[7+has_cardinal_direction]=='.' || format[7+has_cardinal_direction]==',') {
 
                /*
                 * decimal seconds
                 */
 
                int pad_seconds = 0;
 
                if(seconds<10) pad_seconds = (int) strlen(format)-has_cardinal_direction-5;
 
                snprintf(res,buffer_size,"%.3d%.2d%0*.*f",degrees,(int)minutes,pad_seconds,num_decimals,fabs(seconds));
 
        } else {
 
                snprintf(res,buffer_size,"%.3d%.2d%.2d",degrees,(int)minutes,(int)(seconds + 0.5));
 
        }
 
 
        if(decimal_separator==','){
 
                res[strlen(res)-num_decimals-1] = ',';
 
        }
 
        if(has_cardinal_direction){
 
                if (stringbuffer_aprintf(sb, "<subfield code=\"%c\">%c%s</subfield>", subfield, cardinal_direction, res) < 0) return LW_FAILURE;
 
        } else {
 
                if (stringbuffer_aprintf(sb, "<subfield code=\"%c\">%s</subfield>",     subfield, res) < 0)     return LW_FAILURE;
 
        }
 
 
        pfree(res);
        return LW_SUCCESS;
 
 
 
}
 
static int gbox_to_marc21_sb(const GBOX box, const char* format, stringbuffer_t *sb) {
 
        POSTGIS_DEBUG(2, "gbox_to_marc21_sb called");
 
        if (stringbuffer_aprintf(sb, "<datafield tag=\"034\" ind1=\"1\" ind2=\" \">") < 0)      return LW_FAILURE;
        if (stringbuffer_aprintf(sb, "<subfield code=\"a\">a</subfield>") < 0) return LW_FAILURE;
        if (!corner_to_subfield_sb(sb, box.xmin, format, 'd')) return LW_FAILURE;
        if (!corner_to_subfield_sb(sb, box.xmax, format, 'e')) return LW_FAILURE;
        if (!corner_to_subfield_sb(sb, box.ymax, format, 'f')) return LW_FAILURE;
        if (!corner_to_subfield_sb(sb, box.ymin, format, 'g')) return LW_FAILURE;
        if (stringbuffer_aprintf(sb, "</datafield>") < 0) return LW_FAILURE;
 
        POSTGIS_DEBUG(2, "=> gbox_to_marc21_sb returns LW_SUCCESS");
 
        return LW_SUCCESS;
}