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

Definition at line 836 of file lwgeom_functions_basic.c.

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

837 {
839  GSERIALIZED *geom1 = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
840  GSERIALIZED *geom2 = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
841  LWGEOM *theline;
842  LWGEOM *lwgeom1 = lwgeom_from_gserialized(geom1);
843  LWGEOM *lwgeom2 = lwgeom_from_gserialized(geom2);
844 
845  if (lwgeom1->srid != lwgeom2->srid)
846  {
847  elog(ERROR,"Operation on two GEOMETRIES with different SRIDs\n");
848  PG_RETURN_NULL();
849  }
850 
851  theline = lw_dist3d_distanceline(lwgeom1, lwgeom2, lwgeom1->srid, DIST_MIN);
852 
853  if (lwgeom_is_empty(theline))
854  PG_RETURN_NULL();
855 
856  result = geometry_serialize(theline);
857 
858  lwgeom_free(theline);
859  lwgeom_free(lwgeom1);
860  lwgeom_free(lwgeom2);
861 
862  PG_FREE_IF_COPY(geom1, 0);
863  PG_FREE_IF_COPY(geom2, 1);
864  PG_RETURN_POINTER(result);
865 }
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
char ** result
Definition: liblwgeom.h:218
int32_t srid
Definition: liblwgeom.h:355
GSERIALIZED * geometry_serialize(LWGEOM *lwgeom)
LWGEOM * lw_dist3d_distanceline(LWGEOM *lw1, LWGEOM *lw2, int srid, int mode)
Function initializing 3dshortestline and 3dlongestline calculations.
Definition: measures3d.c:30
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: