transform_geom()

Datum transform_geom

(

PG_FUNCTION_ARGS

)

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

{
        GSERIALIZED *geom;
        GSERIALIZED *result=NULL;
        LWGEOM *lwgeom;
        projPJ input_pj, output_pj;
        char *input_proj4, *output_proj4;
        text *input_proj4_text;
        text *output_proj4_text;
        int32 result_srid ;
        char *pj_errstr;
 
        result_srid = PG_GETARG_INT32(3);
        geom = PG_GETARG_GSERIALIZED_P_COPY(0);
 
        /* Set the search path if we haven't already */
        SetPROJ4LibPath();
 
        /* Read the arguments */
        input_proj4_text  = (PG_GETARG_TEXT_P(1));
        output_proj4_text = (PG_GETARG_TEXT_P(2));
 
        /* Convert from text to cstring for libproj */
        input_proj4 = text_to_cstring(input_proj4_text);
        output_proj4 = text_to_cstring(output_proj4_text);
 
        /* make input and output projection objects */
        input_pj = lwproj_from_string(input_proj4);
        if ( input_pj == NULL )
        {
                pj_errstr = pj_strerrno(*pj_get_errno_ref());
                if ( ! pj_errstr ) pj_errstr = "";
 
                /* we need this for error reporting */
                /* pfree(input_proj4); */
                pfree(output_proj4);
                pfree(geom);
 
                elog(ERROR,
                    "transform_geom: could not parse proj4 string '%s' %s",
                    input_proj4, pj_errstr);
                PG_RETURN_NULL();
        }
        pfree(input_proj4);
 
        output_pj = lwproj_from_string(output_proj4);
 
        if ( output_pj == NULL )
        {
                pj_errstr = pj_strerrno(*pj_get_errno_ref());
                if ( ! pj_errstr ) pj_errstr = "";
 
                /* we need this for error reporting */
                /* pfree(output_proj4); */
                pj_free(input_pj);
                pfree(geom);
 
                elog(ERROR,
                        "transform_geom: couldn't parse proj4 output string: '%s': %s",
                        output_proj4, pj_errstr);
                PG_RETURN_NULL();
        }
        pfree(output_proj4);
 
        /* now we have a geometry, and input/output PJ structs. */
        lwgeom = lwgeom_from_gserialized(geom);
        lwgeom_transform(lwgeom, input_pj, output_pj);
        lwgeom->srid = result_srid;
 
        /* clean up */
        pj_free(input_pj);
        pj_free(output_pj);
 
        /* Re-compute bbox if input had one (COMPUTE_BBOX TAINTING) */
        if ( lwgeom->bbox )
        {
                lwgeom_refresh_bbox(lwgeom);
        }
 
        result = geometry_serialize(lwgeom);
 
        lwgeom_free(lwgeom);
        PG_FREE_IF_COPY(geom, 0);
 
        PG_RETURN_POINTER(result); /* new geometry */
}

References LWGEOM::bbox, geometry_serialize(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_refresh_bbox(), lwgeom_transform(), lwproj_from_string(), LWGEOM::srid, and text_to_cstring().

Here is the call graph for this function: