PostGIS  2.1.10dev-r@@SVN_REVISION@@
Datum LWGEOM_closestpoint3d ( PG_FUNCTION_ARGS  )

Definition at line 801 of file lwgeom_functions_basic.c.

References DIST_MIN, geometry_serialize(), lw_dist3d_distancepoint(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_is_empty(), result, and LWGEOM::srid.

802 {
804  GSERIALIZED *geom1 = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
805  GSERIALIZED *geom2 = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
806  LWGEOM *point;
807  LWGEOM *lwgeom1 = lwgeom_from_gserialized(geom1);
808  LWGEOM *lwgeom2 = lwgeom_from_gserialized(geom2);
809 
810  if (lwgeom1->srid != lwgeom2->srid)
811  {
812  elog(ERROR,"Operation on two GEOMETRIES with different SRIDs\n");
813  PG_RETURN_NULL();
814  }
815 
816  point = lw_dist3d_distancepoint(lwgeom1, lwgeom2, lwgeom1->srid, DIST_MIN);
817 
818  if (lwgeom_is_empty(point))
819  PG_RETURN_NULL();
820 
821  result = geometry_serialize(point);
822 
823  lwgeom_free(point);
824  lwgeom_free(lwgeom1);
825  lwgeom_free(lwgeom2);
826 
827  PG_FREE_IF_COPY(geom1, 0);
828  PG_FREE_IF_COPY(geom2, 1);
829  PG_RETURN_POINTER(result);
830 }
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
#define DIST_MIN
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1006
LWGEOM * lw_dist3d_distancepoint(LWGEOM *lw1, LWGEOM *lw2, int srid, int mode)
Function initializing 3dclosestpoint calculations.
Definition: measures3d.c:79
char ** result
Definition: liblwgeom.h:218
int32_t srid
Definition: liblwgeom.h:355
GSERIALIZED * geometry_serialize(LWGEOM *lwgeom)
int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members) ...
Definition: lwgeom.c:1229

Here is the call graph for this function: