gml_is_srs_axis_order_gis_friendly()

static int gml_is_srs_axis_order_gis_friendly ( int32_t  srid )

static

Return 1 if the SRS definition from the authority has a GIS friendly order, that is easting,northing.

This is typically true for most projected CRS (but not all!), and this is false for EPSG geographic CRS.

Definition at line 427 of file lwgeom_in_gml.c.

{
        char *srtext;
        char query[256];
        int is_axis_order_gis_friendly, err;
 
        if (SPI_OK_CONNECT != SPI_connect ())
                lwpgerror("gml_is_srs_axis_order_gis_friendly: could not connect to SPI manager");
 
        sprintf(query, "SELECT srtext \
                        FROM spatial_ref_sys WHERE srid='%d'", srid);
 
        err = SPI_exec(query, 1);
        if (err < 0) lwpgerror("gml_is_srs_axis_order_gis_friendly: error executing query %d", err);
 
        /* No entry in spatial_ref_sys */
        if (SPI_processed <= 0)
        {
                SPI_finish();
                return -1;
        }
 
        srtext = SPI_getvalue(SPI_tuptable->vals[0], SPI_tuptable->tupdesc, 1);
 
        is_axis_order_gis_friendly = 1;
        if (srtext && srtext[0] != '\0')
        {
                char* ptr;
                char* srtext_horizontal = (char*) malloc(strlen(srtext) + 1);
                strcpy(srtext_horizontal, srtext);
 
                /* Remove the VERT_CS part if we are in a COMPD_CS */
                ptr = strstr(srtext_horizontal, ",VERT_CS[");
                if (ptr)
                        *ptr = '\0';
 
                if( strstr(srtext_horizontal, "AXIS[") == NULL &&
                    strstr(srtext_horizontal, "GEOCCS[") == NULL )
                {
                        /* If there is no axis definition, then due to how GDAL < 3
                        * generated the WKT, this means that the axis order is not
                        * GIS friendly */
                        is_axis_order_gis_friendly = 0;
                }
                else if( strstr(srtext_horizontal,
                           "AXIS[\"Latitude\",NORTH],AXIS[\"Longitude\",EAST]") != NULL )
                {
                        is_axis_order_gis_friendly = 0;
                }
                else if( strstr(srtext_horizontal,
                           "AXIS[\"Northing\",NORTH],AXIS[\"Easting\",EAST]") != NULL )
                {
                        is_axis_order_gis_friendly = 0;
                }
                else if( strstr(srtext_horizontal,
                           "AXIS[\"geodetic latitude (Lat)\",north,ORDER[1]") != NULL )
                {
                        is_axis_order_gis_friendly = 0;
                }
 
                free(srtext_horizontal);
        }
        SPI_finish();
 
        return is_axis_order_gis_friendly;
}

References free(), and malloc().

Referenced by parse_gml_srs().

Here is the call graph for this function:

Here is the caller graph for this function: