get_projinfo.c

Go to the documentation of this file.

/*!
   \file lib/gis/get_projinfo.c
 
   \brief GIS Library - Get projection info
 
   (C) 1999-2014 by the GRASS Development Team
 
   This program is free software under the GNU General Public License
   (>=v2). Read the file COPYING that comes with GRASS for details.
 */
 
#include <string.h>
#include <errno.h>
#include <unistd.h>
#include <stdio.h>
#include <grass/gis.h>
#include <grass/glocale.h>
 
#define PERMANENT "PERMANENT"
 
/*!
   \brief Gets units information for location
 
   Note: Allocated Key_Value structure should be freed by
   G_free_key_value().
 
   Prints a warning if no units information available.
 
   \return pointer to Key_Value structure with key/value pairs
   \return NULL on failure
 */
struct Key_Value *G_get_projunits(void)
{
    struct Key_Value *in_units_keys;
    char path[GPATH_MAX];
 
    G_file_name(path, "", UNIT_FILE, PERMANENT);
    if (access(path, 0) != 0) {
        if (G_projection() != PROJECTION_XY) {
            G_warning(_("<%s> file not found for location <%s>"), UNIT_FILE,
                      G_location());
        }
        return NULL;
    }
    in_units_keys = G_read_key_value_file(path);
 
    return in_units_keys;
}
 
/*!
   \brief Gets projection information for location
 
   Note: Allocated Key_Value structure should be freed by
   G_free_key_value().
 
   Prints a warning if no projection information available.
 
   \return pointer to Key_Value structure with key/value pairs
   \return NULL on failure
 */
struct Key_Value *G_get_projinfo(void)
{
    struct Key_Value *in_proj_keys, *in_epsg_keys;
    char path[GPATH_MAX];
 
    G_file_name(path, "", PROJECTION_FILE, PERMANENT);
    if (access(path, 0) != 0) {
        if (G_projection() != PROJECTION_XY) {
            G_warning(_("<%s> file not found for location <%s>"),
                      PROJECTION_FILE, G_location());
        }
        return NULL;
    }
    in_proj_keys = G_read_key_value_file(path);
 
    /* TODO: do not restrict to EPSG as the only authority */
    if ((in_epsg_keys = G_get_projepsg()) != NULL) {
        const char *epsgstr = G_find_key_value("epsg", in_epsg_keys);
        char buf[4096];
 
        sprintf(buf, "EPSG:%s", epsgstr);
        G_set_key_value("init", buf, in_proj_keys);
        G_free_key_value(in_epsg_keys);
    }
 
    return in_proj_keys;
}
 
/*!
   \brief Gets EPSG information for the current location
 
   DEPRECATED: Use G_get_projsrid() instead.
 
   Note: Allocated Key_Value structure should be freed by
   G_free_key_value().
 
   \return pointer to Key_Value structure with key/value pairs
   \return NULL when EPSG code is not defined for location
 */
 
/* superseded by G_get_projsrid(), keep for backwards compatibility */
struct Key_Value *G_get_projepsg(void)
{
    struct Key_Value *in_epsg_keys;
    char path[GPATH_MAX];
 
    G_file_name(path, "", EPSG_FILE, PERMANENT);
    if (access(path, 0) != 0) {
        if (G_projection() != PROJECTION_XY) {
            G_debug(1, "<%s> file not found for location <%s>", EPSG_FILE,
                    G_location());
        }
        return NULL;
    }
    in_epsg_keys = G_read_key_value_file(path);
 
    return in_epsg_keys;
}
 
/*!
   \brief Get WKT information for the current location
 
   \return pointer to WKT string
   \return NULL when WKT is not available for the current location
 */
 
char *G_get_projwkt(void)
{
    char *wktstring = NULL;
    char path[GPATH_MAX];
    FILE *fp;
    int n, nalloc;
    int c;
 
    G_file_name(path, "", WKT_FILE, "PERMANENT");
    if (access(path, 0) != 0) {
        if (G_projection() != PROJECTION_XY) {
            G_debug(1, "<%s> file not found for location <%s>", WKT_FILE,
                    G_location());
        }
        return NULL;
    }
 
    fp = fopen(path, "r");
    if (!fp)
        G_fatal_error(_("Unable to open input file <%s>: %s"), path,
                      strerror(errno));
 
    wktstring = G_malloc(1024 * sizeof(char));
    nalloc = 1024;
 
    n = 0;
    while (1) {
        c = fgetc(fp);
 
        if (c == EOF) {
            break;
        }
 
        if (c == '\r') { /* DOS or MacOS9 */
            c = fgetc(fp);
            if (c != EOF) {
                if (c !=
                    '\n') { /* MacOS9 - we have to return the char to stream */
                    ungetc(c, fp);
                    c = '\n';
                }
            }
            else { /* MacOS9 - we have to return the char to stream */
                ungetc(c, fp);
                c = '\n';
            }
        }
 
        if (n == nalloc) {
            wktstring = G_realloc(wktstring, nalloc + 1024);
            nalloc += 1024;
        }
 
        wktstring[n] = c;
 
        n++;
    }
 
    if (n > 0) {
        if (n == nalloc) {
            wktstring = G_realloc(wktstring, nalloc + 1);
            nalloc += 1;
        }
        wktstring[n] = '\0';
    }
    else {
        G_free(wktstring);
        wktstring = NULL;
    }
 
    if (fclose(fp) != 0)
        G_fatal_error(_("Error closing output file <%s>: %s"), path,
                      strerror(errno));
 
    if (wktstring && *wktstring)
        G_chop(wktstring);
    if (wktstring && *wktstring == '\0') {
        G_free(wktstring);
        wktstring = NULL;
    }
 
    return wktstring;
}
 
