srs_find_planar()

static void srs_find_planar ( const char *  auth_name, const LWGEOM *  bounds, struct srs_data *  state  )

static

Definition at line 469 of file postgis/lwgeom_transform.c.

{
        int32_t srid_from = lwgeom_get_srid(bounds);
        const int32_t srid_to = 4326;
        GBOX gbox = *(lwgeom_get_bbox(bounds));
        PJ_TYPE types[1] = {PJ_TYPE_PROJECTED_CRS};
        PROJ_CRS_INFO **crs_list_ptr, **crs_list;
        int crs_count;
        PJ_CONTEXT *ctx = NULL;
 
        PROJ_CRS_LIST_PARAMETERS *params = proj_get_crs_list_parameters_create();
        params->types = types;
        params->typesCount = 1;
        params->crs_area_of_use_contains_bbox = true;
        params->bbox_valid = true;
        params->allow_deprecated = false;
 
#if POSTGIS_PROJ_VERSION >= 80100
        params->celestial_body_name = "Earth";
#endif
 
        if (srid_from != srid_to)
        {
                LWPROJ *pj;
                if (lwproj_lookup(srid_from, srid_to, &pj) == LW_FAILURE)
                        elog(ERROR, "%s: Lookup of SRID %u => %u transform failed",
                            __func__, srid_from, srid_to);
 
                box3d_transform(&gbox, pj);
        }
 
        params->west_lon_degree = gbox.xmin;
        params->south_lat_degree = gbox.ymin;
        params->east_lon_degree = gbox.xmax;
        params->north_lat_degree = gbox.ymax;
 
        crs_list = crs_list_ptr = proj_get_crs_info_list_from_database(
                                    ctx, auth_name, params, &crs_count);
 
        /* TODO, throw out really huge / dumb areas? */
 
        while (crs_list && *crs_list)
        {
                /* Read current crs and move forward one entry */
                PROJ_CRS_INFO *crs = *crs_list++;
 
                /* Read the corners of the CRS area of use */
                double area;
                double height = crs->north_lat_degree - crs->south_lat_degree;
                double width = crs->east_lon_degree - crs->west_lon_degree;
                if (width < 0.0)
                        width = 360 - (crs->west_lon_degree - crs->east_lon_degree);
                area = width * height;
 
                /* Ensure there is space in the entry list */
                srs_state_memcheck(state);
 
                /* Write the entry into the entry list and increment */
                state->entries[state->num_entries].auth_name = cstring_to_text(crs->auth_name);
                state->entries[state->num_entries].auth_code = cstring_to_text(crs->code);
                state->entries[state->num_entries].sort = area;
                state->num_entries++;
        }
 
        /* Put the list of entries into order of area size, smallest to largest */
        qsort(state->entries, state->num_entries, sizeof(struct srs_data), srs_entry_cmp);
 
        proj_crs_info_list_destroy(crs_list_ptr);
        proj_get_crs_list_parameters_destroy(params);
}

References srs_entry::auth_code, srs_entry::auth_name, box3d_transform(), srs_data::entries, LW_FAILURE, lwgeom_get_bbox(), lwgeom_get_srid(), srs_data::num_entries, srs_entry::sort, srs_entry_cmp(), srs_state_memcheck(), GBOX::xmax, GBOX::xmin, GBOX::ymax, and GBOX::ymin.

Referenced by postgis_srs_search().

Here is the call graph for this function:

Here is the caller graph for this function: