408 double *xpts, *ypts, *zpts, *mpts;
410 int *shpparts, shppointtotal = 0, shppoint = 0;
413 shpparts =
malloc(
sizeof(
int) * lwmultilinestring->
ngeoms);
416 for (i = 0; i < lwmultilinestring->
ngeoms; i++)
419 LWDEBUGF(3,
"Total number of points: %d", shppointtotal);
422 xpts =
malloc(
sizeof(
double) * shppointtotal);
423 ypts =
malloc(
sizeof(
double) * shppointtotal);
424 zpts =
malloc(
sizeof(
double) * shppointtotal);
425 mpts =
malloc(
sizeof(
double) * shppointtotal);
428 for (i = 0; i < lwmultilinestring->
ngeoms; i++)
431 shpparts[i] = shppoint;
437 xpts[shppoint] = p4d.
x;
438 ypts[shppoint] = p4d.
y;
439 zpts[shppoint] = p4d.
z;
440 mpts[shppoint] = p4d.
m;
442 LWDEBUGF(4,
"Linestring %d - Point: %g %g %g %g", i, xpts[shppoint], ypts[shppoint], zpts[shppoint], mpts[shppoint]);
POINT4D getPoint4d(const POINTARRAY *pa, uint32_t n)
#define LWDEBUGF(level, msg,...)
SHPObject SHPAPI_CALL1 * SHPCreateObject(int nSHPType, int nShapeId, int nParts, const int *panPartStart, const int *panPartType, int nVertices, const double *padfX, const double *padfY, const double *padfZ, const double *padfM);SHPObject SHPAPI_CALL1(*) SHPCreateSimpleObject(int nSHPType, int nVertices, const double *padfX, const double *padfY, const double *padfZ