57 int32 srid_to, srid_from;
59 srid_to = PG_GETARG_INT32(1);
62 elog(ERROR,
"ST_Transform: %d is an invalid target SRID",
SRID_UNKNOWN);
66 geom = PG_GETARG_GSERIALIZED_P_COPY(0);
71 PG_FREE_IF_COPY(geom, 0);
72 elog(ERROR,
"ST_Transform: Input geometry has unknown (%d) SRID",
SRID_UNKNOWN);
77 if ( srid_from == srid_to )
78 PG_RETURN_POINTER(geom);
80 postgis_initialize_cache();
81 if ( lwproj_lookup(srid_from, srid_to, &pj) ==
LW_FAILURE )
83 PG_FREE_IF_COPY(geom, 0);
84 elog(ERROR,
"ST_Transform: Failure reading projections from spatial_ref_sys.");
91 lwgeom->
srid = srid_to;
99 result = geometry_serialize(lwgeom);
101 PG_FREE_IF_COPY(geom, 0);
103 PG_RETURN_POINTER(
result);
char result[OUT_DOUBLE_BUFFER_SIZE]
int32_t gserialized_get_srid(const GSERIALIZED *g)
Extract the SRID from the serialized form (it is packed into three bytes so this is a handy function)...
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
void lwgeom_refresh_bbox(LWGEOM *lwgeom)
Drop current bbox and calculate a fresh one.
int lwgeom_transform(LWGEOM *geom, LWPROJ *pj)
Transform (reproject) a geometry in-place.
void lwgeom_free(LWGEOM *geom)
#define SRID_UNKNOWN
Unknown SRID value.