Datum geometry_distance_spheroid(PG_FUNCTION_ARGS)
{
GSERIALIZED *geom1 = PG_GETARG_GSERIALIZED_P(0);
GSERIALIZED *geom2 = PG_GETARG_GSERIALIZED_P(1);
SPHEROID *sphere = (SPHEROID *)PG_GETARG_POINTER(2);
int type1 = gserialized_get_type(geom1);
int type2 = gserialized_get_type(geom2);
bool use_spheroid = PG_GETARG_BOOL(3);
LWGEOM *lwgeom1, *lwgeom2;
double distance;
/* Calculate some other parameters on the spheroid */
spheroid_init(sphere, sphere->a, sphere->b);
/* Catch sphere special case and re-jig spheroid appropriately */
if ( ! use_spheroid )
{
sphere->a = sphere->b = sphere->radius;
}
if ( ! ( type1 == POINTTYPE || type1 == LINETYPE || type1 == POLYGONTYPE ||
type1 == MULTIPOINTTYPE || type1 == MULTILINETYPE || type1 == MULTIPOLYGONTYPE ))
{
elog(ERROR, "geometry_distance_spheroid: Only point/line/polygon supported.\n");
PG_RETURN_NULL();
}
if ( ! ( type2 == POINTTYPE || type2 == LINETYPE || type2 == POLYGONTYPE ||
type2 == MULTIPOINTTYPE || type2 == MULTILINETYPE || type2 == MULTIPOLYGONTYPE ))
{
elog(ERROR, "geometry_distance_spheroid: Only point/line/polygon supported.\n");
PG_RETURN_NULL();
}
if (gserialized_get_srid(geom1) != gserialized_get_srid(geom2))
{
elog(ERROR, "geometry_distance_spheroid: Operation on two GEOMETRIES with different SRIDs\n");
PG_RETURN_NULL();
}
/* Get #LWGEOM structures */
lwgeom1 = lwgeom_from_gserialized(geom1);
lwgeom2 = lwgeom_from_gserialized(geom2);
/* We are going to be calculating geodetic distances */
lwgeom_set_geodetic(lwgeom1, LW_TRUE);
lwgeom_set_geodetic(lwgeom2, LW_TRUE);
distance = lwgeom_distance_spheroid(lwgeom1, lwgeom2, sphere, 0.0);
PG_RETURN_FLOAT8(distance);
}
Datum geography_from_geometry(PG_FUNCTION_ARGS)
{
GSERIALIZED *geom = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
LWGEOM *lwgeom = NULL;
GSERIALIZED *g_ser = NULL;
geography_valid_type(gserialized_get_type(geom));
lwgeom = lwgeom_from_gserialized(geom);
/* Force default SRID */
if ( (int)lwgeom->srid <= 0 )
{
lwgeom->srid = SRID_DEFAULT;
}
/* Error on any SRID != default */
if ( lwgeom->srid != SRID_DEFAULT )
{
ereport(ERROR, (
errcode(ERRCODE_INVALID_PARAMETER_VALUE),
errmsg("Only SRID %d is currently supported in geography.", SRID_DEFAULT)));
}
/* Check if the geography has valid coordinate range. */
if ( lwgeom_check_geodetic(lwgeom) == LW_FALSE )
{
ereport(ERROR, (
errcode(ERRCODE_INVALID_PARAMETER_VALUE),
errmsg("Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type" )));
}
/*
** Serialize our lwgeom and set the geodetic flag so subsequent
** functions do the right thing.
*/
lwgeom_set_geodetic(lwgeom, true);
/* Recalculate the boxes after re-setting the geodetic bit */
lwgeom_drop_bbox(lwgeom);
lwgeom_add_bbox(lwgeom);
g_ser = geography_serialize(lwgeom);
/*
** Replace the unaligned lwgeom with a new aligned one based on GSERIALIZED.
*/
lwgeom_free(lwgeom);
PG_FREE_IF_COPY(geom, 0);
PG_RETURN_POINTER(g_ser);
}
Datum geography_from_geometry(PG_FUNCTION_ARGS)
{
GSERIALIZED *geom = (GSERIALIZED*)PG_DETOAST_DATUM_COPY(PG_GETARG_DATUM(0));
LWGEOM *lwgeom = NULL;
GSERIALIZED *g_ser = NULL;
geography_valid_type(gserialized_get_type(geom));
lwgeom = lwgeom_from_gserialized(geom);
/* Force default SRID */
if ( (int)lwgeom->srid <= 0 )
{
lwgeom->srid = SRID_DEFAULT;
}
/* Error on any SRID != default */
srid_is_latlong(fcinfo, lwgeom->srid);
/* Force the geometry to have valid geodetic coordinate range. */
lwgeom_nudge_geodetic(lwgeom);
if ( lwgeom_force_geodetic(lwgeom) == LW_TRUE )
{
ereport(NOTICE, (
errmsg_internal("Coordinate values were coerced into range [-180 -90, 180 90] for GEOGRAPHY" ))
);
}
/*
** Serialize our lwgeom and set the geodetic flag so subsequent
** functions do the right thing.
*/
lwgeom_set_geodetic(lwgeom, true);
/* Recalculate the boxes after re-setting the geodetic bit */
lwgeom_drop_bbox(lwgeom);
lwgeom_add_bbox(lwgeom);
g_ser = geography_serialize(lwgeom);
/*
** Replace the unaligned lwgeom with a new aligned one based on GSERIALIZED.
*/
lwgeom_free(lwgeom);
PG_FREE_IF_COPY(geom, 0);
PG_RETURN_POINTER(g_ser);
}
GSERIALIZED* gserialized_geography_from_lwgeom(LWGEOM *lwgeom, int32 geog_typmod)
{
GSERIALIZED *g_ser = NULL;
/* Set geodetic flag */
lwgeom_set_geodetic(lwgeom, true);
/* Check that this is a type we can handle */
geography_valid_type(lwgeom->type);
/* Check that the coordinates are in range */
if ( lwgeom_check_geodetic(lwgeom) == LW_FALSE )
{
if ( (! lwgeom_nudge_geodetic(lwgeom)) || lwgeom_check_geodetic(lwgeom) == LW_FALSE )
{
ereport(ERROR, (
errcode(ERRCODE_INVALID_PARAMETER_VALUE),
errmsg("Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type" )));
}
}
/* Force default SRID to the default */
if ( (int)lwgeom->srid <= 0 )
lwgeom->srid = SRID_DEFAULT;
/*
** Serialize our lwgeom and set the geodetic flag so subsequent
** functions do the right thing.
*/
g_ser = geography_serialize(lwgeom);
/* Check for typmod agreement */
if ( geog_typmod >= 0 )
{
postgis_valid_typmod(g_ser, geog_typmod);
POSTGIS_DEBUG(3, "typmod and geometry were consistent");
}
else
{
POSTGIS_DEBUG(3, "typmod was -1");
}
return g_ser;
}
void
lwgeom_set_geodetic(LWGEOM *geom, int value)
{
LWPOINT *pt;
LWLINE *ln;
LWPOLY *ply;
LWCOLLECTION *col;
int i;
FLAGS_SET_GEODETIC(geom->flags, value);
if ( geom->bbox )
FLAGS_SET_GEODETIC(geom->bbox->flags, value);
switch(geom->type)
{
case POINTTYPE:
pt = (LWPOINT*)geom;
if ( pt->point )
FLAGS_SET_GEODETIC(pt->point->flags, value);
break;
case LINETYPE:
ln = (LWLINE*)geom;
if ( ln->points )
FLAGS_SET_GEODETIC(ln->points->flags, value);
break;
case POLYGONTYPE:
ply = (LWPOLY*)geom;
for ( i = 0; i < ply->nrings; i++ )
FLAGS_SET_GEODETIC(ply->rings[i]->flags, value);
break;
case MULTIPOINTTYPE:
case MULTILINETYPE:
case MULTIPOLYGONTYPE:
case COLLECTIONTYPE:
col = (LWCOLLECTION*)geom;
for ( i = 0; i < col->ngeoms; i++ )
lwgeom_set_geodetic(col->geoms[i], value);
break;
default:
lwerror("lwgeom_set_geodetic: unsupported geom type: %s", lwtype_name(geom->type));
return;
}
}
GSERIALIZED* gserialized_geography_from_lwgeom(LWGEOM *lwgeom, int32 geog_typmod)
{
GSERIALIZED *g_ser = NULL;
/* Set geodetic flag */
lwgeom_set_geodetic(lwgeom, true);
/* Check that this is a type we can handle */
geography_valid_type(lwgeom->type);
/* Force the geometry to have valid geodetic coordinate range. */
lwgeom_nudge_geodetic(lwgeom);
if ( lwgeom_force_geodetic(lwgeom) == LW_TRUE )
{
ereport(NOTICE, (
errmsg_internal("Coordinate values were coerced into range [-180 -90, 180 90] for GEOGRAPHY" ))
);
}
/* Force default SRID to the default */
if ( (int)lwgeom->srid <= 0 )
lwgeom->srid = SRID_DEFAULT;
/*
** Serialize our lwgeom and set the geodetic flag so subsequent
** functions do the right thing.
*/
g_ser = geography_serialize(lwgeom);
/* Check for typmod agreement */
if ( geog_typmod >= 0 )
{
postgis_valid_typmod(g_ser, geog_typmod);
POSTGIS_DEBUG(3, "typmod and geometry were consistent");
}
else
{
POSTGIS_DEBUG(3, "typmod was -1");
}
return g_ser;
}
Datum geometry_from_geography(PG_FUNCTION_ARGS)
{
LWGEOM *lwgeom = NULL;
GSERIALIZED *ret = NULL;
GSERIALIZED *g_ser = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
lwgeom = lwgeom_from_gserialized(g_ser);
/* Recalculate the boxes after re-setting the geodetic bit */
lwgeom_set_geodetic(lwgeom, false);
lwgeom_drop_bbox(lwgeom);
lwgeom_add_bbox(lwgeom);
/* We want "geometry" to think all our "geography" has an SRID, and the
implied SRID is the default, so we fill that in if our SRID is actually unknown. */
if ( (int)lwgeom->srid <= 0 )
lwgeom->srid = SRID_DEFAULT;
ret = geometry_serialize(lwgeom);
lwgeom_free(lwgeom);
PG_RETURN_POINTER(ret);
}
Datum geography_in(PG_FUNCTION_ARGS)
{
char *str = PG_GETARG_CSTRING(0);
/* Datum geog_oid = PG_GETARG_OID(1); Not needed. */
int32 geog_typmod = PG_GETARG_INT32(2);
LWGEOM_PARSER_RESULT lwg_parser_result;
LWGEOM *lwgeom = NULL;
GSERIALIZED *g_ser = NULL;
lwgeom_parser_result_init(&lwg_parser_result);
/* Empty string. */
if ( str[0] == '\0' )
ereport(ERROR,(errmsg("parse error - invalid geometry")));
/* WKB? Let's find out. */
if ( str[0] == '0' )
{
/* TODO: 20101206: No parser checks! This is inline with current 1.5 behavior, but needs discussion */
lwgeom = lwgeom_from_hexwkb(str, LW_PARSER_CHECK_NONE);
/* Error out if something went sideways */
if ( ! lwgeom )
ereport(ERROR,(errmsg("parse error - invalid geometry")));
}
/* WKT then. */
else
{
if ( lwgeom_parse_wkt(&lwg_parser_result, str, LW_PARSER_CHECK_ALL) == LW_FAILURE )
PG_PARSER_ERROR(lwg_parser_result);
lwgeom = lwg_parser_result.geom;
}
/* Set geodetic flag */
lwgeom_set_geodetic(lwgeom, true);
/* Check that this is a type we can handle */
geography_valid_type(lwgeom->type);
/* Check that the coordinates are in range */
if ( lwgeom_check_geodetic(lwgeom) == LW_FALSE )
{
ereport(ERROR, (
errcode(ERRCODE_INVALID_PARAMETER_VALUE),
errmsg("Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type" )));
}
/* Force default SRID to the default */
if ( (int)lwgeom->srid <= 0 )
lwgeom->srid = SRID_DEFAULT;
if ( geog_typmod >= 0 )
{
postgis_valid_typmod(lwgeom, geog_typmod);
POSTGIS_DEBUG(3, "typmod and geometry were consistent");
}
else
{
POSTGIS_DEBUG(3, "typmod was -1");
}
/*
** Serialize our lwgeom and set the geodetic flag so subsequent
** functions do the right thing.
*/
g_ser = geography_serialize(lwgeom);
/* Clean up temporary object */
lwgeom_free(lwgeom);
PG_RETURN_POINTER(g_ser);
}
