static void
lwgeom_collect_endpoints(const LWGEOM* lwg, LWMPOINT* col)
{
int i, n;
LWLINE* l;
switch (lwg->type)
{
case MULTILINETYPE:
for ( i = 0,
n = lwgeom_ngeoms(lwg);
i < n; ++i )
{
lwgeom_collect_endpoints(
lwgeom_subgeom(lwg, i),
col);
}
break;
case LINETYPE:
l = (LWLINE*)lwg;
col = lwmpoint_add_lwpoint(col,
lwline_get_lwpoint(l, 0));
col = lwmpoint_add_lwpoint(col,
lwline_get_lwpoint(l, l->points->npoints-1));
break;
default:
lwerror("lwgeom_collect_endpoints: invalid type %s",
lwtype_name(lwg->type));
break;
}
}
static LWMPOINT*
lwmline_locate_along(const LWMLINE *lwmline, double m, double offset)
{
LWMPOINT *lwmpoint = NULL;
LWGEOM *lwg = lwmline_as_lwgeom(lwmline);
int i, j;
/* Return degenerates upwards */
if ( (!lwmline) || (lwmline->ngeoms < 1) ) return NULL;
/* Construct return */
lwmpoint = lwmpoint_construct_empty(lwgeom_get_srid(lwg), lwgeom_has_z(lwg), lwgeom_has_m(lwg));
/* Locate along each sub-line */
for ( i = 0; i < lwmline->ngeoms; i++ )
{
LWMPOINT *along = lwline_locate_along(lwmline->geoms[i], m, offset);
if ( along != NULL )
{
for ( j = 0; j < along->ngeoms; j++ )
{
lwmpoint_add_lwpoint(lwmpoint, along->geoms[i]);
}
/* Free the containing geometry, but leave the sub-geometries around */
if ( along->bbox ) lwfree(along->bbox);
lwfree(along);
}
}
return lwmpoint;
}
static LWMPOINT*
lwpoint_locate_along(const LWPOINT *lwpoint, double m, double offset)
{
double point_m = lwpoint_get_m(lwpoint);
LWGEOM *lwg = lwpoint_as_lwgeom(lwpoint);
LWMPOINT *r = lwmpoint_construct_empty(lwgeom_get_srid(lwg), lwgeom_has_z(lwg), lwgeom_has_m(lwg));
if ( FP_EQUALS(m, point_m) )
{
lwmpoint_add_lwpoint(r, lwpoint_clone(lwpoint));
}
return r;
}
LWMPOINT*
lwmpoint_from_lwgeom(const LWGEOM *g)
{
LWPOINTITERATOR* it = lwpointiterator_create(g);
int has_z = lwgeom_has_z(g);
int has_m = lwgeom_has_m(g);
LWMPOINT* result = lwmpoint_construct_empty(g->srid, has_z, has_m);
POINT4D p;
while(lwpointiterator_next(it, &p)) {
LWPOINT* lwp = lwpoint_make(g->srid, has_z, has_m, &p);
lwmpoint_add_lwpoint(result, lwp);
}
lwpointiterator_destroy(it);
return result;
}
LWMPOINT *
lwmpoint_construct(int srid, const POINTARRAY *pa)
{
int i;
int hasz = ptarray_has_z(pa);
int hasm = ptarray_has_m(pa);
LWMPOINT *ret = (LWMPOINT*)lwcollection_construct_empty(MULTIPOINTTYPE, srid, hasz, hasm);
for ( i = 0; i < pa->npoints; i++ )
{
LWPOINT *lwp;
POINT4D p;
getPoint4d_p(pa, i, &p);
lwp = lwpoint_make(srid, hasz, hasm, &p);
lwmpoint_add_lwpoint(ret, lwp);
}
return ret;
}
static LWGEOM*
parse_geojson_multipoint(json_object *geojson, int *hasz, int root_srid)
{
LWGEOM *geom;
int i = 0;
json_object* poObjPoints = NULL;
if (!root_srid)
{
geom = (LWGEOM *)lwcollection_construct_empty(MULTIPOINTTYPE, root_srid, 1, 0);
}
else
{
geom = (LWGEOM *)lwcollection_construct_empty(MULTIPOINTTYPE, -1, 1, 0);
}
poObjPoints = findMemberByName( geojson, "coordinates" );
if ( ! poObjPoints )
{
geojson_lwerror("Unable to find 'coordinates' in GeoJSON string", 4);
return NULL;
}
if( json_type_array == json_object_get_type( poObjPoints ) )
{
const int nPoints = json_object_array_length( poObjPoints );
for( i = 0; i < nPoints; ++i)
{
POINTARRAY *pa;
json_object* poObjCoords = NULL;
poObjCoords = json_object_array_get_idx( poObjPoints, i );
pa = ptarray_construct_empty(1, 0, 1);
parse_geojson_coord(poObjCoords, hasz, pa);
geom = (LWGEOM*)lwmpoint_add_lwpoint((LWMPOINT*)geom,
(LWPOINT*)lwpoint_construct(root_srid, NULL, pa));
}
}
return geom;
}
static LWMPOINT*
lwmpoint_locate_along(const LWMPOINT *lwin, double m, double offset)
{
LWGEOM *lwg = lwmpoint_as_lwgeom(lwin);
LWMPOINT *lwout = NULL;
int i;
/* Construct return */
lwout = lwmpoint_construct_empty(lwgeom_get_srid(lwg), lwgeom_has_z(lwg), lwgeom_has_m(lwg));
for ( i = 0; i < lwin->ngeoms; i++ )
{
double point_m = lwpoint_get_m(lwin->geoms[i]);
if ( FP_EQUALS(m, point_m) )
{
lwmpoint_add_lwpoint(lwout, lwpoint_clone(lwin->geoms[i]));
}
}
return lwout;
}