/*!
   \brief Get srid (spatial reference id) for the current location
 
   Typically an srid will be of the form authority NAME:CODE,
   e.g. EPSG:4326
 
   This srid is passed to proj_create() using PROJ or
   OSRSetFromUserInput() using GDAL. Therefore various other forms of
   srid are possible, e.g. in OSRSetFromUserInput():
 
   1. Well Known Text definition - passed on to importFromWkt().
   2. "EPSG:n" - number passed on to importFromEPSG().
   3. "EPSGA:n" - number passed on to importFromEPSGA().
   4. "AUTO:proj_id,unit_id,lon0,lat0" - WMS auto projections.
   5. "urn:ogc:def:crs:EPSG::n" - ogc urns
   6. PROJ.4 definitions - passed on to importFromProj4().
   7. filename - file read for WKT, XML or PROJ.4 definition.
   8. well known name accepted by SetWellKnownGeogCS(), such as NAD27, NAD83,
   WGS84 or WGS72.
   9. "IGNF:xxxx", "ESRI:xxxx", etc. from definitions from the PROJ database;
   10. PROJJSON (PROJ >= 6.2)
 
   \return pointer to srid string
   \return NULL when srid is not available for the current location
 */
 
char *G_get_projsrid(void)
{
    char *sridstring = NULL;
    char path[GPATH_MAX];
    FILE *fp;
    int n, nalloc;
    int c;
 
    G_file_name(path, "", SRID_FILE, "PERMANENT");
    if (access(path, 0) != 0) {
        if (G_projection() != PROJECTION_XY) {
            struct Key_Value *projepsg;
            const char *epsg_num;
 
            G_debug(1, "<%s> file not found for location <%s>", SRID_FILE,
                    G_location());
 
            /* for backwards compatibility, check if PROJ_EPSG exists */
            if ((projepsg = G_get_projepsg()) != NULL) {
                epsg_num = G_find_key_value("epsg", projepsg);
                if (*epsg_num) {
                    G_debug(1, "Using <%s> file instead for location <%s>",
                            EPSG_FILE, G_location());
                    G_asprintf(&sridstring, "EPSG:%s", epsg_num);
                    G_free_key_value(projepsg);
 
                    return sridstring;
                }
            }
        }
        return NULL;
    }
 
    fp = fopen(path, "r");
    if (!fp)
        G_fatal_error(_("Unable to open input file <%s>: %s"), path,
                      strerror(errno));
 
    sridstring = G_malloc(1024 * sizeof(char));
    nalloc = 1024;
 
    n = 0;
    while (1) {
        c = fgetc(fp);
 
        if (c == EOF) {
            break;
        }
 
        if (c == '\r') { /* DOS or MacOS9 */
            c = fgetc(fp);
            if (c != EOF) {
                if (c !=
                    '\n') { /* MacOS9 - we have to return the char to stream */
                    ungetc(c, fp);
                    c = '\n';
                }
            }
            else { /* MacOS9 - we have to return the char to stream */
                ungetc(c, fp);
                c = '\n';
            }
        }
 
        if (n == nalloc) {
            sridstring = G_realloc(sridstring, nalloc + 1024);
            nalloc += 1024;
        }
 
        sridstring[n] = c;
 
        n++;
    }
 
    if (n > 0) {
        if (n == nalloc) {
            sridstring = G_realloc(sridstring, nalloc + 1);
            nalloc += 1;
        }
        sridstring[n] = '\0';
    }
    else {
        G_free(sridstring);
        sridstring = NULL;
    }
 
    if (fclose(fp) != 0)
        G_fatal_error(_("Error closing output file <%s>: %s"), path,
                      strerror(errno));
 
    if (sridstring && *sridstring)
        G_chop(sridstring);
    if (sridstring && *sridstring == '\0') {
        G_free(sridstring);
        sridstring = NULL;
    }
 
    return sridstring;
}