double
lwgeom_interpolate_point(const LWGEOM *lwin, const LWPOINT *lwpt)
{
POINT4D p, p_proj;
double ret = 0.0;
if ( ! lwin )
lwerror("lwgeom_interpolate_point: null input geometry!");
if ( ! lwgeom_has_m(lwin) )
lwerror("Input geometry does not have a measure dimension");
if ( lwgeom_is_empty(lwin) || lwpoint_is_empty(lwpt) )
lwerror("Input geometry is empty");
switch ( lwin->type )
{
case LINETYPE:
{
LWLINE *lwline = lwgeom_as_lwline(lwin);
lwpoint_getPoint4d_p(lwpt, &p);
ret = ptarray_locate_point(lwline->points, &p, NULL, &p_proj);
ret = p_proj.m;
break;
}
default:
lwerror("This function does not accept %s geometries.", lwtype_name(lwin->type));
}
return ret;
}
double
lwpoint_get_y(const LWPOINT *point)
{
POINT4D pt;
if ( lwpoint_is_empty(point) )
lwerror("lwpoint_get_y called with empty geometry");
getPoint4d_p(point->point, 0, &pt);
return pt.y;
}
double
lwpoint_get_m(const LWPOINT *point)
{
POINT4D pt;
if ( lwpoint_is_empty(point) )
lwerror("lwpoint_get_m called with empty geometry");
if ( ! FLAGS_GET_M(point->flags) )
lwerror("lwpoint_get_m called without m dimension");
getPoint4d_p(point->point, 0, &pt);
return pt.m;
}
static size_t
asgeojson_point_buf(const LWPOINT *point, char *srs, char *output, GBOX *bbox, int precision)
{
char *ptr = output;
ptr += sprintf(ptr, "{\"type\":\"Point\",");
if (srs) ptr += asgeojson_srs_buf(ptr, srs);
if (bbox) ptr += asgeojson_bbox_buf(ptr, bbox, FLAGS_GET_Z(point->flags), precision);
ptr += sprintf(ptr, "\"coordinates\":");
if ( lwpoint_is_empty(point) )
ptr += sprintf(ptr, "[]");
ptr += pointArray_to_geojson(point->point, ptr, precision);
ptr += sprintf(ptr, "}");
return (ptr-output);
}
static size_t
asgeojson_point_size(const LWPOINT *point, char *srs, GBOX *bbox, int precision)
{
int size;
size = pointArray_geojson_size(point->point, precision);
size += sizeof("{'type':'Point',");
size += sizeof("'coordinates':}");
if ( lwpoint_is_empty(point) )
size += 2; /* [] */
if (srs) size += asgeojson_srs_size(srs);
if (bbox) size += asgeojson_bbox_size(FLAGS_GET_Z(point->flags), precision);
return size;
}
int lwgeom_is_empty(const LWGEOM *geom)
{
int result = LW_FALSE;
LWDEBUGF(4, "lwgeom_is_empty: got type %s",
lwtype_name(geom->type));
switch (geom->type)
{
case POINTTYPE:
return lwpoint_is_empty((LWPOINT*)geom);
break;
case LINETYPE:
return lwline_is_empty((LWLINE*)geom);
break;
case CIRCSTRINGTYPE:
return lwcircstring_is_empty((LWCIRCSTRING*)geom);
break;
case POLYGONTYPE:
return lwpoly_is_empty((LWPOLY*)geom);
break;
case TRIANGLETYPE:
return lwtriangle_is_empty((LWTRIANGLE*)geom);
break;
case MULTIPOINTTYPE:
case MULTILINETYPE:
case MULTIPOLYGONTYPE:
case COMPOUNDTYPE:
case CURVEPOLYTYPE:
case MULTICURVETYPE:
case MULTISURFACETYPE:
case POLYHEDRALSURFACETYPE:
case TINTYPE:
case COLLECTIONTYPE:
return lwcollection_is_empty((LWCOLLECTION *)geom);
break;
default:
lwerror("lwgeom_is_empty: unsupported input geometry type: %s",
lwtype_name(geom->type));
break;
}
return result;
}
LWPOINT*
lwpoint_force_dims(const LWPOINT *point, int hasz, int hasm)
{
POINTARRAY *pdims = NULL;
LWPOINT *pointout;
/* Return 2D empty */
if( lwpoint_is_empty(point) )
{
pointout = lwpoint_construct_empty(point->srid, hasz, hasm);
}
else
{
/* Always we duplicate the ptarray and return */
pdims = ptarray_force_dims(point->point, hasz, hasm);
pointout = lwpoint_construct(point->srid, NULL, pdims);
}
pointout->type = point->type;
return pointout;
}